title: How to change feature's geometry?

text: Hy, I'm creating my layer like this:               crs = core.QgsCoordinateReferenceSystem(4326, core.QgsCoordinateReferenceSystem.EpsgCrsId)     uri = "Polygon"     uri += '?crs=' + crs.authid()     uri += '&field=id_tsk:integer' # IDX_ID     uri += '&field=name:string' # IDX_NAME     uri += '&index=yes'          layer = core.QgsVectorLayer(uri, 'name', "memory")          core.QgsMapLayerRegistry.instance().addMapLayer(layer)      When I look to the capabilities of the dataProvider (in QGis Python console) i've got this:               >>> layer.dataProvider().capabilitiesString()     PyQt4.QtCore.QString(u'Add Features, Delete Features, Change Attribute Values, Add Attributes, Delete Attributes, Create Spatial Index, Fast Access to Features at ID, Change Geometries')      My data are contains in a QAbstractItemModel of my own. So when the model data changed I want to update only element / geometry that has changed and not the whole layer (witch might be big).               self._layer.setReadOnly(False)     feat_2_del = set()     feat_2_add = []     feat_2_chg = {}     # +1 because we want the bottomRight row to be selected     for row_num in xrange(topLeft.row(), bottomRight.row() + 1):         idx_geo = self._model.index(row_num, self._col_geo)         wkt_v = idx_geo.data(QtCore.Qt.DisplayRole)              obj = self._get_object_from_idx(idx_geo)         feat = self._object_map.get(obj.id)              geom = None         if not wkt_v.isNull():             wkt = wkt_v.toString()             if wkt.size() > 0:                 # some object do not have any geometry defined                 geom = core.QgsGeometry.fromWkt(wkt)              old_exist = feat is not None         new_has_geo = geom is not None              if old_exist and new_has_geo:             # we just change the old feature             feat_2_chg[feat.id()] = geom         elif old_exist and not new_has_geo:             # we remove the old feature             feat_2_del.add(feat.id())             del self._feature_map[feat.id()]         elif not old_exist and new_has_geo:             # we create a new feature             feat = core.QgsFeature(id(obj))             feat.addAttribute(IDX_ID, QtCore.QVariant(obj.id))             feat.addAttribute(IDX_NAME, QtCore.QVariant(obj.name))             feat.setGeometry(geom)                  self._object_map[obj.id] = feat             self._feature_map[feat.id()] = obj             feat_2_add.append(feat)          self._layer.dataProvider().deleteFeatures(feat_2_del)     self._layer.dataProvider().addFeatures(feat_2_add)     self._layer.dataProvider().changeGeometryValues(feat_2_chg)     self._layer.commitChanges()     self._layer.setReadOnly(True)     self._iface.mapCanvas().refresh()      When I add a geometry information to an object (addFeatures), the layer is correctly updated (with the according geometry). But If I change the geometry (changeGeometryValues) or remove a feature (deleteFeatures) nothing is done. But my dataProvider should handle it. What am I doing wrong ?


title: Source for all census tables available at each census geography/summary level

text: I'm familiar with using factfinder2 and generally know my way around census products. But one thing that I've never found is a document that has all the tables from major census products like the dicennial census or ACS and the level of census geography at which it's available. Anyone ever come across such a thing?


title: Nearest Neighbor problem in Postgis 2.0 using GIST Index (<-> function)

text: I'm trying to use Postgis 2.0 new function <-> (Geometry Distance Centroid) in order to calculate, for each row of my table (cosn1), the distance to the nearest polygon of the same class. I was trying to use the following code:               WITH index_query AS (       SELECT g1.gid As ref_gid, ST_Distance(g1.the_geom,g2.the_geom) As ENN             FROM "cosn1" As g1, "cosn1" As g2            WHERE g1.gid <> g2.gid AND g1.class = g2.class         ORDER BY g1.gid, g1.the_geom <-> g2.the_geom)      SELECT DISTINCT ON (ref_gid) ref_gid, ENN          FROM index_query     ORDER BY ref_gid, ENN;      But then I realize the warning: > Note: Index only kicks in if one of the geometries is a constant (not in a > subquery/cte). e.g. 'SRID=3005;POINT(1011102 450541)'::geometry instead of > a.geom Meaning that the Index wont be used at all, and the query will take almost the same time as before using:               SELECT DISTINCT ON(g1.gid)  g1.gid As ref_gid, ST_Distance(g1.the_geom,g2.the_geom) As ENN             FROM "cosn1" As g1, "cosn1" As g2            WHERE g1.gid <> g2.gid AND g1.class = g2.class         ORDER BY g1.gid, ST_Distance(g1.the_geom,g2.the_geom)      Can anyone point me a workaround that allows me to improve performance of my query? Thank you very much.


title: Python Script to Create Table Using Data Driven Pages

text: I am trying to convert some python script to display a table (based on a dbf) of features located within a specific data driven page. So far, I have the script successfully able to refresh the map to the specific table, but it does not update the table. I have it set-up as three text boxes that should be updated with three specific fields when the user runs the script from ArcToolbox. Any suggestions on why my table is not updating?               import arcpy, sys, os          #Reference current MXD     mxd = arcpy.mapping.MapDocument("current")          #Get input parameter     Name = arcpy.GetParameterAsText(0)          #Reference  data frames     mapatlasDF = arcpy.mapping.ListDataFrames(mxd, "Layers")[0]     locatorDF = arcpy.mapping.ListDataFrames(mxd, "Locator Map")[0]          #Reference appropriate layers     atlasLyr = arcpy.mapping.ListLayers(mxd, "PinalCreekMapAtlas_HalfMile", mapatlasDF)[0]     locatorLyr = arcpy.mapping.ListLayers(mxd, "Locator Map", locatorDF)[0]     atlasoutlineLyr = arcpy.mapping.ListLayers(mxd, "Map Atlas Outline", locatorDF)[0]          #Reference layout elements by calling ListLayoutElements      for elm in arcpy.mapping.ListLayoutElements(mxd):       if elm.name =="Table1Column1": tab1Col1Txt = elm       if elm.name =="Table1Column2": tab1Col2Txt = elm       if elm.name =="Table1Column3": tab1Col3Txt = elm          #Reference the Data Driven Page object     ddp = mxd.dataDrivenPages          #Set the current page to be the one selected in the script tool     arcpy.AddMessage(Name)          pageID = mxd.dataDrivenPages.getPageIDFromName(str(Name))     mxd.dataDrivenPages.currentPageID = pageID          #Set the appropriate definition queries     atlasLyr.definitionQuery = "Name = '" + Name +  "'"     locatorLyr.definitionQuery = "Name = '" + Name +  "'"     atlasoutlineLyr.definitionQuery = "Name <> '" + Name +  "'"          #Update Sheet Index data frame     arcpy.SelectLayerByAttribute_management(locatorLyr, "NEW_SELECTION", "\"Name\" = '" + Name + "'")     locatorDF.panToExtent(locatorLyr.getSelectedExtent())          #Reference Affected Parcels table and select appropriate records     parcelTable = arcpy.mapping.ListTableViews(mxd, "AffectedParcels")[0]          #Build query and create search cursor to loop through rows     parcelFieldValue = "Page " + Name     queryExp = "\"MapPage\" = '" + parcelFieldValue + "'"  #e.g., "MapPage" = 'Page 01'     parcelRows = arcpy.SearchCursor(parcelTable.dataSource, queryExp)          #Clear all table text values     tab1Col1Txt.text = " "; tab1Col2Txt.text = " "; tab1Col3Txt.text = " "          #iteate through each row, update appropiate text     count = 0     for row in parcelRows:       if count < 30: #Table1 - static position         tab1Col1Txt.text = tab1Col1Txt.text + row.getValue("OwnerName") +"\n"         tab1Col2Txt.text = tab1Col2Txt.text + row.getValue("APN") + "\n"         tab1Col3Txt.text = tab1Col3Txt.text + row.getValue("LengthTrail") + "\n"       if count ==30:           arcpy.AddMessage("Table Overflow") #The code could be reworked to show the last 90 records       count = count + 1          arcpy.RefreshActiveView()     arcpy.AddMessage("PROCESS COMPLETED")


title: Can QGIS perform graphic clipping?

text: In ArcGIS for Desktop there is an option to graphic clip inside the data frame - this clips the feature you're viewing without actually clipping the feature into a new feature. Does Quantum GIS have a tool to do this?


title: Polyline will not turn into a feature layer via MakeFeatureLayer_management (arcpy)

text: I've got a python script that is connecting a sequence of points to make lines (script is based on answers to my previous question). To attempt to make subsequent processing more efficient, I would like to do analysis/processing to each polyline as the points are connected, rather than storing them cumulatively in a shapefile (or feature class) and then reading them back into the script later. So the revised code (just the functional looping part, for the sake of brevity) is basically this:               for start in start_cursor:         startFeature = start.getValue(shapefieldname)         pnt1 = startFeature.getPart()         for end in end_cursor:             endFeature = end.getValue(shapefieldname)             pnt2 = endFeature.getPart()             point.X = pnt1.X             point.Y = pnt1.Y             array.add(point)             point.X = pnt2.X             point.Y = pnt2.Y             array.add(point)             polyline = arcpy.Polyline(array)             arcpy.MakeFeatureLayer_management(polyline, "sightline_lyr")      However, **MakeFeatureLayer_management** doesn't seem to exactly be working. Any subsequent functions I want to apply (e.g., `FeatureToRaster_conversion`) to that new feature layer that was supposed to be created from that polyline has no data. If I add the new feature layer to the data frame while the script is running, it shows up in the table of contents in the appropriate data frame, but with the little red exclamation mark that indicates ArcMap can't find the data for that feature layer. Is an arcpy.Polyline able to be made into a (temporary) feature layer? If not, how can I make a feature layer such that I can do more spatial analysis in further steps? (Short of "make a shapefile or feature class, then read it back in," preferably...)


title: Weird raster visualisation

text: Using the hillshade tool (Arcmap 10) on a 30m SRTM raster downloaded from NED the result is the one in the picture. Why are there those crossed lines all over the raster? How could I remove them? ![enter image description here](http://i.stack.imgur.com/JpP5W.jpg) Thank you, Michele


title: In MapBasic, how do you place a user selected string from a FileOpenDlg() statement into an EditText field for use by a sub procedure?

text: I have created a custom menu that displays a dialog box when clicked on. I have made the "Browse" button to call a FileOpenDlg() statement. I want the file that the user selects to populate in the EditText space. I would then like to have my app read that string and store it, so that it can be used in a sub procedure when the user clicks run. Any suggestions? For Example: I want it to behave like a box that asks you to pick a file to upload. The browse button is used to browse for the file, and I want the "run" button to behave like an "upload" button. I was unable to post the picture on here of my dialog box, but that is the concept i am going for. Here is the code for the dialog box and the handler:                   Sub Leads         Dialog     Title "Leads"         Width 198 Height 80     Control StaticText         Title "Select Leads File"         Position 4, 16     Control EditText         ID "2"     Control Button         Title "Browse"         Calling PickLeadsFile         ID "22"         Position 152, 15     Control StaticText         Title "Run Leads Update?"         Position 4, 52     Control OkButton         Title "Run"         ID "222"         Position 74, 48     Control CancelButton         Position 118, 48              End Sub                Sub PickLeadsFile         Dim Leads_filename As String     Leads_filename = FileOpenDlg("C:\Users\dykesb01\Desktop\Test","","XLS","Open Leads File")         Alter Control 2         Title From Variable Leads_filename         End Sub


title: How to export a .sid image as high DPI .tiff or .bmp from QGIS?

text: I have 2 50 cm resolution areal photographs in .sid format. How can I export these layers as a regular image format (like tiff or bmp) which I can then open in any simple picture viewer program? I've tried saving a map, but that reduces the image resolution. I need to resulting image to retain it's original resolution.


title: How to Configure Geoserver WFS to work with ArcGIS for Desktop?

text: I am trying to figure out how to get a WFS service working with ArcMap. I have Geoserver up and running and I am able to add a WFS layer in QGIS but not ArcMap. I have tried to use the Data Interoperability Tool in ArcMap to connect to my server but there are no layers displayed in select Feature Types dialog. I am able to connect to a WFS server so I'm assuming that I have overlooked some detail about how ArcMap wants Geoserver to be configured.


title: How to implement Spiral sequence rule in GML grids?

text: I am working with gml:RectifiedGridCoverage and currently I am introducing support for the different sequence rules defined inside gml:coverageFunction. These are the supported sequence rules in GML 3.2:   * Linear   * Boustrophedonic   * Cantor-diagonal   * Spiral   * Morton   * Hilbert I have established how to move along the grid when mapping between the range set and the domain set using the previous rules, anyway there is an exception, I haven't found anything about how the spiral sequence is supposed to work, I think I have an idea about how the spiral sequence should work but that would be only a guess and I would like to know what is the "official" way to do it. Does anybody know how the spiral sequence should be implemented? I have covered the other sequences, but if somebody could point to other resources that could help me to reinforce what I already made then that would be awesome.


title: GRASS module g.mapset can't create new mapsets from outside

text: I use pyscripter to control GRASS from outside. I want to create a new mapset named **water**. Here's my script to generate new mapset,               import os     import sys          #set up GRASS environment variables     sys.path.append(os.path.join(os.environ['GISBASE'], 'etc', 'python'))     import grass.script as g     import grass.script.setup as gsetup     gisbase = os.environ['GISBASE']     gisdb = 'C:\Users\Heinz\Documents\grassdata'     location = 'newLocation'     mapset = 'TC'     gsetup.init(gisbase, gisdb, location, mapset)          #create new mapset in the location     g.run_command('g.mapset', flags = 'c', mapset = 'water', location = 'newLocation', gisdb = 'C:\Users\Heinz\Documents\grassdata')          #check result     print g.read_command('g.mapsets', flags = 'l')      and the output,               *** Remote Interpreter Reinitialized  ***     >>>      PERMANENT TC watershed_analysis      The script runs fine but with no new mapset created, what mistake did I make？ By the way, this script takes much time to run, but if I remove the line which creates new mapsets, the script runs faster, so maybe that line is the reason why GRASS fails to generate new mapset. * * * I tried to type those lines in the python shell of the GRASS layer manager,               >>>import grass.script as g          >>>g.read_command('g.mapsets', flags = 'l')     PERMANENT TC water water_1 watershed_analysis          >>>g.run_command('g.mapset', flags = 'c', mapset = 'water_2', location = 'newLocation', gisdb = 'C:\Users\Heinz\Documents\grassdata')     0     >>>g.read_command('g.mapsets', flags = 'l')     PERMANENT TC water water_1 water_2 watershed_analysis      It can generate a new mapset, and it proves that my script is OK. I can create new mapsets inside the GRASS layer manager, but can't from outside, are there anyone successfully do this from outside? I am using GRASS 6.4.3, python 2.7 under Windows 8.1


title: How to get a the minimum value of a raster layer in qgis with python?

text: I'm very new user of QGis and I work with landsat 5 images. In this framework, I would like to make an atmospheric correction subtracting the minimum value present in each band (1 to 4). My problem is that the minimum value present in the band is zero, but 0 is the background or none data. I think that is possible to use some python code (using Qgis console) for this process, where I can to find the minimum value in the band (ignoring values equal zero). Thanks! Marcelo               (For i=1 To n;          For j=1 To m;             If myraster(i,j) > 0 Then (                ND = myraster(i,j);                  If ND < minND Then                            minND = ND);                   Next j             Next i     )


title: Extracting elevation from a TIN file with a point shapefile overlay

text: I have a TIN file that was built off of contour polyline shapefile (I think). I have a point shapfile with a few points. How can I extract elevation of the TIN file where the points occur? I am using ArcMap with python automation. Any free software would be welcome as well. Perhaps a transforming of the TIN to another file type is needed. Thanks.


title: Styling river charts with TileMill - issue with filling a polygon made of disjoint line segments

text: I am trying to style inland river charts produced by the Army Corps of Engineers. I imported the vectors into Postgres, and I’m using TileMill to style the data. I am having trouble correctly filling the inner contour for the river, because the contour is made of line segments that are broken by features like bridges, and the contour fill spills out. How can I fill the contour properly and avoid this issue? Some sort of processing in Postgres? Is this a solved problem… a good way to make gapped line segments into polygons for filling? You can see a screenshot of the issue here: ![enter image description here](http://i.stack.imgur.com/dt0ev.png)


title: How to handle updates to pgRouting network in webmapping environment?

text: I'm building an application for a client of mine that will need to create routes in a PostgreSQL database. So far so good, but I need to support edits and changes to my network and I need to support bulk imports. I've read the tutorials and the examples available on the web and the routing part is clear for me. My main concern is perfomance. If I try to update the network at each update, this will be incredibly slow. Here's the workflow I'm thinking for creating new edges:   1. User inserts a new edge (let's call that A);   2. System looks for existing edges that touches or crosses this new edge (edge B);   3. System splits edge B (if touches) and/or splits edge A (if crosses); This will trigger another event at the database (recalculating the new nodes for each of B splitted parts B and B');   4. System looks for existing nodes for edge A;   5. If nodes found, update the feature to maintain the network topology (from to node ids);   6. If no nodes are found, create them and update the necessary features; Of course there will be similar triggers for updates to the edges (this may be even more complex than the above) and deletes. For a small number of users it may be acceptable, but for bulk imports this is unlikely to work (I have three of four spatial queries being run for each update). My bulk process will be run using async tasks (Celery). As whuber always says, if it's too slow, blaim the algorithm. I have a similar setup with ESRI software for a customized network scheme and works well, but it's not the fastest kid on the block. I have two ways of doing this: in the application layer or in the database. The database is a good choice IMO but I may leak the domain layer to the database and I REALLY don't want to do that. The application is being developed with GeoDjango. What are your thoughts on this? Can you guys comment on the algorithm? This is a more django-specific question: can I query memcached entities using spatial queries?


title: ArcGIS 10 Add-Ins

text: What are the best Add-Ins for ArcCatalog and ArcMap v10? XRay for ArcGIS 10 is very useful, what else is out there? Thanks, Seth


title: Is there a way to update the attribute table of a feature class with a query?

text: I have a layer with multi-polygon feature classes for states and provinces. The shape area for some features classes is < 0 causing the label to be position in odd locations. It seems that the calculation is generally correct, and that if I could update all shape area in a single query I could fix this quickly. Currently I am querying to identify the negative area features classes, and editing a point, which resets the area and then it labels correctly. However it would be much quicker to be able to run a batch update on the value.


title: What is the modern equivalent to ArcInfo CLEAN?

text: What is the arcgis equivalent to the Arcinfo Workstation CLEAN command? Or perhaps better rephrased as, how does one clean topology without installing arcinfo workstation? Assume an arcinfo license level is present, just not the `Arc:` command line environment. **Update:** A key part of the question is "without installing arcinfo workstation". I blame the confusion on ESRI's decision to use the word "ArcInfo" for two distinct things, and then decouple them. In common usage ArcInfo can refer to a) arcinfo workstation, and b) arcinfo license level. You can't have workstation without an arcinfo license, but you can, and increasingly more and more people do, have an arcinfo license without workstation present -- which means the answers suggesting geoprocessing clean tool can't be used. And Kirk is right, this question can also be thought of as _"What is the modern equivalent to a polygon coverage (and how do you clean it)?"_


title: GeoServer Limit SLD editing functionality to certain users in GeoExplorer

text: How would I go about preventing user X from editing layer styles? I have created a user in GeoServer that if they login they only see the layer preview screen. However as I am using the Boundless Suite (OpenGeo) if they add a new site in the GeoExplorer,they login and add a new point/ However the same use can also edit the style for that layer. I have hidden all the other layers, however this layer must be visible. Is there a way to either remove the style editor from GeoExplorer? or prevent users from editing styles from the GeoServer side? I did try creating a edit feature map using the OpenGeo SDK however found that adding the user login functionality pretty complex with little easy to follow documentation. I had all the functionality I required, as I could turn off the edit style feature, but then anyone who managed to access the layer would be able to add a new point.


title: OpenLayers - Can't remove zoom control buttons, other controls simply lay on top

text:
I have a simple OpenLayers script with a Google map as the base layer. By default the map has a + and - zoom control button, which I can't find how to get rid of. I've added an empty controls array to the map options, but it still appears. If I add any other controls, for example PanZoomBar, it just sits ontop of the + and - buttons, semi-obscuring them. The effect is the same whether I add PanZoomBar in the controls array or with map.addControl. Here is what it looks like by default: ![enter image description here](http://i.stack.imgur.com/zWNZP.jpg) And with added PanZoomBar: ![enter image description here](http://i.stack.imgur.com/byrKN.jpg) And here's my code:               <!DOCTYPE html>                   <title>Map test</title>
Map Test

                                   Can anybody help me get rid of the + and - buttons?


title: Is there any methodology or technology to post process the Mobile/Tablet GPS data

text: We have developed an application which is running on android OS for collecting the Lat/Long value of any feature and it runs in any type of android Tablet/Mobile but the problem is that the accuracy is very poor (because the accuracy is depends upon the GPS receiver of that device). Is there any methodology or technology to post process the GPS data?


title: Methods for geotagging/displaying text documents

text: I'm wondering what methods folks are familiar with for geotagging/displaying text documents, or the text from a document. For example, it's easy enough to geotag a photo and display it in ArcGIS or Google Earth. But what if I have an extended textual description of a place? Is it possible to geotag that, export it as a shapefile, and display that text in a GIS? One possibility seems to create a text field, but that seems limited in characters. Any suggestions would be much appreciated.


title: gdalwarp to EPSG 3857 of a Bing tile

text: I'm quite new to GDAL and to the remote sensing world so sorry in advance if I say something silly :-) What I would like to do is to assign from scratch the EPSG 3857 System Of Reference to a Bing tile (a not georeferenced jpg) downloaded from the Bing server. I've already computed the coordinates [x, y] this tile should have in the 3857 SoR and create a world file from them, but how can I now project my tile to 3857 by using this .tfw? ![enter image description here](http://i.stack.imgur.com/jj53t.jpg) Thank you so much in advance! Martina


title: GWR Result Issues

text: I am working with GWR to identify if there are any relationships between 'extreme' events and disease incidence. The use of GWR has been employed because I am working in a region where local differences occur for for selected variables. It is slightly different to trying to find all the variables that explain the incidence rates, rather I am looking at specific variables to see if they do impact or not. I am getting good results for the St. Residuals, local r2 and coefficient rasters, however I am concerned about the results in the results box. I've read the literature about these results and it says for the residualsquares and sigma the smaller the value the better. My results for these two measures aren't what I would consider small. Is this vital or are the results I actually map more important? The AICc value is lower for GWR than OLS. I feel positive about the st. residuals, local r 2 and coefficient rasters but don't want to put these in results in my work if actually underlying those the model is no good. Also, I know it is advised to start with OLS to find a good model, but is this essential if I am not looking to explain the incidence rates fully, rather I am looking to see if the selected variables do impact the incidence rate or not. Thanks


title: How to calculate elevation differences between points and terrain model?

text: I have a points layer with a Z (elevation) field and a terrain model (points layer or polygons layer -triangles- or eventually raster). I'm looking for a way to select the entities of my points layer that are at least 10m under the terrain model. Any way to do this in QGIS 1.8 ?


title: Creation of a GIS Competency matrix for employers

text: I am in the process of developing a technical GIS competency matrix for new employees. The matrix will be used not only to assess new employees but also to monitor employee development. The matrix should contain everything from general GIS concepts to advanced server-side GIS and Web-GIS development. Has anyone had any experience with such a matrix or can shed any light on how this would look? My first attempt would look this this.   1. General GIS concepts: CRS models and transformations, Formats, spatial analysis, georeferencing, heads-up digitizing, symbology, colour   2. Spatial databases: ESRI, Postgresql/PostGIS, Mysql   3. WebGIS: O-O Programming, Openlayers, Serverside GIS, Gdal, Geoserver/Mapserver/Deegree   4. General programming: Python, JavaScript, HTML/CSS, JQuery, PHP   5. Linux Administration: Shells, scripting, commands, administration, monitoring I am also interested in finding out what kind of weighting these various criteria should carry. My personal bias is that I think I place more value upon the contextual nature of the GIS-Work than the Technical side, i.e that I find it more important to first understand the data before one develops the tools, rather than first to achieve the technical prerequisites and then try to understand what the data or tool actually represents. Due to this, I would rather employ a planner who learnt GIS-Skills, than a programmer who has later developed GIS skills. Obviously this only applies to GIS projects of a strong contextual nature. If I wanted a GIS-Administrator to only monitor JAVA projects on Linux, then i´d obviously prefer an IT professional. This is only a first attempt, so any comments are gratefully received. Thanks, Rob


title: Nullable fields in feature service will not save with null values

text: I created my file gdb schema in Diagrammer, which included setting up several coded value domains, and applying them to fields in the feature class. All the fields have their `Required` property as the default `No` and `IsNullable` as the default `True`. Once the file gdb is created from the xml document and the data loaded, the feature class is copied into SDE. An mxd was set up and published as a feature service with mobile capabilities. During testing the mobile project in ArcGIS mobile, I noticed that when I add a new feature, or edit an existing one's attributes, all the fields with domains applied give an error:`Field value is not defined.You must enter a value to continue.` When I edit the SDE feature class in ArcMap, it allows me to save edits with new features having null fields. I am unsure why the mobile project will not allow me to save it, as I have created various projects in the past using an identical workflow, and I have able to save features without entering a non- null value for fields with attribute domains. I see a bug report was submitted for feature services on ArcGIS Online doing the same thing, so it's possible that this is a feature service issue. I am using ArcGIS Mobile/Desktop/Server 10.


title: Point symbology in SLD file

text: I would like to create a style file for a layer in Geoserver. I have point geometries and the SLD cookbook does give nice examples for the simple ways to create your symbology. However, I would like to combine the attribute dependent way and the scale dependent way. My points need to be colored depending on some attribute value and change size depending on the zoom level. Does anybody have an example on how to do this or even some link where more complicated version of SDL files are available? Cheers, Tom


title: Postgis template1 database missing

text: I want to create a template_postgis for connecting postgis via qgis So I follow the step on following page how to create a new postgis database for usage with QGIS under windows But when I follow the steps there's a problem happened               createdb: database creation failed: ERROR:  template database "template1" does not exist      So I google it for answer but no one can help me to create a template1 database! Any suggestion?


title: Finding Help For Setting Up QGIS Web Client in Windows 7

text: Is there any procedures or installation guide of QGIS Web Client for Windows 7 OS? Because I'm really having a hard time setting up the Web Client in Win7 and I don't know any Linux-based commands that's featured in https://github.com/qgis/qgis-web-client. My QGIS Server is running well and I'm using QGIS 1.8. Thanks in advance


title: Upgrade Postgres/Postgis

text: I have Postgres 8.4 and postgis 1.5.3 installed on centos6. Now I need to upgrade both to postgres 9.x and postgis 2.x. Is there somewhere a good tutorial how to upgrade without removing the old versions? I need to be able to run both versions beside each other. Thanks


title: Select all "NULL" values in the attribute table?

text: I have a quite simple syntax question. I want to select ALL my null values in an attribute table of a feature class, not just the values in a single field. So the question would be: what is the valid sql statement that could get me all the null values? Thank you!


title: How to modify height model values outside of road network

text: I have a height model of a quite big area saved in Esri ASCII grid files and a road network of this same area in a shapefile. I'd like to convert every height model value outside the road network to -999.999 (NODATA_value in my .asc-files), as I won't be needing them. This would reduce amount of data massively as I will later use only real values and forget those NODATA_values. Is there any reasonable way to do this, for example in QGIS or ArcMap? Would it even be possible to have a few meters wide buffer around the roads, just to be sure that there will be every needed height value in the end?


title: How to check map tiles are loaded in OpenLayers?

text: How to check if the map tiles are loaded? I need to know whether all map tiles are loaded on to the OpenLayers map as html contents.... Is there any method to find out or function in openlayers to find out for how much percentage the map tiles are loaded.......


title: Programmatically change layer position in the table of contents (QGIS)

text: I am trying move/position a layer to the bottom of the TOC by using addTopLevelItem on all other layers. I unsure how to reference the layers as QTreeWidgetItem's.               def bg_bottom_level_item(self):         layers = self.canvas.layers()         layerlist = []         for layer in layers:             layerID = str(layer.id())             if "Auto_Background" not in layerID:                 layerlist.append(layer)             else:                 pass         for layer in layerlist:             leg = qgis.utils.iface.mainWindow().findChild(QTreeWidget, 'theMapLegend')             item = QTreeWidgetItem(layer)             leg.addTopLevelItem(item)      Error: `QTreeWidgetItem(QTreeWidgetItem): argument 1 has unexpected type 'QgsRasterLayer'` Or does anyone have another solution for moving layers around in the TOC?


title: How to Convert Interpolated Map to Excel using ArcGIS for Desktop?

text: I have interpolated some SO2 data using IDW and have calculated means on county level using Zonal Statistics. I would like to convert this data now into an excel file. I already got the hint, that I should use the Conversion Tool and Table to Excel. However, I can only convert the two tables I have imported from shape files. That is the map with the counties and the data with SO2 values and measuring stations. Do I have to convert the values I obtained from Zonal Statistics? I hope I was able to describe my problem sufficiently. I am using ArcMap.


title: Reading a raster map cells values in C#?

text: Suppose I have a raster map file (TIFF, TGA, JNG... etc) loaded with geo-info manifested by the cells values. I need to read the values of the raster map cells "programatically" in a C# program. P.S: I don't want to use "heavy" tools or libraries (like ArcGis and ArcObjects) for such a simple task.


title: ArcGIS 10.2 Combine Two Feature Classes in One File Geodatabase with Attachment (Pictures)

text: I am working in ArcGIS 10.2, and have point data collected from 2 separate tablets with the same database schema. Both of these files have attachments (pictures) with them. I am trying to combine these two files and can't seem to find a way to keep the attachments when I do combine them. I have tried creating a new file geodatabase and exporting the files as XML's. I have tried loading the data into a new file geodatabase and nothing seems to work. Has anyone run into this, or know a way to combine these two feature classes?


title: Add to geoprocessing environments when creating ArcGIS .net Add-in

text: I would like to add a new environment to the Geoprocessor Environments from an ArcMap.net Add in. I try adding the env to the GPEnvironmentManager with the following code but when ArcMap opens and I look at the environments, its not there.                           IGPEnvironment secTolEnv = new GPEnvironmentClass();                 secTolEnv.Category = "MyCategory";                 secTolEnv.Name = "_sectiontolerance";                 secTolEnv.DisplayName = "Section Tolerance";                 secTolEnv.DataType = new GPLinearUnitTypeClass();                 secTolEnv.Value = new GPLinearUnitClass                 {                     Value = 360,                     Units = esriUnits.esriFeet                 } as IGPValue;                      IGPEnvironmentManager GPEnvManager = new GPEnvironmentManagerClass();                 GPEnvManager.PersistAll = true;                 GPEnvManager.AddEnvironment(secTolEnv);      I am calling this while initializing the add in extension. Its getting called when ArcMap opens and I can step through it fine. Is this possible? I was thinking maybe I need to get the GPEnvironmentManager from ArcMap itself but I have not figured out how to do this yet. Any ideas are appreciated.


title: Create custom locator file - what license level of ArcDesktop do I need?

text: I would like to create custom locators in ArcDesktop to circumvent this problem: Why doesn't US One Address locator style have city, state or zip code fields? So I am going to buy a copy of ArcDesktop, but I think there are 3 levels, so I'm not sure which one to buy? ArcView, ArcEditor or ArcInfo. Which level allows you to create custom locators, aka .lot files?


title: ESRI Maps with Ravis Visual graph

text: How can we integrate ARCGIS maps with Ravis Birds eye's Visualgraph ? My intension is to place the ravis visual graph over the map and the nodes of the visual graph must be placed on their specific location on the map ... Can anybody help me in this regard ? Thanks in advance .


title: Workflow for DXF interchange with QGIS

text: What is the best method of importing/exporting dxf files with QGIS ? I work in an engineering field and we use CAD and specialist design programs, and often go through several iterations during a project. I am having trouble interchanging data between QGIS and dxf formats. I export several layers from QGIS in DXF format and load them into cad. Essentially I xref the exported DXF files into cad, and then draw over the top of them in cad, as I find it quicker and more acruate. The files line up ok, but need to be scaled and transformed to match the coordinates in QGIS. I then detach the xref's, save the new DXF with line work on a single layer. The DXF file coordinates match the qgis coordinates. I cant seem to get the new DXF to import back into QGIS. I would like to be able to export and import a particular data table several times as the project develops. **Screenshot 1** - I exported the pink polygons and edited them in CAD. I then imported the DXF file back in but nothing came in - see "entities". ![enter image description here](http://i.stack.imgur.com/CaFOK.jpg) **Screenshot 2** - Attribute table of imported layer. ![enter image description here](http://i.stack.imgur.com/pP3JV.jpg) Thanks


title: update one shapefile attribute table with values from another dbf table

text: I try to explain here what I am trying to do: I have one shapefile and one independent dbf table with the same fields. In the dbf table all the fields are populated but in the shapefile attribute table just one, lets name it "OneField". What I want to do is to check that the values from "OneField"(Shapefile) are the same as the values in "OneField"(dbf table) and if so, to populate the remaining empty fields in the shapefile attribute table with the ones in the independent dbf table. At the moment I am trying just to copy the values from the independent dbf table to the shapefile attribute table but I am stuck.(when I run this code I get a message that pythonwin stopped working and nothing happens to the tables). Can you give me a hand please? Here is the code:               import arcpy          table = "link/to/table.dbf"     fc = "link/to/shapefile.shp"          # Create a search cursor          rowsTable = arcpy.SearchCursor(table)          # Create an update cursor          rowsFc = arcpy.UpdateCursor(fc)           for row in rowsTable:         row = row.getValue("OneField")         valueTable = row         for row in rowsFc:             row = row.setValue("OneField", valueTable)             rowsFc.updateRow(row)             row = rowsFc.next()         row = rowsTable.next()              del row, rowsFc, rowsTable      Thank you very much


title: How to handle TIME in GIS?

text: Is there any way or any idea on how to handle Time attribute in GIS? For example, which are the ways to specify that an object is at a place at a given time interval and them it is in another place in another time interval, let's say, like a political boundary that has changed.


title: Would you consider online geocoding a breach of privacy?

text: Suppose I have a bunch of addresses of individuals participating in a certain study (most likely - health related, where privacy and ethical considerations are always important issues). Nowadays, providers like Google or Yahoo offer decent results in terms of positional accuracy. The North American Association of Central Cancer Registries (NAACCR) lists such options in their 'Geocoding Best Practices: Review of Eight Commonly Used Geocoding Systems' and 'A Geocoding Best Practices Guide' guides. Cinnamon and Schuurman (2010) for example used BatchGeocode service as a part of their tool to investigate injuries in low resource setting. Would you consider geocoding such addresses using online services, like Google Maps or OpenStreetMap a breach of privacy? PS1 possibly related question. PS2 recent article in Epidemiology (one of the leading, peer-review journals in the field) published short communication detailing instructions on how to geocode using Google Maps & Places APIs. Interestingly, not a word about security/privacy was mentioned...


title: Editing with auxiliary features (lines, points etc.) in Arcgis 9.3.1 - How?

text: In order to digitize exact features in arcgis 9.3.1, our draftsman would like to use auxiliary features (lines and points mostly) that don't have to be saved later on. I told them the easiest way was to digitize the auxiliary features into a seperate shapefile - is there a simpler or better way? (in case I didn't explain what I want very good I'll give an example. Suppose someone has to digitize a lighting post (point feature) 4 meters from a curb. In Autocad our draftsman would have digitized a line four meters long, using snapping he'd put a lighting post feature in the end of the line (in a different layer), and in the end he'd erase the line. Is there any way to mimck the process without actually digitizing into a shp/gdb?) Thanks in advance.


title: How to copy the contents of an attribute from one vector layer to another one that touches it or intersects it?

text: I have 2 layers that intersect, called layer A and layer B. Layer A has an attribute (X) that I want to copy wholesale (for the whole layer) to layer B, so that when object A1 intersects object B1, attribute X1 from A1 is copied to B1. I am trying to do this in QGIS. I used to do this in Mapinfo with an UPDATE query... then I switched to Mac. I am new to QGIS.


title: How to perform intersect operations in the Field Calculator in QGIS?

text: How to use the Intersect function? How would I specify to               REPLACE layer1.field1 with layer2.field2 WHEN layer1.geometry INTERSECTS layer2.geometry      Is that feasible in QGIS?


title: Partial Map Pixel Projections

text: I have the following:   1. A PNG or JPG of a random width and height, in pixels, that represents a partial map of the world (i.e. a rectangle of the United States).   2. The four corners of this image are represented by _real_ latitude and longitude values that if plotted, would represent 4 corner points of the United States (LatTop, LatBottom, LongLeft, LongRight).   3. A coordinate that, if plotted, would represent a point somewhere inside of this "rectangle" (Latitude, Longitude). My question is, what type of projection or formula would I use to calculate the X and Y coordinates (in pixels based off of the original PNG/JPG image) of the coordinate provided in #3 above? Some example values would be:   * ImageWidth = 1312   * ImageHeight = 576   * LatTop = 30.244381   * LatBottom = 29.529928   * LongLeft = -96.764958   * LongRight = -94.404619   * Coordinate Latitude = 30.100000   * Coordinate Longitude = 95.100000 I found some information regarding Plate Carrée projections but I could not successfully apply it to this problem. Perhaps I am missing something. Thank You


title: OGR merge shapefiles of different geom type into a single GeoJSON feature collection

text: I'd like to merge shapefiles of different geometry types into a single geojson feature collection, but OGR says it doesn't support update of GeoJSON files:               $ ogr2ogr -f GeoJSON -t_srs crs:84 -lco GEOMETRY_AS_COLLECTION merge.geojson points.shp     $ ogr2ogr -f GeoJSON -t_srs crs:84 -update -append merge.geojson -lco GEOMETRY_AS_COLLECTION lines.shp -nln merged     ERROR 4: GeoJSON Driver doesn't support update.     FAILURE:     Unable to open existing output datasource `merge.geojson'.      This seems like a common use case to me. Am I missing something obvious in the ogr2ogr manual? I can put together something in python (and finally use fiona) but I'm rather surprised I don't see anything existing when searching here or github... which makes me thing I'm messing something up. thanks!


title: Qgis Server - Ubuntu 13.04 - Can't create a directory

text: I'm Portuguese (sorry about my English) and this is my first post here. It would be very useful to me, to discover how to get to work qgis server+ qgis web client. I'm a Windows user, so, I followed some tutorials to achieve that mission. I had luck with qgis server, but was impossible to figure out how to put the qgis web client reading my qgis projects. So, even feeling myself like a fish out of water,I installed Ubuntu 13.04, and first of all, I tried to follow Anita Graser tutorial : http://anitagraser.com/2012/03/30/qgis-server-on-ubuntu-step-by-step/ The problem is, that I'm stuck in this step: "Adding a QGIS project file It’s time to add a QGIS project to our server. To do that, we move to the QGIS Server folder cd /usr/lib/cgi-bin where you should find qgis_mapserv.fcgi and wms_metadata.xml. I’ve decided to have one folder for each project file. My first project is “vienna”. sudo mkdir vienna cd vienna" When I put in ubuntu terminal: "`sudo mkdir vienna`" it appears:               "sudo: unable to open /var/liib/sudo/vitorjcj/1: Read-only file system     mkdir: cannot create directory 'vienna': Read-only file system"      How do I solve this problem? PS: Take into consideration that is my first time playing with linux :) **EDIT:**               vitorjcj@ubuntu:~$ cd /usr/lib/cgi-bin          vitorjcj@ubuntu:/usr/lib/cgi-bin$ dir          admin.sld  qgis_mapserv.fcgi  wms_metadata.xml          vitorjcj@ubuntu:/usr/lib/cgi-bin$ sudo mkdir vienna          [sudo] password for vitorjcj:           sudo: unable to open /var/lib/sudo/vitorjcj/1: Read-only file system          mkdir: cannot create directory ‘vienna’: Read-only file system


title: Problem loading raster QGIS 2.0

text: I am very new at GIS. I am trying to load a raster so that I can complete a tutorial on georeferencing. The raster file is a jpeg. However, whenever I try to load the raster, it gives me solid colored lines that are completely distorted from the original. I realize this is probably a very simple fix. Example: ![example](http://i.stack.imgur.com/p7Fb4.jpg)


title: Routing with osm2po - starting from middle of the street

text: I'm currently looking at using osm2po's java library to build a simple application that builds a shortest path between my favorite bars :) However, I'd like to start at any location, not at an intersection. I saw a response that's it's theoretically possible (http://www.rqna.net/qna/nwkrrs-can-i-start-and-end-an-osm2po-route-in-the- middle-of-a-road-segment.html), but not really sure how to do this with the latest java API. Is this functionality exposed? I've seen findClosestVerexId(lat, lon), but what I really want is something like findEdge(lat, lon), and then a way to get source and target from the edge. Thanks!


title: How to bulk change the schema name of layer properties in ArcMap

text: I would like to change the data source location (schema) in layer property in ArcMap. Our data are stocked in an Oracle Spatial database. Recently, the DBA has changed the schema name and then I cannot access the layers within my MXD. So I think it is possible to bulk change the schema using ArcPy functions. But according to the documentation I don't know how. Until now I've tried this:               mxd = arcpy.mapping.MapDocument("CURRENT")     layers = arcpy.mapping.ListBrokenDataSources(mxd)     for layer in layers:         layer.replaceDataSource(...) // I don't know what to put here`      The schema name was before BICCDDN and now it is EDMCDD. All table names has not changed at all. Only the prefix schema name. Thanks to help me. Here is the place where I can manually change the schema name. ![enter image description here](http://i.stack.imgur.com/oiGvG.png)


title: Geostationary Satellite Projection

text: What is the equation needed for forward and inverse transformation of Geostationary Satellite projection.?. How to convert it to Geodetic and vice versa. Geostationary satellite projection is mentioned in Proj4 as here and in wiki


title: Duplicating Features To Create A Pattern

text: I'm currently trying to take the 'Shifting Features' script from ArcPy Cafe and this SE response and apply it to a similar script that will one by one, through all the features in a feature class, copy the first feature to another feature class, then create a copy of that feature X times with a set interval and angle. So with this I should be able to take, say a feature class containing two lines, and create 5 lines, 3 feet apart at a 45 degree angle in a new feature class. I figured this should be easy by nesting cursors, one for searching the original FC, one for inserting rows in the new FC, but I apparently don't know much about cursors in general, or SearchCursor or InsertCursor in general. The InsertCursor command runs the correct number of times, and creates the correct number of features, but the features have no geometry. Can anyone help me figure out what I might be doing wrong? Thanks!               import arcpy, math          # all of these will eventually be made arguments          fc = r"c:/geometry/polylines.shp" # source feature class     fc2 = r"c:/geometry/polylines_shift.shp" # destination feature class     c_shift=3 # desired overall offset (hypotenuse)     no_features=5 # number of features to duplicate between offsets     angle=45 # angle of offsets          # given the angle and the hypotenuse, determine the x and y coordinate shift     x_shift=c_shift * math.cos(angle)      y_shift=c_shift * math.sin(angle)          # I know this will copy the existing features from fc to fc2     #This will eventually be replaced with a command to create a blank feature class once I get the meat of the script fixed     if arcpy.Exists(fc2):         arcpy.Delete_management(fc2)     arcpy.Copy_management(fc,fc2)          i_cursor = arcpy.da.InsertCursor(fc2,("SHAPE@XY"))     with arcpy.da.SearchCursor(fc, ["SHAPE@XY"]) as cursor:         for row in cursor:             for n_feat in xrange(1, no_features+1):                 i_cursor.insertRow([[row[0][0] + (x_shift*n_feat or 0),                                      row[0][1] + (y_shift*n_feat or 0)]])     del row, cursor, i_cursor


title: Browse (add/remove) ArcGIS Server Services (Layers) from the ESRI Javascript API viewer

text: I am quite new to the Javascript API. I am not if there is a widget that does what I would like but I am trying to create a web map that allows the users to choose/add/remove layers located on our GIS Server. I cant figure out how to achieve something similar as the picture attachment. If anyone have any samples/tips, please let me know. Thank you much for your time. Alex ![What I would like to achieve \(from silverlight API\) ](http://i.stack.imgur.com/tKqGd.png)


title: tilestache wms layer in qgis and open layer rendered layer

text: I have setup tilestache and styled a map for using mapnik xml and I can call on the url in a web browser, how can I call the url for a wms layer in QGIS or is there a way to view mapnik xml in qgis. I have tried different combinations and i have experience with wms from geoserver and mapserver. making a request at localhost:7980/sa/1/1/1.png gives me a png. i have tried to add it to open layers but i am not getting it right which i thing is caused by a bad url I am passing. var tilestache = new OpenLayers.Layer.XYZ(   "Track for Africa",   "http://www.mysite.com:7980/sa/${z}/${x}/${y}.png",


title: Prevent parameter defaults in service definition file

text: I have a Python toolbox with a single tool designed for 10.2. The tool is intended to be a geoprocessing service on a server. So I run the tool, and after it succeeds, I right click the result and choose Share As --> Geoprocessing Service. In the wizard, I choose "Save a service definition file" and move through the wizard until it switches to the dialog where I can set service properties. I Analyze, and there are no errors or warnings. At this point, I notice that my parameters are being given default values. In particular, they get the values that I used when I ran the tool to create the results file. ArcGIS's dialog will not let me modify the default parameters; the text boxes are read only. I want to prevent this. Ideally, my parameters would have no default values, but empty string default values would also be acceptable. All my input parameters happen to be strings. How can I prevent ArcGIS from using the result's parameters as the defaults?


title: How to get pixel size/cell size of raster layer?

text: After much searching, I've still not been able to figure out how to perform what seems like a simple function: I'd like to get the pixel size/cell size of a raster layer using PyQgis. Is there any native way to do this through QgsRasterLayer or a helper class? Or will I need to read the file in using GDAL functions, e.g.:               geotransform = dataset.GetGeoTransform()     if not geotransform is None:         print 'Pixel Size = (',geotransform[1], ',',geotransform[5],')'      Thanks for your help. Alex


title: How to batch import GDAL supported raster files into GRASS?

text: I am trying to import some GDAL supported raster files into a binary raster map layer using the r.in.gdal module from QGis Grass tool. As I have multiple files to import, I would like to batch import them. Is is possible or should do I have to repeat the same operation for every single file? I work with QGis 1.8.0 Lisboa.


title: How to buffer in pyQGIS?

text: I'm looking for some examples to do simple python scripts in QGIS. How would I do a buffer analysis on a dataset? I can't seem to find much in the manual and QGIS python that quite matches the Esri documentation.


title: Masked data included in ISODATA classification computations?

text: I am using ENVI to mask NDVI MOD13Q1 images using the "pixel_reliability" raster included in the MOD13Q1 hdf. All instances of cloud cover (pixel_reliability=3) have been masked (=0). I have applied these masks to 23 images and stacked them into a layer (there may have been an easier way to batch apply the pixel_reliability mask to the relevant NDVI raster, feel free to share!). When I look at the spectral signature of the NDVI stack in which each time point has been masked with the respective pixel_reliability raster, the masked values show up as zero, distorting the shape of the spectral curve. Does ENVI include the zero values in its ISODATA classification algorithm? Is there a way to have ENVI recognize 0 as "missing data" and ignore the value it the classification process (I have already changed the "No Data" value using the Edit ENVI Header tool and this doesn't make a difference in the spectral signatures)? Is there a better way to use one data cube (pixel_reliability) to mask another (NDVI) at each time point? Many thanks, Emily


title: I've installed OpenGeo Suite, how do I use it with Leaflet?

text: I've just finished installing OpenGeo Suite on my VM Ubuntu web server. I'd like to be able to use Leaflet from leafletjs.com to view my maps with on my mobile devices. I'm supposed to create an HTML page, I assume, and use different snippets of code to make my map with it, however I don't know how to link the OpenGeo and Leaflet together. It mentions using a tile layer URL, however I cannot figure out how to create this URL or where to find it. Also, do I have to create a new Leaflet HTML page for every map I want or am I way off? I'm sure I am missing some key elements of exactly how OpenGeo works, but I'm having a hard time finding any beginners guides and the documentation isn't very beginner friendly in some areas.


title: Are there tools to directly import QGIS projects into Tilemill?

text: Now that there is a tool to convert an MXD to Tilemill, one must then wonder if the ability to open QGIS documents in Tilemill will ever exist.


title: Qgis Server and WCS service : problem to display raster layers

text: I'm using Qgis Server and I am experiencing a few problem. I would like to display some raster layers ( for basemap) with other vector layers but I only obtain a white backgroung. **My situation :** I have one project ( qjs file) containing all the layer i want to retrieve with this url : localhost/qgis/qgis_mapserv.fcgi.exe?SERVICE=WMS&VERSION=1.3.0&SRS=EPSG:2154&REQUEST=GetMap&map=c:/myprojects/qgis/V3_dev.qgs&BBOX=320000,6300000,550000,6500000&WIDTH=1000&HEIGHT=1000&LAYERS=ORTHOPHOTO&FORMAT=image/png ORTHOPHOTO is the group containing some raster layers : http://nsa34.casimages.com/img/2014/06/12/140612095507246249.png But the only thing i got is a perfect white picture, no errors. There is no problem with Bbox or EPSG because i can easily display vector layer. I tested to directly select ortho1 (included in ORTHOPHOTO) but same result. I'm quite puzzled because i read that since the 2.2 version of qgis, raster layer could be use in qgis server. Source (chapter Qgis Server) : http://docs.qgis.org/2.2/fr/docs/user_manual/preamble/whats_new.html **My configuration :** → Windows 7 x64 → Packages installed with OsGeo4w following this tuto : hub.qgis.org/projects/quantum-gis/wiki/QGIS_Server_Tutorial (qgis server / gdal / fcgid) → Apache 2.4.9 with ''fcgid_module'' (and ''cgi_module'' too) → The V3_dev OWS properties : nsa34.casimages.com/img/2014/06/12/140612102742321017.png **What i tried :** → Some re-installation → With the GetCapabilities request i can see the group ORTHOPHOTO queryable. → Qgis server on a ubuntu VM but same result. → i tried to convert raster layer into vector layer : layer quality fail → I read that someone have the same problem on CentOS. Source : www.forumsig.org/archive/index.php/t-38101.html → I recently discovered the GetCoverage request : localhost/qgis/qgis_mapserv.fcgi.exe?SERVICE=WCS&VERSION=1.1.0&CRS=EPSG:2154&REQUEST=GetCoverage&map=c:/myprojects/qgis/V3_dev.qgs&BBOX=320000,6300000,550000,6500000&WIDTH=1000&HEIGHT=1000&COVERAGE=ORTHOPHOTO&RESPONSE_CRS=EPSG:2154&FORMAT=image/png and i obtain : The layer for the COVERAGE 'ORTHOPHOTO' is not found → i'm also wondering : Is it possible that some librairie (like gdal) is missing in the OsGeo4W installer ? Hoping someone can help me :) PS : Sorry for the french links


title: GIS Survey Line Plans

text: Has anyone ever heard of any tools in either Arcmap, QGIS, Mapwindow GIS etc, that could create a line plan (Survey Lines) across a polygon (Survey Area)? A form of input such as distance between lines and minimum length would hopefully be required. If this tool could be tailored to create the shortest, so most efficient line plan, then that would be even better. Obviously a form of automation rather than just manually offsetting poly-lines is required- as that can get very monotonous over large areas! Any help is greatly appreciated!! ![enter image description here](http://i.stack.imgur.com/dYE5r.jpg) The image attached shows the kind of output I need, however the Survey Areas are usually much more irregular!!


title: Disable customized context menu for some cases using python

text: I create a python add-in tool to add new item to the context menu, but I want the new item enabled only when user right click on map when within a buffer of polyline feature. I haven't figure out how to make it, which means for the area qualified (within a buffer of polyline) can see the new item on context menu, others won't show the new item (or the new item is disabled in grey). I am using ArcGIS 10.2.2 and python.


title: Geometry to Shapefile

text: From importing the shapefile to postgis, its geometry is added to the database postgresql. The question, how can I load the shapefile in the Visual Basic.Net using the geometry in the database ?


title: How to translate tens of images from a directory as a batch with gdal_translate?

text: For example i have a dataset of images in *.cub format like               $path: /work/120614/mg_1164/     f01.img.cub_E.cub     f02.img.cub_E.cub     '     '     f0n.img.cub_E.cub      for converting a single image I use `gdal_translate f01.img.cub_E.cub fo1.tif` for whole the datasets.HOW? batch file?


title: List of Offensive Symbols/ Colours for Cartography

text: I was at a data visualization workshop a few days ago and someone there asked if there are any symbols/colours that might offend people when making charts, the presenter didn't really understand the answer so I helped to explain how in cartography using certain colours to represent certain groups could offend. For example using yellow to represent East Asian's, Red to represent Native Americans, or Red to represent basically any minority (since red has a connotation of being negative/bad/dangerous) The only symbols I can think of that have negative connotations are yellow six pointed stars to identify Jewish people, and pink triangles to identify gay people as these were used in Nazi Germany. I realized I only ever learned this from working with other researchers and was wondering if there is some comprehensive list out there for these types of colours/symbols.


title: Tilemill use custom font

text: I want to use the font `Air Soft Condensed heavy` in Tilemill. I placed the font file `Positype-AirSoftCond-Heavy.otf` into the font directory of my project:               /home/<user>/Documents/MapBox/<project>/fonts/      The `palette.mms` file has the line                Map { font-directory: url(./fonts); }      How can I reference now my font to use it?


title: Can leaflet read OSM xml format?

text: I want something similar to `OpenLayers.Format.OSM` for leaflet. i know how to use a geojson file and i have found support for GPX and WFS in plugins, but not for OSM. Simply saying i want to make this map using leaflet and overpass.


title: Why does Fetch Python Plugins responded with only the installed ones?

text: I use QGIS 1.8. Whenever I demand to install plugins through Fetch Python Plugins it displays only those the already installed. I couldn't find for example; point sampling, group stats and photo2shape plugins. Can one help me?


title: Split Line at intersection and attach attributes

text: I want to slit line at intersection and attach attributes and export I have roads table with columns id,road_name,road_type,geom with 2115 records. I try using this query - `select st_astext((st_dump(st_union(geom))).geom) from roads;` it successfully splits the line gives 5114 records but when i am trying to attach attributes it fails `select id,road_name,road_type,st_astext((st_dump(st_union(geom))).geom) from roads;` ![split at intersection](http://i.stack.imgur.com/XeERQ.png) I need to split line at intersection also keep attributes means line 1 divide in 3 having same id is 1 , road_name is abc similar to line 2 divide in 2 having same id is 2 , road_name is xyz and so on.


title: Open Source Challenge: Automated Mastermap update in Oracle Spatial

text: I've been set a challenge which is a bit out of my comfort zone. I'm an esri developer trying to break into the open source world. I've been asked to automate the update of some oracle SDO datasets. I have a boundaries dataset (X) that needs to match topologically with OS Mastermap data. At the moment, the Mastermap update is a manual process (disk in the post!). If a boundary in the OS data is updated e.g a field is reduced as a result of a new road, the same field in dataset X needs to be updated accordingly. At first, I'm expected to flag up any changes that will need to be made to dataset X. Eventually, dataset X will be updated automatically. I also need to store the changes made for audit purposes. I've been challenged to achive the above with an open source stack. My initial thoughts are that a workflow could look like – download update -> import dataset -> select all new and updated toids in OS dataset -> some kind of spatial query to identify effected features in dataset x-> export selected features to archive -> edit selected to snap to new OS boundaries. Has anyone attempted something similar? What tooling would the community recommend? All suggestions welcome, James


title: Optimal video cards for GIS programs

text: I recently purchased a Dell XPS-8300 (i7-2600|16 GB DDR3 RAM|64-bit Windows7) which I intend to use primarily for GIS programs. I bought a refurbished model because it was much cheaper and because I am able to choose a graphics card of my own. I am currently in the process of purchasing said card, but I haven't really seen much information about what kinds of cards are optimal for ArcMap or QGIS rendering as opposed to Battlefield 3 or Skyrim. I am currently looking at an XFX AMD Radeon HD6870 PCIE 2GB Dual Mini or an XFX AMD Radeon HD 6870 900M 1 GB DDR5 DUAL MINIDP HDMI DUAL DVI PCI-E Video Card The 2 GB model is about $40 more than the 1 GB. I have a 460 watt power supply and do not know if the extra ram affects power consumption. The 1 GB model, from all reviews, should work fine with my system, so I am hoping the 2 GB would work fine as well. However, since I already have 16 GB of onboard ram and have read that GPU RAM isn't really important unless you go under the amount needed for a given program and settings, I am wondering if this extra GB would be worth it for my purposes: For example if I have a parcel dataset with hundreds of thousands of polygons, would that extra GB of GDDR5 RAM significantly improve the time it takes to render the polygons when I scroll/zoom in on ARCmap/QGIS; or does it mostly worthless due to the fact that I already have 16 GB that can be used? Also are there any specific advantages/disadvantages to using a ATI vs NVIDIA card for GIS display programs in general? Thanks in advance!


title: Shapefiles on WMS Server

text: I've been reading about using WMS servers to server shapefiles. When a server does support shapefiles over WMS, what format does the client download the shapes in? Will it typically be an image format (png, jpeg, etc) like it is for CADRG, or is there a way to extract the vector data, or will the shapefile data (geometries, etc) typically be part of the metadata available on the server?


title: Where to start with Web GIS development, or which technology to select?

text: I am very new to GIS development, and to be be frank I have no background about it at all. I searched the web but the tutorials I found seemed to assume the reader has some background information. the thing is that I am confused about what to read or learn, there seems to be lots of technologies, and I feel lost since some speak about openlayers, geoserver, mapserver, google maps, and open street maps. So here is what I am supposed to develop, and I hove you could give me an advice about which technology to use, and where should I start reading - given that I know almost nothing -. Case 1: a closed system for about 20 users only, who can specify locations on the map, and the web application will store the latitude and longitude of the locations and show the markers. I wanted to use google maps api, but I cancelled that since there license requires you to purchase the service if the system is a closed one. so what technology should I use in such case? I need a free option, also I will be only using web server, so if the solution includes using my own geoserver, or something like that I won't be able to do it. Case 2: I am supposed to display the roads and routes between two given points, and probably add some notes on the map. For this I case I can use my own map server/geo server, but again I want your suggestions. of course the solution need to be open source finally, I hope you could tell me what to start reading first, thanks for you in advance.


title: Geocoding or adding region to imported Access table in MapInfo

text: I have imported an Access database of 500+ sites into MapInfo. I've managed to geocode all but 80 of them against other tables that have the sites as regions with the same ID number, and replaced the geocoded objects with regions. What can I do about the 80+ sites that have not been geocoded and need a region manually assigning to them? Selecting and adding a region just creates a new record. There is no object to do 'clear map objects' as there would be if the record was geocoded.


title: Can' t edit anything in Arcgis 9.3

text: My name is Ramona and I'm from Romania. I'm new to GIS. I want to digitize a topographic map in ArcGIS 9.3 and cannot. I tried these steps: I created a Personal geodatabase inside a raster dataset, then I loaded a tiff (mosaic) and a Personal geodatabase feature dataset (with more Personal geodatabase feature classes). However, when I start to edit, the software puts the last vertex. I'm using Windows 7 Professional.


title: Qgis 1.9 on Ubuntu 12.04 constant crashes as soon as I try to do anything

text: I can get a google earth physical layer open but as soon as I try to open a vector file it crashes-just goes away-no error message, nothing


title: OpenLayers - Test Control To See If It Is Active

text: I am trying to get the active state of an OpenLayers SelectFeature control. I create this control and add it to my map but I do NOT activate it yet. This part works fine. Here is the script to create the control and add it to the map:               var highlightControl = new OpenLayers.Control.SelectFeature(lyrstands, {             hover: true,                     highlightOnly: true,                     selectStyle: {        fill: "#0000FF",                strokeColor: "#0000FF",                strokeWidth: 2,        fillOpacity: 0.5               }         });               map.addControl(highlightControl);      Then I have a button that calls a function. Here is the code for the button:               new OpenLayers.Control.Button({title: "Select Polygon", trigger: HiLite})         I know the button is calling the HiLite function because if I put a simple alert statement in the HiLite function (for example: alert("The Function Works");) it works correctly. So what I am trying to do is check to see if highlightControl has been activated by testing the highlightControl.active value, if it is already activated the function should do nothing, but, if highlightControl is NOT activated it should activate it. Unfortunately, the function does not work. Here is the function I have written:                  function HiLite(){           if (highlightControl.active==false)                          {highlightControl.activate();}           else {}            }          If I write the function to simply activate the tool without the If Else stuff it works fine. But, for some reason when written as above I cannot get the status of the Active attribute. The OpenLayers documentations states that the Active attribute is a boolean value and I assumed I could check the status using the highlightControl.active==false statement. I would appreciate any insight on why I cannot get this function to work. Thank you!!


title: Define area of study in Google Earth before Import in Qgis

text: First, I would like to thank you, for your answers. **Here is my problem:** I wish implement a GIS solution to manage a water and wastewater network. In my project I have any map for use it like background. So, I think the only solution for having a layer background, is to define it from Google Earth. I must define de scope of study or area of study in Google Earth, after I import it in Qgis. Can someone help me in this step, and give me advice ? Also, if someone has any documentation or tutorial, it can help me. Thanks and regards


title: QGis CadTools CAD Console

text: I need help with the Qgis CadTools CAD Console commands. According to these help files there are six commands. I just can't get them to work. Can someone help me? Thanks


title: How to use Sextante with PostGIS

text: Is there any way to use analysis from Sextante with PostGIS Geometry Type? Edit: I am trying to create a webapp that make some analysis (like interpolation for example) and returns the results. So I don't want to use a desktop solution. Python could be a solution... But how?


title: Altitude Above channel in GRASS 6.4

text: i'm fairly new to the GIS world and GRASS (come from econometric background), i would like to do some species distribution modeling using R and GRASS. So i started to read the book of Tomislav Hengl : **A practical guide to geostatistical mapping** , and specially **chapter 8** about species modeling (using the bei dataset in the spatstat R package). The same analysis is also on the author wiki page : http://spatial- analyst.net/wiki/index.php?title=Species_Distribution_Modelling Tomislav Hengl, used mainly SAGA-GIS and R, and i would like to do the same analysis using R and GRASS (through spgrass6). The problem is that i don't know the equivalent of SAGA ta_channel module 3 (altitude above channel) in GRASS, so far i tried r.stream.distance to compute elevation above stream but the result are really different. This is the part of original code i want to do (problem is wherer there are **):               library(gstat)     library(spatstat)     library(splancs)     library(rgeos)     library(maptools)     library(adehabitat)     library(rgdal)     library(lattice)     library(RSAGA)     # General settings for plots:     trellis.par.get("fontsize") -> fontsize     fontsize$default<-34     fontsize$points<-22     fontsize$text<-34     # environnement pour SAGA     monenv <- rsaga.env(path = "/usr/bin/", modules = "/usr/lib/saga")     #-------------------------------     # 1. Data import and preparation     #-------------------------------          # Load the data:     data(bei)     str(bei)     str(bei.extra, max.level=1)     plot(bei.extra$elev, main="Beilschmiedia (Tropical rain forest trees)")     plot(bei, add=TRUE, pch=16, cex=0.3)     class(bei.extra[[1]])     # Export to a GIS (for visual exploration):     grids <- as(bei.extra[[1]], "SpatialGridDataFrame")     names(grids)[1] <- "elev"     grids$grad <- as(bei.extra[[2]], "SpatialGridDataFrame")$v     writeGDAL(grids[1], "dem.mpr", "ILWIS")     write.asciigrid(grids["elev"], "dem.asc")     write.asciigrid(grids["grad"], "grad.asc")     bei.pnt <- data.frame(x=bei$x, y=bei$y, no=rep(1, length(bei$x)))     coordinates(bei.pnt) <- ~x+y     writeOGR(bei.pnt, "bei.shp", "bei", "ESRI Shapefile")     # reformat the point pattern so it fits the grids:     bei.ppp <- ppp(coordinates(bei.pnt)[,1], coordinates(bei.pnt)[,2], marks=bei.pnt$no, window=as(grids[1], "owin"))     # Generate the wetness index and vertical distance from the channel network in SAGA and import back to R:     rsaga.esri.to.sgrd(in.grids="dem.asc", out.sgrds="dem.sgrd", in.path=getwd())     rsaga.geoprocessor(lib="ta_hydrology", module=15, param=list(DEM="dem.sgrd", C="catharea.sgrd", GN="catchslope.sgrd", CS="modcatharea.sgrd", SB="twi.sgrd", T=10))     rsaga.geoprocessor(lib="ta_channels", module=0, param=list(ELEVATION="dem.sgrd", CHNLNTWRK="chnlntwrk.sgrd", CHNLROUTE="channel_route.sgrd", SHAPES="channels.shp", INIT_GRID="dem.sgrd", DIV_CELLS=5, MINLEN=10))     **rsaga.geoprocessor(lib="ta_channels", module=3, param=list(ELEVATION="dem.sgrd", CHANNELS="chnlntwrk.sgrd", ALTITUDE="achan.sgrd", THRESHOLD=0.1, NOUNDERGROUND=TRUE))     rsaga.sgrd.to.esri(in.sgrds=c("twi.sgrd","achan.sgrd"), out.grids=c("twi.asc","achan.asc"), prec=1, out.path=getwd())**     grids$achan <- readGDAL("achan.asc")$band1     grids$twi <- readGDAL("twi.asc")$band1      Here is my code (the problem is where i put **):               library(gstat)     library(spatstat)     library(splancs)     library(rgeos)     library(maptools)     library(adehabitat)     library(rgdal)     library(spgrass6)     library(ggExtra)     # environnement pour GRASS     loc <- initGRASS("/usr/lib/grass64", home=tempdir(), override = TRUE)     #-------------------------------     # 1. Data import and preparation     #-------------------------------     # Load the data:          data(bei)     bei.data <- as.data.frame(bei)     bei.elev.data <- as.data.frame(bei.extra$elev)     ggplot(bei.elev.data, aes(x, y, z = value)) + geom_tile(aes(fill = value)) + geom_point(aes(x, y, z = NULL), data = bei.data, size = 2, pch = "+") + scale_fill_gradientn(colour = terrain.colors(4)) + opts(title = "Beilschmiedia (Tropical rain forest trees)")     class(bei.extra[[1]])     # Export to a GIS (for visual exploration):     grids <- as(bei.extra[[1]], "SpatialGridDataFrame")     names(grids)[1] <- "elev"     grids$grad <- as(bei.extra[[2]], "SpatialGridDataFrame")$v     writeRAST6(grids["elev"], "dem", zcol = "elev",  ignore.stderr = TRUE)     writeRAST6(grids["grad"], "grad", zcol = "grad", ignore.stderr = TRUE)     bei.pnt <- data.frame(x=bei$x, y=bei$y, no=rep(1, length(bei$x)))     coordinates(bei.pnt) <- ~x+y     # reformat the point pattern so it fits the grids:     bei.ppp <- ppp(coordinates(bei.pnt)[,1], coordinates(bei.pnt)[,2], marks=bei.pnt$no, window=as(grids[1], "owin"))     # Generate the wetness index and vertical distance from the channel network in GRASS and import back to R:     # Work on  GRASS:     execGRASS("g.region", parameters=list(rast="dem"))     execGRASS("g.region", parameters=list(rast="grad"))     execGRASS("r.info", parameters = list(map = "dem"))     execGRASS("r.info", parameters = list(map = "grad"))     gmeta6()     # wetness index     execGRASS("r.topidx", flags = "overwrite", parameters = list(input = "dem", output = "twi"))          # r.stream.extract     execGRASS("g.list", parameters =  list(type = "rast"))     execGRASS("r.stream.extract", flags =  "overwrite", parameters = list(elevation = "dem", stream_rast = "stream", direction = "stream_dir", threshold = as.integer(10)))     #execGRASS("r.watershed", flags = c("m", "overwrite"), parameters = list(elevation = "dem", stream = "stream", threshold = as.integer(10)))     # altitude above channel     **execGRASS("r.stream.distance", flags = "overwrite", parameters = list(stream = "stream", dir = "stream_dir",  dem = "dem", method = "downstream", elevation = "achan"))**     # spline to fill      execGRASS("r.fillnulls", flags = "overwrite" ,parameters = list(input = "achan", output = "achan"))     # exportation     grids$achan <- readRAST6("achan")@data$achan     grids$twi <- readRAST6("twi")@data$twi      So my question is what is really **Saga altitude above channel** (if possible in plain english) and how to compute this using GRASS ? I'm made several researchs (for weeks now), and may be it's right place here. Thank you in advance and sorry for my poor english


title: set infoTemplate for various graphics in graphics layer in loop

text: I am successful at getting an info window pop up on the event click of a marker on a map. Now, what I'm having trouble with is passing unique information for each graphic in the loop to the individual info window. Here's my loop               for(i=0;i<coords.length;i++){                 xcoord=coords[i][1];                 ycoord=coords[i][2];                 type=SimpleMarkerSymbol.STYLE_SQUARE;                 if(coords[i][3]=='stream'){                     var point = new Point(xcoord, ycoord);                 }                 var graphic = new Graphic(point, new SimpleMarkerSymbol(type).setSize(10).setColor(usgscolor));                 //these next four lines are where I'm having confusion                 var infoTemplate = new InfoTemplate();                 infoTemplate.setTitle("title");                 infoTemplate.setContent("content");                 graphic.setInfoTemplate(infoTemplate.setTitle("title"));                 glayer.add(graphic);             }      Then I add each graphics layer at the end and an event listener for the pop up               map.addLayer(glayer);//add all the individual grahics in the graphics layer to the map                        glayer.on('click', function(evt) {             map.infoWindow.setTitle(graphic.getTitle(infoTemplate));             map.infoWindow.setContent(graphic.getContent(infoTemplate));             map.infoWindow.show(evt.screenPoint,map.getInfoWindowAnchor(evt.screenPoint));         });      These various topics on the arcgis developer site is what I'm using for reference but I still can't seem to put the pieces together https://developers.arcgis.com/javascript/jshelp/intro_formatinfowindow.html https://developers.arcgis.com/javascript/jsapi/graphic-amd.html https://developers.arcgis.com/javascript/jsapi/infotemplate-amd.html https://developers.arcgis.com/javascript/jsapi/infowindow-amd.html


title: how can i symbolyze a layer with graduated colors

text: can any one tell me how to use isimplerender interface to symbolize a polygon layer with graduated colors according to a value of a field such as "population" thanx


title: Load a .sav file in ENVI

text: I have a 23-band .sav file that I need to load in ENVI. The 23 bands are 23 time points of MODIS NDVI data. I ultimately want to run a classification on this data. I have been using Data Manager to load the .sav file, but have only been able to load three bands at a time (RGB). How can I load the ENTIRE 23-band file at once?


title: Which projection is being used in this figure?

text: I would like to use the projection being used in this figure. Does anyone know which one it is? Esp helpful if it is available in ArcGIS.   From plant stress page ![http://www.nasa.gov/images/content/711179main1_drought_ESI-2011-226x.JPG](http://i.stack.imgur.com/wPrzW.jpg) thanks!


title: ModelBuilder error on validating when deleting not yet created gdb

text: I have a parent model with a submodel nested in it and tools outside the submodel but within the parent model that connect to the submodel. In the submodel, the first tool is "create geodatabase". In the parent model, the last tool deletes the gdb. When the model is run (not edit), it looks like the model validates before running and comes up with the error that it cannot find the gdb to delete at the very end. Obviously this is because the gdb has not yet been created at the beginning of the model. Is there any way to get around this where it would not validate before running?


title: How can I insert new records to an ArcSDE feature class from a view in ArcSDE

text: I am wondering how I can both insert new records from a view to a feature class and update data that may be newer in the view than in the feature class. I have been fumbling around in modelbuilder and arcpy and haven't made any significant headway yet. I need to be able to export the model or script so that it can be run as a recurring job. Thanks!


title: How to do selection with a WFS source in OpenLayers?

text: One of my students wants to style a WFS layer and have selections too. I got as far as this code http://ian01.geog.psu.edu/geoserver/www/wfs/wfs-sld2.html which styles the WFS (though with out the text labels) and while the selection control seems to be working (i.e. the box is blue not red) nothing is actually selected.


title: Is there a Postgis SQL Coding Style Guidelines?

text: Most programming languages have Coding styles. Some have defined standards, other have just some code of conduct that the majority of programmers follow. Standard coding styles are have obvious benefits, in gis.stackexchange itself, it would help people better understand both questions and answers. From some "googling", there are some reference to Coding Styles for SQL, although they are not consensual:   * http://stackoverflow.com/questions/118288/sql-coding-style-guide   * http://stackoverflow.com/questions/522356/what-sql-coding-standard-do-you-follow   * http://www.simple-talk.com/sql/t-sql-programming/sql-code-layout-and-beautification/ While this SQL coding styles can be used for Postgis, my question is if there is a Postgis\SQL coding style to handle its functions, or at least some common practices.


title: How to draw shapes with highest possible decimal accuracy

text: I am trying to read input of lat and long to draw the shape based on condition.The input is standard array of lat and long with 15 decimal accuracy. I am realizing - I am loosing some decimal accuracy (the point was somewhat offsite and complex polygons were never drawn). The samples are here: point: [[[-81.00392112682097,46.50255058560673]]] the plotted point is: -81.003906 46.502686 and I am generating PDF maps at 1:3000 scale so it is a big deal. polygon: [[[-79.37485174735863,43.6430783790403],[-79.37361793121188,43.64341222822559],[-79.37323169311376,43.64261253986735],[-79.374186559523,43.64234079868168],[-79.37412218650664,43.64220880966229],[-79.37454061111293,43.64210787668709],[-79.37471227248987,43.64248831547784],[-79.37463717063745,43.64251160757044],[-79.37481956085045,43.64285322389105],[-79.37470154365381,43.642884279823875],[-79.37485174735863,43.6430783790403]]] I used eval to read input as list, but probably that was not a good idea. Not sure how and where should I cast the input as decimal?               OrderCoord =  eval(arcpy.GetParameterAsText(0))     point = arcpy.Point()     array = arcpy.Array()      # A list that will hold each of the Polygon objects      featureList = []     for feature in OrderCoord:         # For each coordinate pair, set the x,y properties and add to the          #  Array object.         #         for coordPair in feature:             point.X = coordPair[0]             point.Y = coordPair[1]             array.add(point)              if OrderType.lower()== 'point':         feaERIS = arcpy.Multipoint(array)         elif OrderType.lower() =='polyline':         feaERIS  = arcpy.Polyline(array)             else :         feaERIS = arcpy.Polygon(array)         # Clear the array for future use         array.removeAll()              # Append to the list of Polygon objects         featureList.append(feaERIS)      I added the screenshot of these points at 2000 scales. We have a Flex application where customer can draw and order a geometry. All the data sits in ORACLE and I am working on a GP Tool to print final PDF. ![I added the screenshot of these points at 2000 scales. We have a Flex application where customer can draw and order a geometry. All the data sits in ORACLE and I am working on a GP Tool to print final PDF.](http://i.stack.imgur.com/LVZzw.png)


title: Problem using ArcGIS API for Silverlight

text: I've install Visual Studio 2010 SP1, Microsoft Expression Blend 4, and ArcGIS API 3.0 for Silverlight. When I create a new project using a standard map application in Visual Studio, I see the error in design tab. When I create a project in Microsoft blend 4, I don't see any template and the appliction shows an error. How can I fix the problem?![visual studio error](http://i.stack.imgur.com/JR9Hx.png) ![microsoft expression blend error](http://i.stack.imgur.com/SeOQv.png)


title: Creating Garmin GPS compatible map using TFW+tiff files

text: I have complete set of 1:50000 maps of Sri Lanka in .tiff format. They have .TFW files accompanied with them. Now I need to create maps to use with my etrex 20 device using those tiff files. How can I do it? Thank You.


title: pgrouting 2 pgr_createTopology not populating source and target field

text: I trying to integrate pgrouting 2.x. I am now using pgrouting 1.x. In version 1 I used `assign_vertex_id`. It already dropped in version 2 and became `pgr_createTopology` function               SELECT pgr_createTopology('roads_r', 0.00001, 'the_geom', 'gid')      It returning query result OK but its not populating source and target field of my road data. Is there any bug in this function?


title: How to set identify tool to identify from visible layers as a default when starting ArcGIS

text: How to set identify tool to identify from visible layers as default when starting ArcGIS without the need to change it from top-most layer to visible layer each time stating arcmap


title: Hovering of a shape with OpenLayers doesn't un-highlight the shape

text: So I've got a list of shapes (zipcodes) which I'm trying to display on an OpenLayers Google map. They render fine. When I hover over them, the shape colour changes. Great! when i -leave- the shape, the colour doesn't revert back the original shape colour. Can someone help me out? Here's my code.. ### 1. Render each shape to the map.               // NOTE: shapeData is an array of objects.     function renderShapes(shapeData) {         for (var i = 0; i <= shapeData.length; i++) {             addShape(shapeData[i]);         }     }      ### 2. Add the shape.               function addShape(data) {         if (!data) return;              if (shapeLayer == null) {             setupShapeLayer();         }              // Reference: http://dev.openlayers.org/docs/files/OpenLayers/Format-js.html#OpenLayers.Format.externalProjection         // This reprojects the WKT from GPS to Google Maps format. Yes, external == the source.         var format = new OpenLayers.Format.WKT({             'externalProjection': new OpenLayers.Projection('EPSG:4326'),             'internalProjection': map.baseLayer.projection         });              var feature = format.read(data.WellKnownText);         shapeLayer.addFeatures(feature);              var highlightControl = new OpenLayers.Control.SelectFeature(shapeLayer, {             hover: true,             highlightOnly: true,             renderIntent: "temporary"         });              map.addControl(highlightControl);         highlightControl.activate();     }      ### 3. Create the layer where all the shapes will exist, in.               function setupShapeLayer() {         shapeLayer = new OpenLayers.Layer.Vector(shapeLayerName, {             styleMap: new OpenLayers.StyleMap({                 "deafult" : new OpenLayers.Style({                     strokeColor: '#999999',                     strokeOpacity: 1,                     strokeWidth: 1,                     fillColor: '#46461F',                     fillOpacity: 0.10                 }),                 "temporary" : new OpenLayers.Style({                     strokeColor: '#2d9908',                     strokeOpacity: 1,                     strokeWidth: 1,                     fillColor: '#3eb216',                     fillOpacity: 0.60                 })             }),             rendererOptions: { zIndexing: true }         });              map.addLayer(shapeLayer);         shapeLayer.setZIndex(1);     }      As I said, the highlight works .. it goes from yellowish to green. but it doesn't revert back to the original default style. :(


title: Grass r.mapcalc generates 1 cell map

text: I'm a GIS student and in my project I need to generate vegetation indexes. I've generated it succesfully yesterday using both python2.7 and grass gui .I'm using WinGRASS-6.5.svn-r54428-453 on Win7. Today -I don't know why- grass started to generate 1 pixel NDVI map. The formula is as simple as it could be: float(b4-b3)/(b4+b3) as most of you already know. My data is not corrupted, I'd saved original maps somewhere else and it was my first check that whether there is a corruption. Also using QGis 1.8.1 I can generate correct NDVI map with same data. You can check it here:               r.info map=ndvi4@PERMANENT                                                            +----------------------------------------------------------------------------+      | Layer:    ndvi4@PERMANENT                Date: Tue Jan 01 04:06:35 2013    |      | Mapset:   PERMANENT                      Login of Creator: murat           |      | Location: newLocation                                                      |      | DataBase: D:\gis\grassWS                                                   |      | Title:     ( ndvi4 )                                                       |      | Timestamp: none                                                            |      |----------------------------------------------------------------------------|      |                                                                            |      |   Type of Map:  raster               Number of Categories: 255             |      |   Data Type:    FCELL                                                      |      |   Rows:         1                                                          |      |   Columns:      1                                                          |      |   Total Cells:  1                                                          |      |        Projection: UTM (zone 36)                                           |      |            N:          1    S:          0   Res:     1                     |      |            E:          1    W:          0   Res:     1                     |      |   Range of data:    min = nan  max = nan                                   |      |                                                                            |      |   Data Description:                                                        |      |    generated by r.mapcalc                                                  |      |                                                                            |      |   Comments:                                                                |      |    float(L71175034_03420050627_B40@PERMANENT -                             |      |    L71175034_03420050627_B30@PERMANENT) /                                  |      |    (L71175034_03420050627_B40@PERMANENT +                                  |      |    L71175034_03420050627_B30@PERMANENT)                                    |      |                                                                            |      +----------------------------------------------------------------------------+     (Tue Jan 01 04:06:51 2013) Command finished (0 sec)                                   I uninstalled WinGRASS-6.5.svn-r54428-453 and installed WinGRASS-7.0.svn-r54453-453, I've deleted mapset and location and the recrated them also. I couldn't think anything else please give me some information.


title: How to install the latest build of QGIS?

text: On Dec. 18, Radim Blazek developed the 'Improved raster singleband pseudocolor classification GUI'. Being new at QGIS, I have no idea how to go about updating to this version. Could someone possibly assist me with this? Thanks in advance.


title: Why is the curl syntax unable to create datastore?

text: I am implementing the following code to upload an existing shapefile as per the documentation,               curl -u admi PUT 'Content-type: application/zip' \           --data-binary @/var/www/geo/shapefile/zip/Poa_horridula.zip \           http://localhost:8080/geoserver/rest/workspaces/acme/datastores/Poa_horridula/Poa_horridula.zip      But I am getting following error, > curl: (6) Couldn't resolve host 'Content-type' > >> No such datastore: acme,Poa_horridula The document says that the curl code will create the datastore then why am I getting this error? **UPDATE: I changed the code to upload a zip file but I am still not able to upload any shapefile.** Where am I going wrong? **UPDATE 2: Ok I am almost there** So I changed the command to,               curl -u admin:geoserver -v -XPUT -H 'Content-type: application/zip' --data-binary @/var/www/geo/shapefile/zip/Poa_horridula.zip http://localhost:8080/geoserver/rest/workspaces/acme/datastores/Poa/file.shp      The only problem is that I am able to add the shapefile but the layer is not being published. Is that normal? Do I have to manually publish it?


title: How to add points in degrees + decimal minutes to QGIS?

text: I try to add geocaches with x,y coordinates in a xls-table to qgis and convert to shp.file. Usally geocach-coordinates are like this: N 51° 29.440' E 011° 56.430' (WGS84) (in german its called "bogenminute) . But it seems that qgis only know geographic coordinate system (WGS 84) with coordinates like this one: 51.754683 12.000967 ("dezimalgrad"). So i can charge them from one to another but i think its sensless. There must be a easier way. Is there another geographic coordinate system ,qgis has, to add the points directly from the first format (bogenminute)?


title: How can I add an index to a CartoDB table?

text: First of all, do I need to do this or does CartoDB do this based on the query on the fly? NYC's 311 database is pretty big and it'd be nice to be able to have an index on created_date and complaint_type.


title: Python libraries for correlating Photo to GPX file?

text: Can anyone recommend a Python library(s) that can geotag a photograph with the lat/long coordinates by correlating the timestamp of the photograph to the timestamp of a waypoint within a GPX file? I am aware of desktop applications like GeoSetter which will do this for you, but I would like to write an automation procedure that I could run on a nightly basis that emulates this process:   1. User takes a GPS and a camera into the field.   2. User drops a folder onto a server with a GPX file and the photos taken that day.   3. A Python process geotags the metadata each photo in the folder with the lat/long coordinates by comparing the photo timestamp to the nearest equivalent timestamp (+/- 30 seconds?) of a waypoint in the GPX file.   4. Convert the coordinates found in the photo's metadata to a point in a GIS system. I have procedures in place for steps 1, 2, & 4 but either need to find a library that already accomplishes step 3 or write my own. Any suggestions?


title: How to achieve Gradient Fills in ArcMap?

text: Does anyone know of a native way to achieve the gradient fills like the ones in the attached picture (map_1)? The circular and rectangular fills that I see already in ArcGIS' symbology choices don't seem to take the whole polygon into consideration (map_2). Is this type of thing done using another software package like Photoshop? ![map_1](http://i.stack.imgur.com/Ucxz6.jpg) ![map_2](http://i.stack.imgur.com/nJb7S.jpg)