Identification of consecutive points within a specified buffer

I have a point file, hourly relocations of an animal, and I want to be able to place a buffer around each point and calculate the number of subsequent points which are within the buffer zone. I am looking for a method which will work along the point file, like a moving window, which will only count those subsequent points which are within the buffer zone. For instance at point 10 I place a buffer of 1500m and I want to know whether point 11 is within the buffer and if so then whether point 12 is within the buffer and so forth. I don't want to know whether point 52 is within the buffer zone unless all of the previous points have been within the buffer. I also don't want to know whether points 9 or 8 etc. are within the buffer. I have found and tried an additional toolbox called "moving window point analysis" which works as a moving window on the point file. This works well, but very slowly, and includes all points that are within the buffer zone even if they are not consecutive points. I can't find a way to just make it look at consecutive points. I would like a method which will provide an output table as I have a lot of data points to look at in this way. I am using ArcGIS 10. Any help that anyone can provide would be greatly appreciated. Thanks


how to mosaic rasters in grass gis?

I am new user of GRASS GIS. I need to know how can i mosaic different satellite imageries using GRASS GUI? I was trying to use i.image.mosaic, but could not succeed. I have been looking for online tutorial but could not found one. Reference to any online resource will be really helpful. Thank you.


Accessing CRS info with pyqgis in a standalone application

I'm currently having trouble accessing and defining CRS information in a standalone application. I've read the API documentation and have tried it in the python console in QGIS 2.2.0. The python console works fine. I'm able to get the information I need. But the layers in my stand alone application do not have have valid CRS's. What is the reason for this? Code:               import qgis     from qgis.core import *     from qgis.utils import *     from qgis.gui import *          import pdb          QgsApplication.setPrefixPath( r"C:\OSGeo4W64\apps\qgis", True )     QgsApplication.initQgis()     app = QgsApplication([], True)            print "Testing Json crs:"     layer = QgsVectorLayer("E:/GIS_Data/GeoJsonRendering_Test/test.json" ,  'musym', "ogr")     QgsMapLayerRegistry.instance().addMapLayer(layer)     print layer.isValid()     print layer.crs().isValid()     print layer.crs().authid()          print "Testing QgsCoordinateReferenceSystem:"          t = QgsCoordinateTransform()     projSrc = QgsCoordinateReferenceSystem(3857, QgsCoordinateReferenceSystem.EpsgCrsId)     projDest = QgsCoordinateReferenceSystem(4326, QgsCoordinateReferenceSystem.EpsgCrsId)          print projSrc.isValid()     print projDest.isValid()          print "Testing memory layer CRS:"     vTemp = QgsVectorLayer("Point?crs=EPSG:3857", "temporary_points", "memory")     print vTemp.isValid()     print vTemp.crs().isValid()      Output:               Testing Json crs:     True     False          Testing QgsCoordinateReferenceSystem:     False     False     Testing memory layer CRS:     True     False      So I have two valid layers without valid projections. Here is the output for that code in the python console in QGIS:               Testing Json crs:     True     True     EPSG:3857     Testing QgsCoordinateReferenceSystem:     True     True     Testing memory layer CRS:     True     True      So my code only works in the QGIS python console. I'm assuming that I am not defining settings that need to be defined, but I"m not quite sure what I'm missing.


Python QC field calculator

I am new to Python and am trying to run a simple QA/QC on a particular set of fields. I basically want to verify that fields x and y add up to z, I may have to incorporate rounding since a small difference doesn't make a difference and will throw off my results. This is what I have currently and it is not processing:               def Qcontrol(check):       if (!Forest_Acres! + !NonForest_Acres!) = !ACRES!:          return 0       else:         return 1          Qcontrol(!check!)      all fields are double. **Also need to compensate for Nulls in Forest and Non_forest...


Correcting Classification for LIDAR - ArcGIS

Is there any procedure to automatically correct lidar classification? My LAS Dataset has water listed as ground and buildings listed as vegetation etc... Just wondering if there is any option to fix this.


Open a table for editing in a Console outside of ArcMap

I'm trying to write a C# console application to convert a CSV of Lat,Lon data into a Table using insert/update cursors. However, I am having the darndest time figuring out how to get a reference to an ITable without first opening an instance of ArcMap. How would it be done? The purpose is to open a table inside a geodatabase that _exists_. After I have a reference to it, I can create an insert cursor and turn the CSV into a proper table inside the geodatabase directly without bothering with Access (the devil).


Scripting in GRASS

Can anyone suggest me web-links to learn scripting in GRASS GIS (its utilities, examples etc.) other than those available in GRASS official Wiki/Sites?


Statistics field in the dissolve tool not always accurate

When I run the Dissolve Tool using a Statistics Field (SUM being the statistic type) and proceed to check the SUM values for accuracy by comparing against the data in the original feature class, I mostly get accurate results however, for certain records the SUM will not equal the true total of the values in the field I specified. I tried running the tool again but the second time I had just the records which were not calculating properly selected and it yielded desirable results. Has anyone experienced this and if so do you know what causes it and how to remedy the issue. I should note this is processing on point features. For example :               Numbers in original fc  Dissolve field      Numbers in output field OutputField             1013                      A                         1737          A             1013                      A                         899           B             410                       B                      327                       B                      122                       B               The numbers simply come from a pre existing feature class.


Problem reprojecting data from personal geodatabase using QGIS

I have a dataset in the form of a personal geodatabase (mdb) and am interested in using certain feature classes for further processing. My problem is - I am using USGS data orthoimagery which comes in UTM coordinate system and am unable to reproject my geodatabase feature class (layer) to UTM. The same task is easily done in ArcGIS, but I want to use QGIS as there are other features of QGIS that I intend to use later on. What I have tried: 1) Save geodatabase layer (feature class) as shp and try to reproject that shp --> the reproject task does not proceed beyond 2% completion and I end up with a 1KB shp file 2) Reproject mdb layer (Feature class) directly to UTM --> It does not allow me to do this either. 3) Create shapefile from mdb layer (feature class). Define projection for the shp file. Reproject this shp file. --> this also results similar to (1) So, for some reason, projecting layers that have their origin from an ESRI mdb do not seem to be reprojected in QGIS. Any help/tips on this would be welcome. Thanks in advance. ![enter image description here](http://i.stack.imgur.com/uPys4.jpg)


Generate a report or table for buffer zone

I have a polygon feature class. I create a buffer for this feature with a distance being 2 km.Now the buffer output intersect with features of the original polygon.So say for instance from the orginal polygon id 1 it creates a buffer polygon with id 1a. Now I want a report or table of all the features the 1a intersects with the original polygon. Any suggestions?? Thank you.


return url for arcgisoutput

I have a server setup that has port 8080 assigned to ArcGIS Server installation and port 80 to our custom WebADF application. From one of our application, we make a call that generates an image in the output directory. The issue that I have is that the return url doesnt contain 8080 in the url string. It returns this string: http://sfs012/arcgisoutput/_ags_mapa...8cb868696e.jpg I am expecting to get: http://sfs012:8080/arcgisoutput/_ags...8cb868696e.jpg I have tried lot of things but can't seem to figure out how and where to make thing change. Looking for some feedback. Thanks Jay


Collecting street data to populate a graph stucture for routing

I'm implementing a basic digraph with weighted edges for street routing in a small geo area (600mx600m). How and where can I get the required data to populate such a structure? I've been reading about OpenStreetMap, but am unsure as to how I can export the required data into my application. The kind of data I would need is - all the streets (nodes) in the specified area, their connecting streets (subsequent nodes), and the distances between them.


GDAL Brightness RGB Contrast Gamma Saturation

I've not had much experience with modifying imagery, I've recently been playing with GDAL for some batch processing of geo-referenced images, I've been supplied with some imagery in ECW format that is rather dull and could do with some RGB, Contrast Gamma and Saturation adjustments, So finally my question... Is it possible to adjust these values using GDAL or something similar I have managed to find the following values would most benefit my imagery. Brightness -5 R 0 G +5 B -5 Contrast +75 Gamma +2.0 Saturation +190 This was found using irfanview which when used to batch process the imagery seems to strip the georeferenced inforamtion from it :(. Cheers Jamo


How do I run the nvis tool in QGIS GRASS?

I am trying to follow this workflow: [How to generate a DEM from a countour Shapefile?] and it's running fine, and when I hit view output or view results I end up with a flat raster with shading associated with the heights of my contours on screen. However, I want to see the 3D terrain with hillshading etc, and potentially export a 3d mesh, and from research, think that the nvis module in GRASS can do this for me, but when I try to run this I get the error "module nvis not found"? And QGIS crashes and dies. When I go to Manage and Install Plugins, Grass is there, and it has let me perform all the other modules I've needed? So not sure what to do with this?? Any ideas or help much appreciated! Edit: Link to sample data here (.shp+.gml): https://www.dropbox.com/sh/mwxlr3ij78d8gkr/lKwgu3Pafu


Is it possible to store data on Amazon S3 for use in OpenGeo Suite?

We've got about 500GB of geospatial data we'd like to "web-enable" to various degrees and OpenGeo Suite seems like a good route to go down given our circumstances. We can sort out a server (probably VPS) for a good price, but storing all that data locally on the same server will bump up the price considerably. So is it possible to have the VPS running the OpenGeo Suite applications, while the data is stored elsewhere (e.g., on S3)?


Android API - how get the car position (fixed to the road)?

All the android.location classes receives best REAL position of the GPS. Is that any way to get from Google navigation API the CAR position x,y - fixed (tighten, interpolated) to the closest road on the Google map? thanks in advance for any advice


DLL As Embedded Resource ArcMap Add-In? VB.NET

I tried using StackOverflow, but I don't think anyone can answer my question in context(ArcMap Addins). I am currently developing an add-in that needs to load a custom control DLL as a dependency. Unfortunately, ArcMap uses DLL's as their add-ins. So in turn, I must be able to either import an embedded resource to the Visual Studio Toolbox(for use with the forms), or figure out a way to make the control reference the embedded DLL as its delegate. Let me know if you happen to have any ideas. Thanks, Logan


How to create DEM from contours without ArcGIS Spatial Analyst licence?

I'm trying to create a 3d DEM/DTM from a flat contour shapefile which has height data attached to each contour, but I'm not sure how to do this? I've found this post, on creating a dem from contours using the Topo to Raster tool, (using the ANUDEM technology,) which looks great: [see: creating a dem from contours] but I don't have a Spatial Analyst licence, so wondering if there is another way to do this, using a standard ArcView licence (version 9.2)? I'm currently investigating QGIS as an alternative, [via this process: How to generate a DEM from a countour Shapefile?] though I've only just downloaded QGIS, so while I'm a GIS newbie I am SLIGHTLY more familiar with Arc! Any alternatives, help, or tips would be much appreciated!


GDAL_translate: converting ESRI GRID to Geotiff in batch

**The Challenge:** To batch convert 1000 rasters in ESRI Grid binary to Geotiff using either using a batch file or Python 2.7 without ArcGIS/Arcpy module. **The Tools:** Win7Pro, GDAL 1.9.2 (via OSGeo4W), python 2.6/2.7, 1000 binary ESRI Grid rasters (see image) **File structure:** ![File structure](http://i.stack.imgur.com/WgKRF.png) **Single file conversion:** (piece of cake)               gdal_translate GRID1 -of GTiff GRID1.tif      **Batch Route:**               for %i in (*.img) do gdal_translate -of GTiff %i %~SourceFilename.tif      **The Problem:** ESRI Grids don't have a file extension and are instead treated as file directories by the operating system. I need to find an alternative way to iterate. My only option it seems is to iterate by directory as (to my knowledge) there is no method (to my knowledge) to read the GRIDs into a variable like using the ArcPy method: `rasters = arcpy.listRasters("*", "GRID")` Does anyone have any suggestions about how to modify the (*.img) portion of the code so that the loop will read each GRID directory as a file object so I can iterate using the for loop posted above? Note: I'm also aware that I will need to add the parameter for recursion to the batch file and possibly ignore all directories with the name "info". I know I could just use arcpy and python to accomplish this task, but I wasn't satisfied with the results on google of how to batch convert ESRI GRIDs with only GDAL and maybe python. I appreciate any and all advice!


Why won't GeoServer render my SVG symbols?

I'm trying to use GeoServer to render a point layer with custom symbols. This worked perfectly well when my symbols were PNGs, but we haven't decided exactly what size our symbols should be, so I tried to switch to SVGs instead, to take advantage of their cleaner resizing. Now GeoServer won't render anything at all. The relevant bits from the log seem to be:               2013-03-18 23:00:51,655 INFO [geoserver.wms] -      Request: getServiceInfo     2013-03-18 23:00:51,814 WARN [geotools.rendering] - Assuming rendering buffer = 0, but estimation is not accurate, you may want to set a buffer manually     2013-03-18 23:00:51,822 ERROR [lite.StyledShapePainter] - ShapePainter has been asked to paint a null style!!     2013-03-18 23:00:51,826 ERROR [lite.StyledShapePainter] - ShapePainter has been asked to paint a null style!!     2013-03-18 23:00:51,830 ERROR [lite.StyledShapePainter] - ShapePainter has been asked to paint a null style!!     2013-03-18 23:00:51,835 ERROR [lite.StyledShapePainter] - ShapePainter has been asked to paint a null style!!      And so on (I believe once for each feature, but I haven't counted). My SLD is as follows, but it's a near-perfect match for some of the official examples, so I doubt the problem lies there:               <?xml version="1.0" encoding="ISO-8859-1"?>     <StyledLayerDescriptor version="1.0.0" 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"       xsi:schemaLocation="http://www.opengis.net/sld http://schemas.opengis.net/sld/1.0.0/StyledLayerDescriptor.xsd">       <NamedLayer>         <Name>Hydrant</Name>         <UserStyle>           <Title>Hydrant</Title>           <FeatureTypeStyle>             <Rule>               <Title>Red Square</Title>               <PointSymbolizer>                 <Graphic>                   <ExternalGraphic>                     <OnlineResource xlink:type="simple" xlink:href="http://www.shastacsd.org/mapport/img/map-symbols/hydrant.svg" />                     <Format>image/svg+xml</Format>                   </ExternalGraphic>                   <Size>12</Size>                 </Graphic>               </PointSymbolizer>             </Rule>           </FeatureTypeStyle>         </UserStyle>       </NamedLayer>     </StyledLayerDescriptor>      Can anyone think what I'm doing wrong? Thanks in advance if you can!


OpenLayers WFS protocol query after vector layer creation

I want to be able to create a layer on init, then load features from a WFS protocol afterward. I'm having trouble doing this. Using WFS protocol in vector constructor **works fine and is not my problem** : http://dev.openlayers.org/releases/OpenLayers-2.10/examples/wfs-filter.js What i want to do is Add a layer like such:               var layer = new OpenLayers.Layer.Vector("Results", {         projection: new OpenLayers.Projection("EPSG:4326")     });      Then Query it using the WFS Protocol.... this is what i can't get to work... what i have so far is:               var protocol = new OpenLayers.Protocol.WFS({         version: "1.1.0",         url: MapBase.Configuration.GeoServer.WFS_Service.EndPointURL,         featurePrefix: MapBase.Configuration.GeoServer.WFS_Service.featurePrefix,         featureType: MapBase.Configuration.GeoServer.WFS_Service.featureType,         featureNS: MapBase.Configuration.GeoServer.WFS_Service.featureNS,         geometryName: MapBase.Configuration.GeoServer.WFS_Service.geometryName     });               var response = protocol.read({         maxFeatures: Search.Properties.maxResults,         callback: function (resp) {             console.log(resp);             layer.addFeatures(resp.features);         }     });      resp.features gets populated fine. My main problem is I can't find where to specify my filter and projection as i can in the Vector constructor. Even though i specified the projection in the vector constructor, when i addFeatures from my WFS response, they come back in a different projection because they end up really small spooning the west coast of Africa. Same thing with the filter. I don't see that WFS protocol alone supports applying a filter? It only seems to work in conjunction with a vector layer constructor. This seems wrong though because the WFS protocol as described by the OGC support the filter and when openlayers submits it's WFS request via the vector constructor, it does indeed include it just fine. There should be a way to specify projection and filter in the WFS protocol alone? **EDIT #1:** What i have now:               var protocol = new OpenLayers.Protocol.WFS({         version: "1.1.0",         url: MapBase.Configuration.GeoServer.WFS_Service.EndPointURL,         featurePrefix: MapBase.Configuration.GeoServer.WFS_Service.featurePrefix,         featureType: MapBase.Configuration.GeoServer.WFS_Service.featureType,         featureNS: MapBase.Configuration.GeoServer.WFS_Service.featureNS,         geometryName: MapBase.Configuration.GeoServer.WFS_Service.geometryName,         srsName: "EPSG:4326"     //    ,     //    defaultFilter: myFilter     });          protocol.mergeWithDefaultFilter(Search.GetWFSSpatialFilters());


How to scale map appropriately?

I am suppose to make a map. The specifications are 1:50k scale centered on a certain point. They want to view a distance of 30 miles from that centered point, but when I scale to 1:50k the map is only about an 8 mile range. How do I do this in ArcGIS 10?


OpenLayers.Map is not a constructor

I am trying to enter in my maps application (Geoserver, postgis, openlayer) but it is empty and I get an error with Firebug: `TypeError: OpenLayers.Map is not a constructor` in my map.js script in line `var map = new OpenLayers.Map({projection: "EPSG:900913",displayProjection: "EPSG:4326",numZoomLevels: 7});` What is wrong? EDIT: This is my header html               <script type="text/javascript" src="http://extjs.cachefly.net/ext-3.4.0/adapter/ext/ext-base.js"></script>     <script type="text/javascript" src="http://extjs.cachefly.net/ext-3.4.0/ext-all.js"></script>     <link rel="stylesheet" type="text/css" href="http://extjs.cachefly.net/ext-3.4.0/resources/css/ext-all.css" />     <link rel="stylesheet" type="text/css" href="http://extjs.cachefly.net/ext-3.4.0/examples/shared/examples.css" />     <script src="http://openlayers.org/dev/OpenLayers.js"></script>     <script type="text/javascript" src="./script/GeoExt.js"></script>     <script src="http://maps.google.com/maps/api/js?v=3&sensor=false"></script>          <link rel="stylesheet" href="http://openlayers.org/dev/theme/default/style.css" type="text/css">     <link rel="stylesheet" href="http://openlayers.org/dev/theme/default/google.css" type="text/css">          <script type="text/javascript" src="map.js"></script>      The problem occurred suddenly while I looking the layers before added. I clicked in option `Reload feature type` in a particular layer page, Could this have something to do?


SVG symbols for alphabet characters e.g. A, B, C, ... etc

Can anyone tell me whether there is a SVG plugin or source (or how to do) of letters of the alphabet e.g. A, B, C, D etc. I want to create a symbology with colour e.g. gray but also with letters in the fill pattern e.g. SI to denote "Semi-improved" grassland. Ideally trying to create JNCC phase 1 habitat symbols. Thanks Will


How to add image base tool control to a drop-down list in ArcGIS 9.3.1

I have to retrograde an ESRI add-in application from ArcGIS 10.1 to ArcGIS 9.3.1. My challenge is how in the drop down list can I add base tool control( will be shown as icons inside the drop down list : red and green icons )![enter image description here](http://i.stack.imgur.com/mx8FV.jpg), In ArcGIS 9.3.1 is easier to create a drop down list containing text: ![enter image description here](http://i.stack.imgur.com/qlYj7.jpg) By default a drop down list can't contain image , how can I overcome that?


Is it possible to set a tool as optional within ArcGIS Model Builder?

I have a model with two output files - one is useful always and the other only under certain circumstances. The first (always wanted) is a raster dataset, the other (optional) is a vector feature layer derived from the raster. **Is it possible to make the last step (including the corresponding tool) optional** , or do I have to make two separate models, one with and one without the final step? EDIT: I use ArcGIS 10.1. The output feature layer can help, but it's hard to interpret, so I usually hid it just after creation. It was originally intended as a first step of subsequent model and it will be moved to it, but there are some problems in the next model's creation. Having this layer for each set of parameters is quite handy for me to work on the subsequent model, but almost useless for other users. Plus, there's curiosity behind that - quite a lot of tools have optional outputs, so even if I didn't need it I'd be curious whether it's possible in model builder too.


How to use classified styles with joined layers?

I am trying to set a color range based on this tutorial In LAYER PROPERTIES for "new", in JOINS, I have joinlayer - data-join joinfield - GEO_ID2 targetfield - GEOID10 Where "new" is my shapefile and "data-join" is my DBF file with the data I want to use for the graduated color. I want to join column "geo_id2" in "data- join" with "geoID10" in "new." After I am set up, I click "Apply," then "OK." When I toggle over to Style, I see the column I wish to use --- "total-hous" I have selected a color ramp "blues" and 5 classes at equal intervals. But I'm not getting a range like in the Knight example, just five categories all with a range from 0 to 0. Suggestions?


Model builder: How can I use count points in polygon tool in model builder?

I have two layers: one point and one polygon grid layer. Trying to count the number of points in each grid polygon and have it become an attribute in the polygon layer. Currently using ArcGIS 10 with HawthsTools>Analysis tools>Count points in polygon to accomplish this task manually, but I want to bake it into a model. I can't figure out how to add HawthsTools into model builder. Why is this? I have not yet installed Geospatial Modeling Environment, which is the new iteration of HawthsTools. Does GME play nice with model builder? If not, what solutions can you suggest? The closest I've gotten is finding a python script that will count points in polygons. http://forums.arcgis.com/threads/23838-Using-Point-in-Polygon-in- Modelbuilder. Can I use this in a model? Hand holding on how to accomplish this would be very, very appreciated.


Is there any 'Auto-complete polygon' tool in QGIS?

Is there any 'Auto-complete polygon' tool in QGIS, similar to the one existing in ArcGIS?


How to create perpendicular line from point on one line to an intersection with the nearest polygon layer?

I am trying to create laterals from a layer of sewer mainlines and have yet to Google a thoroughly helpful answer. Hoping to find some help here. I have a layer of mainlines that have a direction associated with them pointing from upstream to downstream. I have tabular data that represents the intersection of a lateral line with the mainline. It has the ID of the mainline, the distance from the upstream manhole where it intersects and the side of the mainline that it is on (numbers 0-12, directly downstream being 0/12 or directly upstream being 6, 7-11 being to the left of the downstream direction of the line). I need to create perpendicular lines from these intersection points to the nearest polygon in my parcel layer and extend that lateral line 2 feet into the parcel. I also need to add that parcel number to the lateral in my attribute table and label the laterals ID_1, ID_2, ID_3, etc. on each mainline. So far in ArcMap 10.0 sp4 ArcInfo license, I have created a set of points representing those intersections on the mainline. With Linear Referencing, I created routes out of the mainlines and then made an event layer of lateral lines using the distance from the "zero" measure of the mainline. Now I'm stuck. How do I make perpendicular lines on the correct side of the mainline for these 45k laterals to cross 2 feet into the nearest parcel polygon?


Heatmap.js is not displaying on OpenLayers

I tried to follow the example on the heatmap.js but for some reason it is not displaying the heat map layer, and the layer control is not working. Maybe I am missing something stupid. Can someone point me out what goes wrong?               $(function() {      var testData={     max: 46,     data: [{lat: 33.5363, lng:-117.044, count: 1},{lat: 33.5608, lng:-117.24, count: 1},{lat: 38, lng:-97, count: 1},{lat: 38.9358, lng:-77.1621, count: 1}]     };     var transformedTestData = { max: testData.max , data: [] },     data = testData.data,     datalen = data.length,     nudata = [];          // in order to use the OpenLayers Heatmap Layer we have to transform our data into     // { max: , data: [{lonlat: , count: },...]}     while(datalen--){     nudata.push({     lonlat: new OpenLayers.LonLat(data[datalen].lon, data[datalen].lat),     count: data[datalen].count     });     }     transformedTestData.data = nudata;     var map = new OpenLayers.Map('map_canvas',                     {                         allOverlays: true,                         projection: new OpenLayers.Projection("EPSG:900913"),                         displayProjection: new OpenLayers.Projection("EPSG:900913"),                         units: "m",                         maxResolution: 156543.0339,                         maxExtent: new OpenLayers.Bounds(-20037508, -20037508, 20037508, 20037508),                         controls: [                             new OpenLayers.Control.Navigation(),                             new OpenLayers.Control.PanZoomBar(),                             new OpenLayers.Control.LayerSwitcher(),                             new OpenLayers.Control.ScaleLine(),                             new OpenLayers.Control.MousePosition()                         ]                     });      var layer = new OpenLayers.Layer.OSM();     // create our heatmap layer     var heatmap = new OpenLayers.Layer.Heatmap("Heatmap Layer", map, testData,                     {visible: true, radius: 15},                      {isBaseLayer: false, opacity: 0.3, projection: new OpenLayers.Projection("EPSG:4326")});     //map.addLayer(layer);     map.addLayers([layer,heatmap]);          map.zoomToMaxExtent();     heatmap.setDataSet(transformedTestData);     // Handler for .ready() called.     });      ![heatmap example](http://i.stack.imgur.com/k8ogU.jpg)


How to pansharpen LANDSAT 8 in GRASS?

Having used QGIS for quite some time for mapping, I was 'forced' to dig into GRASS as pansharpening is not directly executable in QGIS. After understanding more or less the structure, I did some trials with i.fusion.brovey. I was able to generate the rasters with 15m resolution, but when loading into QGIS I can not achieve/recuperate the nice colour contrast as is present in the original 30m rasters. In an earlier post someone confirmed that and proposed to use another function for landsat: 'i.his.rgb' (after using i.landsat.toar, which I don't use). He starts to convert min/max to 1-255 before using i.his.rgb...My original rasters have a range of 0-65535. It did not work out for me... I quit after some days of struggling to just perform a pansharp, I feel stuck... Why should just applying i.pansharp not work for Landsat as suggested by others? Can someone suggest a more straightforward workflow to execute pansharpening of Landsat 8 images? Looking forward to some advice, thanks in advance...


Does the Load Data function in ArcCatalog automatically re-project the data?

If I have a feature class within a dataset in one coordinate system and using the Load Data function, load that data into a feature class of a different dataset that's in a different coordinate system ... does the data automatically re-project? I'm thinking it must because it ends up in the correct location. How does this work? I'm a little confused.


Creating columns in PostGIS using osm2pgsql hstore tags

While I feel this may come down to a database question, it feels GIS-related enough that I feel comfortable posting it here. I'll be happy to take it elsewhere if it belongs there, though! I've started the arduous journey of mangling OSM data to work in an RDBMS format that is useful to me. One of the major problems I am running into is the hstore data type. While it is quite useful for being able to hold an extremely large dictionary of values, they are hard to "get at", in my experience. Specifically, the osm2pgsql tool creates `addr:flats`, `addr:housenumber`, and `addr:interpolation` fields, but leaves out some that are quite obvious to me, such as city, postcode, and street. I would like to parse this data to populate new columns with it. What methods would you recommend for accessing such data? Thanks


FeatureLayer SaveEdits failing

I've created a layer which I have bound a DataFeatureForm. Everything is working hunky dory, up to the point where I try to call the layers SaveEdits method in the EditEnded event. When I call this method, the SaveEditsFailed event is raised, returning the following message: > ESRI.ArcGIS.Client.Tasks.ServiceException: Unable to complete operation The Feature's layer API is saying it takes all the edit crud stuff when I hit it with a browser. _also_ I don't know if it will help, but I am providing the Feature Layer a subset of all the layers Out Fields. Could that also effect this? Thanks for yoru help! _EDIT:_ I was able to get the HResult from the exception. HResult = -2146233088


How do I write XML metadata with Python?

Some Python software I have written generates raster files via GDAL. I would like the software to create an XML metadata file for each raster so that I can include details about how the raster was made, who made it, and so-on. I understand that ISO 19115 is the international geospatial metadata standard. I really have no idea where to start with generating a metadata file, though. Can GDAL do this from scratch? Or would it be better to work with a dedicated XML library in Python?


How to Solve the error - "Python Plugin is Broken"?

I tried following the steps to build a python plugin for QGIS. (step 3 is where I am stuck at) 1) I used the Plugin Builder plugin to build the initial set of files. I call my plugin "PointWithinRadius". After using the plugin builder I have the files: **init**.py, pointwithinradius.py, resources.qrc, etc..(the usual files) 2) Next I ran pyrcc4.exe and was able to get resource.py from resource.qrc 3) Tried running pyuic4.bat as follows: (but get the errors as seen below)               ![C:\Program Files\Quantum GIS Copiapo\bin>pyuic4.bat -d -o "C:\Users\stalmaki\.qg     is\python\plugins\PointRadius\ui_pointradius.py" -x "C:\Users\stalmaki\.qgis\pyt     hon\plugins\PointRadius\ui_pointradius.ui"     Traceback (most recent call last):       File "C:\Python25\lib\site-packages\PyQt4\uic\port_v2\invoke.py", line 44, in     invoke         exit_status = driver.invoke()       File "C:\Python25\lib\site-packages\PyQt4\uic\driver.py", line 69, in invoke         self._generate()       File "C:\Python25\lib\site-packages\PyQt4\uic\driver.py", line 103, in _genera     te         self._opts.pyqt3_wrapper, self._opts.from_imports)     AttributeError: Values instance has no attribute 'from_imports'][1]      I am unable to figure out the reason for this and searching for this error didn't give me much success either. Finally when I try opening my plugin in QGIS, I get the following (see attached screenshot jpeg) I want to build a sort of hello world plugin at first but these errors are not allowing me to get out of the starting blocks. Any help on solving would be great. My development environment: Windows 7 32 bit QGIS 1.6.0 Python 2.5 Qt 4.7 Thanks. Sincerely, ST Image attached below:


How can I load .SID files in ArcEngine?

How can I load .SID files in ArcEngine? I think its possible, but I can't find any links. With Kirk's help I got it here's the code:                 public void add_SID_file(string dir)              {                 RasterWorkspaceFactory raster = new RasterWorkspaceFactoryClass();                 IRasterWorkspace rasterWorkspace = raster.OpenFromFile(dir, 0) as IRasterWorkspace;                      if (rasterWorkspace != null)                 {                     IWorkspace w = rasterWorkspace as IWorkspace;                     IEnumDatasetName dsNames = w.get_DatasetNames(esriDatasetType.esriDTRasterDataset);                     IDatasetName datasetName = dsNames.Next();                     List<string> names = new List<string>();                      while (datasetName != null)                     {                         names.Add(datasetName.Name);                         datasetName = dsNames.Next();                     }                     if (names.Count > 0)                     {                         IRasterDataset pRasterDataset = rasterWorkspace.OpenRasterDataset(names[0]);                         IRasterLayer pRasterLayer = new RasterLayerClass();                         pRasterLayer.CreateFromDataset(pRasterDataset);                         mapControl.Map.AddLayer(pRasterLayer);                      }                 }


How to replicate an unversioned SDE geodatabase to a file geodatabase?

I need to implement a process for creating a file geodatabase copy (a read- only snapshot, essentially) of a production SDE geodatabase without using versioning, adding global ID columns, etc. I cannot change any aspect of the SDE geodatabase. All of the documentation I have read on geodatabase replication says that the parent replica data must be versioned. Indeed, attempting to create a simple check-out replica with either the Create Replica wizard in ArcMap or the Create Replica geoprocessing tool fails with the following error message: > Can't create empty replica. Read only data or unversioned data cannot be > replicated. Data versioned with the option to move edits to base cannot be > replicated. GlobalIDs are required for two way and one way replica data. > Creating one way archiving replica requires archiving to be enabled on the > parent. Is there a way around this? I need to be able to run the process once a month so a scripted or programmatic solution would be preferred (arcpy, ArcObjects., etc.). The amount of data to be synchronized is large, probably 500 datasets/tables/feature classes. One solution that comes to mind is exporting and importing ESRI XML Workspace documents, but this seems heavy-handed and error-prone. Another would be a more fine-grained approach of listing and copying individual datasets with arcpy or ArcObjects but if there is an easier, existing solution (i.e. geodatabase replication, if it would let me), that would be preferable.


Graphic legend GeoEXT

UPDATE: The problem is that I got .x-hide-display for Legend graphic. SO I have the picture which is not displayed.. Can somebody help me with the graphic legend request. I have the following code, and I can get the Value with the legend picture, but the legend does not appear in the legend panel for some reasons.. Thank you in advance.                       mapPanel.layers.each(function(record){                  if(!record.getLayer().isBaseLayer && record.id!=highlightedLayer.id){                      var layName = record.getLayer().params['LAYERS'];                 var legendStyle;                           if (legendStyle.indexOf('sld') != -1) {                     var sldURL = encodeURI(config.geoserverURL + "/styles/" + legendStyle);                     record.set("legendURL", config.geoserverURL+"/"+config.geoserverWorkspace + "/wms?REQUEST=GetLegendGraphic&VERSION=1.1.1&FORMAT=image/gif&legend_options=fontSize:12" + "&LAYER=" + layName);                 } else if (legendStyle.indexOf('png') != -1) {                     record.set("legendURL", legendStyle);                 } else {                     record.set("legendURL", config.geoserverURL+"/"+config.geoserverWorkspace + "/wms?REQUEST=GetLegendGraphic&VERSION=1.1.1&FORMAT=image/gif&WIDTH=20&HEIGHT=20&legend_options=fontSize:11&LAYER=" + layName);                 }             }              });//end each-loop                   legendPanel = new GeoExt.LegendPanel({             cls : 'settings',             title : "Map Legend",             itemId : 'legendPanel',             autoScroll : true,             collapsed : true,             border : false,             filter : function(record) {                 return !record.getLayer().isBaseLayer;             },             preferredTypes : ['gx_urllegend']         });


repost: how to show context menu in Leaflet

I asked similar question earlier but didn't get suitable answer. I am now asking again with more clarification. Please see the following google right click menu. I want to create something like this with leaflet `contextmenu` event. I want to get the latlng of that clicking point by contextmenu event and want to do some custom query with this. So the question is simple, How to create right click menu like Google Maps with leaflet API? The context menu will store the latlng event of the clicking position which will use for further query. Thanks ![enter image description here](http://i.stack.imgur.com/9ZEX4.png)


ArcGIS - Custom, rule-based borders for polygons?

I've been working with a map that looks something like this: ![enter image description here](http://i.stack.imgur.com/D5hTN.png) Wherever two polygons of different types (shown above with different colors) meet there needs to be a thick dashed line between them. The solution so far has been to keep a separate feature class for the lines, which of course needs to be updated every time the polygons change. This has gotten to be pretty tedious and is creating a lot of redundancy in our data. Is there any type of map layer that could render the lines on the fly based on spatial relationships and attribute data? I'm guessing a geometric network would be involved - does anyone have any insights? **EDIT** : Sorry everybody, I made a mistake with the screenshot and accidentally had another layer on. The polys are already dissolved to begin with. Please see above for a new version.


Free download of Soviet topographic maps (also known as Genshtab maps)

In the past, http://poehali.org/maps/ offered a free-of-charge download of Soviet military topographic maps of many areas in the world. Those maps are in the public domain. Now, this site redirects to http://mapstor.com , a website offering a low-cost archived download of those maps, but not free-of-charge. http://www.topomapper.com/ allows for interactive looking at Soviet military maps, but only for parts of the world and only up to a scale of 1:100.000, even for areas mapped at 1:50.000. Where can I download free-of-charge public domain Soviet military maps?


Black marks on edge of geotiffs

Just a quick question, are you aware of any way to crop a geotiff while maintaining the co-ordinate infomation. Im trying to crop a series of geotiffs of multi-beam data which have a black frame on one or two edges (in arcmap 10) which i can't remove, even setting the background colour to no colour does't work (it is set so that only the multibeam image is showing though). Im not actually a 100% sure what the marks are, the data was worked in fledermaus to create DTM's and exported as geotiffs then to arc . Any thoughts? I've read the other entries about this problem but i not really confident on how to use gdal or code. I no longer have access to fledermaus so cannot re-edit etc.


Creating relationship class in ArcSDE?

I can not seem to create a relationship class in ArcSDE, or at least one that works. I am succseful in creating the class however ArcMap does not recongnize it (there are no related tables shown) and if I click on relationship classes properties in the context menu in ArcCatalog I get the following error message: ![enter image description here](http://i.stack.imgur.com/kffuS.png) Again the relationship class is created succesfully the results are included. I get the same if I use the catalog tree, toolbox or python script. ![enter image description here](http://i.stack.imgur.com/b4JOa.png)


How to pass in a database created on the fly in code to Quick_Export in arcpy?

So I am trying to programmatically use quick_export through the arcpy package in an arcmap tool, and so I referenced the quick_export tool's dialog box that comes up when it is executing for the parameters. So in my code, prior to using the quick_export function I am creating a database in postgres using some psql commands, and then checking to see if it exists. however, quick_export has failed. So my only thought process on why the quick_export didn't work is because of the latency of recognizing that my database has been created. *When I created the tool I created a field to name the database as a string field.


GeoServer - Serving html files

Encountering a Geoserver - Openlayer WFS layer not visible I found this interesting answer . Solution about html page on the server surprised me very much as I tried previously to load an html page from the jetty geoserver without any success, encountering always a "HTTP ERROR: 404 NOT_FOUND RequestURI=/web/ws3.html" Of course, I tried several other places on my http://localhost:8080/geoserver. Can someone confirm that its possible to load any html page from geoserver, and if a special configuration is required to reach this status ? Thanks very much, Olivier


Free api to reverse geocode latitude, longitude to census tract?

I have a bunch of latitude,longitude coordinates (mostly city/town centroids) in the United States for which I'd like to find the containing census tract. For instance,               input => nn.nnn, -nn.nnn      output => 010300      Is there a free or low-cost API out there that can help? (Essentially Yahoo Placefinder, with the additional response of census tract.) (I understand that many cities include multiple census tracts and I'd expect that some small towns contain only parts of one or more census tracts. I'm comfortable ignoring that issue for now.)


Strip html tags from ArcMap table

I am importing a KMZ we receive as output from an iOS app into ArcMap using the KML to layer tool. The resulting table includes a "pop-up info" column that contians cells with marked up tags like: `

Description	Stakes for monitoring 2 stakes

Time	N/A

Description	N/A

` Is there a simple way to strip these html tags with the field calculator and VB Script / Python within ArcMap, or from the KMZ itself somehow prior to importing? Thanks! would a period at the end of a sentence interrupt the script? That's the only thing I noticed between the stopping point cells... `

Description	View northwest toward point between alignment and high bank. Standing on alignment on left (left) bank of stream channel.

Time	N/A

Photos

`


QGIS Zonal Statistics just gives 'NULL'

Seems like such a simple plugin, but it just won't work. I'm trying to run Zonal Statistics to find the mean value of a raster's pixels within the area of a vector polygon, but all it will give me is 'NULL'. The polygon and the raster are both in the same CRS. There is one null data point within the area, it has a no data value of -9 and is identified as "null (no data). But I don't believe that that is the issue because choosing an area that does not contain a null still doesn't work. In order to see the new columns in the attribute table I find it necessary to toggle editing on the polygon on and off. I've tried running it with editing both on and off, one time it worked with editing on! But never again. Is it supposed to require editing to be on? I'm using QGIS 1.8.0 on a Mac running OS 10.7.5. ![The dashed line marks the bounds of "LowDustRegion". "Gale_Crater_1_tes_seasonal_Daytime_median" is the course resolution raster behind.](http://i.stack.imgur.com/UgIvl.png) ![Only 'NULL' is recorder in the Attribute table.](http://i.stack.imgur.com/Lmg4q.png)


How do I set layer transparency in QGIS 2.0 with Python?

I have a raster layer loaded in QGIS and I'd like to be able to set its transparency via the Python console. I'm new to the API and have not been able to find out how to do this. Any help would be greatly appreciated! _Update_ : Here's what I see from the console:               >>> l = qgis.utils.iface.activeLayer()     >>> l     <qgis.core.QgsRasterLayer object at 0x07144DF8>     >>> rt = l.rasterTransparency()     Traceback (most recent call last):       File "<input>", line 1, in <module>     AttributeError: 'QgsRasterLayer' object has no attribute 'rasterTransparency'      I'm just loading a single band raster in GeoTiff format with a palleted color scheme... nothing too fancy. Am I getting the layer incorrectly?


ArcMap layout graphics

When I try to copy graphics and text between .mxd documents, size of the elements (width, height, font size) change randomly. Anyone know how to transfer LayoutElements between .mxd without any distruptions?


Adding Basemap and changing projection - still not aligned

I am using the basemap in ArcMap10 to create a shapefile of hand digitised polygons based on building outlines. However, when I change the basemap projection to Brit Nat Grid, it does not line up with the OS building polygons that are also Brit Nat Grid. No matter what course of transform/project/tool I use, I cannot get the two to concur


How to create a radius search in Leaflet

My objective is to create a radius search tools using Leaflet API. It will behave like how many POIs within the radius from a given point. Any tutorial or any blog forum link will be helpful to me. Thanks


How can I convert geoJSON data to GPX format

I want to know how I can convert geoJSON data to GPX format.


Get the shortest path between a start point and end point in pgrouting

I want to find the shortest path for a start point and end point chosen by the user. How can I accomplish this in pgrouting. I have the latitude and longitude of the start point as well as for the end point. How can I accomplish this ?


How to install QGIS Server Ubuntu 12.04?

I am new user to Ubuntu (12.04) and am trying to set up QGIS server. I have QGIS desktop installed and working. In regards to QGIS server I have used               apt-get install qgis-mapserver libapache2-mod-fcgid      and was successful but when I try the **GetCapabilities request it fails** and also when I try start/stopping apache2 I get the error **"No Apache MPM package installed"**. Do I need to configure the apache2.conf file for this to work? Thanks for your time. Below is example of error.               ben@ubuntu:~$ aptitude search apache2-mpm-     p   apache2-mpm-event               - Apache HTTP Server - event driven model        p   apache2-mpm-itk                 - multiuser MPM for Apache 2.2                   p   apache2-mpm-prefork             - Apache HTTP Server - traditional non-threa     p   apache2-mpm-worker              - Apache HTTP Server - high speed threaded m     ben@ubuntu:~$ /etc/init.d/apache2 start     No apache MPM package installed


How to test if Osmosis is installed on openSuse?

update; thanks to the answer i have done some more i teste like you recomended;               linux-wyee:/home/martin/gis/test/test2/a_repo # :~> cnf osmosis     linux-wyee:/home/martin/gis/test/test2/a_repo # cnf osmosis          Das Programm 'osmosis' was found here:       * osmosis [ Pfad: /usr/bin/osmosis, Repository: zypp (gis_-_osmosis) ]          Try installing with:         zypper install osmosis      well just wondering - does this mean i have to install it!!? love to hear from your **btw-** do i need postgis and - yes i guess i have to install postgresql too. **Note;** i have a apache / lamp allready up and running on my opensuse 12.1. I am pretty new to GiS and to Linux. I try to install Osmosis on my openSuse 12.1 System. Here what i have in YAST:               linux-wyee:/home/martin # zypper lr -d     #  | Alias                                       | Name                                        | Aktiviert | Aktualisieren | Priorität | Typ    | URI                                                                                | Dienst     ---+---------------------------------------------+---------------------------------------------+-----------+---------------+-----------+--------+---------------------------------------------------------------------------------------+-------      1 | Aktualisierungen-für-openSUSE-12.1-12.1-1.4 | Aktualisierungen für openSUSE 12.1 12.1-1.4 | Ja        | Ja            |   99      | rpm-md | http://download.opensuse.org/update/12.1/                                             |             2 | UpdatedApps                                 | KDE UpdatedApps Repo                        | Ja        | Ja            |   99      | rpm-md | http://download.opensuse.org/repositories/KDE:/UpdatedApps/openSUSE_12.1/             |             3 | gis_-_osmosis                               | gis - osmosis                               | Ja        | Ja            |   99      | rpm-md | http://download.opensuse.org/repositories/Application:/Geo/openSUSE_12.1/             |             4 | openSUSE-12.1-12.1-1.4                      | openSUSE-12.1-12.1-1.4                      | Ja        | Nein          |   99      | yast2  | cd:///?devices=/dev/disk/by-id/ata-TSSTcorp_CDDVDW_SN-S083C_m1L3456789KLMNOP,/dev/sr0 |             5 | packman-x86_64                              | packman-x86_64                              | Ja        | Nein          |   99      | rpm-md | http://packman.inode.at/suse/12.1/                                                    |             6 | repo-debug                                  | openSUSE-12.1-Debug                         | Nein      | Ja            |   99      | NONE   | http://download.opensuse.org/debug/distribution/12.1/repo/oss/                        |             7 | repo-debug-update                           | openSUSE-12.1-Update-Debug                  | Nein      | Ja            |   99      | NONE   | http://download.opensuse.org/debug/update/12.1/                                       |             8 | repo-non-oss                                | openSUSE-12.1-Non-Oss                       | Ja        | Ja            |   99      | yast2  | http://download.opensuse.org/distribution/12.1/repo/non-oss/                          |             9 | repo-oss                                    | openSUSE-12.1-Oss                           | Ja        | Ja            |   99      | yast2  | http://download.opensuse.org/distribution/12.1/repo/oss/                              |            10 | repo-source                                 | openSUSE-12.1-Source                        | Nein      | Ja            |   99      | NONE   | http://download.opensuse.org/source/distribution/12.1/repo/oss/                       |            linux-wyee:/home/martin #       Regarding the above mentioned listing of repos - can we say that I have osmosis installed yet? How can i test if osmosis is installed? Can this be tested on commandline - in terminal? How?


Programmatically modify MapTip Expression

If I load a layer via a shapefile and set               ((ILayer)myLayer).ShowTips = true;       How can I change the tip content (I want to point it to a different column or have a expression of multiple columns) programatically with ArcEngine 10 and C#?


Remove ArcSDE repository and associated files from PostGresql DB

A bit of background first. I have a PostGres DB with PostGIS, I also have ArcSDE installed in this DB and all feature classes created through ArcSDE are using the PG_Geometry. Everything works well! Over the next twelve months I will be looking at getting all of my desktop and web clients away from esri products to open source equivalents which for me means i'll be free of ArcSDE at last (I'm not an esri basher, I enjoy working with esri products but I can't justify the existance of ArcSDE doing something that should be free). **Anyway I want to know if there is a way which I can "uninstall" ArcSDE from PostGres?** I know that if the ArcSDE tables are left there users wont notice a difference but from an admin perspective it would be much nicer and tidier. Has anyone done this in the past? Any tips or experiences? Thanks Ando


GPS to map workflow, for non-GIS folk

**What applications and workflow(s) do you use for gathering, displaying and storing GPS data, with your _centralised_ corporate data, for people who are not GIS folk?** We're in the process of standardizing across all branch offices and I'd like to hear what other folks are doing. While we have GIS users there is a significant and growing number of staff who are not and will never be GIS practitioners but are nevertheless called upon to work with spatial data. Putting Arcmap on the desktop is too expensive in time and training let alone dollars, which also knocks Quantum GIS and friends out of the running. Google Earth is certainly easy to use and the latest versions import GPS data with a minimum of fuss. GE doesn't work for us because their base imagery in our area is so poor it doesn't line up with GPS data by tens of metres and there is no facility to replace their base with our own. Also the GE elevation model is poor and the GPS tracks disappear beneath the ground unless one discards the GPS elevation on import. The free ArcReader solves admirably the problem of fast and easy geo browsing with custom base data -- provided the central agency has ArcPublisher, which we do -- but there is no means of easily dropping GPS data on top. Garmin now has the free Basecamp, but there doesn't seem to be an easy means of getting the GPS data out of Basecamp and into our central data repository (unless I've missed something?). So how about it, what do you do? UPDATE:At the moment the primary use case is for field staff (water license inspectors, conservation officers, wildlife biologists) to collect points and tracks in the field, bring them back to the office, overlay on our in-house base data, pan/zoom/browse/etc., do the odd Identify query, export as bitmap or wmf to include in a report, and save the GPS data somewhere in an ArcGIS accessible format. We, the real GIS folk, pick up the GPS data from the drop box and stuff it into a proper geodatabase for shared corporate use, after the appropriate data massage is complete.


ArcGIS 9.3 Mapbook

I'm looking for a product that will produce multiple pdf plots of mapsheets specified by a shapefile key (say, 11 x 17) from ArcGIS using a single .mxd template, include and orient a common keymap, etc. The DS Mapbook extension from the ESRI developer samples is one of these, but it's an old VB6 product and am wondering if there's anything else similar, newer and maybe commercially available. thanks in advance,


GDAL Pixel Functions question

I am writing a program in python that converts raster files using gdal, and using VRT as an intermediate format. One of the things that i need to do is convert bands between different formats. I am pretty sure that you do this using Pixel Functions. Since python doesn't have this functionality, i have to use C++ for this part. I used some code from the VRT tutorial website, but it doesn't seem to work.               GDALAddDerivedBandPixelFunc("MyFirstFunction", TestFunction);          typedef CPLErr     (*GDALDerivedPixelFunc)(void **papoSources, int nSources, void *pData,                             int nXSize, int nYSize,                             GDALDataType eSrcType, GDALDataType eBufType,                             int nPixelSpace, int nLineSpace);      I implement that function with the code from http://www.gdal.org/gdal_vrttut.html, but it doesn't seem to be working. I am calling gdal_translate to do the conversion process, but it seems like the pixel function does not stay. The addpixelfunc returns a result of 0. Is there some way to write a GDAL extension to have this function be permanent? **Update:** I have been able to correctly create a pixel function and get it to work if i manually convert the file in the same C++ program as the GDALAddDerivedBandPixelFunc resides, however i am still interested in how i could make this permanent. Although it may not be needed since i might want to let the user change some numbers in the function.


Kartograph Open Street Maps

Firstly, congratulations to Gregor Aisch for creating such a wonderful library, thank you! I am having a little trouble rendering an svg from some osm data. I downloaded the Sweden transportation data from here and this is my json config file:               {         "layers": {             "roads": {                 "src": "roads.shp",                 "simplify": false             }         },         "export": {             "width": 800,             "round": 2         },         "bounds": {             "padding": 0.02,             "mode": "bbox",             "data": [377.9582, 59.3031, 378.1233, 59.3731]          }     }      The problem is that it's just creating a blank svg file. When using the polygons mode and if filtered to say a building it renders fine. Any help would be greatly appreciated. Thanks, Tom


How can i select many rows from a table?

I want to create an action that select many rows from PostGIS table to show many lines in QGIS. How can i select many rows from table with python to add these rows in a layer? Thank you


Uploading Polyline feature in ArcGIS Online

I'm seeing some snapping/topological errors in ArcGIS Online (trial version) when I upload my well snapped polyline layer. Please see the screen shot below. ![Polyline Screenshots](http://i.stack.imgur.com/5iCDv.jpg) What may be the cause of this?


Is it possible to show our custom message in QGIS Status bar using python?

Is there any way to show our custom message in Qgis Status bar using python? Just like in arcgis "IApplication.statusbar.message(0) = "Please wait..." like that is there any option to show progressbar in Qgis like "IApplication.progressbar.show()"


Does Identify tool of ArcGIS for Desktop have Python script behind it?

Can anyone let me know the Python script for "Identify" tool in arcmap10. I need to create a customized script tool for that and add in a toolbar.. Pls help...


Create custom Authentication filter in GeoServer 2.3.0

## Context In my current project, I've the requirement to validate that requests coming to GeoServer (2.3.0) are allowed. The project holds those facts:   * the GS client cannot provide the principal information (the password for instance), GS itself has no connection with a user/role repo So we took the opportunity to use the auth filter mechanism in order to check that:   * a valid request (to a specific WFS' layer) contains a special HTTP Header (let's say X-CUSTOM-VALID)   * This header is a JSON encoded message holding enough information to validate the fact that the request was initiated by a client that was connected to a valid third system (a username, a secret, stuffs like that) ## Status The documentation is telling us that we should be able to do so... However, the documentation is not clear how to create such components and how they should be configured. Debugging GeoServer I managed to find that to configure such a filter, it requires a dedicated Authentication Provider. That, in order to have a panel in the web admin interface (under authentications, in the Authentication Filters list) # Panel Thus my code is composed of those files:   * ProducteurAuthFilterPanel.java   * ProducteurAuthFilterPanelInfo.java   * ProducteurAuthenticationFilterConfig.java   * ProducteurAuthenticationFilterPanel.html These are require to add a panel in the Web Admin Interface. `ProducteurAuthFilterPanelInfo` is glueing the two other along with the `ProducteurAuthenticationFilter` here-after (THE **filter** ^^). The `ProducteurAuthenticationFilterConfig` declares that in its constructor:               setClassName(ProducteurAnonymousAuthenticationProvider.class.getName());     setName("producteur");      # Filter (and Provider) Now, the classes needed to create a filter to be included in a chain (I guess):   * **ProducteurAuthenticationFilter** : the filter implementation extending `GeoServerSecurityFilter` and implementing `GeoServerAuthenticationFilter`   * ProducteurAnonymousAuthenticationProvider: somehow required by the Panel (above) to define the new filter   * ProducteurAuthenticationException: used in the AuthenticationEntryPoint (only Http403ForbiddenEntryPoint for now) Finally, the beans are defined like so:               <bean id="yaanonymousFilterProvider" class="dgarne.java.geoserver.security.ProducteurAnonymousAuthenticationProvider"/>          <bean id="producteurAuthPanelInfo" class="dgarne.java.geoserver.security.ProducteurAuthFilterPanelInfo">         <property name="id" value="security.producteurAuthFilter" />         <property name="shortTitleKey" value="ProducteurAuthFilterPanel.short"/>         <property name="titleKey" value="ProducteurAuthFilterPanel.title"/>         <property name="descriptionKey" value="ProducteurAuthFilterPanel.description"/>     </bean>      At the end of the game, in the Web Admin Interface I've a new Item in the filter panel, and I used it in the default mapping (see image below for references): ![enter image description here](http://i.stack.imgur.com/sXWkX.png) # Problem description Here we are... None of my WFS request issued by a client (OpenLayers) which are matching the default mapping (/**) is going through the defined Filter. While debugging I found that the filter chains defined in the Spring Context are never including my definition, but rather are always including the classical one using either anonymous, digest or basic... # Question So is there anyone able to point me out with a (much ^^) more complete documentation about how I have to do it?


Why are KML/KMZ parts missing after upload to Google Maps?

I have a problem with a KMZ file uploaded to my **Google Maps places** : As you can see I have a polygon KMZ file (converted from a shape file in ArcMap) that has 12 parts but the Google Maps shows only 4 of them at once (i.e., on 3 separate views). How can I fix it so they're all visible? Any idea? P.S: The KML is 3,9Mb that's why I made a KMZ that is 0,9Mb. Thanks in advance!


Openlayers Issue seeing Raster Images on both sides of the dateline

I am using OpenLayers V2.13.1 and I am displaying Raster Layers on both sides of the date line When I zoom Out, I cannot see raster Images on both sides of the date line, if I pan across the date line (to the right in this case) the images east of the date line appear and the Rasters to the left disappear. If I pan westward, I can see all images that are east of the date line correctly, but not those west of the date line. Is there any way around this issue as I am displaying Weather and Typhoon tracks that don't care about the date line?


Adding XY data on QGIS

Does QGIS menu have any XY tool like in ArcGIS to add XY data? I have used the other method to capture XY coordinate. Can anyone help?


merge rasters with nodata cells

I want to create single attribute table base on multiple overlapping raster. Some raster do not overlap completely but have common areas. In that situation I need the value from other raster while keeping the non overlapping raster value as missing. I am using arcmap 10.


Problems with versioning

I have 5 PostgreSQL_PostGIS database in my company and i conneting on them through ArcSDE. Four of database work correctly, only with one i have i realy big problem. When i make two versions of database, and i edit data in one version, reconcile and post on DEFAULT version, database do some kinde of automaticly reconcile on second version, and i see all the changes in second version before i do reconcile. Does anyone know what could be the problem?


How to convert raster to point in QGIS

I wanted to convert Raster into Point Vector in QGIS. Is it possible? QGIS have option to convert it into polygon but i didn't found any tool to convert it into point. can anyone help me?


Convex Quadrilateral Polygons subdivision in equal parts using Python / OGR

Consider an irregular polygon (P) which has four vertices. I would like to subdivide P into smaller polygons as shown here http://i.stack.imgur.com/4SKhB.jpg In the old days I was using ArcGIS 9.3 and a VBA macro written by Miles Hitchen to achieve this. The VBA code is following:               Option Explicit          Const xSplit As Long = 10     Const ySplit As Long = 5     Dim dCoords(xSplit, ySplit, 1) As Double               Public Sub GridSelectedQuadrilaterals()     Dim pMxDoc As IMxDocument     Dim pInFtrLyr As IFeatureLayer     Dim pFtrSel As IFeatureSelection     Dim pPolygon As IPolygon     Dim pEnumIDs As IEnumIDs     Dim lID As Long              ' Get the first selected polygon on the first layer         Set pMxDoc = ThisDocument         Set pInFtrLyr = pMxDoc.FocusMap.Layer(0)         Set pFtrSel = pInFtrLyr         Set pEnumIDs = pFtrSel.SelectionSet.IDs         pEnumIDs.Reset              lID = pEnumIDs.Next         While lID >= 0             Set pPolygon = pInFtrLyr.FeatureClass.GetFeature(lID).Shape             GridQuadrilateral pPolygon             lID = pEnumIDs.Next         Wend              MsgBox "Finished"          End Sub               Private Sub GridQuadrilateral(pPolygon As IPolygon)     Dim pSegColl As ISegmentCollection     Dim lIdx As Long     Dim cx(3) As Double, cy(3) As Double     Dim dx(3) As Double, dy(3) As Double     Dim dx2 As Double, dy2 As Double     Dim x1 As Double, x2 As Double, x3 As Double     Dim y1 As Double, y2 As Double, y3 As Double     Dim l As Long, xs As Long, ys As Long              Set pSegColl = pPolygon              ' Get the corner coords of the quad         lIdx = 0         For l = 0 To 3             lIdx = GetIndexOfNextCornerSegment(lIdx, pPolygon)             cx(l) = pSegColl.Segment(lIdx).FromPoint.X             cy(l) = pSegColl.Segment(lIdx).FromPoint.Y         Next l              dx(0) = (cx(1) - cx(0)) / xSplit         dx(1) = (cx(1) - cx(2)) / ySplit         dx(2) = (cx(2) - cx(3)) / xSplit         dx(3) = (cx(0) - cx(3)) / ySplit              dy(0) = (cy(1) - cy(0)) / xSplit         dy(1) = (cy(1) - cy(2)) / ySplit         dy(2) = (cy(2) - cy(3)) / xSplit         dy(3) = (cy(0) - cy(3)) / ySplit              For ys = 0 To ySplit           x1 = cx(3) + dx(3) * ys           y1 = cy(3) + dy(3) * ys           x2 = cx(2) + dx(1) * ys           y2 = cy(2) + dy(1) * ys           dx2 = (x2 - x1) / xSplit           dy2 = (y2 - y1) / xSplit                For xs = 0 To xSplit             x3 = x1 + dx2 * xs             y3 = y1 + dy2 * xs             dCoords(xs, ys, 0) = x3             dCoords(xs, ys, 1) = y3           Next xs              Next ys              BuildGrid          End Sub               Private Function GetIndexOfNextCornerSegment(lStartIdx As Long, pPolygon As IPolygon) As Long     Dim PI As Double     Dim pSegColl As ISegmentCollection     Dim pLine1 As ILine, pLine2 As ILine     Dim l As Long     Dim lNxtIdx As Long     Dim dAng As Double              PI = Atn(1) * 4         Set pSegColl = pPolygon         For l = 0 To pSegColl.SegmentCount - 2             lNxtIdx = lStartIdx + l             If lNxtIdx = pSegColl.SegmentCount Then lNxtIdx = 0             Set pLine1 = pSegColl.Segment(lNxtIdx)             lNxtIdx = lNxtIdx + 1             If lNxtIdx = pSegColl.SegmentCount Then lNxtIdx = 0             Set pLine2 = pSegColl.Segment(lNxtIdx)             dAng = Abs(pLine1.Angle - pLine2.Angle) * 180 / PI             ' *** Ensure we only deal with angles upto 180 ***             If dAng >= 180 Then dAng = 360 - dAng             If dAng > 20 Then                 ' The start point of this segment is a corner point                 GetIndexOfNextCornerSegment = lNxtIdx                 Exit Function             End If         Next l          GetIndexOfNextCornerSegment = -1     End Function               Private Sub BuildGrid()     ' Now create the polygons on 2nd layer     Dim pMxDoc As IMxDocument     Dim pOutFtrLyr As IFeatureLayer     Dim i As Long, j As Long     Dim pFtrCls As IFeatureClass     Dim pFtrCsr As IFeatureCursor     Dim pFtrBfr As IFeatureBuffer     Dim pPtColl As IPointCollection     Dim pPt As IPoint              Set pMxDoc = ThisDocument         Set pOutFtrLyr = pMxDoc.FocusMap.Layer(1)         Set pFtrCls = pOutFtrLyr.FeatureClass         Set pFtrBfr = pFtrCls.CreateFeatureBuffer         Set pFtrCsr = pFtrCls.Insert(True)              For i = 0 To ySplit - 1             For j = 0 To xSplit - 1                      Set pPtColl = New Polygon                 Set pPt = New Point                 pPt.PutCoords dCoords(j, i, 0), dCoords(j, i, 1)                 pPtColl.AddPoint pPt                 pPt.PutCoords dCoords(j, i + 1, 0), dCoords(j, i + 1, 1)                 pPtColl.AddPoint pPt                 pPt.PutCoords dCoords(j + 1, i + 1, 0), dCoords(j + 1, i + 1, 1)                 pPtColl.AddPoint pPt                 pPt.PutCoords dCoords(j + 1, i, 0), dCoords(j + 1, i, 1)                 pPtColl.AddPoint pPt                 pPt.PutCoords dCoords(j, i, 0), dCoords(j, i, 1)                 pPtColl.AddPoint pPt                      Set pFtrBfr.Shape = pPtColl                 pFtrCsr.InsertFeature pFtrBfr                  Next j         Next i              pMxDoc.ActiveView.Refresh          End Sub      The user input at the top of the VBA code tells the script how many tessellations the user requires for each axis               Const xSplit As Long = 10     Const ySplit As Long = 5      Could anyone help me translate this code in Python? I want to achieve the same thing using Python and OGR.


In Quantum GIS, how do I specify the units for a shapefile?

I'm trying to create a shapefile of fields, drawing polygons around each field, and having the acreage calculated in the attribute table. I believe I have the projections correct. It's not coming up with accurate figures. How do I correct this? Thank you!


How to cut a SHP based on another SHP using R

I have a Shapefile of points with certain information. I have another Shapefile with polygons corresponding to the states of the country I am working with (in this case Mexico). Since I only want to work specifically with 3 states, I wanted to know if there is a way to get rid of the points that dont fall into these 3 states.


CRS Austrian Lambert/WGS84

I am totally new to Qgis (or GIS in total): My problem is to get two layers in same CRS to do an intersection and further analysis described in the blog of underdark (Groupstat tutorial). My is the CORINE-AT clip from here: http://www.umweltbundesamt.at/umweltsituation/umweltinfo/opendata/oed_landbedeckung/?cgiproxy_url=http%3A%2F%2Fhttpapp5.umweltbundesamt.at%2Fdata%2Fdataset.jsf%3Bjsessionid%3D2E72E8BB3EB1C39EEBDC124E516B445C%3Fid%3D1 and the country data from DIVA-GIS: http://gadm.org/data/shp/AUT_adm.zip As far as I could see the CORINE data is in Austrian Lambert and the country data is in WGS84. I tried to change the CRS via right click to layer, via vector-> define actual projection, change the WGS84 to Austrian Lambert and then save the shape file... There were always two outcomes: Error in CRS:... or The layer did not lay on each other (country data is flattened). However the result when doing the intersecting (or clipping, or anything else) is an empty attribute table. After two days of trying this task I would be very greatful for a step by step help. (German is highly welcomed) Thanks a lot from desperate!


Manually add a new point to PostGIS table

I did accidently delete a table row - or well, in fact a whole subway station. My attempt to solve the problem was this PostGIS function:               ST_GeomFromText('POINT(coordinates)')      But then PostgreSQL started to complain about some SRID connection with the `spatial_ref_sys` table. What did I do wrong? Is it possible to add this through QGIS instead?


How to make a group of polygons larger?

I need to edit the vertices of about 30 polygons in an equal way, how would I go about this? Alternatively is there some method to just increase the size of a group of polygons 200%?


Preserving origin IDs when creating Origin-Destination (OD) Cost Matrix?

I have created an Origin-Destination cost matrix showing the distance in feet between several thousand home addresses (origins) and about 24 rental vehicles (destinations) at various points throughout the city. Each address point that was input into the OD matrix analysis also originally contained information on the number of times the resident at each address had reserved a particular vehicle. My goal is to create a table that will show, for all home addresses, the distance from that address to each of the 24 rental cars and the number of times the resident reserved each of the 24 cars. I could then say, for example, "40% of reservations made by residents of Neighborhood A are for cars within 1/2 mile" or "Residents who reserved car X lived an average of 2000 feet away". Here is the problem. In conducting the OD matrix analysis, it seems that all the unique identifiers of the data that I input have been obliterated. It would appear that there is no way to associate the distances between origins and destinations back to the other attributes associated with the home address points except for some sort of spatial join. Yet, between the inaccuracies introduced by geocoding thousands of addresses and the fact that many home locations are apartments with multiple residents in the same building, it is impossible to distinguish which resident made which trips by using a spatial join. How have others typically handled the origin-destination problem when origin points overlap and still need to be associated back to other attributes?


geotools oracle spatial plugin samples / documentation

I plan to use geotools and I saw that there is an oracle plugin available, but the documentation is not very detailed. Are there any examples or documentations available to access the functionality of the plugin? Best regards D3


Geoserver SLD style order not applied

I use a .sld to style a road layer according to 4 classes. Normally, the Painter algorithm should be applied: which come first is drawn first. But this does not really apply. The .sld are saved from QGIS. This is the style with the black line on bottom:                   <?xml version="1.0" encoding="UTF-8"?>     <StyledLayerDescriptor xmlns="http://www.opengis.net/sld" xmlns:ogc="http://www.opengis.net/ogc" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" version="1.1.0" xmlns:xlink="http://www.w3.org/1999/xlink" xsi:schemaLocation="http://www.opengis.net/sld http://schemas.opengis.net/sld/1.1.0/StyledLayerDescriptor.xsd" xmlns:se="http://www.opengis.net/se">       <NamedLayer>         <se:Name>Roads</se:Name>         <UserStyle>           <se:Name>Roads</se:Name>           <se:FeatureTypeStyle>             <se:Rule>               <se:Name>M</se:Name>               <ogc:Filter>                 <ogc:PropertyIsEqualTo>                   <ogc:PropertyName>featurecod</ogc:PropertyName>                   <ogc:Literal>1571</ogc:Literal>                 </ogc:PropertyIsEqualTo>               </ogc:Filter>               <se:LineSymbolizer>                 <se:Stroke>                   <se:SvgParameter name="stroke">#bee8ff</se:SvgParameter>                   <se:SvgParameter name="stroke-width">20</se:SvgParameter>                   <se:SvgParameter name="stroke-linejoin">round</se:SvgParameter>                   <se:SvgParameter name="stroke-linecap">round</se:SvgParameter>                 </se:Stroke>               </se:LineSymbolizer>             </se:Rule>             <se:Rule>               <se:Name>A</se:Name>               <ogc:Filter>                 <ogc:PropertyIsEqualTo>                   <ogc:PropertyName>featurecod</ogc:PropertyName>                   <ogc:Literal>1572</ogc:Literal>                 </ogc:PropertyIsEqualTo>               </ogc:Filter>               <se:LineSymbolizer>                 <se:Stroke>                   <se:SvgParameter name="stroke">#d79e9e</se:SvgParameter>                   <se:SvgParameter name="stroke-width">10</se:SvgParameter>                   <se:SvgParameter name="stroke-linejoin">round</se:SvgParameter>                   <se:SvgParameter name="stroke-linecap">round</se:SvgParameter>                 </se:Stroke>               </se:LineSymbolizer>             </se:Rule>             <se:Rule>               <se:Name>B</se:Name>               <ogc:Filter>                 <ogc:PropertyIsEqualTo>                   <ogc:PropertyName>featurecod</ogc:PropertyName>                   <ogc:Literal>1574</ogc:Literal>                 </ogc:PropertyIsEqualTo>               </ogc:Filter>               <se:LineSymbolizer>                 <se:Stroke>                   <se:SvgParameter name="stroke">#ffebaf</se:SvgParameter>                   <se:SvgParameter name="stroke-width">5</se:SvgParameter>                   <se:SvgParameter name="stroke-linejoin">round</se:SvgParameter>                   <se:SvgParameter name="stroke-linecap">round</se:SvgParameter>                 </se:Stroke>               </se:LineSymbolizer>             </se:Rule>             <se:Rule>               <se:Name>streets</se:Name>               <ogc:Filter>                 <ogc:Or>                   <ogc:Or>                     <ogc:Or>                       <ogc:Or>                         <ogc:PropertyIsEqualTo>                           <ogc:PropertyName>featurecod</ogc:PropertyName>                           <ogc:Literal>1575</ogc:Literal>                         </ogc:PropertyIsEqualTo>                         <ogc:PropertyIsEqualTo>                           <ogc:PropertyName>featurecod</ogc:PropertyName>                           <ogc:Literal>1576</ogc:Literal>                         </ogc:PropertyIsEqualTo>                       </ogc:Or>                       <ogc:PropertyIsEqualTo>                         <ogc:PropertyName>featurecod</ogc:PropertyName>                         <ogc:Literal>1577</ogc:Literal>                       </ogc:PropertyIsEqualTo>                     </ogc:Or>                     <ogc:PropertyIsEqualTo>                       <ogc:PropertyName>featurecod</ogc:PropertyName>                       <ogc:Literal>1578</ogc:Literal>                     </ogc:PropertyIsEqualTo>                   </ogc:Or>                   <ogc:PropertyIsEqualTo>                     <ogc:PropertyName>featurecod</ogc:PropertyName>                     <ogc:Literal>1579</ogc:Literal>                   </ogc:PropertyIsEqualTo>                 </ogc:Or>               </ogc:Filter>               <se:LineSymbolizer>                 <se:Stroke>                   <se:SvgParameter name="stroke">#000000</se:SvgParameter>                   <se:SvgParameter name="stroke-width">1</se:SvgParameter>                   <se:SvgParameter name="stroke-linejoin">round</se:SvgParameter>                   <se:SvgParameter name="stroke-linecap">round</se:SvgParameter>                 </se:Stroke>               </se:LineSymbolizer>             </se:Rule>           </se:FeatureTypeStyle>         </UserStyle>       </NamedLayer>     </StyledLayerDescriptor>      and i get this ![enter image description here](http://i.stack.imgur.com/tn3SQ.png) while with the style with the black line on top                   <?xml version="1.0" encoding="UTF-8"?>     <StyledLayerDescriptor xmlns="http://www.opengis.net/sld" xmlns:ogc="http://www.opengis.net/ogc" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" version="1.1.0" xmlns:xlink="http://www.w3.org/1999/xlink" xsi:schemaLocation="http://www.opengis.net/sld http://schemas.opengis.net/sld/1.1.0/StyledLayerDescriptor.xsd" xmlns:se="http://www.opengis.net/se">       <NamedLayer>         <se:Name>Roads</se:Name>         <UserStyle>           <se:Name>Roads</se:Name>           <se:FeatureTypeStyle>                     <se:Rule>               <se:Name>streets</se:Name>               <ogc:Filter>                 <ogc:Or>                   <ogc:Or>                     <ogc:Or>                       <ogc:Or>                         <ogc:PropertyIsEqualTo>                           <ogc:PropertyName>featurecod</ogc:PropertyName>                           <ogc:Literal>1575</ogc:Literal>                         </ogc:PropertyIsEqualTo>                         <ogc:PropertyIsEqualTo>                           <ogc:PropertyName>featurecod</ogc:PropertyName>                           <ogc:Literal>1576</ogc:Literal>                         </ogc:PropertyIsEqualTo>                       </ogc:Or>                       <ogc:PropertyIsEqualTo>                         <ogc:PropertyName>featurecod</ogc:PropertyName>                         <ogc:Literal>1577</ogc:Literal>                       </ogc:PropertyIsEqualTo>                     </ogc:Or>                     <ogc:PropertyIsEqualTo>                       <ogc:PropertyName>featurecod</ogc:PropertyName>                       <ogc:Literal>1578</ogc:Literal>                     </ogc:PropertyIsEqualTo>                   </ogc:Or>                   <ogc:PropertyIsEqualTo>                     <ogc:PropertyName>featurecod</ogc:PropertyName>                     <ogc:Literal>1579</ogc:Literal>                   </ogc:PropertyIsEqualTo>                 </ogc:Or>               </ogc:Filter>               <se:LineSymbolizer>                 <se:Stroke>                   <se:SvgParameter name="stroke">#000000</se:SvgParameter>                   <se:SvgParameter name="stroke-width">1</se:SvgParameter>                   <se:SvgParameter name="stroke-linejoin">round</se:SvgParameter>                   <se:SvgParameter name="stroke-linecap">round</se:SvgParameter>                 </se:Stroke>               </se:LineSymbolizer>             </se:Rule>             <se:Rule>               <se:Name>M</se:Name>               <ogc:Filter>                 <ogc:PropertyIsEqualTo>                   <ogc:PropertyName>featurecod</ogc:PropertyName>                   <ogc:Literal>1571</ogc:Literal>                 </ogc:PropertyIsEqualTo>               </ogc:Filter>               <se:LineSymbolizer>                 <se:Stroke>                   <se:SvgParameter name="stroke">#bee8ff</se:SvgParameter>                   <se:SvgParameter name="stroke-width">20</se:SvgParameter>                   <se:SvgParameter name="stroke-linejoin">round</se:SvgParameter>                   <se:SvgParameter name="stroke-linecap">round</se:SvgParameter>                 </se:Stroke>               </se:LineSymbolizer>             </se:Rule>             <se:Rule>               <se:Name>A</se:Name>               <ogc:Filter>                 <ogc:PropertyIsEqualTo>                   <ogc:PropertyName>featurecod</ogc:PropertyName>                   <ogc:Literal>1572</ogc:Literal>                 </ogc:PropertyIsEqualTo>               </ogc:Filter>               <se:LineSymbolizer>                 <se:Stroke>                   <se:SvgParameter name="stroke">#d79e9e</se:SvgParameter>                   <se:SvgParameter name="stroke-width">10</se:SvgParameter>                   <se:SvgParameter name="stroke-linejoin">round</se:SvgParameter>                   <se:SvgParameter name="stroke-linecap">round</se:SvgParameter>                 </se:Stroke>               </se:LineSymbolizer>             </se:Rule>             <se:Rule>               <se:Name>B</se:Name>               <ogc:Filter>                 <ogc:PropertyIsEqualTo>                   <ogc:PropertyName>featurecod</ogc:PropertyName>                   <ogc:Literal>1574</ogc:Literal>                 </ogc:PropertyIsEqualTo>               </ogc:Filter>               <se:LineSymbolizer>                 <se:Stroke>                   <se:SvgParameter name="stroke">#ffebaf</se:SvgParameter>                   <se:SvgParameter name="stroke-width">5</se:SvgParameter>                   <se:SvgParameter name="stroke-linejoin">round</se:SvgParameter>                   <se:SvgParameter name="stroke-linecap">round</se:SvgParameter>                 </se:Stroke>               </se:LineSymbolizer>             </se:Rule>           </se:FeatureTypeStyle>         </UserStyle>       </NamedLayer>     </StyledLayerDescriptor>      i get this: ![enter image description here](http://i.stack.imgur.com/yq7xz.png) which is the same. Thanks!


Hardware configuration needed for geoserver to perform faster

I Build an application for enabling simple map layer switcher using EXT JS ![enter image description here](http://i.stack.imgur.com/hxrIS.jpg) in this when user press refresh button( which present in the bottom of Layer Tab the application ) read all checked nodes in coma separate value and send as 'dynamic LAYER GROUP' WMS request to the geoserver and served as a single layer to this application How ever the **map rendering is too much slow for 5 user** to increase the map accessing speed I added 3 instance of geoserver in same system with different port with 1024MB java inital and maximum memory pool for each instance and by using Apache HTTP Server i also done Load balancing . where the 'getmap' request is load balanced and transfers to these three geoserver I tuned all the geoserver instance for production environment But still am facing slow rendering speed 1)Is the slow speed is because of using same system for the multiple instance? 2)I am using same database for the three instance did this have anything to do with MAP availability criteria I need to tune my application to work for 30 user simultaneously presently i having an HP server with XEON Processor + 4 GB RAM using stable geoserver 2.4 and openlayer 2.12 what is the best way to tune my application to speed map rendering speed which is the minimum configuration (hardware+software) needed to Speed up map rendering for 30 simultaneous user


Increase Map Availability in Geoserver

I have an application which access same map from two geoserver using loadbalancing concept the map is accessed using WMS service. where the two geoserver connected to a single same POSTGIS database. even the both tomcat is alloted with 3GB of RAM but the map rendering is still too slow Is this have anything to do with the performance of POSTGreSQL?? is there any way to increase map availability ? Please help me


How to setup and use GRASS?

Previously I've asked an question and Micha gave me a helpful advice:How to find the watershed area for each point in the river network using DEM？ After getting the useful DEM file, I began to setup GRASS and use it to solve my problem, I downloaded it from its official site:Index of /grass64/binary/mswindows/native, choosing this one: > WinGRASS-6.4.2-2-Setup.exe to download. I have downloaded the main program and two datasets, and then launch GRASS after fininshing installation: ![enter image description here](http://i.stack.imgur.com/3e51M.png) I create an new folder and put downloaded datasets in it, in my case, > C:\grassdata But as I choose this folder as my GIS data directory, the window soon vanished and nothing came up. I tried many times but the result was just all the same, so I want to setup GRASS correctly, and then use these two useful modules:   * r.watershed elev=DEM thresh= drain=flow_dir_grid   * r.water.outlet drain=flow_dir_grid basin=basin1 easting=X1 northing=Y1 to solve my problem. How to setup GRASS and could someone provide step by step guidance?Thanks.


Qgis: Select multiple layers in sql-Filter

I have a Shape-File with lines of streets. Now I want qgis only to show the categories: Type: primary and Type: secondary and Type: motorway. How do I do this? This doesn't work:               "type" = 'primary' AND "type"'secondary' AND "type"'motorway'      What do I do wrong?


Creating a new shapefile based on an existing problematic shapefile

I have a large floodplain shapefile and would like to create a new shapefile based on coverage. I'm having major problems with the original (intersections and other errors) which I've been unable to fix in either QGIS/Grass or ArcGIS. Is there any way of easily creating a new, error free version of that shapefile without the errors/intersections, preferably in QGIS or Arc? EDIT: I've tried Repair Geometry - Arc didn't find any issues with the file. I've also tried QGIS's geometry check on the same shapefile. It found a number of undefined errors (about 15-20) and around 30-40 intersections. I tried fixing these in grass, but when I export the file, open it up again in QGIS and run the check again I get the same errors.


Color different in Maps exported from ArcMap to PDF and other formats?

When exporting the same map in ArcMap into different formats the color is not consistent. The map is a roadway network and I use grey 30% solid line as the symbol in ArcMap. However, when i export it to PDF and all the other formats, only the PDF resembles the color in ArcMap. All the rest appear like a color of grey 10%. My problem is I have to generate PDF and picture format (which can be directly inserted into Word document) maps at the same time, but the colors appear different. Instead of spending time to make two different map document with adjusted colors, what can I do? The first image below is a jpeg format map, you can hardly see the grey color roadways in the background. The second one is the pdf format map, the grey color shows up exactly as the same as what I see in the ArcGIS. ![enter image description here](http://i.stack.imgur.com/wQsK0.jpg) ![enter image description here](http://i.stack.imgur.com/NTaS1.jpg)


Is there a way to show the direction of a route?

I use pgrouting and openstreetmaps .I would like to know if there is a way to show the direction of the route on the map (eg something like this ---->---- or any other way,idea) Thanks in advance.


IFeature.Store seems to bypass edit sessions and cannot be rolled back. Why?

I'm working with ArcObjects edit sessions for the first time and am having some trouble with it. The C# code I am about to discuss is part of an ArcMap 9.3 extension DLL. My problem is that calls to `IFeature.Store` inside edit operations seem to simply bypass the edit operation (and edit session) in which they execute. This is my code:               IFeatureWorkspace workspace = …;     IWorkspaceEdit workspaceEdit = (IWorkspaceEdit)workspace;     if (!workspaceEdit.IsBeingEdited()) return;          // at this point we know that we are inside an edit session.          IFeatureClass featureClass = workspace.OpenFeatureClass(…);          workspaceEdit.StartEditOperation();     try     {         IFeature feature;         IFeatureCursor features = featureClass.Search(null, false);         // (ESRI recommends search cursors inside ArcMap, even when updating inside edit sessions.)         while ((feature = features.NextFeature()) != null)         {             feature.set_Value(…); // only as an example; assume that the feature is modified somehow.             feature.Store();         }         workspaceEdit.StopEditOperation();     }     catch     {         workspaceEdit.AbortEditOperation();         throw;     }      As you can see, the edit operation in this code modifies all features in some feature class. Let's say that after a few successful feature updates, one update throws an exception. The exception handler in the above code would then call `workspaceEdit.AbordEditOperation();`. As I understand it, this should rollback all modifications done inside the edit operation (including the previous, successful feature updates). I am however observing that the effects of `feature.Store();` never get rolled back. I have also tried this with `IEditor.StartOperation`/`.StopOperation`/`.AbortOperation` instead of the `IWorkspaceEdit` methods, but that doesn't change the observed behaviour. Can anyone explain to me what I am doing wrong? How can I rollback the effects of `IFeature.Store` in an edit session?