title: HSLayers WPS GUI problem

text: I have a GeoServer running, I have made a process that simply returns the requested tiff. Now here's the deal : when I send my request from the integrated WPS demo in GeoServer, the GUI gives me a list of available maps to choose from. However, when I try the same thing with my (working) HSLayers client (a simple HTML page with an embedded js WPS client), the GUI asks me for an URL and I have no list of maps. I would prefer having them shown. Additionnaly, if I take any random process, the GUI will work fine. If I switch between processes, everything works fine (updating input fields, etc.). But once I switch to my AddGeoTIFF process, the GUI generator bugs and doesn't switch anymore (I have to reload the page). Is this due to a small error in the HSLayers client or is it a problem in the prototyping of my process ? Here's my process' code :               @DescribeProcess(title = "AddGeoTIFF", description = "Returns a GeoTIFF\n")     public class AddGeoTIFF implements GSProcess {          private static final CoverageProcessor PROCESSOR = CoverageProcessor.getInstance();          @DescribeResult(name = "result", description = "The resulting GeoTIFF")     public GridCoverage2D execute(             @DescribeParameter(name = "coverage", description = "The GeoTIFF you wish to obtain") GridCoverage2D coverageA,             ProgressListener progressListener) throws ProcessException {                 return (GridCoverage2D) coverageA;         }     }      I am not good with Java and based it on an existing process, so there might be some horrors. I have been looking for the cause of the issue, and it is (what I believe to be) a parsing error. The (on the fly) GUI generator's part which troubles me can be found in `[HSLayers install file]\source\HSLayers\WPSClient\complexField.js`, on lines 77 to 82. If someone experienced could have a look at it, It'd help. Thanks


title: GeoNetwork map viewer not displaying title name in Thai language?

text: I am using **GeoNetwork 2.10.3** to **Add WMS** server from **GeoServer** and am unable to display its title using the **Thai** language. I have tested an edit string with the Thai language and it worked. It is only the map viewer which cannot display it. Can you help me to fix this problem? ![enter image description here](http://i.stack.imgur.com/voy8g.png)


title: What is a good method for serving many static map layers through ArcGIS Server 10?

text: I'm writing a Silverlight app that consumes some map services from ArcGIS Server 10.0 SP2 and I have about 15 layers of historical imagery I would like to give the user access to. This imagery doesn't change and I'd prefer to use a cached map service so the user gets a smooth panning/loading experience. The catch is the user needs to be able to have any two layers overlaid on each other so they can blend between the two. I can see two options to do this:   * Create 15 cached map services, one for each layer   * Create one multi-layer cached map service I've read several threads where it has been stated that multi-layer caches are not for web applications (and it appears they use the DynamicMapService from the API so they may not even satisfy all my requirements). I'm sure this isn't an unusual situation to be in and I'm wondering what the typical solution might be. Right now the layers are dynamic layers and it doesn't provide as smooth a user experience as I'd like.


title: What is ArcGIS Market place?

text: Does any one know (OR using OR have experience) exactly what is ArcGIS Market place? According to Google search its related with selling apps and GIS data built using ESRI software's. So any additional information/help will be great !! thanks in advance !!


title: How to enforce labels for all polygons?

text: I am trying to label polygons but it does not always work. Either the polygon is quite large, then the label appears, or the polygon is small, then the label is not displayed. I've already tried to lower label size, disable box about overlap labels, disable scale limits, yet, nothing works. If I want the label being displayed, I have to use the node tool in order to make the polygon larger ... this is not the good solution. Has someone already met this problem? (I have looked in other questions, but only polygons geometry validity is mentioned, and I have already checked it)


title: Equivalent function to shapely's envelope in ogr?

text: I wrote a very simple python function that gives me the envelope/bounding box/...based on two coordinate pairs (ower left and upper right corner):               # ...     def get_bounding_box(self, ll_corner, ur_corner):         try:             from shapely.geometry import MultiPoint             mp_env = MultiPoint([ll_corner, ur_corner]).envelope             return ogr.CreateGeometryFromWkt( mp_env.to_wkt() )          print mp_env     POLYGON ((4.8439699999999997 52.3661099999999990, 4.8469199999999999 52.3661099999999990, 4.8469199999999999 52.3706000000000031, 4.8439699999999997 52.3706000000000031, 4.8439699999999997 52.3661099999999990))      As you can see I'm returning an ogr geometry which I use, for instance, in this manner:               bbox = geo.get_bounding_box(lowerleft,upper_right)      ...     p=ogr.Geometry(ogr.wkbPoint)     p.AddPoint(x,y)     if p.Within(bbox):         ....      Isn't there an equivalent function to shapely's envelope in ogr? Both `.ConvexHull()` and `GetEnvelope()` return only the two corner coordinate pairs.


title: GRASS 6.4.2 define Geographic boundaries using MSYS

text: I am working with daily gridded rainfall data of 50 years. The initial data had the extension .nc.gz.tar which I have unzipped using MSYS in GRASS 6.4.2. Now I have 365*50 days of daily data. My problem is that my data covers a large region of which I am only interested in a small region. I have the co-ordinates of my region of interest xx°, yy° N, xx°,yy°E. I want to write a shell script which could extract only this region from each daily file and give me the output. Is there any way I can do this for all the 365*50 days. I would really appreciate any help.                     tar -xvf APHRO_MA_050deg_V1101.1951-2007.nc.gz.tar                  i=1951             while [ $i != 2008]     do        gunzip APHRO_MA_050deg_V1101.$i.nc.gz         r.in.gdal -o input=NETCDF:"APHRO_MA_050deg_V1101.$i.nc":precip output=prec$i         rm APHRO_MA_050deg_V1101.$i.nc     done      Thanks a lot


title: Access Shapefiles from PostgreSQL/Postgis in R

text: I successfully setup a connection to my PostgreSQL/PostGIS database on Ubuntu Server from my Laptop (Ubuntu 12.04). Therefore I installed RODBC in R to setup a connection and I made a correct entry in the /etc/odbc.ini file to setup the DSN as described for example here. Now I'm able to connect to my database with               ch<-odbcConnect("database")      After connecting, I'm am able query all the tables inside the database with SQL queries and even use the spatial functionalities of PostGIS inside the db. Now I would like to work with the spatial data inside R so I downloaded the RGDAL package. The driver for PostGIS should be automatically included in Unixbased systems so I execute               ogrDrivers()      to check this. It gives me back that PostgreSQL = TRUE so I assume this includes the PostGIS driver?! Otherwise it isn't listed. Now, in theory I should be able to read-in a layer like this:               test<-readOGR(ch,"state")      but unfortunately I gives back:               Error in ogrInfo(dsn = dsn, layer = layer, encoding = encoding, use_iconv = use_iconv) : Cannot open file      So did I miss anything? What can I do to make it work?


title: How do you create a CartoDB viz with a dynamic SQL statement?

text: I want to create a CartoDB viz where the where parameter in my SQL statement changes dynamically at runtime. When I share my CartoDB viz from the dashboard, I can't help but notice the the SQL statement is passed as a URL var. I thought that was an easy way to dynamically change the where clause. After playing around though, I found out that the parameter is ignored (and for good reason, to guard against SQL injection), though other parms like showing the title/desc and the ne/sw latlongs are processed. So I'm now trying to create a viz using the CartoDB API. But for some reason, I can't get the creating visualizations at runtime example to work.


title: Installing Geoserver on Ubuntu 13.04 - " trying to overwrite '/usr/bin/shp2pgsql-gui' "

text: I am trying to upgrade Opengeo Suite from 2.0 to 4.0, specifically for the upgrade of PostGis 2.0 to 2.1, I followed the instructions outlined on this page: http://suite.opengeo.org/opengeo-docs/installation/ubuntu/upgrade.html But had no success. I realized I hadn't fully installed opengeo-client or postgis 2.1, so I used apt-get install and received this:               The following packages have unmet dependencies:     opengeo-client : Depends: postgis-2.1 (>= 2.1.0) but it is not installed     postgresql-9.3-postgis-2.1 : Depends: postgis-2.1 but it is not installed     E: Unmet dependencies. Try using -f.      So I used -f and received..               The following extra packages will be installed:       postgis-2.1     The following NEW packages will be installed:       postgis-2.1     0 upgraded, 1 newly installed, 0 to remove and 3 not upgraded.     28 not fully installed or removed.     Need to get 0 B/701 kB of archives.     After this operation, 1,839 kB of additional disk space will be used.     Do you want to continue [Y/n]? Y     dpkg: warning: files list file for package 'postgresql-9.2-postgis2' missing; assuming                     package has no files currently installed     (Reading database ... 384786 files and directories currently installed.)     Unpacking postgis-2.1 (from .../postgis-2.1_2.1.0-1+opengeo_amd64.deb) ...     dpkg: error processing /var/cache/apt/archives/postgis-2.1_2.1.0-1+opengeo_amd64.deb (--unpack):      trying to overwrite '/usr/bin/shp2pgsql-gui', which is also in package postgis 1.5.3-2ubuntu1     dpkg-deb: error: subprocess paste was killed by signal (Broken pipe)     Errors were encountered while processing:      /var/cache/apt/archives/postgis-2.1_2.1.0-1+opengeo_amd64.deb     E: Sub-process /usr/bin/dpkg returned an error code (1)      I was recommended to go back and forth with **apt-get upgrade** and **apt-get -f install** to isolate the issue. But I am assuming that Postgis 1.5 is still in use, but I cannot find where to remove it. (If it is possible to do a complete wipe of the Opengeo suite then I can also do that if it would be simpler, all of my data is backed up) Please be as specific as possible in what to do, I'm not experienced in Linux. Most of my work is on the GIS/Web side. My Ubuntu release is: Ubuntu 13.04 Raring Let me know if I need to provide any more information (and how to provide it) and I'll be as hasteful as possible. Thank you.


title: QGIS 2.0.1: Changing layer properties with python console

text: I'm working on a plugin that will aim to ensure that the following steps either "transparent" to the user with one click :   1. Creating a table with predefined fields   2. Record in this table in postgis, then open in Qgis   3. Changing the properties (especially the properties of the fields tab) this layer to:                - Set "Edit widget" of fields       - Associate a form created with Qt Creator       For steps 1 and 2 is OK. However despite researches in the qgis' Cookbook and forums, I haven't find functions that allow me to carry out step 3 using the python console ... Have you any idea ? Thanks for the answers in advance! Guillaume


title: How to split a raster into multiple rasters using QGIS?

text: I would like to know how to split a raster into multiple rasters in a single operation using QGIS. Although I can clip the raster image, it would be better to split it using a single operation. I am using QGIS version 1.8.0.


title: How to build topology without an ArcEditor or ArcInfo license?

text: I have a huge dataset that I am capturing that was full of errors previously and I want to make sure that the new version is error free - which can be very time consuming with only an ArcView standard license. Therefore, I was wondering if anybody knew of any free/low cost (<£200/$350) alternatives that can be used to build topology without having to shell out for ArEditor? Many thanks Mark


title: ArcGIS JavaScript API Select Features

text: I'm having some trouble using the "selectFeatures" method with a feature layer. The code below builds a query using input text. If anyone has any pointers on what I could be doing wrong here, it'd be a huge help! I added a function to get some feedback after the selection is complete. It appears to be working fine when I test the search, but the results are not being displayed correctly.               //Query                     var query = new Query();                     searchButton.onclick=function(){                         var value = dom.byId('searchText').value;                         query.returnGeometry = true                          query.where = '"Type"' + " LIKE '%" + value + "%'"                         console.log('Query Registered');                         roundabouts.selectFeatures(query,FeatureLayer.SELECTION_NEW, function(selection){                                 console.log('Selection Query: ' + query.where);                                 console.log(selection.length + ' Results Returned');                         });                     };


title: STGeomFromText not working on large polygons

text: I'm using **Microsoft.SqlServer.Types** to do intersections on polygons.               SqlChars scFullSwath = new SqlChars(new SqlString(fullSwath));     SqlGeography fullSwathGeom = SqlGeography.STGeomFromText(scFullSwath, 4326); //Error out on big polygons      According to the documentation: > This **geography** data type method supports **FullGlobe** instances or > spatial instances that are larger than a hemisphere. This is an example of a valid polygon that **works** :               POLYGON((14.1977984372468 49.7086179545186,14.0975789966719 49.490931050796,12.5900898047007 45.9899553708373,11.2416180142584 42.4725018717541,10.0189326479206 38.9414519482468,8.89701484835269 35.3990096092695,7.85664665475729 31.8469172968655,6.88278698469561 28.2865986360566,6.79081552980762 27.9386602738175,13.5662986089156 26.3473709096527,13.6780927958767 26.689514921038,14.8759019115514 30.1864100426377,16.1814448369554 33.6666019429789,17.6150918852402 37.1268498599495,19.2025294573402 40.5631314224601,20.9765131252558 43.9703521754743,22.979305680181 47.3419272187238,23.1129487986833 47.550807797626,14.1977984372468 49.7086179545186))      ![enter image description here](http://i.stack.imgur.com/WLcxl.png) ![enter image description here](http://i.stack.imgur.com/W8NAF.png) ![enter image description here](http://i.stack.imgur.com/beney.jpg) As you can see it has correct **"ring orientation"** as mentioned by others on how to solve this problem. # Problematic Polygon:               POLYGON((14.1977984372468 49.7086179545186,14.0975789966719 49.490931050796,12.5900898047007 45.9899553708373,11.2416180142584 42.4725018717541,10.0189326479206 38.9414519482468,8.89701484835269 35.3990096092695,7.85664665475729 31.8469172968655,6.88278698469561 28.2865986360566,5.9634518930894 24.7192549336294,5.08892530812995 21.1459314523495,4.25119070563764 17.5675630126313,3.44351321077389 13.985005924225,2.66012537113857 10.3990599134859,1.89598509474054 6.81048335641945,1.1465837262999 3.22000326571774,0.407788590200467 -0.371677988700824,-0.324291851935639 -3.96387852745089,-1.05342955105917 -7.55593436365933,-1.78339318076703 -11.1471976352728,-2.51806681006348 -14.7370348965234,-3.26157909469963 -18.3248253652926,-4.01845238057626 -21.9099575448643,-4.79378397261009 -25.4918242201927,-5.59347616002234 -29.0698144501277,-6.42453876771587 -32.6433010407778,-7.29549969159991 -36.211621758752,-8.21697775205904 -39.774051024648,-9.20250398028347 -43.3297576050588,-10.2697317463825 -46.8777402574098,-11.4422726186086 -50.4167283166232,-12.7525722806525 -53.9450245293382,-14.2465807422689 -57.460247128251,-15.9916555981513 -60.9588899763672,-18.0905888714589 -64.4355326552956,-20.7079207471342 -67.8813344075985,-24.1225654627664 -71.2809433475941,-28.8408188355059 -74.6055692800629,-35.8463883696685 -77.7924382401655,-47.2063856590074 -80.7065857738326,-66.9785058803618 -83.0052691165418,-97.9750390100537 -83.9384177606558,-128.643556044951 -82.9375576334944,-148.00391894514 -80.6054307189309,-159.139108412103 -77.6779603190499,-166.029981462794 -74.4857156018706,-170.684840043845 -71.1595774540993,-174.061735811063 -67.7603847042112,-176.654974376848 -64.3161893679373,-178.737627762043 -60.8419881952857,179.528828146465 -57.3464295623487,178.043316526901 -53.8348168444905,176.739494352815 -50.3105766671909,175.572040277034 -46.7760308036764,174.508913762508 -43.2328270400433,173.526773241969 -39.6821924739992,172.608150515157 -36.1250880589781,171.739644782635 -32.5623067842119,170.910730829305 -28.9945370682842,170.112949405406 -25.4224048384334,169.339341980169 -21.846501773014,168.58404528169 -18.2674039867933,167.841992255889 -14.6856847908136,167.10868448838 -11.1019226955191,166.380012729174 -7.51670676722884,165.941432353483 -5.35533773778457,159.776665988585 -6.57902683369825,160.177433883186 -8.72380367556899,160.822079257621 -12.2824960195651,161.443414585553 -15.8408809387369,162.043430549844 -19.398296263682,162.623815450699 -22.9541494522971,163.185982448494 -26.5079178675654,163.73108551617 -30.0591498688561,164.26002344831 -33.6074656910582,164.773430063532 -37.1525586793143,165.271646700603 -40.6941960473808,165.754670150563 -44.2322194161788,166.222063966029 -47.7665451784619,166.672811789969 -51.2971640454233,167.105074008517 -54.8241408276084,167.515774078978 -58.3476136686714,167.899867192761 -61.8677941701833,168.248975202488 -65.3849687010856,168.548650592718 -68.899502591151,168.772353767069 -72.4118497117796,168.86640163402 -75.9225678964884,168.706491743022 -79.4152462878817,167.927467291569 -82.8961624959035,164.732805578133 -86.3800283841654,90.7183887585475 -89.5255796637184,0.869476212320378 -86.5107511675832,-2.53511419326111 -83.0263801264859,-3.35347936173547 -79.5441534398354,-3.52714539202692 -76.0510262182095,-3.43943131235962 -72.5383029542796,-3.21903289869946 -69.0235405384683,-2.9212128486995 -65.5061918147878,-2.57316966777165 -61.9858166317077,-2.18966563058102 -58.4620663215801,-1.77924490961034 -54.9346738854397,-1.34704681853976 -51.4034493155988,-0.896205077108871 -47.8682781762759,-0.428595677041346 -44.329121290399,0.0547408943611908 -40.786015213663,0.553351609541117 -37.2390723089532,1.06722155761446 -33.6884808125902,1.59668364420096 -30.1345044451656,2.14236556311945 -26.5774814836359,2.70516170545828 -23.0178236671284,3.28622299594447 -19.4560147963763,3.88696070360683 -15.8926094397698,4.50906229046143 -12.3282317178813,5.15451864742359 -8.76357458301829,5.82566311373916 -5.19940031270823,6.5252236911824 -1.63654187123703,7.25639066630586 1.92409340837588,8.02290309530139 5.48151764590984,8.82915889854057 9.03465271716534,9.6803551482695 12.5823158658826,10.5826675131422 16.1231989255446,11.5434812953803 19.6558394558106,12.5716910643249 23.1785804132512,13.6780927958767 26.689514921038,14.8759019115514 30.1864100426377,16.1814448369554 33.6666019429789,17.6150918852402 37.1268498599495,19.2025294573402 40.5631314224601,20.9765131252558 43.9703521754743,22.979305680181 47.3419272187238,23.1129487986833 47.550807797626,14.1977984372468 49.7086179545186))      ![enter image description here](http://i.stack.imgur.com/wCLpd.png) KML Version for google earth. When I give it this polygon that is longer than a hemisphere it fails and errors out: **Error:** > 24205: The specified input does not represent a valid geography instance > because it **exceeds a single hemisphere**. Each geography instance **must > fit inside a single hemisphere**. A common reason for this error is that a > polygon has the wrong ring orientation. # Documentation said it can handle **FullGlobe** so what's the problem? ![enter image description here](http://i.stack.imgur.com/b1tYv.png)


title: How should I manage PostGIS Raster data with different projections?

text: I have a requirement to store and manage archaeological geophysics data that is collected as a rectangular array of samples -- a raster image.   * Each raster will usually 20x20 or 30x30 floating-point samples, typically sampled at 1m intervals.   * A survey will consist of one or more of these images in a given location.   * It is possible that two different surveys may take place in different countries, or areas that use different projections, but each survey will use one and only one projection.    * They're never likely to be viewed together, each survey will usually sit by itself.   * The data will only be accessed by a custom front-end, so there will be no users getting direct control of it through `psql` or similar.   * Each sample needs to be stored as it was collected, so I can't reproject it into a common CRS such as Web Mercator because one sample could end up covering more or less area than in the original projection, and analysis will need to be performed on the data. How should I best store the data in a PostGIS Raster database? The options I have come up with are:   1. Ignore SRID constraints and store all the data in one table, writing my front-end code to deal with manipulating the data in a consistent manner.   2. Store all the data in one table, and rewrite the SRID constraint as a compound of SRID and survey ID.   3. Through table inheritance, create a new table for each new SRID.   4. Through table inheritance, create a new table for each survey. 1 and 2 break some of the nice automated parts of PostGIS, but will be otherwise hidden in front-end code. But queries will probably take slightly longer. 3 and 4 could end up with an explosion of tables that would make it harder to manage FK constraints and so on. Practically, the number of rasters per survey is anywhere from 1 to 100 or more, and the number of surveys is likely to run into the hundreds. But the number of distinct projections is likely to remain very low, which favours 3.


title: How much math does a GIS Analyst need to know?

text: For someone studying to pursue a career as a GIS Analyst, what math courses should he/she take? **Update** Here's a long list of free Math courses from MIT to serve as a frame of reference. Which are essential, useful, useless?


title: WFS service is not OGC compliant

text: I am quite surprised to discover several WFS services (v1.0.0) running on the web which are not compliant with the OGC WFS-1.0.0 standard. In the framework of my project, I use a JAXB binding built with the XSD schemas from the OGC, and notably the schemas for WFS v1.0.0. As you imagine, not OGC-compliant services are a problem in such a project. The problem of compliance is for example the XML element: > /WFS_Capabilities/ServiceIndentification which is replaced by: > /WFS_Capabilities/Service For example in the capabilities: http://www.disasterschartercatalog.org/ogc/cecec4wfs?REQUEST=GetCapabilities&VERSION=1.0.0&SERVICE=WFS Is anyone know the reasons of this? is it something common?


title: What reliable GeoCoding service would you recommend me?

text: I'm looking for a web service providing a full address based on coordinates (and reverse). I'm ready to pay for that service. I don't want to use Google's API.


title: OpenLayers - How to show UK Map with OSGB Projection

text: We are using OpenLayers for our GIS Application. We use Google and BingVE as source for the BaseMap. We use Projection in the OpenLayers with following code. We want to show UK Map with OSGB.               if(source=='bingve')     {         zoomlength=17-minZoomLevel;         map=new OpenLayers.Map('map',{                  projection:new OpenLayers.Projection('EPSG:900913'),             displayProjection:new OpenLayers.Projection('EPSG:27700'),             units:'m'              });     }      Does this will Automatically show Map in UK format OR we need to explicitly call Transform Method to convert Lag / Lon.


title: How to create a line in a given distance to an existing one?

text: I am looking for a possibility to create a second line in a given distance to an existing line. The existing line is curved, so simply create parallel line segments with CAD tools will not work. Any solutions? Thanks in advance. Addendum: I mapped vegetation types in the surroundings of small streams. The vegetation types have to be displayed on the map as stream-accompanying lines, prallel to the water bodies. So the second line (vegetation) has different properties compared to the first line (stream) and is a feauture of a different shape file. The solution I found yet is to create buffers around the streams and to trace their borders, but it is time-consuming.


title: random point through raster

text: I've some issue in some task pls help if you can as I've not good python skills. I've a raster with values from 5 to 97 and I want a python script which can scan the raster and generate random points (on some interval) based on raster values. Script/tool in QGIS will only generate points only where cell value is equal or greater than 40. ![enter image description here](http://i.stack.imgur.com/91Bax.png) If any script will work do share pls. or some other way of doing so. Regards


title: What GIS software is used by web maps sites (Google, Bing Maps, etc)?

text: We all have seen beatiful (imo) maps produced by Google and Microsoft. There are many others; some of them are of global use (such as MapQuest, TomTom or HERE), and some are just of country scale. I assume that these companies have obtained/collected/bought raw vector data and then have used a GIS software package and some custom cartography workflows to produce maps with good symbology/layout and then cached them into multiple tiles. Rasterizing vectors into tiles is not that particularly interesting, but I was often wondering what kind of Desktop GIS software have those large companies used to make those maps? Esri, for instance, uses ArcGIS platform for making the Esri basemaps such as Esri World Street Map. What desktop GIS software do Google and Microsoft use for creating maps? I am curious about that because I want to know what other GIS packages apart from Esri ArcGIS are capable of doing this work on such a large scale.


I have a potential project to utilize census data within a GIS and the hard part appears to be the consolidating of all the census data for the state of Florida from the below listed US Census website. Does anyone know if someone has already consolidated all of this data within a database like access or within a shapefile. Thanks, John Mitchell US Census Data: http://www2.census.gov/census_2010/0...ile_1/Florida/ US Census explanation of Data: http://www2.census.gov/census_2010/0...DME_SF1_v2.pdf


title: Integrate MarkerClusterer with Google Map V3

text: I'm having trouble getting the MarkerClusterer into my current Google Map (which has taken a long time to get this far!!). How can I combine the two? I'm using V3 of the api. Here's the example I'm working on http://www.mediwales.com/mapping Here's the MarkerClusterer code:               var center = new google.maps.LatLng(37.4419, -122.1419);     var options = {       'zoom': 13,       'center': center,       'mapTypeId': google.maps.MapTypeId.ROADMAP     };          var map = new google.maps.Map(document.getElementById("map"), options);          var markers = [];     for (var i = 0; i < 100; i++) {       var latLng = new google.maps.LatLng(data.photos[i].latitude,           data.photos[i].longitude);       var marker = new google.maps.Marker({'position': latLng});       markers.push(marker);     }     var markerCluster = new MarkerClusterer(map, markers);      Update: I've attempted to add the clusterer to my current code but it doesn't seem to work. `Places[i]` doesn't seem to feed into the clusterer.


title: Running focal statistic with iterator in ModelBuilder?

text: I need to extent the Value of a grid (water mirror elevation) within the nodata area until reaching a specific barrier (the limit of the watershed), in ArcGIS 10.2 for Desktop. I am using this:               Con(IsNull("grid"),FocalStatistics("grid",NbrRectangle(3,3),"MEAN"),"grid")      So far I have applied it successively to the last output, increasing the size of the rectangle slowly. This is run in a ModelBuilder. ![enter image description here](http://i.stack.imgur.com/MJBjY.jpg) I was getting bored and I have tried to run an raster iterator to apply a focal statistics on each Output. ![enter image description here](http://i.stack.imgur.com/ddomO.jpg) The model runs it only one time and stops. I was expecting to have to force to stop once it reachs the barrier . And the pre-existing value are also interpolated, which is not the purpose here. I would greatly appreciate any tips to improve this process. And nope, I still didn't have time to start learning python ;-) Thanks


title: Remember selected items on server round trip

text: Stack: GeoServer, GeoExt, PostGIS, etc etc (all the standard stuff). I have a map displayed in my browser and performing a box select adds relevant items to my selection vector layer. In addition to this a popup box appears with a summary of selected items and against each layer summary is a hyperlink which when clicked causes a server round trip (this round trip is necessary as it displays the main detail of each selected item). The server round trip results in the map being re-displayed but now my selected items are no longer selected. My question is, can I maintain the list of selected items and keep them highlighted after the map re-draw. I could probably cobble something together using cookies but was wondering if there is an easier solution!? Thank you!


title: Network Analysis: How to restrict certain road types and allow them to be traversable at intersections

text: I am faced with a problem for an assignment. I have tried restricting roads by using a restriction attribute with the following pre-logic script code with value = restricted               restricted = false     if [REST_STREET]= "Y" then     restricted = true     end if      The problem I am faced with when solving the closest facility is that these restricted roads do not allow contiguous roads to be traversable due to the restricted road intersecting with it. I have based the restriction on this question How to restrict portions of road network from the analysis in ArcGIS 10 Network Analyst?. Maybe I should just delete the roads I do not need?


title: Shapefiles for US city and Zip level

text: I already have the information for US states and Counties. I have the information in the Json files. I want the information for **ALL** US city and Zip Code information. Any leads appreciated. I cant really navigate the US consensus sites so if anyone has the information, I'd be really thankful.   * I found this site where I had to pay for Zip codes for different regions. Other places I have found the info for major cities i.e 10 odd cities only. I am looking for the entire US and would apprciate help.


title: Develop offline GIS program/annotator

text: I'm looking for a package to develop a simple offline GIS desktop application that works on a tablet/laptop for using it on the terrain (namely for annotating existing maps). The program has following requirements:   * Able to show/load vector and image data offline (formats doesn't matter much)   * Able to edit basic data. So able to save vector data.   * Desktop or Android   * We want to create our own routines for editing the data. I have looked a bit to ArcGIS SDKs, but offline is in beta there so I was wondering if there are other (maybe open source) toolkits that some of you have experience with. I'm not interested in things like buffering, geoprocessing,... We want to be able to show vector and image data and provide our interface on top of that.


title: I want to export from wikimapia into kml all sugar refinery category points

text: I want to select a country in wikimapia and then export all sugar refinery points into a kml file. Also how do I find the category id for sugar refinery?


title: How to convert WKT to SVG

text: I have some Well-known text (WKT) for representing geometry object such as Point, MultiPoint, LineString, Polygon, MultiPolygon etc. I have a multipolygon with in total 40000 Points. I have found this plugin which convert svg to wkt https://github.com/davidmcclure/svg-to-wkt Is there any plugin which covert wkt to svg?


title: How to create efficient database for routing?

text: I created a database of many important features and locations of my city, and I want to find shortest/fastest path between two locations. Now, I'm calculating the route by using A* algorithm in Java for android device, but it takes too much time to calculate and draw a simple route. Is there any efficient way to store and retrieve the road data from database for route calculation? Thank in advance.


title: Getting the name of the mxd using python

text: I have a python script which reads through all the mxd's in a folder and changes their properties from relative file paths to absolute file paths. I now want to add a bit to the script to change the version of the mxd from 10.1 to 10.0. I want the 10.0 mxd to be saved into a pre-defined folder. The new mxd should have the same name as the original mxd. This is the bit I'm struggling with. How do I get the name of the original mxd? (defined as basename in the code below). I'm assuming I need to use an arcpy.mapping function, but I'm completely stumped.               import arcpy, os     #workspace to search for MXDs     Workspace = r"C:\Testing\ImageMapperReplacement\LRM157"     arcpy.env.workspace = Workspace     #list map documents in folder     mxdList = arcpy.ListFiles("*.mxd")     #set relative path setting for each MXD in list.     for file in mxdList:     #set map document to change     filePath = os.path.join(Workspace, file)     mxd = arcpy.mapping.MapDocument(filePath)     #Get the file name     basename = ??????     #set relative paths property     mxd.relativePaths = False     #save map doucment change     mxd.saveACopy (r"C:\Testing\ImageMapperReplacement\LRM157\mxd2\\" + basename + ".mxd", {10.0})


title: Where has the cartogram plugin for QGIS gone?

text: This plugin used to be downloadable from Carson Farmer's plugin repository page. However, the repository is broken. An old version of the plugin source code is still available at https://code.google.com/p/ftools-qgis/ but using QGIS 1.8 the plugin can not be loaded ("corrupted"). Any help will be greatly appreciated...


title: A tool or way to Extend line by specified distance?

text: I have a purely aesthetic layer which has arrow symbols. Some aren't showing up right because the line is too small. I have selected maybe 50 records where I need to extend this line by a given number (ex. 2 meters). The extend line tool only extends lines to a specified intersection, so this tool isn't what I'm looking for. I've tried editing the shape length field but it won't let me. Is there a simple way of doing this via field calculator or within the editor tool bar? Thanks


title: How to avoid the faint display of tracks when using the time slider of Google Earth?

text: After playing a track with the Time Slider, I often get a faint but clearly visible track remaining. Anyone know why and how to fix it please? Here's an example. The GPX track has been played once and now I'm stuck with this faint remainder.   ![timeslider](http://i.stack.imgur.com/uqTFi.jpg) https://dl.dropboxusercontent.com/u/4019461/GE-TimeSliderFaintPuzzle.jpg


title: QGis 2.0.1 (german version): Transect with square polygons

text: I would like to create with QGis 2.0.1 ( Dufour , **german version** ) several square polygons (1km ²) that seamlessly connect to each other, so that a strip ( transect ) of 9km times 1km length results . ( I once attached a picture of how it should look ) -  http://www.myimg.de/?img=TransektArcMap75b51.jpg So far I know only the corresponding functions in ArcMap from ESRI. QGis seems equipped in the basic version with less features. Although I have the plugin " Rectangles ovals digitizing " installed, but can not choose here the edge length . In addition, the following square is to connect seamlessly . In ArcMap the cursor shows me when digitizing contacts with neighboring polygons or their vertexes or lines on , so you could work " seamlessly " . Furthermore, each of the squares , partially each quadrant , then an attribute (here, an air quality index) to be assigned , in accordance with which the square to be dyed. Are different colors within a shape based on various values (how it works in ArcMap ) possible? About tips on how I can tackle this project , I'd be very happy . PS: For information on tools , etc. ..... I use the German user interface. Best regards and thank you in advance Daniel


title: Save GM directions as KML - with timestamps?

text: In Google Maps I use Get Directions and after saving the result as a 'map' and opening it in My Places I am able to save it as a KML. All fine up to a point, allowing me to display the route in Google Earth. But I'd also like to **animate** it with GE's Time Slider. Given that the GM directions assumes speeds and provides an estimated journey time, is it possible to get GM to save the KML with timestamps please? Otherwise I must edit the 'raw' KML elsewhere and add arbitrary times myself, a relatively tedious task. -- Terry, East Grinstead, UK


title: MODIS HDF to Corrected Reflectance Bands 7-2-1. Change Reflectance

text: Convert HDF to Tiff easily. But I can not get the images to 7-2-1 bands. Like this : ![Ice on the Great Lakes. 26 February, 2014](http://i.stack.imgur.com/GcRCN.jpg) How can I get this image from HDF and what product should I use (MYD09GQ and etc)? Maybe using a ArcGIS? Thanks


title: Is it possible to get all projections form OGR?

text: I'm working on a project where the user can give my program so GIS data which may or may not have which projection it is using. So is it possible to get a list of all the projections from OGR or Proj4? I am quite new to using this GIS stuff so I'm sorry if I'm using the wrong terms. Thanks for all the help. **Edit** The program I am working will take a GIS file from the user as an input and will show all the points. But sometimes the files do not say which projection they are using so I dont know what to do with the data. Instead of just throwing an error message I am trying to show a tree view of all possible projection the user can choose from. Kind of like ArcGIS, where you can change which projection is being used. Also I'm using c++ for this project.


title: No placement inside feature dataset using arcpy XYtoLine?

text: When I run the following code with the outFC creation path set to a feature dataset inside a file geodatabase...               arcpy.XYToLine_management(tempFC, outFC, "POINT_X", "POINT_Y", "NEAR_X", "NEAR_Y",     "GEODESIC", "", formatSpatialRef)      ...I get the following message:               ERROR 000582: Error occurred during execution.      The script runs without errors when I set the outFC creation path to the root level of the gdb. Spatial references for the feature dataset and the output FC are identical. I'm wondering if this is a case of a built-in behavior similar to what is described in this posting: http://forums.arcgis.com/threads/39448-Using-Python-can-t-get-Batch-Project- outputs-to-save-in-a-feature-dataset


title: Is there anything obviously wrong with this hit test code?

text: I want to know if the IPoint of the mouse position is hitting an IPolyline, I wrote this code, and was optimistic that it would work, but it always returns false.                 IPolyline poly1 = pc1 as IPolyline;                  IHitTest test = poly1 as IHitTest;             if (test != null)             {                 int a = 0;                 int b = 0;                 bool c = false;                 double distance = 0;                 IPoint hitPoint = new PointClass();                 bool isHit = test.HitTest(mousepoint,                                         50000,                                         esriGeometryHitPartType.esriGeometryPartNone,                                         hitPoint,                                         ref distance,                                         ref a, ref b, ref c);                       }


title: What are OpenStreetMap license requirements?

text: We are planning to mix OpenStreetMap (OSM) layers with our own layers and afterwards to build an application which will perform common GIS tasks. Assuming we are allowed to mix our layer and OSM layers in one application, where is the most appropriate place to mention OSM copyright info? Perhaps including this info in the layer metadata is enough or maybe we have to include this info in the application as well.


title: How do I use ranges for raster colour map in QGIS

text: I have a raster layer in GeoTIFF format, with continuous data values between -1 and +1. For the colour map, I would like...   * >0 to be green   * = 0 to be grey    * <0 to -0.25 to be pink    * <-0.25 to -0.50 to be dark pink ...etc to 1. Is this possible with QGIS? I tried putting ranges in the colour map box, but it didn't work. If not, I'll do it in R.


title: ArcMap Find Tool Customization

text: I am looking to create a customized version of the find tool where the layer and field tools are pre-configured. Ideally I could add this to ArcMap 10.1 as a button for end users to who use same search using the find button repeatedly. Any help would be appreciated.


title: Could not create datasource

text: I'm fighting with an error I can't figure out where it is coming from. Tried to compile all the things and used the latest unstable/stable ppa repositories as well, can't find the issue where it is comeing from. Servers Setup is:   * Postgres 9.1    * Postgis 2.0 (2.1)   * Mapnik 2.0    * TileCache latest An error occurred: Could not create datasource. No plugin found for type ' postgis ' (searched in: /usr/local/lib/mapnik/input) (encountered during parsing of layer 'landuse_gen0' in map '/var/mapnik/osm.xml') File "/usr/local/lib/python2.7/dist- packages/TileCache-2.11-py2.7.egg/TileCache/Service.py", line 343, in cgiHandler format, image = service.dispatchRequest( params, path_info, req_method, host ) File "/usr/local/lib/python2.7/dist- packages/TileCache-2.11-py2.7.egg/TileCache/Service.py", line 208, in dispatchRequest return self.renderTile(tile, params.has_key('FORCE')) File "/usr/local/lib/python2.7/dist- packages/TileCache-2.11-py2.7.egg/TileCache/Service.py", line 138, in renderTile data = layer.render(tile, force=force) File "/usr/local/lib/python2.7/dist- packages/TileCache-2.11-py2.7.egg/TileCache/Layer.py", line 436, in render image = self.renderMetaTile(metatile, tile) File "/usr/local/lib/python2.7/dist- packages/TileCache-2.11-py2.7.egg/TileCache/Layer.py", line 396, in renderMetaTile data = self.renderTile(metatile) File "/usr/local/lib/python2.7/dist- packages/TileCache-2.11-py2.7.egg/TileCache/Layers/Mapnik.py", line 39, in renderTile mapnik.load_map(m,self.mapfile)


title: GDAL HDF not recognized as supported file

text: I am having the same issues as in this thread; however, I'm new to programming and need more advice on how to get gdal to open hdf files. I downloaded hdf from MODIS, which I believe are hdf 4. I have a mac and would like to work with the data in python. I installed HDF4.2r4 and szip-2.1 with homebrew. I know hdf is not available, in terminal I ran:               gdalinfo --format hdf4      which gave:               ERROR 1: --format option given with format 'hdf4', but that format not     recognised.  Use the --formats option to get a list of available formats,     and use the short code (ie. GTiff or HFA) as the format identifier.      Further when I attempt to open the file:               gdalinfo MOD04_L2.A2003001.1755.051.2010313012055.hdf     ERROR 4: `MOD04_L2.A2003001.1755.051.2010313012055.hdf' not recognised as a supported   file format      I would greatly appreciate any help on how to get gdal to open my MODIS hdf files. Please let me know how can can clarify this question. Many thanks.


title: Is there any GRASS command to measure the distance between a point in rastermap and another in vector map?

text: I have a question about WinGRASS 6.4.3. Is there any GRASS script command to measure the distance between a point in rastermap and another in vector map? I have a raster map and a vector map. Raster is mapped fire points, and vector is mapped houses as point data in same area. Now I can see the raster map superimposed with the vector map on the Map View. I would like to know which fire points are most close to house. So, I would like to know distance between each fire and house. I know there is GUI tool to measure distance, however, I have to measure many distances for many maps. So, I need to know the GRASS command to measure distance.


title: Gamma correcting ECW

text: we would like to modify an existing ECW file with a gamma correction in order to make the image less colourful. Is there a way to do this using e.g. GDAL?


title: Projections when using google maps and rest service on openlayers

text: I am trying to use a rest service on an openlayers map that has 2 google maps layers. The google map layer is set to it's own custom projection (which I had to code into the library by hand, but I believe is the standard EPSG of 900913).. However the layer that I would like to sit on the google map is in Oklahoma State Plane North.               var arcrest = new OpenLayers.Layer.ArcGIS93Rest("Rest Layer", arcserver_rest, {                                 LAYERS: 'show:18,19,20,21,29,30,31',                                 IMAGESR: '2267',                                 BBOXSR: '2267',                                 FORMAT: 'png',                                 TRANSPARENT: true                             },                             {                             isBaseLayer: false,                             units: 'ft',                             maxExtent: new OpenLayers.Bounds(2034036.00510709,  136689.343284666, 2069119.14841628, 148575.857284666)                             }                             );      However, this is projecting into the Congo, and not Oklahoma. Any suggestions.. Thank you.


title: projection conversion from mercator to google maps format

text: Is this               spTransform(Finland3, CRS("+proj=longlat +ellps=WGS84 +no_defs+ellps=WGS84 +towgs84=0,0,0 "))      the right way to project a shapefile with this projection               coord. ref. : +proj=merc +lon_0=0 +lat_ts=0 +x_0=0 +y_0=0 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0       to plot it in `ggmap` on top of a Google map raster? therefore looking like               +proj=longlat +ellps=WGS84 +towgs84=0,0,0


title: Is there a way to use some sort of "instancing" when storing polygons in a spatial database?

text: I am quite new to GIS database addons, and I am still exploring the possible uses for them. In general the idea of storing objects indexed by their position which enables additional operations is pretty cool. However, it's fit for data gathered from real world scenarios, where every park, building etc. will have its own set of points. The simplest way is to simply store them as points and execute regular queries on them. However, imagine we're putting some imaginary objects on the map, belonging to a few categories, but in large quantities. Let's take a star shape for example. Now, I can simply add the set of polygons every time I want to place my star in my "world", but this creates a huge duplication of data and requires me to modify each and every record shall I want to update the shape somehow. The obvious solution to that is to store only the position of said shape, and a key to some other table holding shapes with 0-based coordinates. I can perform things like intersection tests by using `ST_Translate` on the shape and its position; so far so good. Will this approach scale to a bigger number of instances? Will I be able to perform every query I would be able to using regular polygon storage? And lastly, but not leastly, doesn't PostGIS provide a built-in support for that? I was unable to find that in the documentation.


title: What does pgr_createTopology() expect for a table parameter?

text: In this sample http://docs.pgrouting.org/2.0/en/doc/src/developer/sampledata.html#sampledata there is this function:               SELECT pgr_createTopology('edge_table',0.001);      I understand that this creates a topology from edge_table. My question is what is inside the edge_table? To be more specific: my table has ID, latitude, longitude, geom(POINT), and the_geom(LINESTRING) fields, with ten records. This is the client address. Then I transform latitude and longitude into geom(POINT) with PostGIS functions. Then I use st_makeline() to make a line between one point (client) to another, and update the the_geom(LINESTRING) field. One client connects just one time with the other (just the last one returns to the first one, because it needs a target to make a line). That gives me one network. Then I use pgr_createTopology to generate a vertex table, then use the pgr_dijkstra() function, but it returns me the same path that I built. What am I doing wrong before executing the pgr_createTopology? What do I have to put into my table to make this function work? I don't have any map loaded on my table, just points from latitude and longitude. I appreciate any help.


title: How do I edit an existing PMF file using ArcGIS Publisher?

text: I downloaded a trial version of ArcGIS Publisher and the second question on this page (http://www.esri.com/software/arcgis/arcreader/common-questions) suggests that editing is possible in addition to creating PMFs. Ultimately I would like to add a layer to an existing PMF and republish the whole thing.


title: How to calculate weighted mean coordinates within polygons?

text: I'm using QGIS 1.8 on OSX 10.8.3. I have a polygon shapefile of all of the US state legislative districts. I also have a layer of points that are centroids of US Census Tracts, with a field for the population of the tract. How do I generate a layer of points that are the weighted mean of each census tract point within each legislative district? My hunch is to play with spatial joins and use fTools to do mean coordinates, but I'm not sure how to calculate the mean for each polygon instead of the whole layer itself. Any thoughts?


title: How do you create a legend for the bucket divisions in a fusion table map?

text: I am creating a map using Google fusion tables. I have divided the data into "buckets" using the Configure Map Styles window. Is there a way to add a legend to the map telling the readers what the colors (buckets) relate to?


title: Problem with setValue() function in ArcGIS 10

text: I'm trying to add data to an attribute table in ArcGIS 10 using the following code:               def make_floor_no( shapefile ):         "Makes header for number of Floors (FLO) and calculates value"              fieldName = "FLO"         try:             ARCPY.AddField_management(shapefile, fieldName, "DOUBLE")         except:             print "Field already there"                # loop through attribute table             Rows = ARCPY.SearchCursor( shapefile )               for row in Rows:             floors = round( row.getValue( 'HGT' ) / 3.0, 0)             print str(floors)             row.setValue( fieldName, floors )             Rows.updateRow(row)                  However, I keep getting an error along the line               row.setValue( fieldName, floors )      I can't spot anything wrong with this, and have tried a few different options. Is there something wrong with my syntax? The error message is as follows:                File "Z:\ConstructionMaterial.py", line 100, in <module> make_floor_no( InputFile )      File "Z:\ConstructionMaterial.py", line 71, in make_floor_no          row.setValue( fieldName, floors )      File "C:\Program Files (x86)\ArcGIS\Desktop10.0\arcpy\arcpy\arcobjects\arcobjects.py", line 941, in setValue         return convertArcObjectToPythonObject(self._arc_object.SetValue(*gp_fixargs(args)))     RuntimeError: ERROR 999999: Error executing function.


title: What is ArcGISCacheMgr.exe and why does it have a lock on my folder?

text: I've got a 10.0 gp service running on Windows Server 2008. The gp service writes to a folder on the server. The parent of this folder is shared. When it comes time to delete folders on the server, occassionally I find that it is locked by ArcGISCacheMgr.exe. Does anyone know what this process does, or how I can get it to let go of my folder? **Update:** The ArcGISCacheMgr.exe process is running under the ArcSOC user account, so this makes me think it has something to do with ArcGIS server. (The folder is created by the gp service.)


title: Coordinates format in ArcGis 9.3

text: Regarding to Jay suggestion i open a new topic. Do you know how i can set in arcgis 9.3 coordinates only with one locus after comma, i mean like that format: xx xx xx,x not like i have now xx xx xx,xxx (screen below) ![alt text](http://i.stack.imgur.com/PGuxY.jpg)


title: Fusion table styling

text: I am importing data from fusion table data as a layer. but the colour i used in fusion tables is not reflected in my map. it by default remains red small icon. Please let me know if I am doing something wrong. layer1 = new google.maps.FusionTablesLayer({ query: { select: 'geometry', from: '1f1XfTx7J52PZJv3mwnBq3rCiPwb6x1li67ca6lg' } });


title: Derive WGS84 Longitude & Latitude from British Northings & Eastings

text: I'm looking for a fairly simple way to append columns with WGS84 longitude and latitude to a table containing northings and eastings. The source data is Ordnance Survey's CodePoint Open, which I is I think OSBGB36. I want to use this with OpenStreetMap. I have the CodePoint data in a csv file and loaded into an MS Access database, from which I can load it into QGIS. From there I've tried exporting it with a transform operation into various file formats but this doesn't seem to do the trick. I do have access to PostGIS and SQL Server 08R2 but little experience of using sql spatial.


title: Why is my longitude conversion from NY State Plane Coordinate System to WGS84 consistently off by 3.665 degrees?

text: I'm trying to convert coordinates from the NY State Plane Coordinate System for Long Island (zone #3104) to WGS84 using the GeoTools Java API. The latitude value is correct, but the longitude is consistently 3.665 degrees off to the west:               public static void main(String[] args) throws FactoryException, TransformException {             CoordinateReferenceSystem sourceCRS = CRS.decode("EPSG:4456"); //http://www.geotoolkit.org/modules/referencing/supported-codes.html             CoordinateReferenceSystem targetCRS = CRS.decode("EPSG:4326");             MathTransform tr = CRS.findMathTransform(sourceCRS, targetCRS);             DirectPosition sourcePt = new GeneralDirectPosition(982394.0, 199298.0);              DirectPosition targetPt = tr.transform(sourcePt, null);             System.out.println("Source point: " + sourcePt);             System.out.println("Target point: " + targetPt);      }      The code above, earthpoint.us, and Google all return similar values for latitude: ~40.714. For longitude, earthpoint.us and Google show values of ~-74.006, but the code above produces -77.671, putting the point in the middle of Pennsylvania. Using different source NY SPCS values produces similar errors, with the difference in longitude always being 3.665 degrees to the west and latitude being correct. GeoToolkit produces the same results. Am I doing something wrong? Is there a 3.665 constant I should know about? EDIT: the source values are NAD27.


title: Can I mix versions of ArcGIS server and geocoding locator styles?

text: I have ArcGIS server running and a navteq map installed as a map service. I'm trying to add a geocoding service to the server so I can do reverse geocoding. My map came with a few locator files and I was able to create the geocoding service with these files but when I try to start the service, it tells me `Unable to Start service. Error Server Object instance creation failed Participating locator PointAddress [USA_PointAddress]`. The "verbose" log doesn't give me any additional information and I'm at a loss for what to do. The only thing I can possibly think of that might be the issue is that the locator is set as `Style release: 10.1` and the ArcGIS server is version 10.0. Would this be the issue?


title: How to connect to ArcGIS Server on the other machine?

text: There is the ArcGis Server on the machine 192.168.**.22 Machine name is MAPSERVER. User name and password for ArcgIS Server are "test" and "testtest" (for example) ArcGis Server options are: clusters - name "default", computer "MAPSEVER", protocol "TCP port 4004". I want to connect to this ArcGis Server. I used code               ESRI.ArcGIS.ADF.Identity identity = new ESRI.ArcGIS.ADF.Identity("test", "testtest", "default");     ESRI.ArcGIS.ADF.Connection.AGS.AGSServerConnection agsconnection;     agsconnection = new ESRI.ArcGIS.ADF.Connection.AGS.AGSServerConnection("MAPSERVER", identity);     agsconnection.Connect();     IServerObjectManager SOM = agsconnection.ServerObjectManager;       I get error ` > "Identity error"


title: Which projection: Map from South Tyrol (Italy)

text: I need to find out what projection this base-map has for digitizing the information that was drawn on it. It's from South Tyrol (South Tirol) in Italy. I have the following information: Title: **Provincia di Bolzano / Provinz Bozen** Year: It was published **1980** , so it could be anything before that It has this overview: ![enter image description here](http://i.stack.imgur.com/1vRbw.jpg) Next to it is has this code: **4 III NO-B** And the page number: **Foglio 74 / Blatt 74** The x-coordinates are given like this: **E=1650000** And the y-coordinates: **N=5183000** At the bottom of the map it says (some sort of copyright notice):                 Questo supporto consente solo copie a luce fredda            Cartogr. Autorizz. A norma di legge      And then it says:                 E.I.R.A.            Ente Italiano rilievi aerofotogrammetrici            Rilievo aerofotogrammetrico eseguito con strumenti "Galileo-Santoni"      Edit: I did some research and it might be a Gauß-Boaga (Gauss-Boaga) or Roma 40 projection, but I'm not sure yet. Link: https://it.wikipedia.org/wiki/Proiezione_di_Gauss-Boaga


title: Good place for starters to develop ArcGIS map

text: I am new to ArcGis 10.2. I want to develop a java application which displays map. Please guide me , where to start looking thanks


title: ListFeature classes in SDE

text: I have this code:               def backSDE(ActiveSDE):         arcpy.env.overwriteOutput = True         arcpy.env.workspace = ActiveSDE         fcs = arcpy.ListFeatureClasses()         print fcs          #Databases     ActiveSDE = r"C:\Users\me\AppData\Roaming\ESRI\Desktop10.1\ArcCatalog\SERVER_SQLEXPRESS.gds\Client_Survey"     backSDE(ActiveSDE)      The database only has feature classes inside of it. No datasets. When I run this I get "None" in the python interpreter. There should be 10 feature classes. What am I doing wrong? I'm running ArcMap 10.1 * * * **UPDATED (2)** I now have this code               def listFcsInGDB():         ''' set your arcpy.env.workspace to a gdb before calling '''         for fds in arcpy.ListDatasets('','feature') + ['']:             for fc in arcpy.ListFeatureClasses('','',fds):                 yield os.path.join(arcpy.env.workspace, fds, fc)          #Set workspace environment to geodatabase     arcpy.env.workspace = r"C:\Users\me\AppData\Roaming\ESRI\Desktop10.1\ArcCatalog\PETROPAM_SQLEXPRESS.gds\Client_Surve"          #Get list of feature classes in geodatabase     #FCs = arcpy.ListFeatureClasses()     FCs = listFcsInGDB()           #Loop through feature classes in list     for FC in FCs:               #Print the feature class name          arcpy.AddMessage(FC)      This is the error message it gives me on this line from above: for fds in arcpy.ListDatasets('','feature') + ['']: * * * **UPDATED (3):** Here is my updated code that lists the feature classes from the database connection file.               import arcpy     import os          arcpy.env.workspace = r"C:\Users\me\AppData\Roaming\ESRI\Desktop10.1\ArcCatalog\PFS_Survey.sde"     def listFcsInGDB():         ''' set your arcpy.env.workspace to a gdb before calling '''         for fds in arcpy.ListDatasets('','feature') + ['']:             for fc in arcpy.ListFeatureClasses('','',fds):                 yield os.path.join(arcpy.env.workspace, fds, fc)          #Set workspace environment to geodatabase     arcpy.env.workspace = r"C:\Users\me\AppData\Roaming\ESRI\Desktop10.1\ArcCatalog\PFS_Survey.sde"          #Get list of feature classes in geodatabase     #FCs = arcpy.ListFeatureClasses()     FCs = listFcsInGDB()           #Loop through feature classes in list     for FC in FCs:               #Print the feature class name          arcpy.AddMessage(FC)          print FC          arcpy.CopyFeatures_management(FC, "C:/Users/me/Desktop/" + os.path.basename(FC).split('.')[0] + ".shp")      It outputs this:               C:\Users\me\AppData\Roaming\ESRI\Desktop10.1\ArcCatalog\Survey.sde\Client_Survey.DBO.CO_Arapahoe     C:\Users\me\AppData\Roaming\ESRI\Desktop10.1\ArcCatalog\Survey.sde\Client_Survey.DBO.CO_Montezuma     C:\Users\me\AppData\Roaming\ESRI\Desktop10.1\ArcCatalog\Survey.sde\Client_Survey.DBO.CO_Weld     C:\Users\me\AppData\Roaming\ESRI\Desktop10.1\ArcCatalog\Survey.sde\Client_Survey.DBO.Wy_Converse      The problem now is trying to copy it out of the list. I've updated my code to how I want to copy it, but the ListFeature classes has two decimals in the code. **UPDATE 4** I'm starting a new question that is somewhat based off of this one, but different enough to be in a new page. I will link once it's up. Here's the link to the question about copying several sde connections.   Backup several SDE database connections


title: Flow modeling and FOSS GIS; lahars and HEC-RAS

text: **What are the options for FOSS flow modeling of dilute debris flows in QGIS or GRASS?   Have you, or can you refer me to studies where the authors, tested these models/plugins/functions/etc. for use in lahar modeling?**   (Lahars are muddy volcanic flows, I'm interested in ones that will have a very small solids fraction of <10-20%, making some flood models applicable to the study; other add'l info: my field site has poor DEM coverage and is highly channelized, aka not a flood plain) I am working on developing a FOSS workflow for disaster risk reduction, specifically for using DEMs in FOSS GIS that can either be used as input in separate flow models or with flow models that run in GRASS or QGIS. I know of a few functions that can do this (r.inund.fluv, r.hazard.flood, r.drain, etc) but they seem to be untested for this function. I was interested in using HEC-RAS but my collaborator feels that it only runs effectively in Arc as HEC-geoRAS. The other model I was interested in (again, Arc) is the D8-algo based LAHARZ model, which (it seems?) is kind of similar to r.drain, at least in theory. Any information at all on the matter is greatly appreciated!


title: For Better response time-Map loading

text: I am trying to develop an application using OL. In my code I am using base maps of google , yahoo etc. For this I am including the link in the script,this means that I am connecting to respective server(google,yahoo) etc every time through internet connection(I think this is what including the links in scripts mean). But sometimes when I select Yahoo maps , it takes a lot of time to load and in some cases it doesn't load at all. So for a better response time my idea is to use caching concept and store these maps in my memory using cookies etc.... But I donno whether Its possible to do this or not.So will this caching concept work? Are there any other better ways to get good response time with loading of maps???? or How can I load maps better? Your suggestion are much valuable.


title: How to combine data from overlapping polygons?

text: Here's what I'm working with:   * Shape file of the current 32 Congressional districts for all states   * Shape file of the tract boundaries for Texas   * A DBF format table of Hispanic counts and proportions for Texas census tracts (joined to the tract boundaries) Using this data set, I'm trying to determine the following based upon the centers of the census tracts:   1. Total Hispanic population for each Texas Congressional district   2. Hispanic proportion of total population for each Texas Congressional district   3. Count of majority Hispanic census tracts for each Texas Congressional district   4. The shape index value for Texas Congressional district I'm also trying to create a map applying shading symbology to display the number of majority Hispanic census tracts for each Texas Congressional district. My problem arises when I try to combine the attribute data from the census tracts with the congressional districts. I need to combine all of the tracts within individual congressional districts and sum up their Hispanic populations. Dissolve isn't working, because the congressional districts and census tracts don’t share any common attributes. When I join the shape files spatially, most of the attribute data gets lost. If anyone knows how I could go about doing this, your help is appreciated!


title: How to get 'IGraphicsContainer' from the map?

text: I tried to get `IGraphicsContainer` from Data View (active map) using `FocusMap`, `ActiveView`               IGraphicsContainer graphicsContainer = (IGraphicsContainer)neApp.Document.MxDocument.FocusMap;      but the variable always null, is there any different way to do this? * * * I already get from page layout view in this question right here: Get IGraphicsContainer from page layout but i need it from data view (map).


title: How to get empty gdb using python

text: I have one gdb with lots of layer (~500) and to receive empty gdb I delete all records from layers. I Use deleterows / TruncateTable but for one layers take abaut 1-3 minuts. Is better way (faster) to get a empty gdb without records usnig python?


title: Activating and deactivating of button control not working properly

text: I have got two buttons, one for editing purposes and one just to see all the information related to the polygon I selected. Also I added these two buttons on my panel using OpenLayers control. I am working on a single layer. Now, my problem is that the first time I click on one of the buttons it works properly : a popup comes out properly and it also deactivates properly when I click on the button again. However, when I click a second time, it does not work anymore : the popup does not destroy or get cancelled when I click on different features, all the popup remain there. I tried lots of ways but I still couldn't solve the problem. Any help would be appreciated.   Thanks in advance. Code is posted below.                  <!DOCTYPE html>               <meta charset='utf-8' />     <title>My OpenLayers Map</title>


title: What are the options for publishing maps online, with support for high-res aerial imagery?

text: The place I am working wants to publish their maps online - however the only system they have is mapshare, which isn't supported by ESRI anymore. I've suggested DEKHO but I think they are looking at not spending much money. Atm we are going with layered pdfs however I can't export any map with aerial imagery as its 15cm resolution and the pdfs end up blank or only with part of the imagery on it. Could anyone suggest any other options? I was thinking along the lines of a python application but I am not that good with python and I wouldn't know how to do things server side. The issue I have with layered pdfs is that it doesn't allow people to choose their own extents.


title: Extracting Earthquake Depth from the USGS GeoJson Earthquake Feed for use with Openlayers

text: I have some GeoJson from the USGS as follows:               {"type":"Feature","properties":{"mag":4.3,"place":"30km ESE of Jarm, Afghanistan","time":1382069641720,"updated":1382071013039,"tz":270,"url":"http://earthquake.usgs.gov/earthquakes/eventpage/usb000kg6q","detail":"http://earthquake.usgs.gov/earthquakes/feed/v1.0/detail/usb000kg6q.geojson","felt":0,"cdi":1,"mmi":null,"alert":null,"status":"reviewed","tsunami":null,"sig":284,"net":"us","code":"b000kg6q","ids":",usb000kg6q,","sources":",us,","types":",cap,dyfi,general-link,geoserve,nearby-cities,origin,phase-data,tectonic-summary,","nst":null,"dmin":2.372,"rms":0.49,"gap":112,"magType":"mb","type":"earthquake","title":"M 4.3 - 30km ESE of Jarm, Afghanistan"},"geometry":{"type":"Point","coordinates":[71.1215,36.7111,215.63]},"id":"usb000kg6q"},      The geometry contains 3 dimensional data, where it is necessary to set the GeoJSON.ignoreExtraDims parameter to true to ignore that third dimension to load the Geojson into a Vector Layer This third dimension appears to be the depth of the Earthquake! How can i get access to this if it does not appear in the data in the Vector Layer (It does not, I have checked, it is just ignored)


title: How to calculate the optimal zoom-level to display two or more points on a map

text: We want to display several markers on a static map and want to calculate the optimal zoom-level like Google Maps does. We already calculated the bounding rectangle and the centerpoint of the map but now we have a hard time to calculate the correct zoomlevel to display the whole bounding rectangle. Can someone please point us into the right direction?


title: MapInfo "Error during Read: XML Library Load Error Required XML version unavailable"

text: I am using MapExtreme 7.0, I am running into the following error on an asp.net website (.net 3.5) running on a server with windows 2003 32bit. "Error during Read: XML Library Load Error Required XML version unavailable" when trying to access some feature data. I managed to track down some obscure reference in the 6.7.1 release notes and the fix was to install the latest "XML", I'm assuming that refers to MSXML, I have checked and the server has MSXML 6 installed. Is there something else this may be referring to? Or is there another cause of the problem. I managed to reflect the MapInfo dlls and have found the resource that contains the string but I cannot find where the error actually originates from


title: Filtered SimpleFeatureCollection really slow on Map and other operations

text: I'm working with a shapefile having about 85.000 features. I need to select some of these using a filter and then showing these on a map. Now... If I draw the entire shapefile's features the creation of the map, the zooming and moving around is fast but if I select, even only 100, features from the SimpleFeatureSource, these steps are really slow. Here is how I select, for example, the first 1000 features and then draw: The execution stuks at Layer layer = new FeatureLayer (features, style); and at JMapFrame.showMap(map);               public static void main(String[] args) {         File file = JFileDataStoreChooser.showOpenFile("shp", null);          if (file == null) return;      try{                  FileDataStore store = FileDataStoreFinder.getDataStore(file);             SimpleFeatureSource featureSource = store.getFeatureSource();                  FilterFactory2 ff = CommonFactoryFinder.getFilterFactory2();                  String[] selection = new String[1000];                  for(int i=0; i<selection.length; i++){                      selection[i] = i+"";                  }                  ArrayList<Filter> match = new ArrayList<Filter>();                  for (String name : selection) {                 Filter aMatch = ff.equal(ff.property("ID"), ff.literal(name), false);                 match.add(aMatch);             }             Filter filter = ff.or(match);                  SimpleFeatureCollection features = featureSource.getFeatures(filter);                            // Create a map content and add our shapefile to it             MapContent map = new MapContent();             map.setTitle("Quickstart");                       Style style = SLD.createSimpleStyle(features.getSchema());                       Layer layer = new FeatureLayer(features, style);                  map.addLayer(layer);                  // Now display the map             JMapFrame.showMap(map);                  }catch (Exception e)             {                 e.printStackTrace();     }      Are there some faster way to do this? ......


title: How to Select By Location with matching attribute?

text: I have a whole bunch of overlapping polygons (buffers), I want to select with points only those poly's which:   1. the point falls within; and   2. the poly attribute matches the point attribute What would be the best way to do this selection, perhaps with ModelBuilder? (using ArcGIS Desktop 9.3.1)


title: How to apply a color gradient to a shapefile from QGIS console?

text: Anyone knows how to do a color gradient from the console to a shape file in Qgis (orin Python)? Like the image below. ![enter image description here](http://i.stack.imgur.com/km9PU.png)


title: Looking for free Mexican transportation datasets

text: I have been looking for data of Mexico. What I am looking for is the data of a Transportation in Mexico, specifically Chihuahua province. I have checked out on ArcGIS online and Google but came up nothing good. Is there one that I overlooked? I am looking for free gis data. Thank you.


title: Free high-resolution Orthophotography for Quebec

text: Is anyone aware of free high-resolution orthophotography for Quebec? Something 1m resolution or greater would be ideal. I'm looking for imagery to use in a standard GIS as files (so not as WMS or other service, or via the web.)


title: ArcGIS Server Manager keeps getting COM errors

text: I am using ArcGIS Server 10 on IIS6. While publishing out mapping services/restarting services/deleting services, I keep getting various COM errors like the one below..after viewing the Event log it looks like the ServiceWizard.aspx file is generally the root cause. However, usually if I wait thirty minutes or go to the restroom or other seemingly non-related activity and then try again it that I was trying to do will work. The errors are a little to scattered to describe fully. The question is..Is there something I may be doing wrong or should I restart the SOM every few days or another option? FYI, there are a lot of mapping services on this instance. Thanks. "An error has occured on the server. For details please check the Event (Application) log on the web server. Error HRESULT E_FAIL has been returned from a call to a COM component."


title: How to run two PostGIS versions in parallel on Windows?

text: I'm trying to understand how to install PostGIS2.0 on Windows7. The PostGIS site has the following suggestions: > PostGIS Experimental Binaries PostGIS 2.0.0SVN Binaries with Raster, > Topology support + Extras (Tiger Geocoder 2010) As of November 25th, 2011 > (r8242) includes GEOS 3.3.2dev r3516, Proj 4.6.1, JSON-c 0.9, libgdal 1.9 > SVN r23320 (VERSIONS prior to (r7923) require a dump reload since operators > have changed. We now have the 3D spatial indexes, 3D box operators, KNN GIST > for 9.1, and 2 raster map algebra) > > More details about what has been fixed/added can be found at PostGIS 2.0.0 > and PostGIS 1.5.3/1.5.4 fixes , What is new in PostGIS 2.0, PostGIS Raster > support, and improved Topology support These are 32-bit binaries so will not > work on 64-bit PostgreSQL. We are working on getting 64-bit binaries out as > well, but are running into some stumbling blocks. For the time being -- if > you are on Windows 64-bit, you must install the 32-bit PostgreSQL to use > these. >   >   >     For PostgreSQL 8.4 w/ Loader/Dumper including GUI >     For PostgreSQL 9.0 (compiled against PostgreSQL 9.0.3) w/ Loader/Dumper > including GUI >     For PostgreSQL 9.1 (compiled against PostgreSQL 9.1.0 - note 9.1 > structure changed in beta 1 so this will not work with alpha releases >   > > but should be fine for all beta, release candidate releases and final 9.1.0 > release) w/ Loader/Dumper including GUI. Also includes extensions files to > install using new PostgreSQL 9.1+ CREATE EXTENSION syntax. To install: >   >   >     First copy the share/extension files into your PostgreSQL > share/extension folder. >   >     Then in any database you'd like to spatially enable, run at least > postgis (which contains both geometry and raster support) and >   > > optionally postgis_topology commands. >   >   >     CREATE EXTENSION postgis; >     CREATE EXTENSION postgis_topology; >   I currently have PostGIS 1.5 installed on PostgreSQL 8.4. This is what my folder structure looks like: ![enter image description here](http://i.stack.imgur.com/K3ZaG.jpg) I don't fully understand these instructions: > First copy the share/extension files into your PostgreSQL share/extension > folder. > > Then in any database you'd like to spatially enable, run at least postgis > (which contains both geometry and raster support) and optionally > postgis_topology commands. I've downloaded `postgisgui_pgadmin-2.0.0svn` and `postgis- pg84-binaries-2.0.0svn`. In my original folder structure I don't have a folder `extension` within `share`. Do I just copy the contents of `postgis- pg84-binaries-2.0.0svn` into this folder? What does it mean to `run at least postgis`? I've also tried to follow the detailed instructions at the following blog. * * * EDIT: I have installed Postgresql 9.1 with PostGIS 1.5.3. I am trying to figure out where to copy the PostGIS 2.0 binaries. It seems like I should be copying them into the folders here: ![enter image description here](http://i.stack.imgur.com/lq9Qn.jpg) but the advice suggests placing them in the share/extension folder. I tried copying the entire folder into `share/extension/` and then editing the `makepostgis.bat` file to make it relevant to my file system (the port is 5434 as I have two instances of Postgresql running).               set PGPORT=5434      set PGHOST=localhost     set PGUSER=postgres     set PGPASSWORD=****     set THEDB=template_postgis20     set PGBIN=C:\Program Files (x86)\PostgreSQL\9.1\bin\     set PGLIB=C:\Program Files (x86)\PostgreSQL\9.1\lib\      When I run this file as an administrator, I get the following error:               C:\Windows\system32>set PGPORT=5434         C:\Windows\system32>set PGHOST=localhost         C:\Windows\system32>set PGUSER=postgres         C:\Windows\system32>set PGPASSWORD=****         C:\Windows\system32>set THEDB=template_postgis20         C:\Windows\system32>set PGBIN=C:\Program Files (x86)\PostgreSQL\9.1\bin\         C:\Windows\system32>set PGLIB=C:\Program Files (x86)\PostgreSQL\9.1\lib\         C:\Windows\system32>set POSTGISVER=2.0         C:\Windows\system32>xcopy bin\*.* "C:\Program Files (x86)\PostgreSQL\9.1\bin\"     File not found - *.*     0 File(s) copied          C:\Windows\system32>xcopy /I /S bin\postgisgui\* "C:\Program Files (x86)\PostgreSQL\9.1\bin\\postgisgui"     File not found - *     0 File(s) copied          C:\Windows\system32>xcopy lib\*.* "C:\Program Files (x86)\PostgreSQL\9.1\lib\"     File not found - *.*     0 File(s) copied          C:\Windows\system32>"C:\Program Files (x86)\PostgreSQL\9.1\bin\\psql"  -c "CREATE DATABASE template_postgis20"     ERROR:  database "template_postgis20" already exists          C:\Windows\system32>"C:\Program Files (x86)\PostgreSQL\9.1\bin\\psql"  -d "template_postgis20" -c "CREATE LANGUAGE plpgs     ql"     ERROR:  language "plpgsql" already exists          C:\Windows\system32>"C:\Program Files (x86)\PostgreSQL\9.1\bin\\psql"  -d "template_postgis20" -f "share\contrib\postgis     -2.0\postgis.sql"     share/contrib/postgis-2.0/postgis.sql: No such file or directory          ....the errors continue....


title: Field Calculation - Python - If/Then over multiple data columns

text: I have some VB code saved that looks at two columns of data, and populates a third with the low value. Pre-logic Script Code:               Dim output as double     if [MIN_LEFT] < [MIN_RIGHT] Then       output = [MIN_LEFT]     else       output = [MIN_RIGHT]     end if          low = output      I've never used Python before, but I see that I'm going to have to if I want to do this calculation! But I'm stumped by the comparison of the two columns, and how to define that. Can anyone help me? (And I apologize if this is overly basic!)


title: How to intersect widget using script python arcpy and geoprocessing service?

text: Is it possible to develop an intersection widget with geoprocessing service using script python arcpy? I try to do that with ModelBuilder but after publishing, in geoprocessing service i don't find my inputs among the parameters.


title: Cache on demand feature for cached map services

text: I have a cached map service that I am using in a Silverlight based application. For this cached service, I have cached a core area and checked the option for On Demand Cache. I find that the areas within On Demand Cache are very slow. I am trying to understand how this process works. If I am zoomed in at the scale of 4000, it should be generating cache for a very small area, not sure why it takes much longer. I am using ArcGIS online tiling scheme. Does it generate the cache for all scale levels or just the one level. I am looking for some feedback. Thanks Jay


title: Summing attributes intersected by under a radius

text: Apologies - new user. I am working with QGIS on Mac OS X and with shapefiles downloaded from the US Census Bureau that show population data by state and by county. How do I:   1. select certain points (e.g. cities)   2. draw a defined distance radius around those points (e.g. 500km - At the moment I can only get the units to be in degrees)   3. Sum the populations in the counties that are intersected by the circle? Any assistance is greatly appreciated!


title: How to debug 'Layer is not valid error' from Python plugin?

text: I'm trying to understand what goes wrong with a `QGIS` `python` plugin. Basically I'm trying to render a series of `POINT`. The query works fine, at least giving it by hand from inside PostGIS DB:               "SELECT  gid, CAST (gid as text)|| '_' || (ST_DumpPoints(the_geom)).path[2] as key,  (ST_DumpPoints(the_geom)).geom as points FROM ways"      I concat gid and path[2] because after splitting the_geom (`MULTILINE`) field into several `POINT`, gid are duplicated and I think they are no longer usable as key (with this trick key should be unique, but it seems not solve my problem). I tried to create the layer this way:               def createPointLayer(self):         uri = self.db.getURI()         query = "SELECT  gid, CAST (gid as text)|| '_' || (ST_DumpPoints(the_geom)).path[2] as key,  (ST_DumpPoints(the_geom)).geom as points FROM ways"              uri.setDataSource("", "(" + query + ")", "points","","key")         layerName = "NodeLayer"         aLayer = self.iface.addVectorLayer(uri.uri(), layerName, self.db.getProviderName())      I think the query works fine (QGIS takes a long time before display the error. I guess something goes wrong after the query has been executed). When the plugins run QGIS output the following error: > Layer is not valid: > > The layer dbname='pgrouting' user='user' password='user' key='gid' > table="(SELECT gid, CAST (gid as text)|| '_' || > (ST_DumpPoints(the_geom)).path[2] as key, (ST_DumpPoints(the_geom)).geom as > points FROM ways)" (points) sql= is not a valid layer and can not be added > to the map I already had to struggle against the **Layer is not valid error** , but I don't know how to debug it. So I have 2 questions:   1. Could anyone tell me what's my mistake in this specific case?   2. How can I debug 'Layer not valid error'? How can I get more information on which is the problem? EDIT: QGIS log console display the following message: > No key field for query given


title: GeoExplorer dependencies are missing after git clone

text: When I set up GeoExplorer in development mode today I had to manually clone the external js libs that are referenced in git (buildkit, stick, geoext, gxp, openlayers, printpreview, proj4js). They were not included after running "git clone git://github.com/opengeo/suite.git", which prevented me from building GeoExplorer with maven. After adding them manually everything works. What am I missing in my git/maven configuration to skip this step?


title: Postgresql (soln would be basic SQL I think): ORDER by distance to a point and return those distances in array pairs w/ the full records

text: so I know there have been some questions asked about how to calculate distance using postgis and how to order results by their distance from a respective point. I'm not very talented at SQL (at all) and so was wondering if I could get some help returning: a.) a list of records sorted by their distance to a point b.) with each record in a small array containing the record itself and the distance from the record to the point passed as the second geometric var I've already learned how to order my results from this question: Order by nearest - PostGIS, GeoRuby, spatial_adapter So I'm guessing maybe my solution might look something like:               matches =ReferencePoint.where(custom_sql).order("ST_Distance(reference_points.geom, ST_GeomFromText('POINT #{other_loc.x #{other_loc.y}', 3361))")      Which I assume works...haven't actually tested it.... ...that would be perhaps amended like:               #aborted code      OK, I don't want to run a second select command, I want to add the sql to the code in my `order` command. Which means I probably have to use #select....don't know how to do that simultaneously after an `ORDER` command. So yeah, I need some help. Thanks EDIT: Example data: - lets say I had 100 boats scattered about around a buoy. I want sql to return those boats in order of their distance to the buoy, closest first. I'm going to use pseudo-sql code to represent the last part of this (text inside of select method):               buoyshape = Buoy.find(2341).projected_geometry_3361 #I document srid codes in my columns     Boat.order{st_distance(:proj_shape_3361, buoyshape)}.select("* + st_distance(proj_shape_3361, ?)", buoy_shape)      OK so I don't know how to write sql so the last #select component doesn't do what its supposed to do. But basically I want the distance between the buoy and boat to be returned as an additional column for each record.


title: how to limit zoom level in OpenLayers

text: I would like to have only 3 levels of zoom available (from level 6 to level 8).               demolayer2 = new OpenLayers.Layer.WMS(     "abc...","http://localhost:8080/geoserver/gwc/service/wms",     {layers: 'def...', transparent:"true", format: 'image/png', numZoomLevels: 3, minZoomLevel: 6 },{isBaseLayer:false},          { tileSize: new OpenLayers.Size(256,256)});          map.addLayer(demolayer2);           map.addControl(new OpenLayers.Control.LayerSwitcher({'div':OpenLayers.Util.getElement('layerswitcher')}));          map.zoomToExtent(new OpenLayers.Bounds(-4.601615515076983,39.8769407866263,-3.0527184873294764,41.16710732525929));


title: TIFF to MBTILES conversion on a mac

text: I have a large number of tiff files (map tiles) along with tab and tfw file types as part of an Ordnance Survey download. This is for an offline iOS app. I managed to get the tiffs into QGIS with correct projections and save as a virtual raster. What I need to do now is somehow convert or save it to an mbtiles format. Really appreciate any good explanations on how this can be done for a mapping novice, on a Mac! Thanks B.


title: Geodatabase feature to KMZ - display issue with stacked lines

text: I am working with a line feature in ArcGIS. My end product is a KMZ for use in Google Earth. However, there are lines that are stacked on top of each other in the feature. When I convert my feature to a KMZ only one feature is clickable in Google Earth, meaning clicking on stacked lines only displays the attributes of the topmost line. Is there any sort of method to deal with stacked lines? I know that point features that are stacked expand when you click on them to display all stacked points, but this is not a feature with lines apparently.


title: Server Side Thumbnails

text: I am using a combination of base maps from Google, Bing and OpenSteetMap, along with GeoServer and some geo-spatial data stored on a SQL server, as well as OpenLayers to display maps on a website. I am persisting the current center, zoom level, layer configuration and filters in a database. What I am struggling with, is generating a server side thumbnail that can be used to help differentiate saved map views on a "maps" page. Is there a utility, preferably .net, that will help me get an image, or images, back from the basemap sources I named, as well as geoserver WMS? Or at least one that could help generate the URLs for a given center and zoom point? Querying GeoServer WMS seems pretty straight-forward, the basemap sources, not so much.


title: Rendering geoJson polygons in google maps android apiv2

text: I'm in the process of publishing spatial data stored in a postgis database to a google maps based android app and am trying to put together a plan for rendering my data in the most efficient manner possible. The portion of my dataset I'm concerned about consists of about 200,000 polygons that will be filtered by the user name, bounding box etc. I can see situations where I would like to display up to 1000 heavily simplified polygons when zoomed out, but I doubt that I will be displaying more than 20 on a screen without simplifying the data. My current strategy is to query my database based on the filtering factors mentioned above, return the query as a geoJson and render the resulting data directly into google maps. Is this an appropriate strategy or should I be looking at using another service to help render my data before it is loaded into google maps. If it is an appropriate strategy I could use some advice on how to render geoJson polygons in the google maps android api. I'm quite comfortable with the SQL end of this project but I'm not as experienced with javascript or the google maps android apiv2. Thank you for the help.


title: Select features within a country or bounding box effectively using PostGIS

text: I imported the planet OSM file into a PostgreSQL 9.2/PostGIS 1.5 server using `osm2pgsql`. I now want to select certain map features which lie within the bounding box formed by a countries borders, e.g. all rivers and lakes within Germany. I've come up with a working query, but it's horribly slow:               osm=# EXPLAIN ANALYZE     WITH     bbox      AS        (SELECT ST_Envelope(ST_Union(way)) AS way             FROM    planet_osm_polygon             WHERE     boundary = 'administrative' AND admin_level = '2' AND name = 'Deutschland')     SELECT    *                              FROM    (                                     (SELECT     way, tags, width, water, waterway                    FROM        planet_osm_line                 WHERE waterway='river')             UNION                       (SELECT     way, tags, width, water, waterway                 FROM        planet_osm_polygon                 WHERE        water IS NOT NULL)             ) AS q, bbox     WHERE    ST_Within(q.way, bbox.way);                                                                                QUERY PLAN                                                                                ----------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------      Nested Loop  (cost=1558037.52..1676764.41 rows=1 width=192) (actual time=305373.041..311176.612 rows=20435 loops=1)        Join Filter: ((planet_osm_line.way && bbox.way) AND _st_within(planet_osm_line.way, bbox.way))        Rows Removed by Join Filter: 851602        CTE bbox          ->  Aggregate  (cost=18.27..18.28 rows=1 width=1005) (actual time=374.222..374.222 rows=1 loops=1)                ->  Index Scan using planet_osm_polygon_boundary_admin_level_name_idx on planet_osm_polygon  (cost=0.00..18.26 rows=1 width=1005) (actual time=27.989..34.350 rows=8 loops=1)                      Index Cond: ((boundary = 'administrative'::text) AND (admin_level = '2'::text) AND (name = 'Deutschland'::text))        ->  CTE Scan on bbox  (cost=0.00..0.02 rows=1 width=32) (actual time=374.227..374.228 rows=1 loops=1)        ->  Unique  (cost=1558019.24..1564213.69 rows=412963 width=1131) (actual time=299869.930..309268.592 rows=872037 loops=1)              ->  Sort  (cost=1558019.24..1559051.65 rows=412963 width=1131) (actual time=299869.929..300244.843 rows=874853 loops=1)                    Sort Key: planet_osm_line.way, planet_osm_line.tags, planet_osm_line.width, planet_osm_line.water, planet_osm_line.waterway                    Sort Method: external sort  Disk: 693392kB                    ->  Append  (cost=0.00..1314829.27 rows=412963 width=1131) (actual time=64.278..204999.538 rows=874853 loops=1)                          ->  Index Scan using planet_osm_line_waterway_idx on planet_osm_line  (cost=0.00..440480.43 rows=132130 width=1245) (actual time=64.278..142495.404 rows=595834 loops=1)                                Index Cond: (waterway = 'river'::text)                          ->  Bitmap Heap Scan on planet_osm_polygon  (cost=5271.85..870219.21 rows=280833 width=1078) (actual time=186.014..62347.033 rows=279019 loops=1)                                Recheck Cond: (water IS NOT NULL)                                ->  Bitmap Index Scan on planet_osm_polygon_water_idx  (cost=0.00..5201.64 rows=280833 width=0) (actual time=123.545..123.545 rows=279019 loops=1)                                      Index Cond: (water IS NOT NULL)      Total runtime: 311350.392 ms      I've already touched some PostgreSQL settings:               shared_buffers = 2GB     temp_buffers = 16MB     work_mem = 32MB      I tried to get rid of the external disk sort with a high value for `work_mem` and set `work_mem = 768MB`, but the disk sort is still in the query plan. Is there a way to speed things up?


title: QGIS' georeferencer layer doesn't appear?

text: I use QGIS 1.8.0 for Linux, and I want to draw roads manually on the map. ![enter image description here](http://i.stack.imgur.com/DKzFn.jpg) So according some tutorials, I use `Raster>Geoferencer>Geoferencer` to import a bitmap image including roads, and set Geoferencer like this screenshot. However, after "Start geoferencing", the Geoferencer layer didn't appear on the map. Just like this: ![enter image description here](http://i.stack.imgur.com/CtVVL.jpg) I've tried move the reference layer to top and toggle "update drawing order" but it doesn't work. What should I do?


title: Frequency, summarize tool options of ArcGIS for Desktop?

text: I'm looking for some general frequency, summarize tool. What I basically want is some tool like the summarize tool in ArcGIS (right click in the attribute table> summarize) that runs on all attributes and objects that I select (it would be great if that tool could output to excel for example but that is totally optional). I have way to many objects and attributes to use do this manually... I know this should be possible to create with ArcPy (although I don't find summarize in Arcpy), but this looks like something that can be useful to a lot of people so I believe something like this must already exist? ArcGIS is no requirement, if it runs on shape files or anything I can export to, that's fine. If this does not exist, then I will try to develop it myself and maybe post it somewhere online cause it looks like something useful. EDIT: I've written a small tool in python to basically do summarize. It works on shape files and summarizes to Excel (but other formats are easy to use too) and should be OS independent. If anyone would like to see it I will clean the code up and publish it. EDIT2: The tool can be found here: https://bitbucket.org/tcoopman/shapesummarize/


title: Using Google API 3 with Openlayers Brings POPUP

text: > **Possible Duplicate:**   >  Google Maps layer copyright popup every time map updated / user input I recently upgraded from google map api v2 to v3. It was working quite fine previously but now it does work but each time I zoom,pan I get a popup. I had look at many thred and https://getsatisfaction.com/opengeo/topics/how_to_get_rid_of_goggle_layers_map_data_windows one especially that explains that openlayes 2.11 has no problem (2.10 has problem and 2.11 being the latest I assume this). I am using openlayers 2.11 with geoext. How can I solve this problem? ![POPUP](http://i.stack.imgur.com/JV2Dl.jpg)


title: How to generate Tiles with QuantumGIS or Mapnik

text: I would like to know how to generate Map Tiles using QGis. I have heard that it is possible to create tiles using QGis, only, I don't know how to do that. Googling did not work for me since I am not really expertised with programming (Python etc.) and GIS. However, I am experienced with accomplishing tasks that I am not experienced. So, if somebody can direct me to a tutorial/guide (detailed tutorial is preferred), I can manage to generate tiles. I can also use Mapnik but as far as I understand, Mapnik does not have a GUI so I prefer QGis.


title: How to create svg symbols that have modifiable fill color, stroke color and stroke width?

text: I want to create some svg symbols for Qgis using Inkscape, those symbols must have attributes as described in the question. For the last two days, I have been experimenting as instructed by Sourcepole and mimicked some symbols came bundled with qgis that have desired attributes without any success. Finally I experimented with the simplest form of symbol: I created an svg that contains only a circle in Inkscape and tried to modify it. The original file (circle.svg) have this line:               <path        sodipodi:type="arc"        style="fill:#ffffff;fill-opacity:1;stroke:#000000;stroke-width:1.01988637;stroke-miterlimit:4;stroke-opacity:1;stroke-dasharray:none"      I modified it into:               <path        sodipodi:type="arc"        style="fill:param(fill) #ffffff;fill-opacity:1;stroke:param(outline) #000000;stroke-width:param(stroke-width) 1.01988637;stroke-miterlimit:4;stroke-opacity:1;stroke-dasharray:none"      Saved as circle_modified.svg and selected it as point symbol, but Qgis couldn't change all three parameters at the same time. What did I do wrong or should do differently? My system: Qgis 1.8.0 on Ubuntu 12.04 64 bit, Inkscape 0.48, Gedit 3.4.1 text editor. Thanks.