title: How to smooth / "reinterpolate" a raster with GDAL?

text: Short:   I want to change the resolution of a raster and smooth the gray ramp like shown in the images bellow. The preference is to use GDAL, PIL or Numpy. Description:   I'm kriging points into rasters with an output resolution of 20 meters with the High Performance Geostatistical Library. I don't want to change the output resolution because the interpolation time increase exponentially.   With this resolution the output image is ugly (pixelated and aliased). I don't know if it is conceptually correct but I want the image to be smoother like in the example bellow. It's something like 'reinterpolating' the image into a better resolution one. I'm using python so my preferences are GDAL, Python Imaging Library or Numpy. The answer could be theoretical, like pointing out the algorithm name or the concept of this kind of operation. Source:   ![enter image description here](http://i.stack.imgur.com/PP0wg.png) Destination:   ![enter image description here](http://i.stack.imgur.com/T5g3w.png) EDIT Results with gdalwarp cubic spline:   ![enter image description here](http://i.stack.imgur.com/S7ZDn.png)


title: How can I properly add PostGIS data into MapGuide Maestro?

text: I am having trouble connecting a PostGIS database to MapGuide Maestro. Actually, I have added the PostGIS data into Maestro, but when I click to validate the layer it shows an error message: `Validating resource: Session:9754d026-d86b-11e1-8000-00270e208411_en_7F0000010AFC0AFB0AFA//259f852d-fc68-4475-9de1-a5e446fb6e82.FeatureSource Warning - OSGeo.MapGuide.ObjectModels.FeatureSource_1_0_0.FeatureSourceType: Spatial context extent appears to be invalid (or is the default extent returned by FDO provider. Meaning it very likely does not truly represent the extent of your data) Information - OSGeo.MapGuide.ObjectModels.FeatureSource_1_0_0.FeatureSourceType: No primary key defined for class: public:geography_columns, features will not be selectable and feature joins on this class may produce unexpected results` Can anyone please suggest where I went wrong? **Is there any documentation for how to add PostGIS data into MapGuide Maestro?** Also, if I click to see the preview of the layer it gives me an error report of `Fusion Error`. I am using PostgreSQL 8.4, PostGIS 1.5, MapGuide 2.2 and MapGuide Maestro 4.01 Thanks in advance..


title: Overlaping graphics not visible in featureLayer

text: I have constructed featureCollection from Json and feature collection has 44 features. I have given this feature collection as input to the Featurelayer which is showing only 32 graphics objects of the featurecollection.I found out that overlapping features are not rendered at all in the feature layer. How could I resolve this issue and what is the possible solution to this problem.? Thank you in advance. Thank you in advance.


title: What are the equivalents of ESRI products in open source for Web Mapping?

text: I have been using the following licensed products to publish my GIS data on web.   1. web server: windows 2008 R2,   2. GIS desktop: ArcGIS desktop,   3. GIS web server: ArcGIS server,   4. API: ArcGIS API for JS,   5. Database: Oracle with ArcSDE,   6. IDE: Eclipse. But now I want to switch altogether to open source world mainly because of freedom. Fortunately there is plethora information about open source on web but unfortunately are confusing me a lot. There are soo many products like geoserver, mapserver, openlayer, geotool, vaadin, geomajas, uDig, GRASS etc etc. I am not able to related these products with the one I am already using. Can someone please relate the equivalents?


title: InfoWindow vs. Popup

text: Can somebody provide a good explanation/comparison between Info Window vs. Popup Window. It seems both can be used with Feature layers. What are the pros and cons of either and when is more appropriate to use one than the other?


title: How to setup GeoWebCache to produce same grid as GoogleMaps?

text: Is it possible to create gridSet in GeoWebCache which has same number of zoomlevels with scaleDenominators compatible with google maps? There is EPSG:4326 default, but no mention of number of zoomlevels/scaleDenominators in documentation.


title: Use of GIS in fixed telecommunication network management

text: I want to get an idea how GIS can be used in fixed line telecommunication management? there are open source GIS tools as well as Proprietary tools as well. what is the best tool can be used in above ? thanks. Chandima


title: Activity Dashboard for LocalHost

text: ESRI recently put out a new tool called the 'Activity Dashboard' that checks for new maps created, map uses, etc. Sadly, this tool only works for those maps created through the ArcGIS online system. One of the problems I have as the new 'Web Map' guy for my organization is that multiple maps have been created, all of them as (basically) one-offs to fix a specific problem, and no one has a master list of them or the services they are running on. This left me with a very annoying situation. I 'thought' I knew about them all, but I find that I am regularly learning about old maps as they break. The recent address locator that ESRI killed off brought SEVERAL maps to my attention that I wasn't even aware of. I am wondering if there are any tools that have been developed that track maps and/or services that have been created to give someone in my position a better chance of managing this system.


title: Gaussian Plume Model

text: I'm looking for code to implement a gaussian plume model in 3D. Any suggestions on where to start? Are there any python scripts out there? Thanks


title: applying SLD to WMS raster in OpenLayers

text: I am attempting to adjust the symbology of a loaded WMS raster layer (img file through MapServer). It is the NLCD dataset. The WMS layer loads find with               nlcd = new OpenLayers.Layer.WMS(         'nlcd',         'http://server /fswms/html/rlayers',         {             layers: 'nlcd-2006', transparent: true,             projection: new OpenLayers.Projection("EPSG:4326")         },         {             isBaseLayer: false         }     );      And I have a SLD file that I have created in which I am attempting to make all of the pixels with the value of 11 opacity of 0.               <?xml version="1.0" encoding="ISO-8859-1"?>     <StyledLayerDescriptor version="1.0.0"         xsi:schemaLocation="http://www.opengis.net/sld StyledLayerDescriptor.xsd"         xmlns="http://www.opengis.net/sld"         xmlns:ogc="http://www.opengis.net/ogc"         xmlns:xlink="http://www.w3.org/1999/xlink"         xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance">     <NamedLayer>         <Name>nlcd</Name>         <UserStyle>           <Title>ncld test</Title>           <FeatureTypeStyle>             <Rule>               <RasterSymbolizer>                 <ColorMap>                   <ColorMapEntry opacity="0" color="#000000" quantity="11" />                 </ColorMap>               </RasterSymbolizer>             </Rule>           </FeatureTypeStyle>         </UserStyle>       </NamedLayer>     </StyledLayerDescriptor>      The WML layer ends up being map.layers[3] in my test application so I attempting to apply the SLD by the following code but it is not working.               //this is in the init:        OpenLayers.Request.GET({         url: "sldcode.sld",         success: complete     });            //function that attempts to apply SLD     function complete(req) {         mySLD = format.read(req.responseXML || req.responseText);         map.layers[3].mergeNewParams({sld_body:mySLD});     }         Here a link to the dev site: http://rain.nemac.org/~derek/ol/sldTest.html Any help would be greatly appreciated. Thanks, Derek


title: Visualizing KMZ network link in QGIS

text: Is it possible to open and visualize and kmz network link in QGIS. If so what are the setps involved? Thanks in advance Cheers


title: GIS mapping software for visualizing and publishing maps

text: Currently I am writing my Bachelor thesis and for testing purposes I need software where I can visualize and publish maps and put vector data (point, line, polygon) on them. Can anybody recommend such software? I have already founded GIS Cloud. I need something like this.


title: Where can DEM data for all of British Columbia be found?

text: Is there a website or server from which a DEM covering all of BC can be downloaded - preferably at 1:50K scale? GeoBase has tiles for 1:50K and 1:250K, but they do not provide it as a mosaic. It's not very time efficient to download 1000+ tiles and mosaic them.


title: How to get an extent for Raster using python

text: I have to get an extent of a given raster so that I can save it and use the same to clip another raster with that extracted extent. As of now I am stuck in the extent part.               import arcpy     from arcpy import env               desc= arcpy.env.extent("D:\GIS @ UTD\Sem 1\GISC 6317\Lab\Lab 10\temp\Tahoe\Tahoe\Emer\erelev.grid")          print desc      The following gives an output 0 0 0 0 NaN NaN NaN NaN. Any suggestions.


title: How to improve cartographic visualization of overlapping polygons?

text: I need some advice about how to improve the visualization of cartographic information. User can select different species and the webmapping app shows its geographical distribution (polygonal degree cells), each specie with a range of color (e.g darker orange color where we find more info, lighter orange where less info). The problem is when more than one specie overlaps. What I am currently doing is just to calculate the additive color mix of two colors using http://www.xarg.org/project/jquery-color-plugin-xcolor/ As you can see in this image, the resulting color where two species overlap (mixed blue and yellow) is not intuitive at all. ![enter image description here](http://i.stack.imgur.com/2Zmhu.jpg) Someone has any idea or knows similar tools where to get inspiration? for creating the polygons I use d3.js, so if more complex SVG features have to be created I can give a try. Some ideas I had are... 1) The more data on a polygon, the thicker the border (or each part of the border with its corresponding color) 2) add a label at the center of polygon saying how many species overlap. 3) Divide polygon in different parts, each one with corresponding species color. thanks in advance, Pere


title: QGIS Advanced String Filter Expression returns too many results

text: `"GNR"` is a string column and I expect the filter expression `"GNR"='.11'` to only return the feature with the value `'0.11'`. Why are the other features with values of `'.110'`and `'.1100'` included? Is there something wrong with my query? ![string filter expression returns too many results](http://i.stack.imgur.com/enfUf.png)


title: Conversion of diagonally adjacent pixels into one polygon?

text: I am trying to vectorize a raster whose digital numbers indicate different classes they belong to. ArcGIS and QGIS yield similar results: separate shapes out of pixels sharing the class that are adjacent only diagonally. In this sense, this question is related with this one and also this one. I tried the suggestions posted in those questions, but none of them worked for me. I can not use the buffering-approach because my polygons refer to networks of areas; the connectivity/topology of those networks must remain intact, and I have nodes very close to each other. Dissolve is not useful either; it apparently does not work on polygons sharing a single vertex. ## ILWIS could solve the issue: I found that the vectorization tool of ILWIS (see explanation here) has an option for scanning for the 8 neighboring pixels during the targeting of each raster pixel to a polygon. ILWIS calls it "Output map 8-connected with smoothing". This might solve the issue by avoiding the polygon breakage during the vectorization when input pixels are placed diagonally. I downloaded and installed the ILWIS version available here. However, after importing the raster I want to vectorize (I used the built-in GDAL option), the drop-down menu to select the input layer in the 'Raster to Polygon' tool does not allow me to select anything and thus I can not use the tool. Does somebody know how to use this ILWIS' tool, or how to solve this problem?


title: QGIS snapping not working for me

text: I am new to QGIS, and this is a very basic question. I made a new polygon layer, to trace cadastral boundaries off a raster map. I want the polygon nodes on adjoining acreages to snap together. Some do because I have 'avoid intersections' on. But if I have to manually move a node on one polygon onto a node on the adjoining polygon, I cannot get it to snap to it. I've read the manual, tried to sort it out, but I cannot make it work. How do I make it work? Refer to the screenshot. ![Cadastral boundaries won't snap together](http://i.stack.imgur.com/wgoA2.jpg) I have used MapInfo for 14 years. It does not have an inbuilt way to avoid polygon overlaps, but its snapping is very straightforward.


title: Why I see different data from PostGIS and ArcSDE

text: I have installed ArcSDE 10 on PostgreSQL 8.4.3 with PostGIS 1.4.2. I did some tests.   1. I imported one feature layer with 8500 features into ArcSDE through ArcCatalog. and I deleted 8499 features in ArcMap and just kept one. I saved the edits. I checked the data using PostgresSQL client tool, like pgAdmin III. I still see 8500 rows in the table.   2. I imported one feature layer with 3 point features into ArcSDE through ArcCatalog. and I add 2 new points in ArcMap. I saved the edits. So I can see 5 points on this layer. But when checked the data pgAdmin III. I still only 3 rows in the table. I closed and opened ArcMap and pgAdmin several times. and still see they have different data. Is anything I did wrong?


title: Error when trying to import a raster file into a geodatabase

text: I am using ArcGIS 9.3.1. For some time I have had no problem importing a raster file into my geodatabase. That's recently changed and I have absolutely no idea why. Every time I try to import a raster file, I receive an error message that says: "Problem generating Output Raster Name (ERROR 999999: Error executing function.)" I am at a loss...any help would be greatly appreciated.


title: Openlayers error after removing selectFeature control

text: I am using a selectFeature control for a point layer in OpenLayers 2.12. I have a problem after removing the selectFeature control from the map. The reason why I want to remove the control is because the target layer of the control is destroyed and a new one is added every time the map is moved. In this process I also remove the selectFeature control from the map and add a new one for the updated layer. But that does not work properly. After the first move/zoom event of the map get the following error: _Error: getFeatureFromEvent called on layer with no renderer. This usually means you destroyed a layer, but not some handler which is associated with it._ Here is a version of the function `refresh`, which is called on any pan/zoom event of the map:               var mylayer;     var mycontrol;     refresh = function(){               // Remove previous layers and controls         if(mylayer != null){             mylayer.destroy();             map.removeControl(mycontrol);         }              // Get mylayer         mylayer = new OpenLayers.Layer.Vector("GeoJSON",              ...         );                    // Add layer and selectFeature control         map.addLayer(mylayer);            mycontrol = new OpenLayers.Control.SelectFeature(mylayer,{               onSelect: function(e) {             $("#tmpOutput1").html(e.attributes.name + "\n" + e.attributes.id);             }         });         map.addControl(mycontrol);           mycontrol.activate();             }      **UPDATE** To solve the problem I had to change the layer/control removal to this:               if(mylayer != null){         map.removeControl(mycontrol);         mycontrol.destroy();              map.removeLayer(mylayer);         mylayer.destroy();     }      Note that here the order matters (at least in my example): first remove the control, then remove the layer.


title: How to retrieve transport infrastructure data from google map?

text: I imported into postgis osm data using osm2pgsql. My objective was to get a database with roads elements for a specific region of Italy (i.e. traffic lights, driving directions and turning direction at junction). However there are a lot of missing data, because almost all the cells are empty. I would like to know if I did something wrong and if there is an alternative way to obtain these data. For instance I would like to know if there exists a way to retrieve data from google map (transport infrastructure or even real time traffic information)? Thank you Selini


title: How is a geographic heat map different from a colored statistical thematic map?

text: From **Wikipedia** , > A heat map is a graphical representation of data where the individual values > contained in a matrix are represented as colors. Apparently, a _geographic_ heat map is only one kind, and this is the only kind relevant to me. Cartography is about geographic maps. In cartography, I often distinguish between topographic and thematic maps. From **Wikipedia** , > A thematic map is a type of map or chart especially designed to show a > particular theme connected with a specific geographic area. In thematic cartography, we can at least distinguish between types of thematic maps: statistical, color, choropleth, isopleth, dasymetric, point-symbol, and so on. **How is a geographic heat map any different from a colored statistical thematic map?**


title: Multiple geometry columns for points, lines, and polygons in spatialite?

text: I'm rather new to the GIS world and even newer to Spatialite and trying to create a new DB for our application. I would like to create spatial index(es) on the geometry column(s), but our geometry BLOB's have always existed in a single column as WKT. After doing some research, it looks like you need a separate column for points, lines, polygons, multipoint, etc. I'm assuming this from Recipe #6 in the Spatialite Cookbook which states that you need to create your table first, then use AddGeometryColumn to actually add the column which will contain the geometry. In creating this column, one of the items you specific is the GeometryType. I would much rather keep a single geometry column that could handle all of the shape types, but if the only way to create spatial indexes is to have separate columns, I'll go that route. Any advice on this would be greatly appreciated. Editing this to add more info. Folks before me on this project created their own RTree tables and index to pull data that was within an MBR. My goal is to get rid of the home-grown stuff and use Spatialite for our queries. We have a number of shape files that contain points, lines, multilines, and polygons and I would rather insert their BLOB (WKT) into a single geometry column vs. 1 for each type. Our main queries are to find all shapes that are within an MBR. I hope this helps explain a bit better the end goal. Thanks.


title: Extend Polyline in ArcGIS without intersecting line

text: I am trying to extend the ends of each line I have in a seismic survey 2 miles. I have tried using the extend line (Editor) tool but it needs an intersecting line to move to. Does anyone know how to extend a previously created line a specified distance on both ends without intersecting another line?


title: How to create Race Route with arrow directions and distance markers using ArcGIS for Desktop?

text: Anyone have any suggestions regarding how I would approach creating a race route with arrow directions and distance markers in ArcGIS? I could do it manually but I thought I recalled seeing something about a plug in that aided in this type of work.


title: Portal Sharing API - "created" search parameter

text: I was wondering if anyone had an example of a working REST string searching using the "created" parameter. I have read **the docs** three or four times and every time I try to use it, it returns 0 results. Even though I know there are items in there that were created in the time frame I am providing. Example:               http://<srv address>/sharing/search?f=pjson&q=created:[0000001104537600000%20TO%200000001388534400000]      That is searching everything from January 1, 2005 January 1, 2014. The output from the above:               {       "query": "created:[0000001104537600000 TO 0000001388534400000]",       "total": 0,       "start": 1,       "num": 10,       "nextStart": -1,       "results": []     }      What am I doing wrong?


title: How to make lines styles with offset join properly?

text: I want to draw the lanes of a two-lane street in different colors. I tried to apply an offset to a symbol layer which works fine for a single line but the ends of two connected lines do not join properly. Depending on the angle, the lines overlap or there is a gap between the lines. Is there any way to solve my use case?


title: PostGIS - improving sql statement

text: I've managed to create the below statement that creates a table based on a point in polygon query. However it creates a column that gives a true/false value which I then delete out all the f values from. Is there a better way of producing this statement?               CREATE TABLE crime_test AS      SELECT a.*, b.town, ST_within(a.the_geom, b.the_geom) "within"     FROM crime_data a, nh_areas b;     DELETE FROM crime_test WHERE within = 'f'


title: Qgis - Labels automatically resize to fit within the polygon

text: The question of my problem is simple, I want to make labels automatically reduce to fit within the polignos. Type the biggest size is 12 and goat decrease the size to be restricted to the limits of Poligno. Exists based on the calculator using polignos option field area, but this is time consuming due to the amount of resources. Can I use plugins labels manualemnte easy, but that would take too long and would be basically a CAD but without assistance. ![enter image description here](http://i.stack.imgur.com/Xroab.png)


title: encoding problem with QGIS 1.8 and CSV import

text: I have following problem with QGIS 1.8. Lisboa (standalone installer, Win 7 Pro 64bit) - I import CSV - text in CP1250 separated with ";" and using "," as decimal separator. I correctly set the separated text importer, get correct point layer which has all text entries displayed correctly. I save this as shapefile, set the CP1250 encoding but when I load this layer in QGIS (choose CP1250 again during SHP loading) but get point layer with damaged text entries. I use this tool for a long time and did not have such problems. Did you have encounter such problems? thanks


title: Couldn't load SIP module (QGIS doesn't load Python)

text: I recently installed QGIS 2.4, but Python does not load. I am running on a Windows 64 bit machine. I receive the error message "Couldn't load SIP module Python support will be disabled". Traceback (most recent call last): File "", line 1, in ImportError: DLL load failed: %1 is not a valid Win32 application. I have another installation of Python at C:/Python27. It is the 32-bit version. When I remove C:/Python27 folder, QGIS works fine, so there is some linkage between the new QGIS installation and the C:/Python27 python interpreter. I can't permanently delete C:/Python because I need it for other applications, so I am wondering if there are any suggestion as to how to fix this. Thanks for any advice.


title: QGIS error when add WFS

text: I'm running QGIS 2.4.0-Chugiak on win7 64bit and am getting the following error while trying to add a wfs layer (using the simple Add WFS layer command): The layer http://copadesvr02.copade.neuquen.gov.ar/geoserver/wfs?request=getCapabilities&SERVICE=WFS&VERSION=1.0.0&REQUEST=GetFeature&TYPENAME=cite:Comunidades_Indigenas&SRSNAME=EPSG:22182 is not a valid layer and can not be added to the map I managed to add some layers just a while ago but now can't even add those same layers, and I just tried the same WFS service on gvSIG and was successful. Can anyone help me out? The WFS url is: http://copadesvr02.copade.neuquen.gov.ar/geoserver/wfs?request=getCapabilities


Recently I've been designing maps without navigation tools like the Zoom In/Out and Pan tools, or arrows to pan north, pan west, etc. Instead, most modern mapping APIs support the direct use of the map itself to handle navigation - drag the map to pan, double-click to zoom in, scroll the mouse wheel, shift-drag to zoom, etc. **Is there a set of standard accessibility guidelines for online maps?** Are there valid reasons for needing a "click to pan" control, or any of the other old-school navigation controls? For example, are they used by screen readers or other accessibility aids? I note that Google Maps has a Pan control, even though it's not strictly necessary. ![enter image description here](http://i.stack.imgur.com/nww0p.png) Thanks.


title: How to export MXD layout view to GeoTiff?

text: It my be old issue, but I can't find right solution. How to produce GeoTiff from layout view not from data view? Any simple solution? Thanks


title: Changing symbology as per scale reference in ArcMap?

text: I am using ArcMap 10 and I want to display my village parcel as per given scalewise symbology. I have one idea that copy and paste same layer. But my question is that in a single layer, is it possible to give different symbology as per the scale? e.g               Layer Name    symbol   scale     xyz            -----   100000     xyz            *****   200000


title: Integrating Maximo with Mapinfo?

text: I am software developer and have been working on database and Java as major platform for development. I am using Mapinfo for GIS analysis. Now, the system needs to show trouble tickets in Mapinfo that are achieved from Maximo software. How is it possible to connect these two application?


title: How to test if ArcGis server 10.1 is running correctly?

text: I've recently installed arcgis server 10.1, but I'm not sure if it installed correctly, or if the proper ports are open. Here is my reasoning and I may be off: 1) When publishing a new service(i.e newMap), the folder "arcgisserver\directories\arcgisoutput\newMap" remains empty, I'm assuming their should be a copy of the service in there now. 2) When trying to access the service through arcgis javascript api, I get a blank page. I copied and pasted the source code from the tab or window that pops up when you press "view" in the services thumbnail into a .html file. Then I replaced the url with the REST URL given to me under the services' capabilities options. I've also tried replacing localhost with my machine name and I still get a blank page.               http://localhost/myorg/rest/services/MyMapService/MapServer      If anyone can provide me with some working test code, that would be great, as I'm learning to use javascript api and server at the moment and I want to make sure its configured correctly. Thank you.


title: MapInfo - programmatically find if layer is hidden

text: I am working on an integrated mapping project. I would like a button that toggles a layer. I could hold the layers status in a variable and then show/hide it as appropriate but this seems a little hacky. Is there a way to find out whether a layer's "display" value is set to on or off?


title: SharpMap UK OS data

text: I've downloaded the UK raster and the corresponding tab file from the ordnance survey website and have displayed the image using the GDARaster class. however I am now stuck as I have no idea how to use the .tab file that came with the map and how I can display points on the map using long and lat. I can't seem to find a working example of how to do this can anyone help? Thanks.


title: Looking for QGIS introduction course in Los Angeles

text: To Whom It May Concern: I am in the process of launching a commercial real estate business and I wish to incorporate QGIS as a tool for keeping and analyzing demographic and market research information. I am new to QGIS and GIS in general and I need to find an instructor, seminar or class in the Los Angeles, California area to learn how to use QGIS. Also, I am an Apple Mac user and compatibility with Mac is important to me. Please advise. L. Carlos Simental


title: Access Trimble GPS on ArcGIS Engine App

text: Hopefully this isn't too obscure. I have a Trimble Yuma tablet running Windows 7. I need to build an application that allows viewing of points on a map and editing of a geodatabase (feature layers and tables). I'm hoping to use ArcGIS Engine SDK. BUT...how do I leverage the GPS of the Trimble within my Arc Engine app? Anybody done this with a Trimble before? Any code samples out there? Help!!!


title: Add or Import Metadata to a shapefile

text: Is there a way to add metadata to a shapefile in Qgis? I have the shape file and an xml metadata file and I want to import in shape file.


title: How to combine a vector and a raster file?

text: I want to combine a raster and a vector file. After combining, can I query on the resulting map to obtain a specific area? How?


title: Using zip file in html

text: I have my map files vdcar.js (about 27.79 MB) and vdcar.js.zip (about 9.79 MB) in same folder in server. Also have TestPopden.html (http://ngiip.gov.np/THEMATIC_MAPS/TestPopden.html) and .htaccess (empty file) in same location. I found that it takes maximum time (more than 30 minutes) to open the unzipped map file. What I need to do (with TestPopden.html and .htaccess files) to use zipped version of my map file?               <!DOCTYPE html>                 <title>Population Density</title>       <meta charset="utf-8" />              <link rel="stylesheet" href="leaflet.css" />       <link rel="stylesheet" href="POPCEN/Popcen.css" />


title: Changing from portrait to landscape layout?

text: I am trying to switch to from a portrait to landscape layout. I've already switched up the buttons in page and print set up and for some reason, the print layout icon shows my page still having portrait features in a non- completely filled landscape icon. I'm a graduate student in urban planning who still has his GIS training wheels on.


title: How to geocode US addresses using gisgraphy (400k)

text: Would someone please direct me in geocoding around 400k addresses using gisgraphy or openstreetmap? Thanks in advance, Aji.


title: Sextante in QGIS - downloaded, but does not show up

text: I downloaded Sextante, and I put it into the directory C:\users[my_user_directory].qgis\python\plugins\sextante-1.0.9\sextant I re-started the computer, but there is still no sign of Sextante in QGIS, neither in plug-ins manager, nor menu items. Do I have one directory layer too many in the file path above? Is that why it doesn't work? Thanks.


title: Converting dataset between coordinate systems

text: I'm working with NOAA's PRECipitation REConstruction over Land dataset, and am looking for an efficient way to sample the data by pixels on NASA's MODIS grid, using open source tools. The PRECL data uses 0.5 degree blocks, so I'll have quite a few repeated values. I'm currently using a mapping generated by ArcGIS that specifies the array index within the PRECL data for every MODIS pixel. This works (and it's fast, as I'm just doing a series of array lookups), but it'd be nice to get away from that extra data set, especially if there were some way to do the transformation in code, without Arc. Thanks in advance, all!


title: How to know contents of column names in GRASS?

text: I am using GRASS to extract a specific area from a vector map. Here is the script that I am running: v.extract PHL_map out=makati where="AREA_NAME = 'Makati" The problem is, I don't know what are the values of the columnnames. I can only view the columnnames but I cant view the values under the columnnames


title: CQl Filter Don't work with WFS layer in Open Layers

text: I used a WFS layer from geoserver in open Layers it's shown on MAP but when I want to use a filter it doesn't work I don't know where the error this is the code :               <script type="text/javascript">     OpenLayers.ProxyHost = "/cgi-bin/proxy.cgi?url=";          function init() {     //set up a save strategy     var saveStrategy = new OpenLayers.Strategy.Save();          //set up projections     // World Geodetic System 1984 projection     var WGS84 = new OpenLayers.Projection("EPSG:4326");        // WGS84 Google Mercator projection     var WGS84_google_mercator = new OpenLayers.Projection("EPSG:900913");             //Initialize the map     //creates a new openlayers map in the <div> html element with id="map"     var map = new OpenLayers.Map ("map", {        controls:[     //allows the user pan ability     new OpenLayers.Control.Navigation(),     //displays the pan/zoom tools                        new OpenLayers.Control.PanZoom(),     //displays a layer switcher     new OpenLayers.Control.LayerSwitcher()        ],     //projection: WGS84_google_mercator,     displayProjection: WGS84     } );          // map extent     var mapextent =     new OpenLayers.Bounds(-7.70912,33.65164,-7.39876,33.52465).transform(WGS84, WGS84_google_mercator);          //base layers     var openstreetmap = new OpenLayers.Layer.OSM();     var google_maps = new OpenLayers.Layer.Google(     "Google Maps", {     numZoomLevels: 20     }     );     var google_satellite = new OpenLayers.Layer.Google(     "Google Satellite", {     type: google.maps.MapTypeId.SATELLITE,     numZoomLevels: 20     }     );          //wfs-t editable overlay     var wfs_layer = new OpenLayers.Layer.Vector("Editable Features", {     strategies: [new OpenLayers.Strategy.BBOX(), saveStrategy],     //projection: new OpenLayers.Projection("EPSG:26910"),     protocol: new OpenLayers.Protocol.WFS({     version: "1.1.0",     // loading data through localhost url path     url: "http:// localhost/geoserver/wfs",     featureNS :  "http://www.opengeospatial.net/cite",     maxExtent: mapextent,     // layer name     featureType: "patrimoine",     // geometry column name     geometryName: "geometry",     schema: "http:// localhost/geoserver/wfs/DescribeFeatureType?version=1.1.0&;typename=cite:patrimoine"     })     });          map.addLayers([openstreetmap, google_maps, google_satellite, wfs_layer]);     map.zoomToExtent(mapextent);          }          function filter(){                      var filterType ="cql";                 //var filter = document.getElementById('zone').value;                 // by default, reset all filters                 var filterParams = {cql_filter: null};                 if (OpenLayers.String.trim(filter) != "") {                     if (filterType == "cql")                          filterParams["cql_filter"] = "batiment='Bastion de la sqala'";          // In this line my parametre it' batiment I made a test on the batiment wich has a value "Bastion de la sqala"                                     }                 // merge the new filter definitions                 mergeNewParams(filterParams);                  }                  function mergeNewParams(params){                 Region.mergeNewParams(params);                  }      and whene I type in console map.layers[2].features it show me an error **Cannot read property '2' of undefined** And Thanks for the help


title: How do I label my map with attributes?

text: help is always appreciated. I have a state map in QGIS 2.0 with all of the county names and some other info listed in the attributes table. I would like to have the name of the counties and one of the attributes visually displayed on the map in each respective county. This should be relatively simple...but I can't seem to figure out how to do this. Thanks very much.


title: Not all features available in OSM exports in TileMill

text: I'm using TileMill and OSM Bright to style an OpenStreetMap for hiking. The problem is that OSM Bright does not import all the OSM features for styling. I believe it's possible to add features in the osm-bright.osm2pgsql.mml file, but I'm not certain how. Does anyone know how to add, for example, path labels? Thank you.


title: finding the index of a field with its name using arcpy

text: I wan to extract the index of a field with its name using arcpy. it seems what I have written is so dumb. is there any simpler method with out using for- loop.                   import arcpy              fieldList = arcpy.ListFields(self.address_ref_layer1)         self.field_name = []         for field in fieldList:             self.field_name.append(field.name)         #self.remove_pseudo_nodes()         print self.field_name         print self.field_name.index(u"NAMN1")      it could also be written:                   fieldList = arcpy.ListFields(self.address_ref_layer1         name_index = [f.name for f in fieldList].index(u"NAMN1")      the ListFiled function returns list of field objects with many properties as well as name but no function to return the index of the field which has a specific name.


title: How do you set sea level elevation to white in terrain analysis?

text: I have ASTER files I am working with to create a shaded relief background for some other vector shapefiles. In experimenting with different color-elevation options, I am unable to have the Terrain Analysis tool make the surrounding ocean white in color. I have tried manually adding a white color for elevation values 0.0 to 0.0, but no dice. Going to the Table of Contents->Properties and selecting a different color scheme just replaces the Terrain Analysis colors and I have to re-do Terrain Analysis (although that is likely because I don't know how to 'undo' the color change when the 'undo' option is not available). Does anyone have a solution?


title: Import data, analyze spatially, and export list of points using Google Earth

text: I am trying to:   1. import data into Google Earth from an outside data source.   2. Then overlay a map and isolate which datapoints fall into that area.   3. I would like to then export a list of the data points to give to a person that works with our company. Can I do this in Google Earth Pro? Can you direct me to a resource that would walk me through this? Is there a better solution for doing this type of work if it can't?


title: Normalization of Kernel Density Estimate rasters in ArcGIS 10

text: I've been told by several GIS instructors that the pixel values of a Kernel Density Estimate raster output by the point KDE tool in ArcGIS 10.1 represents the probability of encountering a point, not the density of points at that location. I've been told that the KDE surface is normalized to a probability density surface by dividing the surface by the total number of points used to make it. My problem is that there is no reference to this normalization procedure in the ArcGIS help files at all, with the possible exception to a reference to the paper by Silverman (1986) from which the quadratic kernel function used in ArcGIS 10.1 is derived. Can any one tell me if this is the case? In addition, if this is the case, if you select a field value instead of just using the point location to build the KDE surface, does this normalization procedure still occur, and how would you interpret the surface?


title: Convert latitude and longitude into northings and eastings

text: Hi I have the following UK postcode **dy8 3xt** and know that the latitude and longitude is:-               54.452772     -2.156082      I also know that the Eastings, Northings for the postcode is:-               389490     283880      However I am struggling to find the equation that converts lat/long to northings and Eastings, I would prefer to have the equation in both in jScript and c# (I am being greedy)! Can anyone help? **EDIT** Thanks for your help so far guys, I am starting to learn something here esp. the terminology... Some more info, if you click on this link you can see the results I am looking for. The postcode I entered projects to lat/lng using WG S84 and the grid ref projects to OSGB. So my question is how is this done? **WHAT I LEARNT** Thanks to all that answered, I finally got led to here which I can confirm works great


title: When creating annotation from labels is it possible to feature-link another annotation layer?

text: I have an annotation layer that is currently being used in the flex viewer and I would like to add a secondary annotation layer/line so that it sits directly underneath the first. I am able to link the new annotation layer to the polygon I am creating it from, but I wonder if there is a way to also link it to the annotation layer I already have; as I am having trouble getting the annotation to sit underneath the first one. I tried setting the offset, but as it offsets from either the center or the edge of the block, the new annotation does not line up within the top right corner underneath the original annotation. Editing is off the table as I would have some 30,000 records to manipulate manually. Is there any way to link these labels to the original annotation layer before I turn them into annotation?


title: Spatial clustering with PostGIS

text: I'm looking for spatial clustering algorithm for using it within PostGIS- enabled database for point features. I'm going to write plpgsql function that takes distance between points within the same cluster as input. At the output function returns array of clusters. The most obvious solution is to build buffer zones specified distance around the feature and search for features into this buffer. If such features exist then continue to build a buffer around them, etc. If such features not exist that means cluster building is completed. Maybe there are some clever solutions?


title: Label pie chart in QGIS

text: 1.How can I label the pie charts in QGIS so the labels are written under the diagrams? My pie-charts are representing localities: such as villages, towns, cities. So, I want that their names appear under the diagrams.I want them to look like this image ![I want them to look like this image](http://i.stack.imgur.com/y4YjW.png)   1. How can I move manually the labels of the pie-chart diagrams? When I work with this type of diagrams, the MOVE LABEL icon is not active. Thanks a lot!


title: How to host web map app on ArcGIS.com?

text: I have created a web map app on **arcgis.com**. I have installed ArcGIS 10.2 for Server, and created a services and published it on **arcgis.com**, but it is local server. I want to any one opening my web map app from any other PC not my PC **Here** is my web map app.


title: OpenStreetMap Tile Server with Custom Projection Systems

text: I am trying to install the OSM tile server on a local machine with a custom projection system like epsg:3573. For the server installation I followed this tutorial (http://switch2osm.org/serving-tiles/building-a-tile-server-from- packages/) and it works fine with the default settings. For making changes, I played around with osm.xml and changed the srs information in there but still no luck. It would be appreciated if anyone can give me some hints on customizing the OSM tile server with a custom projection system. Thanks, Ebrahim


title: export spatial information from table that is joined to a shapefile in QGIS 2.0

text: I have a table with data values per district for different years. The table looks like this: district name / year / value. The same district is thus present in the data multiple times for different years. In QGIS I joined this table to a districts polygon Shapefile. THis creates a table with for every year a row with district name and value. Can I export this table in a way the spatial information of the district is stored in the table? Thanks in advance!


title: How to improve performance using spatial indexing with STIntersects()

text: I am working on Spatial database. I need to check intersection of state geometry with geometry inside database. I am presently using the query as                SELECT a.ID, a.OBJECTID, a.MSID, a.FACNAME,Convert(varchar(max),a.ShapeGeo),b.State_Name FROM XFAcilityBackUp26June a with(index(testIndexSpatial)),us_stateTest b WHERE a.shapegeo.STIntersects(b.the_geom.MakeValid()) = 1 and b.State_Name = 'Texas'      It is functioning properly but it takes more time(more than 10 minutes) to execute. I have tried indexing also but no luck with it. Is there any alternative way to improve the performance. Thanks in Advance


title: Installing QGIS 2.0 with the newest gdal and geos libraries

text: I want to install QGIS 2.0 on Ubuntu 12.04 LTS. I followed these steps until the final installation command. Before I compiled gdal 1.10, geos 3.4.2 and grass 6.4.3 on my own. In the synaptic package manager I marked qgis (2.0.1-2+precise1) for installation and I can see that the older libraries gdal 1.7, geos 3.3.8 are marked as well. How can I build QGIS with the new compiled libraries?? I want to avoid that the old libraries will be installed as well.


title: Which value should I assign my digitalized features?

text: I have a scanned an old document that contains geochemical interpolated values (no way to work with the digital file). This is the original image: ![enter image description here](http://i.stack.imgur.com/fDZDA.jpg) And this is the digitalization I made from it: ![enter image description here](http://i.stack.imgur.com/DhDcf.jpg) Now, given that the original image is an interpolation of geochemical measures taken in field, which values should I give to each features of my polygon? For example the crescent moon shape in the upper left corner has a maximum value of 100 (in the original file). Is the half of the original value a good idea (so 50 in this case)? Someone has any suggestion? I hope I made myself clear. p.s. software used: QGIS 2.0


title: How to configure Geoserver WMS to make it printable from QGIS/ArcGIS?

text: I have installed geoserver 2.3.1 - I am using this to create a WMS of some data held in PostGIS. I can view the WMS in QGIS and ArcGIS, but it won't print from ArcGIS (blank sheet) and in QGIS I get the message: Some WMS servers (e.g. UMN mapserver) have a limit for the WIDTH and HEIGHT parameter. Printing layers from such servers may exceed this limit. If this is the case, the WMS layer will not be printed So I assume I have to set something in Geoserver or on the server to allow larger print sizes - where can I do this?


title: Setting Zero values back to NULL

text: I'm working with ArcGIS 10.2 for Desktop and OS Windows 7. I need to set numeric vales from Zero back to NULL, but I can't from the field calculator. I'm working with a vector file and I need to display NULL in one column that shows 0 (Zero) as default value. -------- Edit-------   I am adding the error messages that he received ![Error Message 1](http://i.stack.imgur.com/IJEjy.png)   Translation: An error occurred during the process. Please consult the Geoprocessing Results Window to obtain more details.   ![Geoprocessing Results Window](http://i.stack.imgur.com/Ywhel.png)   Error at bottom, translation: The values of the selected column can not be set to null because the field in the database is defined as "Not Null".


title: Manually building a tile server

text: I use this instruction http://switch2osm.org/serving-tiles/manually-building- a-tile-server-12-04/ I have the problem on the step "Starting your tileserver"               sudo mkdir /var/run/renderd     sudo chown username /var/run/renderd     renderd -f -c /usr/local/etc/renderd.conf      The problem:               renderd[4620]: An error occurred while loading the map layer 'default': Could not create datasource. No plugin found for type 'shape' (searched in: /usr/lib/mapnik/input) (encountered during parsing of layer 'world' in map '/home/oleg/src/mapnik-style/osm.xml')     renderd[4620]: An error occurred while loading the map layer 'default': Could not create datasource. No plugin found for type 'shape' (searched in: /usr/lib/mapnik/input) (encountered during parsing of layer 'world' in map '/home/oleg/src/mapnik-style/osm.xml')     renderd[4620]: An error occurred while loading the map layer 'default': Could not create datasource. No plugin found for type 'shape' (searched in: /usr/lib/mapnik/input) (encountered during parsing of layer 'world' in map '/home/oleg/src/mapnik-style/osm.xml')     renderd[4620]: An error occurred while loading the map layer 'default': Could not create datasource. No plugin found for type 'shape' (searched in: /usr/lib/mapnik/input) (encountered during parsing of layer 'world' in map '/home/oleg/src/mapnik-style/osm.xml')     renderd[4620]: DEBUG: Got incoming connection, fd 8, number 1     renderd[4620]: DEBUG: Got command RenderPrio fd(8) xml(default), z(0), x(0), y(0)     renderd[4620]: Received request for map layer 'default' which failed to load     renderd[4620]: DEBUG: Connection 0, fd 8 closed, now 0 left     renderd[4620]: DEBUG: Got incoming connection, fd 8, number 1     renderd[4620]: DEBUG: Got command RenderPrio fd(8) xml(default), z(0), x(0), y(0)     renderd[4620]: Received request for map layer 'default' which failed to load     renderd[4620]: DEBUG: Connection 0, fd 8 closed, now 0 left     renderd[4620]: DEBUG: Got incoming connection, fd 8, number 1     renderd[4620]: DEBUG: Got command RenderPrio fd(8) xml(default), z(0), x(0), y(0)     renderd[4620]: Received request for map layer 'default' which failed to load     renderd[4620]: DEBUG: Connection 0, fd 8 closed, now 0 left      Update 1: One of the problem was the file "renderd.conf". here is file without error.               [renderd]     socketname=/var/run/renderd/renderd.sock     num_threads=2     tile_dir=/var/lib/mod_tile     stats_file=/var/run/renderd/renderd.stats               [mapnik]     plugins_dir=/usr/local/lib/mapnik/input     font_dir=/usr/share/fonts/truetype/ttf-dejavu     font_dir_recurse=1          [default]     URI=/osm_tiles/     TILEDIR=/var/lib/mod_tile     XML=/home/map/src/mapnik-style/osm.xml     HOST=localhost     TILESIZE=256     CORS=*      But I have another errors:               renderd[8967]: An error occurred while loading the map layer 'default': Shape Plugin: shapefile '/usr/local/share/world)boundaries/shoreline_300.shp' does not exist (encountered during parsing of layer 'world' in map '/home/oleg/src/mapnik-style/osm.xml')     renderd[8967]: An error occurred while loading the map layer 'default': Shape Plugin: shapefile '/usr/local/share/world)boundaries/shoreline_300.shp' does not exist (encountered during parsing of layer 'world' in map '/home/oleg/src/mapnik-style/osm.xml')      Update: The question is solved


title: 2 shape files plots together in a GIS, but not in R

text: I tried plotting these two files and a set of points (dropbox folder) in QGIS and MapInfo. They align pretty good (not perfect though). But in R, only the data from www.naturalearthdata.com can plot along with other things like the points. My 'own' shape files plots nicely, but always alone (not adding on top of or "allowing" points on top of) and... My own shape files are rotated (clipped from a larger file where the center has north up) and I can't figure out how to save the applied rotation in MapInfo (i've tried exporting numerous ways). Maybe this is causing my issue?               library(rgdal) library(raster)     #Loading "my" shape files:     shp.WashingtonLand<-readOGR(".","WLmap2_polyline")     plot(shp.WashingtonLand)          #Loading naturalearthdata.com data for only Washington Land     shp.greenland<-readOGR(".","GRL_adm2")     ext.washington.land <- extent(-68, -58, 79.9, 81.2)     shp.clip.WashingtonLand <- crop(shp.greenland, ext.washington.land)     plot(shp.clip.WashingtonLand)          #data points     SamplePoints <- read.csv2("SamplenumberAndCoordinates",header=TRUE)     locs <- subset(SamplePoints, select = c("Latitude", "Longitude"))     coordinates(locs) <- c("Longitude", "Latitude")     plot(locs, col="red", add=T)      Hope it makes sense, it took a while getting all the things together and simplified appropriately.


title: Openlayers and Geoserver: OSM/Google Maps and WMS Overlay

text:
I am experiencing projection problems with the following set up:   * GeoServer with ESPG:4326 based shapefiles (*.shp)   * OpenJump to edit my shapefiles   * Openlayers with EPSG:900913 base layers   * GeoServer provides the shapefiles as WMS to OpenLayers.  Here is my source code:                             <title>Openplanning mit Geoserver</title>         <meta http-equiv="content-type" content="text/html; charset=ISO-8859-1">

                  It seems that GeoServer is converting from ESPG:4326 to EPSG:900913 internally. Here you can see a sample request for Vienna (which is sent by my OpenLayers application to GeoServer):               http://localhost:8888/geoserver/wms?LAYERS=topp:ma_vienna&SRS=EPSG:900913&TRANSPARENT=TRUE&SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&STYLES=&FORMAT=image/png&BBOX=1800244.8899219,6116147.6705496,1839380.6483984,6155283.4290261&WIDTH=256&HEIGHT=256      My problem: the requested file looks distored, as you can see here: ![http://i.stack.imgur.com/UcWBQ.png \(sorry, I had to post a link again after my edit ...\)](http://i.stack.imgur.com/2F32X.png) In case I open the WMS layer in the OpenLayers preview of GeoServer, it looks ok: ![http://i.stack.imgur.com/P8IMj.png \(sorry, I had to post a link again after my edit ...\)](http://i.stack.imgur.com/dLYus.png) The question is: Why is my shapefile distorted in my OpenLayers application, but not in the GeoServer preview? As you can see, the GeoServer preview uses EPSG:4326 and my OpenLayers application ~~EPSG:900913~~ (Update: EPSG:4326 as well; mouse coordinates are displayed as latitude/longitude). But, however, as GeoServer provides a WMS interface with EPSG:900913 capabilities, shouldn't it correctly reproject my vector sources? Thanks for you replies in advance, Rainer.


title: How to reproject layers (using GeoServer and OpenLayers) in order to overlay on the google maps

text:
> **Possible Duplicate:**   >  Openlayers and Geoserver: OSM/Google Maps and WMS Overlay I am writing a code using OpenLayers. In this code, I want to display different layers in different zoom levels. Up until here, everything works. my problem began when I tried to add Google Maps as baselayers. My layers' projection system is not consistent with Google layers, so they should be reprojected. I tried to use Proj4js.transform method but it didn't work. My layers are in EPSG:3006(Sweden) and Google layers are in EPSG:900913. The following code results in a projected layer shifted to right. It means that the re-projection did not work properly. Please let me know what is wrong with my code.               <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN"     "http://www.w3.org/TR/html4/loose.dtd">                        <title>OpenLayers map preview</title>         <meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">                  <link rel="stylesheet" type="text/css" href="http://gisl.nateko.lu.se:8080/geoserver/openlayers/theme/default/style.css">                                                                                Map viewer using Geoserver and MapLayers

location

                 --> Click on the map to get feature info
            this code is buggy, since I tried to test every possible solution. So maybe it is necessary to remove some parts. for instance, remove the Zoom levels part which is not really my problem.


title: Geoprocessing tool not returning features when used in ArcGIS Viewer for Silverlight

text: I can already guess your first thought: "why the hell is this guy using Silverlight stuff, come on!". Unfortunately, that's what we have decided to go with. Here's my situation. I'm producing a very simple Geoprocessing Service that takes an input feature class (i.e feature set when using the dataset in a geoprocessing service) and uses the MakeFeatureLayer tool and does a query (within the MakeFeatureLayer tool) and produces an output. This tool successfully publishes and works as a service in ArcMap. However, when I use the service in the Silverlight viewer, it produces and empy output (the output is in the TOC, but without features). If I use the MakeFeatureLayer tool without any queries, it works. I've also tried the following to get selected features as an output... 1) MakeFetureLayer -> SelectLayerByAttribute -> CopyFeatures 2) IterateFeatureSelection None if these have worked. It seems that whenever a query of some type is used, it will produce an empy output, but only in the Silverlight viewer. Anythoughts as to why the selection is not working. Model Examples. ![Simple use of the MakeFeatureLayer tool](http://i.stack.imgur.com/g8nIU.png) ![Model Properties](http://i.stack.imgur.com/I5SPK.png) ![MakeFeatureLayer tool with SelectLayerByAttribute](http://i.stack.imgur.com/DF6CF.png) ![enter image description here](http://i.stack.imgur.com/VLPDH.jpg)


title: Converting lat lon into co-ordinates (in ft of m)

text: I have a set of latitude and longitude of a series of points which specific the trajectory of a vehicle. If I consider the first point as the origin (0,0), is it possible to convert the trajectory into a set of points with co- ordinates as (X ft, Y ft)? If yes how? I would prefer to have the mathematical formula if possible. I also have a related question. Does the conversion depend on the map projection I am interested in? If yes, the projection I would be interested in is Mercator projection with datum of WGS84.


title: produce a grid of geographic coordinates within 1 mile of each other with QGis or PostGIS

text: For a given county (bounded box, lets say), I need to get a network of gps points such that each point is 1 mile from its nearest neighbor point in each cardinal direction. The result of this initial calculation would be an array of gps coordinates, which if displayed on a map would look like the vertices on a piece of graph paper. With the one mile radius requirement, this process is arduous by hand for reasonably sized areas.. I'm using an API that takes point data combined with a buffer radius value to find locations. The plan is to make the buffer radii something like 0.7 miles in order to ensure 100% saturation of a given area.


title: WMSGetFeatureInfo Not Returning Results

text: I am having a problem getting the results back from a WMSGetFeatureInfo call to my GeoServer. Using OpenLayers 2.12, GeoServer 2.4.0, Windows Server2008, and Microsoft IIS 7.5 with Application Request Routing Cache and URL Rewrite. I do NOT use a proxy because I use URL Rewrite. The URL Rewrite takes my incoming requests and formats them to avoid the cross domain/same origin problem. I have been successfully querying WFS data using this setup for 2 years but when I try to use the WMSGetFeatureInfo it fails. Here is my process: 1) I serve multiple WMS layers from GeoServer to my web page. 2) The user selects a WMS layer from a dropdown. This sets the variable "selLayer" so my WMSGetFeatureInfo call is pointing at the layer the user wants to get information about. 3) When the user clicks on the layer it fires this code which creates a OpenLayers WMSGetFeatureInfo control with an event listener that should accept the return text & xmlHttp request and display them in a div. :                    var infoControl = new OpenLayers.Control.WMSGetFeatureInfo({              url: 'http://xx.xx.xx.xx:8080/geoserver/wms',              title: 'Identify features by clicking',              layers: [selLayer],              queryVisible: true,              //eventListeners: { getfeatureinfo: function (event) { alert("TEXT: " + event.text + ", XY: " + event.xy + ", REQUEST" + event.request + ", FEATURES: " + event.features); } }               eventListeners: { getfeatureinfo: function (event) {                  document.getElementById("divAttributes").innerHTML = event.text + "Test Text";                  alert("PUT TEXT HERE!!!");              }               }          });                    map.addControl(infoControl);          infoControl.activate();      But, when the script runs I do not get return values for the event.text. When I run this in FireFox (with FireBug) I can copy the GetFeatureInfo request into a new browser window and I get the desired return values. So my question is: Why am I not getting the return values when I run this on my map webpage but I do get results by posting the request into a new browser window. Any help would be appreciated. Thanks!


title: Can QGIS handle .adf attribute data?

text: I have been using QGIS for about 4 months now and its been great to learn a new package. I have been trying to open a source raster (*.adf) files in QGIS. I can display the raster set itself, however, the appropriate attribute data fails to join the spatial dataset. Can QGIS display the full source raster for ArcGIS? Or is the a 'proprietary' format that can only be fully displaying in Arc? Thanks for any help.


title: Why would running LocateFeautresAlongRoutes GP tool return null?

text: I'm currently using the `LocateFeaturesAlongRoutes` tool wrapper for C#/.NET for a custom tool I'm developing for ArcMap 9.3.1. I've been able to run this tool in ArcMap outside of code in a tool with the same files as parameters. Unfortunately when I put it in code I keep getting null for the `IGeoProcessorResult`. Is there any reason why I would get a null object for the result instead of just saying it failed? Also, is there anything wrong with the below code? The relevant block of code is below.               LocateFeaturesAlongRoutes lfar = new LocateFeaturesAlongRoutes();     lfar.in_features = @"C:\GIS_Temp\CONGPROP.mdb\CONGPROP2";     lfar.in_routes = @"C:\GIS_Temp\Routes.mdb\routes";     lfar.route_id_field = "ROUTE_ID";     lfar.out_table = "eventTable.dbf";     lfar.out_event_properties = "ROUTE_ID LINE FMEAS TMEAS";     lfar.in_fields = "NO_FIELDS";          Geoprocessor geoprocessor = new Geoprocessor();     geoprocessor.TemporaryMapLayers = false;     geoprocessor.OverwriteOutput = true;     geoprocessor.AddOutputsToMap = true;     IGeoProcessorResult result = (IGeoProcessorResult) geoprocessor.Execute(lfar, null);      Any suggestions?


title: Feeder bus service

text: First of all, a little background. I work for a regional transit agency. We are making a "diagnostic" about our feeder bus service. We'd like to know what proportion of our users could take the bus to go to the train station instead of taking their car. It's been done several time in the pass, but we are now using gtfs as our main datasource so we have to rethink our methodology. ![typical feeder bus service](http://i.stack.imgur.com/bRZFO.png) To be considered "feeding" the train, a bus route must have a stop within a certain distance from a train station (red buffers). Also, synchronicity with the train service is very important because if your bus arrives at the train station half an hour before the train, the waiting time is too long and you'll want to sleep 20 more minutes in the morning and take your car. Let say you take Line A (Blue) at stop 12. You step off the bus at stop 13. The bus arrives at stop 13, which is the stop to go to the Train Station #1 5 minutes before the train. That's very good. That would mean that everybody that takes that bus route at a stop from 1 to 13 included would arrive 5 minutes before that train. Then, the train, going through a very densely populated area with lots a schools and crossings is forced to reduce its speed a lot. In the meantime, the bus picks up passengers at stop 14 to 17 and arrives at Train Station #2 10 minutes before that train. So passenger taking the bus at stops 14 to 17 would all have a waiting time of 10 minutes once arrived a the train station. So, along that bus line, passenger taking the bus at stop 1 to 13 have a waiting time of 5 minutes when those taking the bus at stops 14 to 17 have a waiting time of 10 minutes. Line B, on the other side of the track, pass near Train Station #1, but its stops are too far to be considering "feeding" Train Station #1. It arrives at Train Station #2 7 minutes before the train (doest it for every train during the morning rush hour; it's very well synchronized). So passengers along Line B, taking the bus everywhere from stop 1 to 59 would have a waiting time of 7 minutes. Now, my question. Once I have determined that stops LineA.13 and LineA.17 are feeding my train (it has been done spatially, in PostGIS), and that the waiting time when taking the bus at a stop before #13 is 5 minutes but those after have a waiting time of 10 minutes, how can I assign the waiting time to all the stop before them? I'd like to do it in Postgres/PostGIS (pl/pgsql or pl/python), but I can use pure python (OS or arcpy) as well. I could, I think, loop backward. So, once I found a stop that fits (here LineA.17), assign the same waiting time to stop 16, then 15... until I found another stop that fits my criteria (LineA.13) and then assign the the rest of the stops, the same waiting time as 13. I have no idea how to create such a loop though. I don't think that I can do it in SQL so I'd have to use a procedural language in PostgreSQL. I had an idea of using pgRouting to find the route between each feeder stops so that way, Line A would be split in two (stops 1 to 13 and then 13 to 17). Would that be easier? The next step, will be to use pgRouting to calculate the driving time from all the stops that have a waiting time (sorry to LineA.18 and over!) and compare that with the schedule of the bus to calculate competitivity (is it taking 5 minutes more in bus that in car?) Any ideas? I normally post a long work-in-progress script to show the effort I've made so far, but I'm stuck! Thanks!


title: trouble with code block execution in python script

text: I am trying to calculate values for a field in a python script, but keep getting an error that I don't understand:               ExecuteError: ERROR 000539: Error running expression: Fron000000F1      Traceback (most recent call last):       File "<expression>", line 1, in <module>     NameError: name 'Fron000000F1' is not defined     Failed to execute (CalculateField).      I can't figure out where the name 'Fron000000F1' is coming from to try and track this problem. I have tested the expression in arcmap using the field calculator with no problems. What is wrong? Code is below:               import arcpy, os, csv     arcpy.env.overwriteOutput = True          out_fc_pm = r'C:\temp\temp.gdb\SignDetect__SOL680_DD_with_PM'          def calc_id(pm,pict_fld):       import re       pm_patt = re.compile('\d*\.*\d*')       pm_match = re.findall(pm_patt,pm[5:])[0]       pm = pm_match.replace('.','')       pm_final = pm.zfill(6)       dir_patt = re.compile('[A-Z]\d*')       dir_match = re.findall(dir_patt,pict_fld)[0]       id = pict_fld[:5]+pm_final+dir_match       if id[-1:].isdigit():         return id       else:         return id+'1'          arcpy.AddField_management(out_fc_pm,'OSMI_ID','TEXT')     arcpy.CalculateField_management(out_fc_pm,'OSMI_ID',calc_id('!POSTMILE!','!Front!'),'PYTHON_9.3')


title: More Efficient Spatial join in Python without QGIS, ArcGIS, PostGIS, etc

text: I'm attempting to do a spatial join much like the example here: Is there a python option to "join attributes by location"?. However, that approach seems really inefficient / slow. Even running this with a modest 250 points takes almost 2 minutes and it fails entirely on shapefiles with > 1,000 points. Is there a better approach? I'd like to do this entirely in Python without using ArcGIS, QGIS, etc. I'd also be interested to know if it's possible to SUM attributes (i.e. population) of all the points that fall within a polygon and join that quantity to the polygon shapefile. Thanks in advance for any help! **Here is the code I'm trying to convert. I get an error on line 9 (poly['properties']['score'] += point['properties']['score']) that says TypeError: unsupported operand type(s) for +=: 'NoneType' and 'float'. If I replace the "+=" with "=" it runs fine but that doesn't sum the fields. I've also tried making these as integers but that fails as well.               with fiona.open(poly_shp, 'r') as n:        with fiona.open(point_shp,'r') as s:         outSchema = {'geometry': 'Polygon','properties':{'region':'str','score':'float'}}         with fiona.open (out_shp, 'w', 'ESRI Shapefile', outSchema, crs) as output:             for point in s:                 for poly in n:                     if shape(point['geometry']).within(shape(poly['geometry'])):                           poly['properties']['score']) += point['properties']['score'])                         output.write({                             'properties':{                                 'region':poly['properties']['NAME'],                                 'score':poly['properties']['score']},                             'geometry':poly['geometry']})


title: How to find the shortest route between two points based on cost surface?

text: I have two locations: (A) starting point and (B) designation. Also I have a slope raster. I want to find the cheapest path between these two locations based on the slope (it is more expensive to travel over steep slope). What tools in which order do I have to use in order to get that (either ArcGIS or Grass GIS)?


title: How to import .csv file to CartoDB using SQL?

text: How to import .csv file to CartoDB using SQL, the data given below : POINT | TOP_dm | BOTTOM_dm | X_RD | Y_RD 35AP7600 | 4418 | 4368 | 262053 | 467223 35AP7601 | 4267 | 4217 | 262000 | 467143 35AP7602 | 4358 | 4308 | 261780 | 467148 35AP7603 | 4209 | 4159 | 262098 | 467243


title: Raster iterator on module builder

text: I'm trying to iterate through a folder of raster via "iterate raster" , then the output of the iterate raster is the input for " extract Sub data-set " since I want to choose the cloud mask Sub data-set. The output of extract Sub data-set is the input for "Make Raster Layer" since I want to have a raster with one band out of six from the original raster. The module isn't working!Any help please?


title: How to create a list and append objects to it through a loop

text: Still struggling with this script. Now I need to create an arpy.AddMessage that says "Projected: with the list of feature classes that were projected in the loop" i.e. Projected: filename, filename, filename etc. I've been at this for 3 days and know I am missing something stupid. Any help would be appreciated. Here is my code again:               # Get the spatial reference of the Feature Class from folder A     # Examine all Feature Classes in folder B and report if their     # spatial references match the Feature Class in folder A          import arcpy     arcpy.env.overwriteOutput = True          #try:     #set up the paths     #targetFC = r"C:\MGIS\geog485\Lesson2\Lesson2\CountyLines.shp"     #folderToExamine = r"C:\MGIS\geog485\Lesson2\Lesson2"          targetFC = arcpy.GetParameterAsText(0)     #print "target FC path is "+targetFC     folderToExamine = arcpy.GetParameterAsText(1)     #print "folder to examine"+folderToExamine          #get spatial reference for the target feature class     targetDescribe = arcpy.Describe(targetFC)     #print targetDescribe     targetSR = targetDescribe.SpatialReference     #print targetSR     targetSRName = targetSR.Name     #print "Target feature class is "+targetSRName          # Get a list of my feature classes     arcpy.env.workspace = folderToExamine     listOfFCs = arcpy.ListFeatureClasses()     #print listOfFCs          #Loop through the list of FCs     for currentFC in listOfFCs:         #print currentFC         #Read the spatial reference of the current one         currentFCDescribe = arcpy.Describe(currentFC)         currentFCSR = currentFCDescribe.SpatialReference         currentFCSRName = currentFCSR.Name         #print currentFCSRName              if currentFCSRName != targetSRName:             print "Spatial references don't match"         else:             print "Spatial references do match"            if currentFCSRName == targetSRName:             #skip            continue         else:            # Determine the new output feature class path and name             outCS = currentFC[:-4] +"_projected.shp"             #Reproject datasets that are different from target file             arcpy.Project_management(currentFC, outCS, targetSR)                  #Create a message stating which files were projected             projFCs = outCS[:-14] +".shp"             #print projFCs             projFCList = arcpy.ListFeatureClasses(projFCs)             x = ", ".join(projFCList) + ", "             arcpy.AddMessage("Projected " + str(x))


title: Getting a square mile bounding box anchored by coordinates?

text: Let's say I have a lat/lng coordinate pair. I want to determine the coordinates that would create a square mile box. The existing coordintae pair would be the top left of the box, and the coordinate pair I'm trying to calculate would be the lower right. What's the best way to handle this? A solution using GeoDjango would be great, but PostGIS or even a mathematical formula would be good too.


title: Setting up drawing order for point shapefile

text: I am having a difficulty with creating polyline feature from a point dataset that represents a boundary. This dataset consists of 5084 points which are randomly ordered. I have no Network Analyst extension available. I tried to use X Tools Pro and 'Make Polylines from Points' tool, but not having the right drawing order set up the output is just a random cross of polylines. I have also tried 'Tracing' in the edit session but I've found that it wasn't possible to trace points (only lines or polygons) - please correct me if I'm wrong.


title: What is the current Web Mercator projection code?

text: The web mercator projection, popularised by Google Maps, seems to be given a new EPSG code every couple of years. EPSG:900913 (Google in calculator text) was an unofficial code EPSG:3785 - is the projection I currently use for my datasets However this blog post suggests that the code is now EPSG:3857. This projection can be found on the EPSG site but it seems to also use the code SR-ORG:6864 and claims _"it is not a recognized geodetic system: see WGS 84 / World Mercator (CRS code 3395)."_ So what is the official code to use?


title: How can I export raster data from Arcmap to an R  *.grd file format?

text: I want to work with raster data in R, but I have problems exporting my raster dataset from ArcMap. By "export data"/GRID I get a file with an *.aux extension. What I need is *.grd. Can anybody help me out?


title: Trouble making heat map - red square

text: I'm a beginner to QGIS and I'm having trouble making my heat map. I'm using QGIS 1.8 and Win7. I imported a CSV file with around a 1000 points (X Y coords - WGS84). When I try to make a heat map (radius 10, decay 0,1), all I get is a red square, even after changing color map to pseudocolor. Any ideas? Thanks


title: How to convert high precision dataset to low precision

text: A high precision dataset was loaded from a file geodatabase into arcsde 9.3.1. The dataload was achieved by copy/paste in ArcCatalog. A custom tool (arcobjects) that extracts from the dataset to a low precision personal geodatabase is now failing when attempting to create the output feature class. It is not feasible to upgrade the personal geodatabase to high precision since all the other datasets it extracts from sde are low precision. Can the hp dataset be downgraded?


title: Transforming geometry coordinates from SRID 4326 to SRID 3011

text: In SQLServer geometry field, I want to change the coordinates of geometry from one spatial reference system to another spatial reference system (from SRID 4326 to SRID 3011) by using ST_Transform. But the destination SRID 3011 is not available in sys.spatial_reference_systems table. How i do this conversion ? Is there any other way ?


title: Call a service web from a ComboBox using JavaScript

text: I want to call a web service and use it in combox for search. Here is the function that i used to call my service web:               var xmlhttp = new XMLHttpRequest();     var url = "http://176.31.230.211/mtmap/MapSrv.svc/getCities";     var params = JSON.stringify({         sessId:"1",         cityName:"ra"     });     xmlhttp.open("POST", url, true);          xmlhttp.setRequestHeader("Content-type", "application/json");          xmlhttp. önreadystatechange = function() {             if(xmlhttp.readyState == 4 && xmlhttp.status == 200) {                     console.log(xmlhttp.responseText);             }     }     xmlhttp.send(params);      I created the ComboBox using the GeoExt bibliography:               var promosCombo = new Ext.form.ComboBox({         id : 'promosCombo',      triggerAction : 'all',      emptyText : "Ville, Avenue, Rue...",      editable : false,      store : promosStore,      mode : 'local',      valueField : 'value',      displayField : 'display',     });      here is my store:               var promosStore = new Ext.data.JsonStore({      url : 'http://176.31.230.211/mtmap/MapSrv.svc/getCities',      fields : ['display', 'value'],      root : 'rows',      autoLoad : true     });      What i want exacty is when i start to tape in my ComboBox the list of the city appear so i can choose the city. It's really urgent. Sorry for my english and Thank you guys for your help and your advice.


title: Adding Distance (KM and Miles) to a Map using ArcGIS 10.1

text: Using ArcGIS 10.1, how do you add and show the distance in mile / km markings alongside the marked trails(KM and miles) while using the aforementioned 10.1? ![enter image description here](http://i.stack.imgur.com/GiSYB.jpg)


title: MapGuide and Google Maps Coordinate System problems

text: I'm trying to do a simple MapGuide OS + Fusion overlaying Google Maps and I can't figure out why the coordinate system are off. I downloaded populated places from natural earth. After Googling, it says that Google Maps uses WGS84 and natural earth is WGS84. Here's what I did: I create a map definition from the natural earth data. Then I created a Flexible Layout with that map definition to view it with fusion. I checked the coordinate system of the map. It says "GCS_WGS_1984" When I display the map, it is incorrect. (The points are way to small and do not correctly align) How do I go about trouble shooting this? I'm brand new to GIS and I think it has to do with the coordinate systems, but as far as I can tell they should match.


title: optimize nearest neighbor query on 70 million extremely high density spatial point cloud on SQL Server 2008

text: I have about 75 million records in a SQL Server 2008 R2 Express database. Each is a lat long corresponding to some value. The table has geography column. I am trying to find one nearest neighbor for a given latitude longitude (point). I already have a query with spatial index in place. But depending on where the record is in the database, say first quarter or last quarter, the query can take about from 3 to 30 seconds to find the nearest neighbor. I feel this can be optimized to give lot faster result by optimizing the query or spatial index. Right now applied some the spatial index with default settings. Here is what my table and query looks like.               CREATE TABLE lidar(         [id] [bigint] IDENTITY(1,1) NOT NULL,         [POINTID] [int] NOT NULL,         [GRID_CODE] [numeric](17, 8) NULL,         [geom] [geography] NULL,      CONSTRAINT [PK_lidar_1] PRIMARY KEY CLUSTERED ([id] ASC)      WITH (PAD_INDEX  = OFF, STATISTICS_NORECOMPUTE  = OFF, IGNORE_DUP_KEY = OFF,       ALLOW_ROW_LOCKS  = ON, ALLOW_PAGE_LOCKS  = ON) ON [PRIMARY]     ) ON [PRIMARY]      The spatial Index i am using:               CREATE SPATIAL INDEX [SPATIAL_lidar] ON [dbo].[lidar] ([geom]) USING  GEOGRAPHY_GRID      WITH (     GRIDS =(LEVEL_1 = MEDIUM,LEVEL_2 = MEDIUM,LEVEL_3 = MEDIUM,LEVEL_4 = MEDIUM),      CELLS_PER_OBJECT = 16, PAD_INDEX  = OFF, SORT_IN_TEMPDB = OFF, DROP_EXISTING = OFF,       ALLOW_ROW_LOCKS  = ON, ALLOW_PAGE_LOCKS  = ON) ON [PRIMARY]      Here is the Query I am using:               declare @ms_at geography = 'POINT (-95.66 30.04)';     select TOP(1) nearPoints.geom.STAsText()as latlon      from     (     select r.geom     from lidar r With(Index(SPATIAL_lidar))     where r.geom.STIntersects(@ms_at.STBuffer(1000)) = 1     ) nearPoints      Here is a sample of lat longs in my database . to give an idea of accuracy and density. All the 70 million records are for one city (Lidar Data)               POINT (-95.669434934023087 30.049513838913736)      Now this query gives me results as i described above, but i want to improve the performance as much as possible. My guess is by tweaking the default values of the spatial index i may be above to better optimize the performance. Any clues this? I tried varying the buffer from 10 to 1000 but with almost same results. Also any other suggestions to improve the performance are welcome. Here is the system i am using right now:               Windows 7 64bit Professional     Intel(R) Core(TM)2 Quad CPU    Q9650  @ 3.00GHz (4 CPUs), ~3.0GHz     Ram: 8 GB     NVIDIA GeForce 9500 GT


title: Using Bounding Box to get area names

text: Question : I have to get the area names of a district. For instance, if my current location is chennai, I have to get the area names of chennai.                I was suggested to use bounding box by google api for iOS, but wanna know is it possible to retrieve all the area names of a district.      Good DAy.


title: Why do most GIS packages need a numeric id?

text: This is a simple yet possibly controversial question: why do most (if not all) GIS packages require that a determined layer have a unique not nullable _numeric_ identifier? Why there is the need for such a surrogate key instead of a natural one? **Examples:**   * ArcGIS enforces OBJECTID (or a GlobalID)   * QGIS does not load layers when they don't have a numeric id.


title: Tilemill (Carto) code for styling multiple features from a single layer

text: How can I more efficiently write this piece of code which aims to style two features in my streets layer (major roads and major road ramps) the exact same way.               #tor_centreline [FCODE=201200] {           line-width: 0.5;           line-color: lightgray;        }          #tor_centreline [FCODE=201201] {           line-width: 0.5;           line-color: lightgray;        }      When I input the code as...               #tor_centreline [FCODE=201200],[FCODE=201201] {      ... I break my map because I believe the second [FCODE=""] is not linked to the #tor_centreline... such that it searches all of my layers for [FCODE=""] and when it finds a layer with no "FCODE" attribute field I receive an error.


title: Split and order lines by polygon intersection

text: I am using PostGIS and QGIS 2.2. I have a shapefile of rivers and a shapefile of state parks and other polygons. I would like to show the user a simplified version of a stream as a straight line, with the edges of the stream line highlighted where it intersects a polygon. Here is a picture of my map in QGIS: ![a normal map of a stream](http://i.imgur.com/sCOLCq5.png) **Please ignore the tributaries** - they are well named and easily culled :) Here is the simplified version I would like: ![A String Map of the stream](http://i.imgur.com/phVehXo.png) In order to create the above image, I would imagine I need to split the stream into several lines. In the above example, I would need six. Blue, Green, Blue, Green, Blue, and Green again. Each of these segments would have a distance and a type (blue, or green). The total distance would be normalized to unit length, and all the segments would be displayed for the user. Another solution would be to simply track the Green segments, each with two offsets: a point, and an end point. In this case, there would only be three segments, and they would show where they started and stopped. This would work well if I wanted multiple types. It was fun to find the distance of the stream inside of polygons vs outside using ST_Intersect and ST_Length. However, that was simply a ratio of inside vs outside. My solution needs both length and order. I think I need to split the line into several smaller lines - both inside and outside the polygons. I think my solution needs to satisfy the following properties:   * Length of each segment must be preserved   * Order of segments must be preserved Thank you for reading.


title: Point Sampling Tool does not preserve set layer order in a qgs project

text: My qgs project contains multiple raster layers in a certain order. It also has a point vector layer at the top position. I used the Plugins -> Analyses -> Point Sampling Tool to extract the raster values at the point positions. This operation created a new point shapefile. I opened the attribute table of the new point shapefile and found the attribute column order in random order, rather than the raster order in the qgs project. Furthermore, I even adjusted the Field order in the panel of Point Sampling Tool to the raster order in qgs project and the attribute table still has a random order. I downloaded the Point Sampling Tool plugin from qgis web. Can I do anything to keep the raster order? I have just started to use QGIS, but I am able to code in C and Python.


title: Create tiles given images and locations; convert markers to tiles

text: **What I need** I need to convert a list of markers to a tile layer I've been googleing without success how to do this. Anyone has faced this problem and has a nice easy solution? **Why I need it** I've developed a HTML application that displays a map. This map has many fixed markers (they do not change in time). Due to that:   1. I want to speed up the application   2. I want to hide the marker position from malicious users (so they can't use the information I've gathered for they're own porpose) I want to convert the markers to a tile layer. Thanks!


title: Error Message Orthorectifying ASTER Imagery in ArcMap 10 with DEM

text: I am trying to orthorectify ASTER imagery (VNIR, 15m resolution) in ArcMap using an ASTER DEM that I have downloaded from the GDEM website. The area is in the Andes, there are very high elevation differences so I guess orthorectification is needed. I was using the tool "Create Ortho Corrected Raster Dataset (Data Management)" and I tried with the DEM as input and also with a derived DHM. I always get an error message very quickly. It's the same coordinate system, I have chosen a very short path (C:/temp), I have tried with file formats *.img and *.tiff. I have never orthorectified before. Any ideas?


title: How to mouse hover on kml layer in arcgis javascript api

text: I am using ArcGIS java script API for map. i want show an popup info window when mouse hover on kml layer.I am trying this code for and using mouse-over event , But this is not working on kml laye. my code is: var kmlUrl = "http://www.arcgis.com/sharing/rest/content/items/a7c2aaaa083741f58825cc3f39c53a22/data";                     var kml = new KMLLayer(kmlUrl);                   map.addLayer(kml);             kml.on("load", function() {               domStyle.set("loading", "display", "none");                  });                  kml.on("mouse-over",function (event) {              alert('test');           });          Please anybody can help me          Thanks


title: Java Exception when enabling WPS for a GeoServer workspace

text: I have installed WPS plugin for GeoServer 2.4.3. I have tried its WPS demo, it works. But when I enable the WPS for a workspace. The following error occurs. Have I reconfigured something?               org.apache.wicket.WicketRuntimeException: Method onFormSubmitted of interface org.apache.wicket.markup.html.form.IFormSubmitListener targeted at component [MarkupContainer [Component id = form]] threw an exception     at org.apache.wicket.RequestListenerInterface.invoke(RequestListenerInterface.java:193)     at org.apache.wicket.request.target.component.listener.ListenerInterfaceRequestTarget.processEvents(ListenerInterfaceRequestTarget.java:73)     at org.apache.wicket.request.AbstractRequestCycleProcessor.processEvents(AbstractRequestCycleProcessor.java:92)     at org.apache.wicket.RequestCycle.processEventsAndRespond(RequestCycle.java:1250)     at org.apache.wicket.RequestCycle.step(RequestCycle.java:1329)     at org.apache.wicket.RequestCycle.steps(RequestCycle.java:1436)     at org.apache.wicket.RequestCycle.request(RequestCycle.java:545)     at org.apache.wicket.protocol.http.WicketFilter.doGet(WicketFilter.java:484)     at org.apache.wicket.protocol.http.WicketServlet.doPost(WicketServlet.java:160)     at javax.servlet.http.HttpServlet.service(HttpServlet.java:727)     at javax.servlet.http.HttpServlet.service(HttpServlet.java:820)     at org.springframework.web.servlet.mvc.ServletWrappingController.handleRequestInternal(ServletWrappingController.java:159)     at org.springframework.web.servlet.mvc.AbstractController.handleRequest(AbstractController.java:153)     at org.springframework.web.servlet.mvc.SimpleControllerHandlerAdapter.handle(SimpleControllerHandlerAdapter.java:48)     at org.springframework.web.servlet.DispatcherServlet.doDispatch(DispatcherServlet.java:923)     at org.springframework.web.servlet.DispatcherServlet.doService(DispatcherServlet.java:852)     at org.springframework.web.servlet.FrameworkServlet.processRequest(FrameworkServlet.java:882)     at org.springframework.web.servlet.FrameworkServlet.doPost(FrameworkServlet.java:789)     at javax.servlet.http.HttpServlet.service(HttpServlet.java:727)     at javax.servlet.http.HttpServlet.service(HttpServlet.java:820)     at org.mortbay.jetty.servlet.ServletHolder.handle(ServletHolder.java:487)     at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1093)     at org.geoserver.filters.ThreadLocalsCleanupFilter.doFilter(ThreadLocalsCleanupFilter.java:27)     at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084)     at org.geoserver.filters.SpringDelegatingFilter$Chain.doFilter(SpringDelegatingFilter.java:74)     at org.geoserver.filters.SpringDelegatingFilter.doFilter(SpringDelegatingFilter.java:45)     at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084)     at org.geoserver.platform.AdvancedDispatchFilter.doFilter(AdvancedDispatchFilter.java:49)     at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084)     at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:311)     at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:68)     at org.springframework.security.web.access.intercept.FilterSecurityInterceptor.invoke(FilterSecurityInterceptor.java:116)     at org.springframework.security.web.access.intercept.FilterSecurityInterceptor.doFilter(FilterSecurityInterceptor.java:83)     at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:72)     at org.geoserver.security.filter.GeoServerCompositeFilter.doFilter(GeoServerCompositeFilter.java:91)     at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323)     at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:68)     at org.springframework.security.web.access.ExceptionTranslationFilter.doFilter(ExceptionTranslationFilter.java:113)     at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:72)     at org.geoserver.security.filter.GeoServerCompositeFilter.doFilter(GeoServerCompositeFilter.java:91)     at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323)     at org.geoserver.security.filter.GeoServerAnonymousAuthenticationFilter.doFilter(GeoServerAnonymousAuthenticationFilter.java:53)     at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323)     at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:68)     at org.springframework.security.web.authentication.AbstractAuthenticationProcessingFilter.doFilter(AbstractAuthenticationProcessingFilter.java:182)     at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:72)     at org.geoserver.security.filter.GeoServerCompositeFilter.doFilter(GeoServerCompositeFilter.java:91)     at org.geoserver.security.filter.GeoServerUserNamePasswordAuthenticationFilter.doFilter(GeoServerUserNamePasswordAuthenticationFilter.java:115)     at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323)     at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:68)     at org.springframework.security.web.authentication.rememberme.RememberMeAuthenticationFilter.doFilter(RememberMeAuthenticationFilter.java:146)     at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:72)     at org.geoserver.security.filter.GeoServerCompositeFilter.doFilter(GeoServerCompositeFilter.java:91)     at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323)     at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:68)     at org.springframework.security.web.context.SecurityContextPersistenceFilter.doFilter(SecurityContextPersistenceFilter.java:87)     at org.geoserver.security.filter.GeoServerSecurityContextPersistenceFilter$1.doFilter(GeoServerSecurityContextPersistenceFilter.java:52)     at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:72)     at org.geoserver.security.filter.GeoServerCompositeFilter.doFilter(GeoServerCompositeFilter.java:91)     at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323)     at org.springframework.security.web.FilterChainProxy.doFilter(FilterChainProxy.java:173)     at org.geoserver.security.GeoServerSecurityFilterChainProxy.doFilter(GeoServerSecurityFilterChainProxy.java:134)     at org.springframework.web.filter.DelegatingFilterProxy.invokeDelegate(DelegatingFilterProxy.java:346)     at org.springframework.web.filter.DelegatingFilterProxy.doFilter(DelegatingFilterProxy.java:259)     at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084)     at org.geoserver.filters.LoggingFilter.doFilter(LoggingFilter.java:75)     at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084)     at org.geoserver.filters.GZIPFilter.doFilter(GZIPFilter.java:42)     at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084)     at org.geoserver.filters.SessionDebugFilter.doFilter(SessionDebugFilter.java:47)     at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084)     at org.geoserver.filters.FlushSafeFilter.doFilter(FlushSafeFilter.java:43)     at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084)     at org.vfny.geoserver.filters.SetCharacterEncodingFilter.doFilter(SetCharacterEncodingFilter.java:109)     at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084)     at org.mortbay.jetty.servlet.ServletHandler.handle(ServletHandler.java:360)     at org.mortbay.jetty.security.SecurityHandler.handle(SecurityHandler.java:216)     at org.mortbay.jetty.servlet.SessionHandler.handle(SessionHandler.java:181)     at org.mortbay.jetty.handler.ContextHandler.handle(ContextHandler.java:726)     at org.mortbay.jetty.webapp.WebAppContext.handle(WebAppContext.java:405)     at org.mortbay.jetty.handler.ContextHandlerCollection.handle(ContextHandlerCollection.java:206)     at org.mortbay.jetty.handler.HandlerCollection.handle(HandlerCollection.java:114)     at org.mortbay.jetty.handler.HandlerWrapper.handle(HandlerWrapper.java:152)     at org.mortbay.jetty.Server.handle(Server.java:324)     at org.mortbay.jetty.HttpConnection.handleRequest(HttpConnection.java:505)     at org.mortbay.jetty.HttpConnection$RequestHandler.content(HttpConnection.java:842)     at org.mortbay.jetty.HttpParser.parseNext(HttpParser.java:648)     at org.mortbay.jetty.HttpParser.parseAvailable(HttpParser.java:211)     at org.mortbay.jetty.HttpConnection.handle(HttpConnection.java:380)     at org.mortbay.io.nio.SelectChannelEndPoint.run(SelectChannelEndPoint.java:395)     at org.mortbay.thread.BoundedThreadPool$PoolThread.run(BoundedThreadPool.java:450)     Caused by: java.lang.reflect.InvocationTargetException     at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method)     at sun.reflect.NativeMethodAccessorImpl.invoke(Unknown Source)     at sun.reflect.DelegatingMethodAccessorImpl.invoke(Unknown Source)     at java.lang.reflect.Method.invoke(Unknown Source)     at org.apache.wicket.RequestListenerInterface.invoke(RequestListenerInterface.java:182)     ... 90 more     Caused by: java.lang.IllegalArgumentException: Parameter message can't be null     at org.apache.wicket.feedback.FeedbackMessage.<init>(FeedbackMessage.java:102)     at org.apache.wicket.feedback.FeedbackMessages.error(FeedbackMessages.java:135)     at org.apache.wicket.Component.error(Component.java:1309)     at org.geoserver.web.data.workspace.WorkspaceEditPage$1.onSubmit(WorkspaceEditPage.java:119)     at org.apache.wicket.markup.html.form.Form.delegateSubmit(Form.java:1565)     at org.apache.wicket.markup.html.form.Form.process(Form.java:958)     at org.apache.wicket.markup.html.form.Form.onFormSubmitted(Form.java:920)     ... 95 more


title: Password does not work on custom install of CartoDB

text: I got an installation of CartoDB working from a recent git pull; however, I cannot login. I am using the domain: development.localhost.lan and username 'development' The password is the one that I gave to create_dev_user. I checked that there is an entry in the database for this user but whatever I do I cannot get it to login. Any suggestions for what to do to get it to work?


title: Create a Dynamic Table in your Map Layout - how to use Graphic table element?

text: Is this Solution (Create a Dynamic Table in your Map Layout using graphic table element) is part of this Extension ESRi production-mapping extension ? how can i use graphic table element ?


title: polylines rotation through the model builder

text: I am searching to add rotation to polylines through the model builder. Is anybody familiar with such option? Thanks Refael