title: How do I perform one-step transformation in python?

text: I would like to perform transformation for this example data set.   I have here four pairs of adjustment points and want to transform given point coordinates from primary system to secondary system according to adjustment.               primary_system1 = (3531820.440, 1174966.736, 5162268.086)     primary_system2 = (3531746.800, 1175275.159, 5162241.325)     primary_system3 = (3532510.182, 1174373.785, 5161954.920)     primary_system4 = (3532495.968, 1175507.195, 5161685.049)          secondary_system1 = (6089665.610, 3591595.470, 148.810)     secondary_system2 = (6089633.900, 3591912.090, 143.120)     secondary_system3 = (6089088.170, 3590826.470, 166.350)     secondary_system4 = (6088672.490, 3591914.630, 147.440)          #transform this point     x = 3532412.323      y = 1175511.432     z = 5161677.111<br>         eventually, how do I compute helmert transformation parameters, shifts, rotations and scale factor?      EDIT   at the moment I try to average translation for x, y and z axis using each of the four pairs of points like:                 #x axis     xt1 =  secondary_system1[0] - primary_system1[0]                xt2 =  secondary_system2[0] - primary_system2[0]     xt3 =  secondary_system3[0] - primary_system3[0]     xt4 =  secondary_system4[0] - primary_system4[0]          xt = (xt1+xt2+xt3+xt4)/4    #averaging      ...and so on for y and z axis               #y axis     yt1 =  secondary_system1[1] - primary_system1[1]                yt2 =  secondary_system2[1] - primary_system2[1]     yt3 =  secondary_system3[1] - primary_system3[1]     yt4 =  secondary_system4[1] - primary_system4[1]          yt = (yt1+yt2+yt3+yt4)/4    #averaging          #z axis     zt1 =  secondary_system1[2] - primary_system1[2]                zt2 =  secondary_system2[2] - primary_system2[2]     zt3 =  secondary_system3[2] - primary_system3[2]     zt4 =  secondary_system4[2] - primary_system4[2]          zt = (zt1+zt2+zt3+zt4)/4    #averaging         Is it right way?


title: High quality public transport trip planner examples? Trip visualization examples?

text: Do you know high-quality online public transport route/trip planners? I'm looking for sites with recognized usability, most carrier-owned planners are bad or very bad from Human-Computer Interaction point of view. Also: do you know any papers or works on trip planner GUI design patterns? Example: instead of expandable list with search results we can use the time graph with bars like in Nokia Transit - I just want to see more of such stuff. Thanks! ![nokia maps example](http://i.stack.imgur.com/Ju9s4.jpg)


title: openlayers toggle visibility of layer

text: EDIT: I have cross-posted this to Drupal Answers. I hope this isn't too much against the rules --- I think the answer will be something to do with both Drupal and OpenLayers. http://drupal.stackexchange.com/questions/50733/openlayers-toggle-visibility- of-layer I'm using OpenLayers in Drupal and am stuck while trying to make a layer toggle on and off (as it does from the right-side layer switcher). See this map. http://dev.lambeth.coop/map Click on food growing projects. EDIT: I've changed the code to show the Drupal-specific elements.               (function ($) {       Drupal.behaviors.openlayers_switcher = {         attach: function(context, settings) {           $(document).ready(function() {             $('[id^=openlayers-switcher-block--]').click(function(e) {               $(this).toggleClass('active');               // 27 is the length of 'openlayers-switcher-block--'               // layername is the string following that.               var layername = e.target.id.substring(27);               var ol = $('.openlayers-map').data('openlayers');               var layers = ol.openlayers.getLayersByName(layername);                    // this doesn't work               // if (layers[0].visibility) {                 // layers[0].setVisibility(false);               // }               // else {                 // layers[0].setVisibility(true);               // }             });           });         }       }     }(jQuery));      When I try to toggle visibility - the visibility doesn't change. EDIT: So it appears that the entire click function is being called twice. I've tried putting these lines (in various combinations) at the beginning of the function, as these have been suggested in google searches for 'click events firing twice' etc.               $('[id^=openlayers-switcher-block--]').click(function(e) {       e.preventDefault();       e.stopPropagation();      ... but these don't help - I still get the function called twice. In any case, I'm new to openlayers and would be grateful for a way forward! Thanks!


title: Add table from SDE to ArcMap TOC via ArcMap python window

text: I would like to add a specific table or shape file to my TOC via the ArcMap python window from the SDE connection I'm currently connected to. Since I'm already connected to the SDE connection in ArcMap I'm assuming I don't have to restate the user name and password in the python script right? How would I going about doing this?


title: ArcGIS and SqlServer.Types: Does feature.Shape.SpatialReference.FactoryCode == srid?

text: I created a test implementation of the ESRI Geometry to WKT recipe provided by @SagebrushGIS. It seems to work great, but I'm wondering if the `SRID` value expected by `Microsoft.SqlServer.Types` (specifically `SqlGeometry.STGeomFromWKB`) corresponds directly to `feature.Shape.SpatialReference.FactoryCode`? I stumbled across this thread in the ArcGIS forums where user zroberthu says (not so difinitively) that the two are equal. And when I debug this line (my feature class is projected to NAD_1983_StatePlane_Missouri_Central_FIPS_2402_Feet):               int srid = feature.Shape.SpatialReference.FactoryCode;     // I get srid = 102697      Next, if I go to SpatialReference.org and search by my ESRI projection name, it takes me here, and I see the value `102697` most everywhere I expect to find it. However, if I select the PostGIS insert statement, it shows `srid` is `9102697`, and then an "`auth_srid`" field corresponds to `102697`. So my questions are these: 1) At the end of the day, does the SRID even matter if I'm not performing spatial operations on two different getdatasets? (Presently, I just want a simple way to expose a feature's WKT in the ArcMap UI.) 2) Can I difinitively say an ArcGIS feature's Shape.SpatialReference.FactoryCode corresponds directly to the SRID expected by SqlServer.Types and/or OGC-compliant utilities? I can apply the SRID using other techniques, but if this is true it is handy! Thanks so much. /E


title: How to Partition Line into Segments using ArcGIS for Desktop?

text: I have a shapefile of U.S. states. I extracted the boundary between two states, resulting in a line segment. For example, the boundary between Georgia and North Carolina. Is there a way to further divide this boundary into, say, units of 25 mile segments? Similarly, is there a way to segment the line into equal length segments. For example, taking this border segment, and dividing it into 3 segments of equal length. I use ArcGIS 10.1 and ArcMap. And if we have several lines to split (a few thousands)?


title: How to divide streets layer into 2 mile segments.

text: I have a streets layer which I want to divide into 2 mile segments. I have searched and read enough for days but there is no information on how to do this. The split command allows you to select and divide the poly-line which is manual. It is practically impossible for me to do it for a streets layer which has tens of thousands of poly-lines. Is there any tool, script or combinations of tools I can use to achieve this? I immensely appreciate any help into this problem. Thank you very much in advance.


title: Generate QGIS maps automatic after updating an Access 2010 or Filemaker Pro 12 database

text: I want to develop a database which consist an interaction with QGIS. After updating the data in the database I want to push on a button and start a process in QGIS. After finishing the process in QGIS I want to import the results in the database. Therefor I need or want to build an executable which starts a process in Python. The part in QGIS is already covered, but is it possible to start the processes in QGIS with an external executable? I have found information about starting an external executable from Python, but I want to do it the other way round. Has anyone experience with starting Python processes from an executable? And which executable can be recommended? Thanks in advance.


title: Geocoding Address without Arc Server?

text: I'm helping geocode addresses, and new at it. My office does not have Arc Server, is this required? I found this conversation Address or reverse geocode a point shapefile that points me to address locators. Can I just load the street map and do this in ArcMap without ArcServer? if not, I have been trying to do this through ArcGIS Explorer for Desktop, but it is suprisingly frustrating to migrate from this to ArcMap. It doesn't provide xy locations as a field, and only preserves one field from the original spreadsheet. In ArcMap I was able to do a join, and then 'add xy coordinates'. but there's got to be an easier way.


title: Terrasync COM Port Unavailable

text: I have a Trimble GeoExplorer Series GeoXM that runs Windows PocketPC version 4.2 (Built March 15 2004 for a bit more resolution there). Lately when the unit has been out on a job it will periodically disconnect from the internal GIS receiver (TISP on COM3) and give an error message that says:               "Failed to connect to the GPS receiver because the COM port in in use by another application"      My field techs tell me that this is happening randomly, sometimes several times a day. So far the only way I have to deal with it is to tell them to just keep trying to connect and after a few tries (generally around 20 minutes of waiting) it comes back online and they carry on with their task. This unit is used exclusively to run Terrasync (v3.01) - it doesn't do anything else so I don't understand what other program could be occupying the COM port when Terrasync needs it. I have tried to uninstall and reinstall Terrasync and the problem has persisted. I'd appreciate any tips or tricks to limit this problem - it's really starting to slow my guys down when their out in the field. Thanks!


title: Interpolation Plugin Not outputting image?

text: I am trying to get terrain profiles from shapefile elevation points by creating a TIN. I have been following this very helpful **tutorial**. At the output a new layer is created but with no information/display within the map itself, but the layer's histogram shows that the data is indeed there. Has anyone encountered this before? Any suggestions are helpful.


title: arcpy clip_managment migration from arcpy 10.1 to 10.2

text: i've created script using arpy with arcgis 10.1. One of the goal is to clip raster via shp file, which contains one feature. it worked perfect.   after migration to arcgis 10.2 i've noticed some problem. result of clipping raster in exact shape of feature. it is just clip to extent of clipping feature.               arcpy.Clip_management(inFC, "#", cut_director + "/" + outFC, shp + "/" + shp_clip, "#", "ClippingGeometry")      i've tried to do clipping via arctoolbox (it does the job) and then copy the process as python snippet but after changing input, output and clipping feature i received error. thanks in advance for help edit: i've found answer. when change nodata parametr from "#" (second one in code above) to "256" it works.


title: Which spatial bundling or grouping algorithm optimally groups elements by attribute?

text: As I delve deeper into GIS, I find myself trying to improve efficiencies of certain procedures and processes. The one process I'm working on now is spatially scheduling construction projects. For example, consider the two construction schedules outlined in the images below. The first one is perhaps the defacto non spatially planned schedule; however, the second, is more optimal and due to bundling of contracts would be cheaper as well. This is easy when the network has 30 links; however, the one I am working on has ~50,000 links. **Question: What should I read upon or learn to implement such algorithms? Further more, are there any application type paper that you can reference to help my understanding on the matter?** ![Non bundled Roads](http://i.stack.imgur.com/Tk4i3.jpg) # Which GIS algorithm can I use to transfer the above into the bottom ![Bundles roads](http://i.stack.imgur.com/ja7ei.jpg) **Note:** I'm not looking for a software tool to do the work for me, I'm looking to understand the process and perhaps implement it as a plugin in QGIS or ArcGIS


title: What is an Orphan point?

text: MY apologizes about the next doubt, but I can't understand what is an **orphan point**. I am using QGIS with a module to apply the Linear Reference System (the module is called "LRS"). During the Calibration step, I select the Linear layer and Point layer. After I go to "Errors" step and the module identifies the next errors: > Not enough points Orphan point The question is: What is an orphan point? **Solved** : The problem was solved. The id's of each layer must have the same ID.


title: How to get the current map from a geoprocessing tool execute function?

text: I'm making a geoprocessing tool with the ArcObjects SDK Java 10.2. I need to modify labeling colors of a layer of the current map. Is it possible from the execute function to access the map currently opened in ArcMap ? I need to get the current IMapDocument or IMap in order to set up parameters of the layers (IGeoFeatureLayer). How can I do that please ?


title: Looking for (nearly) worldwide elevation data

text: I know the raw ASTER and SRTM datasets. Is there any error corrected version available? Or did anybody already combine the data? At least the ASTER data often misses tiles for costal regions (if there is only few elevation and the file would be nearly empty). (I don't mind if I have to pay a fair price for it and it could even be just a web service where I get elevation data for a certain coordinate) Jan


title: How to add Routing Table Query-Result on pgrouting?

text: I've followed the pgrouting latest workshop http://workshop.pgrouting.org. I wanted to know how to add the routing table query-result into my webpage directly after the route showed up? regards


title: Can I create a Model in MapInfo, like in ArcGIS? (Modelbuilder)

text: I have some Forumlas I use quiet often. I want to create a Model like in ArcGIS with it(Drag&Drop put in some Parameters, and so forth), is that possible? Is there a GUI for that? Or is it only possible to write MapBasic code?


title: Finding specific pixel values in ArcGIS using Python

text: I have a raster file containing some integer values. I need to find all pixels whose values are not defined in list like this one: `l = [1,2,4,8,16,32,64,128]`. Is there any easy way to solve that problem without converting raster to points and editing attribute table?


title: How to get a bbox selection of huge shapefile as GeoJSON?

text: I have a 4.29 GB Shapefile from which I want to extract features as GeoJSON, defined by a bounding box. This should be as fast as possible, because I want to use it in a web application. My first attempt at this was to convert to SQL using shp2psql and import it into PostGIS. I created a spatial index (gist) and did a full ANALYZE VACUUM on the table. The VM has 4 GB of memory. The following query works, but took 623 seconds(!):               SELECT         ST_AsGeoJSON(geom)AS geometry,     FROM         shapes     WHERE         geom && ST_GeogFromText(             'POLYGON((5.117705762386322 52.09243663877818,5.118895322084427 52.09243663877818,5.118895322084427 52.09314276719956,5.117705762386322 52.09314276719956,5.117705762386322 52.09243663877818))'         )     AND ST_Intersects(         geom,         ST_GeogFromText(             'POLYGON((5.117705762386322 52.09243663877818,5.118895322084427 52.09243663877818,5.118895322084427 52.09314276719956,5.117705762386322 52.09314276719956,5.117705762386322 52.09243663877818))'         )     )      An equivalent using ST_DWithin took even longer. I decided to try an Amazon RDS PostgreSQL 9.3.1 instance with 68 GB of memory and 26 CPU's. Query time dropped to around 30 seconds. Faster but still not fast enough. Preprocessing will take forever, even on Amazon RDS. And it will cost a lot of money. I also tried Fiona and 'ogr2ogr -spat', but they're about as slow as PostGIS. What other options exist for quickly getting GeoJSON out of a large Shapefile for a certain spatial extent? Thanks!!


title: ArcGIS Server 10.0 too slow to generate cache tiles

text: I want to define the zoom level scale for my arcgis web map, but it seems that I should build the cache tiles first (in ArcMap). I have waited for 2 hours but the cache status is still maintain in 0%. I was wondering if this is normal? Thanks. Details:   * All files in TIFF (elevation map): (28,4 GB)   * ArcGIS Server 10.0 Enterprise Windows 7 32 bit.   * In _Analyze Map_ 10 Warnings and 7 Messages   * 7 scales How to speed it up?


title: Why can I see the data from a second sde database connection in my database?

text: I have an sde 9.3 database which when opened in arccatalog shows the data from another sde database. Under settings > Source this is confirmed. Why and how can this happen? The only thing the 2 databases have in common is that the owner is the same.


title: ArcGIS 10 Rest services regulary stop working

text: I have a problem with a new installation of ArcGIS Server. The installation is ArcGIS 10 on Windows Server 2008r2 with IIS 7. I have followed this guide to configure the IIS ArcGISServicesAppPool: http://support.esri.com/en/knowledgebase/techarticles/detail/36887 I did this because my cache wouldn't work without it. The problem I have is that regularly the REST services stop working. The url to /arcgis/rest/services still lists all my services but when I want to open one, I get an error that the service is not available. Also when I want to go to the ArcGIS/rest/admin/ I cannot log in. I can solve this problem by recycling the ArcGISServicesAppPool in IIS. Then I can log in to the rest admin. But then I still need to clear the ArcGIS rest cache before my services work again. One thing I've changed was to set the Idle time-out in IIS for the ArcGISServicesAppPool to 0, because I noticed that every time after that time out, my rest services didn't work anymore (same problem). So after this change, the services will run much longer, but still after a day or so, I need to recycle and clear the cache. I've done some more checks and I notice that every time the ArcGISServiceAppPool automaticly recycles (via settings in IIS) the services stop being available. The solution again is to manually recycle the application pool and then clear the rest cache again. I've solved my problem by adding ArcGISWebservices user to the Administrators group. Don't find it a clean solution, but it's a solution. Can there be other problems with this? Security?


title: Adding fields to raster attribute table

text: I have a raster attribute table associated with an IMAGINE .img raster that I build and open with:               IGPUtilities utils = new GPUtilities();     IRasterDatasetEdit2 raster = (IRasterDatasetEdit2)         utils->OpenRasterDatasetFromString(path);          raster.BuildAttributeTable();     ITable vat = (raster as IRasterBandCollection).Item(0).AttributeTable;      I need to add new column to the table and populate with values. However, calling `AddField` on the table throws " _The method or operation is not implemented_ " exception.               IFieldEdit field = new FieldClass();     field.Type_2 = esriFieldType::esriFieldTypeString;     field.Name_2 = "Name";     vat.AddField(field as IField);  // Throws      Is there a direct way how to modify the existing attribute table or do I have to create one on my own and assign it to the raster using `IRasterDatasetEdit2::AlterAttributeTable()`? I'm using ArcObjects .NET SDK with ArcGIS 9.3.1. **Edited** Once again ArcMap behaves curiously: When I open the raster in ArcMap after my code failed to add the field, I can open the attribute table but the 'Add Field...' option is disabled. If I remove the raster and open it again, I can add fields from ArcMap and also my code works! Everything works fine with GeoTiff images. I do not think this is a locking issue since I'm releasing all COM instances (via `System::Runtime::InteropServices::Marshal::ReleaseComObject()`)


title: Is it possible to find if a value is present in Array without looping in mapbasic?

text: How to find a value is present in array without looping ?               e.g.     For i=1 TO ubount(Array)       Find if value is present     Next      or some more direct statement like               Array.contains(somevalue) like in java


title: How to add element to Mapbasic Array?

text: Code:                Dim name_array() As String      Redim name_array(10)      How to add value to the array? I know we can add value by specifying index like name_array(i) = "some_value" But I don't want to specify any index _i_. Instead is it possible to Add a value to array to last index used. eg. _If i have added to i=5 Will I be able to add to 6th index without specifying index?_


title: When i'm in editing mode in QGIS, how do I make the selection color transparent?

text: ...because i need to split a polygon and i want to see the layer underneath it. But I have to select the polygon in order to cut it, and when i do that it turns yellow!


title: Viewing current users of map (mxd)?

text: Is anyone aware of a script or tool that would allow me to view what users have a paticular mxd open at one time? We are using ArcGIS Desktop 9.3.1.


title: How to open XML file in ArcGis 10?

text: I have a XML file that would like to use as a base layer in ArcGIS 10. What are the steps that I need to follow in order to open that file in ArcGIS 10? The XML file does not have coordinates. When it's open it shows the following: ![enter image description here](http://i.stack.imgur.com/D67c7.jpg)


title: How to find out what features /tags are included in an .osm or .pbf file?

text: I would like to load osm data into a postgis (postgresql) database. In How to load OSM data into PostGIS preserving all Tags? there has been a solution posted:               Download the OSM files (in .osm format) and import using osm2pgsql.          When using osm2pgsql, there's a file called default.style that allows you to control      which OSM feature goes into which postgresql table, e.g.:          way        boundary     text         polygon     node,way   railway      text         linear          means that              features tagged "boundary" go into the *planet_osm_polygon* table         features tagged "railway" go into the *planet_osm_point* and *planet_osm_line*              because the "railway" tag is common between railway stations (osm nodes) and              railway tracks (osm ways).      Now I wonder how to find out which features are contained in the .osm or .pbf file to configure "default.style" Any input on how to make a good configuration file (default.style) would be appreciated. Thanks!


title: Sextante - SAGA- Shape to grid question?

text: I have a vektor shape (Index 1-5). After using the saga tool (shapes to grid) I'm getting a raster (values 1-5). The problem is that the NO DATA area is defined as 129. Is it possible to set 129 as no data? I cannot find the tool.


title: How to match a shapefile with seconds as units with other layers

text: I have some shapefile in **coordinate seconds** (longitude&latitude multiply 3600). I have open google street layer by openlayers plugin in qgis. how to match my shapefile with google street layer in correct position? thanks.


title: How to export only some feature classes to xml using ArcObjects?

text: I need to export some of the feature classes to an xml file ("export xml workspace") quite often so I'd like to automate the process. I've found an example here how to export a **feature dataset** to an xml file. I haven't been able to modify the code so that it would export only the feature classes I want (for example only the ones that start with a letter 'A'). I also need to export the tables and relationshipclasses (schema only). Is this possible at all? I know there is an interface `IFeatureClassName`, but I don't know how to use it in this case. Should this be done with arcpy instead? I'm using ArcGIS 10.0 SP3.


title: How can I get the GPS info from a PNG?

text: I am trying to get the GPS data from several hundred PNG photographs for display in Q GIS (so need to access the EXIF data). The online and downloadable viewers I have tried haven't yielded much information: ![ ](http://i.stack.imgur.com/wnEfL.png) Do the photographs not contain any GPS info ( I am told that they do) or is there a different way of extracting the information? Many thanks, Laura


title: Huge x,y,z dataset into oracle, possible?

text: I have 200 .xyz files containing 1500000 points each. Is this possible to import this dataset to oracle? It makes a total of 300.000.000 points!!! Background is I want to improve a FME script that cuts out a desired region of elevation data. Today it uses a .shp grid to find the .xyz files intersected. This is a done with two linked fme workspaces, very complicated procedure. If I can get it into oracle, an oracle reader with bounds would do the same job more straightforward.. Any input regarding elevation data in databases or FME retrival of elevation data is much appreciated. Best regards!


title: Georeferencing - Arcmap 10.2: Trouble with deleting control points

text: I'm trying to georeference a raster using a rivers vector (polyline) feature class. I started by shifting and rotating the raster. However when deleting a control point to reduce my RMS, both shift and rotation are "being automatically" canceled and I need to shift and rotate the raster back again. It only occurs in "Auto Adjust" mode.; Does somebody have a clue about the origin of the problem? Best,


title: Dense data set of locations for flow analysis

text: I want to do some analysis on real world data related to crowd flows and now I am looking for a publicly available (or at least available to academia) data set of trackings in a crowd. The data set does not neccesarilly have to be from humans, eg. locations of kettle's would also work. The required features are:   * Dense in the sense of many tracklets on a relatively small area during the same time window, the more the merrier   * Individuals/tracklets have to be identifiable (ie. it must be possible to link all locations added by one individual to some kind of id) at least for some time periods   * Bonus for underlying ground truth (eg. a fair going on) I have searched the webs for about 4 hours now but couldn't find anything that fits at least the first two requirements, any ideas are very welcome.


title: Speed differences between SQL Server 2012 and PostGIS UnionAggregate

text: There seems to be a massive difference in speed between PostGIS and SQL Server 2012 with a union aggregate function.               SELECT geometry::UnionAggregate(the_geom)     FROM loc     GROUP BY LocalityCo      vs               SELECT ST_Multi(ST_Union(the_geom))     FROM dummydata.loc      GROUP BY "LocalityCo"      I get less then a second in SQL Server 2012 but around 25 seconds in PostGIS. PostGIS version is: `1.5 USE_GEOS=1 USE_PROJ=1 USE_STATS=1` Is there anything I can do to speed up the PostGIS query or it's just the way it is? Note: Both have spatial indexes.


title: How to create a multi-story building infrastructure web map?

text: Hi please bear in mind that I've got absolutely **NO** experience working with GIS before you begin laughing at this. During a recent meeting, our client was bouncing around their complaints about how detailed our annual reliability reports were and openly admitted to just reading the covering letter and the introduction before skipping 500-odd pages to the recommendation section. They are looking for a means to visually represent the mechanical and electrical health of plant machinery which is not a major issue as we build a web portal with WPF dashboards, guages etc. which they seemed impressed with. However, because each building is multi-story they wanted a 'Google Maps' type interface control to navigate between building levels and ultimately zoom from floor --> room --> plant machine --> machine component. Again this is achievable through WPF controls but to get it to work we have to build each individual map from scratch. After searching the net... I ended up here. Most of the application resources I have stumbled across discuss working with existing SHP files and I am searching for resource that focuses on creating SHP files from scratch. Thanks in advance for any help/advice as it will be greatly appreciated.


title: Why does QGIS show too many labels per feature?

text: When i check to label polygon feature, one out of a total five polygons shows like 10 labels on top of it. That particular polygon is only one polygon, represented by one row in the table, so i dont see why it shows 10 labels on top of it.....? I havent selected any advanced label options, that could have maybe created a conflict. Any idea anyone? I know that in ArcGIS there is an option; show only one label per feature - BUT there no setting like that in qgis, that i have found. a screen of the problem ![enter image description here](http://i.stack.imgur.com/cApr2.jpg)


title: Inner shadow of a OpenLayers map container?

text: Not sure if this question should be here in the GIS section or if I should move over to a more common html/css place, but here we go: I'd like to have an inner shadow of my OpenLayers map canvas so that it looks like that the map is underneath the rest of the web page. Example: ![example](http://farm9.staticflickr.com/8427/7675433478_ecb2785e95.jpg) The square pattern is where the map will be later. Ideal would be to just use               box-shadow: inset 0 0 5px 5px black;      but the map tiles will be displayed on top of the frame which makes the shadow disappear. Any ideas how to achieve this only with CSS? I really don't want to use lots of divs with PNG files as backgrounds, but I guess I have to?


title: How to get JSON from Geoserver using AJAX request

text: I am building a web mapping application using Leaflet and a sort of customized WFS-T API to place and retrieve features from PostGIS. I have the script for putting features in PostGIS tables down solid, and I've managed to publish my tables as WFS services using Geoserver. My problem--which seems simple enough --is how to import those vector features back onto my map. I figured it would be easiest to go with the JSON option Geoserver provides, but when I make the below request using jQuery:                   $.ajax({         url: "http://localhost:8090/geoserver/cite/ows?service=WFS&version=1.0.0&request=GetFeature&typeName=cite:markers&maxFeatures=50&outputFormat=json&callback=?",         dataType: "jsonp",         success: function() {             //         }     }).done(function ( data ) {         console.log(data);     });      or:                   $.getJSON("http://localhost:8090/geoserver/cite/ows?service=WFS&version=1.0.0&request=GetFeature&typeName=cite:markers&maxFeatures=50&outputFormat=json&jsoncallback=?", function(data){           console.log(data);     });      I get the following error:               SyntaxError: invalid label [Break On This Error]         {"type":"FeatureCollection","features":[{"type":"Feature","id":"markers.1","geom...      Now, I know that these two functions are equivalent; what would seem to be a mismatch between specifying json in the url string and jsonp as the datatype is in fact correct because Javascript interprets the response as Javascript rather than true JSON. I've also figured out from reading other forum posts that the problem is the json script is not being properly contained in an object variable; the syntax error is due to missing parentheses on either side -- ({jsonobject}) -- like so. My question is, **how the heck do I fix the syntax to get JavaScript to work with it?** Obviously everything I want is coming through--it's just breaking because of a stupid syntax error that I can't figure out how to overcome. Any help on this would be much appreciated.


title: Create Zonal Grid in R

text: Suppose I have a n*n raster, and I want to create k square blocks (k can be divided exactly by n*n ) for zonal analysis: for example, when n = 4 and k = 4 a 4*4 raster is create with value               1 1 2 2     1 1 2 2     3 3 4 4     3 3 4 4      How can I do this in R?


title: PostGIS raster support

text: I have the following requirements to manage raster data 1) Convert raster-maps from BSB format to PNG or GeoTIFF   2) Cut the borders of the map   3) Merge individual maps to form a contiguous coverage map   4) Tile them for caching and loading   5) Serve them up on a WMS service to clients. I am using OpenLayers->GeoServer->PostGIS stack to manage all my vector data. I have been thinking of loading, storing & retrieve raster data the same way (in PostGIS). So my questions are 1) What is the benefit of storing rasters in DB over doing this in the file system?   2) What are the steps to achieve this and what tools/ commands I would use?   3) Would I leverage the DB cache or would I still need a tile-cache ?   4) What are the bugs/ gotchas in PostGIS raster handling that I need to watch out for? Thanks,   Ramesh


title: Are there any mirrors for QGIS documentation?

text: I have tried several times today (July 30, 2012) to connect to www.qgis.org (to get the online documentation). The site seems to be down. Do you know if there is a mirror to obtain the Lisboa (1.8) full PDF doc?


title: How to limit viewing to map boundary on Android?

text: I have an offline map of a city displaying on Android. Since my map is just small, when you scroll beyond the the boundaries of the map you see only a gray screen and you will probably be lost and could not be able to view the city map. How can I limit the viewing of the map to the boundaries of the city map? What should I do for the user to only scroll on the map area, and not beyond it? My resources are OSM Mapnik, osmdroid-android-3.0.1.jar library. Thank you for your time.


title: Force St_split() to return geometries in specific pattern

text: I'm Using PostGIS `st_split()` to split a polygon geometry by a linestring. Then I use `st_dump()` to return each of the splitted geometry and it is working fine. But the issue is I want to force st_dump() to return in a specific pattern. For example if the line is vertical, return left geometry first and then right as shown below. ![enter image description here](http://i.stack.imgur.com/MYaI7.png) While if the linestring is horizantal, return upper geometry first then lower geometry. ![enter image description here](http://i.stack.imgur.com/lXZJc.png) Currently, it return geometry randomly. Below is the snippet which return me splitted geometries               SELECT  ((ST_Dump(ST_Split(circle, line))).geom) As wkt     FROM (SELECT      ST_GeomFromText('LINESTRING(73.3864598921383 33.8985381624948,73.3962348946526 33.8646669883339)',4326) As line,     ST_GeomFromText('POLYGON((73.3733940124512 33.8850949267837,73.4058380126953 33.8936448555238,73.41313362121582 33.87818310806669,73.377857208252 33.8708431412792,73.3735656738281 33.8772567393002,73.3733940124512 33.8850949267837))',4326) As circle) As foo;


title: Attribute table in Tilemill - incorrect informations

text: I was editing OSM shapefile in `ArcGIS 10.1.` I have added new column to attribute table. Type of columns is "text/string". So content of the lines is only text (type of polygons). After adding this .shp file and opening of attribute table there in TileMill are incorrect informations in lines. Do you have some advices how to resolve this problem? I was trying "encoding" in advanced (TileMill) but it does not work still. See picture here.


title: How to fix the height in !text using MapFishPrint?

text: How fix height in text block using MapFishPrint config.yaml. For example in A4 size paper, my text change by content of ${content}. My config.yaml               layouts:       #===========================================================================       A4 portrait:       #===========================================================================         mainPage:           pageSize: A4           landscape: true           items:             - !columns               # columns can have an absolute position. In that case, they need a width.               absoluteX: 30               absoluteY: 565               width: 790               config:                 borderWidth: 0.3                 cells:                   - padding: 2               items:                 - !map                   width: 788                   height: 543                  - !columns               absoluteX: 600               absoluteY: 565               width: 220               config:                 borderWidth: 0.3               items:                 - !text                     != FIX HEIGHT OF ITEMS                   text: "${content}"                   backgroundColor: #FFFFFF                   fontSize: 12                   font: Helvetica-bold             - !columns               absoluteX: 30               absoluteY: 100               width: 100               config:                 borderWidth: 0.3               items:                 - !text                   text: "Loc"                   backgroundColor: #FFFFFF                   fontSize: 12                   font: Helvetica-bold             - !columns               absoluteX: 130               absoluteY: 100               width: 80               config:                 borderWidth: 0.3               items:                 - !text                   text: "IMG"                   backgroundColor: #FFFFFF                   fontSize: 12                   font: Helvetica-bold             - !columns               absoluteX: 210               absoluteY: 100               width: 400               config:                 borderWidth: 0.3               items:                 - !text                   text: "Titulo"                   backgroundColor: #FFFFFF                   fontSize: 12                   font: Helvetica-bold             - !columns               absoluteX: 130               absoluteY: 50               width: 150               config:                 borderWidth: 0.3               items:                 - !scalebar                   maxSize: 100                   type: bar_sub                   units: m                   vertAlign: middle                   barBgColor: #FFFFFF


title: How to select only polygons that overlap with another layer?

text: I have two shapefile layers created from different thresholds from a predictive model (Pink=80; Yellow=95). Is there a way to create a new layer of only those pink polygons that contain yellow polygons and effectively erase the pink polygons that do not contain yellow? If I use Geoprocessing - Intersect I will end up with just a copy of the yellow layer. ![enter image description here](http://i.stack.imgur.com/P5kg8m.png)


title: Alternatives to Shapefiles as open-source, cross platform dataset types

text: I am working on software that is very ESRI oriented, but a future version will likely not be using ESRI software. It uses Shapefiles and Geodatabases. I'm planning on getting all of my data to Shapefiles in anticipation for future versions of the software that will likely be on Android and other mobile devices. It appears that Shapefiles are the most common datatype for features in the open-source GIS world, but what are the others, and what benefits do they bring? I'm familiar with GeoJSON and KML, but I'm sure there are others. I would like to know all options, but I am particularly interested in dataset types best suited for storing on mobile devices (the data must be accessible without an internet connection).


title: Apply Database based Rules to Openlayers

text: im trying to apply rules, where the data is fetched from an external database. This is my code: // create a new style and add rules var defStyle = new OpenLayers.Style({ strokeWidth: 1 });                       defStyle.addRules([                      new OpenLayers.Rule({                     symbolizer: {                         fillColor: "#ffeecc"                     },                     filter: new OpenLayers.Filter.Comparison({                         type: OpenLayers.Filter.Comparison.LESS_THAN,                         property: "order07",                         value: "04000"                     })                 }),                      new OpenLayers.Rule({                     symbolizer: {                         fillColor: "#ffcc99"                     },                     filter: new OpenLayers.Filter.Comparison({                         type: OpenLayers.Filter.Comparison.BETWEEN,                         property: "order07",                         lowerBoundary: "04000",                         upperBoundary: "08000"                     })                 }),                 new OpenLayers.Rule({                     symbolizer: {                         fillColor: "#ff9966"                     },                     filter: new OpenLayers.Filter.Comparison({                         type: OpenLayers.Filter.Comparison.BETWEEN,                         property: "order07",                         lowerBoundary: "08000",                         upperBoundary: "11000"                     })                 }),             ]);                            var selStyle = new OpenLayers.Style({                 fillColor: "#ffaa00"             });                  // combine styles in a style map             var myStyleMap = new OpenLayers.StyleMap({                 "default": defStyle,                 "select": selStyle             });      what i want to do is something like this: I have a database with ID| Name * * * 1 | Bob 2 | Sam and i want the rules getting applied according to one of those values. Is it possible to include external database information that is not in the shapefile ?? In spoken words the original rule is "Color all vectors where the property:order07 is less than value:04000 with the color #ffeecc" and i want to change it into something like "Color the vectors where the property:fromDatabase xy: is value:xxyy with the color #xxxxxx"


title: What are spatial and attribute indexes?

text: What is the difference between a spatial index and an attribute index. What are their purposes?


title: Export query to identify what polygon certain points fall in

text: I'm using Arc GIS 10. I have shapefiles consisting of points (street addresses). I also have shapefiles consisting of polygons (ie, congressional districts, city council districts, etc... a separate shapefile for each kind of district) How do I run a query that exports a csv of the address and what district that address falls in?


title: How to query data from a vector map (shapefile)?

text: We want to obtain the specific place of the boundaries of a shapefile. We want to be able to get the points of each polygon corresponding to a specific place in order to use these data in an algorithm. GRASS is the GIS that we are using. We have the elevation map from SRTM and the vector map containing the boundaries from http://www.unsalb.org/. Can we query/retrieve the points contained in a specific boundary for all the polygons in the vector map? Sorry first time to use a GIS


title: How to remove roads from dataset

text: I've got a shape file that I loaded into PostGIS that has both road and parcel multipolygons and I need to remove the roads. I tried loading the shp file in qgis and it loads them both as a single layer. I need to remove roads as best I can without removing parcels. I can remove most of them by testing whether the address is equal to null. However this deletes parcels where that data is unavailable. My idea is to create another filter that tests the complexity of the geometry since most roads are complex polygons and parcels are generally simple (not all). Then I think I will find a "close enough" solution where most roads are removed, and little to no parcels are removed. Here is what I've looked at to no avial: st_isSimple : tests whether a polygon is anomalous (self intersection or tangent) st_numPoints : Does not work because it only looks at the first linestring and there are no linestrings Some others... Anyway, what ideas do you have for classifying a polygon as a road? Is there a query I can do? A feature I can extract offline? -N


title: How to design a diagram/framework to exhibit tool's underlaying structure

text: > **Possible Duplicate:**   >  How to generate schemas of GIS hardware / software / workflows? Can anyone give me some examples of framework diagrams that showcase/exhibit a tool, web-mapping application or similar that use a smorgasbord of interconnected elements, libraries, functions and tools. This is especially useful for communicating how your tool works when speaking publicly (per se) about the tool itself. I am in the process of creating a poster for an academic conference (GeoIDE Geospatial Conference) and need some help putting my idea onto paper. Note that I may be using the incorrect terminology... feel free to edit! This is a rudamentary proposal for what I want to display:                            DECISION   ALTERNATIVES                         ↑   ↑                          MAP <------------- ↰                          ↑↑↑                 |          MAPSERVER ----◇-> OPENLAYERS      DECISION MAKER          JQUERY --↥------------------------ ↵      Is there a best practise to this type of diagram/visual aid? What have you come across that helped you design something similar? Any help is appreciated. I'm about to break out the crayons and construction paper!


title: gsconfig installation in windows 7

text: ![enter image description here](http://i.stack.imgur.com/hqRuO.jpg)How do I install gsconfig in windows 7? I tried with gsconfig.tar.gz but the error is ImportError: No module named setuptools


title: Qgis making a backgound map for topolines on Garmin

text: Hi I am a new user of QGIS, I would like to make a Background map with Topolines in gmap or img format. So I started off with Opendem shape files. I have sucessfully combined the tiles I want into one shape file. I would like to use them as a background to view behind an openstreetmap map. On a Garmin Etrex 20 Now I just do not know where to go from here. Do these files need to be converted into raster files or what would be next? Can someone point me to a good tutorial? Best Regards


title: Distance analysis of 3-D surface in ArcGIS using VB.NET

text: Using VB.NET I created a 3D polygon in ArcGIS to model a fault plane and now I want to find two different distances from a point on earth's surface. How to go about this? Can some one help with the code? I have used following code to create polygon with 4 points. For example Lat1, long1, z1 are latitude, longitude and depth of the first vertex. Now I need to find the following distances:   1. Nearest distance from a point (Defined by its own lat and long) on earth's surface to the 3-D surface   2. Nearest distance from the same point to the polygon formed by projecting this 3-D polygon on earth's surface               Dim pPoint As IPoint = New PointClass()     Dim pPointCollection As IPointCollection = New RingClass()     Dim pGeomCollection As IGeometryCollection = New PolygonClass()     pPoint.PutCoords(lat1, long1)     pPoint.Z = Z1     pPointCollection.AddPoint(pPoint)     pPoint.PutCoords(lat2, long2)     pPoint.Z = Z2     pPointCollection.AddPoint(pPoint)     pPoint.PutCoords(lat3, long3)     pPoint.Z = Z3     pPointCollection.AddPoint(pPoint)     pPoint.PutCoords(lat4, long4)     pPoint.Z = Z4     pPointCollection.AddPoint(pPoint)     pGeomCollection.AddGeometry(pPointCollection)     Dim pOutGeometry As IGeometry = CType(pGeomCollection, IGeometry)     Dim pZaware As IZAware = CType(pOutGeometry, IZAware)     pZaware.ZAware = True      For the first distance, we can use "Near Tool" in ArcGIS but I am not knowing how to initiate that and find the distance.For the other distance we can form a 2D polygon by ignoring the Z values of vertices and use the same Near Tool. Please help me with the code for this idea.


title: MapInfo 12 problems with PostGIS database

text: I am simply trying to connect MapInfo 12 to my existing PostGIS database but after hours of trying various options I have hit a brick wall and need help. So far I have:               1) Downloaded and installed the 32bit ODBC PostrgeSQL driver     2) Open MapInfo 12 and enabled DBMS Catalog via the Tool Manager     3) In the DBMS Catalog I created Create Catalog>ODBC from the drop down     4) Select Data Source clicked the Machine Data Source > New and then User Data Source then scrolled to the bottom of the list to find PostgreSQL ANSI/Unicode options - I picked Unicode     5) In the dialog box I filled in my PostGIS conneciton details and hit test, and it successfully connected     6) Back in MapInfo I choose EasyLoader from the Tools list and and hit the ODBC button in the dialogue, then Machine Data Source and from the list pick PostgreSQL     7) Back in the EasyLoader dialog at the bottom i click Map Catolog' which successfully creates the mapinfo schema and table in my PostGIS database.      Then this is where it goes wrong, it now seems I need to make each table mappable               1) Table>Maintenance>Make DBMS table mappable>ODBC connection and again choose my PostgreSQL from the Machine Data Source     2) The new dialog box opens and it has all of my schemas and tables avaliable     3) I select one of the tables and click Open     4) In the Make Mappable dialog box I have      Index Type: PostGIS     Index Column: the drop down list is empty and this is where I think there is a problem      If I choose Index Type: MapInfo the Index Column comes up with geom which is my geometry column but it also populates the X and Y column as though it is trying to find cooridnates to make the geometry. My PostGIS geometries are standard Point, Lines, Polygons, MultiLineStrings and Multi Polygons Does anyone know how to make a table mappable? I am on MapInfo 12.02, PostgreSQL 9.3 64Bit and PostGIS 2.1 I have hundreds of tables to map, so surely there is a better way Thanks


title: How to copy selected features from one sde feature class to another with different attributes using Arcpy?

text: I am new to Arcpy scripting, I need help on copying the features to SDE feature class. My Scenario: We have 2 Geo databases running on Oracle 11g ArcSDE 10.2 . I have been asked to write the python script for query and copy the features from SDE connection1: Feature class1 to SDE connection2: Feature class2. The feature classes are not identical. I need to copy Shape and only 2 attribute values feature class1 to feature class2 and also populate the remaining attribute values of feature class2 with constant values. Please guide me in right direction to achieve this. My Code till now:               # Import system modules     import arcpy     import os     import sys     try:             ParcelSAMNo = arcpy.GetParameterAsText(0)             ModeofQuery = arcpy.GetParameterAsText(1)             TransactionNo= arcpy.GetParameterAsText(2)             arcpy.env.workspace = r"C:\Users\chenna.kishore\AppData\Roaming\ESRI\Desktop10.2\ArcCatalog\LIC.sde\LIC.LandPlan"             arcpy.env.overwriteOutput = True              fc = "LIC.LPLN_CADASTREPLOT"             where_clause = "PIN in (" + ParcelSAMNo + ")"             s_cursor = arcpy.da.SearchCursor(fc, ("SHAPE@", "PIN", "PDAREA"), where_clause)             arcpy.env.workspace = r"C:\Users\chenna.kishore\AppData\Roaming\ESRI\Desktop10.2\ArcCatalog\SDE@REGEO.sde\SDE.Parcels"             outFc="SDE.RPAR_PD_Parcels"             with arcpy.da.InsertCursor(outFc,["PIN", "PAR_AREA", "SAM_NO", "TRANS_NO", "SHAPE@"]) as rowInserter:                 for s_row in s_cursor:                         # Insert the new row into the Parcel Layer                         geom = s_row[0]                         row_values = [str(s_row[1]), s_row[2], 0, int(TransactionNo), geom]                         print row_values                         rowInserter.insertRow(row_values)                         #rowInserter.insertRow([str(s_row[1]), s_row[2], 0, int(TransactionNo), geom])             # Clean up the cursor             del rowInserter             del s_cursor             del s_row     except Exception as e:         print e.message     print "Success Full..!"      I cannot view the data that has been copied. But there is no errors from the script. Any Suggestions? Thank in advance.


title: How to find out if an ArcGIS Server SOM or SOC is installed on a machine?

text: I need to create a visual map of my current ArcGIS Server system including where software and data are installed. The System comprises of MSSQL, ArcSDE 9.3, ArcIMS, ArcGIS Server 10.0, and various servers which host Server Object Containers, and at least one Server Object Manager. I can´t see any SOC or SOM in the windows services. Can anyone give me clue how to get this information from my server environment?


title: Documentation for Esri's ArcPy Python site-package?

text: Has anyone found any informative documentation about Esri's ArcPy Python module?


title: Popup - a various set of attributes depends of selected layer

text: I would like display in popup a various set of attributes depends of selected layer. I have:                   var selectControl =  new OpenLayers.Control.SelectFeature(          [layer1, layer2],         {          clickout: true, toggle: false,         multiple: false, hover: false,         toggleKey: "ctrlKey", // ctrl key removes from selection         multipleKey: "shiftKey" // shift key adds to selection         }     );      and:                   function onFeatureSelect(event) {         var feature = event.feature;         var selectedFeature = feature;         var text = ""             if(feature.cluster) {                 for(var i=0; i<feature.cluster.length; i++) {     //and here should be something about which layer is selected eg.      //if ( selectControl == layer1 )      //then {first set of attributes} else {second set of attributes}      //... How to do htis?      //This idea doesn't work, I received undefined value of attributes          text += [i+1]+". "+feature.cluster[i].attributes.gis_id +"<br>";     }     else {         text = feature.attributes.gis_id;         }         var popup = new OpenLayers.Popup.FramedCloud("chicken",              feature.geometry.getBounds().getCenterLonLat(),             new OpenLayers.Size(100,100),             text,             null, true,              onPopupClose,             function() { controls['selectControl'].unselectAll(); }//         );         feature.popup = popup;         feature.popup.closeOnMove = true;//         map.addPopup(popup);     };


title: Is the order in which ArcMap loads Extensions, Toolbars and Commands defined?

text: When you extend Esri's ArcMap with both an Extension, a Toolbar and a number of Commands, is the order in which these categories are loaded by ArcMap well- defined? (I.e. are Commands always loaded before Toolbars, Extension always loaded before Commands, etc.)


title: Openlayers vector layer with WFS protocol not refreshed properly after Insert,Update and Delete operations

text: In my web map application, I have added the wfs vector layer for feature editing. Everything is working fine except the refresh function of wfs vector layer. Refresh function with force ture option not update layer properly.I am refreshing the entire map page to see the correct output. Please check my code and screenshot below for refrence. Please anyone help me to fix this. ![enter image description here](http://i.stack.imgur.com/Fdt7m.png) My code is:               var wfsRefresh;     function init(){      format = 'image/png';          bounds = new OpenLayers.Bounds(                 -288161.9289062237, -435.65815425360597,                 949406.755566389, 1300000             );     var options = {     controls: [],     maxExtent: bounds,     maxResolution: 5079.826789665053,     projection: "EPSG:404000",     units: 'm'     };          osMap = new OpenLayers.Map('map', options)          //Setup Base layer     tiled = new OpenLayers.Layer.WMS(                 "Oxford_Map", "http://xxx.xx.x.xx:8080/geoserver/cite/wms",                 {                     LAYERS: 'Oxford_Map',                     STYLES: '',                     format: format,                     tiled: true,                     tilesOrigin: osMap.maxExtent.left + ',' + osMap.maxExtent.bottom                 },                 {                                             buffer: 0,                     displayOutsideMaxExtent: true,                     isBaseLayer: true,                     yx: { 'EPSG:404000': false }                 }                               );            untiled = new OpenLayers.Layer.WMS(                 "Oxford_Map", "http://xxx.xx.x.xx:8080/geoserver/cite/wms",                 {                     LAYERS: 'Oxford_Map',                     STYLES: '',                     format: format                 },                 {                                             singleTile: true,                     ratio: 1,                     isBaseLayer: true,                     yx: { 'EPSG:404000': false }                 }                                 );             osMap.addLayers([untiled, tiled]);                  wfsRefresh = new OpenLayers.Strategy.Refresh({ force: true, active: true });            graphicsWFSlayer = new OpenLayers.Layer.Vector("Graphics Data", {       strategies: [new OpenLayers.Strategy.BBOX(),wfsRefresh ],              projection: osMap.projection,        //EPSG:404000       protocol: new OpenLayers.Protocol.WFS({         version: '1.1.0',         url: 'http://xxx.xxx.x.xx:8080/geoserver/wfs',         featurePrefix: 'cite',         featureNS: 'http://www.opengeospatial.net/cite',         maxExtent: osMap.getExtent,         featureType: 'ODC_Private_Graphics_Data',         geometryName: 'the_geom',         extractAttributes: true,                     srsName: osMap.projection,  //EPSG:404000         isBaseLayer: false,         visibility: true,         schema: 'http://xxx.xxx.x.xx:8080/geoserver/wfs/DescribeFeatureType?version=1.1.0&;typename=cite:ODC_Private_Graphics_Data'                 })     });          osMap.addLayers([graphicsWFSlayer]);      }          // This function will be called when we click save icon on map web page.          function saveFeature(){     var selFeat = graphicsWFSlayer.selectedfeatures[0];     var featCoordinates = '';     var geoType = selFeat.geometry.CLASS_NAME.split('.');     var vertices = selFeat.geometry.getVertices();     for (var nd = 0; nd < vertices.length; nd++) {        if (featCoordinates == '') {           featCoordinates = vertices[nd].y + " " + vertices[nd].x      }      else {          featCoordinates += ', ' + vertices[nd].y + " " + vertices[nd].x      }     }     if (selFeat.geometry.CLASS_NAME == 'OpenLayers.Geometry.Polygon' && vertices.length > 2) {      featCoordinates += ', ' + vertices[0].y + " " + vertices[0].x      featCoordinates = '(' + featCoordinates + ')';     }                 WebService.saveGraphicsFeature(keyId, featCoordinates);     wfsRefresh.refresh();      }          // This function will be called when we click delete icon on map web page.     function saveDeletedGraphicsFeature() {     for (var del = 0; del < graphicsWFSlayer.selectedfeatures.length; del++) {     var idSplit = graphicsWFSlayer.selectedfeatures[del].fid.split('.');     var featId = idSplit[1];         WebService.deleteGraphicsFeatures(featId);     }     wfsRefresh.refresh();     }                     //My webservice function is:      <WebMethod(enablesession:=True)> _      Public Function saveGraphicsFeature(ByVal ObjId As String, ByVal coordinates As      String) As String       Try         If updateGraphicsFeature(ObjId, coordinates) = False Then             Dim qryStr As String = "INSERT INTO ODC_Private_Graphics_Data (the_geom)    VALUES (ST_GeomFromText('POLYGON(" + coordinates + ")'))"             dbObj.Connect()             Dim insertCmd As NpgsqlCommand = New NpgsqlCommand(qryStr, dbObj.dbCon)             Dim affectedRows As Integer = insertCmd.ExecuteNonQuery             dbObj.Disconnect()             insertCmd.Dispose()             insertCmd = Nothing      End if      Catch ex As Exception          Finally         dbObj.Disconnect()     End Try     End Function          Private Function updateGraphicsFeature(ByVal ObjId As String, ByVal coordinates As   String) As Boolean          Try              Dim qryStr As String = "UPDATE ODC_Private_Graphics_Data SET the_geom =    ST_GeomFromText('POLYGON(" + coordinates + ")') WHERE gid = " & ObjId              dbObj.Connect()         Dim updateCmd As NpgsqlCommand = New NpgsqlCommand(qryStr, dbObj.dbCon)         Dim updatedRows As Integer = updateCmd.ExecuteNonQuery         dbObj.Disconnect()         updateCmd.Dispose()         updateCmd = Nothing              If updatedRows = 1 Then Return True              Return False          Catch ex As Exception          Finally         dbObj.Disconnect()     End Try     End Function               <WebMethod(enablesession:=True)> _      Public Function deleteGraphicsFeatures(ByVal ObjId As String) As Boolean         Try                  Dim qryStr As String = "DELETE FROM ODC_Private_Graphics_Data WHERE Sno = " & ObjId             dbObj.Connect()             Dim deleteCmd As NpgsqlCommand = New NpgsqlCommand(qryStr, dbObj.dbCon)             Dim deletedRows As Integer = deleteCmd.ExecuteNonQuery             dbObj.Disconnect()             deleteCmd.Dispose()             deleteCmd = Nothing              Catch ex As Exception              Finally             dbObj.Disconnect()         End Try       End Function


title: Automatic capture (Raster to Vector)

text: I would like to know the automatic process for Scanned hand drawn line Raster to vector conversion. Please can anyone help me.


title: How to remove spikes in polygons with ArcGIS 10/Python?

text: In the picture below you can see a polygon layer, representing steep slopes, extracted from a TIN. The polygons have a lot of spikes in them. I would like to remove these spikes while the rest of the borders stay on the exact same location. I would like to do this in a Python script, using arcpy or Python functions. Further info: ArcGIS 10.0, Python 2.6.5, polygon layer is in a GDB ![picture](http://i.imgur.com/GdE2uZW.png)


title: mapserver tolerance in mapfile or template

text: I set tolerance in a layer of my mapfile. It works fine :               LAYER     NAME emprise     STATUS on     TYPE POLYGON     CONNECTIONTYPE POSTGIS     CONNECTION "host=10.1.18.18 dbname=InrapTours user=PostGIS             password=xxxx port=5432"     DATA "geom from activite.emprise using unique gid using srid = 2154"     TOLERANCE 1 #tampon autour du clic sourise     TOLERANCEUNITS kilometers     CLASSITEM "tpeope" #symbologie suivant le champ tpeope         CLASS             NAME "diag"             EXPRESSION 'Diagnostic' #valeur du champ tpeope             STYLE                 outlinecolor 255 0 0              end #end style             LABEL #etiquette à lintérieur d'une classe                 COLOR 255 0 0                 TYPE TRUETYPE                 FONT sans                 SIZE 5                 ANTIALIAS TRUE                 POSITION auto                 PARTIALS FALSE                 MINDISTANCE 300                 BUFFER 2             END # end label             TEMPLATE "emprise_requete.html"      But I 've seen demo where tolerance can be changed using a field created in the template. I tried to do the same by adding :               <input type="text" name="map_emprise_tolerance" size="5" value=[map_emprise_tolerance]>km      Then i can type a number in the field but the result of the query is still the same i.e. tolerance in the mapfile is used. is there an hidden CGI variable to specifiy ? sort of Thx


title: ArcGIS 10.2 Geoprocessing Service based on Python Fails with IOError: "SDE" does not exist Failed to execute

text: I have a Python Script that inserts X,Y Point in the SDE Point Feature Class. It works fine when executed in ArcMap 10.2, but it fails with IOError when published as Geoprocessing Task service in ArcGIS Server 10.2 ![Create Location Python Script](http://i.stack.imgur.com/RNPob.png) It does not make sense to me why!? ![ArcGIS 10.2 Server Log IOError](http://i.stack.imgur.com/gxuKh.png) Any ideas what might be and how to fix it? Thank You.


title: file type needed to add layer containing contour lines

text: What file format is usable/recognized by qgis to add contour lines to a map? I'm working with Linux Mint 15, and qgis 1.7.5


title: High resolution global land cover

text: I'm searching for world land cover data with higher resolution than the global land cover. I'm working on fire prediction algorithms and I want to do a nice browser based demo were you can select a point anywhere in the world and do a fire front simulation. Land cover with 1km resolution is a bit sparse.


title: problems installing qgis 2.0 Add-ons

text: I have installed qgis 2.0, I need to install add-ons and NEXT I get the error 'The supplement package decompression failed. probably is broken or missing in the repository. might also want to make sure you have write permission to the plugin directory: / home/cordero/.qgis2 / / python / plugins "I searched the file path indicated by the error but apparently I have not the directory, and uninstall and reinstall qgis on my pc but the problem persists![enter image description here](http://i.stack.imgur.com/RIkuq.png)


title: Questions about EPSG:6372

text: I'm using QGIS 1.8.0 and as CRS EPSG:6372 a projection for Mexico, defined by the Instituto Nacional de Estadística y Geografía (INEGI). I have the following questions, I hope you can help me out. 1.- If I measure the minimum distance between points projected with that CRS, what are the units? meters?! 2.- I'm also doing some Vector Grids with maps projected with that CRS, when QGIS asks for the parameters for the grids (as in the attached image) are the units meters? So from the image I will be dividing the map into 10000 x 10000 m. squares, right? Thanks in advance!!! ![enter image description here](http://i.stack.imgur.com/p2F8w.png)


title: Where can I download free or open access census datasets for China, including their respective spatial boundaries?

text: I am interested in the equivalent of US Census Bureau open data and its relevant Tigerline spatial data, albeit for China. In Canada, we also have data released by Statistics Canada along with cartographic reference boundaries. I should make it clear that I am having trouble navigating the Chinese National Bureau of Statistics website, http://www.stats.gov.cn/english/. It may be in fact that this website does not offer any free data, but it is located elsewhere? Most importantly, I am looking for free/open data!


title: How to cluster group of points in Google Maps?

text: I have a lot of points (e.g. 100) in Google Maps and I want to cluster each 3 points near to each other in one group. If it is possible, how can I do that?


title: Mapserver + NDVI styling

text: I'm serving a NDVI coverage via Mapserver. So, one band, float32 values, [-1.0, 1.0]. I want to display the NDVI in greyscale, and generate a legend. I've tried 2 things:   * Create two dummy classes, just for testing: ![Dummy classes](http://i.stack.imgur.com/q3xro.png) So, my legend looks like this ![Legend](http://i.stack.imgur.com/Gu0ti.png) And the image, obviously, is red and blue: http://dl.dropbox.com/u/6599273/capturas/ndvi_en_gvsig.png But I'd like to see the image in greyscale... The I've tried the next thing   * Pre-classify values with PROCESSING "SCALE=-1.0, 1.0". I don't know what Mapserver does internally with this processing order, but I can see the image in greyscale: http://dl.dropbox.com/u/6599273/capturas/ndvi_gris_en_gvsig.png The problem is now I can't see any legend, because I've not created any class. I'd like to see all the colors of the grayscale image in my legend. How could I do it? Many thanks in advance


title: Performance testing of SOC machines

text: I would like to know a way to measure performance of SOC machines. We are in process of moving from physical servrs to VM Servers. Thanks Jay


title: Print out the properties of the attribute table in qgis?

text: I have to merge a lot of shapefiles into one great shapefile.   So I want to compare the properties of the attribute tables to see if there are differences between the tableproperties of each shape (type, length, precision, name, etc). With the intention to standardize all the shapefiles eg: ![enter image description here](http://i.stack.imgur.com/kGOTd.jpg) Is there a way to print out a list of the table properties as shown in the tab 'fields' or to export the properties into an .xls or .dbf? Or is there an automated way to do the standardisation? (Perhaps this is a subject for a new thread...) I am using QGIS 1.8 (QGIS 2.2 is also possible). Thanks


title: How to export a png file from a raster file? r.out.png has an error

text: I used r.out.png to export a png file, but it produced this error:               (Sun Feb 27 00:26:32 2011)                                                           r.out.png input=raster_bulacanElevation output=png_bulacanElevation.png              rows = 382, cols = 819     ERROR: Not able to open file for [png_bulacanElevation.png]     (Sun Feb 27 00:26:33 2011) Command finished (0 sec)


title: Is there a way of cutting a raster map that is displayed in the GRASS GIS Map Display?

text: I have a raster map that I generated from cropping/masking from a bigger raster map. It produced a map with only a portion of not null values. I only want to have a map that is zoomed in to the portion of not null values in order to save it into a PNG file. I only succeeded in setting the region (g.region) to the portion, but not cutting the map that only contains that portion.


title: How can I draw Google Fusion Tables on an Esri Silverlight Map?

text: I've been asked to add drawing support for Google Fusion Tables to an Esri Silverlight map application. For example this one showing Full-Year Stream Temperature Monitoring Sites for the Western US. Ideally, I'd like to be able to handle this just like any other mapservice layer by providing the url and add it to the map. Does anyone know of an easy way to do this? **Update** A related question is: would doing this be a violation of the Google Fusion Table Terms of Service? > 5.3. You agree not to access (or attempt to access) any of the Services by > any means other than through the interface that is provided by Google, > unless you have been specifically allowed to do so in a separate agreement > with Google. Does this mean I have to use Googles javascript API method if I want to show layers from fusion tables?


title: Convert a KML ImageOverlay to World Files with GDAL

text: I have many raster outputs from an Open Source radio transmission topographic path tool called SPLAT, georeferenced in KML image overlay format. I have figured out an easy way to use tools to convert from KML Image Overlays to PGW World Files using this small Java tool (which has fallen off the internet, I rebuilt it from Subversion repo):   * WorldFileTool   * 0.3.7 jar download link, my personal build _Is there a way to use GDAL to make this conversion more script-able?_ _If I were to implement a simple script to write out world files would there be any gotchas?_ I realized I had originally referenced Create world file from kml a few weeks ago, don't miss that one.


title: How to avoid Google map geocode limit?

text: I'm creating a custom google map that has 125 markers plotted via a cms. When loading the map I get this message: > Geocode was not successful for the following reason: OVER_QUERY_LIMIT I'm pretty sure it's the way in which I've geocoded the markers. How can I avoid these warnings and is there a more efficient way to geocode the results? UPDATE: This is my attempt at Casey's answer, I'm just getting a blank page at the moment.                      **UPDATE: It doesn't really matter what the solution is as long as the markers load instantly and it's not breaking any terms & conditions. This is all new to me in the first place so why not learn some more new things!!**


title: In ArcGIS 10 - how do I clip an arc

text: I have an arc that has violated the topology rule that the road feature should not extent outside of my boundary. I tried to use the clip tool in editor but cannot complete the clip. I should be able to click on the boundary, then use the trim tool to clip off the extended line. I selected the boundary, its lit up and shows in the TOC as selected, but nothing happens with the clip tool and I can't do anything with the offending road feature that is over the boundary line.


title: Where download sld styles for OpenStreetMap for GeoServer?

text: I display OpenStreetMap with some sld styles. ![my OSM](http://savepic.su/3172484.jpg) This map is not the same OSM. I want that this map is equal OSM map. I use geoserver 2.3.5 Does anyone display OSM through MapServer, ArcGisServer or Geoserver? \- here are some ideas. **Update:** The map may be not relevant OSM map. But it must be beautiful.


title: Baselayer disappears when I draw a layer

text: I have 3 layers, the base layer and two wms layers. When I load the map, one of them always shows up even when I put _visibility:false_. And when it shows up, the base layer disappears. What can I do? My code in openlayer is:                 map = new OpenLayers.Map('map', {controls: [], displayProjection: new OpenLayers.Projection("EPSG:32632"), numZoomLevels: 31, tileSize: new OpenLayers.Size(256,256)});             wms1 = new OpenLayers.Layer.WMS('base', 'http://172.26.10.126/fcgi-bin/mapserv.exe?map=/ms4w/apps/guinea/guinea.map', {layers: 'base'},{isBaseLayer:true, visibility:true,transparent:true});             wms2 = new OpenLayers.Layer.WMS('guinea','http://172.26.10.126/fcgi-bin/mapserv.exe?map=/ms4w/apps/guinea/guinea.map', {layers: 'guinea'},{isBaseLayer:false, visibility:true,transparent:true});             wms3 = new OpenLayers.Layer.WMS('municipios','http://172.26.10.126/fcgi-bin/mapserv.exe?map=/ms4w/apps/guinea/guinea.map', {layers: 'municipios'},{isBaseLayer:false, visibility:false,transparent:true});             map.addLayers([wms1,wms2,wms3]);             map.zoomToExtent(new OpenLayers.Bounds(123360.47,-162161.43,767730.83,418254.41));             map.setCenter(new OpenLayers.LonLat(8.80899,1.83925), 7);             map.addControl( new OpenLayers.Control.LayerSwitcher() );      And my mapfile is:               LAYER     NAME "base"     STATUS ON     PROCESSING "LABEL_NO_CLIP=True"     CONNECTIONTYPE POSTGIS     CONNECTION "host=localhost user=postgres dbname=Guinea_Ecuatorial password=tph"     DATA "geom from gnq_africa_central using unique gid srid=32632"     METADATA         WMS_NAME "base"         WMS_SRS "EPSG:32632 EPSG:4326"     END     TYPE POLYGON     LABELITEM "na2_descri"     CLASS         STYLE             COLOR 185 122 87             OUTLINECOLOR 0 0 0         END         LABEL             TYPE TRUETYPE             FONT "arial"             SIZE 15             COLOR 0 0 255                         ANTIALIAS TRUE         END     END               LAYER     NAME "provincias"     STATUS ON      PROCESSING "LABEL_NO_CLIP=True"     CONNECTIONTYPE POSTGIS     CONNECTION "host=localhost user=postgres dbname=Guinea_Ecuatorial password=tph"     DATA "geom from gnq_provincia using unique gid srid=32632"     METADATA         WMS_NAME "provincias"         WMS_SRS "EPSG:32632 EPSG:4326"     END     TYPE POLYGON     LABELITEM "provincia"     CLASS         STYLE             COLOR 85 185 107             OUTLINECOLOR 0 0 0         END         LABEL             TYPE TRUETYPE             FONT "arial"             SIZE 15             COLOR 255 0 128             ANTIALIAS TRUE             POSITION CC         END     END     LAYER     NAME "municipios"     STATUS ON     PROCESSING "LABEL_NO_CLIP=True"     CONNECTIONTYPE POSTGIS     CONNECTION "host=localhost user=postgres dbname=Guinea_Ecuatorial password=tph"     DATA "geom from gnq_municipios using unique gid srid=32632"     METADATA         WMS_NAME "municipios"         WMS_SRS "EPSG:32632 EPSG:4326"     END     TYPE POLYGON     LABELITEM "distrito"     CLASS         STYLE             COLOR 85 185 107             OUTLINECOLOR 0 0 0         END         LABEL             TYPE TRUETYPE             FONT "arial"             SIZE 10             ANTIALIAS ON             COLOR 0 0 0             POSITION CC         END     END      The layer that always show up is "distritos".


title: Is There A Term For The Old Maps That Had Lines Connecting Pin-Wheels?

text: Many old maps have lines like below. Is there a term for a map with that type of lines? ![enter image description here](http://i.stack.imgur.com/NqfHP.jpg) ![enter image description here](http://i.stack.imgur.com/TQ6wN.jpg) (from http://r-russell0912-dc2.blogspot.com/2011/01/pirate-maps-and-old-maps.html) My best hypothesis is that the lines may be related to Rhumb lines and therefore Mercator projections. Beautiful Example.


title: Ádvance Python Field Calculator - Global Expressions

text: I am trying to create a global expression that populates a new field based upon existing values in an attribute table. I am using QGIS Dufour and the Sextante toolbox, and the "Advanced Python Field Calculator" algorithm. Specifically, I wish to calculate a new field with the value of "Yes", where the values in another field are equivalent to "A", or "B" for example. I presume in the global expression option I need to include an IF statement, something along the lines of: _Global Expression_ IF FieldA == A OR FieldA = B, N = Yes Else, N= NO _Formula to apply to the new field_ Value = N Just wondering if this is possible? New to writing such expressions, and haven't been able to find much help on the subject matter, or an example script with regards to this algorithm and the correct syntax - would anyone know of a good info source with respect to learning to write such expressions in the advanced python field calculator? Thanks very much for your time.


title: Dynamically rendering polygon, point, and line features

text: I have a Feature Collection consisting of features having id, color_code, etc as attributes. I have to render this feature with the color in the color_code attribute. Is this possible and please suggest some links? Thanks in advance.


title: how to produce a Peirce quincuncial map?

text: As far as I know neither PROJ4 nor ESRI tools can apply the Peirce quincuncial projection. Does anybody know what libraries/softwares can manage it?


title: How can I remove a layer from layerswitcher in openlayers?

text: I have several vectors and WMS layers in layer switcher.But i don't need these layers to be displaying in layer switcher.How can i remove these layers from layer switcher? Thanks in advance..


title: Customizing the time slider in ArcGIS Java Script API 3.1

text: I want to customize the time slider. Currently I am using "esriTimeUnitsMonths" in which I am getting only months not the year (Please check the image![screenshot](http://i.stack.imgur.com/n1iGl.png)) I want to put month name +year ex. Jan 09 Feb 09 like that.. I researched on google but I am not able find suitable information. Some one told me that create the array but we can use array on time slider. Any help will be great..thanks in advance : ) Software : ArcGIS Server 10.1 & ArcGIS Java Script API 3.1 Code: `               var timeExtent = new esri.TimeExtent();             timeExtent.startTime = new Date("1/1/2009 UTC");             timeExtent.endTime = new Date("12/31/2011 UTC");             timeSlider.setLoop(true);             timeSlider.setThumbCount(2);             //timeSlider.setTickCount(dojo.date.difference(timeExtent.startTime,timeExtent.endTime,"Months"));             // original timeSlider.createTimeStopsByTimeInterval(timeExtent,1,'esriTimeUnitsMonths');             timeSlider.createTimeStopsByTimeInterval(timeExtent,1,'esriTimeUnitsMonths');             //timeSlider.createTimeStopsByTimeInterval(timeExtent,1,'esriTimeUnitsYears');             timeSlider.setThumbIndexes([0,1]);            // timeSlider.timeStops(3);             timeSlider.setThumbMovingRate(1000);             timeSlider.startup();`


title: Polyline "breakage" calculation in MapInfo?

text: I am a novice MapInfo user, with a basic to moderate knowledge of queries and SQL. My work entails doing risk assessments for shipping traffic based on AIS- tracks (radio-based vessel information system), which are represented as polylines in the dataset. However, some of the polylines are broken, due to missing AIS-coverage. I use several datasets, each of which have different places with broken tracks, and would like some way to quantify the deficiencies of each, as to make it easier to decide which dataset is best for a given location. I am wondering is if it is possible to use MapInfo to calculate what percentage of polylines enter an area (in this case a 1X1 nautical mile box in a grid of Europe) but does not exit the box afterwards, and then give that box a given color based on the percentage. Any ideas or input in regards to this would be greatly appreciated. If this question falls into the "RTFF", "RTFM" or similar, my apologies in advance.


title: 10.1 printing outside margin

text: I have seen on several maps in the case where a polygon overlaps (or maybe underlaps) the data frame, which is distributed to the margins, and text is turned on, that the text can print outside the data frame in the margin area. It shows this on both the pdf export and the print straight to a plotter.   Has anyone seen this same behavior?   ![margin printing](http://i.stack.imgur.com/okeEq.png)


title: How to use ArcGIS to reverse geocode and get road segment number?

text: I am now in a project which need to make reverse geocoding , namely ,given the latitude and longitude(assumed that the location is in a street), I want to get the road segment it is located in. What I want is not the road name ,but a road segment Id wich is a unique identifier of the **road segment**. As we all know , a road will always consist of several road segments .What I want is the road segment Id, not the road ID, road Id is too coarse for me .So , anybody can give me some suggestions?I am now have the road grid data,including a mxd file and the data files.And I have already installed a arcGIS client and import my data to in(illustrated in the screenshot).Because my project is a java application project, it will be very good if these exists a java api for me to call directly. ![enter image description here](http://i.stack.imgur.com/UKf4G.jpg)


title: Fiona - Preffered method for defining a schema

text: I'm fairly new to the world of python and GIS - so this is definitely a beginner question. I am following Tom Macwrights small guide to Shapely and Fiona, and I noticed that even when opening the previously created shapefile, it doesn't contain the schema I defined using:               schema = {            'geometry': 'Point',            'properties': { 'name': 'str' } }      This is evident because when I try to call it using `schema = input.schema.copy()`, it spits out `AttributeError: 'function' object has no attribute 'schema'` Reading the Fiona documentation, it states that: > the schema of its record type (a vector file has a single type of record, > remember) is accessed via a read-only schema attribute. Does that mean that it is not possible to store a schema in the shapefile,a nd therefore I'm always required set up the variables in the beginning, rather than reading from the file using `input.schema.copy()`? Or do I setup the schema using Records? If so, how do you read from a csv file (as in Tom's example) and place the data into the records to be read from `input.schema.copy()`? is it helpful to do this, or am I just complicating things and should define the schema every time i open the shapefile? Thanks for the help - I'm really interested in developing my _understanding_ of the processes occurring.


title: ArcGIS Server on Windows Server 2008 - .NET or Java Backend

text: I am finding myself having to "defend" why I have chosen to setup ArcGIS Server on .NET backend with IIS (vs Java and Tomcat). The host server is running Windows Server 2008 R2. I thought I remembered reading some type of ESRI slick stating that ArcGIS Server with .NET and IIS offers better performance vs Java and Tomcat, at least on Windows. But after a couple hours of searching, I can't find any definitive evidence. Can anyone state any hard, cold facts as to why ArcGIS Server on .NET and IIS is better on Windows? Links, comments, personal experiences welcome.


title: Clipping raster with vector layer using GDAL

text: I have installed GDAL using the Osgeo installer. How can I clip a raster layer with a vector layer programmatically? Is there any GDAL API available which can help me with this? I am using Python.


title: Volume of water in basin with reference to incremental height of DAM?

text: I have a DEM of my area. I did watershed analysis with Arc Hydro. I want to calculate "Volume of water in a basin with reference to incremental height of DAM" Is there any method, tool to give output with info of at h1-v1, h2-v2,......hn-vn(at h-height , v-volume) in tabular or any other form, if we give input values of "min height" at which dam footing with "incremental height" and do calculations up to "max height". Please let me know if any software/tool is doing above calculation?


title: ArcGIS incredibly slow when opening any toolbox tool

text: I've got a fresh install of ArcGIS 10.0, and whenever I try and open a tool from the toolbox (for example, Project Raster) it takes literally 5 minutes to open the dialog, and then at least 30 seconds for any clicks I make in the dialog (for example, to select input files) to register. Does anyone have any idea why this is? The rest of ArcGIS isn't too slow (at least, no more than normal) but using toolbox tools is just impossible! This doesn't happen on my laptop, and this computer is significantly faster than my laptop (this has 12Gb RAM, quad-core processor etc). Are there any known problems that might cause this?


title: GeoExt Tree Example Modification

text: The example provided in http://api.geoext.org/1.0/examples/tree.html shows all overlayes are checked but how do I make these overlayes unchecked while page is loaded?


title: Going round in Circles..Buffering/ArcObjects/VBA

text: I've just started learning VBA and Arcobjects, and at the moment just playing around tyring to do different things. What i am trying to do is Click on a tool, and allow the user to put a point graphic on the map, but also buffer the point(marker) that has been added. I can get the point on the map, but seem to be going round in circles getting the buffer to work. Can someone give me some hints ( i dont want the exact answer just some guidance on where I am going wrong. Code is here.               Private Sub Buffer_MouseDown(ByVal button As Long, ByVal shift As Long, ByVal x As Long, ByVal y As Long)          Dim pMxDoc As IMxDocument ' Create Object for current document     Set pMxDoc = ThisDocument ' Set object to current Document          Dim pMap As IMap ' Create object for current map     Set pMap = pMxDoc.FocusMap ' Set object to the current active map          Dim pPoint As IPoint 'Create object for point on map     Set pPoint = pMxDoc.CurrentLocation '  Set object to current location          Dim pElement As IElement ' Create object for graphic THAT IS ADDED TO THE MAP!!!!     Set pElement = New MarkerElement ' Set object to a new marker element (point graphic)          pElement.Geometry = pPoint ' Assign the geometry of the element to the Point.          Dim pGraphics As IGraphicsContainer ' Create object to hold the graphic     Set pGraphics = pMap ' Set object to current active map          pGraphics.AddElement pElement, 0 '  Set the graphic that will be added          Dim pActiveView As IActiveView ' Create object for active map display     Set pActiveView = pMxDoc.ActiveView ' Set object to the Active view          '--------------CREATE BUFFER------------------------'           Dim pSpaRef As ISpatialReference3      Set pSpaRef = pMap.SpatialReference           Dim pFCBuffer As IFeatureCursorBuffer2 ' Create Object for buffer      Set pFCBuffer = New FeatureCursorBuffer ' Set object for new buffer           pFCBuffer.Dissolve = True      pFCBuffer.ValueDistance = 100      Set pFCBuffer.BufferSpatialReference = pSpaRef      Set pFCBuffer.DataFrameSpatialReference = pSpaRef      Set pFCBuffer.SourceSpatialReference = pSpaRef      Set pFCBuffer.TargetSpatialReference = pSpaRef                Dim pCGLayer As ICompositeGraphicsLayer2      Set pCGLayer = New CompositeGraphicsLayer           Dim pGLayer As IGraphicsLayer      Set pGLayer = pMap.ActiveGraphicsLayer                pFCBuffer.BufferToGraphics pGLayer                       pActiveView.PartialRefresh esriViewGraphics, pElement, Nothing ' [partial refresh of the graphics element.            End Sub      Thanks Halil


title: How do I hide a feature layer service

text: I'm only just getting into dojo, so I have what most would probably consider a very basic question. I've built my own legend in a div and added divs for the legend items. This is all quite manual, but I just can't get my head around the dojo dijits. I've written the code to toggle the divs on and off (just changing their background colour), but what I now need to add is the code to show or hide the layers from the map. How can I use the ArcGISDynamicMapServiceLayer object I created to show and hide that entire service on the map after it's been constructed? All the examples I've seen look at the service and then list the layers inside it. I just want to add/remove or show/hide the entire service.


title: Is there a way to decrease the auto-hide delay for dockable windows in ArcMap 10?

text: I'd like for my dockable windows to hide a little faster than they do now. Is this a user or application setting within ArcMap or ArcGIS? More details: I'm specifically interested in ArcMap, but if it's one setting for the entire ArcGIS desktop--I'll take it.


title: How to create a polygon feature class from polylines in Python Arcgis?

text: I am trying to create a polygon from polylines. I've created this script but it gives errors of Arrays. Thanks.               import arcpy     polyline=r'C:\temp\temp_lines.shp'     Output_polygon=r'C:\temp\temp_polygon.shp'     polygonArray = arcpy.Array()     polygonArray.add(polyline)     cursor = arcpy.InsertCursor(Output_polygon)     row = cursor.newRow()     polygon = arcpy.Polygon(polygonArray, "")     row.SHAPE = polygon     cursor.insertRow(row)     del row     del cursor


title: Is there any place in which I can find a list of invalid database characters for ArcGIS 10?

text: Relating to my question about the split tool here, as suggested I ask a new question: is there any place with a list of invalid characters to use in a database in ArcGIS 10? edit: It is in relation to the error I got before which is "The name starts with an invalid character". It didn't say where beyond that.


title: Landsat in ArcGIS

text: I need to learn how to download, import and compose Landsat 8 in ArcGIS for natural colors (preferably in a GEOTIFF format). I have tried on my own but it seems that I am doing something wrong... can anyone help me and tell me the exact steps please? Thank you in advance


title: GeoPDF Support in GDAL

text: I installed OSGeo4W using a Windows installer... it contains GDAL 1.7.3. I can't get GDAL to recognize PDF as a file type using commands in MSYS shell. Does this version of GDAL not support PDFs? If not, what do I need to do to add support to PDFs for GDAL?


title: Is there a way to determine the datum used when all you have is a lon and lat?

text: I have a dataset of US intersections from 2011 that I'm going to be using in a project together with data based in WGS84. Included in the CSV file is a variety of fields two of which are the lat and lon. The accompanying documentation makes no mention of the datum that was used for the coordinates. Is there a way to calculate it? If not, because it's US data, would it be safe to assume it's one based on NAD83? I am going import this into a geography field in SQL Server 2012. I'm guessing that the SRID is either 4269 or 4326. I'd like to be able to determine which it is and transform it if necessary.


title: Openlayers getFeature from WFS with Specific ID

text: I'm working with OpenLayers and have the NHDPlus Catchments for the North Atlantic. I have built a web mapping application in Javascript for viewing and delineating basins from these catchments, however, we feel it is no longer necessary to delineate the catchments at run time and can do it with a CSV list of each Catchment ID and their respective list of child catchments upstream while just pulling from a feature request each ID from Geoserver and displaying it as a delineated basin. Can I be pointed to where I might find information on how to do something like this?


title: Using openstreetmap data in ArcMap

text: I posted this on the OSM help page but it was suggested I try here instead: I'm new to OSM data and I'm using the ArcGIS editor for OpenStreetMap within ArcMap. The data downloads correctly but the table doesn't seem to include the 'name' field (which is the field I really need). When I follow the instructions listed here, I receive "ERROR 0008000: The value is not a member of ALL." It works when 'ALL' is there but if I try and add 'name' it fails. Any suggestions?


title: Measure areas, lenghts in specific regions (Australia, Sout East Asia)

text: I'm in need to find a suitable projections to measure correctly areas and lenghts for 2 areas of interest. I want to measure both factors for: 1) First region Australia & New Zealand (1st region), 2) Second region - South East Asia defined as: Bornero, Honkong, Indonesia, Macao, Malyasia, Philipines, Singapoure, Thailand, Taiwan, Vietnam. (minlon: 94.89, minlat:-11.34, maxlong:141.056, maxlat:26.46) In regards to Australia & New Zealand:   1. Equal area projection: I could find two suitable projection to measure areas in Australia GDA94 / Australian Albers (EPSG:3577) and the Australian Albers (EPSG:17365). What about having a proper projection to measure areas at once in Australia and New Zeland? Is there any ready to use projection to do it or should I customize one of currently available, if so the question is how? More general question as nonGIS person: How does such projection measure area of relatively small shapes like shapes of the buildings in country? Does the area of interest the projection is applicable for somehow correspond to the size/ area of the object I would like to measure?    2. Equidistance projection Does the term distance correspond to the length in my understanding? For example I want to measure the length of each road I do have in the Australia's dataset I work with, but reading the definitions of equidstance projections I'm a litlle bit confused while in fact they're saying: Equidistant conic projection - It has the useful property that distances along the meridians are proportionately correct, Equidistant conic projection - It has the useful properties that all points on the map are at proportionately correct distances from the center point. What about measuring the lenght? Does the equdistance projection work for that, if not what would be the best option to go for region1/ region2 from my definition. Thank you in advance for an answer, Dawid


title: Spatialite Table Difference in QGIS

text: I am see differences between tables when I view them in spatialite-gui and in QGIS. I have created a table, added columns, and updated the columns in spatialite-gui. Within spatialite-gui I can view all the added columns, as well as the updated fields. However, when I loaded the same table into QGIS 2.2 and open the attribute table, I do not see the columns I have added. I am able to see all the columns in QspatiaLite, but not the attribute table in QGIS (see graphic below). I am using a Windows 7 64-bit computer with QGIS 2.2 from OpenGeoSuites4W, and spatialite-gui 1.7.1 The following are the queries I have written to create, add columns, and update the table.                Creating table with three columns, then adding geometry, and finally inserting geometry into the table from another using a union statement. This was first performed in spatialite- gui 1.7.1           CREATE TABLE smc_bldg2d (          gid INTEGER NOT NULL PRIMARY KEY AUTOINCREMENT,          city TEXT NULL, state TEXT NULL, country TEXT NULL);                SELECT AddGeometryColumn('smc_bldg2d', 'geom', 4326, 'MULTIPOLYGON', 'XY');                INSERT INTO smc_bldg2d          (city, state, country, geom)      SELECT f.city AS city, STATE AS state, COUNTRY AS country, st_unaryunion(st_collect(f.Geometry)) AS geom      FROM smc_bldg2d_usa_poly AS f      GROUP BY city      ORDER BY city;                ALTER TABLE smc_bldg2d ADD COLUMN type TEXT;                UPDATE smc_bldg2d      SET type = 'vector'      WHERE gid > 0;      ![enter image description here](http://i.stack.imgur.com/GIlD3.jpg)