title: How do I determine programmatically what values to add to the proj4js constructor for a CRS or SRS such that I can iteratively render proper overlays

text: I think this is primarily a GIS based question as far as for example:               var crs = L.CRS.proj4js('EPSG:2154',                         '+proj=lcc +lat_1=49 +lat_2=44 +lat_0=46.5 +lon_0=3 +x_0=700000 +y_0=6600000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs',                         transformation);      So if I read in a wms tile layer through the leaflet api to create an overlay on a wfs base map how do I decide how to adjust those lat/lon values? And what is ellps, towgs, and units?


title: ArcIMS .axl marker symbols in 3d

text: I am using ArcIMS 10. I am fine with the .axl file showing 2d symbology and fill color, but i wonder if it is possible to use **3d marker symbols** in the ArcXML file? Would i have to use raster marker symbol tag to get a 3d point like the 3d pin available in ESRI ArcGIS_Explorer.style file? Cannot find any documentation about this.


title: Geoserver WMS request raster clipped with field boundary

text: I have recently set up GeoServer instance and uploaded raster data and a set of field boundaries. I would like to be able to select a field boundary and have the server return the portion of the raster that falls within the polygon (not the bounding box). Is this possible with a WMS request or how can it be done? Any help would be appreciated!


title: Add vertices costs to postgresql routable network

text: In osm2po it's possible to create custom routers as the sample router TrafficSignalsRouter, but the costs seem to be added on vertices after the routing. Is it possible to add these costs on the edges so that I can export them into my postgresql database (1 minute for each traffic signal in this sample)? (by the way the vertices costs seem to be taken into account because it changes the generated path, but the final cost remains the same in the logs)


title: Using ArcGIS to gain information about species locations

text: I am totally new to ArcGIS so I'm really sorry if this a stupid question. I am a conservation biologist and I've got a GDB database, with a bunch of species ranges on it. The output I desire from this, is a list of countries where the species live. So a list of countries that each ranges covers. Can anyone tell me how I would go about getting this? I've been googling how to use ArcGIS guides and getting nowhere...


title: Custome QGIS application with WPS and GDAL/OGR

text: I have basic question on how to start a project with QGIS. There are two kinds of users, and they are supposed to sign in (Login in) before they start using so that depending on the user type the QGIS interface and plug-ins will be presented to them. So basically there are two different QGIS interfaces. So my questions are:   * Is there any way I can add a login interface before QGIS starts?   * Can we customize the QGIS in different ways so that different users can have access to different QGIS interface and plug-in? After login the first user is supposed to do editing, and can save it on a local disk to continue later. But after some certain button is clicked, the process continue to the second user whose main job is to approve the edited features done by the second user. If approved, the newly edited data is committed to a database. If not the it will be reverted back to the first user with some comments.   * Is there any way that such a process can be implemented in QGIS? Would WPS be any help?   * If there is a need in the process to change one format to another, could we use the GDAL/OGR library?  I heard and most examples that I come across show how to convert formats with GDAL/OGR using command line.   * Is there any possible way you could use the API in your code? Like with in your WPS code or inside QGIS? I am guessing some of my questions may not give sense since I am all new to open source based GIS programming, but I would really be grateful for any kind of help. You be cool!


title: Convert shapefile lines to a polygon in GRASS

text: I have a shapefile of a map that is the outline of a number of islands. I would like to convert the lines into polygons (ie all islands filled in) using GRASS. I searched on how to do this with GRASS but just got some very brief and vague suggestions ... but I gave it a go anyway ....               use v.in.ogr to import the shapefile into grass     And view the result (I'm using QGIS to do this): Looks good.     v.clean input=g_hwp_p1 output=g_hwp_p1_clean tool=break,snap thresh=10,10     v.type input=g_hwp_p1_clean output=g_hwp_p1_clean_type type=line,boundary     v.centroids input=g_hwp_p1_clean_type output=g_hwp_p1_clean_type_centroids       This didnt work which is not surprising as I dont really understand what I'm doing. V.info says that...               g_hwp_p1 has 892 lines     g_hwp_p1_clean has 892 lines     g_hwp_p1_clean_type has 892 boundaries     g_hwp_p1_clean_type_centroids has 892 boundaries (0 centroids, 0 areas, 0 islands)          Also v.type reports: 892 incorrect boundaries.     And v.centroids reports: 0 features modified.     The final output looks much like the input (just an outline).      Can anyone spot what I'm doing wrong? Does anyone know what the correct parameters should be? I also tried this in QGIS using the lines-to-polygon fTool but this didnt do a very good job of it (it diddnt fill in the whole islands but just seemed to add little polygons round the edge). Shapefile is here... https://docs.google.com/leaf?id=0B93lksnTC7_cYWIxNTlmNzEtYTcwYi00YmZlLTlhZmEtZDViYmE1MmYxNDk2&hl=en_US Screenshoy of hwp_p1 .... https://docs.google.com/leaf?id=0B93lksnTC7_cN2Q5YTNiMjEtZjY4OS00ZTYwLWIxMDAtYzBlNmY1M2U1MGRl&hl=en_US And a close-up of a small Island - it forms a closed loop.. https://docs.google.com/leaf?id=0B93lksnTC7_cZjQ1ZDU4ZGQtMWNmMC00ZjNiLWFiZmEtM2MyMDE0MmYzODQ5&hl=en_US Thanks.


title: PostGIS two "circles" intersection

text: I am trying to see if two "circles" intersect. The data I have is two points and two radii (in meters). I am using srid 4326. My main issue (i think) is getting the radius in meters.               SELECT (ST_Buffer('POINT (-104.949517 39.740264)',300) &&                ST_Buffer('POINT (-104.999924 35.000089)',300))      My questions are   1. Is there a more clever way to do this?   2. Or how do I convert the radius from dd to meters to get correct buffers?


title: ESRI Javascript API: Return to the Map Extent you were at when you use browser back button

text: Can anyone give me some pointers (code would be helpful) on how to go back to the extent a user was at after they have navigate away and clicked the back button? So I'm viewing a map at X extent I click on some link and page changes, on clicking the back browser button I want to be able to get back to the Map Extent I was just before the page changed. Any help much appreciated.


title: Labeling classes for small roads on large extent maps in arcmap with sql statement

text: I have many small sub divisions that won't label nice, (lot of overlapping) unreadable. While using the bullet leader symbol. What would be the best statement that would place the labels in a clear spot near the sub division that would point to the road and also look clear and aligned horizontal?


title: What are the rules for releasing ArcObjects from memory in .NET?

text: Do I need to release each object? Is some memory management handled for me? ArcEngine 10, VS2010 C#


title: Combining Two Adjacent Raster Datasets

text: I am trying to figure out how to combine two adjacent rasters into one raster using ArcMAP 9.3.1. The problem I am having is that there appears to be a gap in the data. I would like to be able to combine the Rasters and/or Contour Data created from the Raster Datasets into one so that I can get rid of the gap in information. I am eveutually going to do a watershed ananalysis, so I am seeing the gap in information as a problem. I'd Post a Screenshot but I don't have enough reputation points yet. Thanks in Advance!


title: Deleted layers still appear while loading project QGIS 1.8

text: I'm working on a project where I digitized a few maps (georeferenced) and I also added a couple of csv files. I deleted the raster and cvs layers after I included them into the project as shapefiles but when I re-load the project to continue working on it these phantom layers still appear in the layer panel until the project is actually loaded and ready to work.   It doesn't make any problem other than these layers being listed as options to choose, for instance doing interpolation or other spatial options. The problem occurs when I delete the original csv file from the original location, even though it is not visible in the layer panel. If I do that, QGIS closes with error immediately after opening. I cannot move the csv or raster files from the folder they were originally saved because I cannot load the project. Same thing was discussed in qgis.org but the discussion was closed apparently because it was solved, which is not my case. Any clues? Windows Vista - QGIS Lisboa 1.8.0


title: How do I change from lat/long to a metric coordinate system in GRASS?

text: Yes, it's yet another transformation question. I am working with ASTER DEM data, which is LatLong WGS84. I want to process the DEM map into an STL file for use in a CFD simulation. The script I am using to accomplish this requires metric coordinates, though. Here's what my ASTER DEM mosaic looks like: ![enter image description here](http://i.stack.imgur.com/25Z02.png) I have created a new location in GRASS and chosen UTM with the etrs89 datum and grs80 spheroid. (I don't know if this is the correct choice.) Following these instructions, I created a mosaic of some ASTER tiles using `gdalwarp`:               gdalwarp -t_srs "+proj=utm +zone=32 +datum=etrs89" -r bilinear -dstnodata -9999 -co "TILED=YES" AST_DEM.vrt AST_DEM_Rhein_Main.tif      and managed to successfully import them. (EDIT: It would appear that proj4 does not support ETRS89 as a datum:               $ proj -ld     __datum_id__ __ellipse___ __definition/comments______________________________            WGS84 WGS84        towgs84=0,0,0           GGRS87 GRS80        towgs84=-199.87,74.79,246.62                               Greek_Geodetic_Reference_System_1987            NAD83 GRS80        towgs84=0,0,0                               North_American_Datum_1983            NAD27 clrk66       nadgrids=conus,ntv1_can.dat                               North_American_Datum_1927          potsdam bessel       towgs84=606.0,23.0,413.0                               Potsdam Rauenberg 1950 DHDN         carthage clark80      towgs84=-263.0,6.0,431.0                               Carthage 1934 Tunisia     hermannskogel bessel       towgs84=653.0,-212.0,449.0                               Hermannskogel            ire65 modif_airy   towgs84=482.530,-130.596,564.557,-1.042,-0.214,-0.631,8.15                               Ireland 1965      Apparently `proj` doesn't fail with an error if it doesn't recognize a parameter, it just silently proceeds without warning the user, even if the result is inaccurate. Caveat emptor!) I have a couple of problems with the result. First, when I imported the resulting mosaic, GRASS complained that it didn't recognized the datum:               r.in.gdal input=/home/sfbosch/Documents/Quelle/OpenFOAM/GIS/ASTER/Data/AST_DEM_Rhein_Main.tif output=AST_DEM_Rhein_Main --overwrite     WARNING: Datum <unknown> not recognised by GRASS and no parameters found     Projection of input dataset and current location appear to match     r.in.gdal complete. Raster map <AST_DEM_Rhein_Main> created.      Next, the mosaic obviously looks quite, well, warped: ![UTM](http://i.stack.imgur.com/9F8XR.png) I expected some change with the reprojection, but nothing that stark; could this be due to the non-existent (in proj4) datum? I am having trouble recognizing the topographical landmarks I had picked out for myself. Finally, I had been using OpenStreetMap data, and now it won't import because the projections don't match:               ERROR: Projection of dataset does not appear to match current location.          GRASS LOCATION PROJ_INFO is:     name: Universal Transverse Mercator     proj: utm     zone: 32     no_defs: defined     datum: etrs89     ellps: grs80     towgs84: 0.000,0.000,0.000          Import dataset PROJ_INFO is:     name: Lat/Lon     proj: ll     datum: wgs84     ellps: wgs84     no_defs: defined          You can use the -o flag to v.in.ogr to override this projection check.     Consider generating a new location with 'location' parameter from input data set.      I tried it with the override, but ended up with a layer that wasn't visible. I've since learned that I also need to reproject the vector layers with v.proj. Really, all I want is to work with **metric coordinates** ; it almost doesn't matter which system, the topography is the important part. Was there something wrong with my `gdalwarp` approach? Am I making this too complicated? What's the best way to accomplish what I want?


title: Circular insets in ArcGIS

text: Looking for a way to add an AutoCad-type circular inset in ArcGIS 10.1. Have tried many different methods but the best we have been able to do is shown here ![here](http://i.stack.imgur.com/BBoen.jpg) using "insert picture" after various image/raster processes.


title: Creating an "average" polygon

text: I found this thread in the old ESRI forums with the same type of problem I am facing today. The solution to the problem is explained in little detail and since I am not an expert, I was wondering if somebody could describe it to me in little more detail. We are doing a mental map analysis. We have let's say 50 polygons of a certain city square which are all different in shape, size and location. Our goal is to create a single polygon which would represent the answers from our whole input data. So basically, we need a method for generating a single polygon out of these 50 which would be an "average" polygon in size, shape and location. I hope I made my point and managed to explain what we need to do. If somebody would be kind enough to explain me the process it would be much appreciated. Thank you!


title: What is a good way to cluster point data by item counts in ArcGIS?

text: I have point features that have as an attribute an asset count. For instance, point A may have 3 assets whereas point B may have 300 assets. My goal is to have each point be in two categories: clustered or non- clustered. A clustered point means that, in a buffer with, say, a 5 mile radius, has 100 assets total between all points in that buffer. A non- clustered point is one that doesn't fall in that definition. The points could only be counted once. It is also advantageous to have more clustered than non- clustered points. One solution I had to this problem was to buffer every point and then for every buffer count the number of assets. If the number of assets were below 100, I would call the point associated with that buffer (the center point) non-clustered. Otherwise, it would be clustered. Then, I would go through the clustered buffers and start with the buffer with the most assets. Every point within the buffer would be clustered. Then I'd keep going with the points that were not in any buffer previously checked and then go through that process again until there were no more buffers. The problem with that process is that it automatically assumes that the center for the buffers is a point in the data. It may be better to have buffers centered around points not in the data to get more clustered points. Is there a better way to go about this problem? Programming is not an issue for me. I plan on making this a geoprocessing service if possible, so it's going to use 10 SP5.


title: Python with GDAL and PostGIS as a replacement of ArcGIS

text: Few times ago, I worked with daily land surface temperature satellite images of an area with the size of Alaska. I had to process around one thousand images. I downloaded them, mosaiced, projected, clipped to shapefile boundaries, extracted means out of the clipped rasters, and plotted the means on a time-temperature graph. The time- temperature graph was the final output. I did all these in ArcGIS, and even though I created an ArcPy script which helped automating a part of the work, it was still painstaking. Now I want to build a program that does the same thing, but 100% automatically. I merely want to input to the program the coordinates of the corners in order to define the area of interest and get the final output generated which is the time-temperature graph. Do you think this is possible using Python with the GDAL library and PostGIS as a raster data store?


title: Can QGIS report/calculate areas/perimeters of individual polygon components within multipart shapefile?

text: Can QGIS report/calculate the areas/perimeters of individual polygon components within a multipart shapefile without using the "explode" vector tool to generate a new file? using QGIS 2.2.0 on Windows from OSGEO4W


title: Can I do full SQL expressions through the "WHERE" clause in UpdateCursor?

text: Run as a standalone python script for a shapefile in a toolbox on ArcMap. For example:               psql_url_pop = '''SELECT shapefile.NAME_# + shapefile.WIKI_URL FROM public.shapefile'''          rows = arcpy.UpdateCursor("C:\shapefiles\shapefile.shp", psql_url_pop)      etc... Essentially I am trying to concatenate the values of one field with those of another and decided to use '+' over CONCAT because I want all the values.


Getting the following error when running my app: _Error: Tried to register widget with id==legendDiv but that id is already registered_ Not sure how to get rid of it. My javascript for the legend is:               //Add the Legend           legendLayers.push({               layer: layer,               title: ''           });           dojo.connect(map, 'onLayersAddResult', function (results) {               var legend = new esri.dijit.Legend({                   map: map,                   layerInfos: legendLayers               }, "legendDiv");               legend.startup();           });           map.addLayers([layer]);                if (layer.loaded) {               buildLayerList(layer);           } else {               dojo.connect(layer, "onLoad", buildLayerList);           }      and my HTML is:

      Can anybody suggest a workaround or a fix to the error. I've done some searching and it looks like I might need to do a destroy? Not sure how to do that.


title: Leaflet - Tiles stopped loading. Can't figure out why

text: Suddenly my map stopped loading tiles and I went on the leaflet site and downloaded the latest version of everything. ...still I could only make London visible with its latitude and longitude. I have this following code and I can't figure out what I need to get it working. Can someone please point out what I am missing?                             <title>Copenhagen not showing :(</title>         <meta charset="utf-8" />         <meta name="viewport" content="width=device-width, initial-scale=1.0">         <link rel="stylesheet" href="http://cdn.leafletjs.com/leaflet-0.7.2/leaflet.css" />


title: Issues with building a VRT. from OS map tiles in QGIS

text: I am trying to speed up rendering of OS raster map tiles in QGIS. The tiles are 1:10k Tiffs - about 130 of them. When I attempt to build the VRT it says complete, then it shows up an error message: **Warning 6: gdalbuildvrt does not support rotated geo transforms.** The only settings I had on whilst creating the VRT were: resolution - highest Add to canvas after finishing As you can see from the image some of the tiles were successful. I didn't have a problem with creating a VRT for 1:50k or 1:250k map tiles. Any suggestions or advice on speeding up rendering of OS map tiles (such as pyramids) would be much appreciated. Thanks!


title: How do I insert new-line character into Addin Extension Description?

text: Apparently the text box in Esri's Extension manager doesn't wrap text (10.0 SP2). ![enter image description here](http://i.stack.imgur.com/DDy43.png) As a workaround, I'm trying to insert newline characters into the xml file, but haven't found anything that works. Here's the xml:               <Extensions>       <Extension id="AmberGIS_ExtensionMgrAddin_ExtMgrExtension" class="ExtMgrExtension" productName="Extension State Manager" showInExtensionDialog="true" autoLoad="true">         <Description>Lets you specify extensions that will be \n turned off automatically when document closes.</Description>       </Extension>     </Extensions>      **Does anyone know how to put a newline character into the add-in xml config file?**


title: Inconsistent Identify Task- Explain this bug?

text:
Through some `console.log` debugging I've located the exact line of code that is failing me but I'm not quite sure how to fix the problem. Basically, I have an identify function that queries the parcels layer on click and returns the specified field values in an infoWindow popup. When testing, I noticed that sometimes it will return "no information available" instead of the actual data. If you attempt enough times, it will eventually give you the data. I'd like to prevent this from happening and I should never see "no information available" as the parcels layer has no gaps. Here's the code:               function executeIdentifyTask(evt) {         if (map.infoWindow.count == 0 && map.__LOD.level > 14) {             console.log("debug 1");             var tms;             showLoader();             identifyParams.geometry = evt.mapPoint;             identifyParams.mapExtent = map.extent;             identifyTask.execute(identifyParams, function(results) {                 var features = [];                 var commonContent = '';                 var zoneContent = '';                 var munContent = '';                 var zipContent = '';                 var mCheck = 0;                 console.log("debug 2");                 dojo.forEach(results, function(result) {                     if (result.layerName === "Council Districts" && mCheck == 0) {                         console.log("2.1");                         commonContent += "
Council District Name: 	" + result.feature.attributes['Label'] + "

Council District Population: 	" + result.feature.attributes['Population'] + "
 ";                         mCheck++;                     } else if (result.layerName === "Parcels") {                         features.push(result);                     }                 });                 dojo.forEach(results, function(result) {                     if (result.layerName === "Zoning") {                         console.log("2.2");                         zoneContent += "
Zone: 	" + result.feature.attributes['ZoningCode'] + "
 ";                     }                 });                 dojo.forEach(results, function(result) {                     if (result.layerName === "Zipcodes") {                         console.log("2.3");                         zipContent += "
Zip Code: 	" + result.feature.attributes['ZIPCODE'] + "
 ";                     }                 });                 dojo.forEach(results, function(result) {                     if (result.layerName === "Municipalities") {                         console.log("2.4");                         munContent += "
Municipality: 	" + result.feature.attributes['NAME'] + "
 ";                     }                 });                      results = dojo.map(features, function(result) {                     console.log("2.5");                     var feature = result.feature;                     console.log("2.6");                     var content = '';                     console.log("llllllll 2.7");                     if (result.layerName === "Parcels") {                         console.log("2.8");                         content += '
';                         content += "
Owner Name:	" + result.feature.attributes['ONAME'] + "	";                         content += "
PIN:	" + result.feature.attributes['PIN'] + "
";                         content += "
TMS:	" + result.feature.attributes['TMS'] + "
";                         content += "
Lot:	" + result.feature.attributes['LOT'] + "
";                         content += "
Parcel Type:	" + result.feature.attributes['PARCELTYPE'] + "
";                         content += munContent;                         content += zipContent;                         content += "
District:	" + result.feature.attributes['DIST'] + "
";                         content += zoneContent;                         content += "
Legal Description	" + result.feature.attributes['GRM_LegalDesc'] + "
";                         content += "
Square Footage:	" + result.feature.attributes['SHAPE.area'] + "
";                         content += "
Acreage:	" + result.feature.attributes['RE_acres'] + "
";                         content += commonContent;                         content += "
Owner Address:	" + result.feature.attributes['MADD'] + "
";                         content += "
City:	" + result.feature.attributes['CITY'] + "
";                         content += "
State:	" + result.feature.attributes['STATE'] + "
";                         content += "
ZIP:	" + result.feature.attributes['Zip'] + "
";                         content += "
Map Sub Block:	" + result.feature.attributes['MAPSUBBLOCK'] + "
";                         content += "
Condo Label:	" + result.feature.attributes['CONDOLABEL'] + "
";                         content += "
Lots:	" + result.feature.attributes['RE_lots'] + "
";                         content += "
Deed Book:	" + result.feature.attributes['DeedBook'] + "
";                         content += "
Deed Page:	" + result.feature.attributes['DeedPage'] + "
";                         content += "
Date Sold:	" + result.feature.attributes['SaleDate'] + "
";                         content += "
";                         var template = new esri.InfoTemplate(result.feature.attributes['TMS'], content);                         feature.setInfoTemplate(template);                         tms = result.feature.attributes['TMS'];                         console.log("debug 3");                              if (onBase == 1) {                             console.log("3.1");                             var featPropsArray = new Array();                             var featProp = new Object();                             featProp.Key = "TMS";                             featProp.Value = tms;                             featPropsArray.push(featProp);                                  g_featPropsArray = featPropsArray;                             OB_GenerateHitListForFeature(g_JSON_OB_GIS_URL, g_featPropsArray, "Parcels", "search_results", null);                                  if (!$(this).hasClass("shown"))                                 moveLegend(0);                             $('#tabs').tabs('select', 1);                             hideLoader();                         }                         return feature;                     }                 });                 console.log("debug 4");                 map.infoWindow.setFeatures(results);                 console.log("debug 5");                 map.infoWindow.show(identifyParams.geometry);                 console.log("debug 6");                 hideLoader();             });         } else if (map.infoWindow.count == 0 && map.__LOD.level <= 14) {             console.log("Need to zoom in to identify.")         } else {             console.log("debug 7");             map.infoWindow.clearFeatures();             console.log("debug 8");             //map.infoWindow.destroy();             map.infoWindow.hide();             console.log("debug 9");             executeIdentifyTask(evt);             console.log("debug 10");         }     }      When all executes successfully, the entirety of the code is run. When I get "no information available" the **results** block of code is skipped entirely.


title: How to edit layers which are using via OSM Bright in Tilemill

text: How could i edit a map which is edited via `OSM Bright` in TileMill? Layers are in TM added via PostGis. OSM Bright it's simply built according to informations here. I need delete some roads also delete labels. I have deleted and edited .shp file in `ArcGis`from project folder of TileMill but after saving of the .shp file and starting TM there are no any changes. Thanks for help. * * * I think it could be resolve by editing of .shp but it's still without changes. Do you think that changes could be shown after "re-postGis" to pgsql.


title: How to attribute lines with Ids of touching polygons?

text: I have two shapefiles: one contains streets of a city and one contains blocks (i.e. polygons bounded by the streets). I would like to have the streets attributed with the id of the left and right polygon. Is there a way to do that in ArcGIS?


title: How to add SLD to a wms layer dynamically from javascript

text: How to add SLD style which is written in String format and stored in a 'var' variable of javascript and asign it to a WMS layer while creating it.


title: How to get labels to diplay in Google Earth after saving to KML?

text: When saving a QGIS layer as kml, is it possible to have labels be displayed when opened in Google Earth?


title: What is the 'Layer' box in 'Save Vector Layer' good for?

text: In QGIS, in the 'save vector layer as' window, under OGR creation options, what is the 'Layer' box for?


title: STATISTICS not allowed in the Attribute Table

text: I created a shapefile, calculated geometry then went to grab the stats, the statistics was ghosted, not allowing me to grab those numbers. I am using 10.1 ArcView. I deleted the field and created a new field in ArcCatalog but no luck, same results. Not sure why that is happening. Some help, please.


title: How do I make the scalebar accurate when I am working with a google satellite image in QGIS?

text: I start a project in QGIS. I set the CRS to anything (in this case I use EPSG=26918 a NAD83 datum). I also enable on the fly projection. Then I import a vector. In this case a kml, which I believe is a CRS of WGS84. Then I use openlayers and add a google satellite image. Everything looks good and lines up well. I then open a composer window and add a map and a scalebar. When I check to see if the scalebar is accurate I find that it is not. If I go back to the project and check the measurements with the measurement toolbar I find that the measurement is accurate when the ellipsoid button is checked and inaccurate when the ellipsoid button is not checked. The measurement that I get when the ellipoid button is not checked is the same as I get in the composer window when I use the scalebar to figure out a length. I guess I have a sense of why this happens. But I want to know how I can get an accurate scalebar in my composer window when I am using a google satellite image? Thanks, David


title: How to install shapely for python 2.6 (Mac) NOT 2.7?

text: I'm trying to use **Polygonizer** plugin inside QGIS, but it needs **Shapely** to work. I've installed Shapely (pip install Shapely) but that's for **python 2.7**. QGIS works with **python 2.6** (on Snow Leopard) and I can't find previous installers for it. Maybe I'm missing something? Thanks.


title: Is there any Map event fired when map is loaded completely when we device is oriented from portrait to landscape mode?

text: We are using ESRI ArcGIS JavaScript API in one of our Mobile web app. When map is loaded first time, Onload event is fired. But when we change orientation from Portrait to Landscape, OnLoad event does not fired. We need such event so that when map is completely loaded in Landscape mode we can perform some opearion. We know that there is Orientaion change event in javascript but we need event which should fire when Map is loaded completely in Landscape mode after orientation change event is fired?


title: How to create buffer zone in metres for OSGB1936 layer?

text: I am creating buffer zones in metres on top of UK Ordnance Survey map files that are in OSGB1936 which uses degrees. How do I draw the buffers in meters? Many thanks


title: QGIS on Mavericks

text: I have been trying to install QGIS on Mavericks and I have failed. Prior to installing QGIS, I installed GDAL using homebrew. I also installed Numpy, SciPy and pyqt using the detailed instructions provided at: http://www.thisisthegreenroom.com/2011/installing-python-numpy-scipy- matplotlib-and-ipython-on-lion/ I modified my .bash_profile so that when I run `echo $PYTHONPATH` I get `/usr/local/lib/python2.7/site-packages:` I tried first installing qgis running               brew install qgis --with-grass --with-postgis      and I get the error message               libtool: link: cannot determine absolute directory name of `@@HOMEBREW_CELLAR@@/libpng/1.5.17/lib'     make[3]: *** [libcairo.la] Error 1     make[2]: *** [install] Error 2     make[1]: *** [install-recursive] Error 1     make: *** [install] Error 2          READ THIS: https://github.com/Homebrew/homebrew/wiki/troubleshooting          These open issues may also help:         https://github.com/Homebrew/homebrew/issues/3526         https://github.com/Homebrew/homebrew/pull/23375         https://github.com/Homebrew/homebrew/issues/20522         https://github.com/Homebrew/homebrew/pull/23252         https://github.com/Homebrew/homebrew/issues/16361         https://github.com/Homebrew/homebrew/pull/18528      Then I also tried downloading the QGIS .dmg file from http://www.kyngchaos.com/software/qgis However, I get an error message each time I try opening Install QGIS.pkg saying that I need matplotlib that I already have! Is anyone else also having problems installing QGIS on Mavericks? Does anyone know how to solve it? Thank you


title: C# code to display result of spatial query on arcmap

text: the problem is i have to display those points on arcmap that lies inside the envelope.I am explaining what i have done so far. I have one point feature class say featureclass ,open it and used one ARCGIS Snippet i.e. GetActiveviewFromArcmap() and constructed one envelope as-:               Iactiveview pactiveview=GetActiveviewFromArcmap(m_application);     Ienvelope penvelope=new EnvelopeClass();     penvelope=pactiveview.extent;     penvelope.QueryCoord(out minx, out miny,out maxX,out maxY);          Ispatialfilter sp=New IspatialfilterClass();     sp.geometry=penvelope;     sp.geometryfield=featureclass.shapefieldname;     sp.spatialrel=esrispatialrelenum.esrispatialrelIntersects;                    Ifeaturecursor fc=featureclass.Search(sp,false);     Ifeature feature=featurecursor.nextfeature();          What i have to do is to display the output of this spatial query on arcmap.how i will retrieve the selectionset of spatialfilter query and how to diaplay it on arcmap.


title: How to get cell size of raster shown in ArcMap using ArcGIS Engine and C#?

text: ![Cell Size](http://i.stack.imgur.com/mNVJv.png) As the figure shows above, after classifying the raster, we can get the cell size by opening the raster Properties to its Source tab, and we can calculate the area by mutiplying the cell count and the cell size. How can I get the cell size using C# instead?


title: Are there water/wastewater utility GIS symbols for electronic markers?

text: I hesitate to ask for standard symbols, but is there a de facto standard symbol(s)? Or a popular set of symbols to use for electronic markers? Our utility will soon start marking with electronic ball markers similar to these: http://images.pennnet.com/articles/ww/thm/th_193172.jpg (from Water World Article 'Municipality Marks New Pipelines for Future Access').   ![enter image description here](http://i.stack.imgur.com/kSwq7.jpg) I found this one somewhere, but I can't remember where I grabbed it from:   ![alt text](http://i.stack.imgur.com/0boK6.png) Are there others?


title: How to change a text element's font size using python? (ArcMap 10)

text: I need my script to change the font size of a text element if the string is too long to fit in the title section of the map. I've looked into it using Arcpy, but haven't had any luck. Does anyone know how to do this?


title: How would one geocode historical US census records?

text: I want to geocode census records from the U.S. Here is one example from 1870. This particular census page is from the 6th election district, 8th ward of New York County. Here's one from the 1910 Census. This one mentions the county, township, supervisor district, and enumeration district. I can geocode the county, state. I'm wondering if there are any resources which include latitude/longitude for districts, wards, supervisor districts, and enumeration districts from the past to help me geocode census records. I did find NHGIS, but their data is quite sparse for anything more fine grained than the county level before 1940.


title: IGlobeGraphicsLayer not updating in separate thread (Threading Issue?)

text: I'm currently using ArcGlobe to create a Polygon on a globe. Using examples I've found through arcgis.com I can now allow the user to click on the globe and create a polygon on the map. This works as intended and the visualization on the map is seen. However, I'm now trying to do this in a new thread so I can create a polygon overlay and attempt to return a list of Coordinates from a single point of entry (Since there doesn't seem to be a TrackPolygon function like there is in ArcMap). When putting this code into it's own thread, the event handlers on the globe still work; however the `GlobeGraphicsLayer.AddElement` is no longer displaying this data on the globe. Any ideas? Suggestions?


title: Modelling of rent prices - which interpolation method should I choose?

text: I have national dataset of ~1,4 million households. There I have information about rent, size (number of rooms and m2) and some additional characteristics of each household. I'd like to use this data to create surface of rent prices for the whole country and use this information as a proxy for estimation of values of remaining ~1.5 million households that are owned or do not have rent information. Couple of questions here: Is such approach appropriate for the this kind of problem at all? Which method of interpolation would be most suitable to use here? Also, would it be possible to take information about, for instance the size of the household into account? I'm on ArcGIS 9.3 with ArcInfo license.


title: Which GIS support DXF and/or DWG?

text: Which GIS support DXF and/or DWG (preferably read&write)? Update: Are there any open source GIS supporting these or similar CAD files? Support of CAD files generally seems to be problematic. What are the biggest issues/problems concerning DXF/DWG support in GIS?


title: Why does PDF fail to print in ArcGIS10?

text: Can anyone tell me why ArcGIS10 will not print to my Adobe9 "printer" in a Windows7Pro environment. Multiprocessor system, huge memory. Makes a .ai but crashes when i try a "print". ArcGIS10 also crashes when I try and use the native export to pdf function. Arc is up to date. Thanks ahead of time to those who respond.


title: Wifi hotspot localization

text: Using WiFi 802.11 I'm trying to build WiFi hotspot fingerprints based solely on radio signal strength and mac id. Now I'm wondering: Would it be possible to know the position of hotspot neighbors in line of sight in order to map their position and build more accurate fingerprints? I found, surfing around, that some papers have been written about measuring Angle of Arrival (AoA) of signals, so does anyone know of a library or tool (c/c++/java... other) that determines the hotspot position based on this? thx


title: Problem reprojecting annotation in ArcGIS 10

text: I have a large annotation feature class that I want to reproject from decimal degrees to a specific iteration of Albers. Using the Project tool, it appears to work, all annotation features are present in the output, attributes look fine, but the features don't draw. I can select them in the data frame, but can't see them. What to do?


title: MapServer - Polar Area Being Excluded

text: I am experiencing some unexpected results when trying to use MapServer to create images with a polar projection. I am using the standard BlueMarble tif as the base layer (in wgs84 GCS) and am trying to reproject image into WGS 84 / UPS North (WKID 32661). The polar region seems to be getting 'clipped' out of the image. The attached image was created using the mapfile below along with the shp2img utility included with ms4w. Can anyone explain whats going on here? --Mapfile MAP               NAME BMNG     STATUS ON     SIZE 800 500     IMAGETYPE PNG     UNITS METERS     EXTENT -28567783.609 -28567757.697 32567757.697 32567783.609     IMAGECOLOR 255 255 255     CONFIG "PROJ_LIB" "C:/DevTools/ms4w/proj/nad"     PROJECTION         "proj=stere"         "lat_0=90"         "lat_ts=90"         "lon_0=0"         "k=0.994"         "x_0=2000000"         "y_0=2000000"         "ellps=WGS84"         "datum=WGS84"         "units=m"         "no_defs"     END          LAYER         NAME 'BMNG_UPSN'         TYPE RASTER         STATUS DEFAULT         DATA 'D:\stuff\bmng\bluemarble.tif'         PROJECTION             "init=epsg:4326"         END     END #LAYER      END #MAP ![enter image description here](http://i.stack.imgur.com/ZUrbO.png)


title: CartoDB SQL string as integer

text: Because data is imported as string, I need to query string values as if they were integers, is there something I could use to cast values to string? Where could I find documentation for this SQL in CartoDB? I found this


title: Are any fictional worlds available in some GIS format?

text: In a discussion with colleagues, this question arose: is there any fictional world (e.g. Middle-earth, Discworld etc.) which could be used in a GIS context, e.g. a MBTiles archive and a shapefile or any other format usable with OpenLayers / Polymaps / Leaflet?


title: Strange Behavior in Focal Statistics

text: I am trying to calculate the maximum value in a neighborhood on a raster using the Focal Statistics tool for Spatial Analyst in ArcGIS 10. Using python, I loop though several different radii, with each defining a circular neighborhood, as below:               bandwidth = [1609.344*.125,1609.344*.25,1609.344*.5,1609.344*1,1609.344*2,1609.344*4,1609.344*8,1609.344*16]     for bw in bandwidth:         kd = "copc_kd_06_" + str(count)         nhMax = FocalStatistics(kd, NbrCircle(bw, "MAP"), "MAXIMUM", "NODATA")      Note that the "1609.344" is the conversion from meters to miles, so I'm trying to find the neighborhood maximum at 0.125, 0.25, 0.5, 1.0... miles. The problem is that 0.125, 0.25, 0.5, and 2.0 perform well, while 1.0, 4.0, 8.0, and 16.0 produce extremely odd results, as shown below. ![Odd results in focal statistics](http://i.stack.imgur.com/7IDqz.png) So, specifically, instead of the circular neighborhoods given in the command, it produces irregular neighborhoods of circles split in three. Has anyone else seen this, or do you know of a workaround?


title: Openlayers WFS change featureType

text: I am trying to change wfs featuretype "on the fly" in openlayers. (This method http://dev.openlayers.org/docs/files/OpenLayers/Protocol/WFS/v1-js.html#OpenLayers.Protocol.WFS.v1.setFeatureType) This is how I instantiate WFS vector layer. It works fine, features are displayed.                $WFSProtocol = new OpenLayers.Protocol.WFS({         version: "1.1.0",         featurePrefix: 'geog585',         srsName: "EPSG:3857",         url:  "http://localhost:8080/geoserver/wfs",         featureType: "Neighborhoods",          });          $philadelphiaMarketsWFS = new OpenLayers.Layer.Vector("Vector",      {     strategies: [new OpenLayers.Strategy.BBOX()],     protocol: $WFSProtocol,          });     $map.addLayer($philadelphiaMarketsWFS);      Now I have set up a select control on this vector layer, when a feature is clicked I want to change the featureType attribute to something different               $selectControl = new OpenLayers.Control.SelectFeature($philadelphiaMarketsWFS,      {              hover: true,         clickFeature: function(feature)                 {                       $philadelphiaMarketsWFS.protocol.setFeatureType('FarmersMarkets');                       $WFSProtocol.setFeatureType('FarmersMarkets');                       $philadelphiaMarketsWFS.refresh({force: true});                      }                    });                    $map.addControl($selectControl);     $selectControl.activate();      But after the refresh no features are displayed(vector layer is empty) Could anyone help me out with this, could not find any examples on this method... Thanks


title: Pure html based indoor navigation

text: I would like to locate (position) people inside buildings without forcing them to download a native app. I'm focusing on iOS and Android devices. I therefore thought about wifi triangulation (trilateration), GPS repeater, GSM beamforming or even acoustic background spectrum. Is there any product/tool to achieve this task?


title: How to install MapServer on IIS 7.5 with windows 7 OS

text: I want to use the map server in my asp.net application. For that one, I need to install mapserver on IIS. I am using windows 7 with 64 bit os. I have IIS 7.5 installed on my PC. Can any one help me to on install mapserver on IIS 7.5 with windows 7 OS.


title: How to crop a Raster by Region?

text: I have created a grid using the point density function in VM. I would like to crop this image to a region created from buffering a polyline feature set. Is this possible? Alternatively I have used the set clip region tool, however this only hides the rest of the image where I would like to save the cropped region as a new layer/image. Many thanks!


title: How to change the default color used by openlayers hover?

text: It is possible to change the color of highlightCtrl when I use onselectedHover? Default is light blue. I want to change it to another one eg. red? Does anyone know where this value is defined?


title: GeoServer - Caching with Response Headers Settings

text: I'm trying to find a method to cache one layer that its data is updated every 2 hours in the DB. From what I've found the best option is using Response Cache Headers timeout http://docs.geoserver.org/stable/en/user/webadmin/data/layers.html What I've done is:   1. enable direct integration with GeoServer WMS.   2. ass tiled=true parameter to WMS request   3. replace `http://GEOSERVER_URL/wms` with `http://GEOSERVER_URL/gwc/service/wms`   4. check "Response Cache Headers" with 7200 time specified as in http://docs.geoserver.org/stable/en/user/webadmin/data/layers.html and when i look at the headers I see:               Cache-Control:max-age=7200, must-revalidate     Content-Disposition:inline; filename=geoserver-dispatch.image     Content-Type:image/png     Date:Sun, 04 May 2014 18:38:16 GMT     Expires:Sun, 04 May 2014 20:36:36 GMT     geowebcache-cache-result:HIT     geowebcache-crs:EPSG:900913     geowebcache-gridset:EPSG:900913     geowebcache-tile-bounds:12523442.712499999,2504688.5425000004,15028131.254999999,5009377.085000001     geowebcache-tile-index:[13, 9, 4]     geowebcache-tile-index:[13, 9, 4]     Last-Modified:Sun, 04 May 2014 07:06:05 GMT     Server:Apache-Coyote/1.1     Transfer-Encoding:chunked      I thought it worked but after 2 hours the layer didn't get updated! it's like that the layer is cached and frozen like forever. I think I'm missing something. any idea?


title: How to creat multiple rings between two circular buffers?

text: I am working with mathematical-ecological modeling and I am also a new user of QGIS Software. First I am creating 3 circular buffers with same center (shape center) covering a region with roads (shape roads), so I have initialy two shapes and after I create 3 circular buffers. All mathematical work I use R software to produce the number of intermediary rings and distances with each these intermediary rings that will be between the pair of main circular buffers. My firsts questions are: How could I create automated the process of rings or buffers between the main circulars buffers? Are there scripts or other solution for this? My second and last: I have roads crossing the buffers and in case of creation of rings (or intermediary buffers) how could I calculate the length of the roads between the two intermediary rings? **UPDATE** ![enter image description here](http://i.stack.imgur.com/KfHM1.jpg) The problem is the rings made with a non automate process and also made in atribute table. So I wish make together the automate process with the calculation of the roads that cross the pair of rings. Using R like you said it is possible, but I would like to know the both process because it will be better to know how will be more easy for me. it is my first work with QGIS.


title: Find crossing lines in ArcGIS

text: I have a network of lines (pipes) and I want to be able to show the difference between an intersecting line and a crossing line. I was wondering if there is a tool for generating points on all crossings in ArcGIS. If I use "Intersect" I can find all intersecting and crossing points. I only need the crossing points.


title: How do you open/display and extract NetCDF files info?

text: I am working with NetCDF files using R & Matlab and I would like to know if the GIS community is using them. I know that ArcGIS can open & display, however most of the times the results are not good, since the structure of the netcdf files is always different. I would like to know what kind of software/add-ons are you using to open & display netcdf files? do you extract any geographic+time series info?


title: What's this line in Landsat MSS image?

text: I have a Landsat MSS image file that contains a vertical string of incorrect pixels. Does anybody know what this might be? ![enter image description here](http://i.stack.imgur.com/uVG4t.jpg)


title: Python Search cursor: Error handling

text: I have a search cursor that using a SQL query to match addresses and then returns the associated roll number with that address. However, some addresses do not match the roll number and then I believe the search cursor stops. Is there a bit of code I can throw in that will continue with the next address even if a match isn't found? I think an 'if' statement would work well, but I don't know what to write for it! I'm having trouble generating this. My code is iterating through foldernames to get the address and querying a parcel layer. I can post the code later, if you need an example.               import arcpy     import os     import re     import sys     import traceback     import collections          basedir = r"C:\Test"     fc = r"G:\ParcelsFixNames.shp"     field1 = "ADD1"     field2 = "STREETNAME"               try:     #find all foldernames in directory         for fn in os.listdir(basedir):             string = str(fn)     #split name at hyphen              if re.findall('-', fn):     #right side of hyphen is street address, leftside is street name                 lhs, rhs =  string.split('-')                 street = str.lstrip(lhs)                 street = str.replace(street,'.','')                 street = str.rstrip(street)                 street = str.upper(street)                 add = str.lstrip(rhs)                 expression3 = '"{field1}" = \'{add}\' AND "{field2}" LIKE \'%{street}%\''.format(field1=field1, add=add, field2=field2, street=street)                 rows = arcpy.SearchCursor(fc, fields="STREETNAME; ARN; ADD1", where_clause=expression3)                 for row in rows:                      Roll = row.getValue("ARN")                      newname = str(Roll)                      newpath = os.path.join(basedir, newname)                      #if the name exists already add _1 to end                      if os.path.exists(newpath):                         os.rename(os.path.join(basedir, fn),                             os.path.join(basedir, newname + "_" + "1"))                      else:                         os.rename(os.path.join(basedir, fn),                             os.path.join(basedir, newname))             else:                 print add + " " + street + " " + "Already exists...passing"                 continue     except:         print ("Error occurred")         tb = sys.exc_info()[2]         tbinfo = traceback.format_tb(tb)[0]         pymsg = "PYTHON ERRORS:\nTraceback Info:\n" + tbinfo + "\nError Info:\n    " + \             str(sys.exc_type)+ ": " + str(sys.exc_value) + "\n"         msgs = "GP ERRORS:\n" + arcpy.GetMessages(2 )+ "\n"         print (pymsg)         print (msgs)


title: Why are map tiles being repeated? Using MapServer and Openlayers

text: Problem I have: map tiles from my layer are being repeated. Shapefile is projection 3059. I have specified 3059 is in the mapfile and also when creating map object. I know OpenLayers.Layer.WMS is recommended, but when I use it I only see empty map. When I use OpenLayers.Layer.MapServer, map is displayed (and tiles are repeated), so I'm using that for now. **My layer (EPSG:3059) looks like this** ![My layer](http://i.stack.imgur.com/jYWfg.png) Versions:   Openlayers 2.13.1   MapServer 6.0.3 **WMS Server mapfile:**               MAP         NAME country_map         STATUS ON         SIZE 800 500         EXTENT 313187.4 146161.56 766469.92 470367.64         UNITS METERS         SHAPEPATH "shapefiles"         IMAGECOLOR 155 155 155              WEB             TEMPLATE "template.html"             IMAGEPATH "C:/MS4W/ms4w/apache/htdocs/gis/output/"             IMAGEURL "/gis/output/"                  METADATA                 "wms_title"             "MapServer + Openlayers"                 "wms_onlineresource"    "http://localhost:1025/cgi-bin/mapserv.exe?map=C:\MS4W\ms4w\Apache\htdocs\gis\mapfile2.map"                 "wms_enable_request"    "*"                 "wms_srs"               "EPSG:3059"             END         END              PROJECTION             "init=epsg:3059"         END          # example layer         LAYER         METADATA             "wms_title"             "lat_rob_region"         END         NAME "lat_rob_region"         DATA lat_rob_region         STATUS ON         TYPE POLYGON             CLASS                 NAME "lat_rob_region"                 STYLE                     COLOR 0 255 0                 END             END         END     END      **Javascript:**               var lon = 25;     var lat = 57;     var zoom = 5;     var map, layer1;          function init()     {         map = new OpenLayers.Map( 'map',{displayProjection: new OpenLayers.Projection("EPSG:3059")} );              layer1 = new OpenLayers.Layer.MapServer( "Country map",                  "http://localhost:1025/cgi-bin/mapserv.exe?map=C:/MS4W/ms4w/Apache/htdocs/gis/mapfile2.map", {layers: 'lat_rob_region'},                 {gutter: 15},{srs: 'EPSG:3059'});         map.addLayer(layer1);              map.setCenter(new OpenLayers.LonLat(lon, lat), zoom);         map.addControl( new OpenLayers.Control.LayerSwitcher() );     }


title: How to extract information from a tiff file

text: I have a tiff file that I loaded into ArcGIS. It has roads, rivers and other features including the 600m contour for the region. I want to create a file that only shows the 600m contour. How can I do that?


title: Is there a way to calculate closest distances to a point without using Near Analysis?

text: I frequently get asked to calculate distances from a specific point (A drilling well in the middle of the sea for example) to a number of Environmental layers (Protected Areas, Pipelines, Other wells etc, so would include points, lines and polygons). I currently do it manually which takes time so I thought there could be a way to automate this using python. It seems the best way to do this is use Near Analysis, which I do not have licence for. I was just wondering is there anyway to use python to still generate some sort of table that would give me the closest distances from my input feature to my numerous environmental features, the names of said closest features, and the distance (and also direction if possible, but thats not essential). Any ideas?


title: Is there a maximum length of pdf file QGIS composer can plot?

text: In QGIS 2.4 is there a maximum dimension in milimeters composer can produce while exporting to pdf? Trying to produce a 3000 x 900 mm or 2800 x 900 mm map. The QGIS composer displays all the map area but everytime i export it to pdf, the map gets cut at about 2780 mm displaying white space.


title: Translating hdop, pdop and vdop to metric accuracy from given NMEA strings

text: How can I translate dilution of precision to a metric unit? For example: say I got 1.2 in hdop from NMEA sentence of my GPS receiver, in what should I multiply to get a metric accuracy? Just like the way Android or other platforms do. It may be a duplicate (This link), but it was said: " **All the GPS receiver is doing, is multiplying the dilution of precision figure by the resolution of the device.** " I didn't catch that, what is the "resolution of the device"?


title: How to extract features with specific coordinates from a ShapeFile?

text: I would like to extract a number of features from a ShapeFile which match with **specific latitude/longitude values**. I am interested if/how this can be done with QGIS and ogr2ogr. The output should be a ShapeFile as well; for QGIS a layer is sufficient. I assume the command for ogr2ogr could be something like:               ogr2ogr -f "ESRI Shapefile" -where "Y = 40.752459 AND X = -73.983521" output.shp input.shp      However I am having difficulties finding out what the attributes for lat/lon are named. _Optional_ : Also, is there an option for a fuzzy query where I can configure a radius around the coordinates? I am running Ubuntu 12.04. with QGIS 1.9.0 via ubuntugis-nightly and GDAL 1.10.0.


title: How to use SRTM DEM files for shaded relief in ArcMap?

text: I would like to create a nice shaded relief for France using SRTM3 data. Most of the threads on the subject are from awhile ago, when ArcMap was not able to import these without changing the file extension. Now, it seems, ArcMap can import the .hgt files just fine. The problem is that the value for elevation is set at 32767 to - 32767. If I change the symbology settings and set the stretch value to standard deviation. This sets the values much lower and you can see the elevation differences. ![Image showing difference](http://i.imgur.com/l3XzKO8.jpg) A complete set of images is present at this **imgur album** showing the comparison, stretch setting, and final difference. This creates two more problems: - How might I do this on a mass scale? I've searched for tools and unfortunately most conversion custom toolkits seem to be inaccessible. - Each hgt file seems to be shaded differently, due to the changing of the stretch value to standard deviation, as seen below: I apologize if this is something that has been covered before. I am still relatively new to the world of Cartography and GIS so there may be a simple tool that I should be using. Many thanks for any guidance on the subject.


title: How to add automate values to attribute field?

text: I have parcel, Buildings layers, without any information attribute table, and i want to replace "NULL" values with numbers form 1 to 2000, in "Parcel ID", or "Building ID" field thank you.


title: Vba Cons Numbers

text: I want to make a list of cons numbers(ex. 1,2,3...) in field calculator (VBA). When I am putting for example: (dim d d=[field]+1). The result is all my fields are =1. How to solve this arc map (vba). Thanks.


title: Add sequence numbering based on attributes in a table

text: Is there a way to auto populate a sequence of numbers based on another attribute in a table? For example, if field "Name" has a sequence of numbers from 1 to 10, and I want to number all the rows in the table that associated to "Name = 1" from 1 to 10, and then start this process over when "Name = 2" etc (see attached graphic). Can this be done within the Field Calculator, with a plugin, or does it need to be coded? Thanks...![enter image description here](http://i.stack.imgur.com/4b06f.jpg)


title: Use point shape file to spatial query multi polygon shape files with python?

text: I have an address point layer and would like to spatially query this layer using various polygon shape files to intersect this information so that each address contains the attributes from these polygon layers. Is this possible using non ESRI python? I have used arcpy but don't know if this kind of thing is possible with open source python options? If so is the python script fairly simple? Possible examples?


title: Join Table to QueryFeatureClass

text: I've joined a table to a SQL-QueryLayer like posted here Using IRelQueryTable When I set `geoFeatureLayer.FeatureClass = relQueryTable as IFeatureClass;` the Layer in ArcMap forgot that it is an Query-FeatureClass. So in the Source- Tab of the Layer-Properties-Dialog there is no entry for the DataType and in the Fields-Tab the Fields are listed full-qualified in the Alias-Names. But the Join seems to work because it is correctly listed. Any idea?


title: using mapnik xml file to render tiles in python

text: Is it possible to use an xml file in python? without using the bindings for python for mapnik. basically i'm using django/geodjango, postgis and well i decided i'd try and design the map on tilemill then call the exported xml file. I've tried to look it up. but sadly haven't found anything.


title: Associating location data from PostGIS for a Choropleth map?

text: I'm trying to build my first choropleth map from data that we have stored in our PostGIS enabled Postgres database. So far I've been able to build a map on Mapbox using a topojson set of UK postal area data, but I'm stumped on how to show concentration of our data points per postal area from our database. Basically I'm trying to show we have a concentration of user type X in postal areas. I know enough that querying the db point-in-polygon for each polygon in the postal data set would be slow and cumbersome, but is there a way to search the data by postal area and return a merged data set e.g. x users in area 1, y users in area 2? Should I import the postal areas into a spatial table to do the query against in PostGIS? Currently all the data points are output as JSON but I'm not sure how to associated the relevant points with each postal area.


title: How to script a csv to file geodatabase tool?

text: I have many separate field data files in CSV format and I wish to automatically convert these CSV files into a geodatabase to be able to load them into a schema I created for them in my new ArcSDE enterprise database. The SDE version is 10.1. I am able to do this manually in ArcGIS but I have over 500 files to complete and so I do need a script will automate it for me. I would be glad to have the output in a form of a feature class. I do not have an existing output dataset besides the ones I did manually, it needs to be created from scratch using the CSV files as input dataset. I appreciate any help I can get from here and thank you very much for your consideration and time.


title: GML TimePeriod mapping in MapForce

text: I made SQLServer database structure based on xsd schema I got. In that schema two fields have gml:TimePeriodType. In database I put two columns for each field of this type (begin and end) and the datatye is datetime. I use FME to translate database into GML, and then MapForce to map that GML to xsd schema. I supposed that I can map begin and end column to some atributes of TimePeriod columns in my xsd schema. But (guess), I cannt do it. Can some help me what I need to do and where I am wrong? FME is translating this columns into string (e.g. 20140522000000.000), which is wierd to me (I think I need 2014-05-22T00:00:00 format). What I am missing? P.S. Sorry if question is stupid, but I am new to all of this, and I need pointers for going forward. Thanks!


title: How to remove terracing from Civil3D GeoTIFF export?

text: I've received a GeoTIFF DEM that was exported from `Civil3D` and at every contour there is a distinct terrace created when I process an aspect or slope in `QGIS` with the GRASS plugins. ![terracing](http://i.stack.imgur.com/baXts.jpg) I've looked around a bit and I think this is related to the algorithm used to create the DEM on the `Civil3D` side (Inverse Distance Weighting method). Are there parameters within `Civil3D` to prevent this? or Is there a way to process the terraces out after the file is created in `QGIS`?


title: Chaining Rendering Transformations: Making isotherms using Barnes Surface + Contour process?

text: I'm trying to calculate isotherms following the istructions showed in the post on OSGeo blog (http://blog.opengeo.org/2013/01/10/chaining-rendering- transformations-in-geoserver/) I'm using GeoServer 2.2.3 on a Windows 7 machine (8 Gbyte RAM, biprocessor) I've a shapefile with points and temperature and I'm trying to use the Meteo-1.sld Using the layer in preview http://:8080/geoserver/CesareWorkspace/wms?service=WMS&version=1.1.0&request=GetMap&layers=CesareWorkspace:CentralineTemperature-01-01-2001&styles=&bbox=318287.0,4880823.0,515945.0,5142603.0&width=386&height=512&srs=EPSG:32632&format=application/openlayers it works fine but when I try to use le Meteo-1.sld http://:8080/geoserver/CesareWorkspace/wms?service=WMS&version=1.1.0&request=GetMap&layers=CesareWorkspace:CentralineTemperature-01-01-2001&styles=&bbox=318287.0,4880823.0,515945.0,5142603.0&width=386&height=512&srs=EPSG:32632&format=application/openlayers&STYLES=Meteo-1 the map is white. It seems that there are no errors in the log. When I try to use the Barnes Surface WPS directy (gs.BarnesSurface), I obtain a black tiff .... At the following link https://www.dropbox.com/sh/fmbdhzwp5mobpcq/lf43d3SykM there are:   * data.zip --> the shapefile I'm using    * Meteo-1.sld --> the sld I'm using    * geoserver.log --> geoserver log Could someone tell me which configurations I've to use in the sld for my data? Or give me a sample data and a sld that works? Thank you in advance!!


title: How to create an empty group layer within an existing group layer with arcpy.mapping?

text: Referencing this question, Add a new group layer with python ArcGIS 10, I'm trying to use python to add a empty group within a group layer. So far I'm able to add a empty group within the map document, however I'm getting an error when trying to add another empty group within the existing group layer. Here is the error message: > Runtime error Traceback (most recent call last): File "", line 1, in File > "c:\program files (x86)\arcgis\desktop10.1\arcpy\arcpy\utils.py", line 181, > in fn_ return fn(*args, **kw) File "c:\program files > (x86)\arcgis\desktop10.1\arcpy\arcpy\mapping.py", line 88, in > AddLayerToGroup assert isinstance(target_group_layer, Layer) and > target_group_layer._arc_object.isGroupLayer AssertionError Here is my code:                       mxd = arcpy.mapping.MapDocument(r"C:\Temp\test.mxd")             df = arcpy.mapping.ListDataFrames(mxd, "Layers")[0]             groupLayer = arcpy.mapping.Layer(r"C:\Temp\Group.lyr")             arcpy.mapping.AddLayer(df, groupLayer, "BOTTOM")             targetGroupLayer = arcpy.mapping.ListLayers(mxd, "Group", df)             addLayer = arcpy.mapping.Layer(r"C:\Temp\Group2.lyr")             arcpy.mapping.AddLayerToGroup(df, targetGroupLayer, addLayer, "BOTTOM") # error here             mxd.save()             del mxd,df      Can you add empty group layers within an existing group (using v10.1 SP1)? The ESRI documentation (see image below) states that you can add .lyr to groups (maybe not a .lyr group layer?). ![enter image description here](http://i.stack.imgur.com/cOunt.gif) Thanks **UPDATE** Here is the working code (thanks to Jason Scheirer):                   mxd = arcpy.mapping.MapDocument(r"C:\Temp\test.mxd")         df = arcpy.mapping.ListDataFrames(mxd, "Layers")[0]         groupLayer = arcpy.mapping.Layer(r"C:\Temp\Group.lyr")         arcpy.mapping.AddLayer(df, groupLayer, "BOTTOM")         targetGroupLayer = arcpy.mapping.ListLayers(mxd, "Group", df)         addLayer = arcpy.mapping.Layer(r"C:\Temp\Group2.lyr")         arcpy.mapping.AddLayerToGroup(df, targetGroupLayer, addLayer, "BOTTOM")[0]         mxd.save()         del mxd,df


title: Twitter live feed on ArcGIS desktop

text: I am a complete beginner with ArcGIS and i have been looking for a way to live feed tweets to ArcGIS desktop. A quick google search brings this blogpost up from 2009: http://blogs.esri.com/esri/arcgis/2009/12/15/arcgis-explorer-and- geo-service-integration-twitter-add-in/ However, the link to the add-on is dead. Does this add-on still exist? If so, does anyone know where. If not, does anyone know why?


title: ArcMAP path distance function handling of polygon and raster inputs

text: I am using a path distance function to calculate distances from a set of source cells on a land use raster. Unfortunately, my land use raster and my other inputs have different cell sizes. The land use raster has 300 m pixel sizes and my other inputs have ~90m pixels. I would prefer to avoid resampling. I would also like the entire area of the source pixel to be represented as 0 on the resulting distance raster. I will be pulling out individual pixels from the land use raster to use as source cells, and the easiest way to do this is to convert to points . However, given the size differences in the pixels, I'm not sure that it would be accurate to represent the source cells as points. I was thinking of using polygons instead (which can be expanded to the 300 m size of the pixel), but I'm unsure of how the tool will handle the polygons when calculating distance. If the tool will use the centroid of the polygon rather than the edges, then using a polygon is no improvement over using a point to represent the area. I had also considered using a raster as an input, but, again, I'm not sure how the tool will handle the size differences. **Essentially: Does anyone know how exactly ESRI's distance tools handle different inputs during execution when the source inputs may not match the analysis cell size or may span across multiple pixels in the analysis cell size? How does the program handle resampling and/or source cell definitions in these cases?**


title: Problem with GRASS in QGIS

text: I am trying to import a vectorlayer to GRASS from QGIS. I get an error and hope someone here can help. Been trying to fix it for hours, but come no closer to a solution..... This is the error I get: > Unable to create table: 'create table kulturminnepunkt (cat integer, > x2_0214K_ double precision, x2_0214K_I double precision, KOORDH double > precision, OBJTYPE varchar ( 32 ), BYGGTYP_NB integer, BYGGNR integer, > TILBYGGSNR integer, SEFRAKKOMM integer, REGISTRERI integer, HUSLOEPENR > integer, KM_BETEGN varchar ( 100 ), NAVN varchar ( 60 ), KM_VERNEVE varchar > ( 1 ), ETASJETALL varchar ( 3 ), MAALTALL double precision, LENGDEENHE > varchar ( 20 ), ANTALLSKOR integer, SEFRAK_FUN integer, SEFRAK_FUN varchar ( > 2 ), OPPRINNELI integer, FORHOLDAND varchar ( 2 ), FYSISKMILJ integer, > UNDERBYGNI integer, KJELLER integer, YTTERVEGG integer, FASADE integer, > TAKFORM integer, TAKTEKKING integer, SEFRAK_TIL varchar ( 30 ), TIDSANGIVE > integer, VERN_LOV varchar ( 3 ), VERN_PARA varchar ( 10 ), VERNTYPE varchar > ( 4 ), VERNEDATO date, ENDRINGSGR varchar ( 1 ), VURDERTDAT date, FELTREGIST > varchar ( 50 ), DATFGSTDAT varchar ( 20 ), AJOURFORTD date, AJOURFORTA > varchar ( 50 ), KM_DAT varchar ( 3 ), KM_DATKVAL varchar ( 1 ), MATRIKKELK > integer, MATRIKKEL1 integer, GNR integer, GNR1 integer, BNR integer, BNR1 > integer, FNR integer, FNR1 integer, SNR integer, SNR1 integer, KM_ANTALL > integer, KOMM varchar ( 4 ), KM_SYNLIG varchar ( 3 ), KM_KATEGOR varchar ( 5 > ), KM_HOVEDGR varchar ( 2 ), KM_FUNK_OP varchar ( 4 ), KM_FUNK_NA varchar ( > 4 ), KM_MAT varchar ( 2 ), PLANID varchar ( 16 ), PLANSTAT integer, AREALST > integer, LINK varchar ( 254 ), OMRADEID varchar ( 4 ), ORGDATVERT varchar ( > 50 ), KOPIDATO varchar ( 20 ), kulturid integer, Avstand double precision)' > Finished with error Any suggestions?


title: Geoserver 2.3.1 not finding all columns from PostGIS connection

text: so I am not sure if I have missed something or this a particular problem with this version of Geoserver. I want to style a point layer with a SLD, based on values in a particular column. The column is an enumerated type relating to a type of activity. i.e. type_of_activity: { camping, fishing, hunting } But when I reach the configuration part of publishing the layer the column I base my style on is not present in the feature type details, and when previewing the layer, an exception is caused stating that "type_of_activity" does not exist. Any ideas?


title: Criteria for cartographic capacity

text: For a paper at university i am trying to create a cataloge for cartographic capacity of different GIS programs, especially open source like QGIS and so on The aim of this is to eventually have a cataloge that says: if you want to produce a (just one example) street map, do not take program A because it is not possible to change lines into representations and thus the overlying of lines can not be manually permitted. so my Question is: **What are criteria for cartographic capabilities? Or maybe it is easier to say: by which criteria can i define how good an openGIS program is for producing maps?** So far i have tried to think about a few criteria and came to these: -what kind and amount of different labels can be created -can classes be created that are assigned to different labels -symbols available -production of a map sheet -what can be done with the legend I hope these are reasonable but please tell me anyways so far i have only worked with ArcGIS and thus only know about some of its capabilities. Any answers are highly appreciated **Edit:** Okay i´ll try to put it straight for you: I want to find criteria that validate the capabilities of _openGIS_ to produce maps. What do you think are such criteria? What options are required in a GIS program to produce efficient and clearly represented maps? This is as Barret correctly mentioned NOT about knowledge of a certain user, but solely about the capabilities of the program! The above criteria are only a small part of that, so to make clear what i mean.


title: arcgis server 10, vb.net, non pooled service, web adf

text: i am running arcgis server 10. i created a A.MXD file, then started arcgis server and successfully added non-pooled service using A.MXD. i also installed visual web developer 2010 express and created new web mapping application using vb.net, and pointed mapresourcemanager to the one running my A.MXD. now, i want to be able to generate custom layers when user clicks on a button and add them to the service, i.e. it would display it as part of A.MXD on the website but not on the actual file on the server machine (for this reason i made it non-pooled service). plz can anyone help me with suggestions, samples, short tutorials that probably i have missed searching...? or at least a python script. thanks!


title: arcpy : Add Raster Layer layer without lyr file

text: i wonder how i can add a raster layer to a remote (not the current) mxd file with arcpy. I need to not use a lyr file for that. For vector Layer i have the following lines working fine :               mxdPath = "D:\FormationSpecialisee\IGN_IFI\doc1.mxd"     mxd = arcpy.mapping.MapDocument(mxdPath)     layer = arcpy.mapping.Layer(fcADD)     arcpy.mapping.AddLayer(df, layer , "AUTO_ARRANGE")      but it seems to not work for Raster Layer. Thanks in advance. Guillaume


title: How do I get a better thumbnail from ArcGIS Online using REST interface?

text: The thumbnail field points to a very small image that is (ostensibly) generated by ArcGIS - I would like to either make a call to generate a larger image, or specify that a larger thumbnail is required. I've tried adding size=1024,768 to the thumbnail image URL, but that has no effect. The same tiny image is returned. example:               http://www.arcgis.com/sharing/rest/content/items/         c8fee95d9cce4ffc8e382e710e555fef/info/thumbnail/         ago_downloaded.png?size=1024,768&token=oH7Yu8sYVL8yWP4ZBWZTA7-....      Will one of the other myriad services provide this functionality? Does anyone have experience with this? I can't believe I'm the only one struggling with this.


title: Calculating millions of distances with Road Graph

text: I'm new to GIS but thanks to the answers I got in another question here on `GIS stack exchange` I was able to increase my knowledge significantly. This question is related to an earlier question I asked and got answer, but since this is a new question, I thought better to ask a new question instead of continuing on the older one. I have about 160,000 addresses (points) and a (vector format) map of roads and 10 city limits (polygon). I want to:   * Calculate the distance from each of these addresses to the city borders (or possible the city centers - a point - if that is easier).   * That is, for each city I want 10 distances. I don't just want the possibility to do it one by one - I want the total 1,600,000 distances once and for all! I was earlier advised to use either `Road Graph` for `QGIS` og `pgRouting` for `PostGIS`. I have stated using `Road Graph`. I was able to import all the data into `QGIS`, install the `Road Graph` plugin, and calculate the distance between to points. My question now is:   * How do I get `Road Graph` to calculate all the 1,600,000 distances? And can it at all be done with 'Road Graph'? Thanks


title: Convert an area to meters² in PostGIS

text: I have this query :               SELECT ST_AREA(geom::geography) FROM region      The (geom) contains polygons from a shapefile imported into PostGIS with shp2pgsql. I've appended ::geography, else the result of the function was displayed in degrees. Now, I have this kind of answer :               544163926.22948      How can I know in which units it is displayed ? square/meters ? square/foot ? I want it to be displayed in square/meters (m²) ou square/kilometers (km²) How can I do that ?


title: How do I edit feature attributes, e.g. a text field with a typo

text: A feature class in a file geodatabase has a hand full of features where a text field is typo'd. How do I edit this to fix the typos?


title: Copy two datasets that particpate in the same relationship class

text: I have an ArcSDE GDB that looks something like this: ++ GDB   +++++ Dataset1   ++++++++ FeatureClass1   ++++++++ FeatureClass2   +++++ Dataset2   ++++++++ FeatureClass3   ++++++++ FeatureClass4   +++++ Table   +++++ Table_FeatureClass1_RelationshipClass   +++++ Table_FeatureClass3_RelationshipClass I wanted to create a backup of the two datasets as follows:               arcpy.Copy_management('C:/GDB.sde/Dataset1', 'C:/Extract.gdb/Dataset1')     arcpy.Copy_management('C:/GDB.sde/Dataset2', 'C:/Extract.gdb/Dataset2')      The problem is that when I run the first statement it automatically copies `Table` and all related classes, including `FeatureClass3` from `Dataset2`. Then, when it gets to the second line, the script fails because it's trying to create `Dataset2` which already exists. It doesn't give a specific error, but when I remove the relationship classes everything works fine. So does anyone know how to copy multiple datasets that share a related table like this? **UPDATE** : Following blah238's comment below, I tried using the Export XML Workspace tool with the SDE connection as the input, and a local XML file as the output. It ran fine for most of the datasets, and then hit one it didn't like an threw an error 999999: Failed to execute. I don't see anything peculiar about the dataset, other than it's empty and not registered as versioned. I then tried exporting each of the two datasets individually, which worked but was quite slow and added a few hundred MB to the total file size. Not really a practical solution. Also, I found this ArcGIS Forums post describing the same issue. There's an explanation of how to list FC's that participate in relationship classes, so I could potentially check against that before trying to copy.


title: remove segments of lines in R

text: I have these 3 lines:               line1 <- matrix(c(-1.81209, -1.80470, -1.80616, 55.68024, 55.67947,                         55.67544), ncol=2)     line2 <- matrix(c(-1.80468, -1.80330, -1.79690, 55.68322, 55.68017,                       55.67800), ncol=2)     line3 <- matrix(c(-1.79798, -1.80069, -1.78929, 55.68364, 55.68138,                       55.67966), ncol=2)          line1L <- Line(line1)     line2L <- Line(line2)     line3L <- Line(line3)     my.lines <- Lines(list(line1L, line2L, line3L), ID="my.lines")     myLines <- SpatialLines(list(my.lines))          proj4string(myLines) <- CRS("+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0")     myLinesSpatial <- spTransform(myLines, CRS("+init=epsg:27700 +datum=WGS84"))          plot(myLinesSpatial)      ![enter image description here](http://i.stack.imgur.com/DrtFt.jpg) I need to discard any part of any line, that falls within 100m of any other line. This will result in segments of each of the three lines being removed. How can I do this?


title: Easiest way to check if z values of contour lines are correct

text: What is the easiest way to check if z values of contour lines are correct in ArcGIS? I need just to look at the values, to see if equidistance between them is correct and unchanged.


title: Error while publishing a mxd as service uing python

text: I have created an Excel file with Lat/lon values. I added this data tas a layer in the ArcMap. And I "used Display XY data" otion available on right click. I selected the GCS_WCS_1984 co-ordinate system. I then saved it as ArcMap document(.mxd). I get this error when I am trying to publish my MXD file(created above) as service using python script: Layer's data source has a different projection [GCS_WGS_1984] than the data frame's projection (Code 10001) Any pointers will be helpful. I fixed it. For the benefit of others the problem was due to the wrong option which I chose suring creation of shape file. I changed the option to "the data frame" and bingo the error was gone.Cheers !


title: ArcMap - Symbolism Wildcard

text: This might be one of those slap your forehead so simple questions. In ArcMap, how can I use wildcards to aid with map symbology? For the example below, I have tried adding values such as "Field *" and "Field %" but have not had any luck. This would be convenient with tables of 4000+ records with different ending comments. Field Verify Field Verify - Unclogged Field Verify - Too Much Traffic Field Verify - Under Fence Field Verify - Buried Field Verify - Under Asphalt **Edit:** I am scratching around with this code, except I would like it to first _group_ the results and then _insert_ , instead of replacing, the group into the existing symbology. I use the field name, and then repeat the field name in the query, so that tField = FID and tQuery = FID > 500               import arcpy          tLayer = arcpy.GetParameterAsText(0)     tField = arcpy.GetParameterAsText(1)     tQuery = arcpy.GetParameterAsText(2)          mxd = arcpy.mapping.MapDocument("current")     lyr = arcpy.mapping.ListLayers(mxd, tLayer)[0]     arcpy.SelectLayerByAttribute_management(lyr, "NEW_SELECTION", tQuery)     fList = []     rows = arcpy.da.SearchCursor(lyr, [tField])     for row in rows:       fList.append(row[0])          if lyr.symbologyType == "UNIQUE_VALUES":       lyr.symbology.classValues = fList       lyr.symbology.showOtherValues = False          arcpy.RefreshActiveView()     arcpy.RefreshTOC()     del mxdde here


title: Running a new iteration for each column of a dataset in ArcGIS

text: I'm pretty new to Modelbuilding, having just been turned onto it about a week ago on here. So my problem is, I currently have a dataset of habitats that change through time (from Immature to Mature, or remaining Mature, etc). I would like to create a shapefile for each habitat type for each year. Currently I am able to do this by hand, running a Feature Selection iteration on each field and running the "feature class to feature class" tool to export it as a shapefile. What I was wondering was, is there any way to automate this process across a number of different columns? So, run the iteration and shapefile export for the 1995 field, 1996 field, 1997 field, etc, etc. I'm trying my best right now to sort it out in Python, but today was day 1 of learning and was hoping that there was another method out there that someone may be able to point me in the right direction of. Hopefully this makes sense, if you need any more info, please let me know. Thanks very much in advance. **EDIT** - Example data - The cut off for going from young (or Disturbance) to a mature forest is 40 years. So the field "Year" is the year the forest began growing (or year since it was removed) and under each year is the forest class for that year.               FID    Year    1995           1996          1997         1998     0      1992  Disturbance  Disturbance  Disturbance  Disturbance     1      1854  Black Spruce Black Spruce Black Spruce Black Spruce     2      1960  Disturbance  Disturbance  Disturbance  Disturbance     3      1875  Deciduous    Deciduous    Deciduous    Deciduous     4      1956  Disturbance  Disturbance  Jack Pine    Jack Pine     5      1957  Disturbance  Disturbance  Disturbance  Lowland Conifer      And this is the model I am currently using, I would like this to be re-run across all of my year fields (columns) ![enter image description here](http://i.stack.imgur.com/hzgGT.jpg)


title: Postgis create polygon extent/convex hull from spatial table

text: Please help, I have a polygon layer representing building footprints. I want to create new view in postgis that will contain rectangle representing bounding box and also another view that will have convex hulls polygon of all records in buildings table. I'm creating something like this:               CREATE OR REPLACE VIEW xxx.extent AS SELECT ST_ConvexHull(buildings_2013.geom) as geom from xxx.buildings_2013;      but I won't create spatial view. Could you please help?


title: How to define a Feature Set parameter for an ArcGIS python toolbox tool?

text: I'm trying to create an ArcGIS python toolbox tool which use a polygon for input. I want the user to be able to create the polygon when the tool is run. I know how to do this using both a modelbuilder model A quick tour of using Feature Set and Record Set and for a python script in a normal toolbox. This is done using a Feature Set variable and set the schema to a polygon feature class. I know from here Defining parameter data types in a Python toolbox that the datatype for Feature Set in a python toolbox is GPFeatureRecordSetLayer. If I define the parameter as below it is not possible to create a polygon at run time. I’m almost 100% sure that this is because I have not set a schema. But does anybody know how to set the schema or know if there is some other way to interactively create a polygon and use as input for a python toolbox tool?               InteractiveParameter = arcpy.Parameter(     displayName="Interactive polygon",     name="toolname",     datatype="GPFeatureRecordSetLayer",     parameterType="Required",     direction="Input",     enabled = True)


title: Is it possible to know last update of data behind a WFS service?

text: I need to use several WFS service which are not mine and convert them to shapefiles. I'd like to update the shapefiles whenever the data is updated in the WFS services without downloading the whole data. Is there a way to know the moment when data changed?


title: Is it possible to query via bounding box in the Google Maps Static API?

text: We have a `Qt` program written in `C++` to request `wms` tiles and display them. This works fine. Now we want to display imaginery from `Google Maps` and due to some restrictions in our program, we can't use javascript and we have to directly manage the url. For example, in the wms case:               www.idee.es/wms/PNOA-MR/PNOA-MR?SERVICE=WMS&VERSION=1.3.0&FORMAT=image/jpeg&CRS=epsg:32629&LAYERS=PNOA-MR&STYLES=&REQUEST=GetMap&EXCEPTIONS=INIMAGE&TRANSPARENT=TRUE&WIDTH=256&HEIGHT=256&BBOX=621921.34,4794536.07,622951.58,4801396.70      So I think the right choice for `Google Maps` would be the static API. Now, reading it I didn't find anything regarding querying by bounding box. Only querying via a center and a zoom level. Is that the only way to retrieve an image? Wouldn't it possible to query with a bounding box without specifying a center and the zoom level and retrieve an image representing that extension? Just like with the `wms` case before. **NOTE:** I would like and url example if possible. We can't use `javascript` for this.


title: How to get additional subbasin areas designated by points in the river network？

text: **Update 14/1/9** ：I follow this post：Watershed and Stream Network Delineation, and here's my workflow：   1. use **Batch Point Generation** to create batch points related to designated points(green dots)   2. **Batch Subwatershed Delineation** using these batch points and I got this： ![enter image description here](http://i.stack.imgur.com/rjKBQ.png) It looks great, but I still can't clip existing subbasins using designated points to get smaller subbasins inside. Are there tools in **Arc Hydro Tools** or any other tool like **QGIS** or **GRASS** that can help me get my goal？ * * * I have got the subbasin area of every river segment using Arc Hydro Tools： Thanks Martin! How to get the subbasin area of each specified river segment？ ![enter image description here](http://i.stack.imgur.com/JE3NM.png) Now I have some sample stations(the green dots) in the river network, and I'd like to take them as the segment bases, thus I would get more segments, and I want to split the drainage line based on those points and calculate the subbasin area of these new segments： ![enter image description here](http://i.stack.imgur.com/vXZVG.png) I think the key is how to make ArcGIS take these green dots as confluences of rivers as well. And all I want is like this schematic photo： ![enter image description here](http://i.stack.imgur.com/Y5exJ.png) What can I do？ Thank you for giving me any suggestion！


title: Combining MS Office with GIS

text: I would like to incorporate MS Office into my GIS program. I have been using QGIS and was wondering if this would be possible? I am trying to create a one stop shop for my staff to be able to write reports and print maps from one program instead of going between programs. An example would be I have a complete form with blank fields. The fields would be pulled from the map attributes. Any help would be appreciated. Thanks


title: Problem loading WMS layer from external server in Leaflet

text: I'm trying to use the WMS service documented here in a leaflet map. The documentation shows an example for openlayers:               var wms = new OpenLayers.Layer.WMS(       "Population Density",       "http://sedac.ciesin.columbia.edu/geoserver/gwc/service/wms",       {layers: 'gpw-v3:gpw-v3-population-density_2000'}     );      but I want to use this with leaflet. My own WMS hosted layers work fine with the following syntax:               var wmsRivers = L.tileLayer.wms("http://zzz.zzz.zzz.zz:8080/geoserver/opengeo/wms", {       layers: 'opengeo:rivers',       format: 'image/png',       transparent: true,       version: '1.1.0',       attribution: ""     }).addTo(map);      The following does not:               var wmsFootprint = new L.tileLayer.wms("http://sedac.ciesin.columbia.edu/geoserver/gwc/service/wms", {       layers: 'gpw-v3:gpw-v3-population-density_2000',       format: 'image/png',       version: '1.1.0',       transparent: true,       attribution: ""     }).addTo(map);      After several minutes waiting for sedac.ciesen.colombia.edu, I get a 400 `Failed to load resource: the server responded with a status of 400 (Bad Request)` Is there a way to ping the server to know it's functioning? What about version, I assume this refers to geoserver; how do I know what version they are running? Any other parameters I may be missing?


title: Layer Turned Off But Pop Up Showing

text: I'm having an issue with my popup in my map made using ArcGIS online. One of my group layers contain 4 layers. When I turn all of the layers in this group layer off, the second layer's popup still displays when I click on the map. Does anyone know why this is happening? Also, for some reason, when I have all the layers in this group layer turned on, the second layer's pop up (polygon feature) shows before the 1st layer's popup (a point layer). How do I get the first layer's pop-up to show up first before the second layers?


title: Does OpenLayers support right to left languages?

text: Right to Left languages need to be handled differently. For instance, in popups the text and the close box should be reversed. Google Maps supports this if the page has a `<html dir='rtl'>` element, how about OpenLayers?


title: Auto Generate kml files along the roads for more than 8 points

text: **What I am trying to do:** Auto Generate kml files for bus routes. **Input:** Bus stop list for each route with name and lat/long (every list has more than >30 bus stops.) **Output:** kml file of the route passing through the bus stops and along the roads. **Research so far:**   1. Found this tool from geocodezip.com to plot routes with more than 8 points.    2. Found this excellent tool to write kml files from manually clicked points. It works beautifully for less than 8 way-points.  **Possible direction:** Ideally I want to combine the functionality of both the the above tools and have a tool that writes kml code for routes without restriction on no. of waypoints.


title: Online GIS database building

text: My first post here and sorry it’s a one asking for help. A little background; I work a lot with GIS, mainly Arc but also MapInfo and Quantum, and therefore I am proficient in the use of licensed, desktop based GIS. But, I have little to no experience of web based systems. For obvious reasons I cannot go into too much detail so apologies if things are a little vague. This is in the thoughts phase at the minute so may I’m just trying to pull some ideas together. I’m looking to set something up (online) where users can log in and store a point or two on a live online GIS that other logged in users will be able to see. Basically where a user can zoom to a location and create a point, or 2/3, and then add a few bits of information for each point. I’d then like it if a user can log in and amend/delete their information when needed. Can anyone recommend some good, free (UK mainly), online GIS services that may tick a few of the boxes above? Or, any hints/tips as to what may be involved? Out of interest, does anywhere host MasterMap free?


title: How to achieve labeling in openlayers as Google did?

text: I am tried to find a solution for my question. Please look the images below. **Google** ![enter image description here](http://i.stack.imgur.com/9ialm.jpg) Google are awesome!!.. They are doing fantastic technique in labeling part. Please look the marked labels, the labels are aligned automatically according to the available space. **Openlayers** ![enter image description here](http://i.stack.imgur.com/MVdJW.jpg) This is Openlayers output. Here I just added some vector features to the map and given some styles, but the labels are overlapped together,which look congested. **My Queries**   1. Is there any way to avoid this overlapping like Google done?   2. Any Openlayers renderer (OpenLayers.Renderer) available?   3. Is there any automatic label positioning technique available in Openlayers? Thanks for any reply..


title: How to save raster results of IExtractionOp?

text: I'm using IExtractionOp to clip a raster by polygon:               IExtractionOp operation = (IExtractionOp)(new RasterExtractionOp());     IRasterBandCollection rasterBandCollection = (IRasterBandCollection)operation.Polygon((IGeoDataset)rasterDataset, polygon, true);      I tried saving via               // Fails with exception     IRasterDataset dataset = rasterBandCollection.Item(0).RasterDataset;     dataset.Copy("temp" + extension, (IWorkspace)rasterWorkspace);          // Successful, but does not show any pixels in ArcMap when opened     rasterBandCollection.SaveAs("temp" + extension, (IWorkspace)rasterWorkspace, rasterDataset.Format);          // Same result as SaveAs     ((ITemporaryDataset)dataset).MakePermanentAs("temp" + extension, (IWorkspace)rasterWorkspace, rasterDataset.Format);      None of the three methods above worked... I can use the returned result and add it to the map, so the extraction seems to be valid, but I'm unable to store it anyhow.


title: ArcGIS Error 999999: Feature Compare Tool

text: I'm running the ArcGIS feature compare tool against 2 point feature class layers. I'm comparing 9.1 SDE layer against a 10.1 SDE layer. In the feature compare tool one of the input options is 'continue comparison': I checked off the Continue_Compare option to run in order to check other properties after encountering the first mismatch. However, every time I check off this option, and run the tool is runs for >30 mins and I get the noted error. When I uncheck this option the tool runs smoothly and only takes <2 mins, without error. Just wondering why there is a difference and an error occurs.


title: How to save view so all maps are same when exporting maps?

text: I have been exporting a number of different projections for a study. However, the problem comes when I want to edit features or zoom in or out. When I then return to the Layout view, the map is not in the same position or orientation. This results in slightly different maps and when I put them together the difference in zoom and position of labels is quite obvious. I can of course attempt to match them by moving the map around and zooming in Layout view, but you don't get it exactly right (and the slight difference is so clear when flipping through maps in a presentation). Is it possible to create a template, or something like it, in Layout view so that no matter what zooming or moving around you do in Data view, you can return to the exact same position to export an identical map?