WFS Filtered : The filters that worked for me in 1.8 - gives no results in 2.0

Using Quantum GIS 1.8.0 and 2.0.1 desktop on windows XP 32bit Server: http://vmgeoserver.vd.dk/geoserver/VD/ows?version=1.1.0 Setting filter for VD:vma_admdata:               "BESTYRER"  = '230' AND  "KODE_ADMIN_VINTERVEJSERVI"  = '230_V1' AND  "KODE_ADMIN_VEDLIGEHOLD_NI" = '230_N'      In 1.8 this returns the results I'm looking for. Metadata for the resulting layer states _Source for this layer:_               http://vmgeoserver.vd.dk/geoserver/VD/ows?version=1.1.0&SERVICE=WFS&VERSION=1.0.0&REQUEST=GetFeature&TYPENAME=VD:vma_admdata&SRSNAME=EPSG:25832&FILTER= BESTYRER 230 KODE_ADMIN_VINTERVEJSERVI 230_V1 KODE_ADMIN_VEDLIGEHOLD_NI 230_N       Doing the same in Quantum GIS 2.0.1 gives no error messages but the resulting layer is empty. Checking metadata the _Source for this layer_ is:               http://vmgeoserver.vd.dk/geoserver/VD/ows?version=1.1.0&SERVICE=WFS&VERSION=1.0.0&REQUEST=GetFeature&TYPENAME=VD:vma_admdata&SRSNAME=EPSG:25832&FILTER= BESTYRER 230 KODE_ADMIN_VINTERVEJSERVI 230_V1 KODE_ADMIN_VEDLIGEHOLD_NI 230_N       ie. identical. What am I doing wrong?


Distributed WMS (and TMS, optionally) server architecture, GPU rendering

What if we could use the power of dozens of servers to process a lot of incoming WMS requests? And more, using the power of GPU's to do the job of styling the maps for scalability and speed! The theme of my senior research project is exactly this, distributing the WMS service (and probably, TMS) to reach a maximum horizontal processing of GIS data. I'm very excited about this, and I would like to hear from you, the professionals in this field, ideas and criticisms regarding this approach.


Help with null values in rasters

I need to add different boolean rasters of the distribution of some species, since we need to obtain a single file that shows the locations where there are more than one species occurring. To do that, I rasterized all of the shape files showing the distribution and then reclassified them to obtain the boolean rasters. The problem is, now that I'm trying to add them all together using rastercalc, there are files with different extends even though I clipped them with the same shape file (the final single file -the one that I obtain by adding all of the booleans, is all cut and incomplete with null values)... how can I change them to the same extent in order to obtain the complete map... ? the first picture is the shape file I used to clip the other shape files (this is the polygon of the state I live in). The second picture shows an example of one boolean images with the raster of the shape I clipped them (the state's shape) The third one shows the last raster with all of the boolean rasters together... if you compare images 2 and 3, you will see there is lots of missing data ![enter image description here](http://i.stack.imgur.com/aNZrY.png) ![enter image description here](http://i.stack.imgur.com/eYUsU.png) ![enter image description here](http://i.stack.imgur.com/gxHVt.png)


Postgis 2.0 - How to get a Point from given LinestringM and a measure?

If I have a cable route linear referenced and stored in a LinestringM with M value is the length (in meters) from the start of the route - for example:               LINESTRING M (0 0 0, 10 0 20, 12 0 40, 20 0 50, 21 0 70)      Now I know the cable route is broken at _45 meters from start point_ , that would be in the middle of two points `(12 0 40)` and `(20 0 50)` (see the M value). With Postgis, how can I get the broken point coordinate - which should be `POINT (16 0)` in this case?


How do i draw a route on Quantum GIS between nodes?

I have a requirement to trace path on the points present in the shapefile which is loaded in Quantum GIS. This should enable the user to specify their own routes within the present network. **There will be existing points and lines present in the shape file. Whenever the user clicks on the points there should be a line drawn between the points selected and the route should be saved as a shape file.** **I have searched and found in many places of options like pgRouting and shortest path routing and that is not what i require.My need is simple and due to these options being available i m not able to find a simple solution.Its very much similar to this question :How to route railcars along railroads?. which does not yeild a clear solution.** I am a Newbie in GIS. Any guidance will be very much appreciated. Thanks in advance.


How to set field values as labels in QGIS2 Print Composer

I am mapping the distribution of 325 bird species, displaying circles in the grid cells each species occurs in. The data contains as many duplicates of the NAME as there are cells in which it occurs. I produce a map by filtering for a NAME then want to print the maps using the NAME field in a label as a caption below the map. I would like to have it generated automatically rather than write it in manually - in anticipation of being able to generate all species maps as an Atlas (that's another issue). So far, all I get by using `[% name %]` is `[name]`, not the species name. Any ideas of what I might be missing? Nothing I've found by searching has helped so far. Many thanks Niven


ArcGIS Android MapView with custom Spatial Reference SVY21

I have a MapView that is using a custom spatial reference of SVY21 from the server http://e1.onemap.sg/arcgis/rest/services/SM128/MapServer. The problem is the map's location service's accuracy circle due to the SpatialReference.getUnit causing an exception shown below               Could not update the active symbol     java.lang.RuntimeException: Invalid SpatialReference.          at com.esri.core.geometry.SpatialReference.nativeGetUnit(Native Method)     at com.esri.core.geometry.SpatialReference.getUnit(Unknown Source)     at com.esri.android.map.LocationService$a.k(Unknown Source)     at com.esri.android.map.LocationService$a$a.onLocationChanged(Unknown Source)      I really need the accuracy circle. Is there any solution to this?


Is it possible show the feature attributes dialog when I add a feature in code?

I want to know if it is possible to make the feature attributes dialog in Qgis show up when I add a feature by code using: `QgsVectorLayer::addFeature( tFeature )` This is available when I add features using Qgis' digitizing tools but not when I add features by code!


Image classification for a study site falling over two adjacent Landsat 8 scenes (qGIS or ArcGIS)

This question concerns work-flow rather than specific methodological problems. I will be using GRASS to classify Landsat 8 images to create a vegetation map for a study site that happens to fall over 2 separate Landsat scenes. I would really appreciate some advice about which route to take and I'm not able to find any other information about this elsewhere: 1.) Should I **FIRST** find a way to combine the 2 scenes together (which presumably would require combining the 2 scenes of each channel separately - would that be right?), and then cut out the small area covered by my study site and THEN start the classification process. OR 2.) Should I conduct image classification on the 2 adjacent Landsat scenes separately and then combine right at the end? (Although I'm worried this will take much longer, not to mention use up more memory - the study site I am working with is only 1400 km2) N.B. I am using a Mac (Mavericks OSX) and have access to both QGIS 1.8.0 and ArcView 9.2 to conduct any manipulation on this Geotiff files. Thank you very much for your time and advice!


How to use Plus tool in ModelBuilder?

I want to have the **sum** of at least 2000 raster combined to one raster layer using ModelBuilder ( **Plus tool** ). I think I'm missing something here. This model does not work. ![sample model builder flowchart](http://i.stack.imgur.com/ig9gd.jpg)


Memory layer or QgsVectorFileWriter Error?

Currently this excerpt of code will not save to file (as it is supposed to). The memory layer is visable and editable. Yet even when manually saving the layer, any priorly added features (geometries) are non-present. The variables 'FileDirectory._destination' and 'GetMaps._TOPOcrs' are proper as confirmed via the print command. It should be noted, that the error prompt for the FileWriter indeed reads an error with the value '1'. Any clues?                   try:             self.estabLocal()         except ValueError:             iface.messageBar().pushMessage("Error", "Please select a central folder location to save your model to!", level=QgsMessageBar.CRITICAL, duration=3)             return         try:             self.estabCRS()         except ValueError:             iface.messageBar().pushMessage("Error", "Please select a topographic map for your simulation!", level=QgsMessageBar.CRITICAL, duration=3)             return         pointstress = QgsVectorLayer("Point?crs=" + GetMaps._TOPOcrs + "&field=Well:string(3)&field=Discharge:string(6)&index=yes", "Points" , "memory")         pr = pointstress.dataProvider()         QgsMapLayerRegistry.instance().addMapLayers([pointstress])         PointStressFeat = FindDirectory._destination+"Points.shp"         writer = QgsVectorFileWriter(PointStressFeat, "CP1250", pr.fields(), QGis.WKBPoint, pr.crs(), "ESRI Shapfile" )         if writer.hasError() != QgsVectorFileWriter.NoError:             print "Error when creating shapefile: ", writer.hasError()         iface.setActiveLayer(pointstress)         pointstress.startEditing()         iface.actionAddFeature().trigger()


Why do I lose my printer settings when saving a copy of mxd from ArcPy?

Judging from the answer to this question, one cannot change the page size etc from ArcPy, so I set up my "template" mxds so my script can save a new copy of it, change data sources etc. This methods works perfectly for my A4 maps. However, my A0 mxds are losing the printer and page settings from the template. My `A0_template.mxd` is set to print from the plotter on A0 paper using the default printer settings. After the copy of the mxd is produced and I open it, the printer has somehow changed to the normal small printer with the paper size matching the default tray. I am simply using the `saveACopy(filename)` method which is part of the MapDocument class. Why would the correct settings for the A4 mxd be saved, but not for the A0?


Obtain address from PostGIS OSM data like Nominatim does?

We currently use the Nominatim service to do address lookups. We want to store the OSM data in our own servers as our usage gets higher. How do I calculate the address of a given coordinate, like how Nominatim does?


How can I split an osm file into separated files with poly files using osmosis?

I have got a country.osm file and I want to split it alongside of admin_level=8. I extracted previously ~3.000 poly files into a separated folder. I have batch script, that split them, if I give the input and the poly:               osmosis --read-xml %1 enableDateParsing=no --bounding-polygon file=%2 completeWays=no completeRelations=no clipIncompleteEntities=true --write-xml file=clc_%2.osm      **The question:** Are there any way to do this parallel and the input file read just once?


Ways to Speed Up Python Scripts Running As ArcGIS Tools

This is a pretty general question. I just wondering what tips and tricks GIS programmers have used to speed up arcpy scripts that you import into the toolbox and run. I work most everyday writing little scripts to help non-GIS users at my office process GIS data. I have found that ArcGIS 10.0 processing in general is slower than 9.3.1 and sometimes it gets even slower when running a python script. I'm going to list a particular example of a script that takes over 24 hours to run. It's a loop that tabulates the area of a raster in a buffer for each shape in the buffer. The buffer has about 7000 shapes. I don't believe it should run this long. A               while x <= layerRecords:              arcpy.SetProgressorLabel("Tabulating Row: " + str(x) + " of " + str(ELClayerRecords))         arcpy.SelectLayerByAttribute_management(Buff,"NEW_SELECTION", "Recno = " + str(x))                                  # Selecting the record         TabulateArea(Buff, "Recno", MatGRID, "VALUE", ScratchWS + "/tab" + str(z) +".dbf", nMatGRIDc)                          # Tabulate the area of the single row              arcpy.AddMessage ("          - Row: " + str(x) + " completed")         x = x + 1         z = z + 1      Before anyone says it, I have run tabulate area on the entire buffer, but it produces errors if run on more then 1 record. It's a flawed tool, but I have to use it. Anyways, if anyone has any ideas on how to optimize, or speed up this script, it would be most appreciated. Otherwise, if you have any speed up tricks for python, when used in ArcGIS, I'd love to hear them. Thanks in advance for any replies!


How to test whether two countries are adjacent?

I have an interesting question:   Is there a dataset by which I can use to test whether two countries are adjacent? For example, if inputs are Mexico and U.S., they are adjacent; if inputs are England and Russia, they are not. **Is it easy to produce an adjacency matrix for all the countries?** **Edited:** This question might be extended as, if we get a political boundary map, is there any ways to produce an adjacency matrix for all the political entities (e.g. Given a map of Florida, can a counties matrix be built to indicate the adjacency of all the counties)


How to draw curves in Flex Viewer?

I want to draw curves from a given point to another in Sample Flex Viewer like shown below. Any help? ![Routes](http://i.stack.imgur.com/Nk20s.png) ![enter image description here](http://i.stack.imgur.com/vt6Gm.png)


How to fix symbols/labels in a tiled WMS?

I have a postgis stored dataset that is consumed in a qgis project and served on qgis server. I noticed however that symbols and labels of point features that overlap 2 tiles of the map are only rendered in the tile that actually contains the feature. Is there some kind of buffering option or a configuration option for rendering to fix this? ![enter image description here](http://i.stack.imgur.com/GLGsn.png)


Connect JSON to GeoExt Application

I am returning a JSON to a GeoExt app and this is the code that is collecting the return,               var control = new OpenLayers.Control.Click({                          trigger: function(evt) {                           var lonLat = evt.xy;                           if (popup) popup.close();                           var  request = Ext.Ajax.request({                           url:"http://128.196.142.94:9000/info?lat="+ panel.map.getLonLatFromViewPortPx(lonLat).lat +"&lon=" + panel.map.getLonLatFromViewPortPx(lonLat).lon,                           disableCaching: false,                           success: function(response) {                                   var obj = Ext.decode(response.responseText);                                   console.log(obj);                                   console.dir(obj);                                   addToPopup(lonLat, response.responseText);                               },                           failure: function(response) {                                   console.log('server-side failure with status code ' + response.status);                                   addToPopup(lonLat, "Server side Failure");                               }                                    //callback: myhandler                      });     }});      I wish to obtain the returned JSON and then assign it to the popUp but whenever I try to connect to the App it is returning the following error,               OPTIONS http://128.196.142.94:9000/info?lat=55.44873046875&lon=-118.17626953125 405 (Method Not Allowed)     XMLHttpRequest cannot load http://128.196.142.94:9000/info?lat=55.44873046875&lon=-118.17626953125. Origin null is not allowed by Access-Control-Allow-Origin.     test.html:222server-side failure with status code 0      If you simply pass the url `http://128.196.142.94:9000/info?lat=55.44873046875&lon=-118.17626953125` you will obtain the JSON but for some reason there is failure that keeps occurring in the application. Is there CORS that can be used in the OpenLayers or GeoExt. I have tried using ExtJS JSON-P but I am not sure what should I enter in the 'callback' of,               <script src="http://domainB.com/users?callback=someCallback"></script>


Network Analyst directions window options

I have created a vehicle route problem with multiple route zones with one depot in each zone. Each zone has a pickup route that delivers to the depot in the zone. Some addresses have special instructions for where to pickup the item from their property as well as a phone number. My question is whether or not there is a way to make any special instructions or any other information like a phone number appear in the directions window to be printed out for the driver? I'm using Network Analyst in ArcGIS 9.3.1.


How to use ArcView GIS extension in QGIS or how to convert one to do so?

I really would like to know whether is a way to can work with ArcView GIS extensions at QGIS? or at least a way to be converted to plugins at QGIS?


Force a new AddDataDialog to open at a defined directory path

Is there a way to force the AddDataDialog to open to a specific path? I created a tool that when clicked on a specific feature will open the Add Data dialog, but I need it to open starting at a specified location. AddDataDialog doesn't have a property analogous to OpenFileDialog.InitialDirectory.


Editing with ArcGIS Silveright

I am using the ArcGIS Silverlight API to edit data in an ArcSDE SQL Server geodatabase. I have been told that any editable feature classes must be registered as versioned in ArcCatalog before they can be edited via the Silverlight API. Does anyone know why this is the case? From what I understand, you can not use the Silverlight API to edit versioned data in a feature class, so I'm not clear on why you would have to register the feature classes as versioned. Thanks, Seth


Specify a starting location with IgxDialog

I would like to know what is the easiest way to set the starting location of a data dialog to a IWorkspace? My problem is that basically I do not want a connection (.sde file) for the users. I want to set the starting location to a workspace connected via IPropertySet. Is it possible? This is related: Force a new AddDataDialog to open at a defined directory path but does not help. EDIT: following Kirk's suggestion, here's a code snippet. This does not work.                   public static IDataset OpenDatasetDialog(IWorkspace workspace,int handle)         {             var workspaceName = (IWorkspaceName)((IDataset)workspace).FullName);             var gxDatabase = new GxdatabaseClass();             gxDatabase.WorkspaceName = workspaceName;                  // we need to set a different variable as .NET runtime does not allow ref parameters with casts             object gxObj = (object)gxDatabase;                  IGxDialog openDialog = new GxDialogClass();             openDialog.RememberLocation = false;             openDialog.set_StartingLocation(ref gxObj);             var filterColletion = (IGxObjectFilterCollection)openDialog;             filterColletion.AddFilter(new GxFilterTablesAndFeatureClassesClass(), true);             openDialog.AllowMultiSelect = false;             IEnumGxObject datasets = null;             IDataset dataset = null;                  if (openDialog.DoModalOpen(0, out datasets) && datasets != null)             {                 IGxObject gxObject = datasets.Next();                 if (gxObject is IGxDataset)                 {                     IGxDataset gxDataset = (IGxDataset)gxObject;                     dataset = gxDataset.Dataset;                 }                 else                     throw new ArgumentException();             }             openDialog.InternalCatalog.Close();             return dataset;         }


Flex popup from outside resource weather warnings

I am using the USGS service of NHSS NOAA weather warnings in an application. I have the service displaying in my map using   http://rmgsc.cr.usgs.gov/ArcGIS/rest/services/nhss_weat/MapServer   The service shows that identify is supported.   The question is how do I get a popup from this service with the warning type and other info?   If I try to use the same service I get no popup, If I point it at my own server I get a popup with no information.   I have also tried using the atom described here but not sure of the exact usage in this context. Any help would be appreciated.


How to put GeoServer layers behind authentication?

I want to publish map layers using OpenLayers and GeoServer. While this is easy enough to set up, I'm a bit lost when it comes to authentication possibilities. What I want is to create different users that will, by logging in with their own user name and password, be directed to their own map window. How can this be securely done so that none of my users will be able to access each other's layers? I have understood that it is quite easy to "hack" the setup by, for example, making queries to the GeoServer where one simply asks it to list all the available layers. What are my options? I suppose apache could be of help when setting it up in my ubuntu server? Any tutorials on this subject?


Rotate while in Layout View

I am using ArcGIS 9.3 without any extensions. How do you get your map to rotate? I've done it before but it has been awhile and now the "rotate" does not seem to function now. I've tried using while in edit and non-edit. Any suggestions?


How to find out Why a QGIS plugin is rejected for Version 2

I have a selection of python plugins (one rather large one) that I have developed under Qgis 1.8 and that I am migrating to qgis 2.0. The error that the Plugin Manager gives me when I load these plugins (This plugin was developed under 1.x) is rather unhelpful. It also arises in situations that have nothing to do with the API Version. For instance, a Network error yesterday introduced a stream of random characters into a few of my code files (Fun!) and this led to the Plugin Manager telling me my plugin was old. Eventually, after pain, I found out a tab and a dollar sign had been inserted in front of a comment in my metadata.txt, and that was the Problem. Screwed up, yes. Old, no. Is there any way to see a more detailed error when QGIS 2 rejects a plugin as 1.x? If not a proper trace, then at least what class, or even what file was rejected? Anything? The instructions regards API changes at the Qgis Hub Wiki are lovely and I sleep with them under my pillow, but they are not as lovely as a simple error message. Qgis 2, by the way, is beautiful. Thank you.


How to create .dem file in QGIS?

I have an elevation data in .asc file, but need to convert it to .dem file to later use in COMSOL. However, in Raster--> Analyzes-->DEM I got just error message: ERROR 1: DEM Default Setup: Importing source dataset projection failed What is wrong? regards, L.


Possible to run "Merge shapefiles to one" algorithm from Python script?

Developing a python script in QGIS 2.2, I can run algorithms from the **Processing** plugin fine such as those from _QGIS geoalgorithms_ and _SAGA_. Can you do the same for the algorithm " _Merge shapefiles to one_ " from the **Vector** menu and how can you call it from the Python console? EDIT: The " _Merge shapefiles to one_ " algorithm is from the Ftools plugin. Found a link which might be useful: Is there a way to access QGIS plugins in Python?


Drawing text and line-dasharray

Is it possible with TileMills and CartoCSS? ![enter image description here](http://i.stack.imgur.com/WmKWA.jpg)


ArcGis Flex Viewer Helloworld widget doesn't display

I'm trying to create a new widget, but it doesn't show. If I comment out the other widgets in the widgetcontainer in the config.xml file, they don't go away. I don't understand why the changes I make don't show... config.xml file:               <widgetcontainer layout="float">             <widget label="Bookmarks"         left="430" top="90"                     icon="assets/images/i_bookmark.png"                     config="widgets/Bookmark/BookmarkWidget.xml"                     url="widgets/Bookmark/BookmarkWidget.swf"/>             <widget label="Louisville Police" left="590" top="280"                     icon="assets/images/i_police.png"                     config="widgets/Query/QueryWidget_Louisville_PoliceStations.xml"                     url="widgets/Query/QueryWidget.swf"/>             <widget label="Search" left="80" top="280"                     icon="assets/images/i_search.png"                     config="widgets/Search/SearchWidget_Louisville.xml"                     url="widgets/Search/SearchWidget.swf"/>             <widget label="Earthquakes (GeoRSS)" left="410" top="280"                     icon="assets/images/i_rss.png"                     config="widgets/GeoRSS/GeoRSSWidget.xml"                     url="widgets/GeoRSS/GeoRSSWidget.swf"/>             <widget label="Draw" left="60" top="400"                     icon="assets/images/i_draw2.png"                     config="widgets/Draw/DrawWidget.xml"                     url="widgets/Draw/DrawWidget.swf"/>             <widget label="Print" left="390" top="400"                     icon="assets/images/i_print.png"                     config="widgets/Print/PrintWidget.xml"                     url="widgets/Print/PrintWidget.swf"/>             <widget label="My first widget"                     icon="assets/images/i_widget.png"                     config="widgets/Samples/HelloWorld/HelloWorldWidget.xml"                     url="widgets/Samples/HelloWorld/HelloWorldWidget.swf"/>             <!--             <widget label="My first widget"                     icon="assets/images/i_widget.png"                     config="widgets/Samples/HelloWorld/HelloWorldWidget.xml"                     url="widgets/Samples/HelloWorld/HelloWorldWidget.swf"/>             -->         </widgetcontainer>      I click the "Run" button, the application appears in my browser, but that widget isn't there. I've confirmed that the widget is in that directory, I've added the module. I'm following the directions here: http://video.arcgis.com/watch/539/arcgis-viewer-for-flex-creating-a-custom- widget http://resources.arcgis.com/en/help/flex- viewer/concepts/index.html#/Creating_widgets/01m300000010000000/ ![enter image description here](http://i.stack.imgur.com/GzNX1.png) ![enter image description here](http://i.stack.imgur.com/GI1US.png) ![enter image description here](http://i.stack.imgur.com/yRlVw.png)


How do I use wildcards to work around the GDAL character limit?

I am attempting to use the 'Merge' & 'Create VRT' tools in QGIS 2.0 on a folder with 3000+ raster tiles in. However whenever I try to perform this operation the following error appears: **'The process failed to start. Either the invoked program is missing, or you may have insufficient permissions to invoke the program.'** The process works when I select only a few tiles, so I can only assume it's due to a character restriction due to the number of files? Upon searching the website for a solution I came across the concept of using 'Wildcards' to get around this problem: Having problems creating VRT in QGIS from png files However I cannot find any clear explanation on how to use them, and I don't have much experience with using OSGeo4W outside of the QGIS GUI. I would appreciate any help or advice.


what is Network analysis and how to find shortest path

I don't know about Network Analysis.Please any one help to understand the Network Analysis with example. General Process flow and procedure is enough to learn my self Thanks in Advance.


I am having Problem with my geoserver using URL

I am having problem when I am Requesting following URL in GeoServer /geoserver/wfs?REQUEST=GetFeature&VERSION=1.1.0&typename=Final&propertyname=name I am getting following error http:///geoserver/schemas/ows/1.0.0/owsExceptionReport.xsd"> java.lang.RuntimeException: Unable to obtain connection: Connection refused. Check that the hostname and port are correct and that the postmaster is accepting TCP/IP connections. Unable to obtain connection: Connection refused. Check that the hostname and port are correct and that the postmaster is accepting TCP/IP connections. Connection refused. Check that the hostname and port are correct and that the postmaster is accepting TCP/IP connections. Connection refused: connect Can any one help in Resolving this issue


how to make a line from a set of gps points?

As a part of my research I need to create the trajectory of respondents according to the their collected GPS data (including X and Y). I am familiar with python programming, but I am a beginner in ArcGIS. I want to write a script to automate this process I really appreciate your help and comments.


Where does ArcMap GUI get font size and docking guide positions from?

In the screenshot below, ArcMap 10.1 appears to be making some arbitrary choices about the font to use for:   1. The Layers data frame name in the Table of Contents (larger font than expected) - **NOW FIXED**   2. The Main Menu pulldowns (smaller font than expected) - **_NOW RESOLVED_** Can anyone tell me where to find a setting that will enable me to make these the same size as the font used for the Main Menu items, please? ![enter image description here](http://i.stack.imgur.com/e6AIc.png) UPDATES: I finally resolved this by visiting **Control Panel\Appearance and Personalization\Display** and **resetting the font size I had set there from 150% to 125%**. This resolved all. I was also pursuing the mismatch in font size between Menu items and pulldowns via the ArcGIS Discussion Forums. The problem also affected the docking guides used to help resize windows (see graphic below). ![enter image description here](http://i.stack.imgur.com/FmMwX.png)


Why does QGIS scalebar show a wrong scale?

I'm using QGIS 1.7.1 and have added a scalebar. The scalebar usually displays in degrees, unless I switch the map's projection to meters, in which case the units are display in _cm_. This is not correct. Any thoughts on how I can get this in reasonable units? Thanks!


Leaflet Event Craziness: jQuery click events don't propagate to L.DomEvents

I noticed that calling `.click()` on jQuery selected DOM elements in Leaflet ( such as anchor icons ) will not do anything -- the click event isn't propagated to the events created with L.DomEvents. For example, go over to bboxfinder and try to trigger the zoom out button's click event:               $('a[title="Zoom out"]').click();      Nothing will happen, the map won't zoom out. I jumped down the rabbit hole to figure out the why(s) and ended up with more questions than answers. My two favorite questions are listed below. I'm hoping the answers for each question address two parts:   1. Why is this behavior happening? Be as technical as you can so people can learn!   2. What are the benefits of this behavior? This is a well-written, abstracted library. I assume they made these choices for reasons that are not easy for someone like me to see. Any guesses what those benefits are? For example, maybe these choices work well with mobile devices. **Questions** **1.** Why are the jQuery events not propagating to DOM elements? I can call click on the unwrapped jQuery DOM element ( see below ) and fire the event, but i assume jQuery would pass this event down for me:               // this works fine, notice we're grabbing the raw element     $('a[title="Zoom out"]')[0].click();      **2.** When i was digging around in the Leaflet source I noticed button events were created/bound with this code workflow below ( it also happens in 3rd party packages like L.Draw ). The question is why are we using all these events on a single button and how are they working together? And really, why is it like this? I've pulled the code below from this part of the source:                   L.DomEvent          .on(link, 'click', stop)          .on(link, 'mousedown', stop)          .on(link, 'dblclick', stop)          .on(link, 'click', L.DomEvent.preventDefault)          .on(link, 'click', fn, context)          .on(link, 'click', this._refocusOnMap, context);


select by attribute on shapefiles with open-source python modules

I usually use ArcGis 9.3.1, but when running a python script that includes Select_by_attribute (from gp), the script runs very very slowly. Is there a way I can perform a select-by-attribute on a opensource-python-module (such as shapely?)? Will it be faster?


Find Nearest features from an input location

I wanted to find the nearest features (it can be point, line and polygon layers) from an user input or search location and identify 5 nearest locations along with their distances. Have anybody worked on this link below. If so, can anyone point me to the right sample application, please. https://developers.arcgis.com/javasc...arest-amd.html https://developers.arcgis.com/rest/a...nd-nearest.htm On the above link, there is no sample application or code to show how this works and also it can only query against two feature layers. I wanted to query against multiple layers. Secondly, there should be a way to do without using FindNearest class, as this class has only been introduced with V3.7. If this is the case, then how GIS developers across the world would have created application without this class. I think, some one should really have an answer to this.


Changes in esri config in VS2010 has no effect

I've created a DockableWindow in my Add-In (vs 2010, c#, arcobjects 10 sp2) and the esri config was updated with the following:               <DockableWindows>         <DockableWindow id="CTArcMapAddin_AreasDockWindow" class="AreasDockWindow+AreasDockWindowImpl" caption="Areas" image="Images\AreasDockWindow.png">           <InitialPlacement height="300" width="250" state="unpinned" position="left" />         </DockableWindow>       </DockableWindows>      I run the Add-In and the dockable window shows up on the left side... Then when changing the "position" value to "right" and running the addin again it's still positioned to the left...why does that change not show through...a bug maybe?


Calculate ANOVA using GIS

It's possible to calculate a statistic like ANOVA using GIS. I'm using ArcGis, but there's any GIS software to calculate it? Thanks,


Run a model from python problem

I am using arcGIS 10.0 SP5 I have the following code :               import arcpy          arcpy.ImportToolbox(r"C:\Documents and Settings\Derek\FRAME\FORTEST.tbx", "Derek")          arcpy.TEST_Derek()      The model is names TEST The toolbox alias is Derek When I try and run this I get an error               AttributeError: 'module' object has no attribute 'TEST_Derek"      What I want to do is to run the model TEST from the toolbox FORTEST


Defining SQL 2008 store within GeoServer 2.3

I want to define a SQL server 2008 store within Geoserver 2.3. I installed the sql Extension. After defining the different attributes, I get the following error message :               Error creating data store, check the parameters. Error message: Unable to obtain connection: Cannot create PoolableConnectionFactory (The TCP/IP connection to the host localhost, port 1433 has failed. Error: "Connection refused: connect. Verify the connection properties, check that an instance of SQL Server is running on the host and accepting TCP/IP connections at the port, and that no firewall is blocking TCP connections to the port.".)      The data I entered : host -localhost/sqlr2 port - 1433 database - the database name within SQL scheme - dbo user & password = A geoserver user i defined within sql Thanks in advance Lauren


QGIS Quickly switching overlapping raster images

I'm new to QGIS, and haven't used arcGIS in about 5 years, however I've just seen someone using arcGIS and was able to quickly switch between overlapping rasters when checking registration. The top half of the screen was one raster, and the bottom was another, and the mouse was used to move up and down in real time switching between the two (I'm probably not explaining this very well). The images I'm working with are very large (2x1GB geotiffs) and of the same area. There is a good 5second lag between switching layers which makes it difficult to see if there is good alignment. Does this functionality exist in QGIS, to quickly switch rasters in real time?


Watershed layers with names for each one for Europe or Sweden

I am working on freshwater species and I would like, make a list of each watersheds used for each species. So I am looking for a watershed maps which contain name of the watershed. I found some good layers, but, they don't get names for each basins, only for the bigger ones. I would like to know if someone can give me advices, links or good sources? I need mainly the cover for Sweden at least, or Europe. I searched this site, but I found most of the elements to create watersheds. Regards


Is it possible to reproject KML using ESRI's Silverlight API?

My question is similar to this one, but for Silverlight. I need to reproject an entire KML file to NAD83, add it to the map with all the layers, symbology, etc. in tact. Is there anyway to do this with the Silverlight API? thanks


Does anyone display OpenStreetMap through MapServer, ArcGIS for Server or Geoserver?

I use Windows OS. The data (OpenStreetMap) are stored in a postgres/postgis database. How do I display this data correctly? **Update:** Now I use: geoserver-2.3.5.exe, postgresql-9.1.9-1-windows-x64.exe, postgis-pg91x64-setup-2.0.3-2.exe, osm2pgsql **Update2** I lookd at the source for "osminabox". https://git.hsr.ch/index.html#repositoryBrowser;26067c88-f4a8-46e8-be0d-95c1a5091c94;null;Code/osm2gis/geoserver/geoserver_data/styles   1. But I find many sld styles (I have only 4 layers: planet_osm_polygon, planet_osm_line, planet_osm_point, planet_osm_roads).   2. There are a lot of styles for points: osm_poi.sld, osm_building.sld, osm_pow.sld. I don't have this layer.    3. Also I find strings such: "img/recycling.png"  So, If I use intranet system I can't work? **Update 2 continue:** I try to use https://github.com/fotographs/OSM-Styles- for-GeoServer/tree/master/sld. But I cannot find sld styles for the points. https://help.openstreetmap.org/questions/7521/geoserver-style-files - Would it help?


QGIS 2.01 - raster transparency not working properly

![Trasnsparent raster image with underlying vector layer](http://i.stack.imgur.com/TlqM4.jpg)I have loaded a color raster image (aerial orthoimagery) into a QGIS project. The layer is visible but very transparent (in the enclosed image the vector polygons are below the imagery), but when I check Properties/Transparency the transparency is set to 0%. The imagery is apparently color infrared. In QGIS 1.8 I simply went to the Color Map dialog and set the R-G-B color map to 4-1-2 respectively and the display was in perfect true color. I cannot figure out how to do the same with 2.01. I have looked at the options under multiband and singleband pseudocolor and neither seems to give me the option to reset the color map as I was able to do with 1.8. You can see that the image is displaying in true color, but highly transparent. I did not change any color map settings. When I view an older project that created with 1.8 with 2.01 imagery from this same data set displays perfectly, so it is not a problem with the source data. However, when I look at the Style settings for the older project with 2.01 I cannot find any differences between that project and the Style settings for the new project. I don't konw if this is a color map issue or a trasparency issue, but I would simply like to get the image ot display properly.


How to keep all the unique columns of the shapefiles after merge?

Writing a script where the user is asked to select a directory which contains Folder_1 and Folder_2, both containing shapefiles and run various algorithms. The basic logical process of the script is:   1. Create a standard grid   2. Look into Folder_1 and clips a shapefile onto the grid   3. Create a unique column in the attribute table for the newly clipped shapefile and save into the Output folder   4. Loop through any remaining shapefiles and repeat Step 2 & 3   5. Repeat Step 2-4 for Folder_2 I then want to merge everything in the Output folder using **saga:mergeshapeslayers** , however, this only takes the attribute columns from the first selected layer. I would like all the unique attribute columns to appear. Can this be done easily? I've attached the code.               ##Test=name     ##Select_directory=folder     ##Cellsize=number 1000     ##Result=output vector          import os     import glob          #   Folder path of output shapefiles     path_res = "C:\Users\Test\Output\\"          #   Select required directory     os.chdir(Select_directory)               def rank():          #   Create grid with defined parameters         outputs_1=processing.runalg("qgis:creategrid", Cellsize, Cellsize, 24108, 18351.157175, 258293.802316, 665638.226408, 1, 'EPSG:7405', None)         outputs_2=processing.runalg("qgis:fieldcalculator", outputs_1['SAVENAME'], 'ID', 1, 10, 0, True, '$rownum' , None)         outputs_3=processing.runalg("qgis:deletecolumn", outputs_2['OUTPUT_LAYER'], 'longitude', None)         outputs_4=processing.runalg("qgis:deletecolumn", outputs_3['SAVENAME'], 'latitude', None)          #   Set directory and search for all .shp files in Folder_1         os.chdir("C:\Users\Test\Folder_1\\")         for fname in glob.glob("*.shp"):                #   Clip shapefiles files within Folder_1 to Grid and save files to Output folder                  outputs_5=processing.runalg("qgis:clip", outputs_4['SAVENAME'], fname, None)             outputs_6=processing.runalg("qgis:fieldcalculator", outputs_5['OUTPUT'], 'Rank_' + fname, 1, 10, 0, True, 1 , path_res  + "/"+ fname)          #   Set directory and search for all .shp files in Folder_2         os.chdir("C:\Users\Test\Folder_2\\")         for fname in glob.glob("*.shp"):           #   Clip shapefiles files within Folder_2 to Grid and save files to Output folder                 outputs_7=processing.runalg("qgis:clip", outputs_4['SAVENAME'], fname, None)             outputs_8=processing.runalg("qgis:fieldcalculator", outputs_7['OUTPUT'], 'Rank_' + fname, 1, 10, 0, True, 2 , path_res  + "/"+ fname)          #   Paths of the shapefiles in the Output folder with list comprehension         output = [shp for shp in glob.glob(path_res + "*.shp")]          #   Merge the shapefiles         processing.runalg("saga:mergeshapeslayers", output[0], ";".join(output) , Result)               if Select_directory:         rank()     else:         pass


Best Way to Get a List of Selected Rows for a Specific Layer using Layer Object in ArcGIS

What is the fastest/cleanest way to get a collection of rows that are selected in a certain layer? For some background, I am trying to hide selected features in a specified layer.   * The solution involving the use of **SearchCursor** uses the layer name, but my setup may involve duplicate named layers   * Another solution I found copies the selected rows into a new layer, which may be time costly for my setup. Is there a way to get the selected rows by using a layer object which is provided by a script parameter?


ArcGIS Javascript API LabelLayer - when clicking a label is there a way to know which feature is represented by that label?

When adding labels via addFeatureLayer it is impossible to know which feature was clicked? The "click" callback for a LabelLayer provides the label graphic but no reference to the FeatureLayer graphic this label represents. Is there a way of determining the graphic from the FeatureLayer that this label graphic represents? https://developers.arcgis.com/javascript/jsapi/labellayer- amd.html#addfeaturelayer


Import processing in qgis 2.0

I know that this question has already asked but no concrete answer has given!! My environment is Ubuntu 12.04 and I'm using eclipse-pydev as IDE. Thanks in advance for helping me.


WFS service Server Error in '/ArcGIS/services' Application

I have a SSL secured ArcGIS Server 9.3.1 machine on W2k8. I set up a service with WFS. Every time I try to request the capabilities, I get the error: Server Error in '/ArcGIS/services' Application Runtime Error Description: An application error occurred on the server. The current custom error settings for this application prevent the details of the application error from being viewed remotely (for security reasons). It could, however, be viewed by browsers running on the local server machine. Details: To enable the details of this specific error message to be viewable on remote machines, please create a tag within a "web.config" configuration file located in the root directory of the current web application. This tag should then have its "mode" attribute set to "Off". Notes: The current error page you are seeing can be replaced by a custom error page by modifying the "defaultRedirect" attribute of the application's configuration tag to point to a custom error page URL. In the IIS log I see a 500 type error (Internal server error).


how to create an outline polygon (footprint) from raster data with postgis?

how can I create a polygon, containing the actual outlines (aka footprint) of a raster file with postgis? I just need the outlines of the data containing parts of the raster. ST_PixelAsPolygon creates a complex polygons - but I just need the footprint. There is a similar question here How do I create a shapefile showing footprints of Rasters?, but the answers do not cover postgis. D3


resampling errors

I have to reproject my images, but it contains some resampling errors. Any one has any idea on significance of resampling error. Are they neglectable? Is there any good documentation to refer to? I'll do it by GDAL. In my case source coordinate system is WGS84 Lat,Lon (ESPG=4326) and target is MODIS Sinusoidal.


How to color-balance false-composites of pan-sharpened Landsat images in GRASS?

Whenever I try to pan-sharpen composites of some Landsat images in GRASS using `i.pansharpen`, `i.fusion.brovey` or the IHS sharpening method, the output will have some or all of the following characteristics:   * the composite color is in a different hue compared to the un-sharpened composite   * the brightness level is messed up   * the entire composite went all-white/all-black (when using images pre-processed to top-of-atmosphere reflectance or surface reflectance corrections in `i.landsat.toar`) I've also tried **all** of the following; but the colors/brightness remained the same or turned even worse:   * Applied `i.landsat.rgb`, before-and-after the pan-sharpening process   * Played with the `-f` or `-p` flag in `i.landsat.rgb`   * Tried `r.colors` to edit the color table to grey/grey255/grey.eq   * Tried `i.pansharpen` using all Brovey/IHS/PCA methods   * Played with the `-l` flag in `i.pansharpen` to rebalance the blue-channel The GRASS GIS manual did explained on how to perform pan-sharpening and color- balancing, but I can't figure out how to combine both processes in a concurrent workflow. I suspected that this is due to my poor understanding of color-tables, color-histogram, etc. in GRASS.. So, can someone explain to me - **how do you tackle color-balancing problems when dealing with Landsat images after image-processing in GRASS?** Can you share with me your favorite workflow/methods? Many thanks for any feedback!


How to Create and Rotate Raster Dataset using ArcObjects?

I'm following the "How to create a raster dataset" tutorial: http://help.arcgis.com/en/sdk/10.0/arcobjects_net/conceptualhelp/index.html#//000100000464000000 The questions are: 1) How can I create the raster in a file geodatabase? (The example given is incomplete.) 2) How can I rotate the newly created raster? (Standalone TIFF or a geodatabase raster.) I know there is IRasterGeometryProc3, but I can't make it work. Any samples / help would be appreciated. Cheers.


What are use cases of Esri's CityEngine?

I would like to know the use cases of Esri's CityEngine, and what distinguish this tool from the other standard GIS like ArcMap and QGIS? and in which fields could be used ?


How can I edit lines containing M values with the ArcGIS JS API?

I am working with the ArcGIS JavaScript API and want to move and edit the vertices of a polyline analog to the following sample: http://help.arcgis.com/en/webapi/javascript/arcgis/demos/ed/ed_toolbar.html This works fine as long as the polyline does not contain M Values. If it does contain M Values I get the following error in firebug. Errorcode:-2147221399, description: „Geometry does not have M Value(s), hence cannot update the feature." What am I doing wrong? Here’s my code:               function initEditing(layers) {         var myMap = this;         var firePerimeterFL = myMap.getLayer("firePerimeterFL");              var editToolbar = new esri.toolbars.Edit(myMap);         dojo.connect(editToolbar, "onDeactivate", function(tool,graphic,info) {           if(info.isModified){             firePerimeterFL.applyEdits(null, [graphic], null);           }         });              var editingEnabled = false;         dojo.connect(firePerimeterFL, "onDblClick", function(evt) {           dojo.stopEvent(evt);           if (editingEnabled === false) {             editingEnabled = true;             editToolbar.activate(esri.toolbars.Edit.EDIT_VERTICES , evt.graphic);           } else {             editToolbar.deactivate();             editingEnabled = false;           }         });       }


ESRI Flex Viewer 2.4 - How Can I Remove The Header (or UI Component that holds the HeaderController Widget)

I am trying to make some modifications to the general layout of the Flex Viewer. As part of these modifications, I am trying to remove the Header altogether. I removed the HeaderControllerWidget from the config.xml - thinking the map should now occupy the entire browser. But now I just get empty white space where the HeaderController was, leading me to believe some UI component must be defined elsewhere to hold the HeaderController widget. I have searched through the code trying to find this component with no luck. Does anyone have any suggestions as to how I can remove Header altogether? I want the map to occupy the entire browser.


Problem converting KML file to .lyr

I have a problem in converting KML files into lyr ones with ArcGIS 10. I do all the steps perfectly but the conversion fails and the red check mark appears. It is surprising that when I repeat the steps for more than 10 or 20 times eventually it works and the green check mark appears!   I have the same problem with different KML files. what solution do you suggest to overcome this problem? thanks a lot


How to find line length in QGIS while editing?

I am using QGIS (before, I used ArcView). I would like to know if there is any way during the editing phase of a vector (line or polygon) to know the length of the line (total and partial) that is being created. Thanks.


Bring Feature Selection To Top

I'm using IFeatureSelection.SelectFeatures() to select parcel polygons in my parcel layer. It works fine until I create a ClassBreaksRenderer to apply a color ramp based on a field in the feature. The selection proceeds without error but does not draw on top of the renderer. Is there a way to ensure that the feature selection draws on top?


Properly deal with broken project data source Analyze Errors when publishing Geoprocessing Service (Python / ArcGIS Server 10.2)

When publishing a geoprocessing service, I noticed that arcgis does not like to see:               os.path.join(os.path.dirname(__file__), 'projections')      This line seems to throw an **broken project data source** while analyzing the service for publishing to arcgis server. I'm guessing its because the only the main .py file will be copied into the services directory and Esri is giving me a heads up that the data source is going to break because 'projections' folder will NOT be copied over. In my case, I also have some unit test which reference a test_data folder and shouldn't affect the tool when running off the server. Unfortunately I still can't publish. To publish, I:   1. Moved my relative path references to a separate module named service_config.py   2. imported the path constants from service_config.py   3. Added a .pth file to Lib/site-packages to allow the services directory version of the script to still have references to the projection folder without needing it to be copied to the arcgisinput directory. I feel like my solution is a hack. Any suggests on how to improve this tool configuration? My code can be found here: https://github.com/brendancol/esri-extract-data-where


Why should I use the AMD style deceleration in ArcGIS JS API?

It seems to me like this:               var point;     require(['esri/geometry/Point'], function(Point) {         point = new Point();     }      does the same thing as:               var point = new esri.geometry.Point();      Is there any reason to opt for AMD-loaded constructors?


How to tranform a column field into a list using Field Calculator

I'm trying to make a simple field calculation with ArcGIS 9.3: New field = Old field - Old field(first value) which in Calculate Field with python code should be Expression: !Old field! - first Code Block: list= [AngOriz] first = list[0] The error is ERROR 000539: Runtime error : name 'AngOriz' is not defined Failed to execute (Calculate Field). How could I tranform a column field into a list? Ivano


When should you NOT use a spatial index?

I am asking this because I was mainly working with Oracle but for the past year I've been doubling with PostGIS and SQLServer 2008. Most spatial functions in Oracle won't work without a spatial index returning the ORA-13226 error: > 13226, 00000, "interface not supported without a spatial index" //*Cause: > The geometry table does not have a spatial index. // *Action:Verify that the > geometry table referenced in the spatial operator has a spatial index on it. To me this makes sense. You run a spatial query=you must have a spatial index. But as far as I understand, neither PostGIS not SQL Serve require this. PostGIS even seems to have functions (_* e.g. _STContains) that EXPLICITLY won't use the spatial index. So the question is- are there any cases where you should NOT use a spatial index?. Not necessarily whether is a 'take it or leave it' approach i.e. it won't make any difference, but where NOT using the spatial index will impove performance? To me, the last sentence is a contradiction in terms but otherwise why PostGIS would provide these functions?


Get features WMS under a polygon in layer-vector

I have the buffer of a polygon in a layer-vector, I want to know what features of a WMS layer are covered by the buffer, as the SLDSelect Control but with a polygon drawn. any ideas? thanks in advance


Calculate stop times/distances of transport lines - tools and procedure?

I need to calculate stop times of bus lines (extrapolation) I have:   1. Shape file with all stops   2. shape file with control points where bus time was taken   3. Shape file with all bus routes   4. csv file with times in first stops on the route, final stops and control points where time was taken. All of shp files have an attribute "route number" which identifies to which transport route they refer. I need to find a tool (and procedure) that helps me extrapolate times in each of the stop on the route, using the departure time, arrival time and times at control points on the route (which are not stops). I suppose I actually need to calculate distances first and then times based on those distances. I am unix user so OS tools would be preferable Any suggestion how to crack this? Attributes described below. Thanks a lot for any hint. Karel * * * **Stops shp file** Description: Locations of bus stops Attributes: Cod_ubic_p - stop location Code_variant - bus line (referenced in _Bus_ shp file) Ordinal - order corresponding to the stop in the journey of the variant (bus line) * * * **control_points shp file** Description: Locations of control points Attributes: Cod_ubic_c - control point code Code_variant - bus line type (referenced in Bus_lines) Ordinal - order of given control point in the path of the variant (bus line) Desc_locat - Description of Control Point location * * * **Bus_lines shp file** Description: Routes of bus lines Attributes: Cod_line - Transportation line code Desc_line - Description of the transport line (eg 112, D16, etc) Ordinal - Correlative number of the sub-line of the line Code_subline - Code of the transportation sub-line Subline Desc - Description of the sub-line Cod_variant - Code of the route (to link with _Stops_ shp) Desc_variant - Type (Variant) description (A = Ascending, D = descending) * * * **csv sheet _Timetables_** Description: Timetables Attributes: COD_LINE - Transportation line code CODE_DAY_TYPE - Code for the day of the week LEAVES_FROM - Origin of the bus route DEPARTURE - Depature time DESTINATION - Destination of bus route ARRIVAL - Arrival time DESC_LOCAT - Description of Control Point (Stop) location COD_UBIC_C - control point code (referenced in *control_points* shp file) TIME - Time at control point location ORDINAL - Correlative number of the sub-line of the line DESC_VARIANT - Type (Variant) description (A = Ascending, D = descending)


How to change projection and show latitude / longitude in QGIS?

I am new to QGIS and pretty ignorant with the subject, but I want to be able to see the latitude / longitude of a point (with a mouse click) in Quantum GIS. I imported a shapefile (roads of finland from here : http://www.maanmittauslaitos.fi/en/node/6567 ), and there is a box a the bottom that show the current coordinate in unknow units to me : 521636,6947168, same the coordinate capture don't show lat/long. The coordinate system seems to be WGS84 - 4326. I am sure it is a stupid question, but thanks for your help.


Performing statistical analysis of buffers?

I am doing a buffer analysis of incinerators and looking at the output areas within those buffers, is there a way of doing a statistical analysis of several buffers to see if there is a significant pattern between the data within the buffers?


GeoExt Legend Modification To Include Feature Count

GeoExt has Legend capability but is is quite basic which shows only feature types like in example http://api.geoext.org/1.0/examples/legendpanel.html. How could it be modified display legend with count(no of values)? ![enter image description here](http://i.stack.imgur.com/ujlye.png) I am using **_geoserver_** to serve map and using **_geoext_**. The solution I am looking for is to add this feature into legend panel (may be through modification of javascript source) and not make my own solution. Any path, direction code for this would be appreciated.


Re-enabling a web service API call that has stopped working unexpectedly

so we have are running an instance of ArcGis Server, and have a handful of services running. These were set up by someone who is no longer working with the company. Some of them are built into applications that have been functioning just fine. Recently an application we were running though had one of the API calls that the page uses stop working. It was pulling in lat/long from user entry and sending it to the service using GET, to determine what boundary the point fell in. About a week ago and a half ago, this just stopped functioning. I don't think anything has changed, or at least I don't know who could/would have changed anything, and I can't find any evidence of change either. The funny thing is that all of our other services are working fine, including one that does the same kind of point-in-polygon search, but on a different layer. I can't even find the underlying code to access this lookup (if there is any). I have the URL, and I have the service entry in the arcgis server manager, but I don't know what files configure that URL that I could look at. Any thoughts? Generally I could use help figuring out why this may have stopped working unexpectedly, And perhaps how to either turn it on, or whether we need to build this all over again.


Is it possible to unique value style a raster in QGIS?

In ArcGIS one can load a raster layer and choose an option to style the raster based on unique values. This is helpful if I have a raster of discrete values that are not evenly distributed. For example, I have a land cover map that has values ranging from 1-12, 20-30, and 250. In this case a continuous or equal interval classification will not uniquely classify the raster. Is possible to do in QGIS directly or with a plugin?


Does OpenLayers3 support feature editing?

For OpenLayers2, there are nice examples which show how to create vector features and edit them (adding vertices, moving them, ...), e.g. http://dev.openlayers.org/releases/OpenLayers-2.13.1/examples/modify- feature.html Is it possible to reproduce this functionality in OpenLayers3? On the official website, I was only able to find an example which shows how to create vector features but not how to edit them: http://ol3js.org/en/master/examples/draw- features.html


Spatial ETL or some other tools?

and, first of all, sorry for my mistakes, my English is not so good. I need to read geographic data from an ORACLE spatial database in Gauss Boaga SRS, convert them in WSG84 (EPSG 32632), e put them in another ORACLE spatial data base. I've tried to use GeoKettle 2.5 but I think that there is some bugs because it seems not set SRID on the new geometries. I need a tool that works on server side and I would not to extract data from the source database, convert them in some intermediate format, to apply the srs conversion and then put the data in the new destination database. There is any other open source ETL or some other tool? Thank you very much in advance Cesare


KML file saved from KMZ missing tabular data

I saved a KML file in Google Earth from this KMZ ... https://www.rampp- team.com/documents/pennsylvania/watershed/LowerSusquehanna/lancaster_pa/Lancaster_CSLF.zip . Upon opening the KML file in QGIS, all of the tabular data is the same for all of the polygons, even though when the KMZ was loaded into Google Earth, there was color coding based on, apparently, tabular data. Is there any way to preserve the data so that I can work with it in ArcView (shapefiles)? There seems to be missing information in the KML that is present when displayed in Google Earth.


ArcGIS - ModelBuilider - Collect Values not working as expected

I am trying to iterate through a folder of rasters, reclassify each raster, then combine all of the reclassified rasters using Cell Statistics (to get the overall SUM of the reclassified rasters). I have read countless help forums and my Model should be working, but I am not able to connect the "Cell Statistics" block to the "Output Values" from the "Collect Values." For some reason, the only variables available to be under the "Cell Statistics" input picker box are "Raster" and "Reclassified Rasters" - I would expect "Output Values" to be listed as well, as that is the link I am trying to make. Any help would be greatly appreciated!! ![ModelBuilder](http://i.stack.imgur.com/07BZl.png) ![CellStatsPicker](http://i.stack.imgur.com/wW0vb.png)


How do you symbolize "Unique Values, Many Fields" with more than 3 fields?

I have data that has 5 fields with unique values that I need to symbolize at one time. Is there a work around to allow me to do this or a python code I can run to use more than three fields?


How to fix "Connection failed - Check settings and try again" when connecting to PostGIS?

I have Postgres running via Postgres.app on rails already and am trying to connect QGis to the database. PostGIS is working already. (questions already found online deal with this issue only) I do not have a password or username set up for this database, as indicated in my database.yml file for rails:               development:       adapter: postgis       encoding: unicode       postgis_extension: true       schema_search_path: public,postgis     ffsdfsdfsdf  pool: 5       database: gisdatabase       host: localhost      So I've gone into QGis and am trying to set up a new PostGIS connection:               Service:     Host: localhost     Port: 5432     Database: gisdatabase     SSLMode: allow     username:      password:      I've changed the port around and managed to get this error:               "could not connect to server: Connection refused..."      Changing the port back to 5432, I get a different error, thus indicating that the connection is not being refused:               "Connection failed - Check settings and try again."      I'm thinking that perhaps its the blank "Service" form or perhaps an error with the empty username and password fields that are causing the error, but I don't really know. I've tried to find the correct Service name but have been unsuccessful


Is it possible to offset lines with the same symbology?

I'm working in ArcMap 10 and have been tossed an electrical system map drawn in CAD to recreate in ArcMap. In several sections, different line numbers share the same structures for long stretches. Many times these have the same voltage (which is what they are symbolized on). Is there a way for me to continue symbolizing these lines by voltage while showing 3 separate lines that are not stacked for one section, 2 separate for another, etc. etc? Thanks for your help!


How to convert a line to points?

I'm working with QGIS and I want to convert a line into many points, does someone know how? Thanks.


What is the single most influential book every GIS professional should read?

If you could go back in time and tell yourself to read a specific book/text/guide/tutorial at the beginning of your career as a gis-developer, which book would it be? I expect this list to be varied and to cover a wide range of things. * * * This same question was once asked at stackoverflow for programming topics, and it has been considered by some as "one of the most important question ever asked on stackoverflow". So I decided to bring this question to serve our needs here at gis.stackexchange.


Reconcile Versioned Geodatabase with only DEFAULT version?

Follow up to this: Spatial Database Connection String I am trying to reconcile a versioned geodatabase that has only one version, DEFAULT. I am following the steps from the ESRI documentation: A quick tour of reconciling a version. It seems like there is no reason to have to do this, but the "OLE DB" table and "spatial" table are out of sync until I do this. As per comment from @mwalker in linked question above: > If you're using versioned SDE then you may not see changes right away, since > versions are handled with delta tables and views. The source tables won't > get updated until the version is reconciled.


Select vectors overlapping a raster

I have a raster layer with travel times from a specified point in my city's downtown. What I'd like to do is select all of the census blocks that lie at least partially within 45 minutes of downtown (i.e. I'd like to select all of the vector features that contain at least one raster pixel with a value of less than 45). I'd like to do this in QGIS if at all possible; suggestions for other FOSS software are also welcome.


Nearest Neighbor calculation in PostGIS

I'm learning how to use PostGIS and spatial databases for analysis. What I am trying to do is perform a calculation to get the distance for the nearest polygon in a file, using the edge calculation, rather than vertices. Using this answer from Paul Ramsey to a similar question: > CREATE TABLE mytable_distances AS SELECT a.id, b.id, > ST_Distance(a.geom::geography, b.geom::geography) as distance FROM mytable > a, mytable b; I am attempting to apply it to my spatial database. I don't understand the structure of this query though. I think `CREATE TABLE mytable_distances AS` creates a table to store the result but after this part I'm lost. Are `a` and `b` column names? If so, why would I specify two columns to calculate this? My table is called `TestArea` and I have experimented with some basic queries successfully:               SELECT        "TestArea".hgt     FROM        public."TestArea"     WHERE       "TestArea".area > 100      The structure of the database in PGAdmin III is as follows, with my table called `TestArea`. I'm not sure what the nearest neighbor calculation should look like using my column headers (all of these objects are polygons). ![enter image description here](http://i.stack.imgur.com/YXhdk.png)


Accurate and fast nearest neighbour, which distance function should I use and how to it speed up?

I am solving the following problem: I have a map with many points in it and I need to find the nearest point to given location. Each point carries information about it's latitude and longitude. This is for the routing algorithm to find an origin and destination point closest to the user selected locations. I wanted to use the "<->" operator for finding the nearest neighbour. But when testing I encountered the difficulty. It seems that the <-> is not always working properly for ordering points by distance. When running the following query, for the first two SELECT pairs it seems that point p1 is closer to location L. While from other two query pairs the result shows that in fact p2 is closer to location L. I checked that p1 is closer to p2 using http://www.movable-type.co.uk/scripts/latlong.html. This site return the distance in km similar to the result of the third SELECT pair, which is correct, they are both using spherical model of the Earth. By p1 I mean 14.452601,50.076207. By p2 I mean 14.452789,50.075481. By L I mean 14.451722068277277,50.07564315647511 (longitude,latitude). Is this correct, or am i making some mistake?               SELECT     (SELECT ST_Distance(st_setsrid(st_makepoint(14.452601,50.076207),4326),st_setsrid(st_makepoint(14.451722068277277,50.07564315647511),4326)))AS geometrical_from_p1,     (SELECT ST_Distance(st_setsrid(st_makepoint(14.452789,50.075481),4326),st_setsrid(st_makepoint(14.451722068277277,50.07564315647511),4326)))AS geometrical_from_p2,          (SELECT (st_setsrid(st_makepoint(14.452601,50.076207),4326) <-> st_setsrid(st_makepoint(14.451722068277277,50.07564315647511),4326)))AS geometrical_from_p1,     (SELECT (st_setsrid(st_makepoint(14.452789,50.075481),4326) <-> st_setsrid(st_makepoint(14.451722068277277,50.07564315647511),4326)))AS geometrical_from_p2,          (SELECT ST_Distance_sphere(st_setsrid(st_makepoint(14.452601,50.076207),4326),st_setsrid(st_makepoint(14.451722068277277,50.07564315647511),4326)))AS geomethrical_from_p1,     (SELECT ST_Distance_sphere(st_setsrid(st_makepoint(14.452789,50.075481),4326),st_setsrid(st_makepoint(14.451722068277277,50.07564315647511),4326)))AS geometrical_from_p2,          (SELECT ST_Distance(ST_GeographyFromText('POINT(14.452601 50.076207)'),ST_GeographyFromText('POINT(14.451722068277277 50.07564315647511)')))AS geographical_from_p1,     (SELECT ST_Distance(ST_GeographyFromText('POINT(14.452789 50.075481)'),ST_GeographyFromText('POINT(14.451722068277277 50.07564315647511)')))geographical_from_p2     ;      the 8 results for the 4 query pairs are: 0.00104424158783511; 0.00107918396178276; 0.00104493638333828; 0.00108122276233257; 88.6845186250299; 78.2447351379485; 88.8360501337395; 78.4754154775416 My application is written in java with connection to PostgreSQL 9.2 and PostGIS 2.0.3 . Now I am using simple solution written in java which is fast and accurate enough: I just go through all the points P loaded in RAM, calculate their distance from location L and remember the one with minimal distance. I use following formula:               simple_distance = SQRT[((L.longitude-P.longitude)*LON_DEG_M)^2 + ((L.latitude-P.latitude)*LAT_DEG_M)^2]      LON_DEG_M and LAT_DEG_M represent length(in meters) of one degree of longitude/latitude at a latitude of L. I calculate them only once for L and use these values for all points P. This does not give exact distance but I use it just to find the nearest point. (SQRT is not necessary.) On my laptop it takes about 2 milliseconds to find the nearest node this way. In the future I would like to PostGIS to calculate the nearest neighbour(s). On the internet I found the soulution suggesting using <-> operator and I tried ST_distance(geometry,geometry) function as well. As I described earlier I did not manage to make these work properly. Now I have implemented versions using ST_distance_sphere(geometry,geometry) and ST_distance(geography,geography) which seem to work properly but are slower for me. I am using following queries:               SELECT osmid FROM nearest_node_geom_sphere ORDER BY ST_Distance_sphere(geometry,st_setsrid(st_makepoint(L.longitude,L.latitude),4326)) LIMIT 1;          SELECT osmid FROM nearest_node ORDER BY ST_Distance(geography,ST_GeographyFromText('POINT(L.longitude L.latitude)')) LIMIT 1;      Is there a way to use the <-> operator properly to return the order the distances. If no how can I improve the performance of ST_distance_sphere(geometry,geometry) and ST_distance(geography,geography). Which should be fastest to use? For the current number of nodes I am getting the following times: (2 milliseconds for the simple java algorithm described earlier, 50 ms for St_distance_sphere(geometry,geometry), 90 ms for ST_distance(geography,geography)). I have around 9000 points in db running locally.


How to calculate area in square kilometers of a OSM boundary

I have tried using measurement plugin in JSON. It does not work if the boundary is made of multiple selected lines. That plugin works when it is a single line for the whole boundary. Is there any way to do this?


Where did my points go?

I had an excel file marking points on my map. I then added a base map and my points are gone. Even with the base map layer turned off the points do not show. I tried reordering the layers and tried deleting then re-adding the excel file and nothing seems to work. Any ideas?


Las 1.2 format displays extremely long thin line of points

I am having difficulty properly visualizing aerial LIDAR data from the following source: http://csc.noaa.gov/htdata/lidar1_z/geoid12a/data/1124/ Taking any of the laz files, the point cloud displays as an extremely long thin line of points. Zoomed in, there are cyclic perturbations in the points along that line, which leads me to believe that some sort of transform is necessary to properly visualize the data. For comparison sake, the aerial LIDAR data set found here visualizes as one would expect: http://opentopo.sdsc.edu/gridsphere/gridsphere?gs_action=datasetMetadata&cid=geonlidarframeportlet&otCollectionID=OT.092011.2875.1 Both claim to use NAD83 as the horizontal coordinate system. From what I can tell, the only major difference is one is LAS 1.0 and one is LAS 1.2. That being said, I don't really know what I'm doing. Any help in explaining the difference between the two data sets is greatly appreciated. A little bit of background to help guide an answer. I do not come from a GIS background, but rather am interested in comparing my LIDAR point clouds against a geographically accurate model, say a quality aerial LIDAR collect. For now, while learning, I have been using the free applications LAStools and CloudCompare.


Touch up tool for rasters?

I was wondering if it were at all possible to manipulate/massage rasters (eg Geotiff) using Qgis plugins or Grass?? By manipulate, I mean manually smooth off interpolation artefacts, blended surface disagreement, etc. Like a photoshop 'touch up tool' if you will. I've had a crack at trying with GIMP, but I lose all my geotiff info. If anyone hasn't already developed it, they should. Update: GIMP Works for 8bit Geotiff files, however, all data below 0m and above 255m are lost. Need to find a software package that deals with float64 Geotiffs.


My class using ogr to calculate distance is returning incorrect results away from the equator

In an answer to this question ustroetz kindly provided some sample code for use of ogr to do lat/lon distance and point in polygon calculations. I've created a Python class based on his example to do distance calcs but the results aren't correct. I realize that I could simply use haversine but the class is meant to grow to include other ogr functions, this is just the first step. Allowing for the small difference due to the earth not being a perfect sphere, I'd expect the distance between degrees of latitude to be approximately 60 nautical miles anywhere from the equator to the poles, but the code in my class is returning approx 60 nm at the equator and 89 nm at 47 degrees latitude. I'm geo-ignorant so don't know how to go about fixing this. I assume that the problem is in the class **init** where the coordinate transform is set up. Perhaps EPSG 3857 is not the correct projection to convert to from EPSG 4326? Since distance is not calculating correctly I think that a point in polygon method (to be added to the class) will also be inaccurate unless I perform a repair.               #!/usr/bin/python          from __future__ import print_function     from osgeo import ogr     import osr     import sys          class Geo(object):         def __init__(self):             """             Set up coordinate transform             """             inSpatialRef = osr.SpatialReference()             inSpatialRef.ImportFromEPSG(4326)           # WGS 84             outSpatialRef = osr.SpatialReference()             outSpatialRef.ImportFromEPSG(3857)          #Spherical mercator             self.coordTransform = osr.CoordinateTransformation(inSpatialRef, outSpatialRef)              def distance(self, position_1, position_2, units="km"):             """             Returns distance in 'units' (default km) between two Lat Lon point sets                     position_1 and position_2 are dicts containing "lon" and "lat"                          Units (optional) can be "km", "sm", or "nm"             """             point1 = ogr.Geometry(ogr.wkbPoint)             point1.AddPoint(position_1["lon"], position_1["lat"])             point1.Transform(self.coordTransform)                  point2 = ogr.Geometry(ogr.wkbPoint)             point2.AddPoint(position_1["lon"], position_2["lat"])             point2.Transform(self.coordTransform)                  raw_dist = point2.Distance(point1)                  print("Raw distance from ({0:6.2f}, {1:6.2f}) to ({2:6.2f}, {3:6.2f}) is {4:14.6f} km, "                 .format(position_1["lon"], position_1["lat"],                         position_2["lon"], position_2["lat"],                         raw_dist), end="")                  if units.lower() == "km":                 return raw_dist * 0.001             elif units.lower() == "sm":                 #return raw_dist * 0.00062137                 return (raw_dist / 1000) * 0.621371192             elif units.lower() == "nm":                 #return raw_dist * 0.000539956803                 return (raw_dist / 1000) * 0.539956803          if __name__ == "__main__":             g = Geo()             print("{0} {1}".format(g.distance({"lon":147.0, "lat":0}, {"lon":147.0, "lat":1}, "nm"), "nm"))             print("{0} {1}".format(g.distance({"lon":147.0, "lat":47}, {"lon":147.0, "lat":48}, "nm"), "nm"))      The output produced by the code below is:               Raw distance from (147.00,   0.00) to (147.00,   1.00) is  111325.142866 M, 60.1107682357 nm     Raw distance from (147.00,  47.00) to (147.00,  48.00) is  164780.762454 M, 88.9744936905 nm


Why might Ordnance Survey codes all be out by 500m?

I've been given a set of UK Ordnance Survey geocodes, but I think they're all off by about 500 metres. Here is an example - please could someone confirm to me that the codes are wrong, not my calculations? The OS geocode that I've been given is `SO7567`, for a place called Abberley, in Worcestershire. My script converts this to `52.30136731, -2.36798287` here: http://goo.gl/maps/9a3Yu However, Abberley is here, about 500 metres north: http://goo.gl/maps/OAD9L This conversion site suggests my calculations are correct, and it is indeed the original coders (who looked at Abberley, and decided that its Ordnance Survey code was `SO7567`) who were out of whack. Could anyone suggest why they might all be out by exactly this amount - could it be a projection issue? I'm not an geographer, and the person who gave them to me swears they should be right! Thank you for nay help.


postgis error: when adding features geometry has z dimension but column does not

I have to define a geometry column in PostGIS for 2D points, but when I try to import points from DXF to postgres I get the following message: > "postgis error: when adding features geometry has z dimension but column > does not" Is there any way to change DXF files to be 2D geometry? I have tried this with the comand `Flatten` in autocad but the same error appears.


Help with opentripplanner

I have started working with opentripplanner (i am developing a journey planner app). I am following the two & five minute tutorial given on opentripplanner website. but while starting server using "bin\start-server.bat" giving me errors which are shown in picture, i also tried to google for error messages but didn't find any fair solution. kindly help me in this matter.![C:\\apache- tomcat-7.0.33\\otp\\bin>start-server](http://i.stack.imgur.com/fefw5.png)


OpenLayers - Get features inside a polygon

In my OpenLayers map I allow users to draw an ad-hoc polygon vector. Then I need to get all the features from a WFS vector layer that fall within the ad- hoc polygon. Here is the code I have that allows the user to draw the ad-hoc polygon in the map. This code works fine:               function buildIt() {     //CREATE THE NEW VECTOR      var polygonLayer = new OpenLayers.Layer.Vector("Poly Layer");     //ADD IT TO THE MAP     map.addLayer(polygonLayer);          //CREATE A DRAWFEATURE CONTROL     var polyedit = new OpenLayers.Control.DrawFeature(     polygonLayer,//POINT IT AT THE USER DEFINED POLYGON LAYER     OpenLayers.Handler.Polygon     );          // 'point' is one of the named callbacks specified by OpenLayers.Handler.Polygon     // and is called every time a point is added to the feature being drawn     polyedit.handler.callbacks.point = function (pt) {          console.log(pt)     }     map.addControl(polyedit);//ADD THE DRAWFEATURE CONTROL TO THE MAP     polyedit.activate();//ACTIVATE THE DRAWFEATURE CONTROL      } After this I am clueless. How can I get all the features that fall within the user defined polygon? Specifically, I need to get the attributes from each WFS feature (sometimes points, sometimes polygons) that fall within the user drawn ad-hoc polygon. Any suggestions would be appreciated! Thanks.


QGIS UTM zones - TIN interpolation plugin

I am trying to add a Z-variable to a vector in order to eventually calculate the volume or the zone using TIN. I was told I first need to set my UTM in order for my Z-variable to be recognized in meters. If my vector spans accross various UTM zones, which ones do I use?


Inland / crater watershed delination

Even though watershed delineation is a quite popular matter here, I have what seems to be a somewhat unusual demand. I start with a DEM which is constituted of a series of craters, all centered on one point which is the lowest part of the crater. Here is a picture of my DEM What I want to to, is to obtain the drainage basin of each crater. I know how to delineate watersheds using ArcGIS' Hydrology tools, but they are best suited for 'open' hydrological systems (flows try to reach the exterior of the raster). Here is the result I obtain. How comes some of my drainage basins are obviously wrong? Why is Arcgis so lazy that its scope is limited to the rectangle limited by the four most extreme (N-S-W-E) points? Therefore, my question is the following: is it possible to automatically delineate the watershed basin and set the center of my craters as the targets where my flows should go? Thanks in advance for your help! Damien [Edit : added Arcgis output for clarity]


GRASS i.landsat.toar returns empty raster

I am trying to use the GRASS module i.landsat.toar. Currently the files created are empty and it spits out: ![enter image description here](http://i.stack.imgur.com/SFfDN.png) So it seems that it is not picking up something? The location is set to the extent of the landsat image. Does anyone have an idea where I am going wrong? Becky


Making Use Of Railway Route from Google Maps API

I have a specific requirement where in i have to display or highlight only railway stations in the map rather than the whole map. How can this be achieved.?. Please find the below code that i have tried.

                    Any help will be much appreciated.. Thanks in advance..


Ogr: How to add layer description to dgn?

I have created a dgn-file using Ogr, but can't find a way to fill out the description for this layer programatically. (something like _dgnLayer.Description = "Highways")               opts.Add("SEED=" + seedfile);     opts.Add("COPY_WHOLE_SEED_FILE=YES");          // Prepare     Initialization.PrepareGDAL();          // Initialize output DGN     using (var dgnDriver = Ogr.GetDriverByName("DGN"))     {         _dgnDataSource = dgnDriver.CreateDataSource(filename, opts.ToArray());         _dgnLayer = _dgnDataSource.CreateLayer("", null, wkbGeometryType.wkbUnknown, null);         _dgnFeatureDefinition = _dgnLayer.GetLayerDefn();         _dgnFeatureDefinition.AddFieldDefn(new FieldDefn("ErrorTest", FieldType.OFTString));     }      Regards Lars


calibrating vector line in qgis

I have loaded GPS points into a QGIS vector layer. I also imported a vector line from a dxf file which doesn't have georeferencing. I would like to fit in line geometry between the GPS points - how can I do this?


How to extract numbers from each vectors?

I'm trying to create custom legend for a layer where OpenLayers.Rule objects are created dynamically and added to a layer.               new OpenLayers.Rule({                     symbolizer: {                         strokeColor : color,                         fillColor: color                     },                     filter: new OpenLayers.Filter.Comparison({                         type: OpenLayers.Filter.Comparison.BETWEEN,                         property: p,                         lowerBoundary: from,                         upperBoundary: to                     })                 })); and so on ..      I need to extend functionality of legend where for each rule I need to extract how many vectors affected by the rule so that when user changes interval boundaries as above, drawn legend replaced by new numbers. Is there any built-in functionality for this in OpenLayers ?


Why is GPS a free service?

The Global Positioning System ( **GPS** ) **is a free service** that is owned and operated by the U.S. Government and is always available. When we buy a GPS device, we don't pay a monthly fee or pay a tax for GPS support. We only pay the price of the device. **I don't understand how it's possible**. How does the U.S. make money?


How to convert factfinder2.census.gov map data to KML?

I've downloaded a map from http://factfinder2.census.gov/faces/tableservices/jsf/pages/productview.xhtml?pid=DEC_10_SF1_P1 (limited to San Francisco census tracts), and I'd like to convert it to KML format. However, after unzipping the shapefile .zip file and running               $ ogr2ogr -f KML pop_by_census_tract.kml 140_00.shp      I get the following error message:               ERROR 6: No translation for Mercator_Auxiliary_Sphere to PROJ.4 format is known.     Warning 1: Failed to create coordinate transformation between the     input coordinate system and WGS84.  This may be because they     are not transformable, or because projection services     (PROJ.4 DLL/.so) could not be loaded.     KML geometries may not render correctly.     This message will not be issued any more.           Source:     PROJCS["WGS_1984_Web_Mercator_Auxiliary_Sphere",         GEOGCS["GCS_WGS_1984",             DATUM["WGS_1984",                 SPHEROID["WGS_1984",6378137.0,298.257223563]],             PRIMEM["Greenwich",0.0],             UNIT["Degree",0.0174532925199433]],         PROJECTION["Mercator_Auxiliary_Sphere"],         PARAMETER["False_Easting",0.0],         PARAMETER["False_Northing",0.0],         PARAMETER["Central_Meridian",0.0],         PARAMETER["Standard_Parallel_1",0.0],         PARAMETER["Auxiliary_Sphere_Type",0.0],         UNIT["Meter",1.0]]          ERROR 1: Latitude 4553026.920500 is invalid. Valid range is [-90,90]. This warning will not be issued any more     Warning 1: Longitude -13627851.842700 has been modified to fit into range [-180,180]. This warning will not be issued any more      When I open the resulting KML file in Google Earth, I get randomly scattered lines across the entire planet. The example at code.google.com/apis/kml/articles/vector.html works fine, so the installation of GDAL (1.8.0 from macports) isn't totally broken. http://trac.osgeo.org/gdal/ticket/3962 seems related, but I'm not sure if it's saying that Mercator_Auxiliary_Sphere is just not supported by GDAL, or if it's implying I should pass another command line option to make it work. Edit: I wound up accomplishing my goal by downloading http://www2.census.gov/geo/tiger/GENZ2010/gz_2010_06_140_00_500k.zip, which contains shapefiles that ogr2ogr and shp2pgsql can understand, and joining it with the pure table data from factfinder2. I'm still curious about the format factfinder2 is using, but it's no longer blocking progress for me.


Examples of Google styled Maps? (including JSON)

After watching a great video on creating different map styles with Google Maps, I have been having some fun playing around with the Google Maps Style Wizard. I was hoping to find a gallery of user submitted Google Maps Styles, with accompanying JSON code so that I could pinch their hard work & cartography skills, and use on my own applications. Has anyone come across either a collection of styles (+JSON) or at least some great examples of google maps styles, where the JSON is also available? The closest I could find was this.


Efficiently querying a region-partitioned PostGIS table

I have two PostGIS tables: one of `streets`, and one of `regions`, which, apart from spanning the same spatial extent, are a priori unrelated. The `streets` queries I perform can be restricted to a given region, by supplying a known id. To exploit this, I added a `region_id` column to the `streets` table (using a `st_intersects`-based update query, which was quite expensive, but performed only once). My question: once this is in place, will the use of an index on the `streets.region_id` column be enough to make queries of this type maximally efficient:               select * from streets where region_id = <id> and ...      Or is there an alternative mechanism that would perform better?


A formula error in raster calculator with arcgis 10.1

help me to solve a problem at a formula to calculate the susceptibility to landslide that gives me errors when computing with "raster calculator" in arcgis 10.1, the formula is exact but errors in the expressions: "(", ")" or "," and here is the formula: Con("aspectU0" < 90,0,"aspectU0" > 90 & "aspectU0" < 180 , ("aspectU0" - 90) / 90,"aspectU0" > 180 & "aspectU0" < 270,(270 - "aspectU0") / 90,0,"aspectU0" > 270 & "aspectU0" < 360,0) or help me to convert the following formula which is executable in arcgis 9.3 to be executable in arcgis 10.1: con([aspectdem] le 90, 0, [aspectdem] gt 90 & [aspectdem] le 180, ([aspectdem] - 90) / 90, [aspectdem] gt 180 & [aspectdem] le 270, (270 - [aspectdem]) / 90, 0, [aspectdemf] gt 270 & [aspectdemf] le 360, 0) Thank you in advance


Where can I obtain 3D data for Singapore?

I am looking for providers of 3D maps/models of Singapore. This might be either from photogrammetry/laser techniques or based on a DEM + building footprints/landuse with a height attribute. Is any such data available either in public repositories or from specialized companies?


OpenLayers 2.13.1 - Google layers

I'm following the examples of a free on OpenLayers, in the code below I am calling various layers of google maps api. Have checked with the code book and is correct but nothing appears on the broswer Tested with one at a time and appeared, does anyone know what is wrong in the code? Also I've taken the spec version of the api call to google maps and did not work. Nothing appears.                   <!DOCTYPE html>     <html lang='en'>     <head>      <meta charset='utf-8' />      <title>My OpenLayers Map</title>      <script type='text/javascript' src='OpenLayers.js'></script>           <!-- Chamada a Api Google Maps especificando a versão 3.2 -->      <script src="http://maps.google.com/maps/api/js?sensor=false&v=3.2"></script>           <script type='text/javascript'>           //Variável map      var map;           //Função a ser chamada para criar o mapa na leitura da página      function init() {           // A classe Map espera dois argumentos "map_element" que o ID do elemento HTLML geralmente uma DIV e opções de mapa {key:value}      map = new OpenLayers.Map('map_element', {});            //Função de inicialização  - google.maps.MapTypeId.HYBRID (híbrida)      var google_hybrid = new OpenLayers.Layer.Google(         "Google Hybrid",         {type: google.maps.MapTypeId.HYBRID}         );           //Função de inicialização  - google.maps.MapTypeId.TERRAIN (topológica)      var google_physical = new OpenLayers.Layer.Google(         "Google Physical",         {type: google.maps.MapTypeId.TERRAIN}         );           //Função de inicialização  - google.maps.MapTypeId.TERRAIN (topológica)      var google_satellite = new OpenLayers.Layer.Google(         "Google Satellite",         {type: google.maps.MapTypeId.SATELLITE}         );           //Função de inicialização  - google.maps.MapTypeId.ROADMAP (street layer default)      var google_streets = new OpenLayers.Layer.Google(         "Google Streets",         {             );          //Adiciona as Google Layers ao mapa     map.addLayers([google_hybrid,google_physical,google_satellite,         google_streets]);                //Camada de controle que vai mostrar as camadas no mapa      map.addControl(new OpenLayers.Control.LayerSwitcher());           //Verifica se o mapa tem um ponto central e o extende a sua extensão máxima      if(!map.getCenter()){      map.zoomToMaxExtent();           }      }           </script>     </head>          <body onload='init();'> <!-- Chama a função js init() -->           <!--Elemento HTML onde o mapa é exibido -->      <div id='map_element' style='width: 1000px; height: 900px;'>      </div>          </body>     </html>


Where does MapInfo store its error and log files?

I am having errors opening a table with spatial co-ordinates in Microsoft SQL and my question in gis.stackexchange is here. I want to know where are the error and log files of Mapinfo stored. I searched over the Mapinfo installation but was unable to find them.


Creating Latitude Grid from DEM

I am trying to evaluate an equation in the raster calculator and I need latitude values. I have the DEM for the region and am wondering if it is possible to create a latitude raster based on the DEM? Thanks a lot!


How to reproject spatial data using free libraries?

How can I use free libraries to transform spatial data? For example, I want to change the projection of a Shapefile within the code of my C# web application. How do I do that?


How to solve Projection Problem at client side

I want to discuss a scenario here to all.   1. I have developed a Web mapping application in OpenLayers and GeoServer. The customer give us shapefiles with native projection SWREFF99_1330 (EPSG:3008). i published this shapefile on Geoserver after changing its native projection from EPSG:3008 to EPSG:4326 by using ogr2ogr tool. The layer displayed perfectly over Google map layer in my OpenLayers Web Page.    2. Again, I publish this shapefile on GeoServer with its native projection (this time, without changing projection using ogr2ogr tool) and then i go to my page and figure out that the layer was not visible over Google map layer.  Now, my question to all members: Is there any way to change the native projection from EPSG:3008 to EPSG:4326 at client side (in my OpenLayers web page code) ? so that the layer should display correctly over google base layer as in case 1.


How to select nth row in mapbasic?

There are multiple rows in my .tab file in mapinfo. Is there any way to select nth row in mapbasic window of mapinfo?


Looking for map calibration software and city database

I'm in need of a mapping software which, given a map, lets you calibrate it (choose which geographic coordinates project to which pixels on the map image) and then given a database of cities, you can navigate to the specified city on the map you just calibrated. To reword the question, I need a mapping software that can calibrate my map image and accepts some kind of database format to display cities on that map. Where to find such compatible map software and database? Preferably free and open source software and data. If the question is too broad, I appreciate any points making me near to the answer.


Labels not appearing on any layer in ArcMap?

Not sure what the issue is, but I'm unable to get labels to display in ArcMap 10. The obvious things are accounted for (as far as I know). Labels are set in Properties, layers are checked to label in the TOC. I've even restarted to see if that would work. **UPDATE: 3/3/2011** After trying all of the mentioned solutions, I was still unable to get labels to display. It wasn't specific to a single layer, but rather to any layer in the .mxd. I opened a new document, added the layer I'd been trying to label, and was able to display labels just fine. What may have caused the initial document to globally make labels not display?


SLD for layers present in a table in DB

My layers are present in postgis db which are being served to openlayers JS through geoserver. I have my layers corresponding to polygons,lines etc. whose attributes are stored in the same table itself. Like Halls table contain columns (hid, hall_geom,hall_name) etc. I want hall_name of each hall_geom(polygon) to be placed as a label on it in the map. I have tried SLD- styling form geoserver but couldn't make it successful(I think the styles given are for shape files and not tables in db). Please show me a right way