Export a join table to arcmap

I have shape file (Point) and a table which I join using the following code.               public ITable JoinLayer_Table(IFeatureLayer pFeatLayer, IStandaloneTable pStTable)             {                      IDisplayTable pDispTable;                 IFeatureClass pFCLayer;                 ITable pTLayer;                 pDispTable = (IDisplayTable)pFeatLayer;                 pFCLayer = (IFeatureClass)pDispTable.DisplayTable;                 pTLayer = (ITable)pFCLayer;                      IDisplayTable pDispTable2;                 ITable pTTable;                 pDispTable2 = (IDisplayTable)pStTable;                 pTTable = pDispTable2.DisplayTable;                       IRelationshipClass pRelClass;                 IMemoryRelationshipClassFactory pMemRelFact = new MemoryRelationshipClassFactory();                 pRelClass = pMemRelFact.Open("TabletoLayer", (IObjectClass)pTLayer, "HouseNo", (IObjectClass)pTTable, "UNIT_ID", "forward", "backward", esriRelCardinality.esriRelCardinalityOneToMany);                           IRelQueryTableFactory pRelQueryTableFact = new RelQueryTableFactory();                 ITable pRelQueryTab = (ITable)pRelQueryTableFact.Open(pRelClass, true, null, null, "", true, true);                 return pRelQueryTab;             }      The join results is an ITable. Now I need to make a shape file from the Join results and add it to the TOC in ArcMap. How can I do this? Any code example would be greatly appreciated. PS: Actually I wanted to join 2 shape files ( a point shape file and a polygon shape file) and get a point shape file for the JOIN results, but could not do that . So thats why I have a point shape file and the corresponding dbf of the polygo shape file to create the Join. Any help on this also would be greatly appreciated.


Modelling Sea level rise scenarios for a Coastal county in ArcGIS

I am a graduate student @ Clemson University and am working on a GIS based project which is to model different sea level rise scenarios and estimate its impacts on a coastal county and to suggest measures to reduce the risk involved. I have requested data from FREAC (Florida Resources and Environmental Analysis Center) and got 5ft. and 10ft. DEMs for 2 counties. I also have shapefiles of the counties with the parcel data. Now, I am trying to figure out the best way to model the scenarios. Let me know if there are any methods used previously or any suggested methods. Thank you


Reading National Elevation Dataset (ArcGrid/GridFloat/IMG) with python only tools

I've found high precision elevation (1/3 - 1/9 arcsecond resolution) data from the National Elevation Dataset provided by the USGS. It comes int IMG, Arcgrid, and GridFloat formats. I know ArcGrid corresponds to a paid software package, but I'm trying to stick to using freely available tools. I have GPS data that I'm trying to correlate with ground level. Are there any python libraries which will let me transform gps data into altitude using a datafile in IMG, ArcGrid, or GridFloat formats? If you know of an example that would be extremely helpful. Thanks!


GeoWebCache and OpenLayers horizontal resolution error

I have a mapping application (the Arne G. Brekke Bygdebok Collection) which uses OpenLayers to load map tiles from a GeoWebCache server operated by the Statens Kartverk (Norwegian Mapping Authority). It's worked nicely for two years, but this morning the tiles all load as pink, with errors like "400: Requested horizontal resolution: 0.0439453125 , best match: 0.03333333333333333 exceeds 10% threshold. Perhaps the client is configured with an incorrect set of scales (resolutions), or the DPI setting is off compared to the one in GWC?" if you view the image directly. Here is a test case I worked out demonstrating the problem:               <!DOCTYPE html>                   <meta charset="UTF-8" />         <title>Arne G. Brekke Bygdebok Collection</title>
Arne G. Brekke Bygdebok Collection

                                    That's the simplest test case I can work out for it. Based on that, OpenLayers generates URLs like this one (broken onto multiple lines for clarity):               http://opencache.statkart.no/gatekeeper/gk/gk.open         ?LAYERS=topo2         &FORMAT=image%2Fjpeg         &SERVICE=WMS         &VERSION=1.1.1         &REQUEST=GetMap         &STYLES=         &SRS=EPSG%3A4326         &BBOX=11.25,56.25,22.5,67.5         &WIDTH=256         &HEIGHT=256      How can I get it to work again? Bear in mind that I have zero control over the map server -- I can only adjust the JavaScript on my end.


Working with legend patches/symbology in ArcGIS

I am working with a legend in a map from ArcGIS 10.1. The map is used as a template and will be used many times to create different maps. The legend is used to show a soil type and its descriptor: ![enter image description here](http://i.stack.imgur.com/br6z0.jpg) So, in other words, I have set my layer symbology to unique values, add all values, single field and I have turned off all the symbology so that the patches do not show up in the legend. The solution I am looking for, is to avoid having to manually turn off the symbology all the time. I just want to 'Add All Values' and have the item and Description populate in the legend without the symbology patch. I don't think there is a checkbox/button that does this in the layers properties symbology tab, but I was trying to figure out how to do this in the legend properties dialogue box. So far, I can't seem to figure it out. It seems like any of the styles under the Legend-->Properties-->Items tab requires a patch, label and description. There is no option to just choose the label and description. Am I missing it somewhere? Thanks, Mike


How can I delete dangling paths after merging features

In my previous question I've asked how to can I merge different polygons while converting to GeoJSON and TopoJSON. After a day of researching I've realized that, rather trying to do it command line, I should use some editor which can handle Shapefile data. So, I've choose QGIS. I've toggled to edit mode, chose all areas needed and called "merge selected features". After the merge I'm getting something like this: ![enter image description here](http://i.stack.imgur.com/7PZ8P.png) I'm expecting to get all countries in Africa joined into one region. But, as you can see, I'm getting some dangling paths. My question is - what I'm doing wrong, how I can avoid such dangling paths or get rid of them after the merge? Since selecting all this nodes via "Node tool" does not help much.


Web-based thumbnail image in 'Item Description' for a Python toolbox [ArcGIS 10.1]

I'm using locally stored images as thumbnails in the Item Description for a python toolbox [ArcGIS 10.1]. I'd like to load the images from an online source so other users of the toolbox can still see the images. Is there a way to do this?


How to create arrows lines automatically on ArcGIS, for a thematic map of directed flux?

![Sample Flux Map](http://i.stack.imgur.com/hTXsg.jpg) I need to create a Map with the Import and Export Flux between countries. I have the needed data in my table. But I wish to know if there is any tool or extension which allow me to create those lines or arrows automatically (similar to how we do with labels, but those ones create word, which organize themselves in the `Data Frame`). For that all that I should do is to indicate lineament of origin and destiny, representing the seaports. I need to do something like this, but with arrows:


Network Analyst - OD Cost Matrix snapping distance

I am using the ArcGIS Network Analysis - OD Cost Matrix layer tool and I was wondering if anyone knew of a way to incorporate the snapping distance of the origin and destination to my total distance value. Right now, the Origin snaps to the closest line feature in my network and the same for the destination point, and the tool only gives me a distance between these 2 snapping points, but I would like to include the 2 snap distances as well. Does anyone know if this is possible? I know I could use the Near command ontop of this tool to get the total distance, but I was hoping there was a setting in Network Analyst that could solve my problem.


Is d.mon xN deprecated in GRASS in favor of something in the wxGUI?

I'm trying to learn command-line GRASS (6.4.3) using the KyngChaos binaries on a Mac (10.8.4) from Neteler and Mitasova's book, Open Source GIS: A GRASS GIS Approach (3rd ed.). I've learned from this site that some of the commands in the book have been deprecated with the introduction of the wxGUI. Is d.mon one of those commands that has been deprecated? I've found that I can still invoke a X11 graphics monitor using d.mon, but is it preferable to use the graphics monitor in the wxGUI? If so, then can I access the wxGUI's graphics monitor from a terminal session? I'm finding that the command console in the wxGUI behaves somewhat unpredictably compared to a terminal (e.g. I can't cd to change working directory). Also, am I likely to continue running into inconsistencies between the book and the current version of the software as I work my way through the book?


Populate print template title text from a dojo textbox

I have a print widget that starts up great. It's based on the typical ArcGIS print widget sample. However, I want to tweak this a bit. There is a "titleText" property in the layout options that I want to set from a text box that the user types a title into. So, based on the dojo textbox value, I assign that to the titleText property. A lot of the code I find is in AMD so it's really tough to decipher (as I am still starting out in JSAPI). Thanks in advance for any help!


How to add Excel-sheets with meter coordinates to QGIS?

I have an Excel-sheet where the coordinates are in meters. What is the most efficient way of inputting this file to QGIS? There are thousands of files to put so please recommend some straight-forward solution. Thank you, Kimmo


Styling broad lines with short dashes

I'm looking to style some lines in QGIS to look like the steps in OpenStreetMap's default style. Actually, I'm styling those same steps on my own machine. ![steps from OSM](http://i.stack.imgur.com/w1qd4.png) Basically, unless I'm missing something, I need to make a dashed line that is wide with very short dashes. I've tried setting a custom dash pattern, but qgis doesn't seem to take values of less than 1 because the size of the dash doesn't shrink after that point. So I can get a custom dash pattern with a dash that is as short as it is wide, but no shorter. Is there a way past this, or some other way to style the lines to look like steps?


Problem with access to ArcGIS Server Manager

ArcGIS For Server 10.0 _Firefox can not establish a connection to the server._ (Error in all browsers) ArcGIS Server Manager problem In ArcCatalog i can connect but can't _start_ _Configuration MapServer can not be started._ Also in computer management missing _Local Users and Groups_ What is my fault? Thanks


Why do the components loaded by using WMS and WFS from GeoServer not the same?

I have tried to load shape files into GeoServer, but by loading into QGIS, the WMS and WFS versions are not the same. The WMS version seems like correct. They are both loaded by using WGS84. Do you have any suggestions? ![The green rectangles are from the WFS version, while the blue border ones are from WMS version.](http://i.stack.imgur.com/RrJBw.png)


Splitting lines and assigning two new values

It is my understanding that when a line is manually split using the split tool in ArcGIS 10, the geoprocessing follows a create and copy format, meaning that one new feature is created, while the other retains all of the values of the original line segment. My question is, is it possible to set up the behavior so that when a split is performed, a delete and two creates is executed. A little background to my question: I am working with a sewer dataset and utilitizing the attribute assistant add- in provided by ESRI. I have a rule set up to generate a new facility id upon creation of a new main. This works fine for one half of the split line, but since the other half is not considered a "new" line, it retains the facility id of the original segment. I'm hoping there's a way to modify the split tool to create two new, unique lines upon split. Any guidance would be most appreciated!!


How to document field properties of attribute in shapefile?

I would like to document within ArcGIS how the values for specific shapefile attribute fields were derived or obtained. Is it possible to input and export this information as a report with ArcGIS 10?


Old Windows Installers for python library fiona

Due the the recent bug as presented here for fiona 1.1, i am obliged to use a version of fiona in 1.0.x on windows 8.1. however, i can't find older versions of the installer at here. Are there any sources for the archive installers, specifically for the 32bit fiona (1.0.x) ? thanks in advance! P.S. I think it will apply to ALL other previous libraries, and not just fiona


QGIS Bookmarks not transfering with files

My students send me QGIS map files that when I open them have been stripped of the bookmarks. Is there a setting I'm missing or is the normal case.


QGIS map file to be viewed on both mac and pc platforms

Im teaching a class where some students have QGIS installed on their macs and I'm using a pc. Although the file structure is the same mac reference files differently than the pc and I get data link errors when I open a students map project. Any quick conversion methods? I can do it a link at a time but that seems a bit slow to me.


What format to store/use four band Raster data

I have four band aerial photography - RGBI. I need to share it with my users, a mixture of MapInfo, QGIS, and ArcGIS. What format should I store/serve the data in? GeoTIFF isn't well suited to four bands because some applications assume the fourth band is transparency. ECW and JPEG2000 are proprietary, plus not all applications support four bands. I can split it up into different layers/files but would prefer not to.


How to check a polyline is self intersecting

How to check weather a **single polyline** is self intersecting to itself. using arcobject C#


get layer names from a group within a group

Trying to get the layer names from a group that are inside a group. I have the code the gets the layer from a group, but not group inside group. any ideas?               Dim pDoc As IMxDocument     Set pDoc = ThisDocument     Dim pMap As IMap     Set pMap = pDoc.FocusMap          Dim pCompLayer As ICompositeLayer     Dim pCompLayer_child As ILayer     Dim pGrplyr As IGroupLayer     Dim i As Integer     Dim l As Integer     Dim k As Integer                    For i = 0 To pMap.LayerCount - 1      If TypeOf pMap.Layer(i) Is IGroupLayer Then      Set pCompLayer = pMap.Layer(i)          For l = 0 To pCompLayer.Count - 1      UserForm1.ComboBox1.AddItem pCompLayer.Layer(l).Name           ' detect if group is within group?? stuck here      Set pCompLayer_child = pCompLayer.Layer(l)      If TypeOf pCompLayer_child Is IGroupLayer Then      Set pCompLayer_child = pCompLayer.Layer(l)            UserForm1.ComboBox1.AddItem pCompLayer_child.Name           End If          Next l               End If     Next i          UserForm1.Show          End Sub


Problem with using a circle as polygon for querying PostGIS

I'm trying to query PostGIS and find some object features that intersect with a circular buffer created by the user using OpenLayers, JavaScript and PHP. I use the following code for creating the circle:               var circle = OpenLayers.Geometry.Polygon.createRegularPolygon(         new OpenLayers.Geometry.Point(X, Y),         1,         30     );      and then I'm planning to use the following SELECT Query in my php code:               SELECT b.geom     FROM buildings AS b     WHERE ST_WITHIN (b.geom, ST_MAKEPOLYGON(circle, srid));      Q1: how can I pass the geometry created for the circle from my .JS file to .php file? can I pass it as an object? can I do something like below??:               http://mydomain.com/page.php?buffer=circle      Q2: Is my SELECT query correct? Am I using the correct function ST_MAKEPOLYGON? Any suggestions/help would be highly appreciated.


Tilemill Index problem

I have created my own customized tile. Now I am facing some problem with its indexing or naming convention. I am giving an clear example. I have taken an area where 6 tiles were created. In Google Maps or OSM in the same extent its giving 6 tiles and the same tiles. So there is no referencing or projection problem. Even the X or the folder name is also matching. I am getting trouble with the Y column or Image name. Its not matching. I have taken a close look at it. Interestingly I found in Google or OSM the Y column is incrementing from top to bottom where as tilemill created tiles number those incrementing from bottom to top. That is why the naming of tiles are not matching. What should I do?


Where to get Land Use data for the UK?

We're currently looking for a Land Use dataset for part of England (UK) and were wondering what datasets are out there, either free or commercial. One's we've found so far:   * http://www.infoterra.co.uk/landbase   * http://www.geoinformationgroup.co.uk/products/national-lu   * http://data.gov.uk/dataset/land_use_statistics_generalised_land_use_database (more generic statistics) **Important** : We're not looking to generate this data ourselves, We're only interested in nicely pre-packaged datasets.


Can't find the New label tools in QGIS

I can't find the New label tools in QGIS developed by sourcepole. I'm using QGIS 1.7.4. Is it a QGIS tool or a pluggin? Thanks!


Display google maps with Geotools?

I would like to display Google maps (Raster and Vector) in Java using GeoTools. I am stuck because I am new to this library. Are there any tutorials out there that can help me ?


How to convert a CSV from DMS to WGS84 in QGIS?

I would like to automatically convert from a Degree - minute - second system to WGS84 EPSG:32632. I have a CSV with these kind of values:               LATITUDINE  LONGITUDINE     46°01’24,7” 11°13’47,5”     45°42’07,5” 10°55’11,3”     46°01’37,6” 11°06’41,7”     46°15’03,7” 11°11’00,1”      and I would like to automatically convert them. There is a way? They are a lot, if not I would do with a normal converter... I try to import them in QGIS (but only a few of them works) and export them in some way, but I sincerely do not understand how to. Any suggestion? Thanks!


Add hour and seconds to CartoDB/Torque animation

A while ago I made a CartoDB/Torque animation that displayed hours and minutes while playing. Sujay v. helped me out with the code here: https://groups.google.com/forum/#!topic/cartodb/xpxF4fSV7ZI Essentially I changed some lines in the torque.js file as below:               torque.clock.set('' + date_arry[0] + ' ' + date_arry[1] + ',  ' + date_arry[2] + ' ' + date_arry[3] + '');      However I've come back to needing to do something similar, and the torque.js file has changed alot and I can't figure out how to implement something similar now. Would someone be able to advise? The source code is here: https://github.com/CartoDB/torque/blob/master/dist/torque.uncompressed.js Apologies I can't do this myself. I'm pretty good at SQL and GIS, but lacking the .js skills. :-( Thanks James


issues migrating + updating an older vb.net add-in

Key points to consider: -The original add-in was written in VS 2005 for ArcGIS 9.1, uses Framework 3.5, with now deprecated namespaces, but those don't appear to be problematic -Migration is to 9.3.1 -Got my hands on the source + already compiled DLL and TLB for the project -The DLL runs provided I register with regasm, however a conversion issue (string-to-date) occurs because the code takes for granted that system format is in mm-dd-yyyy format and not dd-mm-yyyy format. The MonthCalendar form object is quite rigid in its format, so a slight alteration of the code is needed to account for this exception. That's the easy part. -I opted to use VS Express 2013 since it's the only one I could find immediately. Opening the source project triggers a single-direction conversion of the project file, but I don't believe any data was lost. I am informed that some namespaces are deprecated. I try to fix up the reference and eventually manage to build. "Register for COM Interop" is checked. I also used the 3.0 version of Microsoft PowerPacks rather than PowerPacks.Vs which was originally referenced. I cannot point to this now that the new powerpacks is installed anyway. -I add the newly compiled tlb to a custom toolbar, run the tool, and am informed that the Microsoft Scripting Runtime COM cannot be found, culture neutral, publickeytoken null.  What now? Does it matter that I'm using VS 2013?


Elevation and point data

On arcmap 10.1, I have elevation data that is shown as contour lines (polylines) i also have point data showing the distribution on certain animals. I want to run some kind of statistical analysis to see at what height the animals are most commonly found at, how would i do that/what tools would i use? Thanks


removing message after wfs-t transaction

Hi is there a way to remove the alert (message) after adding new object to database with WFS-T. I am using code from demo example from Geoserver. The code is shown below                 var drawPoint = new OpenLayers.Control.DrawFeature(                     roads, OpenLayers.Handler.Point,                     {displayClass: 'olControlDrawFeaturePoint'}                 );                 drawPoint.featureAdded = function(feature) {                     feature.layer.eraseFeatures([feature]);                     // cast to multipoint                     feature.geometry = new OpenLayers.Geometry.MultiPoint(                         feature.geometry                     );                     feature.style.strokeColor = "#0000ff";                     feature.state = OpenLayers.State.INSERT;                          feature.layer.drawFeature(feature);                 }      Is there a way to remove this alert. Or maybe just to get the part SUCCESS or FAILD??


Installation of SDA4PP

I have downloaded the plug in SDA4PP for QGIS. I'm trying to install it but I don't see the Add 3rd party repositories under the Repositories tab of the QGIS python plugin installer. This step is mandatory according to the installation guide. I don't know it this is due that the only download link I found was for Windows XP and I´m working with Windows 7. Does anyone have an idea how to carryout the installation? Any hint would be really appreciated. Regards.


Loading a saved pqt (saved composer file) with python in a QT Plugin

I'm working on a plugin for QGIS. The aim is to make the interface easy to use for my customer that never used GIS before. So I wanted to know if there is a easy way of opening the composer and then loading the saved pqt with python. I want to incorporate a printing option in my plugin. But the traditional composer will not work as the person will need to open composer every time an then load the saved file. Sorry if this is a stupid question, but I have no idea where to start with this one. Thanks


How to cut a shapefile using a batch?

I would like to cut a shapefile using a batch (like in the photo): ![enter image description here](http://i.stack.imgur.com/cM1RV.jpg) But I have some problems because I don't know what do I have to put in field "extent" and after putting somthing like "example.shp" to see what happens the shape cut not include the road intersected with the shape. In this case I used the SAGA's batch included in QGIS. So, if anyone knows how to cut a shape in batch or other way I will be so thankful!


pgrouting 2: missing pgr_findnearestnodedwithin function

With the PgRouting 1.5 I was using the find_nearest_node_within_distance function and it works, but now I'm moving the project to PgRouting 2 + Postgis 2 + PostgreSQL 9.1, and now that function, I think, was renamed to pgr_findnearestnodedwithin but when I try to execute the query I get the following error               ERROR:  function pgr_findnearestnodedwithin(unknown, numeric, unknown) does not exist      I was looking if there is a new way to import this function but I didn't have success. please help. there are something I'm missing? Thanks in advance


Russian in fieldname changing with Python in ArcGIS

Is there any way to rename fieldname in the attribute table of a layer/shapefile in ArcGIS 10.2 containing russian letters? If yes, is it possible to do it in a script?


How to preserve ordering of results with the wrapper functions

When running               SELECT * FROM dijkstra_sp('ways',232, 421);      I get the result ordered by gid, rather than by the order for reaching the target by starting at the source. What is the reason for that and is there a way to get the results in the _natural_ routing order? Im using pgrouting _1.05-1~lucid1_ from http://ppa.launchpad.net/georepublic/pgrouting


how to plot gps coordinates on my own raster map?

I have lat log values available in degree format, i donot know how to access raster maps using c#, and how to plot accurately. which libraries should i use and what are the controls avialable ,i am using visual studio 2010 ultimate edition. There are controls available for vector maps in which one can zoom ,pan and plot the given coordinates accurately etc etc, but i have a georeferenced raster map with 4000*5000 resolution, i have done tiling with the map,each tile is of 256*256 resolution i need some source code or a manual to use the raster map to develop a navigation app and to plot the location accuraetly on the map, thanks


problem with control ACTIVEX in conflict with Arctoolbox

> **Possible Duplicate:**   >  Problem with ArcToolbox tools I have this great and unsolvable (till today!) problem: in arcgis 9.1 (in both OS, XP SP3 and 7 ultimate 64bit) when I open Arctoolbox appears a message that it is impossible to open eevry tools (but the Arcgis works for every other performance), because there is a conflict with web browser (?) and control ActiveX. I think there is a problem in/with Windows register, but I am not be able to make changes. I have ask help to ESRI Italy and ESRI US, but with no solutions. Anyone have had the same problem? Anyone can help me? I Thanks and regards. Piero Favino Architect Milan, Italy Polytechnic of Milan


Problem with control ACTIVEX and Arctoolbox (9.3.1) on windows 7 64 bit

> **Possible Duplicate:**   >  Problem with ArcToolbox tools "one or more activex controls could not be displayed when ever you try to open any tools from Arctoolbox it display this errors. then go to start run -->regedit--> HKEY_CURRENT_USER\Software\Microsoft\Windows\Curre ntVersion\Internet Settings\Zones\ and delete the pseudonumeric key(L) or Rectangle shape before zone 0. restart the arcgis , and it works" did anyone have ever tested this command, to resolve the problem in object? What happens to the others programms and to OS? Is it secure operation? Thanks to all for the help, because I have no find a solution till today (also from ESRI Italy and US).


ArcObjects: How to make non-feature-linked annotation features with a simple line callout that remains functional?

I would like to make simple line callout annotation that is not feature-linked but remains functional when it is manually moved, i.e. the anchor point remains on the feature while the text moves. I have tried making graphic text elements with simple line callouts and loading them into an annotation feature; however, the text box and leader line are combined into a static graphic when loaded into annotation. The leader line is no longer functional; the anchor point doesn't stay anchored. How might I do this programmatically in ArcObjects? Will this only work if I automatically generate annotation from the feature class I wish to annotate? Any answers or insights would be greatly appreciated. Thanks, all.


How do I add a hyphen to the Field Calculator to merge 2 columns?

I am trying to build a field calculator expression to merge two columns into one. I am looking for the word that makes it a range. Example: I have 2 columns in elevation, `From` and `To`. I want to make them like "`4100 - 4200`" where I add the hyphen. How do I do that ? ![table showing data being worked with](http://i.stack.imgur.com/6JYRu.png)


QGIS Layer Properties Field Tab

When I setup my layer properties I selected some of the attributes to be hidden on the Fields tab. I would like to change the EDIT WIDGET from hidden to a different function however some of them are no longer available in the layer properties dialog box. How do I make them available in the layer properties dialog box under the Fields tab again. Using version 1.7


What's the most promising application (use) of OpenStreetMap data?

In some areas OSM data is very complete and has a level of detail not seen in closed data sources. Another distinguishing factor is that it contains the knowledge of the people who live in that areas. For example, commonly, colleges are mapped with many details, because students spend special care mapping the place they study. Taking into account these characteristics and others and the OSM licence, What's the most promising application of OpenStreetMap data? -Car GPS may have OSM data?   -Government could make use of OSM data?   -Companies could make use of OSM data?   -Education?   -Other possibilities?   -How?


Replacing attribute table entries using python script. Script needed

I am using ArcMap 10 from ArcGIS 10.2.1 package. I got an attribute table, with few collumns. I need to replace some column "String_A" entries with newer "String_B" entries, in the same row. In both "String_A" and "String_B" columns are some values. Some entries in column "String_A" dont need to be changed, only those, where "String_B" entries not . Both columns are strings. I am doing it by hand for now, but it would be nice to have a script in python, that can do it for me, because from time to time I need to update table with some new entries. I am learning python, but with my knowledge for today I can't do this task. So, I need some help with Field calculator and python :) From "String_B" to "String_A" replace entries, example in the picture: http://tinypic.com/r/2dcd79u/8


What traits identify 2D, 2.5D, 3D geometries, and how can you convert downward, such as from 2.5D to 2D?

In answering this question, posed by @Burton449, I came to realize I don't have a definitive grasp on what distinguishes 2D, 2.5D, and 3D geometries. I certainly realize Z/elevation values are integrated into the data, but I'm unclear if the 2.5D or 3D nature of the data is designated at both the file level (projection, headers, etc?) and the individual feature/geometry level, or just the latter. * * * [ **Update 9.26.2013** ] It occurred to me that in order to express vertical units, for the GIS to be spatially-aware, the projection might/should specifically mention the vertical datum. However if I use OGRINFO to interrogate a tiger line roads layer (tl_2012_45079_roads.shp), which came down the pipe as a 3D-enabled layer, the layer's SRS is just good ol' NAD83 (EPSG:4269). While ogrinfo can tell the layer is 3D, I'm not seeing any evidence of this in the WKT geometry. Here's a snippet of output showing the summary details and a bit of the first feature geometry:               C:\>ogrinfo C:\xGIS\Tiger\tl_2012_45079_roads.shp -geom=yes -sql " SELECT * from      tl_2012_45079_roads WHERE FULLNAME = 'US Hwy 1' " | more     INFO: Open of `C:\xGIS\Tiger\tl_2012_45079_roads.shp'           using driver `ESRI Shapefile' successful.          Layer name: tl_2012_45079_roads     Geometry: 3D Line String     Feature Count: 3     Extent: (-81.341364, 33.753578) - (-80.603160, 34.268694)     Layer SRS WKT:     GEOGCS["GCS_North_American_1983",         DATUM["North_American_Datum_1983",             SPHEROID["GRS_1980",6378137,298.257222101]],         PRIMEM["Greenwich",0],         UNIT["Degree",0.017453292519943295]]     LINEARID: String (22.0)     FULLNAME: String (100.0)     RTTYP: String (1.0)     MTFCC: String (5.0)     OGRFeature(tl_2012_45079_roads):6176       LINEARID (String) = 1102215607452       FULLNAME (String) = US Hwy 1       RTTYP (String) = U       MTFCC (String) = S1200       LINESTRING (-81.050945 33.995744,-81.050455 33.995919,-81.049132 33.996364,-81     .049051 33.996392,-81.048981 33.996415,-81.048852 33.996458,-81.04883 33.996466,     -- More  --      * * * Specifically I'd like to know:   * How are the Z-values/flags integrated in 2.5D and 3D geometries, and what is responsible for "announcing" this when a file is accessed? In other words, how might I know _at a glance_?   * Are Z-values of 2.5D and 3D geometries expressed in WKT? I realize OGR can designate the 25D option, but I'm unclear on where the z-value is expressed. Does WKT for a 2.5D or 3D geometry include the z-value?   * Finally, How might you convert down, from a 2.5D to a 2D geometry? I hope this isn't asking too much from one question. Thanks for any help!


PHP MGRS Conversion Formula/Library

Currently, the formula I'm using to convert MGRS to Lat/Long and vice versa is giving me slightly different results than most other conversion utilities. It's usually between 10 to 30 feet off. This isn't a very big difference, but it matters for the particular application of data. I'm wondering if anyone knows of a simple library that does conversions between MGRS and Lat/Long or UTM.


Bringing information from one shapefile into another using ModelBuilder?

I’m trying to use the ModelBuilder to fasten up my work in ArcGis but I don’t know how far I can use it and what’s possible. Here is my problem: I want to analyse agriculture parcels regarding to the watersheds. I have two shapefiles: one with all my parcels and one with all my watersheds. ![Origin Shapefiles](http://i.stack.imgur.com/4lkSf.png) I now used ModelBuilder to Clip the parcels with each watershed. ![Modelbuilder to clip each watershed](http://i.stack.imgur.com/ysMu7.png) This works perfectly and now I have many shapefiles with the parcels from every watershed. (for example: parcels_watershed01.shp, parcels_watershed02.shp,...) What I would like to do, is use the information I can get from these parcels_watershed shapefiles and bring them into my origin watershed shapefile. For example: I would like to know the area of agriculture parcels inside of each watershed and bring this information into the watershed shapefile. How can I do this in an automatic way? Is it possible to use ModelBuilder in order to extract an information from one shapfile into another? (I now opened each parcel_watershed shapefile, looked up the statistics of the area and wrote them down and had to fill them manually into the Attribute Table of the watershed shapefile... this is long and I’m sure there must be a better way, but i don’t know how.) I'm using ArcGIS 10.1 for Desktop.


Python script not working after exporting it from ModelBuilder

trying to figure out why my script wont work (I should mention I'm still at an intro level). This is a script I exported from ArcGIS model builder which worked fine using that:               # Import arcpy module     import arcpy               # Local variables:     DN700_Less_Than_10m_distance_M2S = "DN700_Less_Than_10m_distance_M2S"     Easement_BND = "C:\\20140211_scripting\\Sample_data.gdb\\Easement_BND"     Output_Feature_Class__3_ = "DN700_Less_Than_10m_distance"     DN700_Less_Than_10m_distance_M2S__2_ = "DN700_Less_Than_10m_distance_M2S"          # Process: Densify     arcpy.Densify_edit(DN700_Less_Than_10m_distance_M2S, "DISTANCE", "1 Meters", "0.1 Meters", "10")          # Process: Feature Vertices To Points     arcpy.FeatureVerticesToPoints_management(DN700_Less_Than_10m_distance_M2S__2_, Output_Feature_Class__3_, "ALL")          # Process: Near     arcpy.Near_analysis(Output_Feature_Class__3_, "Easement_BND", "", "NO_LOCATION", "NO_ANGLE")      Any help would be great!


How to convert a NetCDF file to GeoTIFF using GDAL

I have a NetCDF file which contains data for the number of rain days in the UK. I'd like to convert it to a GeoTIFF, but can't get GDAL correctly georeference it. See dumps from gdalinfo below. I want a single band TIFF showing the "raindays_tmean_abs" variable. The "lat" and "lon" variables seem to have sensible min/max values to put the data over the UK, but when I convert it with gdal_translate it ends up off the west coast of Africa. There is a warning about the first dimension not being a "Time or Vertical" dimension. I read that gdal expects the data to be in a certain order, so maybe this is the problem, but I don't know how to modify it. How do I convert/reproject this data? **gdalinfo NETCDF:RainDays1_1961-1990_LTA_25km.nc**               Warning 1: dimension #0 (meaning_period) is not a Time or Vertical dimension.     Driver: netCDF/Network Common Data Format     Files: RainDays1_1961-1990_LTA_25km.nc            RainDays1_1961-1990_LTA_25km.nc.aux.xml     Size is 39, 52     Coordinate System is `'     Origin = (-17.650008352179277,11.329999312758446)     Pixel Size = (0.219999614514803,-0.219999998807907)     Metadata:       axis_0#units=       meaning_period#comment=Meaning periods are all months (12) beginning January, all seasons (4) beginning Winter (DJF) and Annual (1)       meaning_period#long_name=meaning_period       meaning_period#units=       NC_GLOBAL#Conventions=CF-1.0       NETCDF_DIM_EXTRA={meaning_period}       NETCDF_DIM_meaning_period_DEF={17,6}       NETCDF_DIM_meaning_period_VALUES={0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16}       raindays_tmean_abs#comment=Converted from UKCP09 monthly data to climatologies and then averaged to 25km rotated pole grid.       raindays_tmean_abs#coordinates=lon lat       raindays_tmean_abs#long_name=Number of days per meaning period having a rainfall >= 1 mm (days)       raindays_tmean_abs#missing_value=1e+020       raindays_tmean_abs#name=rain_days       raindays_tmean_abs#title=rain_days       raindays_tmean_abs#units=days       rlat#axis=Y       rlat#bounds=bounds_rlat       rlat#long_name=latitude in rotated pole grid       rlat#standard_name=grid_latitude       rlat#units=degrees       rlon#axis=X       rlon#bounds=bounds_rlon       rlon#long_name=longitude in rotated pole grid       rlon#modulo=360       rlon#standard_name=grid_longitude       rlon#topology=circular       rlon#units=degrees     Subdatasets:       SUBDATASET_1_NAME=NETCDF:"RainDays1_1961-1990_LTA_25km.nc":bounds_rlat       SUBDATASET_1_DESC=[52x2] bounds_rlat (64-bit floating-point)       SUBDATASET_2_NAME=NETCDF:"RainDays1_1961-1990_LTA_25km.nc":bounds_rlon       SUBDATASET_2_DESC=[39x2] bounds_rlon (64-bit floating-point)       SUBDATASET_3_NAME=NETCDF:"RainDays1_1961-1990_LTA_25km.nc":lat       SUBDATASET_3_DESC=[52x39] latitude (32-bit floating-point)       SUBDATASET_4_NAME=NETCDF:"RainDays1_1961-1990_LTA_25km.nc":lon       SUBDATASET_4_DESC=[52x39] longitude (32-bit floating-point)       SUBDATASET_5_NAME=NETCDF:"RainDays1_1961-1990_LTA_25km.nc":raindays_tmean_abs       SUBDATASET_5_DESC=[17x52x39] raindays_tmean_abs (32-bit floating-point)     Geolocation:       LINE_OFFSET=0       LINE_STEP=1       PIXEL_OFFSET=0       PIXEL_STEP=1       SRS=GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9108"]],AUTHORITY["EPSG","4326"]]       X_BAND=1       X_DATASET=NETCDF:"RainDays1_1961-1990_LTA_25km.nc":lon       Y_BAND=1       Y_DATASET=NETCDF:"RainDays1_1961-1990_LTA_25km.nc":lat     Corner Coordinates:     Upper Left  ( -17.6500084,  11.3299993)      Lower Left  ( -17.6500084,  -0.1100006)      Upper Right (  -9.0700234,  11.3299993)      Lower Right (  -9.0700234,  -0.1100006)      Center      ( -13.3600159,   5.6099993)      Band 1 Block=39x1 Type=Float32, ColorInterp=Undefined       Min=10.658 Max=22.949        Minimum=10.658, Maximum=22.949, Mean=15.625, StdDev=3.151       NoData Value=1.0000000200408773e+020       Metadata:         comment=Converted from UKCP09 monthly data to climatologies and then averaged to 25km rotated pole grid.         coordinates=lon lat         long_name=Number of days per meaning period having a rainfall >= 1 mm (days)         missing_value=1e+020         name=rain_days         NETCDF_DIM_meaning_period=0         NETCDF_VARNAME=raindays_tmean_abs         STATISTICS_MAXIMUM=22.949373245239         STATISTICS_MEAN=15.624619410087         STATISTICS_MINIMUM=10.657625198364         STATISTICS_STDDEV=3.1508820627895         title=rain_days         units=days     Band 2 Block=39x1 Type=Float32, ColorInterp=Undefined       Min=7.893 Max=17.172        Minimum=7.893, Maximum=17.172, Mean=11.591, StdDev=2.282       NoData Value=1.0000000200408773e+020       Metadata:         comment=Converted from UKCP09 monthly data to climatologies and then averaged to 25km rotated pole grid.         coordinates=lon lat         long_name=Number of days per meaning period having a rainfall >= 1 mm (days)         missing_value=1e+020         name=rain_days         NETCDF_DIM_meaning_period=1         NETCDF_VARNAME=raindays_tmean_abs         STATISTICS_MAXIMUM=17.172254562378         STATISTICS_MEAN=11.591275386767         STATISTICS_MINIMUM=7.8926663398743         STATISTICS_STDDEV=2.28225356537         title=rain_days         units=days     Band 3 Block=39x1 Type=Float32, ColorInterp=Undefined       Min=9.700 Max=21.762        Minimum=9.700, Maximum=21.762, Mean=13.882, StdDev=2.988       NoData Value=1.0000000200408773e+020       Metadata:         comment=Converted from UKCP09 monthly data to climatologies and then averaged to 25km rotated pole grid.         coordinates=lon lat         long_name=Number of days per meaning period having a rainfall >= 1 mm (days)         missing_value=1e+020         name=rain_days         NETCDF_DIM_meaning_period=2         NETCDF_VARNAME=raindays_tmean_abs         STATISTICS_MAXIMUM=21.761842727661         STATISTICS_MEAN=13.88236379732         STATISTICS_MINIMUM=9.7001943588257         STATISTICS_STDDEV=2.9875935432638         title=rain_days         units=days     Band 4 Block=39x1 Type=Float32, ColorInterp=Undefined       NoData Value=1.0000000200408773e+020       Metadata:         comment=Converted from UKCP09 monthly data to climatologies and then averaged to 25km rotated pole grid.         coordinates=lon lat         long_name=Number of days per meaning period having a rainfall >= 1 mm (days)         missing_value=1e+020         name=rain_days         NETCDF_DIM_meaning_period=3         NETCDF_VARNAME=raindays_tmean_abs         title=rain_days         units=days     Band 5 Block=39x1 Type=Float32, ColorInterp=Undefined       NoData Value=1.0000000200408773e+020       Metadata:         comment=Converted from UKCP09 monthly data to climatologies and then averaged to 25km rotated pole grid.         coordinates=lon lat         long_name=Number of days per meaning period having a rainfall >= 1 mm (days)         missing_value=1e+020         name=rain_days         NETCDF_DIM_meaning_period=4         NETCDF_VARNAME=raindays_tmean_abs         title=rain_days         units=days     Band 6 Block=39x1 Type=Float32, ColorInterp=Undefined       NoData Value=1.0000000200408773e+020       Metadata:         comment=Converted from UKCP09 monthly data to climatologies and then averaged to 25km rotated pole grid.         coordinates=lon lat         long_name=Number of days per meaning period having a rainfall >= 1 mm (days)         missing_value=1e+020         name=rain_days         NETCDF_DIM_meaning_period=5         NETCDF_VARNAME=raindays_tmean_abs         title=rain_days         units=days     Band 7 Block=39x1 Type=Float32, ColorInterp=Undefined       NoData Value=1.0000000200408773e+020       Metadata:         comment=Converted from UKCP09 monthly data to climatologies and then averaged to 25km rotated pole grid.         coordinates=lon lat         long_name=Number of days per meaning period having a rainfall >= 1 mm (days)         missing_value=1e+020         name=rain_days         NETCDF_DIM_meaning_period=6         NETCDF_VARNAME=raindays_tmean_abs         title=rain_days         units=days     Band 8 Block=39x1 Type=Float32, ColorInterp=Undefined       NoData Value=1.0000000200408773e+020       Metadata:         comment=Converted from UKCP09 monthly data to climatologies and then averaged to 25km rotated pole grid.         coordinates=lon lat         long_name=Number of days per meaning period having a rainfall >= 1 mm (days)         missing_value=1e+020         name=rain_days         NETCDF_DIM_meaning_period=7         NETCDF_VARNAME=raindays_tmean_abs         title=rain_days         units=days     Band 9 Block=39x1 Type=Float32, ColorInterp=Undefined       NoData Value=1.0000000200408773e+020       Metadata:         comment=Converted from UKCP09 monthly data to climatologies and then averaged to 25km rotated pole grid.         coordinates=lon lat         long_name=Number of days per meaning period having a rainfall >= 1 mm (days)         missing_value=1e+020         name=rain_days         NETCDF_DIM_meaning_period=8         NETCDF_VARNAME=raindays_tmean_abs         title=rain_days         units=days     Band 10 Block=39x1 Type=Float32, ColorInterp=Undefined       NoData Value=1.0000000200408773e+020       Metadata:         comment=Converted from UKCP09 monthly data to climatologies and then averaged to 25km rotated pole grid.         coordinates=lon lat         long_name=Number of days per meaning period having a rainfall >= 1 mm (days)         missing_value=1e+020         name=rain_days         NETCDF_DIM_meaning_period=9         NETCDF_VARNAME=raindays_tmean_abs         title=rain_days         units=days     Band 11 Block=39x1 Type=Float32, ColorInterp=Undefined       NoData Value=1.0000000200408773e+020       Metadata:         comment=Converted from UKCP09 monthly data to climatologies and then averaged to 25km rotated pole grid.         coordinates=lon lat         long_name=Number of days per meaning period having a rainfall >= 1 mm (days)         missing_value=1e+020         name=rain_days         NETCDF_DIM_meaning_period=10         NETCDF_VARNAME=raindays_tmean_abs         title=rain_days         units=days     Band 12 Block=39x1 Type=Float32, ColorInterp=Undefined       NoData Value=1.0000000200408773e+020       Metadata:         comment=Converted from UKCP09 monthly data to climatologies and then averaged to 25km rotated pole grid.         coordinates=lon lat         long_name=Number of days per meaning period having a rainfall >= 1 mm (days)         missing_value=1e+020         name=rain_days         NETCDF_DIM_meaning_period=11         NETCDF_VARNAME=raindays_tmean_abs         title=rain_days         units=days     Band 13 Block=39x1 Type=Float32, ColorInterp=Undefined       NoData Value=1.0000000200408773e+020       Metadata:         comment=Converted from UKCP09 monthly data to climatologies and then averaged to 25km rotated pole grid.         coordinates=lon lat         long_name=Number of days per meaning period having a rainfall >= 1 mm (days)         missing_value=1e+020         name=rain_days         NETCDF_DIM_meaning_period=12         NETCDF_VARNAME=raindays_tmean_abs         title=rain_days         units=days     Band 14 Block=39x1 Type=Float32, ColorInterp=Undefined       NoData Value=1.0000000200408773e+020       Metadata:         comment=Converted from UKCP09 monthly data to climatologies and then averaged to 25km rotated pole grid.         coordinates=lon lat         long_name=Number of days per meaning period having a rainfall >= 1 mm (days)         missing_value=1e+020         name=rain_days         NETCDF_DIM_meaning_period=13         NETCDF_VARNAME=raindays_tmean_abs         title=rain_days         units=days     Band 15 Block=39x1 Type=Float32, ColorInterp=Undefined       NoData Value=1.0000000200408773e+020       Metadata:         comment=Converted from UKCP09 monthly data to climatologies and then averaged to 25km rotated pole grid.         coordinates=lon lat         long_name=Number of days per meaning period having a rainfall >= 1 mm (days)         missing_value=1e+020         name=rain_days         NETCDF_DIM_meaning_period=14         NETCDF_VARNAME=raindays_tmean_abs         title=rain_days         units=days     Band 16 Block=39x1 Type=Float32, ColorInterp=Undefined       NoData Value=1.0000000200408773e+020       Metadata:         comment=Converted from UKCP09 monthly data to climatologies and then averaged to 25km rotated pole grid.         coordinates=lon lat         long_name=Number of days per meaning period having a rainfall >= 1 mm (days)         missing_value=1e+020         name=rain_days         NETCDF_DIM_meaning_period=15         NETCDF_VARNAME=raindays_tmean_abs         title=rain_days         units=days     Band 17 Block=39x1 Type=Float32, ColorInterp=Undefined       NoData Value=1.0000000200408773e+020       Metadata:         comment=Converted from UKCP09 monthly data to climatologies and then averaged to 25km rotated pole grid.         coordinates=lon lat         long_name=Number of days per meaning period having a rainfall >= 1 mm (days)         missing_value=1e+020         name=rain_days         NETCDF_DIM_meaning_period=16         NETCDF_VARNAME=raindays_tmean_abs         title=rain_days         units=days      **gdalinfo NETCDF:RainDays1_1961-1990_LTA_25km.nc:lon**               Driver: netCDF/Network Common Data Format     Files: RainDays1_1961-1990_LTA_25km.nc            RainDays1_1961-1990_LTA_25km.nc.aux.xml     Size is 39, 52     Coordinate System is `'     Origin = (-17.650008352179277,11.329999312758446)     Pixel Size = (0.219999614514803,-0.219999998807907)     Metadata:       axis_0#units=       lon#long_name=longitude       lon#missing_value=1e+020       lon#name=lon       lon#standard_name=longitude       lon#units=degrees_east       meaning_period#comment=Meaning periods are all months (12) beginning January, all seasons (4) beginning Winter (DJF) and Annual (1)       meaning_period#long_name=meaning_period       meaning_period#units=       NC_GLOBAL#Conventions=CF-1.0       rlat#axis=Y       rlat#bounds=bounds_rlat       rlat#long_name=latitude in rotated pole grid       rlat#standard_name=grid_latitude       rlat#units=degrees       rlon#axis=X       rlon#bounds=bounds_rlon       rlon#long_name=longitude in rotated pole grid       rlon#modulo=360       rlon#standard_name=grid_longitude       rlon#topology=circular       rlon#units=degrees     Corner Coordinates:     Upper Left  ( -17.6500084,  11.3299993)      Lower Left  ( -17.6500084,  -0.1100006)      Upper Right (  -9.0700234,  11.3299993)      Lower Right (  -9.0700234,  -0.1100006)      Center      ( -13.3600159,   5.6099993)      Band 1 Block=39x1 Type=Float32, ColorInterp=Undefined       Min=-15.830 Max=3.672        Minimum=-15.830, Maximum=3.672, Mean=-5.246, StdDev=4.393       NoData Value=1.0000000200408773e+020       Metadata:         long_name=longitude         missing_value=1e+020         name=lon         NETCDF_VARNAME=lon         standard_name=longitude         STATISTICS_MAXIMUM=3.671746969223         STATISTICS_MEAN=-5.2457962955782         STATISTICS_MINIMUM=-15.830010414124         STATISTICS_STDDEV=4.3927366732513         units=degrees_east      **gdalinfo NETCDF:RainDays1_1961-1990_LTA_25km.nc:lat**               Driver: netCDF/Network Common Data Format     Files: RainDays1_1961-1990_LTA_25km.nc            RainDays1_1961-1990_LTA_25km.nc.aux.xml     Size is 39, 52     Coordinate System is `'     Origin = (-17.650008352179277,11.329999312758446)     Pixel Size = (0.219999614514803,-0.219999998807907)     Metadata:       axis_0#units=       lat#long_name=latitude       lat#missing_value=1e+020       lat#name=lat       lat#standard_name=latitude       lat#units=degrees_north       meaning_period#comment=Meaning periods are all months (12) beginning January, all seasons (4) beginning Winter (DJF) and Annual (1)       meaning_period#long_name=meaning_period       meaning_period#units=       NC_GLOBAL#Conventions=CF-1.0       rlat#axis=Y       rlat#bounds=bounds_rlat       rlat#long_name=latitude in rotated pole grid       rlat#standard_name=grid_latitude       rlat#units=degrees       rlon#axis=X       rlon#bounds=bounds_rlon       rlon#long_name=longitude in rotated pole grid       rlon#modulo=360       rlon#standard_name=grid_longitude       rlon#topology=circular       rlon#units=degrees     Corner Coordinates:     Upper Left  ( -17.6500084,  11.3299993)      Lower Left  ( -17.6500084,  -0.1100006)      Upper Right (  -9.0700234,  11.3299993)      Lower Right (  -9.0700234,  -0.1100006)      Center      ( -13.3600159,   5.6099993)      Band 1 Block=39x1 Type=Float32, ColorInterp=Undefined       Min=47.594 Max=60.806        Minimum=47.594, Maximum=60.806, Mean=54.186, StdDev=3.229       NoData Value=1.0000000200408773e+020       Metadata:         long_name=latitude         missing_value=1e+020         name=lat         NETCDF_VARNAME=lat         standard_name=latitude         STATISTICS_MAXIMUM=60.805950164795         STATISTICS_MEAN=54.185836173138         STATISTICS_MINIMUM=47.594329833984         STATISTICS_STDDEV=3.228754901642         units=degrees_north


Label features on hover/click with QGIS?

Simple question: Is there a built in/easy way to show labels for features based on a hover or click action with QGIS 2.0.1?


Cache greater than 500mb ArcGIS Server 10.2.1

Will this affect performance? Why 500? And how can I minimize the size of the cache? 2 service less than 500 MB will be faster than one with more that 500MB? Thanks


Which units are used for the Print Composer grid?

I changed the units on the general tab to inches. When I go to set up the grid under snapping, whats the units. From what I have read it seems as though they are map units, it that correct? I am a attempting to get set up a sheet to fit within en existing title block.


Priority labelling in geoserver

Good day can i ask something about geoserver 2.2.1? I have my maps then i want to see the largest cities before the small cities but i don't have population instead i have my PLACE column its not a numbers its a text example               NAME                PLACE          ADDRESS     Ayusan Norte        Village        VIGAN     Barangay VIII       Village        VIGAN       Then i want to create a big point to see all bigest cities How can i do that thanks.. God Bless !


Android Apps to edit WFS-T

Are there any mobile Android Apps where you can edit the Features served via a WFS-T?


Source for US parcel boundary data?

Where can I find shapefile (or other vector) data for US private parcel boundaries (such as for a plat map)?


How to sum up values of neighbouring polygons?

I hope that you can help me with the following problem: I have a vector (polygon) layer. I would like to add an attribute to the layer that - for every polygon - sums the values of a specific field of all its neighboring polygons. To give a more concrete example: I have a polygon layer of districts containing information on population. Now, for every district I would like to know how many people live in all its neighboring districts. As I have more than 300 districts I cannot do this by hand for each and every district. Is there any way do this more efficiently in QGIS?


OSM printing in QGIS-Web-Client

Printing in QGIS-Web-Client generates pdf with only QGIS project layers. I tried to add osm map in QGIS Composer but it doesn't work. It still prints only overlays. Is there any other way to print a map with osm baselayer?


how to run an aml in python

I am trying to run a aml in a python program I am writing and haven't been having much luck figuring out the correct way to run the aml. The aml requires workstation to be in grid mode. Any tips? I am only a novice in workstation... I'm running arcmap 10.0 Thanks, Matt


How to get the a FeatureLayer full extent?

Is there a way to get the extent of a FeatureLayer (with all its features, not only the ones on the map)?


Geoserver: Georectification via SLD fails

## **Situation** I have a raster image and the coordinates of its four corner pixels. My goal is to warp the image with respect to this coordinates ("georectify") **on the fly** when requested via GeoServer to display in OpenLayers. Inspired by this discussion list, my idea was to do the following in GeoServer:   1. Create a WorldImage layer (raster + text file with geodata)   2. Create a Styled Layer Descriptor (SLD) calling the Web Processing Service (WPS) `GeorectifyCoverage` which "georectifies a raster via Ground Control Points using gdal_warp".   3. Apply the SLD to the layer. ## **Problem** When I request the layer for example via the OpenLayers preview, it throws an error reporting:               Could not convert the value EPSG:32632 into the expected type interface     org.opengis.referencing.crs.CoordinateReferenceSystem for parameter targetCRS      Oddly enough, when I use the "WPS request builder" (in the GeoServer Web Application) to call `GeorectifyCoverage` with the same parameters, it gives the right output. So the error must originate somewhere from my SLD document below. **Can anybody give me a clue what might be wrong with it?** * * * ## **Additional Info** SLD document calling WPS rendering transformation `GeorectifyCoverage`:               <?xml version="1.0" encoding="ISO-8859-1"?>     <StyledLayerDescriptor version="1.0.0" xsi:schemaLocation="http://www.opengis.net/sld StyledLayerDescriptor.xsd" xmlns="http://www.opengis.net/sld" xmlns:ogc="http://www.opengis.net/ogc" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance">         <NamedLayer>             <Name>GeorectifyCoverage</Name>             <UserStyle>                 <Title>Georectify Coverage</Title>                 <Abstract>Warps image to specified GCP</Abstract>                 <FeatureTypeStyle>                     <Transformation>                         <ogc:Function name="gs:GeorectifyCoverage">                             <ogc:Function name="parameter">                                 <ogc:Literal>data</ogc:Literal>                             </ogc:Function>                             <ogc:Function name="parameter">                                 <ogc:Literal>gcp</ogc:Literal>                                 <ogc:Literal>[[[0, 0], [308639.500, 5154180.500]], [[0, 450], [308639.500, 5150929.500]], [[600, 0], [315730.500, 5154180.500]], [[600, 450], [315730.500, 5150929.500]]]</ogc:Literal>                             </ogc:Function>                             <ogc:Function name="parameter">                                 <ogc:Literal>targetCRS</ogc:Literal>                                 <ogc:Literal>EPSG:32632</ogc:Literal>                             </ogc:Function>                            </ogc:Function>                     </Transformation>                     <Rule>                         <RasterSymbolizer>                             <Opacity>1.0</Opacity>                         </RasterSymbolizer>                     </Rule>                 </FeatureTypeStyle>             </UserStyle>         </NamedLayer>     </StyledLayerDescriptor>      where `GeorectifyCoverage` requires the following input:   * `data`: input raster (my layer name)   * `gcp`: List of Ground control points (my corner coordinates)   * `targetCRS`: system to use for the output raster (my CRS corresponding to the GCP) **Notes:**   * When I call a different WPS function like `Contour` using the same document structure and raster image it works perfectly.   * When I change `EPSG:32632` to `EPSG:4326` the same error is reported.   * Changing the raster format (jpeg, TIFF) does not helo.   * Processing the same data offline with the command line tools `gdal_translate` (assign GCPs) and `gdalwarp` results in the desired output. **My System:**   * Ubuntu 12.04   * OpenGeo Suite 3.0.2   * GeoServer 2.2 Alternative approaches using OpenGeo would also be welcome.


How to search text string and select item in QGIS?

I have a table attribute with complete roads name. I have to search items by digiting a part of the road name and I need to colour all the items selected. eg: ROAD NAME: / road 01 / road 02 / road 11 / road 12 I need to colour all items that contain "1" (so all the items except the second one). I would like to set the string operator is this way: ROAD NAME ILIKE ' * 1 * ' but it does not work _edit_ I work in advanced search in table attribute


Is it possible to undertake spatial clustering with QGIS?

Is it possible to undertake spatial clustering with QGIS, and also end up with proportional symbols in a quasi 3-d symbology?


Openlayers : Zoom In and Map Layers Blocked wih WHITE space

I've a problem where something is blocking my map layer when I zoomed the map click to view               function load() {         map = new OpenLayers.Map({             div: "map",             maxExtent: new OpenLayers.Bounds(101.076797485352, 3.60602392425537,101.077194213867, 3.60696087346802)                  });         // setup tiled layer         floorplan7 = new OpenLayers.Layer.WMS(                     "assetmapper:floor_plan_t7 - Tiled", "http://192.168.56.101:8080/geoserver/assetmapper/wms",                     {                         LAYERS: 'assetmapper:floor_plan_t7',                         STYLES: '',                         format: 'image/png',                          },                     {                         buffer: 0,                     //    displayOutsideMaxExtent: true,                         isBaseLayer: true,                         yx : {'EPSG:4326' : true}                     }              );               map.addLayers([floorplan7]);         map.addControl(new OpenLayers.Control.LayerSwitcher());         map.addControl(new OpenLayers.Control.ScaleLine({}));         map.zoomToMaxExtent();          }      What should I do to clear the 'WHITE' area?


Quantum GIS Installation Failure

I installed GSL 1.16 & GDAL 1.10, then I proceed to install QGIS 1.8.0 After running the .pkg(installer) I get the following message. * * * Failure The GDAL 1.9 framework is required. * * * Does this mean I should be using GDAL 1.9 instead of 1.10 or is it just not recognizing 1.10?


Problem changing CRS to enable clip function

I have two layers, one point file and one polygon. The point file is in WGS 84 and polygon in NAD83-UTM10. The project CRS is NAD83-UTM10. They project fine with on-the-fly projection turned on. I want to clip the point file with the polygon, and understand that all layers need to be in the same CRS to enable clipping or other ftools functions. When I change the CRS of either layer to the other (changing the point file to NAD83 for example), it projects far away from where it was (ie. doesn't line up). This is the same if I right click to change the CRS or if I "save as." Same no matter which CRS I select. I'm running QGIS 2.0.1 Dufour on Mac. Any help would be greatly appreciated! Thanks!


How to create a new table for a vector file?

I will use `v.db.connect` but i want to use a new table, and not an existing table to connect it to. What is the GRASS script for it? Thanks!


Restart ArcGIS Server Map Service From Command Line

Is there a way to restart an ArcGIS Server v10 Map Service from the command line in Windows 7?


PyQGIS: How to color raster layer with existing colorBrewerMap

With pyqgis I have added multiple new raster layers to the legendinterface().   Because i like the ColorBrewer style i want to color some of them with the " **Spectral** " gradient.   How can I do this in python? The API-help only gives advice how to color with individual set colorramps.   I found out that I can query the existing Ramps with `QgsColorBrewerPalette.listSchemes()`, but how can I set this style for my layer?                 qshade = QgsColorRampShader()     qshade.setColorRampType(QgsColorRampShader.INTERPOLATED)     qshade.setColorRampItemList("Spectral") # Similar to this?


ArcGIS: recode to different geography

What is the easiest way to recode data from one geography to another? For example, if I have population data for one set of polygons, and I'd like to deduce it for different polygons - a bit like a spatial join but taking account of partial overlaps. For this job I am willing to assume even distribution of population over area. I imagine you could produce a raster then use zonal statistics to convert back to polygons, but is there a quicker way?


Can I avoid decimal conversion calling CostDistance

Our client asked us to incorporate blockages and delays for movement of boats and ships on the local waterways. In our one process we are creating polygons of the distance a boat can move at top speed for up to 90 minutes. The change required us to look at several things and we decided to set the raster to a time-based seconds per meter. In our test sample I have a point located on the waterway. I have time zones established for how far a boat will move with no slowdowns.               ie.  Boat speed = 40 knots          @15 minutes  10 nm          @30 minutes  20 nm          @45 minutes  30 nm          @60 minutes  40 nm          @90 minutes  60 nm     To define the seconds/meter we compute 40 * 1.94384449 rounded and converted to integer.      this yields a default raster value of 78 which is the number of seconds it takes for the boat to move one meter. And when we convert nautical miles to meters the maximum distance a boat can travel at full speed is 8667360 meters. All of these numbers are within a decent realm for storing as integer values in a cost raster, yet when I call CostDistance() it takes 33 hours to run and the result shows the 8667360 in scientific notation. (which means it is being stored as floating point in the raster) So my question is What can I do to force CostDistance to remain integer so it will not run like a dog and take 33 hours?


Create Image Mosaic JDBC datastore in command line

I would like to publish raster from PostGIS database, and I found this tutorial http://docs.geoserver.org/2.1.4/user/tutorials/imagemosaic- jdbc/imagemosaic-jdbc_tutorial.html , everything works fine, but I need to configure GeoServer in command line (curl for example). Is there an opportunity to do this? Thanks for the information. Jb


How to obtain timestamp and user id from .osm file in PostGIS enabled Postgres database

At present, osm2pgsql seems to ignore the timestamp and user id of the elements, a serious drawback for some applications. According to a bug report this is not solved by the addition of the `\--extra-attributes` tag. I've just tested this in osm2pgsql 0.83.0 and am wondering if there is a good way around it? My plan B is to do pre-filtering with osmosis and then load the data into R using osmar. Still, it would be nice to have the PostGIS option working also.


What will GIS be like in 5 years (and in 50)?

How will GIS technology be used 5 years from now ... and in 50 years?


how to add the label for located charts in ArcGIS

I am creating a Location Pie chart for my study area using Male and Female population in ArcGIS software, the software creates for legend of charts but I want to create a label on chart itself, how to do it in ArcGIS please help me ![sample map for locational chart](http://i.stack.imgur.com/CtHqf.png) ![the chart should have label like this on its location](http://i.stack.imgur.com/pwkaY.png)


Given a geocode or geo location in US, how to find exits that are within some distance from the geo location for US?

We have a need to find the geocodes of US highway exits that are within some distance from a given geocode. I do not know how to do this? Can one please direct to the right source or provide a solution to me?


Updating MapGuide Enterprise 2011's Query Tool

I am trying to create a custom theme for my company's website using MapGuide Enterprise 2011, and I am looking for advice on fine-tuning the software's built-in query tool to allow for more advanced querying. The built-in Query tool is (as far as I can tell) only able to perform a single query on a layer, however I am hoping to be able to query multiple layers/multiple queries in the same layer for more precise searching. Being able to save my query results as a selection would also be ideal. Any assistance/feedback would be greatly appreciated. Thank you.


How to open the Shape files and extract the data from these files?

We are in need of US highway data. We found it for Florida at the state GIS website . We downloaded the data for Florida highway intersections from here (this is a zip file) . We do not know how to use these files and extract the geocode values of the intersections (exits) from these files. Any help is greatly appreciated.


How to apply custom Google map style in OpenLayers?

Recently I stumbled upon interesting article showing potential of adding custom styles to Google Maps. Some examples of styled maps can be seen here and you can design your own from scratch here. I'd like to apply 'greyscale' theme to my OpenLayers base map: ![enter image description here](http://i.stack.imgur.com/PZCZU.jpg) Style description according to Google demo would look somehow like:               var styles = {       'Greyscale': [         {                         featureType: 'all',           rules: [             {saturation: -100},             {gamma: 0.50}           ]         }       ]     }      My Google map in Open layers currently looks simply like:               var gmap = new OpenLayers.Layer.Google("Google Streets",         {         'numZoomLevels': 20,         'sphericalMercator': true}     );      How could I apply greyscale theme to my map?


Why is my animation exported from ArcGlobe not playing?

The exporting process has finished and an avi file was created. When I try playing it, it raises an error in WMP as well as VLC. VLC gives me some information: "audio or video format not supported "undf" I am exporting with no compression to avi. Oh yeah, and unfortunately there is no way to fix this.... Any ideas? ![enter image description here](http://i.stack.imgur.com/fXkBu.png)


Why is my Polygon (WGS84) in SpatiaLite 4 not displayed by QGIS?

I made a polygon in a SpatiaLite 4 database as follows:               INSERT INTO test_pg (id, name, geom) VALUES (3, 'test', ST_PolyFromText('POLYGON((8.556824 47.185269, 8.611009 47.185269, 8.611009 47.218568, 8.556824 47.218568, 8.556824 47.185269))', 4326));      I assume that the order of coordinates is always LON LAT - is that correct? When I do `SELECT ST_AsText(geom) FROM test_pg` it correctly outputs the polygon coordinates. Now in QGIS, I add the SpatiaLite layer but the polygon does not get displayed (even after clicking on this zoom to layer thing). Geometry seems to be valid, I checked it with the QGIS tool. _I also found out that when moving the cursor on the map in the east direction, 47.xxx (LAT in this case) gets increased, when moving the cursor in the north direction, 8.xxx gets increased. This should be reversed, shouldn't it?_ Where are my errors? **EDIT** : I also INSERTED a polygon with LON LAT coordinates reversed. This polygon actually gets displayed in QGIS, but again, 47.xxx is displayed as LON and 8.xxx as LAT where it should be the opposite actually. _**EDIT2_** : Consider this question as closed. For some reasons, it suddenly works. Don't ask me why, maybe restarting QGIS several times helped...


How to get ArcGIS Templates in VS2010

I have a work computer which has ArcGIS desktop installed first and then VS2010 (just the way IT does it) and then templates are not available, anyone know what I need to do to see the templates? I am running this under XP.


QGis on Ubuntu: Install Dufour w/ Master already installed?

As per title really: I've got Master installed & updating from repos nightly (I assume). But it has bugs which directly affect the work I'm trying to do, so I also need Dufour to do some processing. Does anyone know if it's possible to have both installed concurrently, as it is in windows? Seem that "qgis" runs master...


GIS Server - evaluation architecture

Building from my previous questions: QGIS Server limitations? I want to evaluate performance of both MapServer and GeoServer. Can both servers exist on the same machine (each serving out WMS) or should I set them up on separate servers? thank you!


Weird behaviour of ST_Project with PostGIS2

somehow I can't understand the results of ST_Project, but maybe this problem is just me being stupid (again): Note:   * I use geography-type here, but the results are the same for geometry.   * I use a rather long distance, because I need, at least in one direction, for my application a large distance. Because of this I found the problem in the first place. So I just want to project a point towards North, East, South and West. With the direction north, everything is as I expect it:               SELECT ST_AsText(ST_Project('POINT(-30 -0.5)'::geography, 1e7, radians(0)));           st_astext               -----------------------------     POINT(-30 89.4874128664049)      But in the other directions, here for example East, I get this:               SELECT ST_AsText(ST_Project('POINT(-30 -0.5)'::geography, 1e7, radians(90)));                   st_astext                       ---------------------------------------------     POINT(59.8315233716443 0.00116689648502556)      Shouldn't this be `POINT(59.8315233716443 -0.5)` ? This behaviour occurs in all directions save north. I'm new to the whole earth coordinate system-stuff so maybe I missed something, but I thought that when one walks directly west (or south) one would, by definition, not leave the current longitude (or latitude)-line but only change the other. Any help would be appreciated! Thank you.


Merging many rasters using Python

I have a list of many many rasters. I want to have a simple and easy to understand script that takes this and merges it into one big mosaic. I found a similar question but the accepted answer adivises against using R for this purpose. I could also use R but are very inclined towards Python. I really need a script as I am going to simply share it with some people that want it in this form. I am inclined to using gdal_merge.py but I wanted some thoughts on the problem first as **I am potentially going to handle massive amounts of images**. Why merge so many images? I also found it impractical but was something I was asked to do.


Unable to load Shapefile using pgShapeloader / shp2pgsql-gui (pgAdmin3's PostGIS Shapefile Import/Export Tool)

I am trying to load a shapefile (http://dds.cr.usgs.gov/pub/data/nationalatlas/statesp020_nt00032.tar.gz) into PostGIS 2.0.1 using pgAdmin3's PostGIS Shapefile Import/Export Tool. ![enter image description here](http://i.stack.imgur.com/yh7bA.png) However its not allowing me to import the shapefile, showing the error messages below. Any idea what happen? **Output from PostGIS Shapefile Import/Export Tool**               ==============================     Importing with configuration: statesp020, staging, geom, C:\gisdata\statesp020_nt00032\statesp020, mode=c, dump=1, simple=0, geography=0, index=0, shape=1, srid=4269     Shapefile type: Polygon     PostGIS type: MULTIPOLYGON[2]     Failed SQL begins: "SET CLIENT_ENCODING TO UTF8;     SET STANDARD_CONFORMING_STRINGS TO ON;     BEGIN;     CREATE TABLE "staging"."statesp020" (gid serial,     "area" float8,     "perimeter" float8,     "statesp020" float8,     "state" varchar(20),     "state_fips" varchar(2),     "order_adm" int4,     "month_adm""     Failed in pgui_exec(): ERROR:  function addgeometrycolumn(unknown, unknown, unknown, unknown, unknown, integer) does not exist     LINE 15: SELECT AddGeometryColumn('staging','statesp020','geom','4269...                     ^     HINT:  No function matches the given name and argument types. You might need to add explicit type casts.          Shapefile import failed.


how to use mapproxy to serve WMS from reprojected google maps tiles, i have searched for docs, but there's no docs which describe how to add google maps as source.


How to automate export in Tilemill?

I have installed Tilemill to export my map as mbtiles,and I used mbutil to export as PNG images & UTFGrid.So forth,I have a Tilemill project with MSS(Carto) files.Now the data source is going to change weekly and I wish to automate this action. So does Tilemill support command-line scripting or something else to export my project with given properties ?? Current System : Windows 7 x64 **EDIT in response to @Sunil:** > TileMill doesn't have any features for automatically manipulating your Carto > stylesheets, but you could use a bash script to manipulate the Carto files > directly. TileMill also has an export command that can be called from bash. On Mac OS X, you could run a command like this in the Terminal for more details about using the export command: /Applications/TileMill.app/Contents/Resources/node /Applications/TileMill.app/Contents/Resources/index.js export --help **The exact command depends on how you installed TileMill** Unfortunately,same command for Mac OS X does not work on my current system.What does the last sentence mean,is there different installations of Tilemill ?


Output of RasterToPolygon not supported as shapefile

I made a python code for arcGis 9.3 which do the following step : > 1) Convert all the raster (.tif) in a folder to polygon shapefile in another > folder > > 2) Intersect all the converted shapefile polygon (previous) with a mask The problem is that the intersect function doesn't work and it returns shapefile with no populate fields. So I try manually to use the intersect Tool with the converted polygon create by my python code. And I have this error when trying to put the converted polygon shape in input : PolygonName # Does not exist or is not supported And when I Copy/paste this same polygon in another folder, then it works. My code is as following :               import sys, string, os, arcgisscripting, glob     import arcgisscripting     gp = arcgisscripting.create()          # Set a default workspace     gp.workspace = "Z:\PPR\PPR76\PPRI_AUSTREBERTHE\Modelisation_Ruissellements\Tentative Macro Amelioration MNT"          # Set a toolbox     gp.toolbox = "management"               pathImages = gp.GetParameterAsText(0)     pathOutputConversion = gp.GetParameterAsText(1)     pathOutputIntersect = gp.GetParameterAsText(2)     masque = gp.GetParameterAsText(3)          # ON CONVERTI TOUTES LES IMAGES DU DOSSIER EN POLYGON          for infile in glob.glob(os.path.join(pathImages, '*.tif')):         flnm = os.path.basename(infile)         # Avec cette fonction : fileName contient le CHEMIN+nom du fichier sans l'estension et fileExtension contient l'extension (.qqch). Ici on voulait huste le nom du fichier sans le chemin donc on l'a applique a flnm          fileName, fileExtension = os.path.splitext(flnm)              # Local variables...         Input_Raster= os.path.join(pathImages, infile)           Output_Polygon_Features = pathOutputConversion + "\ " + fileName+ ".shp "              # Process: RasterToPolygon...         gp.RasterToPolygon_conversion(Input_Raster, Output_Polygon_Features, "SIMPLIFY", "")          gp.AddMessage(" FIN DE CONVERSION RASTER TO POLYGON ")          # ON INTERSECTE TOUS LES POLYGON IMAGES AVEC LE MASQUE DU BV          for infile in glob.glob(os.path.join(pathOutputConversion, '*.shp')):       flnm = os.path.basename(infile)       # Local variables...       Output_Feature_Class = pathOutputIntersect + "\ " +flnm       Input_Features = masque +" #;"+os.path.join(pathOutputConversion, infile)+" #"            # Process: Intersect...       gp.Intersect_analysis(Input_Features, Output_Feature_Class, "ALL", "", "INPUT")      Thanks for your help,


Get selected features on ArcMap in arcMap Addin (eclipse)

I am building a custom arcMap Addin button which onclick event tells info about selected features. Following is my java class ---               public class Break extends Button {          /**      * Called when the button is clicked.      *       * @exception java.io.IOException if there are interop problems.      * @exception com.esri.arcgis.interop.AutomationException if the component throws an ArcObjects exception.      */     @Override     public void onClick() throws IOException, AutomationException {         // TODO Auto-generated method stub                   MapBean mapBean = new MapBean();         try{             IMap map = mapBean.getMap();             ISelection selection = map.getFeatureSelection();             IEnumFeature enumFeature = (IEnumFeature)selection;             enumFeature.reset();             IFeature feature = enumFeature.next();             while (feature != null){                 System.out.println("FEATURE ID: " + feature.getOID());                 feature = enumFeature.next();             }         }         catch (Exception e){             e.printStackTrace();         }         }     }         I implimented the button on ArcMap but nothing seems to work onclick event , might be its not getting the ArcMap at all.Could anybody tell how to resolve this issue. Thanks in advance////


Get latest tiles from TileMill's tile server in Leaflet

My leaflet app gets the tiles from TileMill's tile server in the standard way:               L.tileLayer('http://<myip>/tile/SteveCountryVic/{z}/{x}/{y}.png', { ... });      These tiles are often old, not showing the latest styles. I'm not sure if they're coming from TileMill itself or some intermediary cache. I notice that TileMill includes an 'updated' parameter:               http://<myip>/tile/SteveCountryVic/{z}/{x}/{y}.png?updated=1372394902000      Manually copying that updated timestamp means I get the latest timestamp. But what's a good way to not have to manually update the URL every time I change the style?


Overlay wms + google layer in Openlayers

I want to use OpenLayers, GeoServer layer with the google maps api as a base. In my GeoServer layer is redesigned to EPSG: 900913 when I open the geoserver openlayers with everything ok. ![Layer Configuration](http://i.stack.imgur.com/JkEMa.png) On a map created in OpenLayers projection EPSG: 900913 for working with the Google Maps API, then all my layers took this projection. Not forgetting that the GeoServer natively redesigns layer as configured above. Tested with alloverlay property to see where he was going to stop the WMS layer but she even appears. Another attempt, instead of using the default maxResolution bounds and used them in the openlayers geoserver, open the layer in openlayers with geoserver teclei F12 and copied the properties and bounds maxResolution to my map, it did not work. The question is: why the WMS layer does not appear? My code:                   <!DOCTYPE html>     <html lang='en'>     <head>      <meta charset='utf-8' />      <title>My OpenLayers Map</title>      <script type='text/javascript' src='OpenLayers.js'></script>           <script src="http://maps.google.com/maps/api/js?sensor=false&v=3.2"></script>           <script type='text/javascript'>                var map;      var mercatorProjection = new OpenLayers.Projection('EPSG:900913');      var latLongProjection = new OpenLayers.Projection('EPSG:4326');      var bounds = new OpenLayers.Bounds(                     -51.2809219360352, -30.2445888519287,                     -51.0207977294922, -29.9661273956299                 );           function init() {            map = new OpenLayers.Map('map_element',{         maxExtent: bounds,         maxResolution:  0.0010877400636672,         units: 'm',         //allOverlays: true,         projection: mercatorProjection,         displayProjection: latLongProjection      });                var google_streets = new OpenLayers.Layer.Google(         "Ruas",         {numZoomLevels: 20}         );          var wms_layer = new OpenLayers.Layer.WMS(     'WMS',     'http://localhost:8080/geoserver/Teste/wms',     {layers: 'Teste:onibus', transparent: true},     {isBaseLayer: false,         opacity: 0.7}     );               //Adiciona as camadas ao mapa     map.addLayers([google_streets, wms_layer]);               var point = new OpenLayers.LonLat(-51.01,-30.01);      point.transform(new OpenLayers.Projection("EPSG:4326"),          map.getProjectionObject());      map.setCenter(point, 10);       //Camada de controle que vai mostrar as camadas no mapa      map.addControl(new OpenLayers.Control.LayerSwitcher({}));           //Verifica se o mapa tem um ponto central e o extende a sua extensão máxima      if(!map.getCenter()){      map.zoomToMaxExtent();           }      }           </script>     </head>          <body onload='init();'> <!-- Chama a função js init() -->           <!--Elemento HTML onde o mapa é exibido -->      <div id='map_element' style='width: 800px; height: 800px;'>      </div>          </body>     </html>      Another thing I picked up the WMS Network link on this same layer and put the overlay on Earth was perfect. So my layer is properly designed . My page result: ![page](http://i.stack.imgur.com/aEi59.png) Link geoserver:               http://localhost:8080/geoserver/Teste/wms?service=WMS&version=1.1.0&request=GetMap&layers=Teste:onibus&styles=&bbox=-51.2809219360352,-30.2445888519287,-51.0207977294922,-29.9661273956299&width=478&height=512&srs=EPSG:900913&format=application/openlayers      Code openlayers on geoserver:                                                            <script>             var map;             var untiled;             var tiled;             var pureCoverage = false;             // pink tile avoidance             OpenLayers.IMAGE_RELOAD_ATTEMPTS = 5;             // make OL compute scale according to WMS spec             OpenLayers.DOTS_PER_INCH = 25.4 / 0.28;                  function init(){                 // if this is just a coverage or a group of them, disable a few items,                 // and default to jpeg format                 format = 'image/png';                 if(pureCoverage) {                     document.getElementById('filterType').disabled = true;                     document.getElementById('filter').disabled = true;                     document.getElementById('antialiasSelector').disabled = true;                     document.getElementById('updateFilterButton').disabled = true;                     document.getElementById('resetFilterButton').disabled = true;                     document.getElementById('jpeg').selected = true;                     format = "image/jpeg";                 }                      var bounds = new OpenLayers.Bounds(                     -51.2809219360352, -30.2445888519287,                     -51.0207977294922, -29.9661273956299                 );                 var options = {                     controls: [],                     maxExtent: bounds,                     maxResolution: 0.0010877400636672,                     projection: "EPSG:900913",                     units: 'm'                 };                 map = new OpenLayers.Map('map', options);                      // setup tiled layer                 tiled = new OpenLayers.Layer.WMS(                     "Teste:onibus - Tiled", "http://localhost:8080/geoserver/Teste/wms",                     {                         LAYERS: 'Teste:onibus',                         STYLES: '',                         format: format,                         tiled: true,                         tilesOrigin : map.maxExtent.left + ',' + map.maxExtent.bottom                     },                     {                         buffer: 0,                         displayOutsideMaxExtent: true,                         isBaseLayer: true,                         yx : {'EPSG:900913' : false}                     }                  );                      // setup single tiled layer                 untiled = new OpenLayers.Layer.WMS(                     "Teste:onibus - Untiled", "http://localhost:8080/geoserver/Teste/wms",                     {                         LAYERS: 'Teste:onibus',                         STYLES: '',                         format: format                     },                     {                        singleTile: true,                         ratio: 1,                         isBaseLayer: true,                        yx : {'EPSG:900913' : false}                     }                  );                      map.addLayers([untiled, tiled]);                      // build up all controls                 map.addControl(new OpenLayers.Control.PanZoomBar({                     position: new OpenLayers.Pixel(2, 15)                 }));                 map.addControl(new OpenLayers.Control.Navigation());                 map.addControl(new OpenLayers.Control.Scale($('scale')));                 map.addControl(new OpenLayers.Control.MousePosition({element: $('location')}));                 map.zoomToExtent(bounds);                      // wire up the option button                 var options = document.getElementById("options");                 options.onclick = toggleControlPanel;                      // support GetFeatureInfo                 map.events.register('click', map, function (e) {                     document.getElementById('nodelist').innerHTML = "Loading... please wait...";                     var params = {                         REQUEST: "GetFeatureInfo",                         EXCEPTIONS: "application/vnd.ogc.se_xml",                         BBOX: map.getExtent().toBBOX(),                         SERVICE: "WMS",                         INFO_FORMAT: 'text/html',                         QUERY_LAYERS: map.layers[0].params.LAYERS,                         FEATURE_COUNT: 50,                         Layers: 'Teste:onibus',                         WIDTH: map.size.w,                         HEIGHT: map.size.h,                         format: format,                         styles: map.layers[0].params.STYLES,                         srs: map.layers[0].params.SRS};                          // handle the wms 1.3 vs wms 1.1 madness                     if(map.layers[0].params.VERSION == "1.3.0") {                         params.version = "1.3.0";                         params.j = parseInt(e.xy.x);                         params.i = parseInt(e.xy.y);                     } else {                         params.version = "1.1.1";                         params.x = parseInt(e.xy.x);                         params.y = parseInt(e.xy.y);                     }                          // merge filters                     if(map.layers[0].params.CQL_FILTER != null) {                         params.cql_filter = map.layers[0].params.CQL_FILTER;                     }                      if(map.layers[0].params.FILTER != null) {                         params.filter = map.layers[0].params.FILTER;                     }                     if(map.layers[0].params.FEATUREID) {                         params.featureid = map.layers[0].params.FEATUREID;                     }                     OpenLayers.loadURL("http://localhost:8080/geoserver/Teste/wms", params, this, setHTML, setHTML);                     OpenLayers.Event.stop(e);                 });             }                  // sets the HTML provided into the nodelist element             function setHTML(response){                 document.getElementById('nodelist').innerHTML = response.responseText;             };                  // shows/hide the control panel             function toggleControlPanel(event){                 var toolbar = document.getElementById("toolbar");                 if (toolbar.style.display == "none") {                     toolbar.style.display = "block";                 }                 else {                     toolbar.style.display = "none";                 }                 event.stopPropagation();                 map.updateSize()             }                  // Tiling mode, can be 'tiled' or 'untiled'             function setTileMode(tilingMode){                 if (tilingMode == 'tiled') {                     untiled.setVisibility(false);                     tiled.setVisibility(true);                     map.setBaseLayer(tiled);                 }                 else {                     untiled.setVisibility(true);                     tiled.setVisibility(false);                     map.setBaseLayer(untiled);                 }             }                  // Transition effect, can be null or 'resize'             function setTransitionMode(transitionEffect){                 if (transitionEffect === 'resize') {                     tiled.transitionEffect = transitionEffect;                     untiled.transitionEffect = transitionEffect;                 }                 else {                     tiled.transitionEffect = null;                     untiled.transitionEffect = null;                 }             }                  // changes the current tile format             function setImageFormat(mime){                 // we may be switching format on setup                 if(tiled == null)                   return;                      tiled.mergeNewParams({                     format: mime                 });                 untiled.mergeNewParams({                     format: mime                 });                 /*                 var paletteSelector = document.getElementById('paletteSelector')                 if (mime == 'image/jpeg') {                     paletteSelector.selectedIndex = 0;                     setPalette('');                     paletteSelector.disabled = true;                 }                 else {                     paletteSelector.disabled = false;                 }                 */             }                  // sets the chosen style             function setStyle(style){                 // we may be switching style on setup                 if(tiled == null)                   return;                      tiled.mergeNewParams({                     styles: style                 });                 untiled.mergeNewParams({                     styles: style                 });             }                  // sets the chosen WMS version             function setWMSVersion(wmsVersion){                 // we may be switching style on setup                 if(wmsVersion == null)                   return;                      if(wmsVersion == "1.3.0") {                    origin = map.maxExtent.bottom + ',' + map.maxExtent.left;                 } else {                    origin = map.maxExtent.left + ',' + map.maxExtent.bottom;                 }                      tiled.mergeNewParams({                     version: wmsVersion,                     tilesOrigin : origin                 });                 untiled.mergeNewParams({                     version: wmsVersion                 });             }                  function setAntialiasMode(mode){                 tiled.mergeNewParams({                     format_options: 'antialias:' + mode                 });                 untiled.mergeNewParams({                     format_options: 'antialias:' + mode                 });             }                  function setPalette(mode){                 if (mode == '') {                     tiled.mergeNewParams({                         palette: null                     });                     untiled.mergeNewParams({                         palette: null                     });                 }                 else {                     tiled.mergeNewParams({                         palette: mode                     });                     untiled.mergeNewParams({                         palette: mode                     });                 }             }                  function setWidth(size){                 var mapDiv = document.getElementById('map');                 var wrapper = document.getElementById('wrapper');                      if (size == "auto") {                     // reset back to the default value                     mapDiv.style.width = null;                     wrapper.style.width = null;                 }                 else {                     mapDiv.style.width = size + "px";                     wrapper.style.width = size + "px";                 }                 // notify OL that we changed the size of the map div                 map.updateSize();             }                  function setHeight(size){                 var mapDiv = document.getElementById('map');                      if (size == "auto") {                     // reset back to the default value                     mapDiv.style.height = null;                 }                 else {                     mapDiv.style.height = size + "px";                 }                 // notify OL that we changed the size of the map div                 map.updateSize();             }                  function updateFilter(){                 if(pureCoverage)                   return;                      var filterType = document.getElementById('filterType').value;                 var filter = document.getElementById('filter').value;                      // by default, reset all filters                 var filterParams = {                     filter: null,                     cql_filter: null,                     featureId: null                 };                 if (OpenLayers.String.trim(filter) != "") {                     if (filterType == "cql")                          filterParams["cql_filter"] = filter;                     if (filterType == "ogc")                          filterParams["filter"] = filter;                     if (filterType == "fid")                          filterParams["featureId"] = filter;                 }                 // merge the new filter definitions                 mergeNewParams(filterParams);             }                  function resetFilter() {                 if(pureCoverage)                   return;                      document.getElementById('filter').value = "";                 updateFilter();             }                  function mergeNewParams(params){                 tiled.mergeNewParams(params);                 untiled.mergeNewParams(params);             }         </script>


Using arcpy.Describe Method on Annotation layers

I'm testing layers in the current document using the `Describe` function. My code is:               for lyr in arcpy.mapping.ListLayers(df):          if(lyr.isGroupLayer == False and lyr.isRasterLayer == False):              x = arcpy.Describe(str(lyr))      It's throwing an error on the annotation layer: > "LSDSectionBoundaryAnno\Default" does not exist I'm thinking that the Describe function doesn't like annotation layers or the subclasses? Can I use Describe on annotation layers and if not, how would I filter out annotation layers? There doesn't seem to be a 'isAnnotationLayer' property. Thanks, Mike