title: Limit what fields can be edited GeoExplorer OpenGeoSuite

text: I have an issue with limiting what fields can be edited by the user. If I link my layer to my PostGIS table, the user is able to edit all the fields for a specific point. -If I create a view I can limit what they can see, however views are read only? Any obvious solutions I am not thinking about to allow me to limit which fields they can change or fill in when modifying a layer? Using OpenGeo Suite, GeoExplorer. Thanks


title: geoext popup + geoserver content.ftl

text: I am trying to include a jpg to a geoext popup. There are various shapefiles, but only one needs to load pictures. With only openlayers it was easy to do, I just used the content.ftl from geoserver 2.4. However, we switched to geoext and sadly I am not able to figure out how to do it now. Can somebody give me a hint how to force GeoExt.Popup to use the different content.ftl of geoserver - or is there another way?


title: Polygon to raster without losing any polygons

text: I'm trying to convert polygons to raster. But the output raster dataset doesn't have polygons which went into as input data. I need all the date in the attribute to follow with to the raster dataset as well. What am I doing wrong? I'm using ARCGIS 10.1.


title: Visualise data belonging to an area

text: I am looking for a way in QGIS (1.8.0) to visualise data, which is related to an area (polygon on a shape file). There are no x-y-coordinates with the data, just the info of the belonging area, means at an arbitrary position inside the area. Example: Samples x1,x2 and x3 are located in area a1, y1 to y5 in area a2. As the samples have certain properties on which I intend to do queries on, distinct locations seems somehow necessary (e.g. sample x1 and x3 have a certain property and therefore should be shown/highlighted ). I initially thought of assigning each sample some random coordinates from inside the areas, but what a hassle.. So, how to approach that??


title: What are the best ways to connect two coordinates on mapinfo?

text: I have with the help of this forum drawn a coordinate system with longitude and latitude at different location. I am trying to connect these sites what is the best way to connect. Is it just drawing lines or something else?


title: Geoserver not processing the data from drupal

text: I am using drupal cms to render the maps. I am using geoserver since I need to render more than say 50000 points. I am following the tutorial given in http://drupal.org/node/845852. I am getting stuck up as geoserver is able to read the values present in drupal database through `http://localhost/drupal6/wfspath?request=GetCapabilities` but unable to process them, I came to this conclusion because I am not able to see the point image in layer preview->openlayers and the layer preview->GML file is blank........... I am attaching wfs_data.png - gml at drupal(http://localhost/drupal6/wfspath) wfs_GetCapabilities.png - `http://localhost/drupal6/wfspath?request=GetCapabilities` given in geoserver store test_gml.png - geoserver->layer preview->gml test_OL.png - geoserver->layer preview->openlayers I badly need help on this, Thanx in advanceenter image description here ![enter image description here](http://i.stack.imgur.com/g2MZU.png) ![enter image description here](http://i.stack.imgur.com/cEOZd.png) ![enter image description here](http://i.stack.imgur.com/2HV3U.png) ![enter image description here](http://i.stack.imgur.com/Hq9p9.png) This is my Log file but I cannot understand the error               2011-05-10 04:03:16,599 DEBUG [geoserver.filters] - filtering http://localhost:8080/geoserver/wms     2011-05-10 04:03:16,602 DEBUG [geoserver.ows] - Getting layers and styles from LAYERS and STYLES     2011-05-10 04:03:16,603 DEBUG [geotools.styling] - number of fts set 1     2011-05-10 04:03:16,603 DEBUG [geoserver.ows] - establishing point style for testws:drupal_testwfsview     2011-05-10 04:03:16,607 DEBUG [wms.responses] - setting up map     2011-05-10 04:03:16,624 DEBUG [wms.responses] - asking delegate for write to org.geoserver.ows.DispatcherOutputStream@12e1a27     2011-05-10 04:03:16,918 DEBUG [geoserver.filters] - filtering http://localhost:8080/geoserver/wms     2011-05-10 04:03:16,923 DEBUG [geoserver.ows] - Getting layers and styles from LAYERS and STYLES     2011-05-10 04:03:16,924 DEBUG [geotools.styling] - number of fts set 1     2011-05-10 04:03:16,924 DEBUG [geoserver.ows] - establishing point style for testws:drupal_testwfsview     2011-05-10 04:03:16,927 DEBUG [wms.responses] - setting up map     2011-05-10 04:03:16,932 DEBUG [wms.map] - setting up 660x330 image     2011-05-10 04:03:16,932 DEBUG [geotools.styling] - number of fts set 1     2011-05-10 04:03:16,934 DEBUG [renderer.shape] - Affine Transform is AffineTransform[[1.422222222222222, 0.0, 330.0], [0.0, -1.422222222222222, 165.0]]     2011-05-10 04:03:16,935 DEBUG [geotools.rendering] - Computed scale denominator: 2.795411320143589E8     2011-05-10 04:03:16,936 DEBUG [geotools.styling] - number of fts set 1     2011-05-10 04:03:16,936 DEBUG [geotools.rendering] - processing 1 stylers for http://drupal.org/project/wfs:drupal_testwfsview     2011-05-10 04:03:16,936 DEBUG [geotools.rendering] - creating rules for scale denominator - 279,541,132.014     2011-05-10 04:03:16,937 DEBUG [geotools.rendering] - Expanding rendering area by 6 pixels to consider stroke width     2011-05-10 04:03:16,937 DEBUG [geotools.rendering] - Querying layer drupal_testwfsview with bbox: ReferencedEnvelope[-236.25 : 236.25, -120.234375 : 120.234375]     2011-05-10 04:03:16,945 WARN [data.wfs] - Query does not provide a CRS, using default: Query: [liteRenderer]        feature type: drupal:testwfsview        filter: org.geotools.renderer.lite.FastBBOX@3bb53f        [properties: geometry]     2011-05-10 04:03:16,955 DEBUG [geotools.xml] - [ class org.geotools.geometry.DirectPosition2D ] is not of type class [Ljava.lang.Double;     2011-05-10 04:03:16,956 DEBUG [geotools.xml] - [ class org.geotools.geometry.DirectPosition2D ] is not of type class [Ljava.lang.Double;     2011-05-10 04:03:17,588 DEBUG [data.wfs] - response head: <?xml version="1.0" encoding="utf-8" ?><wfs:FeatureCollection        xmlns:drupal="http://drupal.org/project/wfs"        xmlns:wfs="http://www.opengis.net/wfs"        xmlns:gml="http://www.opengis.net/gml"        xmlns:ogc="http://www.opengis.net/ogc"        xmlns:xsi="http:/          2011-05-10 04:03:17,593 DEBUG [geotools.rendering] - Style cache hit ratio: NaN , hits 0, requests 0     2011-05-10 04:03:17,593 DEBUG [renderer.shape] - Style cache hit ratio: NaN , hits 0, requests 0     2011-05-10 04:03:17,593 DEBUG [wms.responses] - asking delegate for write to org.geoserver.ows.DispatcherOutputStream@9d6ff7     2011-05-10 04:03:17,594 DEBUG [map.png] - Writing png image ...     2011-05-10 04:03:17,594 DEBUG [geotools.image] - Encoded input image for png writer     2011-05-10 04:03:17,594 DEBUG [geotools.image] - Getting a writer     2011-05-10 04:03:17,595 DEBUG [geotools.image] - Setting write parameters for this writer     2011-05-10 04:03:17,595 DEBUG [geotools.image] - Writer is NOT native     2011-05-10 04:03:17,596 DEBUG [geotools.image] - About to write png image     2011-05-10 04:03:17,639 DEBUG [map.png] - Writing png image ... done!     2011-05-10 04:03:18,222 DEBUG [geoserver.filters] - filtering http://localhost:8080/geoserver/ows     2011-05-10 04:03:18,224 INFO [geoserver.wfs] -      Request: getServiceInfo     2011-05-10 04:03:18,241 DEBUG [geoserver.requests] - Query is net.opengis.wfs.impl.QueryTypeImpl@19c9b66 (group: [], propertyName: [], function: null, filter: null, sortBy: [], featureVersion: null, handle: null, srsName: null, typeName: [{http://drupal.org/project/wfs}drupal_testwfsview])      To gt2: Query:        feature type: drupal_testwfsview        filter: Filter.INCLUDE        [properties:  ALL ]     2011-05-10 04:03:18,242 INFO [geoserver.wfs] -      Request: getFeature         handle = null         service = WFS         version = 1.0.0         baseUrl = http://localhost:8080/geoserver/         providedVersion = null         query = [net.opengis.wfs.impl.QueryTypeImpl@19c9b66 (group: null, propertyName: null, function: null, filter: null, sortBy: null, featureVersion: null, handle: null, srsName: null, typeName: [{http://drupal.org/project/wfs}drupal_testwfsview])]         maxFeatures = 50         outputFormat = GML2         resultType = results         traverseXlinkDepth = null         traverseXlinkExpiry = null         formatOptions = {}     2011-05-10 04:03:18,254 WARN [data.wfs] - Query does not provide a CRS, using default: Query:        feature type: drupal:testwfsview        filter: Filter.INCLUDE        [properties: geometry, body, nid, title]     2011-05-10 04:03:18,892 DEBUG [data.wfs] - response head: <?xml version="1.0" encoding="utf-8" ?><wfs:FeatureCollection        xmlns:drupal="http://drupal.org/project/wfs"        xmlns:wfs="http://www.opengis.net/wfs"        xmlns:gml="http://www.opengis.net/gml"        xmlns:ogc="http://www.opengis.net/ogc"        xmlns:xsi="http:/          2011-05-10 04:03:18,896 DEBUG [geotools.gml] - closing reader org.geotools.feature.collection.DelegateFeatureIterator@17c555d     2011-05-10 04:03:23,541 DEBUG [geoserver.filters] - filtering http://localhost:8080/geoserver/web/     2011-05-10 04:03:23,557 DEBUG [geoserver.web] - cannot convert url: jar:file:/C:/Program%20Files%20(x86)/GeoServer%202.0.3/webapps/geoserver/WEB-INF/lib/web-core-2.0.3.jar!/org/geoserver/web/css/blueprint/screen.css to file (URI is not hierarchical), falling back to the inputstream for polling     2011-05-10 04:03:23,564 DEBUG [geoserver.web] - cannot convert url: jar:file:/C:/Program%20Files%20(x86)/GeoServer%202.0.3/webapps/geoserver/WEB-INF/lib/web-core-2.0.3.jar!/org/geoserver/web/css/blueprint/print.css to file (URI is not hierarchical), falling back to the inputstream for polling     2011-05-10 04:03:23,567 DEBUG [geoserver.web] - cannot convert url: jar:file:/C:/Program%20Files%20(x86)/GeoServer%202.0.3/webapps/geoserver/WEB-INF/lib/web-core-2.0.3.jar!/org/geoserver/web/css/geoserver.css to file (URI is not hierarchical), falling back to the inputstream for polling     2011-05-10 04:03:23,571 DEBUG [geoserver.web] - cannot convert url: jar:file:/C:/Program%20Files%20(x86)/GeoServer%202.0.3/webapps/geoserver/WEB-INF/lib/web-core-2.0.3.jar!/org/geoserver/web/css/blueprint/ie.css to file (URI is not hierarchical), falling back to the inputstream for polling     2011-05-10 04:03:23,578 DEBUG [geoserver.web] - cannot convert url: jar:file:/C:/Program%20Files%20(x86)/GeoServer%202.0.3/webapps/geoserver/WEB-INF/lib/web-core-2.0.3.jar!/org/geoserver/web/css/geoserver-ie.css to file (URI is not hierarchical), falling back to the inputstream for polling     2011-05-10 04:03:23,584 DEBUG [geoserver.web] - cannot convert url: jar:file:/C:/Program%20Files%20(x86)/GeoServer%202.0.3/webapps/geoserver/WEB-INF/lib/web-core-2.0.3.jar!/org/geoserver/web/js/jquery-1.2.6.min.js to file (URI is not hierarchical), falling back to the inputstream for polling     2011-05-10 04:03:23,596 DEBUG [geoserver.web] - cannot convert url: jar:file:/C:/Program%20Files%20(x86)/GeoServer%202.0.3/webapps/geoserver/WEB-INF/lib/web-core-2.0.3.jar!/org/geoserver/web/js/jquery.inline-info.js to file (URI is not hierarchical), falling back to the inputstream for polling


title: How to install PostGIS Manager and RT Sql Layer plugins in QGIS 2.4?

text: This might be a rather simple question, but both plugins are not visible in my QGIS menu. I'm running 2.4 on Ubuntu 14.04. Thanks for every piece of advice.


title: ESRI projection to SRID?

text: I am looking at MapX and MapY coordinates. The docs say they are "truncated State Plane coordinates: NAD 1983 State Plane Louisiana South FIPS 1702 (US Feet)." I am trying to reproject into a common lat/long like wgs-84 by using PostGIS so I need the srid for this projection. I found the spatial reference page for the state plane coordinates http://spatialreference.org/ref/esri/102682/ -- but it does not list the SRID. I found this page of common SRIDs -- but if I plug in the SRID the projection does not output correct lat/longs. How can I find an authoritative SRID?


title: How can I place a couple of html elements on the map?

text:
I need to place html elements for drawing in it some text and graphic by Raphaeljs. Is it possible? I'am trying to create an empty marker for each visible feature and then erase icon container and place into it raphael canvas:               var markers = new OpenLayers.Layer.Markers('Markers');     map.addLayer(markers);          var size = new OpenLayers.Size(1, 1);     var offset = new OpenLayers.Pixel(0, 0);     var icon = new OpenLayers.Icon('/i/b.gif', size, offset);          var marker_items = {}, current_items = {}, centroid, x, y, key, marker, $marker;          layer4.events.register('moveend', layer4, function(e){         current_items = {};         for (var i=0, n=layer4.features.length; i<n; i++) {             centroid = layer4.features[i].geometry.getCentroid();              x = centroid.x, y = centroid.y;             key = 'll-'+x+'-'+y;             current_items[key] = true;             if (marker_items[key])                 continue;                  // create new marker and store it in hash by key             marker_items[key] = new OpenLayers.Marker(new OpenLayers.LonLat(x, y), icon.clone());             markers.addMarker(marker_items[key]);                  // erase icon container, add canvas, and draw raphael graphics             $marker = $(marker_items[key].icon.imageDiv).empty();             $canvas = $('

').appendTo($marker);             Raphael($canvas.get(0), 70, 70).piechart(35, 35, 35, [55, 20, 13, 32, 5, 1, 2, 10]);         }              // remove markers that not visible         for (i in marker_items) {             if (!current_items[i]) {                 markers.removeMarker(marker_items[i]);                 delete marker_items[i];             }         }     });      I have successfully draw graphics, but maybe there is more simple way?


title: How to create Extended Data for KML exported from shapefile by ArcGIS for Desktop?

text: I have a GIS shapefile of 300+ points that I am **converting to KML using ArcGIS for Desktop**. Once in KML I want to create custom balloons by editing the KML file. However, the HTML feature in ArcGIS does not set up the Schema and Extended Data that I need (e.g. if I open a shapefile in Google Earth Pro and save it as KML, it creates the Schema and Extended Data automatically). I don't have access to Google Earth Pro right now. Is there another way to create this in ArcGIS or is there a 3rd party program that will do this? I don't have much programming experience, so I can't really go that route. Thanks in advance!


title: Spatialite: Convert polygon to multi-polygons based on common field

text: I have a list of polygons which should be multipolygons:               PK_UID | PolygonID | Geometry     ------------------------------     1      | 1         | polygon     2      | 2         | polygon     3      | 1         | polygon      How do I convert them to multipolygons in Spatialite based on the common field (PolygonID)? (e.g So they look like this):               PK_UID | PolygonID | Geometry     -----------------------------------     1      | 1         | multi-polygon     2      | 2         | polygon


title: How to improve poor performance times for PostGIS 2.1 Tiger Geocoder?

text: I've recently upgraded my PostGIS DB from 2.0 to 2.1. I initially followed the instructions from Waiting for PostGIS 2.1 - Install PostGIS Tiger Geocoder as an Extension and 2.7. Installing, Upgrading Tiger Geocoder and loading data to load (for the first time) the 2012 Tiger Geocoder into the DB by grabbing the updated geocoder files from the postgis-pg92-binaries-2.1.1devw32.zip file. Then I downloaded and installed the data for 5 western states. (OR, WA, ID, CA, MT). Finally, I ran the `SELECT install_missing_indexes();` command after loading all of the data. I also ran `vacuum analyze` on the whole DB afterwards. That all went fine, and the geocoding process is working. However, while I am admittedly running this DB on a not-so-speedy laptop, my results of 2-8 sec results versus the 61ms results noted here seem a way out of whack. **Any ideas on what I can do to improve performance?** Guides I worked through to try to improve performance:   * OpenGeo "Section 29: Tuning PostgreSQL for Spatial"   * OpenStreetMap "Tune the database "   * Part 1: Getting Started With PostGIS: An almost Idiot's Guide (PostGIS 2.0)   * Best practices for PostGIS production environments? and I've updated all of my postgresql.conf settings to match the highest settings from any of the documents. **System specs** : Win 7 64 bit - Intel 2.3Ghz, 4GB RAM, Paging space 2000MB to 8000MB (fixed), PostgreSQL 9.2 (32 bit), Postgis 2.1 **Additional Details** : When I run the following query using "explain analyze" I get the following results:               SELECT g.rating, ST_X(g.geomout) As lon, ST_Y(g.geomout) As lat,          (addy).address As stno, (addy).streetname As street,          (addy).streettypeabbrev As styp, (addy).location As city, (addy).stateabbrev As st,(addy).zip          FROM geocode('<house #> <street name>, <city> <state> <zip>') As g;           Function Scan on geocode g  (cost=0.25..15.25 rows=1000 width=68) (actual time=185.606..185.606 rows=1 loops=1)


title: CartoDB Vis - Multiple layers mouse interaction?

text: My map visualization has two layers: 1) A city layer where the cities are represented as standard markers. 2) A county layer underlying the cities which has county level aggregate data. The layer is a choropleth. -- Can we get mouseover and / or mouseclick events working on both layers simultaneously? I need to be able to click/hover a city or a county and get the popup data response. (and probably highlight the county polygon) At this point I can access one or the other depending on their order in the vis. After looking around a bit it looks like at least at some point this was handled by a variety of workarounds. I was hoping to get it going somewhat automagically. If not, is this template the "right" way of doing it? I could stuff all the data into one table but that feels like it will make more complication. Users will be entering new marker points into the DB and it would be preferable to keep that data on its own. https://github.com/mhkeller/cartodb-templates/blob/master/basic/point-polygon- hover.html Appreciate your help, I'm in the GIS deepend again.


title: smoothing results from scipy griddata interpolatioin

text: I am using the griddata interpolation package in scipy, and an extrapolation function pulled from fatiando:               import numpy as np     import scipy     from scipy.interpolate import griddata     import matplotlib.pyplot as plt          def extrapolate_nans(x, y, v):         '''           Extrapolate the NaNs or masked values in a grid INPLACE using nearest         value.              .. warning:: Replaces the NaN or masked values of the original array!              Parameters:              * x, y : 1D arrays             Arrays with the x and y coordinates of the data points.         * v : 1D array             Array with the scalar value assigned to the data points.              Returns:              * v : 1D array             The array with NaNs or masked values extrapolated.         '''              if np.ma.is_masked(v):             nans = v.mask         else:             nans = np.isnan(v)         notnans = np.logical_not(nans)         v[nans] = scipy.interpolate.griddata((x[notnans], y[notnans]), v[notnans],             (x[nans], y[nans]), method='nearest').ravel()         return v               grid_x, grid_y = np.mgrid[0:1.5:50j, 0:1.2:50j]     points = np.random.rand(50, 2)     values = np.random.random_integers(1,10,50)          x = []      y = []     for i in points:         x.append(i[0])         y.append(i[1])     n = plt.Normalize(values.min(), values.max())               grid_z = griddata(points, values, (grid_x, grid_y), method='cubic')     extrapolate_nans(grid_x,grid_y,grid_z)     plt.contourf(grid_x,grid_y,grid_z)     plt.scatter(x,y,c=values)     plt.colorbar()          plt.show()      I think the results are decent, but the problem is that the extrapolated values have a distinct zigzag pattern: ![ZigZagging contour lines](http://i.stack.imgur.com/L0K79.png) **How do i smooth these lines? Or is there an alternative interpolation method in python that will generate smoother results?**


title: QGIS 1.8.0. How to resolve map request error. Import of WMS layer

text: When attempting to import a WMS layer into an existing QGIS project, the following error message appears: Response: msDrawMap(): Image handling error. Failed to draw layer named 'WTR_RSTRCTN_PNTS'. msSDELayerOpen(): SDE error. msSDEGetLayerInfo(): . (1) msSDEGetLayerInfo(): SDE error. SE_layer_get_info(): Layer not found. (-20) msSDELayerOpen(): SDE error. msSDEGetLayerInfo(): . (1) msSDEGetLayerInfo(): SDE error. SE_layer_get_info(): Layer not found. (-20) Previous maps imported from the same catalogue have worked without a problem. Is there a fix for this error? Thanks


title: How to configure osm2po to use the specified maxCost in TspDefaultMatrix

text: **The problem** I am using _osm2po_ 's `TspDefaultMatrix` to create a cost (time) matrix:               new TspDefaultMatrix(graph, vertexIds, maxCost, Log.stdout(Log.LEVEL_LOG), properties)      Where `properties` contains:               findShortestPath: false     ignoreRestrictions: false     ignoreOneWays: false     heuristicFactor: 1.0      I tried many value for `maxCost` (`0`, `1`, `2`, `Float.MAX_VALUE`). Except for `maxCost = 0`, I always get this in the log:               DEBUG  Heuristic MaxCost set to 0.2134047      followed by some:               WARN   Matrix [8][9] unreachable!     WARN   Matrix [8][3] unreachable!      because, indeed, some of the vertexId are more than `0.2134047` hours appart. **The question** I want to know how this value (`0.2134047`) is calculated and how I can effectively tell `TspDefaultMatrix` to use my `maxCost`. By the ways, I did not found (was not searching a lot) the osm2po source code. Is it available ? I would have been able to answer my question myself. **EDIT:** By the way, I was using version `4.8.8` of _osm2po_ to build the `.gph` file and at runtime.


title: Fill tool creating large output files from sections of a Raster

text: ![Image of data from one ofthe clipped rasters, uncompressed size does not match what the size in my drive says](http://i.stack.imgur.com/PrdSg.jpg) (This raster information for one of the clipped rasters) I have a raster image representing a landscape that is about 1.3Gb. I clipped the raster into 31 sections which vary in size from 10Mb to 40Mb. Whenever I use the Fill tool on a section, the output raster is 1.3Gb (like the whole original raster). This is too much memory to process. Why is the Fill tool creating such a large output file for just a section of the original raster? Is there a way to reduce this size? Edit: Additional Information The Rasters were clipping using the methods in this previous question Clip a raster using a shapefile into multiple rasters? Edit 2: Python Script changing the mask environment before using Fill (It works, but still provides the 1.3Gb "Fill" Raster)               import arcpy     from arcpy import env     from arcpy.sa import *     arcpy.env.workspace = (r"Workspace\Path")              for i in range (31):                  arcpy.MakeFeatureLayer_management("shapefile","featurelayer" + str(i),'"FID"='+str(i))                   arcpy.env.mask = ("featurelayer"+str(i))                  arcpy.ClipManagement(“fullextentraster”, “#”, raster +“clip”+str(i)+ “.tif”, “file_name” + str(i), “0”, “ClippingGeometry”)                   outFill = Fill(“rasterclip”+str(i)+“.tif”)             outFill.save (r“Save\Location\Folder\rasterfill”+str(i)+ “.tif”)       Edit 3: Python Script changing the processing extent environment before using Fill               import arcpy     from arcpy import env     from arcpy.sa import *     arcpy.env.workspace = (r"Workspace\Path")              for i in range (31):                  arcpy.MakeFeatureLayer_management("shapefile","featurelayer" + str(i),'"FID"='+str(i))                   arcpy.env.extent= ("featurelayer"+str(i))                  arcpy.ClipManagement(“fullextentraster”, “#”, raster +“clip”+str(i)+ “.tif”, “file_name” + str(i), “0”, “ClippingGeometry”)                   outFill = Fill(“rasterclip”+str(i)+“.tif”)             outFill.save (r“Save\Location\Folder\rasterfill”+str(i)+ “.tif”)


title: How to convert Lat/Long to UTM using Proj(4j) similar to JScience UTM.latLongToUTM?

text: I'd like to use Proj(4j) to convert between UTM and Lat/Long coordinates. I used to use JScience for UTM <-> LatLong conversions like this:               UTM.latLongToUtm(latlong, ReferenceEllipsoid.WGS84);     UTM.utmToLatLong(coordinate, ReferenceEllipsoid.WGS84);      But I would like to accomplish this using Proj4j. So I'm currently having set up two `CoordinateReferenceSystems`:               CRS_FACTORY.createFromParameters("EPSG:4326", "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs");     CRS_FACTORY.createFromParameters("?", "+proj=utm +? +zone=?);      Where I don't know how to fill up the ? gaps. Why do I have to specify a zone-parameter? Is there a way to accomplish the same above in JScience using proj4j?


title: Convert Labels to Annotation greyed out

text: This is probably easy, but haven't found a solution. Convert Labels to Annotation is greyed out on right clicking a feature in the TOC. I have no reference scales set. Thanks.


title: Processing SQLite files using ArcGIS Desktop?

text: I'm using qgis and grass for my work. I've produced a shapefile which is over 100mbs (s a result of some statistical analysis). If i convert it to sqlite format (spatially enabled) the data shrinks down to 15 mbs. Which is acceptable for sharing it with email or dropbox. BUT my colleague is using arcgis 9.3 - will he be able to read the data without any problems?


title: How to configure QGIS Sextante to use SAGA?

text: I have been trying to set up SAGA to work using Sextante in QGIS. I have followed the instructions to configure SAGA in Sextante, with the following input in the path for the SAGA folder C:\Program Files\SAGA-GIS\saga_gui. When I then try to use a SAGA algorithm I get this message:Saga folder is not configured. Please configure it before running SAGA algorithms. When I go back to the Sextante configuration options, the path for the Saga folder is empty. There is no information under the Sextante history log. Anyone got any ideas about how I can resolve this issue?


title: gdal2tiles colorpallete error

text: merged few files using gdal_merge gdal_merge.py -pct -o 713.tif sa713.tif na713.tif wa713.tif ea713.tif When i try to make tiles it says gdal2tiles.py: error: Please convert this file to RGB/RGBA and run gdal2tiles on the result. From paletted file you can create RGBA file (temp.vrt) by: gdal_translate -of vrt -expand rgba 713.tif temp.vrt then run: gdal2tiles temp.vrt . How to make tiles of this raster with out loosing the colorpallete data , i can do this using tilemill , but i want to do this on command line . Any help ?


title: How to delete duplicate lines?

text: I was wondering if anyone could help. I have a tab file in MapInfo consisting of just polylines. Each is meant to have its own reference number, however some are duplicated (e.g two lines occupying the same position with the same ID's) I can identify these with a query, but what I need to know is if it's possible to delete duplicates and leave one remaining (it doesn't matter which one as they are identical). Thanks, Alan (ps. I also have access to QGIS if it is easier in that)


title: Is there a programmatic way to remove a user from a license

text: I'm new to the Enterprise level of the GIS world and working at this level I find that there are many instances where another user has taken an ArcEditor or ArcInfo license when it's needed for a more important situation that has arisen. My question has two parts. * * *   1. Is there a way to remove a user from a license programmatically?   2. Is there a way to transfer a user from one license-type to another? (i.e., ArcEdit to ArcView) _Programmatically would include Python, VBA, VB, or other._ I am familiar with "Where Have All The Licenses Gone" but it's very inconvenient as sometimes we have users that leave their computers on and then leave for the day. Rather than unplugging their systems and potentially losing their work it would be a lot more efficient to just strip a specific license from them or transfer their current license to another license.


title: How to bypass 10 character limit of field name in ShapeFiles?

text: I'm exporting geometry with attached text attributes from oracle database to esri shapefile format (.shp) with Java and Geotools library. Attribute columns in our database have names with more than 10 characters, and Geotools forces to truncate them. I understand that it's due to .shp or .dbf files specification. I can workaround this, by creating some simple txt file with "shrtname"="The full and long name", but obviously it won't be understand and imported by any other software than ours. Is there official way to do the mapping from short field names to long full- text names? For example xml file next to all other .shp .dbf .shx files? What are the rules then? greetz


title: QGIS QgsPluginLayer broken on 2.0.1?

text: Is there a generic issue with the QgsPluginLayer class in QGIS 2.0.1? When I take the basic Watermark sample from here: https://github.com/sourcepole/qgis- watermark-plugin and make a couple of mods so that the plugin runs on 2.0.1 (all that appeared to be required was to add a metadata.txt file) it fails with an error "Could not draw Watermark plugin layer because:" and thats it. It seems like there is a similar issue active with an OpenLayers plugin as reported here: http://hub.qgis.org/issues/8679 but I am working on a plugin using the QgsPluginLayer class which is also now broken so I am wondering whether anyone can suggest a workaround or have any idea about when a fix may be available? Thanks in advance!


title: How to switch Map Units on Data Frame Properties back to Feet in ArcMap?

text: Using ArcView 10.1, I can not switch the Map (Units) nder Data Frame Properties to Feet. It is currently at Decimal Degrees. Display can switch. How do you switch it back to Feet?


title: Analyzing Raster Data From Vector Layer Buffer

text: I am using QGIS 2.0. I have satellite images of night-time lights all over earth. I also have a vector shapefile of country boundaries including the largest cities on earth. I want to use buffers to capture a certain amount of light pollution emanating from cities that is viewable from the raster image. My question is, how do I use buffers to analyze pixels (or some aspect of the raster data) within the buffer area? I assume I will want to assign some value to a certain pixel and then total those values by using the buffer to analyze all data within the buffer. I have attached a snip of what I am working with. I am very new to GIS still, so bear with me. Thanks in advance. ![Example](http://i.stack.imgur.com/LLvuB.jpg)


title: Creating a layer in arcmap by dragging a table from an OracleDB

text: I am working on a project for a client, they gave us some data from their oracleDB, and when I connect to it in arc map I can just drag and drop the tables from the catalog window into arcmap and they are drawn. Now we have a bunch of legal deeds that I then extract the data from using OCR, then from the text I programmatically create a shape file and insert it into a table in the same database. However when I drag that table(which is in the same DB with the data they gave me, which I can just drag and drop to draw) I get the warning that says "Could not add the specified data object to the map", and it is not drawn. However when I right click on the table from the arc catalog window and hit create layer, then save the layer, then add the layer from that file that was just created, everything works fine and it is drawn. Am I missing something? Or is something not set up properly. I would like to have the ability to drag and drop for the client as their data does that they already have. I admit I am very new to acrGIS and the esri suite, this is the first gig I have been exposed to it, and I feel like I'm missing something simple and dumb. I've tried googling but haven't found a solution that works.


title: Arcgis Server name include domain

text: When I create a web app using the template it points to `cos-t420:8399/arcgisoutput` but I need it to be `cos-t420.web.com.sg:8399/arcgisoutput`   How can it be done? am using the java version.


title: Choose what zoomlevel to read from for mbtiles file in GDAL

text: GDAL is supposed to read MBTiles files and to do so, it selects the highest zoomlevel available. However, is there a way to choose another zoomlevel? i.e. override the default and go for a lower one? I'm currently accessing the file like this (and it defaults to the highest zoomlevel):               gdal_translate tiles.mbtiles export.tiff


title: Need advice on creating polyline templates (polyline vs polyline ZM) for future network dataset

text: I am about to request large amounts of polylines from a variety of sources (municipalities). The idea is to collect bike path data and build a regional bike network (I work for the region administrative board). To help municipalities (and myself once i need to assemble all the data), I plan to create template files. I now realize I do not know enough about network datasets and 2D to 3D conversions. I basically dont know what I dont know, hence this very open question. I have:   * Arcgis 9.3 (hopefully 10.2.1 in the fall but not for certain).    * High resolution DEM of the region.   * limited bike path shapefiles (segmented polyline ZM, not a network dataset) I need:   * advice in creating template shapefiles that are easy and flexible enough for the municipalities to use and still be fairly straight forward for me to put together for a network. The templates will contain 3 attributes (short integers) that denote road type, traffic situation and recreational value.  specifically I would like to know more about polyline vs polyline ZM? - bear in mind that the municipalities are using a variety of software (CArcGIS, MapInfo, CAD(?) etc. etc.). Thank you for any suggestions.


title: ArcGIS, clip issue when using a PGDB

text: When I clip PolylineMZ feature class inside a Personal Geodatabase it gives odd results for vertices. Vertices are horizontally shifted and do not corespond with original. When I do this in .shp, results are fine. Any ideas what is wrong?


title: How to avoid ImportError: No module named arcpy causing Failed to execute (Script) error?

text: I have a proplem with arc hydro in arcpy,i want to create a simple script where automated some tools from archydro in arcpy.. I use arcmap 10.1 and python 1.7,my import is one dem and my export is a 6 raster file (fill,fdr,fac,cat,str,strlink) and run ok without proplem but if I add one more tool cachmented to polygon that converts a raster file to polygon feature class then show me error i import this python script on arcmap               import arcpy          arcpy.ImportToolbox("C:\Program Files (x86)\ArcGIS\Desktop10.1\ArcToolbox\Toolboxes\Arc Hydro Tools.tbx", "archydrotools")     #Set Ouput Names     fill=arcpy.GetParameterAsText(0) + '\\fill'     Fdr=arcpy.GetParameterAsText(0) + '\\fdr'     Fac=arcpy.GetParameterAsText(0) + '\\fac'     Str=arcpy.GetParameterAsText(0) + '\\Str'     StrLnk=arcpy.GetParameterAsText(0) + '\\StrLnk'     Cat=arcpy.GetParameterAsText(0) + '\\Cat'     Catchment=arcpy.GetParameterAsText(0) + '\\Catchment'          # Local variables:     arcpy.CheckOutExtension("spatial")     arcpy.env.overwriteOutput = True     arcpy.env.workspace = arcpy.GetParameterAsText(0)     dem = arcpy.GetParameterAsText(1)     # Process: Fill Sinks     arcpy.FillSinks_archydrotools(dem, fill, "", "", "ISSINK_YES")     # Process: Flow Direction     arcpy.FlowDirection_archydrotools(fill, Fdr, "")     # Process: Flow Accumulation     arcpy.FlowAccumulation_archydrotools(Fdr, Fac)     # Process: Stream Definition     arcpy.StreamDefinition_archydrotools(Fac, "", Str, "1")     # Process: Stream Segmentation     arcpy.StreamSegmentation_archydrotools(Str, Fdr, StrLnk, "", "")     # Process: Catchment Grid Delineation     arcpy.CatchmentGridDelineation_archydrotools(Fdr, StrLnk, Cat)     # Process: Catchment Polygon Processing     arcpy.CatchmentPolyProcessing_archydrotools(Cat, Catchment)      ![my error](http://i.stack.imgur.com/0aCLU.jpg) but my export file in the folder they become export the all files regular my proplem is the last command because is a feature class I think any idea ?


title: ArcObjects vector projection using VB

text: I am creating a ArcGIS console application using ArcGIS 10.2 and Visual Studio 2010. I need to project a shapefile to a new location (.gdb) using console inputs and the Project tool within the Data Management toolbox. Code below:               Imports ESRI.ArcGIS.esriSystem     Imports ESRI.ArcGIS.Geoprocessing     Imports ESRI.ArcGIS.DataManagementTools     Module Module1     Private m_AOLicenseInitializer As LicenseInitializer = New preInGen.LicenseInitializer()          <STAThread()> _     Sub Main()         'load product code and version to version manager         ESRI.ArcGIS.RuntimeManager.Bind(ESRI.ArcGIS.ProductCode.Desktop)              'ESRI License Initializer generated code.         m_AOLicenseInitializer.InitializeApplication(New esriLicenseProductCode() {esriLicenseProductCode.esriLicenseProductCodeAdvanced}, _         New esriLicenseExtensionCode() {})              'get state abbreviation from user         Console.WriteLine("State abbreviation: ")         Console.Write(">")         Dim StAbb As String = Console.ReadLine()              'get FIPS from user         Console.WriteLine("FIPS: ")         Console.Write(">")         Dim FIPS As String = Console.ReadLine()              'create new geoprocesor         Dim gp As IGeoProcessor = New GeoProcessor         gp.OverwriteOutput = True              'create parameter inputs         Dim param As IVariantArray = New VarArray         param.Add("\\\\prhqgis01\\f\\arcgisserver\\MapMS2_Data\\Parcel_Boundaries\\" & StAbb & "\\p" & FIPS & ".shp")         param.Add("GCD_WGC_1984")         param.Add("\\\\dvhqgis02\\d\\David\\natdb\\" & FIPS & "\\" & FIPS & "_1.gdb\\p" & FIPS & ".shp")         param.Add("GCS_North_American_1983")         param.Add("WGS_1984_(ITRF00)_To_NAD_1983")              'project parcel file         gp.Execute("Project_management", param, Nothing)              'ESRI License Initializer generated code.         'Do not make any call to ArcObjects after ShutDownApplication()         m_AOLicenseInitializer.ShutdownApplication()     End Sub          End Module      When I do a test run, I get a generic COM exception: "Error HRESULT E_FAIL has been returned from a call to a COM component." Can anyone help me find the exact source of the error?


title: Multiple FeatureType in WFS Protocol

text: I would like to pass multiple featureTypes to WFS protocol in order to accomplish a search feature in my app. the idea worked just fine with one featureType, but now I need to add another one to expand my search scope. Here is what I did:               var featureTypes = ["Layer_1", "Layer_2"];            var wfsProtocol = new OpenLayers.Protocol.WFS.v1_1_0({                   url: URL,                                   featureType: featureTypes,                  featureNS: "http://mapmap.org",                   srsName: "EPSG:3857",                   version: "1.1.0",                   extractAttributes: true,                   isBaseLayer: false,                   visibility: true                     });               wfsProtocol.read({                     filter: filter,                               callback: processTheQuery,          scope: [new OpenLayers.Strategy.Fixed({preload: true})],              })      But still didn't work, Any idea?


title: Join multiple spatial tables by distance

text: I'm trying to join two tables, one which contains road names and the other just a point. After exhausting research I'm now stuck. The idea is to join up the two tables using only distance, ie, select all road names which is within 500 meteres of a postcode My current query does nothing except waste time and I'm wondering where I'm going wrong?               SELECT "R"."name", "P"."Postcode"          FROM "planet_osm_roads" AS "R" JOIN "gb_postcodes" AS "P"           ON st_dwithin(st_transform("P"."Geom", 900913), "R"."way" , 500)          ORDER BY "P"."Geom" <-> st_setsrid(st_makepoint(0.06703376770019531,  51.51646108300825), 900913)               LIMIT 10;


title: Adding FeatureLayer on a thread crashes ArcMap

text: I asked this question on the ESRI forums but got nothing so I thought I would try it here (hope that is consider ok). I am creating an extension for ArcMap 10.2 using Visual Studio 2012 C# WPF. Consider this code snippet:               //Open the CSV file     IWorkspaceFactory WorkspaceFactory = new TextFileWorkspaceFactoryClass();     IFeatureWorkspace FeatureWorkspace = (IFeatureWorkspace)WorkspaceFactory.OpenFromFile(FilePath, 0);     ITable Table = FeatureWorkspace.OpenTable(FileName);          IStandaloneTable StandaloneTable = new StandaloneTableClass();     StandaloneTable.Table = Table;          IMap Map = MxDocument.FocusMap;     IStandaloneTableCollection StandaloneTableCollection = (IStandaloneTableCollection)Map;     StandaloneTableCollection.AddStandaloneTable(StandaloneTable);          //Create the Event Layer of the data source     IDataset Dataset = (IDataset)StandaloneTable;     IName TableName = Dataset.FullName;          IXYEvent2FieldsProperties XYEvent2FieldsProperties = new XYEvent2FieldsPropertiesClass();     XYEvent2FieldsProperties.XFieldName = ColumnX;     XYEvent2FieldsProperties.YFieldName = ColumnY;     XYEvent2FieldsProperties.ZFieldName = String.Empty;          IXYEventSourceName XYEventSourceName = new XYEventSourceNameClass();     XYEventSourceName.EventProperties = XYEvent2FieldsProperties;     XYEventSourceName.EventTableName = TableName;     XYEventSourceName.SpatialReference = SpatialReference;          IName XYEventName = (IName)XYEventSourceName;     IXYEventSource XYEventSource = (IXYEventSource)XYEventName.Open();          //Map the Event Layer     IFeatureLayer FeatureLayer = new FeatureLayerClass();     FeatureLayer.FeatureClass = (IFeatureClass)XYEventSource;     FeatureLayer.Name = FileName;          //Add the layer extension     ILayerExtensions LayerExtensions = (ILayerExtensions)FeatureLayer;     LayerExtensions.AddExtension(new XYDataSourcePageExtensionClass());          //Add it to the MXD     Map.AddLayer(FeatureLayer);     FeatureLayer.Selectable = true;     MxDocument.UpdateContents();     MxDocument.ActiveView.Refresh();      This is from a function in a class that is called from the UI Form. What I am doing is creating a CSV file from a .Net DataTable (not in the above code), adding the CSV as a Table in Sources of the TOC, and then mapping the Table as Layer (the last 4 lines above). This all works fine if I do it on the main UI thread of the Form when the user clicks a button. But what I would like to do is it put it on a TPL thread since it can be a time consuming process and here is when I run into problems. The above code will CRASH arcmap completely. It cannot even be trapped. I noticed two things. One, the crash happens on the line "Map.AddLayer(FeatureLayer);" so when it tries to add to the layer collection. So I removed that line just to see what the Table in sources looks like. Second, looking at that Table I notice that if I right click it the options to Open and Display are lit-up for a second and then grey out. Again, this only happens when on a Thread. Another thing to note is that when I break on any of the lines of code above in the Thread and inspect the objects I get an error "The function evaluation requires all threads to run". I would think this is more an issue with Visual Studio and the implicit evaluation settings so it should not affect arcmap and it certainly should not cause it to crash. But I thought I would mention it. Thanks for any help.


title: Python loop Polygon to raster not working

text: I want to convert a large number of different shp polygons to raster in ArcGIS 10.1, using an existing field column "class" for the new raster classes. I am not very fluent in python but I thought this would work, and it doesn´t:                # Import arcpy module      import arcpy     ... # Local variables:     ... project_input = "D:\\RESEARCH_PROJECTS\\NORD-STAR\\AMAZONIA\\1-RAW_DATA\\INPE\\PRODES DIGITAL\\2004\\"     ... raster_output = "D:\\RESEARCH_PROJECTS\\NORD-STAR\\AMAZONIA\\1-RAW_DATA\\INPE\\PRODES DIGITAL\\2004\\convert\\"     ... # List features     ... arcpy.env.workspace = project_input     ... fcList = arcpy.ListFeatureClasses()      ... # Loop     ... for featureClass in fcList:     ... # Output     ...     raster_outputpath = raster_output + featureClass     ...      ... # Process: Polygon to Raster     ...     arcpy.PolygonToRaster_conversion(featureClass, "CLASS", raster_output, "CELL_CENTER", "NONE", "30")     ...     print "finished polygon to raster"      ...          Runtime error  Traceback (most recent call last):   File "", line 20, in    File "c:\program files (x86)\arcgis\desktop10.1\arcpy\arcpy\conversion.py", line 2436, in PolygonToRaster     raise e ExecuteError: ERROR 999999: Error executing function. Failed to execute (PolygonToRaster).            The script actually calls the first polygon, but that is it: > Executing: PolygonToRaster "D:\RESEARCH_PROJECTS\NORD- > STAR\AMAZONIA\1-RAW_DATA\INPE\PRODES > DIGITAL\2004\PDigital2004_00157_pol.shp" CLASS "D:\RESEARCH_PROJECTS\NORD- > STAR\AMAZONIA\1-RAW_DATA\INPE\PRODES DIGITAL\2004\convert\" CELL_CENTER NONE > 30 Start Time: Wed Nov 21 15:17:11 2012 ERROR 999999: Error executing > function. Failed to execute (PolygonToRaster). Failed at Wed Nov 21 15:17:11 > 2012 (Elapsed Time: 0,00 seconds) Thanks for your help Javier


title: Leaflet: How to move the layer-control menu?

text: This might be a silly question, but I could not find a documented way to achieve this. I'd like to freely position the layer-control menu, probably to the top left near the default zoom-in/out button. My layer control looks like this:               // Group layers as overlay pane     overlayPane = {       "Endpoints" : endpointMarkerLayer,       "Links" : linkLineLayer     };          // Add a layer control element to the map     layerControl = L.control.layers(null, overlayPane);     layerControl.addTo(map);      Where `endpointMarkerLayer` and `linkLineLayer` are layers containing markers and polylines respectively. Is there an option to specify where the menu should appear? Or alternatively, how can I get a reference to the DOM-objcet of the control-menu, such that I could assign it a custom class and override positioning in CSS?


title: Is there a way to control which ArcMap AddIn is loaded first?

text: I have a shared DLL that use across multiple ArcMap AddIns. It rarely changes, but recently I've had to modify it and I started getting MissingMethodExceptions whenever I made calls to classes in that shared DLL. After some investigation, I realized that there are other ArcMap AddIns that carried an older version of the shared DLL that were being loaded first. And if I clear out all of the other AddIns and blow away the assembly cache, I run by AddIn with new Shared DLL with no issues. Is there a specific order of AddIn loading from the assembly cache? In my case, it appears to be oldest first. Is there a way to alter the order of loading? Related Links:   http://forums.arcgis.com/threads/29870-Add-ins-that-share-same-custom-library


title: How to find most common classification within buffers?

text: I have a point shapefile that identifies the location of incinerators in England and an output area classification shapefile. I have created a buffer for all the points in the shapefile because i want to see what area classifications fall within the buffer of all the incinerators i have mapped. Next, I intersected the buffer i have created with the output area classification layer. What I want to know is how to i analyse the data within the buffers of the incinerators? I want to to know what is the most common classification throughout all the buffers....


title: How to select all dangling polyline features with Arcview license?

text: I have a city roads shp file and would like to identify (select) all the locations where a roads ends (in other words dangles). I am working in ArcView 10. I do not want to use a topological tool such as ETgeowizards because it can detect and dangles but fixes them automatically with no further options to identify exeptions or even visualize where fixes occured (in my case I will have many exeptions). ArcEditor provides a good interactive environment for identifying which errors are expections - but I am trying to see if I can do this with just ArcView. The approach I have been experimenting with is to try to identify/select segments in the network that only touch one other segment. I thought the dissolve tool might be an option, but it doesnt look like it has the settings I want in the "unsplit lines" option...


title: In ArcMap, can polygon labels always be displayed in a fixed area of the screen?

text: If I'm inside an area, I want its properties to be displayed in a fixed area of the screen (top right, let's say).


title: Open dataset(s) with centroids or other geometry for WOEIDs?

text: The Yahoo! GeoPlanet data is available for download but this data only includes the object hierarchies without any associated geometry. Do anyone know if there is an open source for geometry (centroids or otherwise) for associated WOEIDs? I wonder if there's a lesser-known link somewhere on the GeoPlanet site, or perhaps another party has done some work to correlate open data sets (census, geonames, etc) with WOEIDs. A WOEID (Where on Earth IDentifier) is a unique reference identifier assigned by Yahoo! to identify any feature on Earth.


title: ArcGIS Server 10 Rest Admin login

text: I am trying to login to the http://myserver/arcgis/rest/admin page and when I log in with credentials I get no error message but I am returned to the same login page. If I put in bad credentials I get an error message saying bad login. I can login to the rest/services page and the ArcGIS Server Manager Site. I found this but I am not sure if this applies http://forums.esri.com/Thread.asp?c=204&f=2388&t=257171 It looks like I am using a reverse proxy..from looking in the rest.config file under the rest directory. Does anyone have a clue on what I should investigate?


title: Looking for ArcGIS API for javascript and Jquery samples

text: I'm new in ArcGIS javascript API .I found tons of samples of ArcGIS javascript API and DOJO . it was very interesting for me but i know Jquery and i'd like to use Jquery to develop webgis.Is there any resources and samples about using ArcGIS javascript API and Jquery ? Thanks for advance.


title: Can QGIS preserve layers when exporting a PDF?

text: I saw this posted about a year ago (Can QGIS preserve layers when exporting to PDF?) but figured it was worth asking again -- is this possible? My usual workflow with ArcGIS is to export to AI then edit the layers that way, is there some way to do that from QGIS? Now I just get everything in one huge layer. Thanks.


title: Arguments for opening a GIS area in an important agricultural-related bank?

text: the managers of the bank where I work at have been considering creating a "GIS-unit" that would help with economic-related research / business inteligence. They asked me to raise a list of reasons of why/how such group would benefit the bank businesses. So, my friends, I pass this question to you all. Which arguments can you think of for making the case for the bank to create such group? Just a note, this bank is strong in agricultural-related investments. TL/DR: How can GIS benefit an agricultural-bank businesses??? Thank you.


title: What is the difference between Vector and Raster data models?

text: From: http://support.esri.com/en/knowledgebase/GISDictionary/term/vector%20data%20model > _**vector data model_** : [data models] A representation of the world using > points, lines, and polygons. Vector models are useful for storing data that > has discrete boundaries, such as country borders, land parcels, and streets. > > _**raster data model_** : [data models] A representation of the world as a > surface divided into a regular grid of cells. Raster models are useful for > storing data that varies continuously, as in an aerial photograph, a > satellite image, a surface of chemical concentrations, or an elevation > surface. All I have understood from the above is that both vector and raster data constitute of "latitudes and longitudes", only. The difference is in the way they are displayed. Latitudes and Longitudes in Vector data are displayed in the form of lines, points, etc. Latitudes and Longitudes in Raster data are displayed in the form of closed shapes where each pixel has a particular latitude and longitude associated with it. Is my understanding correct?


title: Convert a shapefile with northing/easting to ASCII raster with latitude/longitude

text: Linked below is a polygon shapefile whose projection is NAD83 and whose values are LU/LC codes. Coordinates are meters of northing and easting. I would like to convert it to an ASCII raster with coordinates measured in latitude and longitude. I'm sure it is a simple process, but I can't figure out the re- projection. Tools I have are ArcGIS or SAGA or R. I can even open the ASCII text and modify the header if no recalculation is necessary. http://hawaii.gov/dbedt/gis/lulc.htm EDIT: Thanks to Dan Patterson for pointing this out: the projection is actually UTM 4, and NAD83 is the DATUM. The output raster should have the same datum. If I want output to have lat/long coordinates, does that mean it has no projection?


title: QGIS- How to create an "area of interest" polygon mask

text: I've created a GP tool using ArcGIS that creates an "area of interest" mask from selected features. Effectively it creates a new donut polygon layer where the hole was the previously selected features. Can I replicate this task in QGIS? Example: For a specific map, I want to highlight a county within the state and mask the rest of the state and features outside the county with a gray, semi- transparent polygon layer above. Conceptually I would do it similar to this.   1. User loads state, county, and river layer   2. User selects county of interest   3. Run a utility (if it exists) that exports all of the unselected features to a donut shapefile without the originally selected polygon. If this layer mask was an "in memory" layer only available to QGIS, that would be just fine, too, maybe even better. Does this type of tool or capability already exists within QGIS or will I have to write this as a custom function? Here is how I did it in ArcGIS with ModelBuilder and the result I would be looking for, where the yellow layer and river layer exists statewide (all the way across the image), and the gray layer is transparent with a "donut hole", allowing the bright yellow portion to display through and draw the viewers attention the the area of interest: ![enter image description here](http://i.stack.imgur.com/vXNn8.png) ![enter image description here](http://i.stack.imgur.com/HADvK.jpg)


title: WMS projection doesn't match basemap

text: I have a flex viewer 3.0 map in wgs 84 (4326) (ESRI basemap) with a wms which is 4326 as well. for some reason the wms doesn't overlay correctly. it is way off. appears on the maps but it is flattened quite a bit.


title: How to represent overlapping lines in ArcGIS 10.1?

text: I have lots of lines (in one feature class) which share the same location but I want to represent them as lines that are next to each other because they have different properties (different OWNERS). So, my first try was transforming them in "Cartographic Line Symbol' -> 'Line Properties' set my offset different for each one. Which is fine except: ![enter image description here](http://i.stack.imgur.com/mSkJH.jpg) As you can see when I choose to represent the lines according to their owner (color= owner). The offset works great, but I need all the lines to be connected and sink in a junction. Something like this: ![enter image description here](http://i.stack.imgur.com/Xt1hy.jpg) So I tought the solution is to create a representation, but not sure how to use a text field (OWNERS) in order to create those parallel lines (the equivalent of Symbology->Categories-> Unique Values-> Value field: OWNER). I know I should use the offset geometric effect. but I'm not sure if this is the solution and how should I do that. Thanks for help!


title: What is the Market share for open source web-GIS systems?

text: I would like to know the market share of 3 web-GIS systems, deegree, GeoServer, Geomajas). Which one is more popular? Are there any evaluations regarding the market share of open source web-GIS systems?


title: How to integrated GIS functionality into an existing web-based business process management software?

text: in a project I am doing, there is a need to combine data from different sources into a GIS and show it in one map. The data could be from ERP system like SAP or other sources. which tool could be used for combining different data sources in GIS? which software would be the best? I am thinking about Geoserver or Geomajas? Thank you for any response


title: Counting number of points in a point shapefile - ArcGIS 10.0 / arcpy

text: I assumed this would be an easy thing to do. I'm loading a point shapefile into a map like this:               layer = arcpy.mapping.Layer(file_name)      Is there a way to get the number of points in the shapefile? I see using the 'Describe' command that I can get some info about things... For instance,               arcpy.Describe(layer).table.fields      Gives me some info about the fields in the underlying table but I can't count the rows.


title: How to only show labels for an arbitrary selection of items?

text: I'm curious how others go about solving this problem: You've created a map for something with a large number of features that are labeled. The customer/client asks that you only show labels for X, Y and Z, based on some seemingly arbitrary decision (eg. what they deem important features). How would you go about doing this? Some ideas:   * Create a new string column for this special label and only fill in a value for the features they want to see (could result in duplicate information)   * Create a new boolean column and flag the features they want to see with true, then use the conditional labeling in QGIS 1.8 to only display the label when the boolean is true


title: Adding a new Coveragestore in a certain Workspace by REST

text: So, simple thing as it seems but I can't seem to make it work! I'm trying to add a raster data file (GeoTIFF) to geoserver via REST API, using curl on linux. REST API _is_ working, since I created a new workspace using it:               curl -u admin:geoserver -v -XPOST -H 'Content-type: text/xml' \      -d '<workspace><name>restProba</name></workspace>' \      http://localhost:8080/geoserver/rest/workspaces      It's correctly created as shown by Geoserver's admin ui, but when I try to add a coveragestore to that workspace by executing               curl -u admin:geoserver -v -XPOST -H 'Content-type: text/xml' -d     '<coverageStore><name>int_dec</name><enabled>true</enabled>     <type>GeoTIFF</type> <url>$home/int_dec.tif</url></coverageStore>'          "http://localhost:8080/geoserver/rest/workspaces/restProba/coveragestores?configure=all"      even though it is correctly loaded, it goes directly into the default workspace, which is not the one I want it in. Is there any way of defining which workspace i want it in? I thought it'd enough by mentioning it in the URL where I point to restProba's workspace URI but it seems to be omiting it. Thanks for your help :-) pd: newbie with geoserver, I've tried searching for this but either i'm not using the correct search criteria, im a fool or it's just too simple for anyone to ask that i'm even fooler for asking it :p


title: ArcGIS/ArcPy Python, how to get latitude and longitude of a projected point?

text: I have a point feature in a feature class that is being accessed by ArcPy. The point is projected but I need to find an efficient means to get the unprojected latitude and longitude for this point. Is there a method other than reprojecting (unprojecting), getting a search cursor on the new feature class, finding the feature, then getting the lat/lon off the feature's shape?


title: Statewide (and more) Address Point data source

text: I am seeking a data source for a statewide address point layer. I do not need resident or business information at these addresses; I just need the validated structured addresses (preferably with standardized elements) and the coordinate. This is for an address suggestion application on clients without internet access. We could do server side queries to an internet service, but not client side. (So, for example, this is why we cannot use Google for queries.) We will be querying _partial_ addresses to build suggestions, so generally geocode services will not work. Partial addresses will have a zip code. Ideally though, we would obtain a database of address points, allowing a partioned table query server side against that database. The main State we want is Missouri, followed by Illinois, but could go as high as a nationwide database. We have interest at this time in countries other than the United States. Do not worry about budget issues on this. I need to gather information on available sources before working up a budget/RFP.


title: List of WMS Servers

text: > **Possible Duplicate:**   >  GIS Map Service Software I'm experimenting with some applications and I would like to know what are our options for wms servers, such as MapServer and Geoserver. Do you guys know any other implementations? OS preferred, but private is ok too. Can we create a comprehensive list of WMS server software?


title: Are there any known problems using QGIS on Mac OSX Mountain Lion?

text: A few days ago Apple released a new version of the OSX (10.8) and I'm thinking about upgrading. So, I was wondering if anyone had any problems with GIS applications (QGIS) and other similar software after the upgrade? Thanks,


title: Heads up digitizing - Polygons overlay differently with two different Satellite images

text: I have been digitizing building footprints from satellite image using ArcGIS 10.1 . Field data collections is done with the use of four corner point GPS way point. Due to multipath errors, we are using satellite image to complement it for better accuracy. I have two different satellite images which one of them is a recent one but lower resolution, while the other one is an old one with a better resolution. I have tried to use the old one with better resolution in my heads up digitizing. It went well . My problem is once overlayed with recent image, polygons bounce off few meters to the road. **Polygon Coordinates and Projection:** Projected Coordinate System: WGS_84_UTM_zone_38N Projection: Transverse_Mercator false_easting: 500000.00000000 false_northing: 0.00000000 central_meridian: 45.00000000 scale_factor: 0.99960000 latitude_of_origin: 0.00000000 Linear Unit: Meter Geographic Coordinate System: GCS_WGS_1984 Datum: D_WGS_1984 Prime Meridian: Greenwich Angular Unit: Degree **Old Satelite image Coordinates** Spatial Reference: GCS_WGS_1984 Angular Unit: Degree (0.0174532925199433) Datum: D_WGS_1984 **New Satellite image Coordinates** Spatial Reference: WGS_84_UTM_zone_38N Angular Unit:Degree (0.0174532925199433) Datum: D_WGS_1984 What could be the problem and how would you go around it? Thanks


title: Thiessen polygon input help

text: I am new to 10.0 and have never used the thiessen polygon tool -- i have a set of points and would like to create a thiessen polygon; what do i need to input (i.e. what type of input (.shp etc)?


title: relationship between sp package in R; mapwindow; ArcGIS/Quantum GIS

text: I am using ArcGIS at office for some visualization work, but I feel it kinda dump to performing the same tasks several times for visualizing some data on different maps. I am think of automating it using R. I have come across several things, `sp` package for R, ArcGIS for visualization, MapWindow for drawing maps, and another opensource Quantum GIS. But I puzzled on the relationship between. For example, if I use MapWindow to draw the maps, can I use only R for visualization work? And what's difference between R+`sp` package and Quantum GIS + python + R? I am definitely a newer here, and thanks for the help. And a little bit more details here: I am trying to visualize the location of beds in an elderly home which the elderly is infected with some bacteria, the bed will be shown red if the patient gets the bug, and shown green if the patient does not. It is perfectly OK for me to use only ArcGIS if I only need to visualize the clustering of bacterial infection in one elderly home. But the problem is: I have 33 of them. And I don't like repetitive tasks to be done manually. So can I use R plus some kind of GIS system to do the repetitive task, and give me the 33 charts?


title: Is it possible to omit the headings of the CSV output of a GetFeature request?

text: I am making a web service which will consume the result of a WFS GetFeature request made against Geoserver. It would tidy up the code on my side if I was able to specify that the CSV output from Geoserver should not include column headers. Is this possible? To me, this page suggests that it probably isn't possible, but I thought I would check.


title: How to find nearest Lat/Lon to the given Lat/Lon from set of Lat/Lon values?

text: I have a set of lat/lon values and gps lat/lon value. I want to find nearest lat/lon value to the gps point from that lat/lon set. Is there any some formula or some other way so that i could implement that in my c++ program.


title: Ecw not supported with qgis-mapserv-fcgi but supported in QGIS Desktop - resolved

text: I have problem with qgis-mapserver version 2 and ECW. Qgis Desktop works well with ECW (gdal support for ecw with the library ERDAS_ECWJP2_SDK-5.0) but not qgis-mapserver (not supported file in Apache error logs). qgis-server must use the same gdal librarie that qgis desktop (QGIS desktop and server are all 2 installed on my computer (ubuntu 12.10). This is what I do not understand. Thank you Note: > I just watched the license to use the ERDAS_ECWJP2_SDK-5.0, including this > paragraph: 1.9. "Server Application" means clustering software products That > You Develop or maintain maintenance and are able of running on a Server To > Provide an image decoding or processing capability. This License Agreement > Does not Authorize You to use the SDK to create Server Applications. > > With this SDK for a "free" installation, or it is installed with a "Desktop > Read-Only" license or with a license "Server Read-Only End User" but not the > 2. > > Next time, I would read longer uses licenses, something I rarely do. Thank you for your answers


title: Building a Polygon from an array of points in Arcobjects and C#

text: I have seen this code in Arcpy which is used to create a Polygon from a string of Points(X , Y).Can you please let me know how I can us something like this which can generate polygon from an array of X,y values in C# and Arcobjects? I already take alook at Arcobjects 10 .net framework sample code sample of Building a Polygon Using Points . but it was too confusing for me. I really appreciate your time if you can point me to a tutorial or a sample which contains a process of polygon creation same as what was in arcpy code.


title: How to run VB scripts in ArcGIS?

text: Some questions about VB scipts:   1. Can I read and modify the code of a VB script? If the answer is yes, how can I do it?   2. How can I run a VB script?   3. Can I use it as a VBA code eg.using variables instead of certain data as defined in a model created using the Model Builder? Thanks


title: How to organize a kml created from a shapefile in QGIS?

text: I have a shapefile that I would like to export into a kml to share it with my colleagues. This shapefile contains a few hundred polygons. I have QGis 2.2.0 on Windows 7. I would like to create a kml that somehow keeps the same style as my shapefile. Indeed, in Qgis, I have a simple category style applied to a certain field ("Soil Type"). When I export the shapefile as a kml (after renaming the "Soil Type" field as "Name" to have it as a title of my kml polygons), the kml created doesn't have any sub-folders per Soil Type. As a result, if I want to set a color per Soil Type, I need to manually create folders in Google Earth, move every single polygons into the different new folders, then apply a style to the folders. It works, but it takes a long time when you have a lot of polygons. Another solution I found is to create a different kml of each "Soil Type" in QGis. But again, it's not very efficient. Is there any solution more efficient to organize the kml created from a shapefile in QGIS? Thank you


title: Why and How does copying and pasting a "problem" model with errors in it into a new empty model fix the model (ArcGIS, ModelBuilder)

text: I often work with ModelBuilder and Models that have previously worked seem to become unstable or develop issues after lots of editing. They return error messages for parts of the model that did not previously return errors. A solution often recommended is to simply copy and paste the contents of the model into a new empty model. Does anyone have expertise or advice to share that can explain why this fixes models and if this is likely to be a permanent fix? I am developing a set of models that I will be sharing and I need them to be robust enough that they work for the model users. One of the reasons I am asking this is that often when contacting ESRI support they will focus entirely on the reported error messages within ModelBuilder even when you tell them the Model has previously worked OK, or that it runs OK on a different computer or from the Edit vs open window. It seems to me that very often the error messages returned by models are of no use at all, especially when the problem for example is fixed by copying and pasting into a new empty model. It could just be that these error messages are a symptom of a model that is becoming unstable? Any thoughts or advice would be very gratefully received.


title: Is there a OSM data server that returns vector data based on a bounding box?

text: I would like to know if exist some OSM data server, not tiled, just data. I know it is possible to import OSM data in your own server, but because you need robust server architecture, it wont be possible for now. I am looking for a server that returns OSM data using a bounding box. A server that can be used with Mapnik OSM Plugin.


title: How to select polygons by distance to a single polygon in another layer?

text: I have a shapefile A with a group of polygons, lakes.   Now I want to select the polygons(lakes) that are within a specific distance to a polygon in a second shapefile, B, with the polygon a city. In other words, I am looking for something like ArcGIS "Select By Location" and "are within a distance of the source layer feature"   I tried Select by Location, and also Select by Radius, but didn't manage to get any of them to work.


title: How to get the file paths of the loaded layers in a QGIS plugin?

text: I am very new to GIS, QGIS and Python so please bear with me. I am using QGIS 2.0 and I would like to develop a QGIS plugin that would display the raster layers of the main window in order to allow for some specific interaction. I have a very hard time getting this to work! I generated the skeleton of my plugin through plugin builder, and changed the [moduleName]2dialog.py file to try to display the main window layers. So I have two questions : 1/ Is it possible to get the layers? 2/ Can I draw it inside a GqsMapCanvas inside my plugin GUI, or do I need another object? This is the code I have for now (generates an error):               from PyQt4 import QtCore, QtGui     from qgis import core, gui, utils     from qgis.gui import QgsMapCanvasLayer          mc = iface.mapCanvas()     layers = mc.mapRenderer().layerSet()     if not layers:         self.ui.displayFileNameTextEdit.setText("No layer to display.")     else:         for layer in layers: #layer is one id             self.ui.qgsMapCanvas.setExtent(layer.extent())             self.ui.qgsMapCanvas.setLayerSet( [ QgsMapCanvasLayer(layer) ] )             self.ui.qgsMapCanvas.refresh()      Thank you very much! I'll appreciate any idea.


title: How to install QGIS 2.0 on Ubuntu 13.10?

text: How can I install new version of QGIS (2.0)? When I try with the indication of the official web site of qgis, terminal write me that the last version of qgis (1.7.5 wrowclaw) has several plugins that not permise it. I am a new user of ubuntu 13.10. Thank you very much,


title: Hard to get my head around, generating average price for houses in each quadrant

text: I'm able to determine the average price of houses sold in a particular area (like the squares on image below) dependent upon postcode. Although I can generate the average sale price for a particular postcode, I'm having difficulty utilizing these to the generations to the average house price for each quadrant. As in image. Assuming I know the average house sale price for each of the letters noted on the image below, **how would you generate the average price for each colored quadrant?** ![enter image description here](http://i.stack.imgur.com/iKuhh.png) **Further information:**   1. The information contained within each quadrant (red, orange, green) should only impact that quadrant. (It's these values I hope to utilize to form an average for each separate quadrant)   2. These quadrants are intended to be eventually presented within Google Maps using circles (which I know how to do)


title: Get all points with in some range of a point (or range of point)

text: I have got a list of points and a fixed range (10m). So what I have to do is for each point I have to find all the points present in a table with the 10m range of the given point. For example: Let say I have following list of points. P1, P2, P3, P4, P5, P6 Now I have to loop through each point. First I'll get P1 and I need all points (in a table) present inside the circle where radius is 10m and center is P1. Similarly I'll do for all other points from P2 to P6. My result will be all the outcome points in each loop. I am not sure whether I have explained well or not. But I can clarify if you can comment. I also need to better solution, where I want to calculate the outcome points at time by creating a polygon may be and then somehow considering the range. Thanks


title: Is there a way to get the extent of a query and not each individual feature?

text: I'm using this solution to get the extent of a query. It brings every feature geometry back to the client where I then calculate the extent of all the individual geometries. Since these can be polyline and polygon features this is a very expensive process. Is there a way to form the query so only the extent is returned?                   app.filter.on("zoom-to-features", () => {         console.log("set map extent to include all annotations matching this criteria");         all(common.GeoUtil             .getFeatureLayers(app.map)             .map(l => {                 var query = new EsriQuery();                 query.returnGeometry = true;                 query.outFields = [];                 query.where = where;                 return l.queryFeatures(query);             })).then(results => {                 var features = <Array<{ features: Array<esri.Graphic> }>><any>results;                 var extents = features.map(result =>                     result.features                         .map(f => common.GeoUtil.getExtent(f.geometry))                         .reduce((a, b) => a ? a.union(b) : b, null));                      var extent = extents.reduce((a, b) => a ? b ? a.union(b) : a : b, null);                 extent && app.map.setExtent(extent);             });          });


title: Is it possible to add a timeline using cartodb?

text: I am trying to add a timeline in one of my visualization and I don't really know how to or even if you can add one in CartoDB. Does anyone know if you can and how?


Is there an average size for British postcodes? Or a method of determining the size for a particular British postcode? Thanks. Considering creating a related program and this would be valuable information. I'm referring to postcodes like:   1. BT28 1AW   2. LS9 8AH   3. L15 3DE By size, I'm referring to the amount of area they cover (1 miles - 20miles)


title: How can I correct a .pgw file for a PNG image?

text: I have a PNG file and I need its .pgw file to be set with NAD 83 as the projection. I have created a .pgw file: 799   0   0   -799   -122.562   37.1946 However when I open the PNG in QGIS or ArcMap 10.2 it is not being placed in the right location. I have looked around on how to form a .pgw file and I think I am doing it right, but something is obviously missing. My PNG file has an associated KML which is where I got the coordinates for my .pgw file. My PNG and.pwg files have the same file name as well.


title: How to find nearest city using Geonames features codes

text: I've got about 300,000 lat/lng coordinates of events that I'm combining into a large map. The coordinates come from all over the world. For each point, I want to find the nearest city of at least XYZ residents. I've written a Node script that uses the geoNames database of world cities and does all the calculations. It all works fine but for the fact that there seems to be a far amount of overlap in some of these cities. For example, `cities15000.txt` lists both "North Atlanta" and "Atlanta," and there is both a "Manhattan" and a "New York City." I was hoping I could limit the search to the highest level in the city hierarchy using geoNames' feature codes, but the documentation is not terribly clear. North Atlanta is "PPL" (populated place), Atlanta is PPLA (seat of a first-order administrative division). I thought I could limit to PPLA, but Baltimore is PPLA2 (seat of a second-order administrative division). There does not appear to be anything near Baltimore (say, "City of Baltimore" or something) that I'm missing. Where do these codes come from? Is there a way to eliminate subsections of cities? I also don't mind using a different source, though I do find the geoNames data clean and easy to manipulate.


title: how to load a raster into the Map Canvas using QGIS Action

text: Is it possible to write an Action in QGIS that will load a raster into the current map canvas based on a url/path of the raster as specificed in the attribute table of a shape file? The shape-file is a index-tile.


title: Plotting pixels on a map image

text: I've got a list of some hundreds of latitude/longitude pairs close to each other (within a country). I want to create an image with a map, and circles overlaid around these points. Basically, what I need is:   * an image of a map   * a way of translating latitude/longitude pairs to pixels Where can I obtain a fitting map (e.g. roadmap or satellite map) with the parameters needed to plot the points? I suppose there is some FOSS tool that does most of this already?


title: Problem connecting QGIS to MySQL

text: I'm trying to create a new vector layer and store it in a MySql-database, but when I test the connection I get an error with no detailed information at all. If I give an erroneus password I get a message that says: > MySQL connect failed for: > db_qgisgeodata,host=localhost,port=3306,user=xxxx,password=ppppppp Access > denied for user 'uuuu'@'localhost' (using password: YES) Because I get a detailed error message if i give the wrong login credentials I think QGIS is able to connect to the database, but there is something more to it. The user has full permissions on the database. The problem is the same if I try with the root-user.   1. The user has full permissions on the database. The problem is the same if I try with the root-user.    2. The database is empty.   3. Windows 7x64   4. QGIS 1.8.0   5. MySql 5.5.28   6. it's no problem to connect with the database from other applications. I guess I'm missing something quite obvious here, but I can't figure out what. Any ideas anyone?


title: OpenLayers rules and context, what is the difference?

text: OpenLayers.Style employs two properties, "rule" and "context" which accordingly contain an array of OpenLayers.Rule and a set of methods. I have noticed that "context" can somehow substitute "rule". For example, in a case where there is a vector layer displaying temperatures as features, with "temp" in its properties for the values.               var temperatureStyle1 = new OpenLayers.Style({         //common symbolizer styles         fontFamily: "Arial",         label: "${temp}",          fontColor: "${coltem}",         fontWeight: "bold"         },         {         context: {             coltem: function(feature) {                  var temperature=parseInt(feature.attributes.temp);                 var col="#000000";                 if(temperature>35) { col="#CF2900"; } //Hot                 if(temperature<0) { col="#0000FF"; } //Cold                 return col;              }         }     });      is somehow equivalent to                   var temperatureStyle2 = new OpenLayers.Style({             //common symbolizer styles             fontFamily: "Arial",             label: "${temp}",              //fontColor: "${coltem}", not really needed             fontWeight: "bold"             },             {             rules: [                 new OpenLayers.Rule({                     filter: new OpenLayers.Filter.Comparison({                         type: new OpenLayers.Filter.Comparison.LESS_THAN,                         property: "temp",                         value: 0                         }),                     symbolizer: { fontColor: "#0000FF" }                     }),                 new OpenLayers.Rule({                     filter: new OpenLayers.Filter.Comparison({                         type: new OpenLayers.Filter.Comparison.GREATER_THAN,                         property: "temp",                         value: 35                         }),                     symbolizer: { fontColor: "#CF2900" }                     }),                 new OpenLayers.Rule({                     elseFilter: true,                     symbolizer: { fontColor: "#000000" }                     })                 ]             });      What are the benefits in each of these two approaches? I can only see that with "context" I can change only one style (of the symbolizer) at a time, while with "rules" I can redefine a whole symbolizer -- instead of making many context methods for each style doing the same comparison every time. Is there anything more that I cannot see? What's the reason for the existence of the "context" object?


title: QGIS 1.8 Did Not Install Properly-- Will not recognize modules or DLLs-- conflicting with existing PyQt?

text: I know this has been asked and re-asked, but I can't seem to get any of the proposed solutions to work. Upon clicking the start menu shortcuts to QGIS desktop, nothing happens. When I run from command prompt, I first received the error that "qwt5.dll" was missing. This was installed and present in the Quantum GIS Lisboa/bin folder, and I copied it into the QGIS app/bin folder, along with the next series of dll's it said were missing. It then told me it could not find PyQt4, and then that the SIP module was missing and that python would be disabled. I ensured that the paths were included in my environment variables, but no change. I was curious if the existing versions of PyQt4 (and the accompanying SIP module) that were already installed on my computer could be the cause of this? This is on my work machine which has heavy security measures in place. On my home machine, which also had PyQt4 and SIP installed previously, I had no issues with installation or launch. Both machines are 32-bit Windows 7. Does anyone have any ideas? EDIT: After much, much screwing around with paths and bat files, I got the desktop, browser, GRASS and MSYS to launch. I still cannot launch OSGeo4W.... when I run the bat file in command prompt, I first got "findstr not recognized".... copied, pasted, re-ran... got "_path=C:/Apps/QGIS//apps/Python27/scripts find: .exe no such file or directory, find: .bat no such file or directory "........hmm..


title: Python GDAL package missing header file when installing via pip

text: I asked this question on Stack Overflow before I realised this site existed. As this site it more appropriate, I'll ask it here too. I'm trying to install gdal from pip `pip install gdal` inside a virtual environment (Ubuntu). It fails because it cannot find `cpl_port.h`               extensions/gdal_wrap.cpp:2853:22: fatal error: cpl_port.h: No such file or directory     compilation terminated      However GDAL is installed correctly and the header file is located at `/usr/include/gdal/cpl_port.h`. Is there some environment variable for GDAL that needs to be set in order for pip to find the header files?


title: Why is export of joined table half empty (all zeroes) even though joined table was fine?

text: I have joined two tables (excel 10 but saved as 03 to import it into ArcGIS) in ArcGIS 10 based on a column. This column has some empty records in both tables. I need the agreeing records to be added to the right and the others to be on the bottom of the table. After fiddling with it I got it to work; we filled in the empty cells of that column with dummy numbers. Anyway, the problem now is that when I export the joined table, the second (right) half of it is empty (all zeroes) even though the joined table was just fine. I tried exporting it as dbf and txt and neither worked. Does anybody have an idea how to fix this? Thanks a bunch!


title: Automating an ArcGIS Project that is a display of map

text: I'm relatively new to GIS and looking to automate a project. I have my 'base map' created in ArcGIS, but now need to create several hundred more so looking for an automated method to accomplish this. I have a desktop version of ArcGIS 9.3. I'm currently using VB, Python isn't installed and will require going through IT and business case... Things I need to figure out: 1. Select the Area of Interest (Polygon) 2. Zoom to the select area 3. Change title 4. Change legend to add some details (can be either the layer name or the description for the layer). 5. Export to PDF. I have figured out #1, but stuck on #2, because although it does zoom to the polygon of interest it isn't choosing the same zoom I'd get if I manually clicked "Zoom to Selected Features". Any suggestions on things I should look into or resources? Thanks, Reese


title: Making a Button in mapinfo with conditional function

text: I have following code :               Sub create_ButtonFromVendor           Dim i As Integer           i=1           Create ButtonPad "allVendor"  As             ToggleButton             Icon MI_ICON_ZOOM_QUESTION               Calling get_Vendor               HelpMsg vendors(i)             Show            For i = 2 TO Ubound(vendors)               Alter ButtonPad "allVendor" Add               ToggleButton               Icon MI_ICON_ZOOM_QUESTION               Calling vendor_ButtonAction(vendors(i))               HelpMsg vendors(i)             Show           Next         End Sub      It is small sub procedure. Here I wanted to create multiple buttons with different function in a loop. How would I add function in a loop? In the code I added above does not work as handler must be sub procedure without parameter.


title: How to draw points in MapBasic?

text: how can we draw a point object in mapinfo using mapbasic, my objective is to generate an array and then map it. Thanks


title: Looking for software to interface GPS data with PC

text: I have GPS Unit from Garmin (GPS 18xPC). I'm looking for software I can install on my PC that will allow me to interface with the data I collect on my GPS unit. Any help will be great


title: arcgis 10.1 account manager, authorizeSoftware script

text: I am trying to run "authorizeSoftware" file to complete creating new account manager in Arcgis 10.1, but after running               ./authorizeSoftware      I got this message: > -------------------------------------------------------------------------- > Starting the ArcGIS Software Authorization Wizard > > Run this script with --help for additional information. > -------------------------------------------------------------------------- > Unable to access X server (:0.0). Please check your DISPLAY variable or run > this script in silent mode. Run with the --help option for more information Any idea how I can fix this?


title: Download Google Map tiles, images

text: I would like to use the google maps downloader to save satellite images from google maps (http://www.allallsoft.com/gmd/index.html), but is it legal to use? In the terms of service you can read that: "Unless you have received prior written authorization from Google (or, as applicable, from the provider of particular Content), you must not: (a) copy, translate, modify, or make derivative works of the Content or any part thereof; (b) redistribute, sublicense, rent, publish, sell, assign, lease, market, transfer, or otherwise make the Products or Content available to third parties...(d) use the Products in a manner that gives you or any other person access to **mass downloads or bulk feeds of any Content** , including but not limited to numerical latitude or longitude coordinates, imagery, and visible map data..." This program downloads the tiles and match them. The internet explorer does the same. When I view google maps in internet explorer, it downloads the tiles to "C:\Users\username\AppData\Local\Microsoft\Windows\Temporary Internet Files" and put them together. So it copies imagery to folder and mass downloads them.


title: How can I change the precision in the attribute table?

text: I'm working with QGIS, v. 2.0.1. I opened a csv-table (from an Excel-file) and now the precision is very high but I don't need this. At the moment the values are for example 15.00000000000000, how can I change this to 2 decimals? And is it possible to set the precision before I open the table in QGIS? The methods suggested by simo just worked partially... I opened the table and it looked good but as I wanted to create a shapefile of the layer the precision switched back to many decimals. What can I do? There are only two decimal places in the csv file and there are also only two decimal places when I load the csv file into QGIS. But as soon as I create a shapefile, there are 14 decimal places.


title: Raw State Data on Health and Social Outcomes

text: I am looking for raw health and social outcome data (poverty, obesity, education levels, etc.) by state. I would like to download it as a csv or excel file. It needs to have the state fips listed in one of the columns. There is a lot of county data out there, but can anyone point me in the right directions for state data?


title: Nearest road in the vicinity of a point without queries (openstreetmap)

text: I know that many similar questions about map fitting have been asked. Everywhere there were answers about sql queries like this:               SELECT ST_AsText (ST_ClosestPoint (pt, line)) AS cp_pt_line,      ST_AsText (ST_ClosestPoint (line, pt)) As cp_line_pt      FROM (SELECT 'POINT (100100)' :: geometry As pt,      'LINESTRING (20 80 98 190 110 180 50 75)' :: geometry As line      ) As foo;      I'm looking fot any library or map fitting methods without using queries. Preferably in java or python. Maybe some has a little example like reading data about nearest roads (I don't need code for map fitting algorithm, I can implement it myself) Thanks for any answer, help or link.


title: Casting a Geography type in PostGIS

text: I have the following table defined, which uses PostGIS geography.               CREATE TABLE test_geog (         id SERIAL PRIMARY KEY,         boundary GEOGRAPHY(Polygon)     );      I added the following test polygon to the table:               INSERT INTO test_geog VALUES (ST_GeographyFromText('SRID=4326;POLYGON((0 0, 0 5, 5 5, 5 0, 0 0))'));      I am trying to determine whether a point lies within any of the polygons in this table. I have this query:               SELECT ST_ContainsProperly(ST_GeographyFromText('Point(2 2)'), area)     FROM (SELECT boundary FROM test_geog) AS area;      This yields the following error:               ERROR:  function st_containsproperly(geography, record) does not exist     LINE 1: SELECT ST_ContainsProperly(ST_GeographyFromText('Point(2 2)'...                 ^     HINT:  No function matches the given name and argument types. You might need to add explicit type casts.      How do I convert this "record" into a `POLYGON`? I'm confused because it seems like the column has already been declared to only hold `POLYGON` types, but for some reason that is not what I am pulling out of the database. I attempted to cast the record to a `POLYGON` like this:               SELECT ST_ContainsProperly(ST_GeographyFromText('Point(2 2)'), CAST (nitf_area AS POLYGON))     FROM (SELECT boundary FROM source_imagery) AS nitf_area;      But that gives me this error:               ERROR:  cannot cast type record to polygon     LINE 1: ...tainsProperly(ST_GeographyFromText('Point(2 2)'), CAST (nitf...      What am I not understanding here?


title: Toggle Identify with Button

text:
I am attempting to set up an Identify button that will only allow the identify functionality to be performed when the button is pressed. I have gotten to the point that the identify only works when the button has been pushed the first time but I can't seem to turn the button off. I have pasted a few bits of code below or you can follow the link to see the entire thing. http://geoville.org/viewers/PopUp/ **Funtion to activate Identify**               function activateIdentify(){         if (dijit.byId("tool_identify").checked) {             dojo.connect(map, "onClick", executeIdentifyTask);         }       }      **Button to Activate function**
Identify
      **Function to toggle buttons**               function toggleButtonIcon(tool) {              //only the tools in the toolbar are dijit togglebuttons so can iterate thru them         dijit.registry.byClass("dijit.form.ToggleButton").forEach(function(togbtn) {           if (togbtn == tool) {             togbtn.attr("checked", true);           }           else {             togbtn.attr("checked", false);           }         });       }  //end toggleButtonIcon      Thanks for your help. -Mike


title: How to join attributes from nearest point in QGIS?

text: I have two layers with points, call them A and B, each with different attributes. For each point in A I want to join all the attributes from the nearest point in B. Ideally this would also include a field measuring the distance between both points. I know how to so this in ArcGIS: right click layer A>joins and relates>joins, then join data based on spatial location and choose to join attributes from nearest point. I could not find similar functionality in QGIS. Any pointers?


title: matplotlib does not work in ArcGIS 10 and Eclipse (Python)

text: I have installed matplotlib (v. 1.0.1) which works with iPython quite good. But if I insert the same script in ArcGIS 10 or Eclipse 3.7.0 (with PyDev) I get always errors. As I consider to use matplotlib to create graphs with ArcGIS Python Scripts developed in Eclipse thats both pretty unfortunate - so far eclipse works quite well for ArcGIS 10. F.i. I tested this one: http://matplotlib.sourceforge.net/examples/api/unicode_minus.html In iPython it works perfect. In Eclipse I get this error: import matplotlib.pyplot as plt ImportError: No module named pyplot And in ArcGIS I get this one: Running script Script1... : bad screen distance "640.0" Failed to execute (Script1). Do I still have something to configure? Best Werner