title: How to add WFS layer from geoserver to openlayers

text: Ok, I know that this question has been asked here many times and I apologize for asking it again but I've looked all answers and I still can't fix my code to work. I would really appreciate if someone could tell me where I'm going wrong. I'm trying to add from geoserver vector layer with this attributes:   * layer name:marine,    * workspace:marine,    * workspace url: "http://www.geoserver.org/marine"   * epsg of layer on geoserver: 3765, epsg of map used in openlayers: 900913 This is part of my code that could help you see where I'm wrong:               <script>                     var geographic = new OpenLayers.Projection("EPSG:4326");             var mercator = new OpenLayers.Projection("EPSG:900913");                                     var world = new OpenLayers.Bounds(-180, -89, 180, 89).transform(                 geographic, mercator             );             var center = new OpenLayers.LonLat(16.00, 44.00).transform(                 geographic, mercator             );                         var map = new OpenLayers.Map("map-id", {                 projection: mercator,                 displayProjection: geographic,                 units: "m",                 numZoomLevels: 16,                 maxExtent: world});             map.addControl(new OpenLayers.Control.Permalink());             map.addControl(new OpenLayers.Control.MousePosition());                  var osm=new OpenLayers.Layer.OSM();             map.addLayer(osm);                       var bing = new OpenLayers.Layer.Bing({                 key: "AqTGBsziZHIJYYxgivLBf0hVdrAk9mWO5cQcb8Yux8sW5M8c8opEC2lZqKR1ZZXf",                 type: "Aerial",             });             map.addLayer(bing);                  var marine = new OpenLayers.Layer.Vector("marine", {                 strategies: [new OpenLayers.Strategy.Fixed()],                 protocol: new OpenLayers.Protocol.WFS({                     version: "1.1.0",                     url:"http://localhost:8080/geoserver/wfs",                     featureNS: "http://www.geoserver.org/marine",                     featureType: "marine",                     geometryName: "the_geom",                     srsName: "EPSG:900913"                        })             });             mapp.addLayer(marine);      somebody told on one topic that if I put my html file to geoserver-->data_dir-->www could help but it didn't. I spent more then half day trying to figure out where I wrong and I would really appreciate help. I just can't understand why is adding wfs layer such a big problem for everyone??


title: We sent Client Esri Files and Now they want to See Polygons

text: We are a voter list company. Every ten years new political boundaries are written for federal and state offices. All of the corresponding shape files are ESRI ready. One of our clients is getting voter data from us and wanted to supplement it with shape files. We do not use ESRI ourselves nor do we understand the raw data and polygons as we are quite illiterate in this area---we simply are giving the raw ESRI data to an end user. We sent the folder of shape files to a client who is also unschooled in ESRI and now they want to "see the Polygons". See sample data that we set the client at: http://www.srrproject.org/resources/redistricting-shapefiles/#shapefiles- content. My questions is the same as I posed to them: can they not "see the polygons" themselves from the shape data that we provided? Educate me on how these shape files are are to be employed. NEW as of Oct 10: Hi @Kirk, @RyanDalton, @Rayner, @PROBERT, @Aaron, I understand this is all free advice here in this wonderful exchange but I want to say that if what I am asking is above and beyond _free_ , we can discuss offline a small consultative gig for you by contacting us at: anne.horrigan@eMerges.com. But here is my specific challenge: The client has 50,000 complete addresses across the US. They want to link each address with its corresponding congressional, state senate and state house district. I am looking for step by step---aka, For Dummies--- instructions on how the client may accomplish this task on a regular basis using the shape files mentioned earlier. For example, what software, what to load, how to import & export.


title: How to fix "The selected file already seems to have a world file!" when image hasn't been georeferenced before?

text: I am georeferencing historical aerials from 1988 in QGIS 2.2.0. Sometimes it works perfectly and other times after I have places my points and start the georeferencer a window pops up saying: "The selected file already seems to have a world file! Do you want to replace it with the new world file?" Yes/No. I click yes thinking it will now take me to the transformation setting window to save my image, but nothing happens. This happens for images I haven't georeferenced before. Please let me know what I am doing wrong.


title: Shrinking the co-ordinates into a unique number

text: I have a bunch of data's in my gpx files. I would like to shrink those bulk data's(co-ordinates) into a single ID or sort of something. any suggestions?


title: How to convert OSM data to a projected Coordinate system

text: What is the easiest way to get OSM data into a shapefile projected data format for QGIS? Sounds simple but I've run into problems with this before.


title: What are some really cool web/mobile GIS UI examples?

text: I need some ideas on eye candy for our product. We have to display maps with lots of data on accessible through the web. Most sites I have seen have a Google, or OSM map with markers and layers overlayed. While most of them are functional - I feel like they have not not caught up with times. What are some cool ideas out there? Any examples of slides with custom maps, sliders, pan controls? Also highly interactive maps that are not just for the browser but for mobile apps as well.


title: Cast arcpy result as an integer instead (arcpy.GetCount_management)

text: I am trying to get a number by counting how many points are in a shapefile. And this works (I think), except I then am running into trouble using that number somewhere else. Eventually, I'll be using that count in some math (field calculator), but while debugging I'm running into an error that will end up causing me trouble later on. This code:               TotalPoints = arcpy.GetCount_management(Path_Pts)     arcpy.AddMessage(">>>> PROCESS: COUNT PATH POINTS {" + TotalPoints + "}")      gives this error:               TypeError: cannot concatenate 'str' and 'Result' objects      I tried casting it as a int, which it ALSO doesn't like:               TypeError: int() argument must be a string or a number, not 'Result'      So I've got a 'Result' object and need to turn it into a number. How can I do that -- or is using the arcpy function unnecessary or overly complicated here? (I'm learning arcpy and python as I go.)


title: How to show script tool dialog from a button in ArcMap 9.3.1?

text: Thanks to Jason I was able to properly set up the parameters of a script tool. I would like to get the same dialog box from the script tool to appear when used from a button in ArcMap. I know how to create a button and getting working, but only if it does not contain parameters. I'm assuming that this customization needs to be carried out using VB, or whatever other language. So, is there some generic code that can be used to create the same dialog box as the one that appears in the script tool? I'm not a programmer (although I can hack may way through stuff in Python) so I don't know where to start.


title: Is it possible to "strip" unwanted/excess data from a shapefile?

text: So I have a dataset that contains a unique identifier, and a .dbf file that I need to join to my shapefile. This is the only field that I can do a join on. The .dbf field is a 8-digit unique number, but the matching field in my shapefile has extra data both in front of and behind a number of rows. When I run the join, it understandably comes up blank as none of the rows match exactly. If the problem was reversed, I would just strip the extra data from the .dbf file. Is there a way I can do that in Arc? (using vrs 10, if that matters). Thank you!!


title: How can a browser argument be used to position a web-map

text: How can I enable the unique location and zoom level of a user on a web-map to be shared? For example, if a user zooms to `lat: 75, lon: 35, zoom: 8`, I would like to be able for the user to share that location through a link (for example, `www.mymap.com/75/35/8`). As I'm using OpenLayers, I understand how to record where the lat/lon/zoom of what the user is viewing, and also how to store this in a database if necessary, but I do not have a clear idea how (or if) this argument can be shared via the browser. My `OpenLayers` code is usually something of this form:               var mapCenter = new OpenLayers.LonLat(-35, 75);                  map.setCenter(mapCenter, 3);        For example, I would like `www.mymap.com/23/15/10` to work, and all other permutations of `lat/lon/zoom`.


title: Update, Insert Cursors don't work on Features that have annotation?

text: I've been working on a script to update my feature classes by grabbing new data from one feature class, and applying it to another feature class. It does this with several cursors, basically get value from the new data and set value to the old data. It works just fine on feature classes without feature-linked annotation. However, if the feature class has feature-linked annotation it will run through the script and not save the updates or appends. Is this a known limitation, am I doing something wrong, or is there a work around to handle this? My other post shows the layout of my code for reference.


title: Generating a vector contour map of fictional terrain

text: Originally I had asked this on the Graphics Design SE but after some realization of the question I was asking, it's probably better to ask it here. I'm making a fantasy world map in illustrator that I'd like to take to the next level by adding terrain contours. Currently, it's a flat political map like this map of Canada for example. I'd like to add an effect like what was used in this map here to show the relative depths of the water: ![enter image description here](http://i.stack.imgur.com/N9oSKl.jpg) This is obviously not done by hand. They would have used a program that compiled real terrain data into vector contours which would get exported to Illustrator/Indesign/editor of choice. I don't have the luxury of real data, so how would I go about generating my own using the map I've already drawn as a starting point? I have a good idea of what it is I need to do -- I'd like to export my existing map to a tool that will generate/edit terrain and then export it back as a vector for Illustrator. I'm just a total newb at GIS and don't have any idea of the software or tools of the trade. I've looked glancingly at things like R and Artificial Terrain Generation, and have used tools like QGIS in the past, but I could really use some direction on what's the best way to go about doing this and what software programs or libraries could help me. I'm no stranger to coding and know C and Python pretty well. I'd appreciate any help I can get. Thanks!


title: Remove a Layer using Arcpy

text: > **Possible Duplicate:**   >  How can I remove a feature layer in an Arcpy Script This question could seem too stupid, but i'm struggling with it, and it cannot work. I have an mxd file containing raster layers, and I would like to remove a layer. I don't know why it does not remove, it doesn't pop up any errors. I'm trying this code:               mxd = arcpy.mapping.MapDocument("photo.mxd")     df = arcpy.mapping.ListDataFrames(mxd)[0]          for i in arcpy.mapping.ListLayers(mxd , "TerraColor_SanFrancisco_US_15m.tif"):         arcpy.mapping.RemoveLayer(df , i)          arcpy.RefreshActiveView()      Any help please.


title: How to select multiple value in select by attributes in ArcGIS 10

text: Does anybody can help me with select multiple values from the attribute table? I try this command: "Classes"='14'AND'07' but doesn't work. Thank you.


title: How can I make sure that tiles display properly when you first go to the page?

text: This page is not showing my map tiles properly when you first go to the page. http://bear.hol.low.li/Pocket_Docent/map.htm If I resize the window it centers and distributes the tiles properly. How can I make sure it tiles and refreshes correctly when the page first opens?


title: How to get "p value" for time series data to detect significant change

text: I have yearly maximum NDVI data (1982-2008) and extracted all value to Excel. Also calculated slope value for each point during the 27 years. All of you know the slope value represents changes per unit (per year). Now I would like to look for significant changes in NDVI over the 27 years. So I have to get "p value" for each point. How I calculate it? I tried it in Excel. Firstly I have to estimate t-test, but i don't know which value should be in place of 'mu' for this case! After estimating it I can find p value using TDIST function in MS excel.


title: Customizing map.on("click") while using default editor widget

text: My question is whether or not it is possible to use the default editor widget in my case given the needs of my application. My map consists of points displayed from four different FeatureLayers. They represent four different types of traffic counts, and my application allows users to associate/dissociate them with traffic count data requests. (Count data requests are non-spatial records stored in a data table.) They are displayed on my page below the map in a Dojo DataGrid. If no traffic count request is selected in the DataGrid, the Editor and TemplatePicker should behave as normal (i.e. they should select points, display the Attribute Inspector, and allow the points to be dragged around). If a request is selected, the Editor should be **temporarily disabled** so that clicking the point toggles it as being associated/dissociated with the traffic count request clicked in the DataGrid. The AttributeInspector should not be displayed As it stands right now, if there is a request selected in the DataGrid, clicking a point on the map still brings up the AttributeInspector. Is it possible to customize the Editor in this way? The definition of my TemplatePicker and the logic of what I want map.on("click") to do is below. Thanks               function generateTemplatePicker(layers) {         widget = new TemplatePicker({             featureLayers: layers,             rows: "auto",             columns: 2,             grouping: false,             style: "width:98%;"         }, "templatePickerDiv");              widget.startup();              widget.on("selection-change", function() {             unselectRequest();             selected = widget.getSelected();         });     }          map.on("click", function(evt) {             if (selected) {                 var newPoint = new Graphic(evt.mapPoint, selected.symbol, null);                 selected.featureLayer.applyEdits([newPoint], null, null);                 widget.clearSelection();             } else {                 if (!evt.graphic) {                     if (requestEditingMode) {                          unselectRequest();                      }                      else {                          //*HIDE ATTRIBUTE INSPECTOR*//                      }                 } else if (!requestEditingMode) {                     //*SHOW ATTRIBUTE INSPECTOR*//                 } else {                     var pointAlreadyRequested = pointIsRequested(evt.graphic._graphicsLayer.id, evt.graphic.attributes.OBJECTID);                     if (pointAlreadyRequested) { removePointFromRequest(evt); }                      else { addPointToRequest(evt); }                 }             }         });


title: fail to create relationship class

text: I am trying to make a relationship class between two tables, but I get this message at the end: " Failed to create the specified relationship class". Can you tell me why? Thanks, Ligia Sorry I wasn't more precise:- I use ArcEditor licence level, ArcGIS 10.0 version, and we're talking about a file geodatabase. I have a table for "ducts" and a table for "cables". A "duct feature" may have several "cable features", the forein key is the ID of the ducts. It's a one-to-many simple relationship.


title: Georeferenced 3D Model Import

text: How can we open a 3D MAX File(3D Building model with exterior and interior elements)in Autodesk Infrastructure/LandXplorer or ESRI CityEngine? How can we georeference a 3D vector building model?


title: How to convert a vector layer to raster?

text: In qgis in the raster menu there is an item " _convert vector to raster_ ". This behaves strangely - I do not understand how it is supposed to work. I can select a vector file, an attribute and an output raster. But, for some reason the output file is only accepted if it exists and even if I select an existing raster nothing seems to happen. Also I do not understand how the program knows what kind of raster to create, what cell-size, projection etc. How can I convert vector to raster?


title: geotools parsing kml with multigeometry only gives me first geometry

text: I'm trying to parse the states.kml file (http://docs.geotools.org/latest/userguide/_downloads/states.kml) on the geotools website using the following code:               Parser parser = new Parser(new KMLConfiguration());     SimpleFeature f = (SimpleFeature) parser.parse( is );     @SuppressWarnings("unchecked")     Collection<SimpleFeature> placemarks = (Collection<SimpleFeature>) f.getAttribute("Feature");          for (SimpleFeature fe : placemarks) {         Geometry g = (Geometry) fe.getAttribute("Geometry");         for (int i = 0; i < g.getNumGeometries(); i++) {             System.out.println(g.getGeometryN(i));         }     }      Each of the geometry fields for each placemark in that file is a MultiGeometry made up of a Point and a Polygon. When I run through my for loop and print out though, only the points are printed. Does anyone know how I can access the Polygon also? Is it accessible, or did it get dropped in parsing? Thanks


title: editing existing lines

text: I have a quick editing question. I have a main line (street) from which I want to extend other small lines to reach some points (manholes), but I want to keep just one feature (the main line). I can't use the extent tool, because it works only with line-line case. ps: I am using ArcGis 10 (unfortunately;P ) Thanks, Ligia


title: percent of total population is x distance from nearest point (ArcGIS 10.2)

text: I'm getting my but kicked by some statistical analyses in ArcGIS 10.2 and need a bit of assistance. I have two data sets. One is a raster which represents the population of a country at 100m resolution, the other is 35k related points. Both population and the arrangement of points vary widely over space. From these two data sets I need to derive the following: 10th, 25th, 50th, 75th, 90th percentile distance between population and the nearest point. I have virtually no background in statistics - It appears as though I'm in a bit over my head at the moment and need a hand staying afloat.


title: Is there a web mapping server works on a linux shared hosting account?

text: **Question:** Is there a web mapping server software out there that works on a typical linux shared hosting account? **Background:** I have been trying find a web mapping server that would allow me to serve WMS services, WFS would be a bonus. I do not have root access and this is a shared linux account. It is a typical shared hosting account that can use php, ruby, python, cgi etc. I don't have root access. **Additional Information:**     * As with most shared hosts java isn't an option.   * I can install binaries locally, I just found out that I can request SSH access which I would assume would open up a few more options. Most shared hosting accounts wouldn't have this option.


title: Moving ArcGIS JSAPI Basemap Gallery Widget down 10px when opened

text: How can I move the Basemap gallery 10px down from its original position when it's opened. The basemap gallery sample, I am referring can be found here: http://developers.arcgis.com/en/javascript/sandbox/sandbox.html?sample=widget_basemapdynamic Also, how can I change the text "Switch Basemap" to an icon image? Any thoughts?


title: Mixing coodinate systems

text: I'm planning to create an online map with openlayers (2 or 3) or leaflet. I have to show some points/markers (~30.000), polygons (~500) and lines in an area of ~10x10 kilometers in several layers. This data comes in local x and y coordinates. The background will be made of raster tiles, perhaps OSM. Now I'm wondering, whats the best strategy to do this.   * Is it possible to make the local coordinate system the default for the whole map in one of those map libraries?   * If not, should the data be converted into the coordinate system of the map (WGS 84) in the database or a cache or doesn't that matter?


title: How to reclassify a raster (tiff) with multiple bands?

text: i have a raster file containing information about the windspeed average for the year. My purpose is to reclassify the 15 classes based on their rgb- composite. The color is the only information given. So i want something like this. A is my inputraster.               A.1 == 255 && A.2 == 175 && A.3 == 176 = 1      And that i want to do with all different rgb-composite. Any suggestions? I tried r.reclass r.recode and the mapcalc from QGis. Thx for your help. Update 18.11.2013: I'am still struggling with this issue. It seems like there is a erorr that has nothing to do with the mapcalc syntax. I just tried a easy expression, like "A*100". I get the following error: Problem loading output layer:               Oooops! The following output layers could not be open     Output raster layer: H:/reclassify/test2.tif     The above files could not be opened, which probably indicates that they were not correctly produced by the executed algorithm          Checking the log information might help you see why those layers were not created as expected          This algorithm requires GRASS to be run. A test to check if GRASS is correctly installed and configured in your system has been performed, with the following result:          GRASS seems to be correctly installed and configured      I want to look at the log tab from the algorithm, but when the error appears i can't get access to the log tab. I only can close the error message and after that the r.mapcalc windows are closed too. So i'am not able to mark any lines from the log tab. I really need help.


title: ArcMap Raster to NetCDF - lat/long

text: I need to convert a few files which are projected in WGS84 into netCDF. NetCDF requires lat/long. I am a bit confused on how the process is handled by ArcMap especially after looking at the online online Help for Arc10.2. Should I convert my raster to lat long projection before converting the file to netcdf or is the lat long conversion is automatically carried out when using the Raster to netCDF tool?


title: How to define correctly the Grass region in Qgis

text: I know that my question could be pretty basic for some of you. I am trying to learn GRASS. I have some raster layers and a vector file in EPSG 32614. Using the grass plugin I defined a new map set, defined a region using the limits of my raster layers. ![enter image description here](http://i.stack.imgur.com/CKVvE.png) ![enter image description here](http://i.stack.imgur.com/dL8mX.png) I succesfully imported both the vectors and the raster in to grass. When I try to run r. cost I have the following error: Raster map is in different zone (14) than current region (0) ￼Unable to open raster map Nevertheless, the metadada of my raster and vector layers are in 32614 (UTM zone 14N). The grass region settings tool tell me this : ![enter image description here](http://i.stack.imgur.com/oU5q4.png) I looked in to this question Raster map in different zone () than region () But It does not tell how ti fix the problem. I would appreciate any hint on the matter. I have the same problem working with 2.1 64 bits standalone and also in Lisboa 32 bits standalone working under win 7 64 bits. Thanks in advance Gerardo


title: WFS Feature Caching in QGIS

text: I am using QGIS 2 and using a WFS with a 5,000 feature limit. When I load a dataset without caching the features and it hits the 5,000 feature limit QGIS seems to give up on reloading data. When I zoom in on the dataset, it continues to request the WFS from the original extent, reoloading the data, but not in the new extent. The upshot of it is that when I load a dataset, hit the limit and zoom in, I still can't see the entire dataset without removing and re-adding the dataset Is this something wrong with the way I have QGIS set up, or is it something to do with the WFS settings?


title: converting a python list to arcmap table

text: I have values stored in a list that need processing in the costDistance tool in ArcMap 10.2 (desktop). For this they need to be converted into a table or file readable by arcmap's tools. Anyone with experience doing this?


title: Scale and Z factor have no effect on hillshade analysis in QGIS

text: I used an SRTM DEM to create a hillshade in QGIS 2.2.0. However, the resultant hillshade seems to be vertically exaggerated. Also, changing the settings for scale and z-factor does not have any effect on the resultant hillshade. I have tried both options in QGIS, Analysis-> DEM (Terrain models) and Terrain analysis -> Hillshade. Both methods gave the same results. Is this an isolated incident, a bug, or am I missing something?


title: How to recreate unique reference number with QGIS or PostGIS?

text: I have a shp file with about 500 odd records in it. It has an ID field but I want to add in a new field with a unique reference number. Can I use the field calculator to generate a unique reference? I’m looking for the reference to be alpha numeric, so something like UG_1 or UG1. Alternatively I could load the data into PostGIS and create the field and reference here but I’m still learning the ropes of SQL statements so I’m presuming that this is possible?


title: How to transform FileGDB featureclasses from custom local coordinate system?

text: Background: Started working for company that is in middle of system-wide IT upgrade including GIS data. Existing in-house data is in a FileGDB, is in an assumed coordinate system and represents 1000 sq\mi of service area in about 150 feature classes. I do have very good common data between existing dataset and 'real world' data (PLSS section data mostly). System: Working with ArcGIS Desktop 10.0 with ArcInfo license level. In old school command line ArcInfo I would have written an AML to run a Transform on the data coverages using the section corners for tic references. Not sure what the process is for File GDB. I looked into the Transform command, and it only appears to work at the Workstation software install level and not at the Desktop level, even though I have the Info license. Not sure if the Transform would work on a GDB instead of a coverage in any case. Question(s):   * Is there an equivalent Transform command for File GDB featureclasses?   * What are any other options for converting? UPDATE: My dataset is currently not in any coordinate system - it's floating\assumed\in-house coords. The data is good within itself but not at all related to any real world cood system. The old ArcInfo Transform command would use ID matched tics to move\warp\stretch data into a new spatial location. I do have data segregated in Feature Datasets with the File GDB. Thanks!


title: How to make Field Value constant size(length) by appending zeros towards the end?

text: I have been racking my brain, but I can't seem to find a solution to appending my field values in my .dbf to correspond with the values of my feature's attribute table. All I need to do is to append 0's to the field value if the field value is less than 11 characters in length. It needs to be done in the field calculator through VB Script or Python.


title: How to create polygon from xy data?

text: I have huge number of x and y data . these are left longitude,right longitude,buttom latitude and up latitude. each row is an extent of a polygon.i'd like to create a polygon using these data without programming in ArcGIS. is there any tool to create polygon automatically? ![enter image description here](http://i.stack.imgur.com/SQlyz.png)


title: rtree python polygon index

text: Sorry for a very basic questions about rtree. I would like to check whether the point is inside the polygon and the distance (point to the nearest point in polygon) as well. I was thinking about using No-SQL database, but both MongoDb and GeoCouch support Polygon records, but don't provide searches like: \- what polygons are around the point (for example, I want to find near what polygons a user is currently located; and after that check what in the precise distance from the user to the nearest polygon point) So, I want to figure out how to make these kind of spatial searches using python and related libraries only (shapely, geopy, rtree). Here are the steps I am following:   * I have shapefile with the targeted polygons (saved from QGis)   * I want to make Rtree Index, to improve performance when working with Shapely and geopy I am wondering, whether I can:   * make Rtree Index for shapefile directly? If not. I can use shapefile lib to collect polygon coordinates or bounding boxes:               >>> shapes[33].points     [[120.9937913, 24.7902669], [120.9939797, 24.790369], [120.994035, 24.7902849], [120.9938467, 24.7901828], [120.9937913, 24.7902669]]          >>> shapes[33].bbox     [120.9937913, 24.7901828, 120.994035, 24.790369]      Now I am a bit confused, and trying to understand how to create rtree index file for my polygons:               >>> idx = index.Index()  //creating index     >>> left, bottom, right, top = (24.7901828, 24.790369, 120.9937913, 120.994035) // can I put bbox coordinates here? I assume that coordinates are in the form [xmin, xmax, ymin, ymax]     >>> idx.insert(0, (left, bottom, right, top))       Would it be a right way for me to make Rtree Index?:               import shapefile     from rtree import index      sf = shapefile.Reader("polygons.shp")     idx = index.Index()     shapes = sf.shapes()     id_num = 0     for shape in shapes:         sh = shape.bbox         xmin = sh[1]         ymin = sh[0]         xmax = sh[3]         ymax = sh[2]         idx.insert(id_num, (xmin,ymin,xmax,ymax))         id_num = id_num + 1       After I execute the script above, I can get the following result from idx.bounds:               >>> idx.bounds     [24.7796055, 120.9839184, 24.7986571, 121.0025152]      To what bounds these coordinates are related to? Is it the bbox for the all the polygons? Is this a right way to check whether there are polygons near the Point? Assuming that the point coordinates are 24.789077, 120.999792               >>> list(idx.nearest((24.789077, 24.789077, 120.999792, 120.999792)))     [1, 2, 13, 29, 61, 62, 59, 60, 27, 28, 55, 5, 26, 54, 78, 49, 50, 3, 9, 10, 22, 72, 43, 39, 42, 85, 40, 8, 17, 76, 77, 36, 73, 71, 33, 34, 69, 63, 66, 64]      I am getting these results and not quite sure what they mean.. What does "nearest" exactly mean here? I cannot measure the distance here, like 5 meters.. so what is the "nearest" output here and how it's calculated?


title: How to fix "Error processing specified tolerance!" in Polygon Centroids tool?

text: I am using QGIS 1.8.0 to find the centroids of 92 features (i.e. smaller polygons) in a polygon layer. However, while using (Vector -> Geometry Tools -> Polygon centroids), I encountered the following error: "Error processing specified tolerance! Please choose larger tolerance..." Can someone please explain to me what went wrong and how I can fix it? Thanks!


title: Moving a shapefile point to specific coordinates

text: In QGIS is there a simple way to move a shapefile point in a shapefile point file to specific coordinates? Without exporting the points to WKT changing the x,y values and then importing the WKT via delimited text. I can move the point easily by hand through the Move Features button on the Digitizing toolbar and my current workaround is to import a simple KML file with the point and then by hand move the point feature of the shapefile on top of the KML point.


title: Input data to Maxent niche modelling tool?

text: I have following question regarding the input format of files to **Maxent**. From **here** I downloaded the environmental layers but they are in tiff format and I guess the tool needs in asc format where number of columns,rows etc and something like a matrix is there. How do I get the exact format required by the tool? How do I solve it.


title: Help with Parameters for MCC-LiDAR

text: I am attempting to use MCC-LiDAR to reclassify points (.las) originally classified using TIFFS. However, I seem to either be misunderstanding what the parameters really mean or where to get the information I need. Two parameters are needed. First, scale parameter which I have set to 0.75 since the point spacing (which is the same as post spacing?) is 0.75cm. Second, the threshold parameter which I have set to 0.5. Additionally, I set MCC-LiDAR to output an .asc grid at 1m resolution. The environment is forested with a mixed canopy, in leaf off conditions and the ground environment is karstic. Can anyone recommend what I should input for the parameters? I'm really having a hard time understanding them. ## Hillshade -- Azimuth: 135, Altitude: 15, dpi: 100 ![Hillshade -- Azimuth: 135, Altitude: 15, dpi: 100](http://i.stack.imgur.com/zvaCn.png)


title: Mapinfo Runtime is throwing error

text: when i am trying to call MapInfo from my PowerBuilder application, With map info pro(full) version "Open table" command is working fine but with mapinfo runtime(proviewer) "Open table" is crashing for the below syntax saying error calling external object. MIOleObject.application.Do("Open Table 'C:\Jail\DataDir\GcuPerm.tab' as Zones Interactive") Please suggest me what could be the problem?


title: arcmap envelope coordinates from webbrowser tool

text: So I have my own php file with textfields which I run in my add in. Now I want to get the variables from the textfields and automatically give the value of the textfields to my `env.PutCoords(x,y,x,y)` How do I do this?


title: Looking for free / opensource Lidar web viewer

text: I have lidar data which I want deploy my users using web viewer. I will probably hold the data in PostGIS (not sure yet). Is there good opensource/free lidar viewer that I can integrate with Openlayers site?


title: How can i uncheck leader lines show in piechart using ArcObject?

text: I'd like to uncheck the leader lines of my piecharts programmatically . I surf in the internet and found nothing about this.It's simple to uncheck the leader lines in ArcMap. Just go to properties of the piechart and Uncheck Show Leader lines but How can i uncheck the capability using ArcObject? What Interface and Property should i use to unchek th leader lines ? I'm using ArcGIS 10, ArcEngine 10,c# language Thanks for Advance.


title: How to call processing toolbox in Qgis' python console

text: I have been trying to use processing toolbox via python console in Qgis 2.0.1. I am able to write the following comamnds: ![enter image description here](http://i.stack.imgur.com/6nP4A.png) I am particularly interested to run 'r.statistics' from GRASS module. ![enter image description here](http://i.stack.imgur.com/qo1bU.png) My question is how to arrange the input parameters sequentially and what dfoes each parameter mean? For example, in the above screenshot, when I type the help, I can see the parameter inputs, processing.runalg('grass:r.statistics', base, cover, method, -c, GRASS region parameter, Grass region cell size, output) I can understand that 'base = input raster', 'method = avg, mean, std' and 'output=outputfile to eb declared'. But I am not sure what are 'cover=', '-c','GRASS region means'? How can I get further help on that? if I were to run the function 'r.statistics' via processing toolbox,how should the exact/correct command look like? processing.runalg('grass:r.statistics', 'C:\mydata\dem.tif', cover, '1', -c, GRASS region parameter, Grass region cell size, 'C:\mydata\outputdem.tif') Can anyone help me what to put in place of rest of the parameters? Thanks in advance.


title: Does arcpy.SearchCursor function support field indexing?

text: Does the arcpy.SearchCursor function support field indexing? In the sample script below, can I sent my cursor to loop through the values in the first field by setting an index to search the first field?               for x in arcpy.SearchCursor(filepath):         print x[0]


title: Extract polyline of border between differing polygons

text: I have a shape file with polygons with an associated variable that can take different values. I would like to extract the borders between polygons of differing value (red and blue) as a polyline (dark black line). I would prefer solutions in R but will accept other open source solutions. ![enter image description here](http://i.stack.imgur.com/aHkFN.png)


title: MongoDB: Search for polygons which have a certain location

text: Good day :). I am trying to implement a search feature on MongoDB using Mongoose.js , the short version of which is this: **I tag locations on a map as polygons. and give them names. Then when an user searches for a place using coordinates, I check which polygon the location falls within, and return the polygon.** I know that the reverse of what I am trying to do works by specifying a polygon, and then searching for locations which fall within the polygon, but is the opposite possible? So instead of saying:               var polyA = [ [ 10, 20 ], [ 10, 40 ], [ 30, 40 ], [ 30, 20 ] ]     query.where('loc').within.polygon(polyA)      I want to say               var 'loc' = [ 11, 24 ]     query.where(loc).within.polygon('polygon')      where `'polygon'` is a polygon schema. I believe it might be possible to do this with mongoDB 2.4 as it'll support GeoJSON, but is it currently possible to do it with mongoDB 2.2? Appreciate the help. **What I've thought of so far** I think it might not be possible even with MongoDB 2.4. So I'm thinking of using some linear algebra to solve this.   1. I compute the center of the polygon, and store it with the polygon area.   2. Then when I search the database I find a polygon with the center/area-point within a certain distance of the location.   3. Loop through each two points in the polygon, and set an `y `<?>` m(x) + c` type constraint, store the constraints in a JavaScript object.   4. Check if location coordinates satisfy the constraints, and return false if any of them is not satisfied. Else return true if satisfied. In effect I'm drawing a polygon's restrictions, and checking if the said location meets all of these. The problem here's that if I have complex polygons (which I'll more likely have) the test would take longer to return an answer. **On reinventing the wheel** _I sometimes get reprimanded for 'reinventing the wheel' by using MongoDB instead of a proper GIS database, which normally has some of the features I ask about. I am aware that what I am trying to do is what has probably been solved years ago by RMDBS, but I want to do it in Mongo. One prevalent way which this feature is likely already implemented is in geocoding, because to find the town/state/country which a point lies in, one would do the above. I thus sincerely want to try to do this in MongoDB as it complements the other computations that I'm already performing in Mongo. Thanks again_


title: How to map Traffic Count Points to Streets?

text: I've got point data representing traffic count locations and I am trying to figure out the best way to represent them on a map.I thought about kriging to provide a continuous surface but I don't think that is a good option. I could do a spatial join to join the points to the streets but I'm not sure how I could represent the actual counts. Any ideas?


title: default arcgis account credentials for 10.2 AMI

text: I'm installing the GeoEvent Processor on an ArcGIS 10.2 AMI Server and the install dialog is asking for the password for the arcgis account. Because the server was created with Cloud Builder, I never set a password for the account. Is there a default password for the arcgis account?


title: GPS Tracking by SMS

text: I have a GPS Tracker which sends the GPRMC string of the NMEA 0183 protocol, at fixed periodicity, to a predefined number as SMS. I am looking for a software which can port these incoming SMS to a COM port to plot them dynamically to Google Earth or a custom GIS application, and use the name of the sending GPS as a label on the map. THe software should be capable of simultaneously tracking multiple GPS trackers simultaneously.


title: Has anyone successfully installed HAZUS MR4?

text: I am attempting to install HAZUS-MH MR4 (v1.4) on a Windows XP SP3 computer with ArcGIS 9.3 SP1 installed. The installation seems to be running fine up until it attempts to start services. Then it displays this error message: ![alt text](http://i.stack.imgur.com/aULHW.jpg) I realize it's a bit of a long shot, but has anyone else installed HAZUS, recieved this error and found a successful solution for it? I've pretty much exhausted the solutions on the FEMA site, and am hoping someone may have a pointer in the right direction.


title: Add Cache as a Raster?

text: According to step 4 of the MapCruncher sample, http://help.arcgis.com/en/sdk/10.0/arcobjects_net/conceptualhelp/index.html#/Sample_Multithreaded_MapCruncher/000100000mmm000000/ _Once the cache is finished processing, open Engine or Desktop, navigate to the cache folder, and add the cache as a raster._ Basically I want to use the precooked maps in a dynamic display environment, but I'm not clear what they mean by "add the cache as a raster". How can I do this and do you think it would even work?


title: Ordnance Survey grids - download or generate?

text: Does anyone know a source from which to download vector UK Ordnance Survey grids at varying resolutions? **Or** failing that, a tool that can create an OS grid at a specified resolution? There are plenty of ways to generate a vector grid, but I've found no way of populating the attributes with OS grid references.


title: Where to download shapefile of 2km Tetrad squares for UK?

text: does anyone know where I can download a shapefile of 2km Tetrad squares for the UK please? A google search wasn't very helpful so I thought I'd try here. A Tetrad is derived from the British National Grid, but is not the same as the British National Grid ( http://en.wikipedia.org/wiki/Tetrad_%28area%29 ), so this is not a duplicate question. Thanks.


title: How to build a LAS dataset for simple viewing in Arcgis

text: I've been given a few LIDAR files to show my boss how useful the data can be. The problem I have is that I simply need to be able to view it using LASTools but I need to build a LAS Dataset first. I only have a Basic level license without Spatial/3D Analyst at the moment so I don't think I can build the dataset through ArcGIS. Is there another program where I can build the dataset and use it in ArcGIS?


title: What are some GIS outlets for my open source MapBasic utility module and library?

text: I'm actively developing a MapBasic module that provides additional MBX programming functionality by wrapping .NET 4.0 and Win32 API calls in easy to use MapBasic functions and subs. MB_Utils acts as an abstraction layer so a user can easily do things that are not directly supported in MapBasic: get Windows user name, display custom progress bars and dialogs, perform file i/o, etc. I've sort of been filling in the frustrating shortcomings of MapBasic I've come across. I'm looking fro some suggestions on where to post MB_Utils so as to gain exposure to other MapBasic developers. MB_Utils is hosted on GitHub: https://github.com/mfcallahan/MB_Utils https://github.com/mfcallahan/MB_Utils/releases


title: How to determine which polygons, if any, points fall within?

text: I am a amateur at GIS and need some help. I have a group of shapefiles of various areas and an Excel spreadsheet of addresses I have geocoded to points. I want to determine whether each point falls within the boundary of a shapefile, and which feature within the shapefile it falls within. For example I have:               Shapefile Zone: polygons Zone #1, Zone #2, Zone #3, etc.     Shapefile Area: polygons Area #1, Area #2, Area #3, etc.     Shapefile Address: points Address #1, Address #2, Address #3, etc.      For my results I would like it to list:                Address           Falls within feature     Address #1     Zone #2, Area #1     Address #2     Does not fall within any feature     Address #3     Zone #3      I know how to incorporate all the polygon features into a shapefile but not how to have it identify what polygon features of the shapefile each point falls into. I am using ArcGIS and Quantum GIS. So in summary my issue is taking coordinates and determining if they fall within certain polygon shapefiles I have. The goal in the end is to start with a list of geocoded addresses (points) and end with a determination of which (if any) polygon's each point falls in.


title: Is the geoxml3 javascript API compatible with the mapstraction API?

text: I am getting ready to revamp a google maps application and I was considering utilizing mapstraction. I am currently using all google maps js, with a kml layer overlay, but obviously by using the google maps js library I do not have as much control over the markers. In any case, I would like to utilize the geoxml3 to parse the kml but if I could also use mapstraction that would be an added bonus.


title: Set initial extent programmatically Android

text: I need to dynamically set the extent when the map loads. I have graphics that are being added to the map layer and want the extent to allow all the graphics to be shown. Currently I am getting a min x and y value and max x and y value based off the points that are being used to add the graphics. Then I try to use set extent with an envelope that contains the min and max values of x and y. So am I missing something to set the bounds so that the map will show the correct portion of the map?


title: LEAFLET - store and retrieve user markers & diagrams from database

text: I am intending to use leaflet in order to build a "microscope laboratory" for my institute's medical school. I am quite new to javascript, ajax and json technologies so I would like to have some general instructions on how to setup the whole architecture. I have my tile server (with microscope image tiles ready to use). I have set up my leaflet client application which composes tiles to a virtual map. I am also using some leaflet plugins (like PanZoom, Zoomslider, Draw, Fullscreen) in order to give users the option to mark and to draw on the picture. There are two further things I would like to do:   1. I would like to **save** (in a database? in an XML?) the markers and the drawings (circles, polylines, polygones) that a user has made on an image   2. Later, each time an image is requested and loaded into the browser I would like to show up those saved markers and drawings. Is there anyone who could help me to select the right tools, technologies and methods for implementing such a virtual laboratory?


title: Getting time dimension from a GRIB file

text: I have a GRIB file from NOAA that contains several variables for a specific subregion in Norht-East Atlantic. I'm using ECMWF grip-api Python library to parse the file. While I'm able to get the values for each {latitude, longitude} pair for each GRIB message with this code:               def iterate_latitude_longitude():     f = open(INPUT)          while 1:         gid = grib_new_from_file(f)         if gid is None: break              iterid = grib_iterator_new(gid,0)              #Get the missingValue for this GRIB message (e.g. NaN)         missingValue = grib_get_double(gid,"missingValue")              i=0         while 1:             #Get the value from the geo iterator obtained previously, value is a tuple with [latitude, longitude, value]             result = grib_iterator_next(iterid)             if not result: break             [lat,lon,value] = result                  sys.stdout.write("-Point # %d - lat=%.6f lon=%.6f value=" % (i,lat,lon))             #Check if missing value is present or not             if value == missingValue:                 print "missing"             else:                 print "%.6f" % value             #increase iterator             i += 1              grib_iterator_delete(iterid)         grib_release(gid)          f.close()      While this is great, this file contains forecast data, it was extracted from this place, more precisely, this file: multi_1.nww3.t00z.grib2 Which contains world data. There are other files for different times: 00-60-12-18-00, thus, every 6 hours. What I'm missing is the relation between the value obtained for each latitude,longitude and the time dimension. Does anybody know how to extract this info from a grib file, whether using Python or Java? I know I can use things like wgrib2 or pygrib to convert to CSV and do the other way around...but it just looks so inefficient Thanks! Alejandro


title: Mosaic (ArcGIS 10)

text: I would like to merge two raster into one image. Both rasters have value range different from each other. On one hand, the values range from (Image 1) 21.5 to 15.7. One the other hand are 18.15 to 13.3 (Image 2). Image 1 is found outside of Image 2 and they do not have overlapping areas. But I want both images to be displayed as one with the same value range and color ramp. I am using ArcGIS 10. Thank you.


title: Compare workspace schemas in Modelbuilder

text: My model has two model parameters:   1. A file gdb workspace or sde workspace containing the data   2. An empty "template" file gdb As the very first process in my model, I want to check if the schema of the data workspace matches the schema of the template gdb. If they don't match, the model should stop and preferably give some meaningful error message (I suppose by using some Python in a `Calculate Value` tool). I know this can be done simply by using the `Export XML Workspace Document` tool now in `10.1`, but I'm still on `10` for now, so I need my models to work properly on 10.


title: Looking for formula for the Robinson projection

text: I'm looking for a source that describes how the Robinson projection is calculated. I'm going to use this for converting long/lat into Robinson coordinates in both R and in javascript.


title: OpenLayers Popup relative to a feature's position - prevent overlaps - only one popup in one time

text: There are many questions related about this topic. However there is no complete example for this widely used pop-up type: Basically i need a popup which will appear the top-right of a Feature (Point) when user clicked on that feature. When i am able to do this, unfortunately pop-up overlaps part of the feature. So probably i should set an offset to the feature. Additionally when map is clicked, the popup should be removed and also when another feature is clicked, the old popup should be removed as well as new popup has been created. Which event listeners should i use FeatureSelected or FeatureClicked? Which popup type should i use Anchored or sth else? Could be projection issiues problem? I just started to use openlayers and javascript. Any help would be appreciated!


title: ArcGIS Publisher PMF Version Compatibility

text: Does anyone know of version compatibility issues of creating a PMF? Or links to that documentation? I have not seen anything from Esri. I have asked Esri Support, and will update with their response. Scenario 1:   * ArcMap 10.2 with publisher   * Microsoft SQL Sever Enterprise Geodatabase   * 10.1 Geodatabase   * Contractors use ArcReader 10.1 64 Bit Scenario 2:   * ArcMap 10.2 with Publisher   * 9.3 File Geodatabase   * Contractors use ArcReader 10.1 64 Bit Scenario 3:   * ArcMap 10.2 with Publisher   * 10.2 File Geodatabase   * Contractors use ArcReader 10.1 64 Bit I have the above 3 scenarios available. What compatibility issues would I have (or have to consider) when creating a PMF in my ArcMAP 10.2 for users with ArcReader 10.1? Obviously it would be appropriate for the users to upgrade their ArcReader to the latest version or for whatever version the creators of the PMF are using, but I do not have that control (contract limitations). Thoughts?


title: Problem with interpolation of 3D points in GRASS

text: I have a sample of 3D points (x,y,z) with tepmerature attribute. I can import my data as a 3D vector map. I want to create a volume of voxels by using v.vol.rst. When I use this module, this error happens WARNING: strip exists with insufficient data Processing all selected output files will require 0 bytes of disk space for temp files can somebody help me how to solve this problem? Thanks


title: Creating a DEM from CAD contour polylines

text: Is it at all possible generate a digital elevation model raster from CAD contour polylines with elevations? I am using straight AutoCAD 2014, not Map3D or Civil3D. I was able to bring the dxf of the contours into QGIS, but was not able to get the elevations into the attribute table. If I could get the elevations to transfer as an attribute, I think I could create the DEM from that, but as of now, I'm stuck at the polyline layer with no elevation attribute. Any thoughts on how to get attributes to transfer from the CAD file, and convert the lines?


title: How to maintain GlobalIDs while copying shape file data into Esri FGDB feature class?

text: I want to copy a shape file data into an Esri file geodatabase feature class. That feature class having a GlobalID column and it is being auto generated by the system while copying. The GlobalID values maintained as string in the shape file and i want to retain the same GlobalIDs in copied feature class too.


title: Administrative boundaries missing in OSM and/or Maperitive?

text: I'm lost trying to render maps with administrative boundaries using OSM and Maperitive. For some cities, e.g. Mexico City, there are boundaries. But, e.g. for Belo Horizonte, there are none, at least in the SVG I export for Illustrator (at zoom level 16) from Maperitive. I find references in the OSM XML for Belo Horizonte to `<tag k="boundary" v="administrative"/>` and `<tag k="admin_level" v="7"/>`. Also, on the OSM website, the map of Belo Horizonte has an admin boundary: http://www.openstreetmap.org/relation/368782#map=11/-19.9185/-43.9604 ...although sometimes you may have to switch to MapQuest view and then back to Standard for it to appear. So, it could be that no boundaries have been included in the OSM data for Belo Horizonte, or they have been and Maperitive isn't reading them. Advice very welcome... =====update==== The answer here: finding administrative boundaries in Guinea using OSM suggested using http://overpass-turbo.eu/ to get admin boundary data, but what I downloaded for BH (7MB) wasn't the boundary of BH. ===update 2=== Found this: http://global.mapit.mysociety.org/ but the search found no boundaries for BH. I guess there are none in oSM data.


title: Best way to combine 7 different shapefiles into one geodatabase

text: I have 7 different point shapefiles, and want to join them into one geodatabase. What would be the best method to do this in ArcGIS10?


title: Microsoft Access mdb "is not a valid or recognized data source" in QGIS

text: The answer to RyanDalton's question "Can QGIS read an ODBC connection?" explains how to establish a link in QGIS with an MS-Access table. Following his steps, I can establish an ODBC link (I successfully complete up to step 7), but when the .mdb is imported, I get an error that the file "is not a valid or recognized data source" in QGIS. I assumed QGIS would read each table as a non-spatial layer, as mentioned in the OGR_OGBC specs, but this is not happening. Can anyone tell me why the tables are not being read? Thanks in advance. I am using a 64-bit version of Windows with 32-bit QGIS and OGDB drivers.


title: Case insensitive CQL-OpenLayers filter?

text: I have the following filter in OpenLayers 2.13               var filter_c = new OpenLayers.Filter.Comparison({             type: OpenLayers.Filter.Comparison.LIKE,             property: 'p_name',              value: '%'+document.getElementById("nm").value+'%'     })      I also use GeoServer 2.1.3 and PostgreSQL 9.1 / PostGIS 2.0. This filter is case sensitive and is applied to a vector layer. I can provide more code if you want. I try to use `ILIKE` instead of `LIKE` but no luck. What am I missing? How can I make this filter case insensitive? Here is the complete, original code               var filter_c = new OpenLayers.Filter.Comparison({         type: OpenLayers.Filter.Comparison.LIKE,         property: 'p_name',          value: '%'+document.getElementById("nm").value+'%'     });          //pass the filter to the layer     var prot =  new OpenLayers.Protocol.WFS({         url:  "/geoserver/wfs",         featureType: "pins",         featureNS: "http://www.mysite.gr",         defaultFilter: filter_cl     });          var _CallBack = function(resp) {         pins.addFeatures(resp.features)         var cb = pins.features.length;         if (cb == 0){alert("Nothing Found");}     };          var response = prot.read({callback: _CallBack});               //refresh to render POIs     pins.refresh({force:true});      The `pins` layer takes data from a table in postgreSQL/PostGIS. This table has `id, name, geometry, category, etc.`Names are in `UTF8` and are in Greek, such as `Ακρόπολη, Παρθενώνας`etc. I simply set it like `pins = new OpenLayers.Layer.Vector("LayerTitle", {renderers: ["Canvas", "SVG", "VML"]})` and the I have the aforementioned code to search it


title: How to import AutoCAD .dwg files as overlay on map in Drupal (using OpenLayers)

text: I am using OpenLayers on my Drupal site. The request is to import AutoCAD .dwg files as an overlay to maps on the site. Is this possible and how? I found the OpenLayers Editor for Drupal, am I on the right path? If this is possible by converting the .dwg file to another format, it would be an acceptable solution, too.


title: Spatial join -different in model builder and right-click join

text: I am trying to incorporate a simple spatial join in a model but the spatial join options in model builder are different from the right click-->join options. How can I do the join as in the image below in the data management--> ![enter image description here](http://i.stack.imgur.com/MALyT.png) Above process works fine but I have to do it manually. The process below fails/doesn't give the result I want and I need to incorporate this in the model as otherwise I have to do this manually and then complete the process in a second model. ![enter image description here](http://i.stack.imgur.com/wL4WB.jpg)


title: How to read coordinates with variable Gauß-Krüger zone (meridian) in QGIS?

text: Is there a sensible way to visualize coordinates like this as point layer in QGIS? X values starting with 2 are EPSG 31466, 3 is 31467 and 4 is 31468. I didn't try, but i think the same problem will arise with UTM zones?               hauptstadt.txt:     ---------------     ID  Name    x   y     1   Kiel    3574000 6021000     2   Hamburg 3567000 5936000     3   Hannover    3551000 5805000     4   Bremen  3487000 5884000     5   Düsseldorf  2556000 5676000     6   Wiesbaden   3446000 5548000     7   Mainz   3447000 5541000     8   Stuttgart   3513000 5405000     9   München 4467000 5334000     10  Saarbrücken 2572000 5457000     11  Berlin  4593000 5822000     12  Schwerin    4461000 5945000     13  Potsdam 4573000 5807000     14  Dresden 5411000 5659000     15  Magdeburg   4475000 5777000     16  Erfurt  4433000 5649000


title: ArcObjects 10.1 Error: No such interface supported (Exception from HRESULT: 0x80004002 (E_NOINTERFACE))

text: I am trying to execute the below sample code Copying or loading data while preserving GlobalID values from ArcObjects 10.1 help. I am creating IWorkspace from FileGDB-Path using `workspaceFactory.OpenFromFile` and passing to this function.               // Converts a GlobalID field to a GUID field.     public static void ConvertGlobalIdToGuid(IWorkspace workspace, String datasetName)     {         // Open the table.         IFeatureWorkspace featureWorkspace = (IFeatureWorkspace)workspace;         ITable table = featureWorkspace.OpenTable(datasetName);              // Get the GlobalID field.         IClassEx classEx = (IClassEx)table;         if (!classEx.HasGlobalID)         {             throw new Exception(String.Format("No GlobalID column in table: {0}.",                 datasetName));         }         String globalIDFieldName = classEx.GlobalIDFieldName;              // Convert the GlobalID column to a GUID column.         IClassSchemaEditEx classSchemaEditEx = (IClassSchemaEditEx)table;         classSchemaEditEx.UnregisterGlobalIDColumn(globalIDFieldName);     }      But unfortunately, i am getting exception on the line               IClassSchemaEditEx classSchemaEditEx = (IClassSchemaEditEx)table;      Exception trace:               System.InvalidCastException was unhandled       Message=Unable to cast COM object of type 'System.__ComObject' to interface type 'ESRI.ArcGIS.Geodatabase.IClassSchemaEditEx'. This operation failed because the QueryInterface call on the COM component for the interface with IID '{F14139C3-B498-4B42-AA6C-5960028554E0}' failed due to the following error: No such interface supported (Exception from HRESULT: 0x80004002 (E_NOINTERFACE)).       Source=CMSS.Windows.FGDB.Utilities       StackTrace:            at CMSS.Windows.FGDB.Utilities.frmUtilities.ConvertGlobalIdToGuid(IWorkspace workspace, String datasetName) in C:\Dev Temp\CMSS.Windows.FGDB.Utilities\CMSS.Windows.FGDB.Utilities\Form1.cs:line 40            at CMSS.Windows.FGDB.Utilities.frmUtilities.btnConvertGlobalId_Click(Object sender, EventArgs e) in C:\Dev Temp\CMSS.Windows.FGDB.Utilities\CMSS.Windows.FGDB.Utilities\Form1.cs:line 86            at System.Windows.Forms.Control.OnClick(EventArgs e)            at System.Windows.Forms.Button.OnClick(EventArgs e)            at System.Windows.Forms.Button.OnMouseUp(MouseEventArgs mevent)            at System.Windows.Forms.Control.WmMouseUp(Message& m, MouseButtons button, Int32 clicks)            at System.Windows.Forms.Control.WndProc(Message& m)            at System.Windows.Forms.ButtonBase.WndProc(Message& m)            at System.Windows.Forms.Button.WndProc(Message& m)            at System.Windows.Forms.Control.ControlNativeWindow.OnMessage(Message& m)            at System.Windows.Forms.Control.ControlNativeWindow.WndProc(Message& m)            at System.Windows.Forms.NativeWindow.DebuggableCallback(IntPtr hWnd, Int32 msg, IntPtr wparam, IntPtr lparam)            at System.Windows.Forms.UnsafeNativeMethods.DispatchMessageW(MSG& msg)            at System.Windows.Forms.Application.ComponentManager.System.Windows.Forms.UnsafeNativeMethods.IMsoComponentManager.FPushMessageLoop(IntPtr dwComponentID, Int32 reason, Int32 pvLoopData)            at System.Windows.Forms.Application.ThreadContext.RunMessageLoopInner(Int32 reason, ApplicationContext context)            at System.Windows.Forms.Application.ThreadContext.RunMessageLoop(Int32 reason, ApplicationContext context)            at System.Windows.Forms.Application.Run(Form mainForm)            at CMSS.Windows.FGDB.Utilities.Program.Main() in C:\Dev Temp\CMSS.Windows.FGDB.Utilities\CMSS.Windows.FGDB.Utilities\Program.cs:line 22            at System.AppDomain._nExecuteAssembly(RuntimeAssembly assembly, String[] args)            at System.AppDomain.nExecuteAssembly(RuntimeAssembly assembly, String[] args)            at System.Runtime.Hosting.ManifestRunner.Run(Boolean checkAptModel)            at System.Runtime.Hosting.ManifestRunner.ExecuteAsAssembly()            at System.Runtime.Hosting.ApplicationActivator.CreateInstance(ActivationContext activationContext, String[] activationCustomData)            at System.Runtime.Hosting.ApplicationActivator.CreateInstance(ActivationContext activationContext)            at System.Activator.CreateInstance(ActivationContext activationContext)            at Microsoft.VisualStudio.HostingProcess.HostProc.RunUsersAssemblyDebugInZone()            at System.Threading.ThreadHelper.ThreadStart_Context(Object state)            at System.Threading.ExecutionContext.Run(ExecutionContext executionContext, ContextCallback callback, Object state, Boolean ignoreSyncCtx)            at System.Threading.ExecutionContext.Run(ExecutionContext executionContext, ContextCallback callback, Object state)            at System.Threading.ThreadHelper.ThreadStart()       InnerException:       I have tried re-install/new-install the ArcObjects SDK. But the exception still occurs. It would be helpful, if someone give some idea about this problem. Thanks.


title: Using cursors to update field that has relationships

text: I am still learning python and arcpy. I have spent too much time trying to figure out this code. I need help. What I have is two feature classes that have the same column names. One has older data than the other. FC1 has data that has 8 relationships (feature linked annotation). The FC2 is a master layer that gets update so attributes change. I cannot just swap fc's because the feature linked annotation is on a live web service. I tried to use Calculate Field in an Arcmap session (found out it has to be in an edit session). I started the Calculate Field and it ran for over 20 hours and I had to kill the session. I have to update 10 columns. I am tried to do this in Python with arcpy.CalculateField_management. Doesn't work and I posted a question. It was recommended to use cursors. The problem is I keep getting an error message that "Objects in this class cannot be updated outside an edit session [ttp.ttp.tintp_1]" Here is my code:               import traceback     import arcpy     from arcpy import env     import os          anno_conn = "\\\\path\\Connection to database.sde"     env.workspace = anno_conn          #Variables     Table_A = "tintp_1"     Table_B = "tintp_r3"     fieldlist = ["link_id"]     workspace = os.path.dirname(anno_conn)          try:              with arcpy.da.Editor(workspace) as edit:             prows = arcpy.SearchCursor(Table_B)                  for prow in prows:                 prior = prow.getValue("prior_visi")                 field = prow.getValue(fieldlist[0])                 print (str(prior) +  "    " + str(field))                      rows = arcpy.UpdateCursor(Table_A, "visi_link = " + str(prior))                      for row in rows:                          row.setValue(fieldlist[0], field)                     rows.updateRow(row)                     print "updated row"                      del row                 del rows                  del prow             del prows               except Exception, err:             traceback.print_exc()      Any help would be greatly appreciated. This has my experienced colleagues stumped as well.


title: What is the OGC way of representing Related features?

text: I have a point shapefile of Persons (The point is located at the registered address of the person) and a polygon shapefile of parcels. There is an owner_id attribute in the parcels shapefile, which links the owner in the persons shapefile to the the parcel ( For the sake of argument lets assume that one parcel can have only one owner, while one person can own multiple parcels) If I was working with ArcGIS, I would have created a relationship class to link these two featureclasses in the Geo-database, and then Published the data using ArcGIS Server. I could then use Query Related Records (Operation) on the Rest End point in my Web App. But I'm working with the Open source stack (Qgis+PostGIS+Geoserver+OpenLayers). What is the way of modelling these relationships in this stack? Right now I have the primary-key foreign-key relationship set up in the PostGIS Database, and have published these two layers as two different services. Whenever I need to find related records, I have to manually run WFS requests on the other layer. Is there a better, more standard compliant or best-practice way of doing this in the Open Source Stack?


title: Asp code for distance in SQL sever

text: I've looked around but I can't find the answer anywhere and im getting kinda frustrated now. I'm using sql server and what im trying to do is to find apartments within a 0.5km radious of a said apartment,I have both lat and lon in a varchar type of data. Ive read a few solutions but I cant seem to find how to put all that within the asp code. Thanks in advanced.


title: Best choice for building static maps: PostGIS, SpatiaLite or Shapefile?

text: I usually build static topographic maps, which are mostly printed. I'm often using OSM and rasterfiles (mainly different DEMs). The areas I'm mapping seldom exceeds 20.000 km². I'm using exclusively open source software like QGIS, SAGA and GRASS. I know the problems with shapefiles, but they are easy to handle. And Im not very familiar with postgis or SpatiaLite. I tried postgis, but got a lot of problems handling it, SpatiaLite seem to me more easy. What do you think: Should I stay with shapefiles or switch over to databases like SpatiaLite or postgis. Whats the advantages and disadvantages in my case?


title: How to calculate position in a picture based on latitude/longitude

text: I have got a picture on which a user manually marks some points. Let's say for simplicity's sake these points always create a 4 point convex polygon. For each one of these points I know latitude and longitude. Now the user enters latitude/ longitude for the 5th point and I have to show it on the same picture. What is the best way to calculate its screen position? ![enter image description here](http://i.stack.imgur.com/qRzqQ.jpg) I have tried the following algorithm:               1) to calculate a geometrical center of the polygon on screen          2) to find an average longitude for the leftmost and  the rightmost     points and an average latitude for the topmost and the lowest points          3) to place a point at the same relation to the geometrical center     as its' geocoordinates are to this average point.      I have a test image, where I know exact geocoordinates for each point on the 10px grid, and I can say that this algorithm doesn't work. What would be the better way?


title: OpenLayers - Resize a point on zoom out?

text: I need to resize a point on zoom. When you deep zoom the point get larger and on zoom out the point is smaller. In fact, the post here - How to resize a point on zoom out? seems to have exactly the answer. @CaptDragon posted an image of exactly what I need and then a demo, HOWEVER the demo is now dead. I've looked into the API that @Mapperz mentions but for whatever reason cannot get my point to resize larger when zooming in. If anyone can provide a working example similar to @CaptDragon - I would greatly appreciate it.


title: ArcMap throws NullReferenceException for OnClick event of Add-In

text: I have this below code which is a toolbar button to open a dock able window in ArcMap 10.2 Desktop application and what happens is that I am able to install the toolbar button but when I a click that button it shows Null reference exception in my try catch. **CODE:**               namespace DNRToolBar.FieldOffice     {     public class WaterRightsInfoButton : ESRI.ArcGIS.Desktop.AddIns.Button     {              public WaterRightsInfoButton()         {         }              protected override void OnClick()         {             try             {                 UID dockWinID = new UIDClass();                 dockWinID.Value = ThisAddIn.IDs.DNRToolBar_FieldOffice_WaterRightsInfoWindow;                 // Use GetDockableWindow directly as we want the client IDockableWindow not the internal class                 IDockableWindow dockWindow = ArcMap.DockableWindowManager.GetDockableWindow(dockWinID);                      if (dockWindow.IsVisible())                 {                     dockWindow.Show(false);                     this.Checked = false;                 }                 else                 {                     dockWindow.Show(true);                     this.Checked = true;                 }             }             catch(Exception ex)             {                 MessageBox.Show("This is the exception" + ex);             }              }              protected override void OnUpdate()         {         }     }     }      **The exception is thrown on on the second line after my try phrase at "dockWinID.Value". Even when I give it a value it still throws the exception no matter what.** **Config Code**               internal static class ArcMap     {      private static IApplication s_app = null;       private static IDocumentEvents_Event s_docEvent;            public static IApplication Application      {        get        {          if (s_app == null)             s_app = Internal.AddInStartupObject.GetHook<IMxApplication>() as IApplication;                 return s_app;        }      }           public static IMxDocument Document      {       get       {        if (Application != null)           return Application.Document as IMxDocument;               return null;        }        }      public static IMxApplication ThisApplication      {        get { return Application as IMxApplication; }      }      public static IDockableWindowManager DockableWindowManager       {        get { return Application as IDockableWindowManager; }      }      public static IDocumentEvents_Event Events       {         get        {           s_docEvent = Document as IDocumentEvents_Event;       return s_docEvent;        }      }      }      Thank you in advance!!


title: IZ.InterpolateZsBetween, how to use it for multiple parts in polyline

text: I want to interpolate Z values for the vertexes in the polyline. Shape file contains few polylines. Polylines do have certain vertexes with correct Z value, others do not have Z value (set to ie. 0). As on the picture below. How do I interpolate (linear interpolation) vertexes with 0 value? I figured that IZ.InterpolateZsBetween will do that. However, with my programming skills at this point, I am able to use basic geoprocessing tools in Python and realy have no clue how to implement ArcObject in it. Any ideas? Any known examples of that kind of code would be appreciated. ![enter image description here](http://i.stack.imgur.com/02JY0.png)


title: Example using sharpmap in mvc asp.net

text: I am wanting to open shapefile files in a ASP.NET MVC application, was researching and found sharpmap but have not found examples of how to implement the MVC. Does anyone have any examples to help me develop.


title: Can I change a bunch of address points to labeled points without having to do each one separately?

text: Hi I am trying to make a map with all of our customers in a certain area, I got the addresses on the map so now i have 300 dots. I was wondering if there is a way to make those labeled points without having to do all 300 separately? if so could someone please tell me. thanks


title: How to calculate perimeter of polygons in Quantum GIS?

text: Is there a way to calculate perimeter of polygons with the field calculator of qgis, like there is for area calculation?![enter image description here](http://i.stack.imgur.com/PgPwR.png)


title: Why is Spatialite delivering "Invalid SRID" error message

text: Qspatialite in QGIS 1.8 (and also Spatialite-gui) delivers an "Invalid SRID" error message when I try to load shapefiles with a SRID value of 3814 (Mississippi Transverse Mercatur). The SRID is within the database file spatial_ref_sys table. For troubleshooting purposes I tried loading the shapefiles directly into a blank QGIS session. QGIS builds its own custom CRS proj4 code and correctly reprojects the shapefiles on the fly against WGS84 Orthophoto (Bing Maps), but it doesn't attribute the projection as EPSG 3813, instead it's called "Custom 100003". I figured it should find EPSG 3814. Also, using the QGIS plugin, Qspatialite has a dialog that allows the user to scroll through a list of available projections, using both common and EPSG names, with SRID value listed as well. SRID 3814 is not listed in dialog, but I double checked and it actually is in the database's spatial_ref_sys table. Thanks for the user who had helped me discover that the SRID is already in the table, and doesn't need to be added. Can someone help identify the problem here? It's not getting read. EDIT: Here is the data you requested regarding my prior question (I couldn't copy it into a comment field) This shapefile correctly reprojects on the fly to WGS 84. EsriShapefile loaded into QGIS 1.8 (Ubuntu 10.04 -- UbuntuGIS Repositories version): USER:100003 - * Generated CRS (`+proj=tmerc +lat_0=32.5 +lon_0=-89.75 +k=0.9998335 +x_0=500000 +y_0=1300000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs`) myEsriShapefiles.prj ------ `PROJCS["mstm",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",1300000.0],PARAMETER["Central_Meridian",-89.75],PARAMETER["Scale_Factor",0.9998335],PARAMETER["Latitude_Of_Origin",32.5],UNIT["Meter",1.0]]` compare this to http://spatialreference.org/ref/epsg/3814/ 3814.prj -----               PROJCS["NAD83 / Mississippi TM",GEOGCS["GCS_North_American_1983",DATUM["D_North American Datum 1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.017453292519943295]],PROJECTION["Transverse Mercator"],PARAMETER["central_meridian",-89.75],PARAMETER["latitude_of_origin",32.5],PARAMETER["scale_factor",0.9998335],PARAMETER["false_easting",500000.0],PARAMETER["false_northing",1300000.0],UNIT["m",1.0]]      comparison of prj files below: 1>Esri.prj 2>3814.prj (the order of the parameters is different)               PROJCS["mstm",GEOGCS PROJCS["NAD83 / Mississippi TM",GEOGCS          ["GCS_North_American_1983",DATUM ["GCS_North_American_1983",DATUM          ["D_North_American_1983",SPHEROID ["D_North American Datum 1983",SPHEROID          ["GRS_1980",6378137.0,298.257222101]], ["GRS_1980",6378137.0,298.257222101]],          PROJECTION["Transverse_Mercator"], PROJECTION["Transverse Mercator"],          PARAMETER["False_Easting",500000.0], PARAMETER["false_easting",500000.0],          PARAMETER["False_Northing",1300000.0], PARAMETER["false_northing",1300000.0],          PARAMETER["Central_Meridian",-89.75], PARAMETER["central_meridian",-89.75],          PARAMETER["Scale_Factor",0.9998335], PARAMETER["scale_factor",0.9998335],          PARAMETER["Latitude_Of_Origin",32.5], PARAMETER["latitude_of_origin",32.5],          UNIT["Meter",1.0]] UNIT["m",1.0]]


title: Open Source GIS Problem

text: I am new at open source GIS platform. I am using Ubuntu 12.10. I have tried to install the followings but without 2-3 tools I could not install anyone. Can any body help me how can I install the following software from scratch. I mean dependencies on each other and installation priorities.   1. PostgreSQL   2. PostGIS    3. GDAL   4. GEOS   5. PROJ4   6. PGRouting (cgal & gaul)   7. Tilemill   8. Mapnik   9. Mod_tiles   10. QGIS


title: Should the IRelationalOperator.Equals return true regardless of vertex orientation?

text: I don't recall every having an issue with this method, but I am now. As I understand IRelationalOperator.Equals, vertex orientation shouldn't affect equality and if I take a polyline and reverse orientation, IRelationalOperator.Equals should still return true. However, I am not usually seeing this: IRelationalOperator.Equals usually returns false although occasionally it will return true after a `ReverseOrientation()` (very strange--I'm seeing this maybe 1 out of 20 calls). Is my understanding of this method correct? If so, is this a bug or is there something wrong with my implementation?               IPolyline polyline1;      IPolyline polyline2;           //...other code          polyline2.ReverseOrientation();     IRelationalOperator relop = polyline1 as IRelationalOperator;     Trace.WriteLine(relop.Equals(polyline2 as IGeometry));


title: How to programmatically access an Add-ins image?

text: From this thread: Can you programmatically change the button image for a ESRI.ArcGIS.Desktop.AddIns.Button? I've learned to change a command image programmatically, but... how can I programmatically get an image from the "Images" folder of the Addins, avoiding to write the physical file path? Is there a way to quickly get access to the Addins resources?


title: Contours generated from merged TIFs huge - much bigger than individually

text: I'm generating contours from four 5x5 degree SRTM DEMs in southeastern Australia. I'm using gdalwarp to reproject the DEM, then gdal_contour to generate contours. First, I generated contours for each block separately. The shapefiles were 9, 45, 73 and 270 MB. Then I merged the source TIFFs, and generated a single contour file. It's 2.0GB (that is, more than 4x the size of the combined individual contours). Why is this? Here are the two scripts: Individual contours:               for f in srtm_66_19 srtm_65_19 srtm_66_20 srtm_65_20     do     gdalwarp -s_srs EPSG:4326 -t_srs EPSG:3785 -r bilinear $f.tif $f-3785.tif           # contours     gdal_contour -a elev -i 20 $f-3785.tif $f-3785-contour.shp          done      Merged contours:               gdal_merge.py srtm_66_19.tif srtm_65_19.tif srtm_66_20.tif srtm_65_20.tif -o srtm.tif     f=srtm     gdalwarp -s_srs EPSG:4326 -t_srs EPSG:3785 -r bilinear $f.tif $f-3785.tif      gdal_contour -a elev -i 20 $f-3785.tif $f-3785-contour.shp      And the complete file listing:               -rw-rw-r-- 1 ubuntu ubuntu  20M Sep  9 17:36 new/srtm-3785-contour.dbf     -rw-rw-r-- 1 ubuntu ubuntu  399 Sep  9 17:08 new/srtm-3785-contour.prj     -rw-rw-r-- 1 ubuntu ubuntu 2.0G Sep  9 17:36 new/srtm-3785-contour.shp     -rw-rw-r-- 1 ubuntu ubuntu 7.6M Sep  9 17:36 new/srtm-3785-contour.shx     -rw-rw-r-- 1 ubuntu ubuntu 404M Sep  9 17:07 new/srtm-3785.tif     -rw-rw-r-- 1 ubuntu ubuntu 8.9M Sep  9 13:34 srtm_65_19-3785-contour.shp     -rw-rw-r-- 1 ubuntu ubuntu  68M Sep  9 12:28 srtm_65_19-3785.tif     -rw-rw-r-- 1 ubuntu ubuntu  69M Sep  9 12:28 srtm_65_19.tif     -rw-rw-r-- 1 ubuntu ubuntu  73M Sep  9 13:38 srtm_65_20-3785-contour.shp     -rw-rw-r-- 1 ubuntu ubuntu  67M Sep  9 12:28 srtm_65_20-3785.tif     -rw-rw-r-- 1 ubuntu ubuntu  69M Sep  9 12:28 srtm_65_20.tif     -rw-rw-r-- 1 ubuntu ubuntu  45M Sep  9 12:28 srtm_66_19-3785-contour.shp     -rw-rw-r-- 1 ubuntu ubuntu  68M Sep  9 12:28 srtm_66_19-3785.tif     -rw-rw-r-- 1 ubuntu ubuntu  69M Sep  9 12:28 srtm_66_19.tif     -rw-rw-r-- 1 ubuntu ubuntu 270M Sep  9 13:37 srtm_66_20-3785-contour.shp     -rw-rw-r-- 1 ubuntu ubuntu  67M Sep  9 12:28 srtm_66_20-3785.tif     -rw-rw-r-- 1 ubuntu ubuntu  69M Sep  9 12:28 srtm_66_20.tif      Why so big? Is there anything I can do about this? I'd prefer to work with merged files as they're easier to manage in TileMill.


title: Is it possible to export dynamic labels from ArcMap to KMZ? (Or alternate approach)

text: I am preparing increasingly more maps and layers for export to KML/KMZ as it seems to be the most convenient way to quickly analyze spatial data by my superiors. In order to export labels I seem to have to convert them to annotations. Even this approach is not without limitations; the annotations seem to be rasterized and almost always illegible unless I make them very large and adjust the pixels and DPI values of the output image. Is there a way to export labels directly, or as vector? If not is there another simple way to export labels to KML/KMZ? COMMENT: Thank you all for your suggestions. It seems that the existing workflow seems to be the best option. I exported my annotations to polygons and this workaround seems to work pretty well. Arc2Earth produces erroneous KMZ files that are riddled with errors. The recommended ArcGIS script is only available for legacy ArcGIS 9.x versions. Dynamic Labels DO NOT export at all into KMZ using the default tool.


title: Lidar data in ArcGIS 10.1

text: Problem: I have a set of LAS files that doesn't display in Arcmap. Details: I try to use LAS to tin in Arcmap to convert my LAS points to a TIN. However, when browsing to the folder with the files, it doesn't recognize them as LAS files: they do not show up even if I choose *All. I have tried a couple of conversion tools available in LAStools, but without success. I have also tried to use FME to convert them to shapefiles, but the only thing I get is the bounding rectangle. I seem to remember hearing something about this being a problem in Arcmap 10.1 with data that has been extracted using FME (which I'm not sure these are, but still), but now I can't find anything about it online. Has anyone else encountered this? Any possible solutions, workarounds, or at least a bug report? To answer @Andre Silva's comment: I'm trying to convert the point cloud to a raster to do hydrological analyses with it, e.g. calculate catchments etc. I have used other las files without issues before, so my interpretation is that this problem is specific to these files. As to the exact details of the files, they're <50 MB .las files that's been filtered for ground points only.


title: How to iteratively pass the values of particular fields to a convert polygon to raster process in modelbuilder?

text: I want to create a series of rasters that are representative, in both size and value, of particular individual fields of a vector feature class. I know how to use the 'polygon to raster' tool in ArcGIS to this effect but I am encountering difficulty in creating a process in modelbuilder that will automate a conversion process for each field. The specific problem is that I don't know how to pass the value of each field to the 'polygon to raster' tool in an iterative manner. I have used the 'iterate feature selection' in order to group the fields I want, illustrated below, however there seems to be nothing configurable within the "Value Field" of the 'polygon to raster' tool that will permit an iteration of the value of each field I have selected within 'iterate feature selection.' How can I iterate this process successfully? ![iterating feature selection polygon to raster ArcGIS Modelbuilder](http://i.stack.imgur.com/0lziW.png) The Value bubble in the model can be passed as a precondition, cellsize or environment:Output configuration keyword. When I run the model I encounter 2 errors: Failed to execute (Make Feature Layer)and Failed to execute (Iterate Feature Selection). Thank you in advance.


title: PostGIS. Obtain the rotation between two equal Geometries

text: I've got two features with the same drawing (poliline) but one is rotated relative to the other. I would like no know in PostGIS the rotation between the two geometries to apply this rotation to a third geometry (poliline). I know how to do the last step with the st_rotate function but i don't know how to do the first one. Regards,


title: Area precipitation using IDW

text: I'm having the following problem with IDW automation. I have a catchment and several gauging stations measuring daily precipitation amounts. I would like to interpolate the precipitation amounts from these stations (points) to the whole catchment. In the next step, I would like to calculate average precipitation value for the catchment. I can do it for one day (1 value for each station). The problem is that I have 30 years of daily precipitation measurements in each station, and so I need to do this process thousands of times. I have created a simple model in model builder which could do this process once, but I am not able to force it to do it repeatedly. I need help how to automate this process. Thanks Here is the modelbuilder scheme: ![enter image description here](http://i.stack.imgur.com/UAts3.png)


title: The International Date Line wrap around

text: Using OpenLayers, I added a WFS layer (on GeoServer) with a filter that returns all features (black) that intersect my polygon (yellow) placed over some Latin American countries within certain dates. ![enter image description here](http://i.stack.imgur.com/et7fs.png) However, the feature that crosses horizontally across the map does NOT actually intersect my polygon. This feature is somewhere in the pacific ocean between Hawaii and Fiji and NOT in Latin America. The problem is that instead of crossing the International Date Line it is being rendered on the map by wrapping around the whole world. The problamatic feature is defined: > POLYGON((-179.700417 14.202717,-178.687422 13.992875,179.024138 > 8.24716,-179.98241 8.035567,-179.700417 14.202717)) I have many of problematic date line features like this but have narrowed it down to this one for this example. I can't just ignore it in my application because i have many of them. I've tried using "wrapDateLine: true" in base layer and WFS layer with same results. Not sure if this would be a GeoServer problem or an OpenLayer problem. **Does anyone know a solution to my international date line problem?**


title: Process a raster with multiple geoprocessing tools in arcgis 10 using python - flow between processes

text: I'm writing a short python script that includes processing rasters with different geoprocessing tools in ArcGis 10, for example, the three processes below, is this best way to use multiple tools successively, i.e. create temporary rasters in the scratch workspace and then delete them all at the end of the script? By "best", I guess I mean is this the fastest/most efficient way to do this?               # calculate distance (inR is the input raster)     tmp1 = EucDistance(inR, "", "1000", "")          # convert to integer raster     tmp2 = Int(tmp1)          # Reclassify again     arcpy.CalculateStatistics_management(tmp2, "1", "1", "")     tmp3 = Reclassify(tmp2, "VALUE", "0 NODATA", "DATA")


title: How to cut one ESRI Shapefile with another ESRI Shapefile?

text: I have one ESRI Shapefile as a map with many objects in it. I have another ESRI Shapefile only with borders in it. This shapefile is bigger in size (covering a bigger area) with many border lines. The border lines form many mosaic pieces (tiles). How can I cut a part out of the map shapefile with the border shapefile as one mosaic piece? Where is the answer described in the Manual (en or dt) or in the answer to another question?


title: How do I fire a function after all features have been selected by a box in OpenLayers?

text: I have a map rendering with Polygons on a Vector layer set up as features. I set up the features so they fire off a function when the feature is selected and this works fine. The issue I'm having is when I try to allow users to select multiple features using a box. This causes the function to fire for each feature that is selected. I understand why this is, but what I really need is to have a function that fires once all of the features are selected so I can run a function using all selected features instead of having it run once for each. Here is my selectFeature code:               window[mapName].addControl(control);                             control.activate();                             var selectCtrl = new OpenLayers.Control.SelectFeature(layerVector, {                                 clickout: true,                                 onSelect: onFeatureSelect,                                 onUnselect: onFeatureUnselect,                                 box: true                             });                             window[mapName].addControl(selectCtrl);                             selectCtrl.activate();                         }                     });      It looks like there's away to attach Handlers to Boxes but so far I've been unsuccessful. Also, I'm afraid that if the function runs after the box selection is finished, it might not be getting the full list of selectedFeatures. Any help would be GREATLY appreciated. Thanks


title: How to compare two DEMs and show their difference?

text: I want to compare my DEMs (ASTER GDEM2 and SRTM4.1). I know, that the GDEM has some extreme "holes" with totaly wrong elevations (I found 3 of them in my working area). Now I need to know, if there are more of them to eventually correct them with the SRTM. I use SAGA, GRASS and QGIS. Thanks!


title: Publish SQL Views from Postgres/Postgis to Geoserver

text: I've trying and reading how to publish a Postgres/Postgis View in Geoserver. I know that I can create an SQL View in Geoserver, but I'd like to define the View in Postgres, not in Geoserver. Does anyone has an example of how it's done? So far I've tried this, but I can't see the content of the View, not even in Postgres.               CREATE OR REPLACE VIEW camino_view AS      SELECT c.id, c.nombre, t.id AS tipo_id, t.nombre AS tipo, l.gid AS geo_id, st_transform(l.the_geom, (-1)) AS the_geom     FROM camino c, camino_tipo t, camino_layer l, camino_link e     WHERE c.tipo_id = t.id AND e.camino_id = c.id AND e.camino_geo_id = l.gid;          CREATE TABLE camino_layer     (     gid serial NOT NULL,     "Tipo" character varying(50),     "Nombre" character varying(50),     the_geom geometry,     CONSTRAINT camino_layer_pkey PRIMARY KEY (gid),     CONSTRAINT enforce_dims_the_geom CHECK (st_ndims(the_geom) = 2),     CONSTRAINT enforce_geotype_the_geom CHECK (geometrytype(the_geom) = 'LINESTRING'::text OR the_geom IS NULL),     CONSTRAINT enforce_srid_the_geom CHECK (st_srid(the_geom) = (-1))     )     WITH (     OIDS=FALSE     );          CREATE TABLE camino_link     (      camino_id integer NOT NULL,      camino_geo_id integer NOT NULL,      CONSTRAINT camino_link_pkey PRIMARY KEY (camino_id, camino_geo_id),     CONSTRAINT camino_link_camino_fkey FOREIGN KEY (camino_id)       REFERENCES camino (id) MATCH SIMPLE       ON UPDATE NO ACTION ON DELETE NO ACTION,     CONSTRAINT camino_link_camino_layer_fkey FOREIGN KEY (camino_geo_id)       REFERENCES camino_layer (gid) MATCH SIMPLE       ON UPDATE NO ACTION ON DELETE NO ACTION     )     WITH (     OIDS=FALSE     );          CREATE TABLE camino_tipo     (     id serial NOT NULL,     nombre character varying(48) NOT NULL,     CONSTRAINT camino_tipo_pkey PRIMARY KEY (id)      )          CREATE TABLE camino     (     id serial NOT NULL,     nombre character varying(40) NOT NULL,     tipo_id integer NOT NULL,     CONSTRAINT camino_pkey PRIMARY KEY (id),     CONSTRAINT camino_tipo_fkey FOREIGN KEY (tipo_id)       REFERENCES camino_tipo (id) MATCH SIMPLE       ON UPDATE NO ACTION ON DELETE NO ACTION,     CONSTRAINT camino_name_uq UNIQUE (nombre)     )     WITH (     OIDS=FALSE     );      Any help will be appreciated. Thanks,


title: How to add a geoJSON layer to leaflet map?

text: i am newbie in leaflet & its geojson features... I am trying to test leaflet's functionalities. In fact i have built a map and its layer:               var map = new L.map('map', {scrollWheelZoom:true, doubleClickZoom:true, boxZoom:true, zoomControl:true});     var osmTiles = new L.TileLayer('http://localhost/googleapi/images/image1/tiles/{z}_{x}_{y}.png', {     maxZoom: 7,     minZoom: 3     }).addTo(map);         and now i want to draw a second layer (geoJson) upon it with a marker, so following the tutorila (cloudmate):               var geojsonLayer = new L.GeoJSON();     var geojsonFeature = {     "type": "Feature",     "properties": {         "name": "Coors Field",         "amenity": "Baseball Stadium",         "popupContent": "This is where the Rockies play!"     },     "geometry": {         "type": "Point",         "coordinates": [83.06877, -156.35742]     }     };     geojsonLayer.addData(geojsonFeature);      I am testing in my local PC (have not installed node + built leaflet) so maybe i am missing something!!! Any help?


title: Loading features incrementally in a vector layer

text: I am new to OpenLayers. In one experiment I am doing I want to load, incrementally, features in a vector layer. In my setting, these features should be loaded from distinct OSM files (generated dynamically) in a HTTP server. I know how to create a vector layer from a single OSM file, I have done it with a code similar to this:               ...     map = new OpenLayers.Map('map_element', {         maxResolution: 156543.0399,         numZoomLevels: 16,         units: 'm',          projection: new OpenLayers.Projection('EPSG:900913'),         displayProjection: new OpenLayers.Projection('EPSG:4326')     });          vectorLayer = new OpenLayers.Layer.Vector("POIs", {         strategies: [new OpenLayers.Strategy.Fixed()],         protocol: new OpenLayers.Protocol.HTTP({         url: "my_dynamically_generated_osm_file.osm",         format: new OpenLayers.Format.OSM()         }),         projection: new OpenLayers.Projection("EPSG:4326")     });          map.addLayers([<other_layers_here>, vectorLayer]);     ...      However, I am wondering how I can create first the vector layer without features, and afterwards adding new features defined in different OSM files according to the requests of the user. I know vector layer objects have a method `addFeatures()`, so I just need to know how could I obtain these features from the distinct OSM files provided by my http server ?


title: The shapes of the USA map

text: I'm newbie in GIS. I want use the shapefiles of USA map for my Geoserver. Where can I download shapes for USA? I know about the TIGER\Line shapes but it is separate shapes of counties and thematic maps (e.g. american indian areas). I want something like OpenStreetMap. P.S. Sorry for my english )


title: Can you change the add-in filename without changing the Visual Studio project name?

text: My solution has several projects in it, one of which is an Add-in project, and I'd like to keep the naming convention of all of the projects intact but change the filename of the generated .esriaddin file. Can the generated .esriaddin filename be changed automatically without without changing the Visual Studio project name or do I need to make a post-build event to rename it? Thanks