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

text: 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.


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

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


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

text: 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?


title: Heatmap.js is not displaying on OpenLayers

text: 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)


title: How to pansharpen LANDSAT 8 in GRASS?

text: 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...


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

text: 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.


title: Creating columns in PostGIS using osm2pgsql hstore tags

text: 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


title: FeatureLayer SaveEdits failing

text: 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


title: How do I write XML metadata with Python?

text: 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?


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

text: 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:


title: How can I load .SID files in ArcEngine?

text: 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);                      }                 }


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

text: 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.


title: Graphic legend GeoEXT

text: 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']         });


title: repost: how to show context menu in Leaflet

text: 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)


title: ArcGIS - Custom, rule-based borders for polygons?

text: 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.


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

text: 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?


title: Black marks on edge of geotiffs

text: 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.


title: Creating relationship class in ArcSDE?

text: 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)


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

text: 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.


title: GeoServer - Serving html files

text: 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


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

text: 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.)


title: Strip html tags from ArcMap table

text:
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

`


title: QGIS Zonal Statistics just gives 'NULL'

text: 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)


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

text: 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?


title: ArcMap layout graphics

text: 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?


title: Adding Basemap and changing projection - still not aligned

text: 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


title: How to create a radius search in Leaflet

text: 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


title: How can I convert geoJSON data to GPX format

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


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

text: 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 ?


title: How to install QGIS Server Ubuntu 12.04?

text: 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


title: How to test if Osmosis is installed on openSuse?

text: 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?


title: Programmatically modify MapTip Expression

text: 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#?


title: Remove ArcSDE repository and associated files from PostGresql DB

text: 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


title: GPS to map workflow, for non-GIS folk

text: **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.


title: ArcGIS 9.3 Mapbook

text: 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,


title: GDAL Pixel Functions question

text: 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.


title: Kartograph Open Street Maps

text: 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


title: How can i select many rows from a table?

text: 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


title: Uploading Polyline feature in ArcGIS Online

text: 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?


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

text: 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()"


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

text: 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...


title: Create custom Authentication filter in GeoServer 2.3.0

text: ## 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?


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

text: 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!


title: Openlayers Issue seeing Raster Images on both sides of the dateline

text: 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?


title: Adding XY data on QGIS

text: 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?


title: merge rasters with nodata cells

text: 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.


title: Problems with versioning

text: 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?


title: How to convert raster to point in QGIS

text: 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?


title: Convex Quadrilateral Polygons subdivision in equal parts using Python / OGR

text: 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.


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

text: 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!


title: How to cut a SHP based on another SHP using R

text: 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.


title: CRS Austrian Lambert/WGS84

text: 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!


title: Manually add a new point to PostGIS table

text: 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?


title: How to make a group of polygons larger?

text: 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%?


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

text: 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?


title: geotools oracle spatial plugin samples / documentation

text: 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


title: Geoserver SLD style order not applied

text: 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!


title: Hardware configuration needed for geoserver to perform faster

text: 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


title: Increase Map Availability in Geoserver

text: 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


title: How to setup and use GRASS?

text: 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.


title: Qgis: Select multiple layers in sql-Filter

text: 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?


title: Creating a new shapefile based on an existing problematic shapefile

text: 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.


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

text: 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)


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

text: 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.


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

text: 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?


title: Calculating tiles extents using some level data

text: I have some data looking like this:               <?xml version="1.0" encoding="utf-8"?>     <levels>       <level nr="0" tileWidth="512" tileHeight="512" envSizeX="5000" envSizeY="5000" minX="465000" maxX="470000" minY="105000" maxY="110000" scaleFactor="1" tilesCountX="10" tilesCountY="10" />       <level nr="1" tileWidth="512" tileHeight="512" envSizeX="4613" envSizeY="4613" minX="465000" maxX="469999" minY="105001" maxY="110000" scaleFactor="1.0837333333333332" tilesCountX="10" tilesCountY="10" />       <level nr="2" tileWidth="512" tileHeight="512" envSizeX="2286" envSizeY="2286" minX="465000" maxX="469955" minY="105045" maxY="110000" scaleFactor="2.1674666666666664" tilesCountX="5" tilesCountY="5" />       <level nr="3" tileWidth="512" tileHeight="512" envSizeX="1134" envSizeY="1134" minX="465000" maxX="469915" minY="105085" maxY="110000" scaleFactor="4.3349333333333329" tilesCountX="3" tilesCountY="3" />       <level nr="4" tileWidth="512" tileHeight="512" envSizeX="563" envSizeY="563" minX="465000" maxX="469879" minY="105121" maxY="110000" scaleFactor="8.6698666666666657" tilesCountX="2" tilesCountY="2" />       <level nr="5" tileWidth="512" tileHeight="512" envSizeX="275" envSizeY="275" minX="465000" maxX="469777" minY="105223" maxY="110000" scaleFactor="17.339733333333331" tilesCountX="1" tilesCountY="1" />       <level nr="6" tileWidth="512" tileHeight="512" envSizeX="134" envSizeY="134" minX="465000" maxX="469658" minY="105342" maxY="110000" scaleFactor="34.679466666666663" tilesCountX="1" tilesCountY="1" />       <level nr="7" tileWidth="512" tileHeight="512" envSizeX="67" envSizeY="67" minX="465000" maxX="469623" minY="105377" maxY="110000" scaleFactor="69.358933333333326" tilesCountX="1" tilesCountY="1" />       <level nr="8" tileWidth="512" tileHeight="512" envSizeX="33" envSizeY="33" minX="465000" maxX="469554" minY="105446" maxY="110000" scaleFactor="138.71786666666665" tilesCountX="1" tilesCountY="1" />       <level nr="9" tileWidth="512" tileHeight="512" envSizeX="16" envSizeY="16" minX="465000" maxX="469432" minY="105568" maxY="110000" scaleFactor="277.4357333333333" tilesCountX="1" tilesCountY="1" />       <level nr="10" tileWidth="512" tileHeight="512" envSizeX="8" envSizeY="8" minX="465000" maxX="469432" minY="105568" maxY="110000" scaleFactor="554.87146666666661" tilesCountX="1" tilesCountY="1" />     </levels>      I have to use this data to designate tiles' extents (`minX`, `maxX`, `minY`, `maxY`). I wrote a method that gets a level and calculates its tiles. Then it gets the tiles and chooses these which belong to a viewport area. I thought that everything works fine, but I noticed the differences between data (when I compare `GdalInfo` data from metadata of `TIFF` file and my counted values there are some differences... **e.g.**   `GDAL_INFO = 465000, 465554.871466667, 109445.128533333, 110000`   `MY_EXTENTS = (465000, 465499.9, 105001, 105500.9)` Here's my method:               public static IEnumerable<Envelope> GetTiles(string xmlFileName, int level, Rect viewportData, bool debug = false)     {         // viewport data         double x = viewportData.X, y = viewportData.Y;         double width = viewportData.Width, height = viewportData.Height;              // load level data         XDocument doc = XDocument.Load(xmlFileName);         if (doc == null)         {             throw new FileNotFoundException("File '" + xmlFileName + "' not found.");         }              XElement levelData =             doc.Element("levels")                 .Descendants()                 .First(element => int.Parse(element.Attribute("nr").Value) == level);              double minX = double.Parse(levelData.Attribute("minX").Value);         double maxX = double.Parse(levelData.Attribute("maxX").Value);         double minY = double.Parse(levelData.Attribute("minY").Value);         double maxY = double.Parse(levelData.Attribute("maxY").Value);              int tilesCountX = int.Parse(levelData.Attribute("tilesCountX").Value);         int tilesCountY = int.Parse(levelData.Attribute("tilesCountY").Value);              Envelope[][] extents = new Envelope[tilesCountX][];              Envelope bounds = new Envelope();         bounds.MinX = 0;         bounds.MaxX = 0;         bounds.MinY = 0;         bounds.MaxY = 0;         for (int i = 0; i < extents.Length; i++)         {             extents[i] = new Envelope[tilesCountY];             for (int j = 0; j < extents[i].Length; j++)             {                 // count extents for pieces                 double tileGeoWidth = (maxX - minX)/(double) tilesCountX;                 double tileGeoHeight = (maxY - minY)/(double) tilesCountY;                 double offsetX = (double) j*tileGeoWidth;                 double offsetY = (double) i*tileGeoHeight;                      extents[i][j] = new Envelope();                 extents[i][j].MinX = minX + offsetX;                 extents[i][j].MaxX = minX + tileGeoWidth + offsetX; // poprawic dla ostatnich (niewymiarowych) tilesow                 extents[i][j].MinY = minY + offsetY;                 extents[i][j].MaxY = minY + tileGeoHeight + offsetY;                      // check if a piece matches to the viewport                 // 1. find geo-bounds                 if (offsetX <= x)                     bounds.MinX = offsetX;                 if (offsetX <= x + width)                     bounds.MaxX = offsetX;                 if (offsetY <= y)                     bounds.MinY = offsetY;                 if (offsetY <= y + height)                     bounds.MaxY = offsetY;             }         }              // 2. increase "max" bounds (+ width/height)         bounds.MaxX += width;         bounds.MaxY += height;              var intersectedTiles =             extents             .SelectMany(es => es)             .Where(e => EnvIntersects(e, bounds))             .ToList();              if (debug)         {             Console.WriteLine("GLOBAL_EXTENT = {0}, {1}, {2}, {3}", minX, maxX, minY, maxY);             Console.WriteLine("TILES_COUNT = {0}/{1}", tilesCountX, tilesCountY);             for (int i = 0; i < extents.LongLength; i++)             {                 for (int j = 0; j < extents.Length; j++)                 {                     Console.WriteLine("EXTENTS = ({0}, {1}, {2}, {3})", extents[i][j].MinX, extents[i][j].MaxX, extents[i][j].MinY, extents[i][j].MaxY);                 }             }             foreach (Envelope intersectedTile in intersectedTiles)             {                 Console.WriteLine("INTERSECTION = ({0}, {1}, {2}, {3})", intersectedTile.MinX, intersectedTile.MaxX,                                   intersectedTile.MinY, intersectedTile.MaxY);             }         }              return intersectedTiles;     }      And an example call:               Rect viewportData = new Rect(100000, 200000, 70000, 50000);      GdalRetile.GetTiles(@"D:\#PYRAMID\su60ne.xml", 1, viewportData, true);      Could anyone help me to correct this code?   Thanks in advance!


title: Why am I getting an ERROR 000728: Field *fieldname* does not exist within table even though field exists?

text: I do some script for a shapefile's table. And this script works for 2 fields and do not for other ones, because of:               ERROR 000728: Field FID_158_01 does not exist within table      But this field is fine and EXIST! ![enter image description here](http://i.stack.imgur.com/rB4SR.jpg) Script executed in field calculator. I take field names from Iterator (Iterate MultiValue). Expression:               a(!%field_name%!)      Codeblock:               def a(f):       if f==0:         z=1         return z       elif f==(-1):         z=0         return z       elif f==(1):         z=1         return z      Crashlog:               Executing (EXTENT redo math): Model13 H:\TEST\CHECK.shp FID_154_01;FID_153_01;FID_158_01;FID_155_01;FID_153_02     Start Time: Fri Sep 28 14:12:22 2012     Executing (Iterate Multivalue): IterateMultivalue FID_154_01;FID_153_01;FID_158_01;FID_155_01;FID_153_02     Start Time: Fri Sep 28 14:12:23 2012     Succeeded at Fri Sep 28 14:12:23 2012 (Elapsed Time: 0,00 seconds)     Executing (Calculate Field): CalculateField H:\TEST\CHECK.shp FID_154_01 a(!FID_154_01!) PYTHON_9.3 "def a(f):\n  if f==0:\n    z=1\n    return z\n  elif f==(-1):\n    z=0\n    return z\n  elif f==(1):\n    z=1\n    return z"     Start Time: Fri Sep 28 14:12:23 2012     Succeeded at Fri Sep 28 14:12:23 2012 (Elapsed Time: 0,00 seconds)     Executing (Iterate Multivalue): IterateMultivalue FID_154_01;FID_153_01;FID_158_01;FID_155_01;FID_153_02     Start Time: Fri Sep 28 14:12:24 2012     Succeeded at Fri Sep 28 14:12:24 2012 (Elapsed Time: 0,00 seconds)     Executing (Calculate Field): CalculateField H:\TEST\CHECK.shp FID_153_01 a(!FID_153_01!) PYTHON_9.3 "def a(f):\n  if f==0:\n    z=1\n    return z\n  elif f==(-1):\n    z=0\n    return z\n  elif f==(1):\n    z=1\n    return z"     Start Time: Fri Sep 28 14:12:24 2012     Succeeded at Fri Sep 28 14:12:24 2012 (Elapsed Time: 0,00 seconds)     Executing (Iterate Multivalue): IterateMultivalue FID_154_01;FID_153_01;FID_158_01;FID_155_01;FID_153_02     Start Time: Fri Sep 28 14:12:24 2012     Succeeded at Fri Sep 28 14:12:24 2012 (Elapsed Time: 0,00 seconds)     Executing (Calculate Field): CalculateField H:\TEST\CHECK.shp FID_158_01 a(!FID_158_01!) PYTHON_9.3 "def a(f):\n  if f==0:\n    z=1\n    return z\n  elif f==(-1):\n    z=0\n    return z\n  elif f==(1):\n    z=1\n    return z"     Start Time: Fri Sep 28 14:12:25 2012     Failed to execute. Parameters are not valid.     ERROR 000728: Field FID_158_01 does not exist within table     Failed to execute (Calculate Field).     Failed at Fri Sep 28 14:12:25 2012 (Elapsed Time: 0,00 seconds)     Failed to execute (EXTENT redo math).     Failed at Fri Sep 28 14:12:25 2012 (Elapsed Time: 3,00 seconds)      ## Any ideas? When i do all the steps in one model it gives me an error. But when i just make union output file then stop the model and RUN THE SECOND part with field calculations -> error dissapeared and all works fine!!!


title: Benefits of relationship classes and how they can be accessed in c#

text: At the moment I've a personal geodatabase with two imported feature classes from shape files and an own created table. This table has columns for the IDs of the elements in the two feature classes. With ArcCatalog I've created two relationship classes between each sfeature class and the table. Both are 1:1 relationsships. Within ArcMap the user can open the table and get access to more detailed informations by using this relationsships (menu -> Related Tables). Background: One feature class holds features consisting of Polylines, the other features consisting of Points. With the Table I can create a logical relationship between one Point and one Polyline and can add more informations in this table which are only valid for this special relation. Now I've to work with this table in C# and I've to access the other feature classes too. It's no problem to read the IDs out of the table and start a QueryFilter for the feature class where I need the informations from. Sure, this is one possible way but there must exist a more efficient way where I can use the relationship classes for faster operations. Does anybody know how I can query the relationship classes to get all informations at the same time without querying two different tables "manually"? Thanks a lot! Edit: To be more specific. If the table has the two fields ID_EdgeTo and ID_EdgeFrom and one field has a relationship to the feature class EdgesTo and the other to EdgesFrom. When I get through the rows and see a combination of ID_EdgeTo and ID_EdgeFrom I want to query the relationship classes to get more information for the ID's in this specific row. Its more about whether the relationship classes are retrievable and how this has to be done. The table in the database has the form like a "manually created" relationship class with relationship attributes


title: QGIS "Coordinate capture" plugin - how to turn it on?

text: I have enabled the Coordinate Capture plugin via the Plugins menu. But I cannot seem to "turn on" or activate the plugin and get it running. How do I do that?


title: Is it possible to download GPX another coordinate system except WGS(GCS) with GPS Tool on QGIS1.8.0?

text: I got GPS data with using bessel1841 as coordinate system. When I download GPX data on GPS tool, it is automatically treated as WGS84. Can I use another coordinate system when I download GPX on GPS tool. It is impossible to convert data from WGS84 to Bessel1841 on QGIS as far as I know - but I asked here: Is it possible to convert data from WGS84 to Bessel1841?


When we convert data from bessel1841 to WGS84, parameter setting must be used because it is impossible automatic convert on QGIS. Reference: Converting shp file from bessel1841 to WGS84 Is it possible to convert data from WGS84 to Bessel1841? Your advice will be big support for me.


title: rooftops solar potential - solar radiation

text: I'm trying to calculate the solar radiation for rooftops, but i have a big problem, i don't have LIDAR data. I need to give the buildings the information of the surface too (height). I've tried to do these, but for a building I'm going to get for example tree types of heights when I should only get one. So I want to ask if anyone could suggest another way to calculate the buildings height and the height of surface, so that I can calculate the solar radiation for rooftops without using LIDAR data because I can get it. Thank's. Mónica


title: Using way-points to create pathway - MATLAB

text: I'm currently working on a map routing system and am not sure if anyone has done this before. I have collected various way-point coordinates on the map using Google map. How can I create a pathway using the way-points to generate the shortest path?


title: arcgis python attribute error, object has no attribute

text: So I am very new to python and working through my first online course. Does anyone know how I can fix this error:               arcpy.Clip_analysis("RI_Sewers","Providence","C:/student/PythEveryone10_1/PythonInArcGIS/RhodeIsland.gdb/Sewers2","#")      Runtime error       Traceback (most recent call last):      File "<string>", line 1, in <module>      AttributeError: 'module' object has no attribute 'Clip_analysis'      Yes, I have imported `arcpy`. FYI it is preloaded in ARC 10.1 and you don't need to do that anymore. Some functions work in `arcpy`, but not all functions. For example the `arcpy.mapping.ListLayoutElements` function works fine, but clipping and buffers return the error above. I can run the clip and buffer tools fine. Recently I have installed Pyscripter. I have uninstalled and reinstalled that as well. I have used both 64 bit and 32 bit versions. I am currently using the 32 bit version. I have updated my `PYTHONPATH` and tried importing `arcpy.conversion`. I also completed a repair of ARCGIS Desktop software and that did not fix the problem. I also tried running `arcpy` in Pythonwin after importing `arcpy` and got the following error:               File "C:\ESRI\Python27\ArcGIS10.1\Lib\site-packages\pythonwin\pywin\framework\scriptutils.py", line 326, in RunScript      exec codeObject in __main__.__dict__      File "C:\student\PythEveryone10_1\HandleErrors\Debugging.py", line 11, in <module>      arcpy.Buffer_analysis(streams, streams + str(buff), buff)      AttributeError: 'module' object has no attribute 'Buffer_analysis'      I know the scripts I am running are written correctly because they are straight from the esri online course. Thanks


title: Is GIS a tool or a science?

text: Is GIS a tool or a science? There's a lengthy and dense paper on the matter here. Could someone offer a more digestible summary?


title: Using Google Earth as view with SQL Server or Postgres Database

text: How do I connect Google Earth with any relation database? I have **live database** containing **spatial** co-ordinates i.e latitude and longitude. User should be able to view these live data on google earth.


title: How to add Georeferenced PDF as layer to QGIS 2.0?

text: I have some georeferenced pdf's that I would like to open in Qgis. I have read that GDAL version 1.8 and higher will support pdf, and I have GDAL 1.10 installed with Qgis 2.0 on windows 7, But when I add raster layer then choose my pdf file, I get an error that my file is not a supported source.


title: How to get routes and download them from OpenLayers?

text: I am trying to generate a route in openlayers with a start and end points and some points during track. I read about openroute: http://wiki.openstreetmap.org/wiki/OpenRouteService#RouteLink But I don't know how can I include it on my webpage. My idea is about to have a button which when you press it, it generate the route between all markers in map. I want that possibility of download it would be available. In openroute webpage you can introduce points and it generates the route but I don't find an api to see how I can introduce it on my webpage. Do you know other way to generate route in openlayers?


title: Linebreaks in CartoDB info window?

text: I have a "notes" field in my CartoDB table, that I tried to populate with multi-line texts, but there doesn't seem to be a way to add line breaks. If I use `<br />` I just end up with a literal `<br />` in my info window. Can I force a line break?


title: Use QuantumGIS to bulk-import csv files

text: I have a large set of .csv-files containing geographic information. Every dataset consists out of 2 .csv-files! This is my folder-structure. All files in Tracks are formatted equally and all in Routes too, obviously.               ------- Tracks     -- 1.csv     -- 2.csv     -- 3.csv      -- ...     ------- Routes     -- 1.csv     -- 2.csv     -- 3.csv     -- ...      I can import them all for themselves, but that would take days of time (2000 datasets ...). So what I would like to do is to load all those .csv-files at once and group them. So after the import I would like to have a lot of layer groups named 1,2,3 (like the files before the comma) and every group contains the layers representing the two files. Like this:               ---- Group 1     -- tracks1     -- routes1     ---- Group 2     -- tracks2     -- routes 2     ---- Group 3     -- ...      Is that possible in any way? I cannot write python code, so I'm hoping for a plugin ;) Thanks a lot!


title: How to invert polygons (or other options for exterior styling)?

text: I have a shapefile containing many polygons that I want to use to mask a raster. I have tried `gdalwarp -cutline <shapefile> -crop_to_cutline -srcnodata -9999 -dstnodata -9999 <source_raster> <target_raster>`, but the results are visually unappealing because of the low resolution of the raster (i.e. jagged edges). It would look better if I could simply place the vector layer over the raster as a masking layer in the map. I would need to see the underlying raster _inside_ the polygons, but _not outside_. Two options have occurred to me for accomplishing this.   1. Invert the polygons and save as new shapefile then fill with opaque background color.   2. Style the vector layer such that inside the polygons is transparent and outside is opaque. Being new to GIS, I haven't figured out how to accomplish either of these approaches. I have at my disposal the gdal framework, QGIS, and GRASS. I think that I would prefer a GRASS command-line solution, as I may need to repeat this process more than a few thousand times in the next few days. Also, please feel free to direct me to the appropriate documentation that I have obviously missed. I hope that this question isn't a duplicate of one of the many I've read about masking rasters with vectors. It appears to me that this case is different because **I specifically want the final map to display partial pixels from the raster via a smooth cut by a vector**.


title: How to distinguish between File and Personal Geodb workspace?

text: I am looking for a way - given an IWorkspace object - to check if it was created from a Personal-GDB or a File-GDB. I tried using IWorkspace.WorkspaceFactory to check if it's an instance of e.g. AccessWorkspaceFactory but unfortunately this doesn't work for the fgdb. According to .NET the fgdb workspace was created by an AccessWorkspaceFactory, too. Duh. So far I've only come up with the idea that one could check if it's an pgdb by trying to create the workspace using the according factory. The same goes for the fgdb, obviously. Like so:               try {        IWorkspaceFactory factory = new AccessWorkspaceFactoryClass();        factory.OpenFromFile(workspace.PathName);        // if we got that far it seems to be a personal gdb        pgdb = true;     } catch (Exception) {        pgdb = false; // nope, not a personal gdb (or valid path or whatever)     }      But this solution doesn't seem to be very elegant. Are there any data structures to check where the workspace came from?


title: PostgreSQL, find points near points

text: I'm going to have a database with a 'people' table and a 'shops' table. I'd like to find all people within 500m of a shop, and details from both tables. I'm used to working with mySQL, but I'm getting the idea that this would be better in PostgreSQL especially as this needs to scale to a lot of records. I just can't find a suitable example query or helpful documentation. Both tables will be points (and it's where people live, so not a moving target!). The aim of this is to be able to loop through the results (in PHP) and say > The nearest **butcher** shops to **Tom** are: **Browns Butchers, Smiths > Butchers, and Tescos**. > > The nearest **butcher** shops to **Richard** are: **Browns Butchers, and > Tescos**. If anyone knows a mySQL query that wouldn't be too bad, then that's great too. It looks like this cutting down the area first will be most helpful and so I found this which is a bit more than just an sql query. http://www.movable- type.co.uk/scripts/latlong-db.html


title: Can I use the globe plugin on a Windows OS and where can I find the plugin for QGIS?

text: I am trying to add my layers (health facilities and pertinent infrastructure within a region) to a 3D image of the landscape and someone recommended I use the Globe plugin. I can not find the Globe plugin. Does anyone know where I can find it?/ Will the plugin work for a Windows OS?


title: Open Source solutions for building a Traffic Management System

text: I am trying to set an Open Source alternative to proprietry ITS (Intelligent Transportation Systems) The goal is to build a Traffic Management System for a middle size city center. Main requirements are:   * Visualization of road network (especially junctions) on the city map    * Managing traffic lights, cameras and other equipment on roads    * Reporting real time data from sensors, cameras (big data ?)   Also a traffic simulation tool is required. I couldn't find an open source alternative to it because 2D and 3D simulation is needed. Can you advice me the best FOSS4G components to build up a system like this?


title: Table for feature layer in PostGIS database cannot be published and enabled on GeoServer

text: I have tried to publish a table for feature layer in a PostGIS database to GeoServer 2.2.5. The layer was added and published. However, the layer has got a yellow exclamation mark in the "Enabled?" field in Layers view. It is also missing from the Layer Preview page. Everything seemed OK during the publishing process. I could also add the table to QGIS and it's shown well. Anything I was doing wrong? Can anyone give me some hints about it? Thanks! Alex


title: Missing function ST_Extent in PostGIS db

text: I hope you can help as I am tearing my hair out here. Do any of you elite geo peeps have any idea why some specific functions could be missing from my PostGIS install? I am trying to use Mapnik with PostGIS but when I try to run queries from Python it fails with an error of:               ERROR:  function st_extent(geography) does not exist at character 73     HINT:  No function matches the given name and argument types. You might need to add explicit type casts.     STATEMENT:  SELECT ST_XMin(ext),ST_YMin(ext),ST_XMax(ext),ST_YMax(ext) FROM (SELECT ST_Extent(geom) as ext from reseller_map_core_reseller) as tmp      This is the part of the mapnik code:               shapes = mapnik.Layer('Geometry from PostGIS')     shapes.datasource = mapnik.PostGIS(host='localhost', user='user', password='password', dbname='sales_map', table='sales_map_core_sales', geometry_field='geom', srid='4326')     shapes.styles.append('My Style')     m.layers.append(shapes)      I'm just trying to overlay all of the regions from the selected table onto my mapnik generated image. I have installed PostGIS successfully, created a spatial db using extensions and my other GeoDjango/PostGIS seems to work so far. I have looked in the db using pgAdmin and I can see that the st_extents() is missing from the list for some reason even though the others _seem_ to be there. Any ideas why this specific function (and possibly others) could be missing when I have installed everything? I am using:   * PostGIS: 2.0.3   * Postgresql: 9.2.3   * Mapnik 2.10 Thanks for your help!


title: Converting NEAR results to Meters/Kilometers?

text: I am having trouble figuring out how to convert the units from the results of my NEAR calculations. I want to convert them to meters or kilometers, but cannot figure out how to do so... ![\[This is a screenshot of my workspace\]\[1\]:](http://i.stack.imgur.com/Znezj.png) Explanation: The Linear Units for all of my shapefiles and Data Frame are "Meters" so I am not sure why it is not calculating in meters? The real distance (as measured from Google Maps) in the highlighted row should be about 500 meters, but it is only returning 28?


title: Openlayers simplify linestring on the fly

text: I have an issue on my OpenLayers openstreetmaps where I am loading LineString features from a kml file. Everything works fine except performance due to the sheer complexity of the LineString and the number of linestrings. I wanted to use the simplyfy() function to simplify the linestring geometry/feature. Here is some code I wrote to simplify on the fly. The problem is in the line of code below:                           feature.geometry.components[i].simplify(0.1);      This does not seem to modify the original geometry feature components at all. What am I doing wrong? I think we might need to use removeComponents and then add the simplified geometries using addComponents() but how to do this?               preFeatureInsert: function(feature)          {         if (feature != "undefined" && feature.geometry != "undefined" && feature.geometry.CLASS_NAME == "OpenLayers.Geometry.Collection");         {             if (feature.geometry.components != "undefined" && typeof(feature.geometry.components) != "undefined")             {                 for (var i = 0; i < feature.geometry.components.length; i++)                 {                     if (feature.geometry.components[i].CLASS_NAME ==  "OpenLayers.Geometry.LineString")                     {                         feature.geometry.components[i].simplify(0.1);                     }                 }             }                                    }     }


Is there any way of knowing what Travel Modes are available for directions between two places in Google Maps v3? Of course, requesting directions for each travel mode and checking which ones return ZERO_RESULTS is a possibility, but it implies four http requests and four google api requests, which I'd very much like to avoid.