title: how to rotate PIXMAP based symbols correctly in Mapserver 5.0.4?

text: I'm on mapserver 5.0.4 I want to draw a polyline with a symbol on every 100th pixel. The symbol is a PNG image with size 10 x 7 pixels. I want the symbol to be drawn perpendicular against the polyline The class in my mapfile drawing the pattern is this:                   CLASS             NAME 'warmfront'             EXPRESSION ([type] == 2 )  #warmtefront             STYLE                 SYMBOL 'warm'             END             STYLE                 SYMBOL 'circle' #draw a connecting line                 SIZE 2                 MINSIZE 1                 MAXSIZE 5                 COLOR 255 0 0             END          END  # CLASS      The symbol "warm" is defined as:               SYMBOL         NAME "warm"         TYPE PIXMAP         IMAGE 'symbols/warm.png'         GAP -10000000000     END      The image looks now like this: ![screenshot](http://i.stack.imgur.com/cVg0Z.png) The red triangles are now indeed drawn along the line, but the changing the GAP value doesn't seem to have any effect. What am I missing? EDIT: as suggested in the comments below, I changed the smbol from a picture to a vector:               SYMBOL     NAME 'warm'     TYPE VECTOR     FILLED TRUE     POINTS       0 4       2 0       4 4       0 4     END     GAP -10     END      The image now looks: ![enter image description here](http://i.stack.imgur.com/Ud6ny.png) No symbol is drawn, or the symbol is drawn invisible? Thank you ====================================================== EDIT: As suggested below I added a SIZE in my STYLE object:               CLASS         NAME 'warmfront'         EXPRESSION ([type] == 2 )  #warmtefront         STYLE             SYMBOL 'circle'             SIZE 40             MINSIZE 40             MAXSIZE 40             COLOR 255 0 0          END          STYLE             SYMBOL 'warm'             SIZE 30             MINSIZE 30             MAXSIZE 30             COLOR 0 255 0         END          END  # CLASS      and my SYMBOL def:               NAME 'warm'     TYPE VECTOR     FILLED TRUE     POINTS       0 4       2 0       4 4       0 4     END     GAP -100      The image looks now like this: ![updated screenshot](http://i.stack.imgur.com/z9jTK.png) But still the SYMBOL is drawn for every point on the line. When i use a pattern in my SYMBOL def like `PATTERN 1 60 END` the triangles do have a gap in between but they are not drawn along the direction of the line. What do i miss?


title: WMS in Google Chrome (ArcGIS API JavaScript) not working

text: I am doing a small web with WMS, but I have a problem with Google Chrome (and IE... but with this one I don't mind). Down of this lines there is two website: Galicia: This one contains WMS from Galicia (north of Spain) and it _loads correctly_ on Firefox and Chrome. Right http://signo-geo.com/alejandrowms/mapacorrecto2.html Andalucia (south of Spain): In this second App, the WMS on Chrome is not loading but in Firefox it works fine. http://signo-geo.com/alejandrowms/mapaprueba2.html I was looking for information about that and only find this thread.               http:// forums.arcgis.com/threads/98149-Chrome-WMS-and-content-type-application-vnd.ogc.wms_xml-fails..      Here is an Image which shows the problem: ![enter image description here](http://i.stack.imgur.com/dMctd.jpg) Why is the WMS not appearing in the second App, even though the code is similar to the first one? And why does this problem occur only in Google Chrome and not Mozilla Firefox?


title: How to get replaceDataSource to work with personal geodatabase?

text: I am having trouble using the arcpy.mapping's replaceDataSource with a personal geodatabase. The script runs, but does nothing if I run the script in the python window of ArcMap. If I run it in PyScripter, PyScripter stops working and gives the message: _EOFError: [Errno 10054] An existing connection was forcibly closed by the remote host_. It ought to be a simple fix; the feature classes are the exact same, just from a personal geodatabase that has been copied to a different location.               import arcpy          fcSource = r'c:\users\athom\Documents\GIS\GIS_computer_06_20_2013\GIS_LOCAL_MASTER\DATA\Vector_Data\Toponomy\TH_toponomy_250k.mdb'     fcName = 'sl_anno250k'     mapDoc = r'c:\users\athom\Documents\GIS\GIS_computer_06_20_2013\GIS_LOCAL_MASTER\MAPS\Map_Projects\TEST\AAAJ2011-004-02 TH Traditional Territory.mxd'     MxD = arcpy.mapping.MapDocument (mapDoc)     for lyr in arcpy.mapping.ListBrokenDataSources(MxD):         if lyr.supports ("DATASOURCE"):             if lyr.dataSource == r"N:\DATA\Vector_Data\Toponomy\TH_toponomy_250k.mdb\sl_anno250k":                 lyr.replaceDataSource(fcSource, 'ACCESS_WORKSPACE', fcName)                 arcpy.RefreshTOC()


title: How to transfer styles from QGIS to ArcGIS?

text: We have a QGIS project that has symbols for a lot of different lines. Data comes from a geodatabase that is accessible to both QGIS and ArcGIS. We now want to use the symbology of the QGIS project to style the features in the ArcGIS project. My first approach was to use SLDs. The QGIS style can be exported to SLD easily, but how can you style your features according to a SLD using ArcGIS? Setting up a simple WMS does not suffice, both QGIS and ArcGIS need editing capabilities on the features. If one could provide a SLD with a WFS (in a ESRI-readable way) that would also be a possibility. Have found nothing for that, either.


title: ST_Intersection point of two 3D linestrings - how is the Z value calculated?

text: I have two 3D linestrings (geometry type `LINESTRING Z`) in a PostGIS database that intersect in 2D space, but not in 3D space. Given that there is technically no 'shared portion' of geometry in 3D, how is the Z value returned by `ST_Intersection` calculated? Here is the query I am using:               SELECT feature, str1, (ST_Dump(geom)).geom.ST_Z as z_value     FROM     (SELECT ST_Intersection(linework.geom, align.geom) as geom, linework.feature, linework.str1     FROM align     INNER JOIN linework ON ST_Intersects(align.geom, linework.geom)     WHERE ST_isvalid(align.geom)='t' AND ST_isvalid(linework.geom)='t') q1;      What I would like is a query that returns the z value of the point of intersection but **on one of the linestrings only**. `linework` contains surveyed features that are crossed by an alignment in `align` - I would like to 'lift' elevations off the surveyed features at the points of intersection but I'm not sure how to go about it.


title: How to add colour relief raster legend to QGIS print composer?

text: I created a colour relief GeoTif from a DEM with the Relief function of Raster Terrain Analysis in QGIS 2.2. Now I´m trying to creat a map of that in the print composer and of course I would like to have a colour bar legend of the elevations. What I found so far is that it is possible to add raster legends in QGIS but only If I render the GeoTif as Singleband Pseudocolour with a reclassified colourmap. But that screws around with the colours and classification of my created colour relief. Any suggestions how to accomplish this in QGIS? Martin


I'm looking for some functionality beyond what ESRI already provides in the JSAPI. For example, the ability to display the LOD in a tooltip on the zoom control. It seems like there should be a dojo+jsapi answer to openlayers+geoext but I have not found it. Can anyone point me to Any Open Source Extensions to ESRI JSAPI?


title: 2 legend objects with arcgis javascript api

text: Can i create 2 legend objects in one script? I'm currently working on a web-application in which i already do have a legend/TOC running and it looks like this:               var legendLayers=[];         legendLayers.push({layer:raum, title:'Raumkategorien', slider:true});             map.on("layers-add-result", function(){             var layerInfo = new TOC({                     map: map,                     layerInfos: legendLayers,             }, "legendDiv");             layerInfo.startup();         });         map.addLayers([raum]);      Now I added some geoprocessing functionality and i wanted to add a new TOC when the result map is loaded in the web-app. I tried this by doing the following:               function gpJobComplete(jobinfo){             var mapurl = mapServiceURL + "/" + jobinfo.jobId;             var energy = new ArcGISDynamicMapServiceLayer(mapurl,{                 id:"Beste Energieformen",                 "opacity": 0.85             });             map.addLayers([energy]);                  var outputLayers = [];             outputLayers.push({layer:energy, title:"Beste Energieformen"});             map.on("layers-add-result", function(){                 domUtils.show(dom.byId('outputDiv'));                 var outputInfo = new TOC({                     map: map,                     layerInfos: outputLayers                 }, "outputDiv");             outputInfo.startup();             });         }      I get a result map and it works just nice in the web app, but there's no new TOC. Am i trying it the wrong way?


title: How to improve cartography by processing map exports with GIMP?

text: I am about to make a map with a heavy emphasis on cartography. I would really like to spice it up by processing the image in **GIMP** and I was wondering if anyone had any tips for bringing the cartographic pain to some ArcMap exports. My current plan is to export certain layers and elements individually and process them, but I will be clicking in the blind so I am looking for suggestions. What are some techniques for processing these exports with GIMP to produce a map?


title: How to select objects by N:M join?

text: I have a layer with some geometries and a data layer that defines an M:N relation. Is there any way I could select (on the map) objects that are in relation with specific id? Basically I would like to select by expression like               object_id in (select object_id from relations where relative_id = 12345)      I can do that from the db manager `select * from objects where object_id in`..., but that means creating new layer each time, which is extra work and needs to be deleted or there would soon be hundreds of them. And the qgis joins don't seem to support M:N relations. Is there a way to select (in qgis, not sql sense) by such query? The layers happen to be in spatialite database. They look something like               create table objects (object_id integer primary key, geometry, attributes...);     create table relations (object_id integer references objects (object_id),                             relative_id integer references something);


title: How to output parts of layers (shapefiles) that do not intersect using ModelBuilder?

text: I have two layers (shape files) that I want to select what do not intersect, to after all make a kernel of it. How do I do that in ModelBuilder without scripting? If someone can help me, even directing me to something to read, would be awesome. I'm using ArcGIS 10.0


title: file formats to represent georeferenced 3D data

text: I am looking for an overview of file formats that can store geo-referenced 3D geometries. There are no further requirements such as texture support, etc. I have found the following formats so far and would appreciate extensions to that list:   * Shapefile   * GML   * KML/KMZ   * GeoJSON


title: Delineate polygonal catchment areas from a raster using open-source GIS

text: I'm 90% sure this basic question has been asked before but looking through the other posts I simply do not see a useful answer. I need to calculate polygonal(vector) catchment areas from a raster based on minimum driving distance/cost from a set of points. I've already got the raster, and for that matter, the nodes(intersections with cost attribute) used to interpolate it. What I need now are polygons for each point describing the area nearer it than any other point. This is basically like a voronoi diagram but based on network distance rather than euclidean space. How exactly can I do this using only open-source tools? I'm comfortable with the command line if I just know what I'm looking for. I feel like GDAL is going to be the answer, but if so, what exactly is the command I'm looking for? Thanks!


title: determine the catchment area from dem or hillshade

text: I want to determine the catchment area of river from hillside. Is there a plugin that I can use in qgis ?


title: Postgis: ST_Buffer breaks because of no forward edges

text: I try to buffer a huge road network using ST_Buffer. After 8h of processing I get an error message: `ERROR: GEOSBuffer: TopologyException: No forward edges found in buffer subgraph` Here is my query (adapted from: memory issue when trying to buffer/union large dataset using postgis):               CREATE TABLE Nebenstrassen050 ("geom" geometry(multipolygon,3857), gid serial primary key);     CREATE SEQUENCE bseq;          WITH ordered AS (       SELECT wkb_geometry AS geom       FROM Nebenstrassen       ORDER BY ST_GeoHash(st_transform(wkb_geometry,4326))     ),     grouped AS (       SELECT nextval('bseq') / 100000 AS id, ST_Buffer(ST_Collect(ST_MakeValid(geom)),50) AS geom       FROM ordered       GROUP BY id     ),     groupedfinal AS (       SELECT (ST_Dump(ST_Union(geom))).geom AS geom       FROM grouped     )     INSERT INTO Nebenstrassen050 SELECT * FROM groupedfinal;      I tried to implement ST_MakeValid to avoid this error but it doesn’t show any effect. I’ working on a Win7 64bit machine with Postgres 9.3 and PostGIS 2.1. Any Ideas?


title: Common Web Server for Web Mapping?

text: I'm planning to create interactive routing map system( for small town only ) which will be running online. My server side language will be PHP and PostGIS/PostgreSQL with using pgRouting and client side script OpenLayer/JS. I'll be using XAMPP for development side and going to upload on my host( running Apache server) after the development. The base layer will be OSM. I read about Mapserver and GeoServer but this might not be fit for my hosting. I'm doubting to start since I'm scared for future failures like incompatibility and speed issue. If possible, can anyone give me a successful project using the same components I will use? Can anyone suggest for the perfection of this idea? Please let me know if this idea will work. Thank you in advance.


title: Select Layer by Attribute: selection using loops

text: I'm quite new with python and i'm currently stuck trying to make a loop... I have a shape file of a grid with several fields of different species (more than 150). The values of each record are 0 or 1 depending on the presence of the species on that grid. What i need to do is to make a selection of grids for each specie. So far this is my attempt:               import arcpy      import os, sys, datetime, time         arcpy.CheckOutExtension("Spatial")      from arcpy.sa import *      from arcpy    import env      arcpy.env.overwriteOutput = 1          arcpy.env.workspace = "G:\\Sara_hab_sp\\test"         arcpy.env.scratchWorkspace = "G:\\temp\\default.gdb"                   arcpy.MakeFeatureLayer_management("G:\\Sara_hab_sp\\Confocc_pergrid.shp", "G:\\Sara_hab_sp\\test\\grid50")     gridlist = arcpy.ListFields ("G:\\Sara_hab_sp\\test\\grid50")          for grid in gridlist:         print(grid.name)         query = 'grid = 1'         selectgrid50 = arcpy.SelectLayerByAttribute_management ("G:\\Sara_hab_sp\\test\\grid50", "NEW_SELECTION", query )      I always get an error message about the syntax of the query but i tried different things and it seems to me the error it's more related to the fact that maybe it's not possible to enter a list of fields in a query, although i have seen it's possible to do it with a list of different values of the same field... Any suggestions on how to deal with this?


title: Script to stack multiple landsat bands in different folders

text: I wanted to find out if there is a possible script I could use to stack landsat bands into a single image, and then apply that stack command to multiple folders (I have about 100 landsat scenes that are in different folders, each folder is the landsat image name, and each folder has 8 bands) Any help would be greatly appreciated. Oh and I am new to scripting language. Many thanks in advance!


title: Weekly raster data for temperature and precipitation in US

text: I am looking for weekly raster data for US for temperature and precipitation. All the major datasets PRISM, WorldClim etc. are either on a monthly scale or daily scale. I was wondering if anyone has aggregated data for weekly values from daily data or can point me to links where I can get it. I am looking for data for years 2005-2011. Thanks.


title: Problems with SEXTANTE

text: maybe the answer is somewhere, but i can't find it right now... I'm getting problems with SEXTANTE and algorithms with GRASS, SAGA... Always is the same, they can't open the files, the error is the folowing: > Oooops! The following output layers could not be open Convex Hull: > ***/home/luismi/Escritorio/para_eliminar/interpolaciones.shp*** The above > files could not be opened, which probably indicates that they were not > correctly produced by the executed algorithm Checking the log information > might help you see why those layers were not created as expected This > algorithm requires SAGA to be run. A test to check if SAGA is correctly > installed and configured in your system has been performed, with the > following result: SAGA seems to be correctly installed and configured As you can see is not a problem with spaces between names, neither is a problem of writable directory because I tried to write in Desktop folder. I'm using Ubuntu 12.04 and Qgis is 2.0 version (was updated from repositories from 1.8 and I think i had the same problem...) Please help... I need this correctly working to continue my project. Thanks to everybody!!!


title: style manipulation using Google Earth api

text: I am creating data in QGIS and then serving it using GeoServer and displaying using the Google Earth api. Is there any way to access the styles using a JavaScript api and alter the look and feel of my data, (i.e. changing colors, etc.)?


title: How to include variable in where clause of arcpy.Select_analysis?

text: I am trying to loop through a shapefile, selecting each feature in turn and copying it to a temporary shapefile to by included in a union analysis. I'm using a cursor to find the ID name for each feature which I'm setting to a varible 'Name'. Whenever I try to use this variable as part of the where clause in arcpy.Select_analysis I get an error: ExecuteError: ERROR 999999: Error executing function. An invalid SQL statement was used. An invalid SQL statement was used. Failed to execute (Select). The code I'm using is:               Name = 101     where = "\'\"StudyID\" = \\'"+str(Name)+"\\'\'"     arcpy.Select_analysis("C:\\input.shp", "C:\\output.shp", where)      If I type it out without using the variables:               arcpy.Select_analysis("C:\\input.shp", "C:\\output.shp", '"StudyID" = \'101\'')      it works fine What to I need to do to fit the variable into the sql statement?


title: Unable to add PostGIS table to Geoserver

text: I am able to add individual tables to Geoserver using,               curl -u admin:geoserver -v -XPOST -H 'Content-type: text/xml' -d '<featureType><name>quadrella_indica</name></featureType>' http://localhost:8080/geoserver/rest/workspaces/opengeo/datastores/species/featuretypes      But I need to loop through the entire list of shp file and so,               for f in *.shp        do         a=${f,,};        curl -u admin:geoserver -v -XPOST -H 'Content-type: text/xml' -d    '<featureType><name>${a%.*}</name></featureType>' http://localhost:8080/geoserver/rest/workspaces/opengeo/datastores/species/featuretypes      done      Only returns a 400 Bad request error,               * About to connect() to localhost port 8080 (#0)     *   Trying ::1... connected     * Server auth using Basic with user 'admin'     > POST /geoserver/rest/workspaces/opengeo/datastores/species/featuretypes HTTP/1.1     > Authorization: Basic YWRtaW46Z2Vvc2VydmVy     > User-Agent: curl/7.22.0 (x86_64-pc-linux-gnu) libcurl/7.22.0 OpenSSL/1.0.1 zlib/1.2.3.4 libidn/1.23 librtmp/2.3     > Host: localhost:8080     > Accept: */*     > Content-type: text/xml     > Content-Length: 47     >      * upload completely sent off: 47out of 47 bytes     < HTTP/1.1 400 Bad Request     < Server: Apache-Coyote/1.1     < Content-Type: text/plain     < Transfer-Encoding: chunked     < Date: Sun, 25 Nov 2012 18:19:31 GMT     < Connection: close     <      * Closing connection #0     Trying to create new feature type inside the store, but no attributes were specified      The above error is repeated for all the shp files that are being added to geoserver. Any idea what is wrong?


title: How to divide a feature class into ranges?

text: I have a shapefile consisting of around 4000 polygons. The polygons are created from buffering around a similar number of islands somewhere in the Norwegian archipelago. Due to an error in the dissolve functionality using arcpy dissolve fails. I am assuming this is because arcpy.Dissolve_management can not handle too many overlapping features. To try to work around this I plan on splitting my initial file into ten separate files consisting of around 400 features which are then dissolved separately. These ten files are then merged together and a final dissolve is performed. How can I do such a selection? The arcpy.select_analysis allows for a where clause. How can I use this to select the n first features, then the n+1 to n*2 features, then the n*2+1 to n*3 features etc? My final goal is to make a small arcpy function which could be a preliminary fix for us mere mortals relying on ESRIs bug fixing.


title: Coordinate system of a raster file to be introduced for transformation

text: I am working on raster data using GDAL. When I extract coordinate system information of the files I got this:   `PROJCS["Geostationary_Satellite", GEOGCS["GCS_Unknown datum based upon the custom ellipsoid", DATUM["D_Not specified (based on custom ellipsoid)", SPHEROID["Custom_ellipsoid",6378140,298.252981]], PRIMEM["Greenwich",0], UNIT["Degree",0.017453292519943295]], PROJECTION["Geostationary_Satellite"], PARAMETER["central_meridian",0], PARAMETER["satellite_height",35785831], PARAMETER["false_easting",0], PARAMETER["false_northing",0], UNIT["Meter",1]]` I have to transform the coordinate system of my other files into this coordinate system. But the problem is that I don't know how to introduce this coordinate system to my program. I can't find EPSG or Proj4 for it. Any solution or alternative?


title: How to stack labels

text: How can labels be stacked? The data has multiple points stacked at the same point location. All of the labels from all of the points at one location need to be stacked in a list next to the "group" point. thanks.


title: What are the pros and cons of PostGIS geography and geometry types?

text: My company uses `geometry` (`the_geom`) data type for storing geospatial data. I've recently been acquainted to the concept of `geography` (`the_geog`) data type which, as I understand it, stores the `SRID` along with the geometry. What are the differences between `geography` and `geometry`, and is there any advantage of using one of them in large databases?


title: How to implement ringmaps in ArcGIS?

text: In a recent article by Stewart and colleagues in IJHG I stumbled upon interesting technique of visualizing data using ringmpas. Some more information referenced in the article about this technique here and here. ![enter image description here](http://i.stack.imgur.com/fmPbg.jpg) [Source] From the article: > _Ring maps were created in Adobe Illustrator through the application of a > custom script that dynamically drew, distributed, and symbolized all graphic > map elements. The values for symbolization were read from a Comma Separated > Value (CSV) file that contained all county attribute data. Three county- > level ring maps were developed._ Is there a way to automate implementation of a map like that in ArcGIS 10?


title: Valid feature class, table names in ArcGIS

text: hi What are valid rules for name of feature class or table in ArcGIS. Is there any supporting document some body know which can explain all the limitation in ESRI platform. In ESRI support site i am kind of lost. Thanks


title: Connect a sub-model with multiple outputs to current model

text: I have a model in model builder that is clipping multiple feature classes to my study area. It is running with an iterator and creating 4 new clipped feature classes that are only represented by one output in model builder. Now I am trying to add this model as a sub model to run further tools on the individual outputs from my sub model. My problem is that the sub model only has one output and any tool that I attach to it will run for each feature class that was created from the sub model. Is there a way to have the sub model show each individual clipped feature class to connect them to separate tools?   * For Clarity* I'm currently running a bunch of clips to my study area at the start of my model that looks like this. ![enter image description here](http://i.stack.imgur.com/5U5v4.png) There area 4 feature classes being clipped in my data set. In an attempt to clean my model up, I created a sub model with an iterator to clip all those feature classes and it looks like this. ![enter image description here](http://i.stack.imgur.com/I3EpZ.png) Now all my feature classes are run through the iterator, however, there is still only one output icon in the model. So I can not connect individual outputs from that iterator to the rest of the model. If I connect that final output, it will run the model for every layer that goes through the iterator. ![enter image description here](http://i.stack.imgur.com/miElY.png) The above model is showing how the sub model is represented in my current model, but only has one output icon and therefore I can't connect individual output features from the sub model to the current model. Thanks


title: Table linked to dropbox file won't update

text: I have a table in my cartodb account linked to a file in my public dropbox folder. I have the settings in Cartodb such that the table syncs with every hour. When I look at the table in Cartodb, it says that it is syncing fine. However, it's not actually updating the table. When I make a change in the dropbox file, Cartodb doesn't register it. Example: right now, my table says it synced 14 minutes ago and that my table was last updated 3 hours ago. This is incorrect. I have continually been updating my table over the last three hours, the last time being about 30 min ago. Does anyone know why Cartodb may not be recognizing when my table is updated? Has anyone else had this problem?


title: Accessing PAT,NAT of coverages using arcobjects

text: I want to access the polygon attribute table (PAT) and node attribute table (NAT) of an arcgis coverage using arcobjects VBA, is this possible?


title: 2011 ESRI DevSummit

text: This is probably more of a community wiki question. Who will be attending the ESRI Developer Summit next month in Palm Springs? Which sessions and/or new products look most interesting to you? What tips would you give for attendees who are there for the first time? Thanks. And for the UK there is http://www.esriuk.com/trainingevents/events/dev_conference11/ [Birmingham, UK] for those cannot attend the US ESRI Summit.


title: Migrating from GIS to CAD

text: There seem to plenty of resources out there to help users of CAD software move to using desktop GIS but is the reverse also true.


title: How to get polygon coordinates for Google Map from .poly files (preferably in nodejs)?

text: I have `.poly` files with polygons of countries. I want to use those polygons to make Google Map polygon layers. How may I read `lng` and `lat` points of polygons from those poly files? They look like these:               australia_v          1          0.1446763E+03    -0.3825659E+02          0.1446693E+03    -0.3826255E+02          0.1446627E+03    -0.3825661E+02          0.1446763E+03    -0.3824465E+02          0.1446813E+03    -0.3824343E+02          0.1446824E+03    -0.3824484E+02          0.1446826E+03    -0.3825356E+02          0.1446876E+03    -0.3825210E+02          0.1446919E+03    -0.3824719E+02          0.1447006E+03    -0.3824723E+02          0.1447042E+03    -0.3825078E+02          0.1446758E+03    -0.3826229E+02          0.1446693E+03    -0.3826255E+02     END     !2          0.1422483E+03    -0.3839481E+02          0.1422436E+03    -0.3839315E+02          0.1422496E+03    -0.3839070E+02          0.1422543E+03    -0.3839025E+02          0.1422574E+03    -0.3839155E+02          0.1422467E+03    -0.3840065E+02          0.1422433E+03    -0.3840048E+02          0.1422420E+03    -0.3839857E+02          0.1422436E+03    -0.3839315E+02     END     END      Which software or nodejs module may I use to read `.poly` files to get usual `lat`, `lng` coordinates?


title: How do I prevent 3rd party application from loading custom ".style" files?

text: _Backround:_ I run Geosoft Target on my machine which has some compatibility functionality with ArcGIS. (mxd conversion, export to, etc.) I use the Geosoft Target application "stand-alone" and never use the ArcGIS funcions. _Issue:_ Target installs a ".style" file in ints instalation directory which appears in the Styel manager. If i uncheck the GEosoft Target ".style" file do clear it from the Style Manager it succeeds only in the current local session - the next time i open any other ArcMap doceument including the one in which i disabled it, the style file is automatically re-loaded. The one possible patchwork solution would be to delete the Geosoft ".style" file but I am more interested in doing it properly in ArcGIS - **is there a global setting that prevents 3rd party syle files from loading?** ![enter image description here](http://i.stack.imgur.com/8SNTn.jpg)


title: Publishing MSSQL spatial data as WFS with QGIS 2.0 results in empty layers

text: I am experiencing some very strange behaviour with a WFS layer served by QGIS Server through Apache, with data taken from an MSSQL table. There appears to be no problem at all getting the data into a QGIS project from MSSQL, but when I publish the layer it is empty -- usually. A few times over the last several days the layers have briefly worked, but it never lasts for long. What could be the problem? * * * ## The setup ### SQL Server   * Spatial data stored in a `geography` column in MS SQL Server Express 2012   * SQL View selecting the rows and columns to be exposed by the WFS ### QGIS   * QGIS project with one layer taking everything in the View, and one using a filter to select the same rows from the Table   * Only three points in the result set, and they appear correctly in the QGIS project   * Something I've noticed that seems strange to me is that each layer has two feature counts, one next to the name and one next to the feature symbol. These don't always match, as shown here:![enter image description here](http://i.stack.imgur.com/QHyXO.png)   * Sometimes I see blank/[3], sometimes [0]/[3], and sometimes [3]/[3]. I haven't been able to find out what the difference between the two is supposed to be. Project properties OWS server tab:   * `Service capabilities` is checked, with `Title`, `Organization`, `Person`, and `E-Mail` filled   * `WMS capabilities` has `CRS restrictions` checked with only the used value of `EPSG:4326` listed   * `Published` column is checked for both layers under `WFS capabilities` ### Web server   * Windows Server 2012   * Apache httpd 2.4   * `mod_fastcgi` and `qgis_mapserv.fcgi.exe` work correctly to serve raster and vector (`.TAB`) files as WMS layers ## The result When I connect to the WFS to a QGIS project, either locally on the server or remotely on my desktop, the layers are identified correctly. The field definitions, layer names, etc. appear correctly. However, most of the time there are no features present. If I pass the parameters               MAP=[QGIS Project].qgs&SERVICE=WFS&VERSION=1.0.0&REQUEST=GetFeature&TYPENAME=[Layer name]      to the QGIS Server, the response is               <wfs:FeatureCollection xmlns:wfs="http://www.opengis.net/wfs" xmlns:ogc="http://www.opengis.net/ogc" xmlns:gml="http://www.opengis.net/gml" xmlns:ows="http://www.opengis.net/ows" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:qgs="http://www.qgis.org/gml" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://www.opengis.net/wfs http://schemas.opengis.net/wfs/1.0.0/wfs.xsd http://www.qgis.org/gml http://[server]/[path to qgis fcgi]?MAP=[QGIS project].qgs&SRSNAME=EPSG:4326&SERVICE=WFS&VERSION=1.0.0&REQUEST=DescribeFeatureType&TYPENAME=[Layer name]&OUTPUTFORMAT=XMLSCHEMA">       <gml:boundedBy>         <gml:Box srsName="EPSG:4326">           <gml:coordinates cs="," ts=" ">-0.000001,-0.000001 0.000001,0.000001</gml:coordinates>         </gml:Box>       </gml:boundedBy>     </wfs:FeatureCollection>      This is definitely not the correct bounding box for the layer. The strangest thing is that every now and then, for no apparent reason that I've been able to identify, the points will appear correctly remotely. Unfortunately this always stops after awhile, again, without apparent cause. * * * I haven't been able to find any logs for QGIS other than the very minimal and infrequent messages that occasionally appear in the Log Messages panel. The Apache log is not showing any errors (which makes sense, since the FCGI process does run "just fine" and provides complete, well-formed output very quickly -- it's just demonstrably the wrong output). Are there any tricks to getting MSSQL features into WFS layers? Further debugging checks I can make?


title: How to get the intersecting polygon of two polygons? (using JSTS intersection)

text: I want to intersect two polygons within an OpenLayers layer. I found the JSTS library, which should be able to intersect two polygons, and return the new created polygon of the intersection. The function to use: intersection               polygon1.intersection(polygon2)      Using the OpenLayers Polygons directly, I get errors, resulting of the missing functions, added by JSTS. And if I use JSTS Polygons, there seems still something missing. `this.getFactory().createGeometryCollection is not a function` How do I get an intersection of those two polygons?


title: Conversion to Raster issue in ArcGIS 9.3.1

text: Trying to convert a polygon to raster in ArcGIS gives me an error... I am wondering if there is another way around to do this? I am using ArcGIS 9.3.1. The error i am getting is number 9999 when I try to use the feature class to raster tool. Here is what I did using the Conversion to Raster from polygon:   1. Choose as a feature class from a file geodatabase    2. Choose landcover as the value field    3. Select cell center type   4. 10 as a cellsize. Is the tool I am using a bug? Is there better way to do this conversion? What about a arcscript?


title: Use of Pictometry images in GIS

text: I am wondering if anyone has worked with Pictometry in GIS and their .pmi files. I recently acquired imagery from them for my area and have discovered that I cannot add the pmi files to the webmap, I presume I would have to export them. The problem is that there are literally thousands of images and the export would be unrealistic (I'm not very happy with their deliverables, wasn't involved in the acquisition process). Has anybody worked much with their imagery and have any insights in how to best use it in GIS?


title: GDAL tools for QGIS on OSX - what is the exact directory for gdal binaries?

text: When I try to use the gdal raster tools, eg clipper, I get the message "The process failed to start. Either the invoked program is missing, or you may have insufficient permissions to invoke the program." I have read the forums discussing this issue, and they seem to have solved it by using multiple paths for python and pymods, separated by "/usr/bin:/Library/Frameworks/GDAL.framework/Versions/1.7/Python/" I am not sure of what this means, but I do know that setting the path to gdal binaries in the gdal tools settings to the path above does not work. Is there a dumbed-down version of how to correctly set the paths in qgis? Any help would be greatly appreciated. The OS I am using is OSX 10.5 (Leopard), and I am using the 'GDAL complete' package. I tried versions 1.7 and 1.8, not sure if that makes a difference.


title: How to set up an OSM server with TileCache and Mapnik?

text: I have been trying for two and a half weeks so far to get a local copy of OpenStreetMap running on a server. I have downloaded the planet file and imported it into a PostGIS database called 'osm'. I have used OSM Mapnik tools to generate an XML stylesheet for Mapnik to use. I have used TileLite to prove that Mapnik can render OSM tiles from the database. The tiles even look the way that I want them to look. My problem now is that I cannot get TileCache to work with Mapnik. I have a MapServer instance installed that I am using to serve Shapefiles. This works with TileCache. The default 'basic' layer in the TileCache configuration file works as well. Please help with my OSM layer:               [osm]     type=Mapnik     mapfile=/var/maps/bin/mapnik/osm.xml     spherical_mercator=true     bbox=-16697000,8610000,-16667000,8640000     maxResolution=156543.0339/4     levels=18     srs=EPSG:900913      I have read every last blog post, forum post, and tutorial I can find. Any help would be appreciated. I suspect I have either missed something or I am doing something stupid. The problem may even be in OpenLayers.


title: How can I timestamp an arcreader print?

text: I'm setting up an published arcreader file, and was searching for a way to add a date/time stamp to the printed map. Is there a way to add a dynamic text box to a published pmf file that will show the date that the map was opened?


title: Setting "REAL" work/scratch workspace in ESRI/ERDAS etc

text: I thought the issue was only in ARCGIS but it seems that ERDAS does the same thing. I have my C with 1.7TB and my D with 10TB and my scratch set to d:\scratch but still ERDAS (and ESRI) fills my profile with it's temp data...why? see the image below (in this example I have tried to set the temp space to c:\scratch) In ARCGIS I have tried C:\Program Files (x86)\ArcGIS\Desktop10.0\Utilities\AdvancedArcMapSettings.exe--> run as admin-->System Paths-->last option "Output Temp file path" and set to D:\scratch. I have tried "solutions" from http://forums.arcgis.com/threads/3300-ArcGIS-10-Default-Geodatabase **How can we resolve this?**![enter image description here](http://i.stack.imgur.com/M24Bn.png)


title: What's a good way to display a custom raster basefile in Openlayers

text: I have a land cover raster that I'd like to display as part of an Openlayers based application. I'm able to load the raster directly as a PNG but my Openlayers application becomes VERY slow. I think I want something that can tile the raster and send level of detail tiles that Openlayers can fetch, but I'm not sure since I'm a GIS newbie. Is this something that GeoServer can handle automatically? I'd love to be able to laod rasters and shapefiles onto the server then link to those layers in Openlayers. Is this the correct approach? **Edit** : I ended up using GeoServer and adding the rasters as data sources. I was able to copy and paste a lot of code from the "Layer Preview" section of GeoServer to dump into my OpenLayers application. Thanks for the help.


title: Intermediate points on a small circle

text: I need to interpolate lat/lon points along a small circle, with a given centre and radius. I've successfully used these formulas for great circles. What is the equivalent for small circles?


title: Best practices for visualising clustered features in web maps?

text: I'm looking for the best symbol to display clustered features in the ArcGIS Server JS API. I don't like the "fly-out" animated cluster symbol, as used in the Silverlight API clustering or the JS clustering sample (they feel a bit gimmicky). So far the best option I've found is the Google Maps marker clusterer symbol. I won't be showing a number in the centre of the cluster symbol, but will be varying the size of the symbol to show concentrations of features, as in this mockup: ![enter image description here](http://i.stack.imgur.com/SclGR.png) What other clustering symbol options are there? Can you point me to a map showing a nice implementation of marker clustering? Thanks


title: How to select specific filenames nested within several folders?

text: I have hundreds of shapefiles that are stored in separate folders, which are nested within several other folders. Among the shapefiles I only need the ones that are named 'POINT*.shp' How do I get these specific shapefiles so I can then add fields, calculate fields, etc.? If there is a better option than `os.walk`, please tell. I'm only beginning to learn python, so this is very new to me. Thanks!


title: Algorithm to fix topological errors in GIS software apps

text: Related question here. I notice that ensuring topological correctness is essential for GIS applications, this is because the input from user or polygon boolean operation can have severe topological problems ( even though the polygons look correct) that would compromise on the quality of the subsequent operations. Clean polygon is how Geo Wizards do to ensure the correctness of the topology. Arcgis also has a command to clean up the slivers. My question is **_not_** about how to use the existing software packages to ensure that the polygon input data is topologically correct; rather, my question is about how these software packages implement those cleaning procedures. In other words, what is the algorithm I can use to make sure that I can fix all the topological errors, given a set of polygonal inputs?


title: Converting geographic coordinate system in R

text: I have points in geographic coordinate system and I wanted to convert them to swiss grid (CH1903+). Sample data:                id       lon      lat       2 7.173500 45.86880       3 7.172540 45.86887       4 7.171636 45.86924       5 7.180180 45.87158       6 7.178070 45.87014       7 7.177229 45.86923         8 7.175240 45.86808       9 7.181409 45.87177       10 7.179299 45.87020      Respected Results:               id       E              N           2     2579408.2431  1079721.1499     3     2579333.7158  1079729.1852     4     2579263.6502  1079770.1125     5     2579928.0358  1080028.4605     6     2579763.6471  1079868.9218     7     2579698.0756  1079767.9762     8     2579543.1019  1079640.6512     9     2580023.6226  1080049.2672     10    2579859.1889  1079875.2740


title: Using the ArcGIS Local Routing Tool . . geoprocessing in the runtime

text: When I create a model in ArcGIS using the Network Analyst Tools, I am able to create Locations and make the Route Layer and Solve. However, I need to be able to add barriers as well as locations, before solving for the Route. Is there a way to do this? Thanks! Renee


title: How to stack rasters for ecological analysis?

text: The attached image shows two different forest openings. Each opening has the following layers associated with it:   1. Point shapefile 1 = A red X which represent a large "legacy" tree in the gap   2. Point shapefile 2 = White points that indicate sapling locations   3. Raster Layer = A kernel density estimate of the magnitude per unit area of point shapefile 2 I am interested in creating a graphic that _visually_ depicts the average magnitude per unit area between the two rasters. Each raster has a different spatial extent. My thoughts were to "stack" the rasters based on the location of point shapefile 1 (i.e. the legacy tree) and perform the analysis using raster algebra (e.g. Float("Raster1" + "Raster2")/2). I should point out that the full dataset will have ~ 30 forest openings. **How can I stack these rasters so that they are centered on the legacy tree in order to perform the raster algebra?** ![enter image description here](http://i.stack.imgur.com/BwJRb.jpg) **Edit** : As @Alexandre pointed out, the solution was to georeference the rasters based on X. ![enter image description here](http://i.stack.imgur.com/Ws1PC.jpg)


title: How to configure tools in GeoExplorer composer?

text: I'm running GeoExplorer within the opengeo suite sdk, how do I override the tools configuration in composer.html ? It's only when i edit Composer.js that tools are configured the way I want. For instance, the Add Layers(+) button should not display the Upload layers menu. any changes I make in composer.html are ignored! probably this has to do with object oriented javascript which I havent fully grasped. My aim is to reorganise menus in the composer and include custom tools. Will appreciate any responses/tips. My current composer.html http://pastebin.com/HY5caSWe Thanks, Allan


title: Contour tool in ArcGIS 10.1 giving syntax errors

text: I am trying to run the contour tool on a raster in 10.1 and get the following error messages? Does anyone know if this is a bug and is there a work around? > ERROR 010328: Syntax error at or near symbol SPACE. > > ERROR 010267: Syntax error in parsing grid expression. Failed to execute > (Contour).


title: OS X Mavericks: accessing GDAL from Terminal issue

text: I successfully installed QGIS on my iMac with OS X 10.9.1 Mavericks using the installers provided by Kyngchaos. However, I'm having troubles accessing GDAL / ogr from terminal. It's giving me the following error message: -bash: ogr2ogr: command not found Does anybody have an idea of how to resolve this? Thank you, Ralf


title: Convert a SpatialPolygonsDataFrame to raster using rasterize function

text: I'm trying to convert my map (a SpatialPolygonsDataFrame) into a raster, using rasterize function, here is a bug that I don't understand...               > class(parcel) # this is my map     [1] "SpatialPolygonsDataFrame"          > attr(,"package")     [1] "sp"          > head(parcel)     AREAHA RECNO      AREA AREASQKM        OS_06 GRD_CAT_06     0 0.0519     0  519.1941        0 chemin prive        way     1 0.0127     1  127.0963        0 chemin prive        way     2 0.0596     2  596.2204        0 chemin prive        way     3 0.0248     3  248.0452        0    route D90        way     4 0.0375     4  375.0141        0 route privee        way     5 0.1260     5 1260.3941        0 chemin prive        way          > dim(parcel)     [1] 20069     6          > object.size(parcel)     88090144 bytes          > r <- raster(ncol=4001, nrow=4001)     > extent(r) <- extent(parcel)     > parcel_ras<-rasterize(parcel, r, "GRD_CAT_06", CRS("+init=epsg:27572"))     Erreur dans data.frame(ID = as.integer(putvals), value = lvs) :    arguments imply differing number of rows: 20069, 13      I just don't understand this comment : what is "putvals" ? There are 20069 polygons in my map, so what is "13" ? I can't reproduce a smaller example of my map, so here are my data : http://dl.free.fr/fxdT2oczB


title: Is it possible to show migration arrows on Leaflet map?

text: I need to develop Leaflet map which will show the whole world. Also, I need to draw some arrows which will show paths of migrations in the world. The important thing to say that those arrows need to be curved. Do you have some advice how can that be shown on the map in Leaflet API? My ideas: Is it possible to show some direction arrows with RaphaelJS? I see that Leaflet has Raphael plugin so we can add Raphael object on the map. Is it possible to draw a polygon which will be arrow shape and overlay it on the map?


title: Adding icon in particular places using ArcGIS and show data based on the location

text: I am very new to ArcGIS. Currently I am working on a project where I have to show positions of some sensors(placed in fixed locations) which are frequently(every 30 minutes) sending data to a server. I also need to show the data (2/3 values, like: temperature, PH lavel and so on) in map when anybody click to the icon of the sensor or it will automatically show the updated data table. I have basically almost no idea about ArcGIS. So it would be better if anybody give some step by step instruction or any way so that I can accomplish the task. Thanks


How can I build a network from a latitude, longitude points? I have some points(latitude and longitude) and i wish to build my network based on this points.


title: Use OpenStreetMap in Basemapgallery - Manually

text: I dont know how to add the OpenStreetMap to the basemapgallery. I want to do this manually because when i use ShowARCGIS basemaps:true i get too many basemaps i dont use. Here's my code:                       function createBasemapGallery() {             //add the basemap gallery, in this case we'll display maps from ArcGIS.com including bing maps             var basemapGallery = new esri.dijit.BasemapGallery({             showArcGISBasemaps: false,             bingMapsKey: 'BingKey',              map: map  }, "basemapGallery");              var basemapRoad = new esri.dijit.Basemap({ layers: [new esri.dijit.BasemapLayer({ type: "BingMapsRoad" })], id: "bmRoad", title: "Bing Maps Wegen", thumbnailUrl:"images/basemap_streets.jpg"              });      basemapGallery.add(basemapRoad);           var basemapAerial = new esri.dijit.Basemap({ layers: [new esri.dijit.BasemapLayer({ type: "BingMapsAerial" })], id: "bmAerial", title: "Bing Maps Luchtbeelden", thumbnailUrl:"images/basemap_imagery.jpg"              });      basemapGallery.add(basemapAerial);           var basemapHybrid = new esri.dijit.Basemap({ layers: [new esri.dijit.BasemapLayer({ type: "BingMapsHybrid" })], id: "bmHybrid", title: "Bing Maps Hybride", thumbnailUrl:"images/basemap_imagery_with_labels.png"              });           basemapGallery.add(basemapHybrid);                        var layer5 = new esri.dijit.BasemapLayer({ url:"http://server.arcgisonline.com/ArcGIS/rest/services/Canvas/World_Light_Gray_Base/MapServer" }); var basemap5 = new esri.dijit.Basemap({ layers:[layer5], title:"Lichtgrijs Canvas", thumbnailUrl:"images/light_gray.png" });      basemapGallery.add(basemap5);              var layer6 = new esri.dijit.BasemapLayer({ url:"http://server.arcgisonline.com/ArcGIS/rest/services/World_Topo_Map/MapServer" }); var basemap6 = new esri.dijit.Basemap({ layers:[layer6], title:"Topografisch", thumbnailUrl:"images/basemap_topographic.jpg" });      basemapGallery.add(basemap6);              basemapGallery.startup();                   dojo.connect(basemapGallery, "onError", function(msg) {alert(msg)});              }


Where can I download the exercise data for The Free Quantum GIS Training Manual? The pdf of the manual just mentions the sources of the data, but doesn't give any links.


title: Modifying tables based on title (MapInfo)

text: I have a query regarding the Mapbasic feature of MapInfo, I want to change the line style of an item depending on name (e.g if the .tab name contains the world group). I think I can work this out with a few bits of syntax help. So for example if I change one tab this appears in the MapBasic window               Set Window 376609128 Front     Set Map Window 376609128  Layer 10 Display Global    Global Line (2,2,16711680)    Arrows On      Would I need to write a statement similar to if contains "xx" then do that? I have a very basic knowledge of scripts syntax but is my basic understanding correct, would this be the way to do it? Thanks, Alan


title: QGIS over WatchGuard Mobile VPN with SSL

text: something quite strange occured today. I connected to my office through WatchGuard Mobile VPN with SSL enabled. I got loads of errors about SSL certificates, WMS errors (there were WMS layers in the QGIS project file) and eventually, QGIS 2.01 crashed. I got some colleagues to connect via VPN from one of our other offices and the same happend. QGIS is installed on my work laptop and this happens when I try to open figures on our network servers via VPN. When I am in the office, with the same laptop and same QGIS installation, and I try to open the same file, these errors do not occur. Therefore the problem is somewhere in the VPN connection. Why would QGIS fall over the VPN connection? Thanks Antonio


title: Leaflet: how to draw curved lines?

text: So, I'm trying to draw curved lines in my application which is based on a Leaflet map. As far as I know Leaflet does not support drawing curved lines at the moment, so my question is really more about how to best get some curved line funtionality in there. The problem is made more difficult because my code will have to be IE-8 compatible, which means that pure SVG is no solution. (Leaflet actually does a fall-back to VML if it detects that no SVG support is present...) So one possibility would be to myself write some leaflet extension code based on SVG _with its own fallback_ to VML. This would be a hell of a lot of work. :/ Does anyone have a better proposition?


title: locating latitude-longitude point in a piste-map with matrix transforms

text: I am trying to plot points of a particular latitude longitude on this image of a Swiss mountainous landscape. ![http://skimap.org/data/989/60/1218033025.jpg](http://i.stack.imgur.com/cNkdm.jpg) What I am doing is I have some 54 known location-points plotted on this image (not visible in the linked image), which is a piste map, and is therefore an artist's impression of the landscape. I found out the latitude-longitudes and also height above sea level for these 54 known points. I opened my piste map in the GIMP, and plotted the 54 points, while at the same time reading their co-ordinates in points, in the piste-map. So, I have x-y coordinates (lat- long) and z coordinate (height above sea level) for all those 54 points (so its a 3d co-ordinate system that I have), and I also have X-Y coordinates of those places in my piste-map (a 2d co-ordinate system), and I assume a Z=1 for all those points. To locate a particular lat-long on my piste-map(I have the coordinates for this place, i.e. the lat-long but not elevation), I need to calculate their X-Y on the piste-map. My approach is, that I calculate distances between this point which is to be located and those 54 known points , and then I sort them in ascending order. The first two points in this order will be the nearest points to the one which we want to locate. Now, I find my 3rd known point, and for that, I loop over the rest of the points, to get the nearest possible point such that this point forms a triangle with the first two closest points with the point which is to be located falling inside the formed triangle. Now, I have three known points forming a triangle say triangle ABC. For these known points, I have a triangle A'B'C' in the piste map, and the point to be located lies in the latter. I figure out the transform matrix which when multiplied with the 3X3 matrix of the three known points(in this 3X3 matrix, each row corresponds to a point, with the column values x,y and z in the 3d scape) results in a 3X3 matrix of those points in the piste-map. I multiply this transform matrix calculated with a new 3X3 matrix for three points, the first one being the point to be located, and the remaining two are the last two known points in the first 3X3 matrix. This way, the first two elements of my resulting matrix are the X and Y values for the point to be located, and I proceed ahead to draw the point on the piste-map. Is this approach correct with respect to the way I am using the transform matrix? Thanks! I would really appreciate if any alternative approaches are suggested as well!


title: Future of ArcGIS Mobile

text: I'm not sure we can get a straight answer from ESRI on this question at his point because there is so much uncertainty with Microsoft in this area right now but "What is the future of ArcGIS Mobile for Windows?". The concern is the diminishing future of Windows Mobile 6.x, which is probably the main platform ArcGIS mobile is designed for. When Microsoft recently announced they are changing the name of Windows Mobile to Windows Embedded Handheld, they have said a new version of this will be released in second half of 2011 and will utilize the latest development tools such as Visual Studio 2010 and Silverlight/XNA. They did say there is a clear migration path for 6.x apps but that sounds pretty vague to me. That said, along with the fact you currently need Visual Studio 2008 and WinForms for ArcGIS Mobile development, I don't see it sticking around too much longer. I see ESRI efforts going more toward the iPhone, Android, Blackberry, and WP7 API's. Maybe the WP7 API will be the future for the new Windows Embedded Handheld OS? Any thoughts, predictions?


title: map with Pacific in center

text: I would like to create a basemap with the Pacific centered. I have found this threat But it doesn't make sense to me. The last option looks the easiest. When I have loaded the map with all the countries and I choose the option 'Set project CRS from layer' (before or after shifting the map), nothing happens. When I choose "zoom to layer extent", I get the whole world again with Africa in its center. I tried the first solution as well, but what is meant with "Load NE world boundaries into QGIS" (how do you do that) ? If anybody could help ... or does anybody has a shapefile of the world with Pacific centered ? Many thanks.


title: Make QGIS interface darker in Ubuntu

text: Is there a way to make the interface of QGIS 2.0 darker in Ubuntu. I would prefer a darker theme closer to the Ubuntu Ambience style. Or is there a way to create personal themes?


title: How to install ArcGIS License Manager 9.3 silently?

text: Is there a way to install ArcGIS License Manager 9.3 silently (by using only command line)? (in windows environment xp/7).


title: Unable to expand projected co-ordinate system folder in ArcMap 10.1

text: On trying to re-project any raster from GCS to any PCS, in the "Spatial Reference Properties" window, my "Projected Co-ordinate Systems" folder is not expanding. Where as the GSC folder is expanding. Even I tried to clear the spatial property of the layer and tried to set it to the desired PCS, there also its not expanding. After restarting ArcMap, before clicking on the ArcToolBox, if I go to layer properties, there it expands. Once I clicled on ArcToolBox, it doesn;t expand either from Layer properties or from Project Raster Tool. I have attached the the image of my window. If at all any one can guide me a solution. Please.. please.. please... I'm stuck almost since last 3 months. Thanks ![Without opening ArcToolBox, From Layer Properties](http://i.stack.imgur.com/7xzcf.jpg)![After opening ArcToolBox, from Project Raster tool](http://i.stack.imgur.com/31QlP.jpg)


title: How to compose WFS GetFeature link from WMS capabilities definition?

text: This publication of geoinformation provides download links for `GetCapabilities` documents via WMS. The platform also runs a WFS service as can been seen here. I wonder if it is possbile to compose a **download link** for ShapeFiles/GML or another format which can be converted ... from the given WMS information. I managed to display the street tree information in QGIS and could load details for individual features. However, this is pixel data and I would like to download **vector data** if possible. The second GetCapabilities document defines layers named `pflanzstandorte2011` and `pflanzstandorte_aktuell`. It would be nice if these were accessible via the platform's WFS API.


title: Subsurface data in Google Earth

text: Are there any known hacks to get subsurface data into Google Earth? Sub-ocean-surface works with the release of Google Ocean, but you cannot go below the sea floor. The other major limitations of Google Earth (no vertical data, single placemark balloon) have hacks to get around them (thin buildings, custom icons). I'm hoping someone has a subsurface hack. **Edit** : Possible approaches might be   * Find the default image and replace it with a transparent PNG   * Capture the network requests to download tiles and respond with transparent PNGs   * Convince GE to set the Level of Detail (LOD) for the primary layer to some negative or infinite value. On OS X, I think the following folder/code looks promising:               cd Google Earth.app/Contents/MacOS/shaders     $ ls atmosphere_ground_sun_off*     atmosphere_ground_sun_off.arbfp1     atmosphere_ground_sun_off.asd     atmosphere_ground_sun_off.ps_2_0     atmosphere_ground_sun_off.arbvp1     atmosphere_ground_sun_off.cfg     atmosphere_ground_sun_off.vs_2_0      Those files appear to contain assembly code instructions to the graphics card, and there is mention of "groundTexture", "groundRayleighMap", and other possible hack-able things.


title: Migrating Defense Solutions objects from ArcGIS 9.3.1 to 10.0 gives Unable to cast COM object error?

text: I've got an ArcGIS Engine application that uses the IMADtedLayer interface to access VPF DTED catalog layers and adds them to a map. The code ran fine under v9.3.1. I just updated to v10, and now I get: Unable to cast COM object of type 'ESRI.ArcGIS.DefenseSolutions.MADtedLayerClass' to interface type 'ESRI.ArcGIS.DefenseSolutions.IMADtedLayer'. ... No such interface supported Here's the relevant C# code:               // geodbFWS is the IFeatureWorkspace object that contains the catalog     // sCatName is the catalog name (string)     ITable catTable = geodbFWS.OpenTable(sCatName);          if (catTable != null)     {         IMADtedLayer aMADted = new MADtedLayerClass();         bool bSetup = aMADted.Setup(catTable);     }      The error occurs on the "bool bsetup" line. I've tried substituting new MADtedLayer(); for new MADtedLayerClass(); - with the same result. I've also uninstalled all ArcGIS components, cleaned out the registry and re-installed. I get the same error. I'm using MSVisualStudio 2010, and the intellisense recognizes the IMADtedLayer interface, the MADtedLayerClass, and the .Setup method - as well as the fact that it should return a boolean. I'm running the Windows 7 Professional Operating System. Has anyone else encountered this problem? Any insights will be greatly appreciated!


title: New python GIS Programmer : Fundamentals of the Python Language for ArcGIS

text: I am new to programming and I'm trying to learn "Fundamentals of the Python Language for ArcGIS" and then I've been trying to compile an example code but python GUI keeps throwing throwing this error **"Subprocess startup error"** IDLE's subprocess didnt make connection. either IDLE cant start a subprocess or personal firewall software is blocking the connection. **and the second error message** Socket error: No connection could be made because the target machine actively refused it now I use avast anvirus and I have checked cant seem to find how it may have been blocking the process so I thought maybe someone can help out. Thanks


title: How to make a custom ToolbarControl Item in ArcEngine?

text: I would like to create a custom toolbar item in ArcEngine. I add items like this ((AxToolbarControl)toolbarControl).AddItem( myCustomItem ) I need to know how to create a custom item in ArcEngine 10. I'm using C# and VS2010. Thanks.


title: When is it safe to use ObjectIDs as keys in a relationship class?

text: I've always avoided using ObjectIDs when defining relationship classes. It is my understanding that some actions on a geodatabase will re-set objectIDs (like import/export). However I do see a lot of people using them. Am I just being paranoid, or are there cases when it is safe to use ObjectIDs when defining a relationship class?


title: How to add a Database Server in Catalog?

text: I have ArcSDE on an Oracle database and I want to add a database server connection to ArcCatalog, I've seen several sites but they don't say what I have to write down. Let me explain: From the tree at the left I have several options, from hard drives, database connections to database servers and gis servers. When I right clic the Database servers option and choose add database server I get a dialog box asking me for a Database Server, so, what is what I have to write there? I've tried with:   * servername\user   * servername\oracle11g   * servername\ESRI_sde Thanks. And SDE Administration


title: v.dissolve (and other scripts) not in GRASS path

text: I see that I have the v.dissolve module in my grass64/scripts dir, but it is not found in my path from the GRASS command line. I have the GRASS_ADDON_PATH environment variable set: GRASS_ADDON_PATH=/usr/local/grass-6.4.1/scripts/ And the module itself is located in: /usr/local/grass-6.4.1/scripts/v.dissolve What am I missing? Thanks, Roger


title: Obtaining a contour as a single polygon in ArcGIS?

text: I have a DEM that I want to retrieve a single contour polygon from at a given elevation. I can use the Contour tool to obtain a series of polylines at a given elevation by setting the contour interval to a point where only one contour is produced (e.g. 9999m) at the base contour interval (e.g. 105m). What I want is a polygon representing the specified elevation, which happens to be the average water surface elevation of a river. You can see an example of my data in the image below. Any thoughts on how I can accomplish this? The polyline contours run up to the edge of the DEM and then break off, so there is no continuous line which surrounds the whole river. The second image is essentially what I would like to be output. You can see that I'll also have to solve the problem of islands within the river channel. I am attempting to automate/model this task so ArcObjects are fine. ![Contour1](http://i.stack.imgur.com/gS0hG.png) ![Contour2](http://i.stack.imgur.com/Qr8Qx.png)


title: Speed of various raster data formats

text: I'm having trouble locating any discussion or comparative benchmarking of different raster file formats (e.g., for use in data analysis in R). Does anybody have any insight into why particular formats may be faster or slower? Or should the differences be minimal? Specifically, I'm interested in if converting a raster (e.g., a GEOTIFF file) to a different format (e.g, a netCDF) is ever worthwhile for the purpose of speeding read/write and other operations. Thanks in advance! Best, George


title: How to set up a Revision Control for ESRI File Geodatabase?

text: I'm interested in doing revision control (Mercurial) on some file geodatabases. Does anyone have an idea of which files are important (need to be saved) and which can be ignored? Note that I'm not interested in ESRI database versioning. These are essentially empty (no data), and I want to make sure they are in sync with code changes. For example, if I add a new feature to the database, I add some things in code that uses that feature. If I want to go back to an old revision of the code, I want the database that matches. I don't even need to be able to compare versions, just get back to a particular one if needed. Some of the files include: gdb, timestamp, *.freelist, *.gdbindexes, *.gdbtable, *.gdbtablx, *.spx, *.atx


title: How can I find points inside a polygon in ArcGIS server?

text: I have a feature class of points. My client wants to draw a polygon on the map to select the subset of points inside the drawn polygon. Is there any tools in ArcGIS API, which I can use it for this purpose? if there is a sample for ArcGIS API for Flex, that would be great.


title: ModelBuilder Output Creates .lock file and Cannot Overwrite When Ran a Second Time

text: I am having trouble overwriting an output shapefile from ModelBuilder. The model works great when ran the first time, but I get an ' **000210 : Cannot create output .** ' error if I re-run the model. The error states 'The output cannot be created. Potential reasons include data locking, an incorrect path, and limited access rights..' I've noticed that there is a .lock file associated with the shapefile so I assume this is my problem. I've checked the overwrite outputs of geoprocessing options checkbox in the geoprocessing options. My model is very simple, I am just using a multivalue iteration to buffer a line segment at different buffer widths and then clipping to a polygon. I am using the inline variable subsitution based off of the buffer value for my output shapefile name. Anybody know how to fix this? or possible work arounds? I am using ArcMap 10.2


title: How to merge hundreds of vector layers and sum attributes in QGIS?

text: I'm still at a fairly beginner level with QGIS and I've been banging my head against the wall trying to figure this out, so I would like to turn here for the internetz help. I am trying to construct a map of mean daily rainfall in Africa over the past ten years. I have access to shapefiles for daily rainfall (with identical columns in the attribute table) and I would like to sum up each day for an entire year period. Then, I want to take the mean of the last ten years and create one shapefile. So far, I've been working with only two specific days with the thought being that if I can combine two, I can combine 365. I've tried Vector>Data Management Tools>Merge shapefiles to one, but this just creates one larger shapefile and doesn't sum up overlapping areas. I've also tried Vector>Data Management Tools>Join attributes by location, which gets closer to what I want, but only keeps the geometry of the target vector layer and loses that of the other. Detailed guidance would be sincerely appreciated. Please let me know if any of this isn't clear. Also, let me know if there may be a quicker way to go about creating a shapefile of the ten year average other than what I've proposed. Below are a couple pics: Project with two layers to merge (I added the continent outline for a bit of context) ![enter image description here](http://i.stack.imgur.com/qEBZk.jpg) Sample of layer attribute table (I only need to sum the values from the "Contour" column--this is daily rainfall in mm) ![enter image description here](http://i.stack.imgur.com/ClchQ.jpg)


title: How do I query for geometry inside the Arctic Circle?

text: I have a large PostGIS database of lat/lon points. How can I query for any points that lie above the Arctic Circle?


title: Problem with loading GeoJSON file in OpenLayers generated by PHP

text: I'm trying to read some geo-data from PostGIS database using PHP, creating a GeoJSON file and showing it on a base map using OpenLayers. The PHP code works just fine. The generated data in GeoJSON format is valid (I checked it with jsonlint). When I load the page, the developer tools shows no errors, but also doesn't load any data on the map. I have just started learning OpenLayers and don't know much about it but saw some examples for this and tried to do the same. Am I doing something wrong in the Openlayers code? The code is:               var map = new OpenLayers.Map({         div: "map",         layers: [             new OpenLayers.Layer.WMS(                 "WMS", "http://vmap0.tiles.osgeo.org/wms/vmap0",                 {layers: "basic"}             ),             new OpenLayers.Layer.Vector("GJ", {                 strategies: [new OpenLayers.Strategy.Fixed()],                 protocol: new OpenLayers.Protocol.HTTP({                     url: "gj/DB2GeoJson.php",                     format: new OpenLayers.Format.GeoJSON({                         extractStyles: true,                          extractAttributes: true,                         maxDepth: 2                     })                 })             })         ],         center: new OpenLayers.LonLat(-112.169, 36.099),         zoom: 11     });      and the GeoJson file generated is:               {     "type": "FeatureCollection",     "features": [{         "type": "feature",         "id": "97",         "geometry": null,         "properties": {             "popn_total": ["97"],             "boroname": ["Staten Island"]         }     }, {         "type": "feature",         "id": "66",         "geometry": null,         "properties": {             "popn_total": ["66"],             "boroname": ["Staten Island"]         }     }, {         "geometry": ["{\"type\":\"MultiPolygon\",\"coordinates\":[[[[578620.717363209,4495974.81786636],[578535.357633272,4495915.80999248],[578524.626146389,4495908.44992811],[578475.151658836,4495981.94016235],[578550.53822566,4496033.27562186],[578571.6880873,4496047.75223892],[578620.717363209,4495974.81786636]]]]}"]     }]      }


title: Sextante plugin

text: I'm looking for the sextant plug-in which is supposed to be showing in Plugin Manager under Manage and Install plugins. I currently have 30 plugins installed and enabled in QGIS 2.0 and 121 more showing as available so repository access looks ok. However cannot see sextant in here although have missed the obvious before. Not sure if other settings are ok. Running Windows 7. It has been suggested a reinstall using OSGeo4W Advance mode would possibly sort this. Help appreciated. Thanks UPDATE: So Sextante is called Processing in QGIS 2.0. Is the Processing toolbox a separate standalone within QGIS and not directly related to the base QGIS Vector / Raster, etc tools?


title: Get folder structure in code from arcgis/rest/services

text: I have never worked with maps before, and I got assigned on this project with ArcGIS. Noone to ask, cause noone has worked here either so I am turning to you guys. I have to create a tree of checkboxes, and clicking them turns the visibility of the clicked tile on or off. I have managed to find some methods, `ArcGISTiledMapServiceLayer` and `ArcGISDynamicMapServiceLayer`. but they don't help me because I cannot give them a link like: `http://sampleserver3.arcgisonline.com/ArcGIS/rest/services`. but what I need is to take all the folder names from there, ![enter image description here](http://i.stack.imgur.com/5Zi25.png) make them the parents in the tree, then get the services that interest me (in this case there are 3 names in the first folder, don't know which one would be a tile) ![enter image description here](http://i.stack.imgur.com/mlArE.jpg) and put them as leaves in the tree. An example of the behavior I need to reproduce is this: ![enter image description here](http://i.stack.imgur.com/qjYjj.jpg) in the link I pasted, go to `BloomfieldHillsMichigan` ,then `BloomfieldHillsMichigan/Parcels`, click on `ArcGIS.com Map` from the "View In:" secition, then the "Details" tab, then "Content". And there is the checkbox I need to reproduce, and the text is "Parcels". And all of this has to be in a tree structure, with all the folders we had in the folder section. Is there a way I could accomplish this? Thanks


title: how to choose stratified random sample on multispectral image

text: I'd like to choose random samples from multispectral image according to each Land Cover image. The multiband image has 70 bands representing land surface temperature data from MODIS. Land cover image has 9 classes. I need 50 sample (pixel) values of each land cover type which extracted from each band of multispectral image. How can I choose these random samples in ARCGIS, ERDAS and EVI? ENVI can be used to choose random samples on multiband image as ROI. But statistic output of ROI gave me the descriptive statistics of mean of all sample. What I want is the mean value of each separate samples of each land cover types in each bands. Thanks in advance


title: Calculating surface and lines in 3D projection in R

text: For my internship research I am trying to compare the surfacearea and length of geometric shapes and lines in one coordinate system to another. I have a SpatialPolygonDataFrame in the European equal area projection ETRS-LAEA and I have a function that can calculate the surfacearea and the length of all polygons. I can transform it to another projection and calculate the surfaceareas and the lengths of all polygons in that projection. My next step would be to calculate the surfacearea and the length of the polygons in ETRS89(x,y,z). Since this is a 3D projection the functions I am using for the 2D projections are not working. Are there functions for 3D projections? I can not find them online. I am hoping there has been some smart guy with a similar problem and already fixed it. I am expecting this is not the case and that it is (for me) impossible in R.


title: How to access input parameters (layer name,field name) in rest SOE?

text: I want to access input parameters in my SOE(Server object extension) code which I am writing in arcobjects using C#. I have constructed the complete URL(including layername id and field name) for calling SOE in my silverlight web application. I want to know that how to get these input parameters to use in code. As I am beginner in SOE(Server object extension) development I need a direction for developing SOE. Please help me out to move further as I am badly stuck with this issue.


title: Unable to load 'nviz' through GRASS

text: I am using QGIS 2.0.1-Dufouron on Windows 7. When I try to load nviz via GRASS, I get the message "Module nviz not found". Weirdly enough, I can load it via the toolbox and be able to 'run' it (it reaches 100%). But afterwards, nothing happens. **Problem:** I merely want to rotate and play around with the view points of layers. This is all for aesthetic purposes, don't want to do anything technical (not yet anyway). I have read posts on people suggesting to install the plugin SEXTANTE which may solve the issue, however it doesn't show up on my list of plugins (maybe I don't have the correct URL?). I am considering trying ArcMap if I have to but would rather do it in QGIS. Any suggestions please?


title: How can I translate a GML definition of SRS to a usual format (e.g. WKT or PROJ.4)?

text: From GDAL 1.9.0 a new utility is available: **_gdalsrsinfo_**. Although it can yield an XML (GML-based) definition from a source equivalent format or file, it seems that it can only _read_ `usual GDAL/OGR forms (complete WKT, PROJ.4, EPSG:n or a file containing the SRS)`, so no GML definitions. Indeed:               $ gdalsrsinfo -o xml '+proj=longlat +a=3396190 +b=3396190 +no_defs'     <gml:GeographicCRS...      but:               $ gdalsrsinfo -o wkt '<gml:GeographicCRS...</gml:GeographicCRS>'     ERROR 1: ERROR - failed to load SRS definition from <gml:GeographicCRS...      Any of you is aware of a way I could make this translation? This question is somehow linked to: Can GeoTools parse GML CRS definitions?


title: Click a marker icon to open a pop-up, using OpenLayers

text: # Background Looking to create a pop-up that performs an Ajax request when the user clicks a marker on an OpenLayers map. The markers, however, have different icons. The following is working code that can be used to accomplish the desired task:               var sites new OpenLayers.Layer.Markers( 'Sites' );     map.addLayer( sites );     addSiteMarkers( sites );          function addMarker( layer, lat, lng, img ) {       var size   = new OpenLayers.Size( 32, 37 );       var offset = new OpenLayers.Pixel( -(size.w/2), -size.h );       var icon   = new OpenLayers.Icon( 'img/marker/' + img, size, offset );       var marker = new OpenLayers.Marker( getLocation( lat, lng ), icon );            layer.addMarker( marker );            marker.events.register( 'click', layer, clickMarker );     }          function clickMarker() {       alert( 'click' );            // Pop-up code will go here...     }          function addSiteMarkers( sites ) {       addMarker( sites, -113.3017895, 49.1939798, 'store.png' );       addMarker( sites, -113.3072538, 49.2013626, 'building.png' );     }      The above code adds two markers (of hundreds) each having a different icon (from a set of fifty different icons). Each click will display a pop-up and each pop-up will have information that is retrieved from a database (using a jQuery ajax call). # Problem I suspect this is a terribly inefficient way to use OpenLayers and map markers for opening pop-ups when a user clicks an icon. # Question How do you use OpenLayers to present a variety of different markers that when clicked display a pop-up (with content dependent on the actual marker that was clicked)? ## Update #1 Most of the examples show how to use a Vector layer, rather than a Marker layer, to handle click events. Using a Vector layer makes making markers more mind-boggling. # Related   * http://dev.openlayers.org/releases/OpenLayers-2.11/examples/marker-shadow.html   * http://stackoverflow.com/questions/7702945/adding-popup-info-in-openlayers   * http://stackoverflow.com/questions/9895951/how-to-display-multiple-markers-with-individual-framecloud-popups-in-openlayers   * http://stackoverflow.com/questions/6026420/openlayers-vector-features-instead-of-markers   * http://dev.openlayers.org/docs/files/OpenLayers/Control/SelectFeature-js.html   * http://docs.openlayers.org/library/overlays.html   * http://dev.openlayers.org/docs/files/OpenLayers/Popup-js.html   * http://dev.openlayers.org/docs/files/OpenLayers/Marker-js.html Thank you!


title: Implementing custom (Vincenty/Haversine formula) functions for MySQL vs Using PostGIS?

text: I currently am using MySQL for all of my relational data. I have hit a crossroads with geo calculations. Since, MySQL does not have good support for spatial functions, I am contemplating on using Postgres's PostGIS extention. But, for this I will have top run 2 DB servers (MySQL and Postgres) as I do not want to move away from _InnoDb_ (for all the **goods** that it offers). >   * Does implementing the functions cause any performance issues(lag)? >   * Is the trouble worth it or should I move ahead with PostGIS? >   * Is moving the whole DB from MySQL to Postgres any good (Although I am > not exactly considering this)? > > * * * > >   * Also, is implementing Haversine or Vincenty better? I do understand that > Vincenty formula provides much better accuracy. But is there too much cost > on CPU? > >


title: What does maps.nokia.com use?

text: The maps on maps.nokia.com are pretty amazing. I assume they're generating their own tiles but, does anyone know about the infrastructure behind site: the tile server used? and their JavaScript Library?


title: QGIS Grid Division/Grid Product etc. results are wrong projection

text: I have several raster grids, all 4320 columns, 2160 rows in WGS 84 projection. They are loaded from ESRI ascii files. When I use the Processing Toolbox to carry out arithmetical transformations on the grids (e.g. multiply one by another, divide one by another) the resulting grid is correct, but offset from the original projection by some number of pixels. The layer purports to be in WGS 84, but clearly isn't. Currently I am working round the problem by using the raster package in R, but it seems to me that the QGIS functions are broken (unless I am missing something obvious). Any suggestions gratefully received.


title: Can I use google api or google maps in web application developed with arcgis api or esri technology?

text: I need to know that is it possible to use google map (as a basemap) and its GIS functionalities(like find/search features) in my web application developed with arcgis api for silverlight/javascript. Also I dont want to use arcgis server. I want to use open source GIS technology.


title: displaying combined data values in ArcMap

text: I got some data in an attributes table and I would like to know how I can combine them in order to get an unique value in the map. For instance,   * Attribute: Type   * Value: Street Those values I want to represent in the following format: **Type / Street** Thank you in advance.


title: How to calculate the intersection of 2 Circles?

text: I'm trying to figure out how to mathematically derive the common points of two intersecting circles on earth's surface given a center Lat/Lon and a radius for each point. Ex. Given:   * Lat/Lon (37.673442, -90.234036) Radius 107.5 NM   * Lat/Lon (36.109997, -90.953669) Radius 145 NM I should find two intersection points with one of them being (36.948, -088.158). It would be trivially easy to solve this on a flat plane but I don't have any experience solving equations on an imperfect sphere such as the earth' surface. Any help would be appreciated!


title: How to sort a field according to value range?

text: I have to check an attribute table field (MY ID NO) of polygon feature against a tiff image of parcel map where fields are written as 12-58,65,87,2565-89999,58648. Example: Parcel map tiff(image) contains plot no (MY ID NO in poly feature's attribute table) against each plot, each tiff has a corner box that describe the plot numbers included in this tiff map. Now I have to check is all the opposite plot no (collected from tiff image file, and also mentioned in the corner-box describing included plot no) is populated in the MY ID NO field. Now I want a field calculator script (I think it is easy and better way to check, suggestions also preferred) to populate "Sorted Range" field. **I need only "Sorted Range" field**. My table as below: ![Portion of Attribute Table](http://i.stack.imgur.com/jd9RZ.png)


title: What is the best quick way to populate a field (in ArcGIS 10)? (for rare language/font e.g. Guzrati numerals)

text: * I have polygons generated from shapefile (esri) format    * I have a tiff image of plot/parcel map(this tiff has plot id scribed in somewhat consecutive in manner) of that polygon feature. I would like to fill the `plot number` field by juxtaposing(placing two abreast) the polygon and image file. I am using Arc GIS 10; what is the best way to do this? Example: I have following tiff map, where from a polyline-shapefile has been produced by arc-scan, thereafter a polygon feature has been made from it. ![My dummy tiff image](http://i.stack.imgur.com/wAvHS.png) Now i want to populate "Plot No" field of polygon feature from above tiff image. But now "see-and-fill" method taking much more time (since 1700+ polygon/plot is many). Under this circumstances i seek the best and quick (some what automated) way to populate "Plot No" field. N.B. Plot No in the map is shown English Numerals, but in case of me these numbers are font that parochial (even abbey software does not support this font. e.g. Guzrati "╜") Thanks in advance


title: Openlayers Measure Line and Area at the same time from the same button

text: I have been using separate buttons to measure lines and areas in my Openlayers maps (as per the Openlayers examples) but would ideally like to get the output for both length and area at the same time. The ideal solution would be to have a solid line for the line measurements and a dashed line connecting the start and end point to symbolize the polygon. Then simply output the two answers to a div on the website. As a bonus, store each line length in an array that can be output to a menu as well. My question is two fold, is it possible, and does anyone have / know of an example? I've searched long and hard and can't find an attractive solution to this. Thanks for your help,


title: Calling rest services from WCF service

text: I have a question. I have an AddIn that allows spatial query using map service. I have logic that connects to map service on the client side. I am looking into the possibility of moving this logic to a WCF service so that I could use the same logic for more than one application. I am not sure if I could use rest service from WCF service, based on my search it seems like this is not an option. Is SOAP service is only way to go or if there is a way to use rest service. Appreciate any feedback. Thanks Jay


title: Split a polygon based on density of points in ArcGIS 10.1

text: I have two layers:   1. Polygon layer representing administrative boundaries   2. Polygon layer representing topography (i.e. a map) In some of the administrative areas, the buildings (shown in the map) are grouped in one area while the rest of the area is something else. I would like to split the administrative area into two, one containing all the buildings and another containing no buildings. See diagram. ![Example of how the administrative areas should be cut](http://i.stack.imgur.com/OmnkV.png) I created a union between my two layers and then dissolved to create a third layer which identified the "building/not building" for each administrative area, this is close to what I want except that it doesn't create a single cut, but creates a complex shape that wraps around each building and go along roads etc. I have a few thousand polygons like this so I would like to be able to cut them automatically.


title: Migrating Feature Datasets from ArcSDE to postgis

text: I am currently planning the migration of a database from ArcSDE/ Arcmap to postgis / Qgis. The data is stored in raster catalogues and feature datasets. Is there an equivalent to feature datasets supported by postgis or do I need to group the layers by adding a prefix corresponding to the current feature dataset.


title: Are GIS Analyst salaries going to always be so low?

text: I regularly see GIS Analyst job postings (often in very expensive cities) asking for a plethora of experience paying in the $40k - 60k a year range, many are even below $40k or $30k. Considering the amount of patience it takes to deal with Esri software and all of the work-arounds one must routinely implement GIS professionals are being severely underpaid. Are the majority of GIS Analysts independently wealthy and just do GIS for fun? I do not see any other logical reason for the low pay reality GIS Analysts are in. Has there ever been an attempt to create a GIS Union?


title: Unable to save to in Geojson format with the correct projection

text: I have a shape file (Polygon) in Qgis that I wish to save in Geojson in order to display it with the mapbox javascript api. Original projection is in POSGAR 97 My issue is that it does not reproject when I set projection to WGS84. Even if in the T.O.C I set the projection to WGS84 and **Save As**. Even if I uncheck the on the fly re-projection. This is how coordinate looks like:               "coordinates": [ [ [ 324345.0, -3579435.0 ]      It should looks like:               "coordinates": [ [ [ -32.43450, -35.794350 ]      Curious think is if I **Save As** in KML, it works and the point of the polygon looks like:               <coordinates>-64.866399861691903,-32.338025763844477...


title: Using the feature class name to populate a field in the features attribute table

text: I'm trying to duplicate the model shown in Use Feature Class "Name" from Iterate Feature Classes to Calculate Field with Model Builder 10.1 I have over 700 feature classes in a gdb and I need to use the name of the feature class to populate a new field within the feature class. I tried to duplicate the model above and I cant figure it out. Can someone help? What is The Feature Type Parameter?


title: Categorized style hiding map

text: I am having an issue with the style properties of my map. I imported the US Zip Code Tabulation Area map from a government website, and in a separate excel file I have a list of people connected to various zip codes. My original problem was QGIS was removing the leading zeros from my CSV file, which caused a lot of zip codes to appear blank. I created a .csvt file to declare the zip code column as a string, and this worked and allowed me to keep the leading zeros. However, now when I join this data to the zip code map, and then try to style based off of the 'name' category, the map disappears. I checked the attribute table, and all the names are correctly assigned to the various zip codes, but the map disappears anyway. The map previously worked before I added the .csvt file but obviously with the missing areas due to leading zero issues. I am not sure why the map is disappearing after adding the .csvt file. I've included images to show how everything appears to be correct in the properties and attribute table, but the map just won't show up. Here is the album where I show my process and the error.


title: ArcEngine 10 Editor:  Split/Rotate/Cut Features toolbar?

text: I am currently looking for the _Split/Rotate/Cut/Reshape_ tools in the toolbar control provided by ESRI, however according to ArcGIS Engine commands page there isn't one. Is it called something else or does it only exist in ArcEditor?


title: Drawing Day and Night on a Google Map

text: I'm looking to plot day/night on a Google map, for an arbitrary point in time. I'm familiar with generating map tiles; I'm just looking for an algorithm to tell me whether a particular point on the globe is currently in daylight or darkness, or to otherwise plot the curve of the day/night interface onto the map. I've done some searching, but it's possible I don't even know enough about the problem domain here to know what terms to search for! Any ideas? Doesn't have to be perfect -- basically, I'm comparing Flickr geolocation data of sunrise and sunset photos (and their "date taken" timestamps) with reality, and this is to help me visualise it.


title: How do I randomly alter a raster map of habitat types?

text: I have a raster of habitat types for a specific area in Scotland. I need to create future habitat scenarios with changes in habitat to assess the population viability of a bird species. For example, in the future there might be 10% more forestry in the area. I would like to alter the current map by randomly adding the forestry in blocks of a certain size. I am, so far, thinking along the lines of selecting random points from a raster which identifies areas where the forestry could occur and growing the correct sized blocks using some sort of cellular automata. Does this seem like the best way of going about this? Is there a better method? If this is the best way available, how could I do this in, preferably, R? (I am currently looking at the rpoints function in "spatstat" along with the CellularAutomata package) I also have access to GRASS, QGis and ArcMap 10 if there are simpler ways in any of them. Thanks.


title: How to import ESRI grid into Matlab

text: I have some ESRI grid data that I would like to import into Matlab. As a typical ESRI grid data I have the two folders with the *.adf files and *.dat, _.nit and_.dat files. I would like to know how can I import the information stored in those files into Matlab. Tank you all!


title: How to select features with NULL attributes in QGIS?

text: I need to identify which features have null attributes in the attribute table of a shapefile. When i try [field]= 'NULL' it selects the features with value=0 (see picture)![enter image description here](http://i.stack.imgur.com/gI8vo.png) Suggestions?


title: Invalid filter when using MapServer as a WFS

text:
I'm trying to use Openlayers and Mapserver to display a shapefile. My HTML file is as follows:                           <title>OpenLayers Example</title>

                       This loads a google-map base map, and should overlay a shapefile over all of England. When the page loads the coordinates of the bounding box of the map are sent to MapServer. I can see that it is sending just the lon/lat values. In my mapfile, I have the projection of the data defined as follows, in both the intro part, and also in the layer:                PROJECTION         "init=epsg:27700"      END      However, the response from the server gives me an error `WFS server error. Invalid or Unsupported FILTER in GetFeature`, and I don't know whether it is due to an incorrect bounding box, or whether I have a mistake with my query using `PropertyName`. From searching online this appears to have been a bug, but I don't know whether it is still a problem or not.               <?xml version="1.0" encoding="ISO-8859-1"?>     <ows:ExceptionReport xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"      xmlns:ows="http://www.opengis.net/ows" version="1.1.0" language="en-US"      xsi:schemaLocation="http://www.opengis.net/ows http://schemas.opengis.net/ows/1.0.0/owsExceptionReport.xsd">       <ows:Exception exceptionCode="InvalidParameterValue" locator="filter">         <ows:ExceptionText>msWFSGetFeature(): WFS server error. Invalid or Unsupported FILTER in      GetFeature : <         Filter ogc="http://www.opengis.net/ogc">       <BBOX>         <Envelope gml="http://www.opengis.net/gml" srsName="EPSG:4326">           <lowerCorner>-55.50109375 43.134977942511</lowerCorner>           <upperCorner>55.24109375 59.593781185019</upperCorner>         </Envelope>       </BBOX>     </Filter>     </ows:ExceptionText>       </ows:Exception>     </ows:ExceptionReport>      What do I need to do to fix this? Should I make sure that I am passing a bounding box in `epsg:27700` coordinates? Or can Mapserver accept lon/lat and project them into the required coordinate system? * * * I have asked another question on the same topic recently, but I had a different problem then!


title: How do I find the average distance within a polygon to a point?

text: For my purposes, I can assume a homogeneous distribution of points within the polygon. What I want to find is the average distance from points within a polygon to a specified point. I don't know where to start and have found similar info online but nothing exact. Ideally, I would like to do this for an entire Thiessen polygon layer. I have many polygons and a layer with one point in each polygon, and within each of those I would like to find the average distance to the point. If you are confused by my question, here's an example: Within a country, what is the average distance one must travel to reach the Capitol? It can be assumed that people are evenly distributed throughout the country.