title: False DescribeFeature Request URL in GeoExt QueryPanel

text: I got the GeoExt QueryPanel working so far, it recognizes the layers available on my Geoserver but it is not able to get the attribute columns for the dataset. I configured an osgeo proxy on my Apache Server, I use the override- ext-ajax.js and have the `OpenLayers.ProxyHost = "/cgi-bin/proxy.cgi?url=";` class configured. A DescribeFeatureType request is sent when selecting the layer I want to filter. The same happens when trying to choose the attribute column in the query panel, but there are none displayed. The sent request URLs look like this:               http://localhost/cgi-bin/proxy.cgi?url=     http%3A%2F%2Flocalhost%2Fcgibin%2Fproxy.cgi%2Fwfs%3F     VERSION%3D1.1.0%26REQUEST%3DDescribeFeatureType%26TYPENAME%3Dcases%26query%3D      The strange part here is, that the request is sent via the proxy.cgi twice. The request URL should look like this, going passed geoserver´s WFS:               http://localhost/cgi-bin/proxy.cgi?url=     http%3A%2F%2Flocalhost%3A8080%2Fgeoserver%2Fwfs%3F     VERSION%3D1.1.0%26REQUEST%3DDescribeFeatureType%26TYPENAME%3Dcases%26query%3D      Does anybody have a clue what´s going wrong here and why the request URL is wrong? The strange thing is that even if I remove the Open.Layers.ProxyHost class and the override-ext-ajax.js from my application, the request is sent that way. This leads to the conclusion, that the source of error can be found in the QueryPanel.js file, yet I couldn´t spot it.   Can it, on the other hand be related to my server? Other requests sent in the application (GET & POST) work fine. The (relevant) code looks like this, maybe you guys can spot the error...               OpenLayers.ProxyHost = "/cgi-bin/proxy.cgi?url=";          var owsUrl = "/geoserver/wfs";          var symbolizer = OpenLayers.Feature.Vector.style["default"];          symbolizer.graphicName = "square";          var store;          Ext.onReady(function() {                      store = new GeoExt.data.WFSCapabilitiesStore({                 url: Ext.urlAppend(owsUrl,                    "VERSION=1.0.0&REQUEST=GetCapabilities"),                 protocolOptions: {version: "1.1.0"},                 autoLoad: true,                 listeners: {                     load: initApp                 }                 });             });             function initApp() {             var vecLayer = new OpenLayers.Layer.Vector("Query", {                 eventListeners: {                     featuresadded: function() {                         var extent = vecLayer.getDataExtent()                         mapPanel.map.zoomToExtent(extent)                     }                 }             });                        var queryPanel = new gxp.QueryPanel({                 title: "Query",                 region: "center",                 width: 350,                 bodyStyle: "padding: 10px",                 autoScroll: true,                 layerStore: store,                 bbar: ["->", {                     text: "Query",                     handler: function() {                         queryPanel.query();                     }                 }],                 listeners: {                     query: function(panel, store) {                         var rule = new OpenLayers.Rule({                             filter: panel.getFilter()                         });                         rule.symbolizer[queryPanel.symbolizerType] = symbolizer;                         var style = new OpenLayers.Style("", {rules: [rule]});                         wms.mergeNewParams({                             sld_body: new OpenLayers.Format.SLD().write({                                 namedLayers: [{                                     name: panel.selectedLayer.get("name"),                                     userStyles: [style]                                 }]                             })                         });                              store.bind(vecLayer, {                             initDir: GeoExt.data.FeatureStore.STORE_TO_LAYER                         });                     }                 }             });      Any help is appreciated! Cheers, Philippe


title: Is machine learning used much in GIS?

text: I'm new to GIS. Is machine learning used much in this field, and where? I'm investigating a job opportunity with a mapping company working on traffic awareness. My background is in supervised machine learning (in text classification) and some use of Hadoop, and they assure me that my skills would be appropriate. I've been interested in GPS and GIS for some time, but I want to make sure that my machine learning skills (which I've invested in heavily -- including an MS degree) would be put to good use.


title: Given city and state, how to get the nearest zip code in US

text: We have a situation where user enters a city, state for which we need to find a corresponding zip code. In case there is no corresponding zip code, we are required to find a zip code that close to the city, state. Can you please throw some ideas as to how to go about doing this using Google maps/geocode API. **EDIT:** Is it possible for US that a city entered by user does not have a zip code?


title: QGIS Server Capabilities Response doesn't contain any Layers

text: This question is similar as this one QGIS Server not working (problem with fast CGI), but with more details: I installed the latest QGIS from trunk on Linux (Ubuntu 11.04) in the non- default location '/usr/local/qgis/'. I installed apache2 from the repository. Next,I followed the steps outlined in this post on Linfiniti's blog (http://linfiniti.com/2010/08/qgis-mapserver-a-wms-server-for-the-masses/), but probably I did not understand the part about telling apache about the library path for my CGI; where do I add the ScriptAlias code to? I created a qgis project with two vector (shapefiles). Under WMS server tab of the project properties I checked 'service capabilities' and gave the title 'trial'. I copied the qgs project file to the /usr/lib/cgi-bin/VECEA directory. Trying http://localhost/cgi- bin/VECEA/qgis_mapserv.fcgi?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetCapabilities in the browser gets me the 'capability' xml file. It contains the user information, but there is no information on the data layers (the qgis project has two shapefile vector layers) and the bounding box is not defined (all set at 0), is that right (I did define an 'advertised extent' in the qgis project under the WMS tab) In QGIS I can connect to the server. However, it just gives me one row with the name of the project (ID 0), whereas I would expect it to show the two vector layers that were part of the project. When adding the layer, it is empty. **Edit** : and now it stopped working altogether. When trying the getCapabilities, I am getting a 500 Internal Server error. Looking in the apache2 error log I am finding the following: [Tue Sep 27 20:53:03 2011] [error] [client 127.0.0.1] Premature end of script headers: qgis_mapserv.fcgi /usr/lib/cgi-bin/VECEA/qgis_mapserv.fcgi: error while loading shared libraries: libqgis_core.so.1.8.0: cannot open shared object file: No such file or directory **Edit** : Got the 'getCapabilities' part to work again. In http://lists.osgeo.org/pipermail/qgis-user/2011-March/011549.html, it was suggested to add LD_LIBRARY_PATH to /etc/apache2/mods-> enabled/fcgi.conf, adding the line "DefaultInitEnv LD_LIBRARY_PATH /usr/local/qgis_git/lib". Only, I had to edit '/etc/apache2/mods-available/fcgid.conf' Still stuck with no layers showing when connecting to the server in QGIS


title: How to create a PostGIS layer using QGIS API?

text: I want to create a new PostGIS layer and save changes made to it. Presently, I am connecting Python to PostGIS using psycopg2 but I want to know whether there is any QGIS API function available. Note: I am able to open an existing PostGIS layer, but I don't know how to create a new one.


title: Openlayers transition effect

text: The OpenLayers TransitionEffect='resize' property is not working while using geoWebcache service of geoServer                 var demolayer = new OpenLayers.Layer.WMS(             "Eg", "http://localhost:8080/geoserver/gwc/service/wms",             {                      layers: 'Eg',                      format: 'image/gif' ,                      tiled: true,                      tilesorigin: [29.06213671874999, 28.912181118327]                     },                     {attribution: '©copy rights',                      tileSize: new OpenLayers.Size(256,256)},                     { transitionEffect: 'resize'                     }                        );      Regards, M_mam


title: Error when setting projection of a tif file in QGIS. Any ideas to solve this?

text: I get this error when trying to define a projection for a tif file in QGIS: ERROR 1: Too many points (3249 out of 3249) failed to transform, unable to compute output bounds. How do i work this out? any ideas? Cheers, Vasco


title: Does Web ADF 10 work with ASP.net 4?

text: I already have an ASP.net application using VS 2010 and .NET framework version 4 and ArcGis Server 10, and i need to add just a page, have map and toolbar. When I try to add map control and run the page the map is not displayed. So i just want to add page to that application having map and toolbar controls.So can anyone can suggest a solution? Thank you in advance.


title: Polygonize problem

text: > **Possible Duplicate:**   >  Polygonize function doesn’t overlay the shape file over the raster file When using "Polygonize" in Quantum GIS, I'm getting a polygon shapefile output that is a mirror picture of the (1-band) raster image. Where have I done wrong?


title: color gradient (or color ramp) in Arcgis and Qgis

text: I am currently working on my thesis and I would like to ask you: does anybody know the mathematical formula that calculates the color gradient in Arcgis and Qgis? In other words, which mathematical equation is used to represent the color transition from one color to another in a buffer? And also: is there a different equation that calculates the color classification in other shapes, such as square, circle etc? A good example can be found in How to Calculate Color Gradient (www.ehow.com), but this is rather general. I would prefer something that concerns Qgis and Arcgis. EDIT: 20/01/2013 ![enter image description here](http://i.stack.imgur.com/JewKA.png) The arrow doesn't mean anything...


The overriden function is public override void DrawDynamicLayer(esriDynamicDrawPhase DynamicDrawPhase, IDisplay Display, IDynamicDisplay dynamicDisplay) What's the difference between dynamicDisplay and Display? They both have similar functions like "DrawPolygon" and "DrawText" and then how about if I cast my map to an ActiveView and access ScreenDisplay, how are all of these display types different? How should I know which one to use?


title: Polygon dissolve in DB

text: I have a polygon feature class with more than 100,000 polygons in DB2 table. I would like to dissolve them based on a single field to reduce the row numbers. Is there any commands / SQL in DB2 to dissolve these without using ArcGIS Desktop tool?


title: BaseCustomLayer draw is not hit

text: I create an object derived from BaseCustomLayer which draws a polyline.   I call Map.addLayer(myBaseCustomLayer) and the function               public override void Draw(esriDrawPhase drawPhase, IDisplay Display, ITrackCancel trackCancel)      is hit. BUT if I set DynamicMapEnabled = true the draw function is not hit ! This is driving me crazy. Thanks for any hints.


title: GIS files for US Special Districts?

text: I'm looking for GIS files (shp, geojson, kml...) for US special districts, but I can't find them anywhere. Special districts are single-purpose local governement distinct from counties and municipalities. For instance, a water management, or a mosquito eradication district. Some counties have files, but only for districts inside their boundaries. I'm looking for something for the US as a whole. Thanks.


title: Specifying output coordinate system in ArcPy

text: I'm updating somebody else's old Python scripts, and trying to make sure I understand their logic and desired outputs while streamlining the code. In this particular script, it specifies an output projection in the processing environment settings.               import arcpy          Output_CS = "GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',...(etc, etc, etc)"     arcpy.env.outputCoordinateSystem = Output_CS      I think this means that all the subsequent geoprocessing operations (e.g., `arcpy.Buffer_analysis`) output in WGS84 without needing an additional `arcpy.Project_management` step. If that is the case, is there any accuracy risk to running a buffer analysis directly to a geographic coordinate system? In other words, does it do the projecting after the buffer's been appropriately sized? (The original input data is in a projected coordinate system, one of the NAD83 state plane systems.)


title: OpenLayers Drag&Drop a pin to select position?

text: Is there an OpenLayers example that allows the user to set a position with drag&drop. I would like to have a pin and set the position after the user drags & drops the pin to a new position.


title: GeoPortal - A database Reference has not been configured for connection Tag

text: I've followed the instructions on setting up to a tee. I am using Oracle as a database. My database structure has been setup. I receive the following error message when I try to search for anything:               Could not get a csw search engine: Did not get key = local in repository      com.esri.gpt.framework.context.ConfigurationException: A databaseReference has not been      configured for connectionTag:       When I login, regardless of what credentials I use, it displays "Welcome!" in green, althoug the login hasn't actually happened. In my gpt.xml file, I have made the following changes:               <interactiveMap         jsapiUrl=http://serverapi.arcgisonline.com/jsapi/arcgis/?v=2.3m         apServiceUrl=”http://server.arcgisonline.com/ArcGIS/rest/services/World_Imagery/MapServer”         mapServiceType="dynamic"         geometryServiceUrl="http://tasks.arcgisonline.com/ArcGIS/rest/services/Geometry/GeometryServer"         locatorUrl="http://tasks.arcgisonline.com/ArcGIS/rest/services/Locators/TA_Streets_US/GeocodeServer"         locatorSingleFieldParameter="City"         />      I commented out LDAP authentication, and opted for simple. I changed the U/N and P/W for this example.               <simpleAdapter>             <account username="exampleUseName" password="examplePassword" encrypted="false"/>             <roles>                 <role key="gptRegisteredUser"/>                 <role key="gptPublisher"/>                 <role key="gptAdministrator"/>             </roles>     </simpleAdapter>      My question: What does this error message mean? What haven't I configured? EDIT: Checked the Tomcat logs and found this:               org.apache.tomcat.dbcp.dbcp.SQLNestedException: Cannot create JDBC driver of class 'oracle.jdbc.driver.OracleDriver' for connect URL '*connection string*'


title: postgis - managing single dataset, multiple coordinate systems

text: I am setting up a forest inventory database that aims to store information about forest plot data from various countries. The current project is in Mozambique, others are in Laos, etc. Each project provides their forest data in different coordinate systems. Other than that the plot data (fields) is identical (nearly), so it makes sense to store it all in one table. I have come across the problem today where I want to view multiple project data in Mozambique. The CRS from one project is supplied in WGS84, the other UTM 36S. The geometry table has set the CRS for plots to WGS84, so when I load in QGIS I can't see the plots in UTM. What is best practice?   * To choose a single CRS (WGS84) for all plot data?   * To store plot data from different countries in different tables?   * Or depending on which record I am viewing QGIS can dynamically use the relevant projection?


title: extract labeled feature in ArcMap

text: We use Arcmap to label our features(points), now since there are too many features, so some of them are labeled while some not, now I wonder if we can extract the features labeled?


title: SID and/or ECW support on QGIS for Linux

text: Running QGIS on Ubuntu 11.04 and QGIS for Android on Galaxy Tab 10.1 (which is an awesome project by the way) NEED to figure out a way to get SID and/or ECW files loaded. Saw some posts about compiling GDAL manually but haven't had any luck. Has anyone successfully accomplished this? Thanks in advance!!!!!


title: How to get Geoserver working via IIS7?

text: I have installed GeoServer 2.1.4 on a Win 2008 64bit server that uses IIS as its default web server but there is no obvious documentation for getting it published via IIS. Does anyone know of a good documented procedure for getting a vanilla install of GeoServer accessible through IIS7? I don't really want to have to install Tomcat/Apache as I know my production server Overlords will not be happy about having it on the production machine. If I find the solution I will certainly be writing it up! Thanks people...


title: how to install QGIS 2.2.0 on mac OS 10.6.8?

text: Does someone know how to install the latest QGIS on a mac os 10.6.8? The dmg package available on http://www.kyngchaos.com/software/qgis does not run on my OS version. I was only able to find a developement version of QGIS from Sept 2013 which is not fully stable. Thanks


title: How to serialize a HTML form to GeoJSON data?

text: I'm developing a Cordova/PhoneGap mobile application and using NodeJS for server and MongoDB as database. I need to serialize GPS data from mobile phone to send it to my webservice with latitude and longitude as a 2dsphere index. When I use ajax command _data.serialize_ it brings me something like this:               {         "_id": "536a5d2a2708ab840ee9c8c8",         "lat": "0",         "lon": "0",         "time": "0",         "speed": "0"       }      How can I turn it to a GeoJSON 2dindexed?


title: Direct DBMS sql queries and ArcSDE geodatabase

text: i want to make updates from Arcobjects on several tables stored on a ArcSDE 10.1 versionned geodatabase based on Postgresql 9.0.5 , none of these update require spatial data , and making these updates through ArcSDE is not easy, but i think to make a direct Postgresql connection to make these sql updates. **Question :** what are the consequences and the problems that could happen from this queries , knowing that i wanna change just attribute data , actually i'm using just the DEFAULT version of the geodatabase, making several versions could make a problem toward the use of direct sql queries ?


title: Error open KMZ file from Oregon 650t Garmin

text: This is similar to this Previous post however; it does not mention about dpi or resolution to save a good quality image with very limited kmz files to be use for Garmin? What is the dpi when do you save the kmz file from Google Earth Pro into Oregon 650t Garmin ? Is there a standard file number to save a kmz file with a specific number ? I had a Non- user set up his ArcMap 10.1 and export user's map to jpeg and gave them to me to geoference them in GEP . My question Do I need to know the Map scale from ArcMap exported to jpeg file ? or regarding of what the user map scale set at and export jpeg at a higher resolution to preserve a good quality image ? Garmin will not open kmz files greater than 1 MB ? Am I wrong ? or do I need to resize the pixel size ? My Garmin GPS does not have usb chip in, so I am thinking I might need that for the CustomMaps folder ?


title: How to make an emf into a character marker symbol north arrow?

text: I have an emf that I would like to make into a North Arrow. How do I get the emf to be a Character Marker Symbol in ArcGIS 10?


title: Basic User Questions: Tracing a Border in QGIS

text: **The project:** I'm plotting a dilapidated cemetery, I'm using GPS waypoint data to plot the headstones in QGIS. The map I'm making is to show the orientation of the graves, and who is in them. I installed the Google maps plugin, now I want to trace the cemetery's road and border. Simple as that. WHY IS IT SO F* ***** HARD TO MAKE A LINE IN THIS PROGRAM!? Online tutorials are NO HELP. I don't get why I can't just make a layer and trace the google image with a simple pen, why do I have to do a trig problem to make a line! I uploaded the current .gpx file of the headstone waypoints and they sync up nicely though they will need some post processing. I am not experienced in GIS, but this project has to be done. I am very frustrated with this program, can you help me? -J


title: What simple solutions are available for sharing big data with end users?

text: I have approximately 1700 raster tiles showing classified land cover with a total size of approximately 300GB. I am looking for a simple solution to privately serve these data to non-technical end users. The end user needs to visually assess these tiled raster products. Very little analysis would be required, only the ability to make the raster tiles transparent with a color map and overlay them on aerial imagery. My initial thought was to create a gigantic mosaic of these tiles and serve the file via FTP. Is there a better method to serve these data to the end-user, perhaps one that could elegantly handle tiled data? I have very little web mapping experience, which should weigh into the solution. An 8TB, dual processor server and the ArcGIS family of products are available to me.


title: snapshot of a map with Mapwingis

text: I'm trying to get a snapshot of a map with Mapwingis. My problem is that I tried dozen of ways but my picturebox stays desperately empty ! Thanks for any help, Olivier Mapwingis 4.8.8


title: Can OpenLayers2 render a completely custom font-face in a text or point symbolizer?

text: I am trying to style a feature using a custom font. Geoserver has an extension allowing external graphics in SLD to be defined by a font. OpenLayers doesn't appear to support this particular Geoserver SLD extension. Should OpenLayers.Format.SLD.v1_0_0_GeoServer do this? I did't see anything in that class to indicate it would. If not, is there some other way to get open layers to load a true type font from an external resource? I'm thinking something along the lines of CSS3 @font-face which would then make the font available when styles are created in the the many ways OL provides. As a last resort, I could export all the glyphs in the font to SVG, but I would rather avoid that if possible. **EDIT:** Several OpenLayers examples involve client side feature styling:   * Vector Behavior Example   * Styled Layer Descriptor (SLD) Example Geoserver provides an SLD extension that allows the use of TTF marks to style data. It's not clear to me if OpenLayers allows this or not. For example, could OpenLayers consume an SLD with a point symbolizer like:                <PointSymbolizer>             <Graphic>               <Mark>                 <WellKnownName>ttf://Webdings#0x0064</WellKnownName>                 <Fill>                   <CssParameter name="fill">                     <ogc:PropertyName>#0000ff</CssParameter>                 </Fill>                 <Stroke>                   <CssParameter name="stroke">#ff00ff</CssParameter>                 </Stroke>               </Mark>               <Size> 45 </Size>               <Rotation> 0 </Rotation>             </Graphic>           </PointSymbolizer>      which works fine in Geoserver. If OpenLayers could do this, it would seem that either the font would have to be on a search path somehow.


title: Error with data frame in python script

text: I'm running a python script from arccatalog and it works, however, I am trying to add layers to an mxd and need to declare the data frame. I keep getting this error message. type 'exceptions.AttributeError': 'unicode' object has no attribute '_arc_object' Failed to execute (StepByStep) Selection of my code is below.               # Import arcpy module     import arcpy     import sys              # Licenses     arcpy.CheckOutExtension("3D")         arcpy.CheckOutExtension("spatial")         arcpy.OverwriteOutput = True          # Script arguments     arcpy.env.extent = "E:\\School_Test\\SchoolTest.gdb\\Extent"         arcpy.env.cellSize = "3"          mxd = arcpy.GetParameterAsText(0)         arcpy.env.scratchWorkspace = arcpy.GetParameterAsText(1)         arcpy.env.workspace = arcpy.GetParameterAsText(2)         spRef = arcpy.GetParameterAsText(3)         Survey_Data = arcpy.GetParameterAsText(4)         Survey_X = arcpy.GetParameterAsText(5)         Survey_Y = arcpy.GetParameterAsText(6)         Survey_Z = arcpy.GetParameterAsText(7)         Sample_Data = arcpy.GetParameterAsText(8)     Sample_X = arcpy.GetParameterAsText(9)         Sample_Y = arcpy.GetParameterAsText(10)         Selection_Criteria = arcpy.GetParameterAsText(11)         Site_Grid = arcpy.GetParameterAsText(12)          df = arcpy.mapping.ListDataFrames(mxd, "Layers") [0]      Below this I declare variables then geoprocessing. Any suggestions?


title: ApplySymbologyFromLayer

text: I am a new user of Python on ArcGis10 i am trying to apply symbology settings to multiple files but to no success. I have tried my own ways copied syntaxes from various places still getting error. The following is the one that make sense to me but What am i doing wrong??               import arcpy     ... from arcpy import env     ... env.workspace = r"C:\Experiments"     ... inputLayer = ["b","c","d","e","f","g","h","i","j","k","l"]     ...      ... symbologyLayer = "a"     ... arcpy.ApplySymbologyFromLayer_management (inputLayer, symbologyLayer)     ...      ***Runtime error <type 'exceptions.RuntimeError'>: Object: Error in executing tool***


title: How to print to a specific scale with QGIS?

text: How can I print a map with qgis to a specific scale? E.g. I have a map of Austria (projection: albers equal area) and want a printout on paper at 1: 3000000? Thanks Kurt


title: Features in .shp can't be identified & won't render

text: Using QGIS 2.4.0 on Windows 7 Service pack 1. I have several .shp's, all polygons, that give me the following problem: 293 features are listed in the metadata and shown in the attribute table but the majority of them won't render, and can't be identified or selected in the map. I can select a feature in the attribute table and zoom to it but the feature appears empty and when one attempts to identify or select it on the map it says "no features are found." I've checked geometry validity. I've exported the .shp. This problem did not happen until after we updated a few months ago from Windows XP (yes I know) to Windows 7. There appears to be no pattern I can see that explains why some features are visible and others are not.


title: Point Pattern Analysis Algorithms

text: I'm looking for point pattern analysis algorithms and literature on point pattern analysis, online resources and book titles welcome. The topics may range from general descriptions of algorithms to concrete use cases in any field of research. _Update Jul 31, 15:54:_ I'd be especially interested in local methods to detect linear features in point patterns. e.g. create road/track geometries from gps points.


title: How can we generate GeoJSON from SQL data with latlong

text: Is there any way we can create GeoJSON from SQL server table with Latitude & Longitude Data?


title: GIS Real Time Information Dissertation Topics

text: Does anyone have any good websites for real time information GIS projects, or have any inputs towards a dissertation topic within this field? Where does one for instance find useble data for such a project? Any and all advice will be greatly appriciated.


title: How to customize layer context menu for a specific ArcMap project?

text: Is it possible to customize group layer, feature layer, raster layer and table context menus in the TOC for a specific ArcMap document (VB.NET)? The goal is to add custom tools, specific to the layer type, to that layer's default context menu. I understand how to do this for all documents, but not a specific document only. An alternative would be to create custom layers, but that seems like a ton of effort if a custom class cannot inherit directly from the group/feature/raster layer classes along with their implemented interfaces (want to take advantage of the built-in functionality, like drawing). When I tried to derive from GroupLayerClass, for example, the interpreter says that ComClassAttribute and ComSourceInterfacesAttribute cannot both be applied to the same class (likely showing my ignorance here). Does anyone have suggestions how to achieve this functionality?


title: two symbols for a single label in legend

text: I am working on a spatial plan where I'm supposed to have two states - the existing one and the future one. I have always worked with separate maps - a single map for a single state. But I was wondering is it possible to put those two states together, i.e. to have two symbols in legend for one label - one on the left side of label and one symbol on the right side. So far, I have one solution - to make two legends with appropriate symbols, put those legends close enough each other and turn off one of labels. But, is there a simpler way to do this? I hope you understand my question. Greets! Here's an image of this problem.![Sorry about cyrillic, I did this map a hurry, cyrillic is official language in my country. Title "Existing state" is a part of one legend, and title "Plan" is in another legend in which I turned off label and put it close enough so it looks like it is all one legend.](http://i.stack.imgur.com/9edxD.jpg)


title: What GIS editors work with GIS data that is stored in SQL-Server

text: If I store my vector GIS data in SQL-Server, what editors can I use to edit the GIS data?


title: How to clip a lines vector layer to a polygon layer in QGIS?

text: I've perused through previous questions and haven't been able to find an appropriate solution to my problem yet, so I wanted to post here. I have a lines layer (a contour map of mean annual precipitation in Africa) that I want to clip to a polygon layer (shape of the African continent). However every time I use the clip tool (Vector-->Geoprocessing tools-->Clip), only an empty shapefile without attributes is created. After looking at other questions in the forum, I've already addressed two potential issues: (1) I've put the layers and the project in the same Coordinate Reference System (EPSG: 4326 - WGS 84) and (2) I converted the lines layer from multipart to singleparts after seeing this solved the problem for another user. Something else to note, and I'm not sure if it even makes a difference: the lines vector layer was converted from a raster file using Raster-->Extraction-->Contour. I would really appreciate any troubleshooting advice, because I'm pretty stuck at this point. Thanks! Below are a couple screenshots to add context: Project overview: ![enter image description here](http://i.stack.imgur.com/8eHJm.jpg) Clip tool: ![enter image description here](http://i.stack.imgur.com/2W76M.jpg) * * * **Edit/Update:** I found a solution by clipping the original raster file using the Africa polygon as a mask. I then converted that clipped raster file to a contour map and it solved my problem. As Chris W pointed out below, there was a lot of data off the coasts and I may have not waited long enough for original file to process (i.e., 40+ hours). Thanks to everyone for their suggestions and troubleshooting tips.


title: Edit tilemill legend to add circle symbol

text: I am trying to edit my tilemill legend so that instead of displaying a red rectangle as below, it displays a red circle (I'm making a proportional symbol map). ![enter image description here](http://i.stack.imgur.com/CnaQN.png) What line would I need to add in my legend code (see below) to make this happen? ![enter image description here](http://i.stack.imgur.com/3oefp.png) ![enter image description here](http://i.stack.imgur.com/9XTjK.png) On the mapbox website it says : "The advantage with html/css is that you can quickly make edits to the legend directly in TileMill, and maintain the ability to manipulate the legend styling with css even after the map has been exported to MBTile format. However, you are limited to right angles and solid colors, and may have to write many lines of code to create a relatively simple design." https://www.mapbox.com/tilemill/docs/guides/advanced-legends/ Hopefully "limited to right angles" does not mean that I cannot create a circle?


title: POST to the CartoDB imports API from Google AppEngine / Python

text: Creating tables using the CartoDB SQL api is a known bug/missing feature, CREATE TABLE queries sent via the SQL API create shadow tables that never show up in the dashboard. One potential way around this is using the new imports api to post a small csv file in the desired schema to create a table that then can be updated/inserted to via the sql api. http://blog.cartodb.com/post/54101913823/got-files-weve-got-a-import-api This works well with files sent via curl. It would be nice to know a little more about what parameters the imports api expects and can respond to. Specifically, I am attempting to post using urlfetch in AppEngine, and not having much luck. Like this..               url = 'http://mycartodb.cartodb.com/api/v1/imports?api_key=myapikey'          with open('testtable.txt', 'r') as fh: #testtable is a simple 3 line csv         payload = fh.read()          result = urlfetch.fetch(         url=url,         payload=payload,         method=urlfetch.POST     ).content          logging.info(result)      These requests successfully hit the api...               {         "item_queue_id":"d18799a5-3849-4058-ac23-49ea7d6e194e",         "success":true     }      ...but when I query the job status it reports an error of unknown origin. One hitch may be that the "Content-Length" header cannot be set in AppEngine, and maybe a successful imports api POST request requires this?               {         "id": "blah-blah-blah",         "user_id": "blah-blah-blah",         "table_id": null,         "data_type": "file",         "table_name": null,         "state": "failure",         "error_code": 99999,         "queue_id": "blah-blah-blah",         "get_error_text": {             "title": "Unknown",             "what_about": "Sorry, something went wrong and we're not sure what. Try uploading your file again, or contact us and we'll try to help you quickly."         },         "tables_created_count": null,         "synchronization_id": null,         "success": false     }


title: Failed load delimited layer in QGIS python console

text: I'm trying to follow the map rendering section from the pyqgis cookbook. I am able render shape file using qgs vector layer, however getting error prompt while loading delimited text file through qgs vector layer. ![enter image description here](http://i.stack.imgur.com/OUxl0.jpg) my script is(running in QGIS python console)               from qgis.core import *     from qgis.gui import *     from qgis.utils import *     from PyQt4.QtCore import *     from PyQt4.QtGui import *     from PyQt4.QtXml import *          mapinstance = QgsMapLayerRegistry.instance()     mapinstance.removeAllMapLayers()          layer = QgsVectorLayer("f:\scripts\map\wb.shp", "wb", "ogr")          if layer.isValid():         print "shape layer successfully loaded - adding to map instance"         mapinstance.addMapLayer(layer)     else:         print "shape layer failed to load!"           uri = "f:/scripts/map/points.csv?delimiter=%s&xField=%s&yField=%s" % (",", "longitude", "latitude")     vlayer=QgsVectorLayer(uri,"points","delimitedtext")          if vlayer.isValid():         print "point layer successfully loaded - adding to map instance"         mapinstance.addMapLayer(vlayer)     else:         print "point layer failed to load!"               print mapinstance.mapLayers()     # create image     img = QImage(QSize(800,600), QImage.Format_ARGB32_Premultiplied)     # set image background color     color = QColor(255,255,255)     img.fill(color.rgb())     # create painter     p = QPainter()     p.begin(img)     p.setRenderHint(QPainter.Antialiasing)     render = QgsMapRenderer()     # set layer set     lst = [ layer.id() ] # add ID of every layer     render.setLayerSet(lst)     # set extent     rect = QgsRectangle(render.fullExtent())     #rect.scale(1.39800703)     rect.scale(1.1)     render.setExtent(rect)     # set output size     render.setOutputSize(img.size(), img.logicalDpiX())     # do the rendering     render.render(p)     p.end()          # save image     img.save('f:\scripts\out.png',"png")      and points.csv file contains following data longitude,latitude 88.3697222,22.5697222 87.3194986,22.3243243


title: How do I convert a .dxf to a .shp containing a field with height information in QGIS?

text: I received two .dxf for my project that represent my area of interest. One of them contains buildings, the other one terrain elevation. They both are 3D. I can convert them to shapefiles successfully using the Dxf2Shp plugin. They also display nicely compare to sample below). ![sample of my model area](http://i.stack.imgur.com/oA1M3.png) However, I loose the height information if I do so. The only field I get in the attribute table is "myid", that simply counts the polygons. Is there a way to convert the files to .shp creating a height field in the attribute tables, that takes the z values from the .dxf? This would already be fine, as the projection is metric. Best regards


title: Why would ArcGIS 9.3 ModelBuilder ignore precondition?

text: For one section of a model in ModelBuilder 9.3, the tool doesn't recognise that there is a precondition. I have saved, shut down and restarted, and refreshed the link between the input data and the tool as well as refreshed the precondition link with no luck. The tool is supposed to calculate fields in a feature class (an output from earlier in the model) but it doesn't pick up that the fields were added in the previous step (precondition step). Any advice would be appreciated! Here is a picture: http://www.tiikoni.com/tis/view/?id=a72e160


title: PostGIS - ST_Collect but only include POLYGONs not POINTs?

text: In PostGSI, is there anyway to turn a GEOMETRYCOLLECTION into a MULTIPOLYGON by dropping all non-POLYGONs/non-MULTIPOLYGONs? I'm doing a query that cuts a bunch of geometries (using `ST_Intersection`) and I want to group them together. I'm using `ST_Collect` to turn this into a GEOMETRYCOLLECTION. I could use `ST_Union`, but that's slower. Is there an alternative approach where I can drop all the non-POLYGONs/non-MULTIPOLYGONs from a GEOMETRY collection and just return a MULTIPOLYGON? I know would be less accurate than `ST_Union`, but it should be faster. For what I'm doing a small hit in accuracy is worth it if it speeds it up.


title: Where to find .mxd examples for ArcGIS 10

text: I'm beginner of ArcGIS . I'm trying to learn how to use ArcMap. I would like to find some example MXD files for ArcMap, but couldn't find anything googling. Does anyone have any idea where can i find any MXD samples for ArcGIS 10? Thank you for any help =)


title: Overwrite existing map service in ArcGIS Server 10.1 (through ArcPy)

text: I am using ArcGIS DT and Server 10.1. I want to over write existing map service though ArcPy and going through this link and getting an error`Runtime error Traceback (most recent call last): File "<string>", line 22, in <module> NameError: name 'draftPath' is not defined` and I am searching for this at my end. I just want to know is there any alternative to overwrite map service automatically after 1 hour in ArcGIS Server 10.1 ? Any help will be great !!


title: Adding a WMS - connection timed out?

text: I am trying to add a WMS service to qgis but am having some problems. After adding the url and name of the service and clicking on connect, I am receiving the following error message:- "Could not understand the response. The wms provider said: Download of capabilities failed: Connection timed out" Does anyone know what the issue might be?


title: Nearest neighbor analysis not showing results

text: I'm trying to do a nearest neighbour analysis on fish shoals distribution. I go to Vector > Analysis Tools > Nearest neighbour analysis. I choose the appropriate layer and when I hit ok the black and white wheel is on, and nothing happens; I think it IS doing the analysis because the blue bar is 100%; but just not showing anything. Has anyone had this problem and know how to solve it? I'm using: QGIS 1.8.0-Lisboa GDAL 1.9.1 With Mac OS X 10.7.5 Cheers, Nathalie


title: Intersection with JTS Geometry

text: I want to calculate an intersection between two geometries with the intersection function from JTS Geometry. What happens, if the geometries have different SRIDs?               SimpleFeature sf1=...;      //SRID 31467     SimpleFeature sf2=...;      //SRID 4327          Geometry geom1 = (Geometry) sf1.getDefaultGeometry();     Geometry geom2 = (Geometry) sf2.getDefaultGeometry();          Geometry overlay = geom1.intersection(geom2);      Does the function calculate the intersection in respect to the different SRIDs or does it ignore it? Which SRID does the calculated geometry have? I didn't find any documentation about that.


title: How to get CQL_FILTER parameters from vector polygon

text: map projection is EPSG:900913, displayprojection for map EPSG:4326. vector layer is polygonLayer = new OpenLayers.Layer.Vector("Polygon Layer"); with draw polygoon option with following code I get CQL_FILTER               var myLatLonSquare = polygonLayer.features[0].geometry.transform(new OpenLayers.Projection("EPSG:4326"), new OpenLayers.Projection("EPSG:4326"));          var testArray = "WITHIN(GEOM, ";          testArray = testArray + myLatLonSquare;          testArray = testArray.replace("MULTIPOLYGON(", "POLYGON");          antennaLayer.mergeNewParams({         "CQL_FILTER": testArray     });      as it's shown on the images(the site does not let me to upload images yet) it gets different location coordinates. How to solve it? UPDATED: ![enter image description here](http://i.stack.imgur.com/gdVUr.png) ![enter image description here](http://i.stack.imgur.com/lr74c.png) ![enter image description here](http://i.stack.imgur.com/HN1eJ.png)


title: Working out percentage in QGIS field calculator

text: I have a project in QGIS where one of the things I need to do is work from a police data file with 9 months of crime information, linked to details of county population (age, employment, etc) to report the % of different crimes in each county (7 counties). When I follow the example screenshot given, word for word, to work out a %, I get an message which says: "An error occured while evaluating the calculation string.". I'm not sure where i'm going wrong. Any guidence would be a lifesaver.


title: How to choose first letters of words in an attribute table?

text: I want to write a script (Python or VB) in Field Calculator which can pick up or choose the fist letters of words and put them in a new field. For example:               Lepidium catapycnon..........> Lc     Gold blossom tree............> Gbt      I'm thankful if someone can help me with it.


title: Error with Leaflet Routing Machine

text: I downloaded the Leaflet Routing Machine from https://github.com/perliedman/leaflet-routing-machine. Unzipped and opened the index.html in my browser to see how it looks, and got an error as shown in the image below:               ERR_FILE_NOT_FOUND file://router.project-osrm.org/viaroute?instructions=true&loc=57.74,11.94&loc=57.6792,11.949&jsonp=_l_routing_machine_0      ![enter image description here](http://i.stack.imgur.com/h5BcB.png)               index.html          <!DOCTYPE html>                   <meta charset="utf-8" />         <title>Leaflet OSRM Example</title>         <link rel="stylesheet" href="http://cdn.leafletjs.com/leaflet-0.7.2/leaflet.css" />         <link rel="stylesheet" href="../css/leaflet-routing-machine.css" />         <link rel="stylesheet" href="index.css" />


title: Making a map that can be easily resized or scaled

text: I've recently been dealing with a map requestor who wants maps resized at the very last minute -- from 8.5x11 to 22x34, or from large poster-sized to "something that can fit in a report" (those vague ones are the best). Obviously there are some "people management" things I can do to help him better plan his projects and requests, but beyond that I still need to produce multiple sizes of the same map with a minimum of time investment. It is tedious to manually adjust all the elements -- legend, title, labels, various text fields -- to create a new layout. I did find "Scale Map Elements proportionally to changes in Page Size" in the Page and Print Setup dialog: ![page and print setup dialog](http://i.stack.imgur.com/ee0Uj.png) However, it still leaves me with things to do:   1. It doesn't work on all map elements (particularly text rectangles)   2. It does nothing to the labels within the map; the 16-point font labels that are great for large poster maps are much too large for small figures (and similar scale issue if making the map bigger) What are additional strategies for making an ArcMap layout more easily scalable? (My goal is to have as few font-size tweaks as possible.)


title: Offline satellite map of a specific area for an iOS app

text: I have a client that wants an offline satellite map as part of an iOS application. The map is of a specific area and needs to be able to be loaded without a connection. From my research I'm pretty sure I can do what I need to do using TileMill and route-me. There are two things that I don't understand:   1. Where can I get a satellite image that I can use in TileMill? (Paid is fine)   2. How do I get it into TileMill?


title: Hosting geoprocessing tool on geoserver

text: Like ArcGIS Server, is there a way to host, geoprocessing tools onto Geoserver?


title: How to convert NetCDF file to Grib2 or Grib files

text: I'm in trouble to read `NetCDF` file data. Basically NetCDF file contains Ocean current data. How do you convert a NetCDF file to a `Grib2 or Grib file`? I already have procedure to read data from Grib2 file using `OpenGrads version 2.0.` What I need overall is to generate image for ocean current components available in netCDF and display on map. Can anyone direct/correct me how to accomplish this task.


title: Wrong pixel sums within buffers in Esri ArcGIS 10.0

text: I have a raster of Zambia, 100m² resolution (about 650mb). Buffers (single) at 1,2,3,4,5,10,15,20,25,30,50 km. I'm running "zonal statistics as table" in the model builder for all buffers (actually, the buffers are run in the model builder as well). How do I know the sums are wrong? Some of the larger buffers have smaller pixel sums than the smaller buffers. I've cut down the dataset so that there are about 150 points to buffer. I've run the problematic points/buffers individually and they returned identical results.   This might be related or suggest where the problem is coming from: When I tried to clip the buffers to the country outline (most of the problematic buffers are around the edges or where the points are in close proximity to each other). I get a `Invalid Topology [Topoengine error.]` message. How can I resolve this?


title: Problem with nested query in PhP and PostGIS

text: I have a simple problem regarding a nested query. I would like to select all buildings in a given circular buffer and create a kml file of those selected objects. The following query does the job very well and returns the results when I run it in PgAdmin.               SELECT total_pop, ST_AsKML(ST_Transform(geom,26918)) AS kmlgm FROM h_buildings WHERE ST_DWithin(ST_Transform(geom, 26918), (SELECT ST_Transform(ST_GeomFromText('POINT(1116858.6062789 7086290.6680572)', 900913), 26918)), 1000)      the "POINT(1116858.6062789 7086290.6680572)" is the center of the circle, and the value 1000 is the radius. As you can see the query has another SELECT inside it (nested). The following PHP code is able to generate kml file based on the results of query. I have already tested it with a simple SELECT query and it works fine. However, when I use the above-metioned query it does not work, and it does not give an error. Does this mean that PHP does not support nested query?! I'm a beginner, Can anybody please help me to find where the problem is and how it can be solved? PHP code:               <?php          $dbconn = pg_connect("dbname='mydb' user='postgres' password='****' host='localhost'");     if (!$dbconn) {         echo "Not connected : " . pg_error();         exit;     }     $srid = '4326';          $sql = 'SELECT name, ST_AsKML(ST_Transform(geom,26918)) AS kmlgm FROM Buildings WHERE ST_DWithin(ST_Transform(geom, 26918), (SELECT ST_Transform(ST_GeomFromText('POINT(1116858.6062789 7086290.6680572)', 900913), 26918)), 1000)';            # Try query or error     $query_result = pg_query($dbconn, $sql);     if (!$query_result) {         echo "An SQL error occured.\n";         echo pg_last_error($dbconn);         exit;     }          $kml ="<?xml version='1.0' encoding='UTF-8'?><kml xmlns='http://earth.google.com/kml/2.1'><Document>";          for ($i = 0; $i < pg_numrows($query_result); $i++) {              $townname = pg_result($query_result, $i, 0);         $townkml =  pg_result($query_result, $i, 1);              $kml .= "<Placemark><name>Town: ".$townname."</name><description>".$townname."</description>".$townkml."</Placemark>\n";     }          $kml .= "</Document></kml>";          header('Content-Type: application/vnd.google-earth.kml+xml');     header("Content-Disposition: attachment; filename=sample75.kml");          echo $kml;          pg_close($dbconn);     ?>


title: Finding largest number of observations within a given radius

text: I have the following problem: Given the location of insured properties (say houses with fire insurance), i have to find the largest sum of insured values which fit into a 300m radius. Now, my first approach would be to divide the map into grids of appropriate size(e.g. 55m), assign the properties to the grids and calculate the sum of the grids in the first step, and then aggregate blockwise (10x10 grids), and finally search for the largest block. But I have a fealing that there must be a more elegant approach.


title: ArcMap Draw Camera Range with ArcPy

text: i am working on some camera data. I have some points which consist of azimuth, angle, distance, and of course coordinate field attributes. In arcmap I want to draw shapes like this: ![Sketch of angular sector](http://i.stack.imgur.com/iah8i.jpg) The shape is pink. Can I get the shape with postgres-postgis? Thanks.


title: Can a folder connection in ArcMap Catalog window be added using Python?

text: I have created a set of general scripting tools for map automation, etc. In one of my tools a geodatabase is created in the workspace folder the script points to. Just out of curiosity, is there any way I can add a folder connection in the ArcMap catalog window once the GDB is created so the user doesn't have to manually create one (assuming they haven't already done so), or is that out of ArcPy's scope?


title: Define the "system" in Geographic Information System

text: When speaking or writing about GIS all too often the concept in peoples' minds for Geographic Information Systems can be rephrased as "the software for map data". While software is integral to GIS, thinking of spatial software as _the GIS_ is a narrow view, unnecessarily limiting our conversation and understanding. I have some ideas about what the "system" in geographic information systems really means, which I'll throw into the punch bowl below, but let us hear yours first. Secondly, is it too late to bring into public consciousness an expanded truer definition which will arise in answer to the GIS moniker, or do we need to find a new term to encompass it? if the latter, what do you suggest?


title: How Do I Disable Clustering in SharePoint?

text: I'd like to disable or remove Clustering points within Sharepoint. How do you do it? Thanks all


title: Why does exported map from print composer differ from map canvas? QGIS 1.8.0

text: I have a small project consisting of a Google maps satellite open layer with two new shapefiles overlaid on it. I have then used print composer to add a grid and some labels. I have also selected lock layers for map items and draw map canvas items in the map tab. However, when I come to import the exported map in to another document, the layers have shifted and the map has zoomed in to an offset position. This occurs in both pdf and image formats. How do I ensure the map I produce in print composer does not alter when saved for export? I hope someone can help as I have noticed another user asked this same question back in February without any response (User 15214 on Feb 15th). Thanks in advance


title: Find which grid cells contain a point

text: I have a dataframe with x,y coordinates. I've sampled 10 x,y coordinates from my dataframe:               y <-c(24.637032,21.569813,21.687,19.677183,24.72438,21.433136,24.637032,22.21065,21.443762,25.273882)          x <-c(-156.02597,-112.00287,-112.00287,-105.27803,-159.47324,-157.858,-156.026,-109.68626,-112.12681,-112.15393)      I've converted these points to a `SpatialPoints` objects:               library(sp)     myPoints <- data.frame(x,y)     coordinates(myPoints) <- c('x', 'y')     proj4string(myPoints) <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84")      I've used these points to make a `SpatialGrid`object:               bb <- bbox(myPoints)     myPointsGrid <- GridTopology(cellcentre.offset =  bb[, 1] + (cs/2), cellsize = c(1, 1), cells.dim = ceiling(diff(t(bb))/cs))     myPointsGrid <- SpatialGrid(myPointsGrid, proj4string = "+proj=longlat +ellps=WGS84 +datum=WGS84")      I need two things 1. find which grid cells contain a point 2. Find the number of points in each grid cell. I tried...               over(myPointsGrid, myPoints)      ...but this has thrown an error. How can I complete this task?


title: How to convert shapefile to vector emf?

text: I want to convert a shapefile to an EMF file. The EMF file must be a vector image, not a raster image. Please suggest to me some good packages or tools. Thanks in advance.


title: ogr2ogr convert all shapefiles in a directory

text: I have a directory with several shapefiles. How can convert all these shapefiles to MapInfo with ogr2ogr. I know how I can convert one file. And I can make a batch script writing a line for each file. But isn't there an easier way to convert all the files in a directory (and subdirectory). Hope someone can help.


title: Batch Convert MapInfo .TAB Files to .KML

text: > **Possible Duplicate:**   >  ogr2ogr convert all shapefiles in a directory As the title suggests, i have a few hundred .TAB files to convert to .KML. This is likely to be a fairly irregular thing so FME (while brilliant!) is not really an option, and FW Tools doesn't appear to offer batch conversion, without possibly scripting the solution. There there isnt a way in MapInfo/MapBasic to do this either. Does anyone have any ideas? Thanks in advance AM


title: How to use ogr2ogr for batch conversion?

text: I have a thousand of arc coverage would like to covert to shapefile. How can I use the ogr2ogr to convert it in batch rather than convert it one by one? It takes a lot of time! :( Thanks for any idea!


title: Making multiple separate selections from attributes

text: I am using ArcGIS 10 with ArcEditor license. I have a 10 year composite of month-to-month temperature values for 17 weather stations. That makes 120 months. I have to select each and every month separately and create a layer out of that selection. Is there an easy way to do that so that I can avoid having to make the repetitive 120 selections?


title: OpenLayers weird problem with GoogleMaps base layer and WMS overlay

text: I'm trying to overlay a WMS layer over GoogleMaps base one. The WMS layer is defined like this:                var catastro_layer = new OpenLayers.Layer.WMS(                "WMS Catastro",                http://ovc.catastro.meh.es/Cartografia/WMS/ServidorWMS.aspx?SERVICE=WMS&VERSION=1.1.1,                {                        layers: "catastro",                        format:"image/png",                        srs: "EPSG:3785",                        transparent: false                }        );      GoogleMaps layers are defined just like here http://openlayers.org/blog/2010/07/10/google-maps-v3-for-openlayers/. As result, my WMS layer appears as one base layer more. So, I can select it _instead of_ the other base layers, not _over_ them. This screenshots shows what I mean: http://dl.dropbox.com/u/6599273/capa_wms1.png The WMS requests are ok: http://ovc.catastro.meh.es/Cartografia/WMS/ServidorWMS.aspx?SERVICE=WMS&VERSION=1.1.1&LAYERS=catastro&FORMAT=image%2Fpng&SRS=EPSG%3A3785&TRANSPARENT=false&REQUEST=GetMap&STYLES=&EXCEPTIONS=application%2Fvnd.ogc.se_inimage&BBOX=-1721973.3737,3228700.0733875,-1702405.4944625,3248267.952625&WIDTH=256&HEIGHT=256 But this isn't what I need. I want a GoogleMaps base layer and a transparent WMS overlay. Now, if I change the WMS layer defintion replacing transparent: false by transparent: true:                var catastro_layer = new OpenLayers.Layer.WMS(                "WMS Catastro",                http://ovc.catastro.meh.es/Cartografia/WMS/ServidorWMS.aspx?SERVICE=WMS&VERSION=1.1.1,                {                        layers: "catastro",                        format:"image/png",                        srs: "EPSG:3785",                        transparent: true                }        );      My WMS layer appears as a overlay: http://dl.dropbox.com/u/6599273/capa_wms2.png, just what I want. But now, the requests are incorrect: http://ovc.catastro.meh.es/Cartografia/WMS/ServidorWMS.aspx?SERVICE=WMS&VERSION=1.1.1&LAYERS=catastro&FORMAT=image%2Fpng&SRS=EPSG%3A900913&TRANSPARENT=true&REQUEST=GetMap&STYLES=&EXCEPTIONS=application%2Fvnd.ogc.se_inimage&BBOX=-1721973.3737,3228700.0733875,-1702405.4944625,3248267.952625&WIDTH=256&HEIGHT=256 As you can see, the epsg:3785 has been replaced by epsg:900913. But the only change I've made is "transparent: false" by "transparent: true". Why is this happening? Many thanks in advance, and best regards


title: Is it possible to produce a decent looking shaded-relief map at 50K-scale-or-larger using SRTM3 DEM?

text: I'm trying to make **a decent-looking shaded-relief map** of a particular area in **50K-scale-or-larger** using open-source tools. The best (and free) option I have for DEM data for that area right now is the **global 90m resolution SRTM3 data**. So far I have tried using **GDAL** tools (`gdalbuildvrt`,`gdalwarp`,`gdal_translate`,`gdaldem`) and **GRASS** tools (`r.resamp.rst`,`r.shaded.relief`) to produce the shaded-relief. **But the output always retained these characteristics** :   * **Image is too pixelated at 50K-scale-or-larger**. I've tried re-sampling them to 30m and even 5m using RST/bilinear/cubic methods, but the images are somehow still pixelated   * **The landform is not clear or too generalized.** I have to increase the z-factor up-to 4.0 to obtain better landform. I don't know if this is normal.. I've read most of the work done by Tom Patterson in http://www.reliefshading.com and http://www.shadedrelief.com, and I'm trying to emulate his work - but I realized that most of the pointers he gave revolved around small-scale maps, and the use of propriety softwares to get the job done (I can only afford to use open-source stuff for now). So I'm at a lost right now.. Perhaps the enhancement for the SRTM3 DEM can be done in GIMP somehow? **So my ultimate question: is it possible to make use of SRTM3 data to make a decent 50K-scale-or-larger shaded-relief maps? If I can, how should I go about processing it?** Appreciate any leads on this matter, thanks!


title: GDAL ECW Support on OSX

text: I have a large number of ECW format raster images and would like to use with GDAL (and ultimately in QGIS) on Mac OSX. Having seen that there is no official ESRI support for the ECW library on any other platform other than Windows, I managed to get hold of the old version of the libecwj2-3.3 library. There are a few versions knocking around, but none of them will successfully compile. Anyone had any success compiling the ECW library on OSX (10.6.7) and then compiling GDAL with ECW support. Thanks. Jon.


title: Online/offline tool for collection of coordinates for simulation purposes

text: Is there any good free online tool to collect coordinate positions given, say a trajectory; or a polygon. For example, retrieve all coordinates between two points, every 500 meters. Many thanks.


title: FME assigning sequential start and end for each line segment

text: Hi hope someone can help me. To simplify, I have a line that is 500 meters long. As part of its attibutes it has a attribute called 'Start Chainage' of 0 and an attribute called 'end chainage' of 500. In FME I have used snipper to cut the line into lengths of 20 meters. What I need is to give each one of those segments their own sequential start and end chanages. So the first segment will have a start chainage of 0 and an end chainage of 20, the next segment will have a start chainage of 20 and end chainage of 40 and so on..... The reason for this is so that I can identify at what point (chainage) along the road some other feature crosses or intesects it. Thanks in Advance.


title: Error when displaying .shp from postgis with mapserver

text: I have postgis database with lines in .shp format in 3D. I made connection in map file:                    LAYER       NAME paths_3D      CONNECTIONTYPE POSTGIS      CONNECTION "host=www.*****.com port=5433 dbname=mountain user=**** password=****"      DATA "the_geom from paths_3D"      STATUS ON      TYPE LINE      CLASS         NAME "all_paths"         STYLE            COLOR 188 143 143         END     END     END      Then, I added checkbox in html file pointing on that .shp in postgis database:               <input type="checkbox" onClick="UpdateMap()" name="layer" value="paths_3D" [paths_3D_check] > Paths<br/>      So, when I open my web page and click on checkbox named Paths, it gives me this error:               writeLog(): Unable to access file. mslog.log msDrawMap(): Image handling error. Failed to draw layer named 'paths_3D'. msPostGISLayerWhichShapes(): Query error. Error (ERROR: column "oid" does not exist LINE 1: ...ary(ST_Force_2D("the_geom"),'NDR'),'hex') as geom,"oid" from... ^ ) executing query: select encode(ST_AsBinary(ST_Force_2D("the_geom"),'NDR'),'hex') as geom,"oid" from paths_3D where the_geom && GeomFromText('POLYGON((561889.240401714 5072468.31385521,561889.240401714 5093238.31385536,589594.131887715 5093238.31385536,589594.131887715 5072468.31385521,561889.240401714 5072468.31385521))',find_srid('','paths_3D','the_geom'))      I don't know what is wrong, please help.


title: solved adress tag: applying them in the right way

text: i wanted to run a osmconvert & filter with the following goal - to get the following tags out of the datas: addr:street addr:housenumber addr:city so i decided to do the following command               ./osmconvert hamburg-latest.osm.pbf --all-to-nodes -o=hamburgrestaurants_2.o5m     ./osmfilter hamburgrestaurants_2.o5m --keep="amenity=restaurant" -o=restaurant_2.o5m     ./osmconvert restaurant_2.o5m --csv="@id @lon @lat shop name  addr:street addr:housenumber addr:city website email" --csv-headline -o=restaurant_2.csv      but i failed - there were no streets no housenumbers and cities. Now i decidet to use overpass-API which seems to be less complicated. Can someone advice me here to go the first step. how to apply the correct adress tags for getting   * adress like    * street , housnomber and .... eg. mailadress and website.... **update:** hello and good evening - a little update: i found out what was wrong here: **Issue:** there's an extra space between "name" and "addr:street" in the above mentioned osmconvert csv string. when i remove it and all will be set. just want to share this solution with all others greetings


title: How does one export the list of stored WMS-servers in QGis?

text: Could anyone point me as to where QGis stores the list of WMS-servers you added? I have built a fairly large list of WMS-server on one computer, but the QGis- project file which gets accessed by my collegues doesn't load ALL of the layers (to keep startup quick). Occasionally though, a wms layer is required for which I have the needed url info saved, but they cannot acces it. Therefor, I'd like to make local duplicates of the WMS server list I have in my setup onto the other computers. **So far the idea, the problem is I'm clueless as to where to find the stored list, and if I can just copy/paste that into the matching folder on another computer** (Using QGis 2.0.1 still on most computers, fiddling around with 2.2 locally before I go and upgrade on all (~12) PC's.


title: Alternative satellite and hybrid layers for Europe published on WMS

text: For an application I am building I require to show satellite (high resolution) and street view layers, preferably also a hybrid view. Unfortunately, due to the application not running on a typical GUI, restrictions of the Google Terms prohibit me from using Google Maps. Looking for alternatives, I have had very little success in finding public WMS hosting alternative datasets. Hosting the data ourselves would require a redefinition of the project's scope which is, at this point, not an option. Looking at questions asking for lists of WMS, I have unfortunately found many servers that are no longer online or providing the services I need. Do you know any WMS hosting the layers I need or, in case you don't, what WMS list would you consider most up to data and regularly maintained?


title: QGIS custom form and autocompleted data entering with "validation check"

text: In Qt-Designer I created a custom form where I can choose serveral species names in two different combo boxes. This combo boxes are linked via "value relation" to the accordant non spatial tables (gattung and art) in my SpatiaLite database. The whole species name consists of "gattung" and "art". Because these tables ("gattung" and "art") contain a lot of species names (latin) wich one can choose I think about an eased editing feature.      1) I'd like to generate an autocompletion when entering the first 3 characters into a line edit (instead of a combo box) for "gattung" so that on can choose only these species that beginns for example with "Que" from a list.      2) When I've chosen a name that matches with these 3 characters for example "Que" for "Quercus", in the second field "art" only valide names can be chosen. For example "petrea" ore "robur" as a subset of the whole "art-names" which one can choose. Valide names are "Quercus robur" or "Quercus petrea". This could be done via line edit when entering 4 characters or with a combo box. I thought about python and a lettercode (Quepetr for Quercus petrea) to solve this but doesn't know where to start and also where to enter the code for the widgets in QT-Designer. I'm grateful for any help or idea!


title: Gis Exam Books and Exam Study Material

text: Do any one have books or study material for ARC GIS Desktop Associate , i can also buy used books .


title: Python GDAL: Write new raster using projection from old

text: If I read in a raster image as an array, then make some changes to the values in the array, how do I then save the array as a raster with the same projection information as the original array? In particular I am conducting processing on some ISIS3 cubes from Mars. These are not projected in any of the nice SetWellKnownGeogCS options. Perhaps this makes my problem somewhat unusual, but I thought it worth documenting my solution all the same.


title: calculating st_simplify tolerance for openlayers

text: I want to display a vector layer in openlayers. Users can scale the map before they load the vector layer and I'd like to simplify the geometry with st_simplify before sending it to the client. Does anyone have a simple formula I can use to calculate an appropriate st_simplify tolerance based on the map's properties (eg zoom, scale or extent?)?   * http://postgis.refractions.net/documentation/manual-1.4/ST_Simplify.html   * http://dev.openlayers.org/docs/files/OpenLayers/Map-js.html (I don't think this is relevant but have a php middle layer and am passing the layer as kml) I've googled around and maybe I'm looking the wrong way because I can't believe i'm the first one with this question! thanks in advance!


title: What format is lat/long stored in OSM PostGIS?

text: **Question:**   How do I to get the lat, long for all the toll booths stored in OSM for Virginia? **Current Query:**               > gis=# select * from planet_osm_nodes where ('toll_booth' = ANY(tags));                id     |    lat    |    lon     |                            tags                                 ------------+-----------+------------+------------------------------------------------------------       496668427 | 448545276 | -881565606 | {name,"Airport Parking",barrier,toll_booth}       525372041 | 463020204 | -874380697 | {name,"Swift Run Gap Entrance Station",barrier,toll_booth}       526241181 | 470778112 | -870435097 | {name,"Front Royal Entrance Station",barrier,toll_booth}       534508045 | 458513289 | -877777774 | {name,"Rockfish Gap Entrance Station",barrier,toll_booth}       563894077 | 470183057 | -857852817 | {barrier,toll_booth}       ...        Rest of the table **Additional Info**   1. OSM Import Log   2. Failed Attempts   3. Tutorial -- manually building a tile server   4. Data source -- Cloudmade/virginia.osm.bz2.


title: Programmatically modify layer datasource on project save event (QGIS)

text: I have a QGIS plugin which allows users to add WFS/WMS layers to the map. I would like to be able to scramble the user ID part (token) of the uri when the project is saved/closed and then use my own instance of QgsProjectBadLayerHandler to fix the uri when a project is opened. The result being that only users with the plugin and correct access levels are able to to view the WMS/WFS layers, so if a qgis project file is opened by someone without the plugin they will be unable to loads the layers. I am having no luck changing the uri when the project is saved. Here is what I have so far:               self.project.projectSaved.connect(self.project_saved)     self.project = QgsProject.instance()          ...          def project_saved(self):         reg = self.LoginAndStore.check_reg()         token = str(reg[1])         toc = self.iface.legendInterface()         toc_layers = toc.layers()         for layer in toc_layers:             try:                 dp = layer.dataProvider()                 src = str(dp.dataSourceUri())                 if token in src:                     no_token = src.replace(token, "NOTOKEN")                     dp.setDataSourceUri(no_token)                     #dp.reloadData()             except AttributeError:                 pass         f = QFileInfo(r"C:\Users\me\Documents\testing.qgs")         self.project.write(f)      Any suggestions on how to get this working or other methods to prevent unauthorized access are greatly appreciated.


title: How to map the frequency of overlapping points within an urban area?

text: I've been asked to do a favor for someone since I'm "the mapping guy" and could use some theoretical help by way of examples since my actual cartography skills are a little rusty. I have a series of points (~200) that represent events within a city over the course of a year. These events occur all over, but it's common that some points will overlap (e.g. events occurred at the same location). I'd like to display all of the points, but in the case of overlapping points would like to highlight the fact that the frequency of events is higher at those locations. I can think of a few ways to do this (size, color), but could use some examples or information on techniques that I could use to turn this data into an informative map. I'm currently using ArcGIS 10 to project and display the points.


title: How to configure osm2po to ignore certain class types?

text: I'm trying to import data using osm2po but am only interested in driving distances, so I'd like to omit all footway, cycleway, steps, etc. in the imported dataset. How can I configure osm2po to do this ? The documentation isn't completely clear... I tried leaving the config of wtr.finalMask as only car, which I interpreted it to mean it should only create ways which include the car type.               # only convert ways containing one of these flags     wtr.finalMask = car      But this is not the case in the resulting dataset. The ways with class types for just bike, foot, and ferry are all present. Do I need to just remove these lines from osm2po.config ?               wtr.tag.highway.pedestrian =     1, 62, 5,   bike|foot     wtr.tag.highway.track =          1, 71, 10,  bike|foot     wtr.tag.highway.path =           1, 72, 10,  bike|foot     wtr.tag.highway.cycleway =       1, 81, 15,  bike     wtr.tag.highway.footway =        1, 91, 5,   foot     wtr.tag.highway.steps =          1, 92, 5,   foot     wtr.tag.route.ferry =            2,  1, 10,  ferry     wtr.tag.route.shuttle_train =    2,  2, 50,  rail|car     wtr.tag.railway.rail =           3,  3, 50,  rail      Or is there some other mechanism for excluding them ?


title: using Z and M feature geometry to create line graph in Python?

text: I am a Python (v2.6.5) noob trying to plot a series of 3D polylines into graphs. So far I've come up with:                   # Import system modules         import arcpy, os, sys, string         from arcpy import env              print "Modules imported."              # Set environment settings         arcpy.env.overwriteOutput = True         arcpy.env.workspace = r"X:\HaulRoad\SpatialData"         arcpy.env.scratchWorkspace = r"X:\HaulRoad\SpatialData"              print "Workspaces set."              # Set local variables         infc = "\HaulRd.gdb\Treatment\XS_2009BEBRklnTerr_interpZ_Rte_isect"         graph_template = r"\XS\XS_Seg02.55.grf"  # template graph originally created with 3d Analyst profile tool, customized in Advanced Properties, then exported to grf format              print "Local variables set."                   #create list of fields of interest from FC or table of interest         fields = ['SegID','STA_calc','shape.Z','shape.M']                   #create empty list for cursor to populate         list = []              print "Table and empty list created."              # use search cursor to get field value for a given record          rows = arcpy.SearchCursor(infc, "", "", "SegID; STA_calc")         for row in rows:             list.append(row.SegID)             list.append(row.STA_calc)              del row, rows              print "Rows appended to temporary table.  List deleted."              #remove duplicates from list         list = dict.fromkeys(list)         list = list.keys()              print "Duplicates removed."              #create a temporary memory variable         memoryFC = "in_memory" + "\\" + "virtualFC"         print "Memory Feature created."              for n in list:             arcpy.TableSelect_analysis(infc, memoryFC, "STA_calc = " + str(n))             out_Seg = fields[0]             out_STA = fields[1]             out_graph_name = n             out_graph_pdf = r"\XS\Script\XS_chart_Seg" + str(out_Seg) + str(out_STA) + ".pdf"             graph_data = memoryFC                  # Create temporary feature layer              arcpy.MakeFeatureLayer_management(infc, "XS_lyr")                  # Create the graph from temporary feature layer             graph = arcpy.Graph()                  # Specify the title of the Graph             graph.graphPropsGeneral.title = "Segment " + out_Seg                  # Specify the subtitle of the Graph             graph.graphPropsGeneral.subtitle = "Station " + out_STA                  # Specify the title of the left axis             graph.graphAxis[0].title = "Elevation (ft)"                  # Specify the title of the bottom axis             graph.graphAxis[2].title = "Station (ft)"                  # Add a vertical bar series to the graph             graph.addSeriesLineVertical("XS_lyr", shape.Z, shape.M)                  # Output a graph, which is created in-memory             arcpy.MakeGraph_management(graph_template, graph, out_graph_name)                  # Save the graph as an PDF             arcpy.SaveGraph_management(out_graph_name, out_graph_pdf, "MAINTAIN_ASPECT_RATIO", 800, 600)              #clean in-memory         arcpy.Delete_management("in_memory")              print "PDFs created and memory cleaned."      I get an error in line 82 where I am adding the vertical line series. I assume the problem is that I am trying to plot values based on feature geometry as opposed to values in the attribute table. Is this possible without some intermediate conversion (e.g. converting vertices to points)? The source data is a polyline feature class with Z and M values in an Arc 10.0 SP4 file geodatabase. Thank-you,   Will PS. My overall intent is communciated in an earlier post: Arc 10.0 dynamic chart titles? I was unsuccessful in my attempts to perform from within Arc. -----EDIT 6/18/2014 error message I receive in IDLE follows:      Traceback (most recent call last): File "C:\temp\MakeGraphTest_v100_20140616.py", line 82, in graph.addSeriesLineVertical("XS_lyr", shape.Z, shape.M) NameError: name 'shape' is not defined


title: How to convert state name abbreviations to full names?

text: I have a data table with columns that includes city names (ex: San Jose) and state name abbreviations (ex: CA). Is there a way to define the state names as their full names (i.e., CA = California) and then combine city name and state name (ex: San Jose, California)? I have approximately 120 records. So what I think I have to do is define each state as their full names (ex: NY = New York, FL = Florida, and so on) first. I guess I could create a new field and add these to it, and then combine this new field with full state names with city names? Does that make sense or is there another way to do it? Need to achieve this using Python. Something like the following, but it doesn't seem right:               def stateNames(): # Add this pre-logic script in the code block       CA = 'California'       FL = 'Florida'       GA = 'Georgia'       HI = 'Hawaii'       MN = 'Minnesota'       NY = 'New York'       # and so on...     return stateNames()          stateNames + !City_Name! # Add this in the box below pre-logic script


title: Can't install ArcGIS License Manager 10.1 on Win Server 2003

text: i try to install a License Manager on a Windows 2003 Server. The setup starts, everything looks fine. When i press "next" the setup crashes, the setup window closes. No error. No entry in the event viewer... I already tried to create a installation log. > "C:\MyPackage\Setup.exe" /L*V "example.log" But nothing is created... Whats the problem? :/


title: How to reproject all featureclasses in a database to another database?

text: I am very new to arcpy and have ArcGIS 10.0 I would like to make a script where I would only have to change two variables in a script (database_input and database_output) so that I can reproject all vector featureclasses within the database and save them into another database(which exists already). Would anyone be able to put me on the right track? My code so far. It is basically an exported modelbuilder model using the batch reproject tool.               # Import arcpy module     import arcpy               # Local variables:     shp1 = "Path_to_Geodatabase_old\\old.gdb\\shp1"     shp2 = "Path_to_Geodatabase_new\\new.gdb\\shp2"     test_etrs_gdb__2_ = "G:\\Tifosy_Serverinstall\\test_etrs.gdb"          # Process: Batch-Projizierung     arcpy.BatchProject_management("Path_to_Geodatabase_old\\old.gdb\\shp1;Path_to_Geodatabase_new\\new.gdb\\shp2", test_etrs_gdb__2_, "PROJCS['ETRS_1989_UTM_Zone_N32',GEOGCS['GCS_ETRS_1989',DATUM['D_ETRS_1989',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',32500000.0],PARAMETER['False_Northing',0.0],PARAMETER['Central_Meridian',9.0],PARAMETER['Scale_Factor',0.9996],PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]", "", "DHDN_ETRS_1989_8_NTv2")


title: Partially updating raster (DEM)?

text: How to update an .img DEM in ArcView 9.3? I do have DEM which need to be updated only if few parts. The area to be updated has been produced in separate DEM with Spline Interpolation with Barriers. The area to be updated itself is S-shape, but I received larger- rectangular DEM (btw. I was hoping to receive -9999 value outside barriers, but it seems that outside area was also interpolated). My aim is to to update only area within the polyline/barriers. How to do that?


title: Adding feature object to PostGIS feature type from GeoServer

text: I am trying to set up a system that will post feature objects periodically to a PostGIS db and display them via GeoServer and OpenLayers on a web site. Is it possible to add the new feature objects to the db via GeoServer using the REST API or geoserver-manager (my system is Java based) or will I have to add the features directly into the PostGIS db?


title: Cut polygon and retain data

text: I have these two polygons:               library(sp); library(rgeos); library(maptools)          coords1 <- matrix(c(-1.841960, -1.823464, -1.838623, -1.841960, 55.663696,                         55.659178, 55.650841, 55.663696), ncol=2)     coords2 <- matrix(c(-1.822606, -1.816790, -1.832712, -1.822606, 55.657887,                         55.646806, 55.650679, 55.657887), ncol=2)     p1 <- Polygon(coords1)     p2 <- Polygon(coords2)     p1 <- Polygons(list(p1), ID = "p1")     p2 <- Polygons(list(p2), ID = "p2")     myPolys <- SpatialPolygons(list(p1, p2))     spdf1 = SpatialPolygonsDataFrame(myPolys, data.frame(variable1 = c(232,                                                                        242), variable2 = c(235, 464), row.names = c("p1", "p2")))     proj4string(spdf1) <- CRS("+proj=longlat +datum=WGS84 +ellps=WGS84     +towgs84=0,0,0")     plot(spdf1, col="red")                    coords1a <- matrix(c(-1.830219, -1.833753, -1.821154, -1.830219, 55.647353,                          55.656629, 55.652122, 55.647353), ncol=2)     p1a <- Polygon(coords1a)     p1a <- Polygons(list(p1a), ID = "p1a")     myPolys1 <- SpatialPolygons(list(p1a))     spdf2 = SpatialPolygonsDataFrame(myPolys1, data.frame(variable1 = c(2),                                                           variable2 = c(3), row.names = c("p1a")))     proj4string(spdf2) <- CRS("+proj=longlat +datum=WGS84 +ellps=WGS84     +towgs84=0,0,0")     plot(spdf2, col="yellow", add=T)      ![enter image description here](http://i.stack.imgur.com/i0lG6.jpg) I want to cut out parts of `spdf1` that are intersected by `spdf2`. However, I want `spdf1` to remain as a SpatialPolygonsDataFrame and to retain any information contained within `spdf1@data`. I¹ve tried gDifference as follows, which cuts out parts of `spdf1` that are intersected by `spdf2`, but then converts `spdf1` to SpatialPolygons (i.e. discarding information contained in `spdf1@data`).               gDifference(spdf1, spdf2, byid=T)      How can I cut into `spdf1` with `spdf2` but retain data contained in `spdf1@data`? I have checked and tried this similar question without how to overlay a polygon over SpatialPointsDataFrame and preserving the SPDF data?


title: How to host my GeoServer+Openlayers web mapping application on the internet?

text: I have created a web mapping application prototype using GeoServer (running on Tomcat as the servlet), Openlayers and PostGIS. Since it was a prototype, at that time, i did not bother about the hosting aspects of it. It was working fine on my (laptop) localhost: 8080. Now as part of my job, i am requested to develop a similar application but this time it has to go on the web, such that person x can access it via www. The agency i work for has a website running however they want it to be a separate thing of its own, hosted separately from the website. These are my questions:   1. what are the steps/options to follow in order to host a web mapping application on the web such that someone outside of the agency's intranet or for e.g. in my case, someone other than the localhost can access it as well? (Considering that i have one that is functioning and accessible via localhost only, and now i want someone else to access it also on their machine).   2. Is it possible to have the application accessible via the agency's website? If so, how can i implement this? In the event that the agency decides to have it as part of their website in future. I have read through some post and responses here that is somehow related however I'm getting confused so I am hoping someone can help respond to my specific query.


title: ArcGIS Online Map using URL does not work with PhoneGap

text: I wanted to show my public accessible ArcGIS Online map with URL on a simple Hello World PhoneGap project with Eclipse. Eclipse threw out the following error in LogCat. Web Console:RequestError: Unable to load file: //www.arcgis.com/sharing/rest/content/items/{MapId}?f=json status: 0 at http://serverapi.arcgisonline.com/js...3.4compact/:15 Just wondering anyone have had any successful experience with using ArcGIS Online web map with PhoneGap?


title: Google maps wms layer with 3857

text: I am adding wms layers to my google maps application. As Google Maps projection is a web mercator, I thought it would be better to request wms layers with srs=3857 (or equivalent) instead of the 4326. I searched a bit for fellow implementation and I saw that most (if not all) use 4326 instead. I wanted to experiment by myself, I request the 3857 to the map server and do some reprojections for the BBOX coordinates. I have the two running so I can compare them. The result is: - 4326, the tiles look OK acvtually, I can't notice much of a shift. - 3857 the wms tiles are completely off, huge lat shift My question is: I am right to think that I should be querying for a 3857 projected wms layer or should I stick to 4326? and is it normal that I have a big shift in my 3857 tiles? **edit** This is the code I use for creating the layer                   layer = new google.maps.ImageMapType({             getTileUrl: function (coord, zoom) {                 var proj = _map.getProjection();                 var zfactor = Math.pow(2, zoom);                 // get tile coordinates                 var top = proj.fromPointToLatLng(new google.maps.Point(coord.x * 256 / zfactor, (coord.y + 1) * 256 / zfactor));                 var bot = proj.fromPointToLatLng(new google.maps.Point((coord.x + 1) * 256 / zfactor, coord.y * 256 / zfactor));                 // convert from wgs84 lat,lng to new srs coordinates                 var topConverted = srsConversion(top);                 var botConverted = srsConversion(bot);                 //create the Bounding box string                 // handles 1.3.0 wms by ordering lat, lng instead of lng, lat                 var bbox;                 if(version == "1.3.0") {                     bbox = botConverted.lng + "," + botConverted.lat + "," + topConverted.lng + "," + topConverted.lat;                 } else {                     bbox = topConverted.lng + "," + topConverted.lat + "," + botConverted.lng + "," + botConverted.lat;                 }                 //base WMS URL                 var url = lURL + "&BBOX=" + bbox; // set bounding box                 return url; // return URL for the tile             },             isPng: true,             tileSize: new google.maps.Size(256, 256)         });


I have published an Image service with ArcGIS Server 10, containing the results of an interpolation operation. The Image which is returned from the service is in greyscale. I wish to render it with a stretched color ramp, with the following properties: ![enter image description here](http://i.stack.imgur.com/4nCsp.jpg) The RestEnd Point API, mentions that a rasterFunctions Parameter has to be given as a renderingRule, for the Export Image Operation. The given samples for the rasterFunctions are very limited, And I haven't been able to create a rasterFunctions which causes any noticeable change in the exported output. Where can I find more details about rasterFunctions? Or build the rasterFunctions Object to get the same result as shown in the window?


title: How to listen for the start of a Protocol event?

text: I have a service that sometimes takes awhile to provide a response. I'm trying to capture the initial call to the service so that I can put a "Loading...please wait" message until I get the response. I can't seem to find a way to get at the ajax portion of the HTTP Protocol request. Essentially, I'm looking for something analogous to jquery's ajaxstart. Has anyone done this or know how to do this? I looked into Openlayers event listeners, but couldn't find much there that helped.               var myVectorLayer = new OpenLayers.Layer.Vector("My Vector Layer", {           projection: "EPSG:4326",           strategies: [new OpenLayers.Strategy.BBOX({resFactor: 1})],           protocol: new OpenLayers.Protocol.HTTP({               url: "./resource/getData.php",               params:{                                   zoom: me.map.getZoom(),                   density: 1               },               format: new OpenLayers.Format.JSON(),            // Over-riding handle Response to get at feature attributes               handleResponse : function(resp, options)               {                 var request = resp.priv;                                 if (options.callback)                 {                   if (request.status >= 200 && request.status < 300)                   {                     // success                     if (resp.requestType != "delete") {                       resp.features = this.parseFeatures(request);                     }                     console.log('success');                     resp.code = OpenLayers.Protocol.Response.SUCCESS;                   } else                   {                     // failure                     resp.code = OpenLayers.Protocol.Response.FAILURE;                     console.log('failed');                   }                   options.callback.call(options.scope, resp);                 }                      me.updateData(resp.features);                    }             }),


title: GDS Workspace Path resovles as owner not user

text: I am having trouble retrieving the correct path to the current map layer dataSource path. I am using ArcGIS Desktop 10.2 (10.2.0.3348) The workspace path that gets returned is correct apart from the current users name is not used but replaced by the owner of the Geodatabase Server. Which does not exist on the users machine and so returns empty dictionary on commands like **arcpy.da.ListSubtypes(lyr.dataSource)**. e.g.               mxd = arcpy.mapping.MapDocument('CURRENT')     df = arcpy.mapping.ListDataFrames(mxd)[0]     for lyr in arcpy.mapping.ListLayers(mxd, '*.Couplings'):         print lyr.dataSource      Which returns;               C:\Users\<owner>\AppData\Roaming\ESRI\Desktop10.2\ArcCatalog\GWMAP_SQLEXPRESS.gds      What I expected would be;               C:\Users\<user>\AppData\Roaming\ESRI\Desktop10.2\ArcCatalog\GWMAP_SQLEXPRESS.gds      The user **owner** is the owner of the Geodatabase Server as listed in the **Properties->Administration** tab. I have tried setting the default environment settings in ArcMap & ArcCatalog. Plus setting the **arcpy.env.workspace** to the correct path. What am I doing wrong? Obviously I could easily use string replace to correct this but would prefer if there is a better approach. Any help most appreciated


title: Faster routing algorithm on pgrouting

text: I have downloaded data for india as a dbf file then converted it using osm2po to make it routable but the problem is when i make any routing query it takes a lot of time around 300 seconds and these results are after indexing source, target and id columns. So anyone can please tell how to enhance speed of pgrouting algorithm, or which algorithm should be used, time under 5 seconds would be OK for me, thanks in advance


title: How to set zoom level beyond 15 in Google Maps?

text: How do we set zoom level in Google Maps? We can set zoom level to 15 but how can we set zoom level to more than 15?


title: gdal_contour creates an empty shapefile

text: I am running: gdal_contour -a elev hazard_interp4.asc haz_ak_contour_4_10 -i 10 the hazard_interp4.asc file has values in band 1 (singel band gray). With the above command, I get a shapefile that has no values. If I burn a color schema into the raster with gdaldem color-relief, then I can successfully create a contour shapefile. However I would prefer to not have to do this step. **Update** : the output of hazard_interp4.asc is:               Driver: AAIGrid/Arc/Info ASCII Grid     Files: hazard_interp4.asc            hazard_interp4.asc.aux.xml     Size is 300, 300     Coordinate System is `'     Origin = (-19003400.000000000000000,11484799.899999998509884)     Pixel Size = (10698.000000000000000,-15086.032999999999447)     Corner Coordinates:     Upper Left  (-19003400.000,11484799.900)      Lower Left  (-19003400.000, 6958990.000)      Upper Right (-15794000.000,11484799.900)      Lower Right (-15794000.000, 6958990.000)      Center      (-17398700.000, 9221894.950)      Band 1 Block=300x1 Type=Float32, ColorInterp=Undefined       Min=0.005 Max=0.595        Minimum=0.005, Maximum=0.595, Mean=0.141, StdDev=0.141       NoData Value=-9999       Metadata:         STATISTICS_MAXIMUM=0.59498530626297         STATISTICS_MEAN=0.14081028623394         STATISTICS_MINIMUM=0.0050000008195639         STATISTICS_STDDEV=0.141377304974


title: Is there a way to test if a feature class is registered with geodatabase?

text: I've got feature classes inside an Oracle 11g ArcSDE geodatabase. I'd like to perform the arcpy method RegisterWithGeodatabase_management on all those feature classes. But before that, I'd like to test if the feature class is already registered, as the registration method will return an error in that case. How can I test if a feature class is registered with geodatabase or not?


title: Validation of Regression Kriging

text: I am interested in exploring the idea of generating a predictive model for a target variable in a spatial context. I would then like to correct the residuals of this model using kriging, i.e. in a nutshell you sum up the results from both models and take that as the final prediction: ![enter image description here](http://i.stack.imgur.com/bMfEF.jpg) The classical approach is to fit a GLM and then krige (is that how you turn it into a verb?) its residuals. I would like to fit something like a Random Forest. Which leads me to my question. How do you validate these models? If I wanted to do cross-validation do I separate the data (e.g. 10-folds) and reapeatedly fit the model, krige the corresponding residuals -> predict on test set, etc. Or do I have to separate the data into 3 sets, one to fit the model, see how it predicts on another data set, krige THOSE errors and then test on the final data set, repeat. What is the usual apporach to do this correctly? Thank you in advanced. I'm using R by the way.


title: Extract information from two layers

text: I have two shape files: road layer and a district layer. Is it possible to know which roads belongs to which districts? I'm currently using QGIS.


title: Unable to open project

text: I recently got a virus in my computer and, after having that fixed, was unable to open my Quantum GIS project. I received an error message:               Microsoft Visual C++ Runtime Library     Runtime Error!     Program: C:\PROGRA~2\Quantum GIS Lisboa\apps\qgis\bin\qgis.exe     This application has requested the Runtime to terminate it in an unusual way.     Please contact the application's support team for more information.      I then clicked "okay" and it kicked me out of qgis. I uninstalled qgis and reloaded it (version 1.8.0). When I tried to open my project, I got the same message. Then I opened qgis desktop and tried to open the project. I got a message that said `Exception - Bad allocation`. I need help. I have to recover my project.


title: Which free/cheap tools can open ESRI Geodatabase Files?

text: Ok I don't need anything advanced, I just don't know _anything_ about this GIS business. I'm trying to get local parcel information and the government site states it is in "ESRI geodatabase format." I've tried Google and get a wiki article that seems to be an advertisement and an esri.com that won't load a single page... so I'm stuck for information. I don't know what software I would need, and I'm really not hoping to plunk down some ridiculous amount of cash for a proprietary format our lovely gov't reps decided to use... Really, I just want to be able to correlate parcel listings with a location on a human (me) readable map. An added [super helpful] bonus would be the capability to feed in a large listing of parcels and see their locations at the same time, as an overview. Any ideas/help?


title: How to use the ESRI JS API to map SQL Server Spatial Types without ArcGIS Server

text: I have heard of people using the ESRI Javascript API to create webmaps with spatial data from SQL Server or other RDBMS's. From what I can tell, they are not using ArcGIS Server as a middleware. In particular, I recently heard of a company that uses a SQL Server database with spatial types and Entity Framework to link the database to the app, but then they are using the ESRI JS API on the front end to render the data. This is something that really interests me as I have worked with SQL Server, Entity Framework, ASP.Net, and the JS API separately and I do not have access to ArcGIS Server. Does anyone know of a particular architecture or workflow for simply using the JS API as a front end for a non ESRI backend?


title: Offset shapefile to start from prime meridian

text: I have Natural Earth shapefile (for example this file: link) covering whole globe in -180:180 range. I want to offset this file, so that it starts from Greenwich prime meridian instead -180. How can I do this with gdal/ogr?