Rotating Map Elements in ArcGIS 10.0

I have imported a shapefile whole polygon is diagonally oriented. Rather than rotate the polygon, I would like to adjust the scale bar/legend/title to fit this unique orientation. Would like to know it both through GUI and python.


Are there any alternatives to ArcScene?

I think I've had it. I've been struggling to create polished scenes and flyover animations using the ESRI tools. SketchUp is great but it does not handle large datasets and most 3rd party models (3D DXF, COLLADA) do not retain spatial reference when imported into SketchUp. I've been able to build fairly large scenes in SketchUp in the past after I prepared my GIS data in ArcGIS 9.x and then used the Sketchup Plugin to export my features. Adding new features was very easy simply export them to SketchUp then copy the model from a new SketchUp project and paste into place in my primary model. This is no longer possible in ArcGIS 10 and the "Collada" method; files exported from ArcGIS are centered at 0,0,0 when opened in SketchUp. The existing "ESRI" workflow seems to be designed only around improving the 3D characteristics of simple extruded multipatch features in SketchUp but the ArcGIS --> SketchUp workflow via the legacy Google plugin has been obliterated. I am not looking for a freeware but rather **a solid, polished software that is capable of producing large 3D scenes from GIS data**. My **budget is up to $8,000** or less. What are alternatives are there? I have fairly large datasets; Usually LiDAR surface and my scenes need to be very large and detailed so this software needs to handle large datasets and also need to be able to retain spatial reference in files imported from 3D DXF, Collada, etc.


Which GIS Software is Right for small manufacturing trends analysis?

I'm a graphic designer who has been tasked with choosing a mapping platform for the manufacturing company I work for. One of my colleagues began our project by using a trial version of ArcGIS to develop some relatively simple maps - just a few shape layers using vectors to describe trade regions, zones associated with sales reps and the like - using light gray basemap with streets and political boundaries plus a population density map as a layer (from ESRI's map server). Although a fair amount of work has been done in ArcGIS, being able to use it commercially is going to incur a licensing fee we would prefer to avoid. Ultimately we'd like to have a customer facing map viewable online through our website. I'm thinking it would be cool to have an app featuring the map tied into a gps locator, so that customers could automatically be fed contact and product info particular to their location. I've downloaded QGIS and played around with that a bit, making maps is fairly easy and straightforward, but integrating map projects into a website seems to be an involved affair compared to ArcGIS. The other option I'm considering is Google Maps Engine, which seems to be a rather attractive choice - map making is fairly simple and straightforward, website integration should be well supported considering its Google, just waiting to hear back from Google on pricing/trials. If anyone could give me some advice on which program would be right for my situation, it would be greatly appreciated. Thus far, the choices are: 1) Stay in ArcGIS and pay the license fee - which would be a fairly seamless transition from the trial run, but I'm concerned the program is way too comprehensive for our needs. 2) Go to QGIS and avoid fees, but pay more in time spent recreating vector layers, setting up an online browser for customers - and again I'm concerned this is possibly too much program for our needs. 3) Go to Google Maps Engine, assuming they get back to us.


ogr2ogr fails to convert GeoJson file > 200MB

As a followup to my previous question (Converting SimpleGeo Places GeoJSON file to shapefile), I am now struggling with what I perceive might be an ogr2ogr "memory" limitation, but I cannot be sure. I have successfully converted the SimpleGeo dump files into proper geojson files (verified that are directly readable in QGIS). Now I want to convert this data to SHP or other formats (SQLite, PostGIS). I am running the following command with success on smaller files (where "places_fr.gejson" has a file size of ~118MB):               ogr2ogr -f "ESRI Shapefile" places_fr.shp places_fr.geojson      When I attempt to run this same exact command with a geojson file that is 340MB, ogr2ogr crashes in one of two ways:   1. "lh_table_new: calloc failed" ( _first attempt_ )   2. Windows throws an "ogr2ogr.exe has stopped working. Please close the program." ( _subsequent attempts_ ) I am running the ogr command via FWTools in a Windows 7 environment dual core with 8GB of installed RAM (4GB available, 12GB virtual memory available). Is there a memory limit to ogr2ogr converting these datasets? I have also tried converting into SQLite with the same results, so I don't think it is an SHP conversion issue. Thanks.


.NET library to convert KML to WKT?

Anyone know of a free OpenSource .NET library i could use to convert KML to WKT? I know i can use openlayers and there are many JavaScript possibilities, but i need to do this on the back-end.


Generating Custom Map image using Python with small installation footprint?

As a part of reporting requirements there is a need to generate region maps. These highlight various regions (countries or states/provinces) with different colors (or shades of a color) depending on the input data. Usually those are maps of the whole world with the country granularity. For some countries, however, the requirement is to have state/province granularity and highlight regions (states, provinces) rather than the whole country. The output should be similar to the image shown here, but with all countries by default not colored at all (or white) and other countries/states/provinces with a (shade of) color depending on the value of statistics. Given input below, what would be the (c)leanest way to generate such maps? I would prefer a **`Python`** solution, for this **`desktop`** requirement, but it does not have to be.               data = [ # Country, (optional) State, Data Value         ('CD', 'AB', 0.05), # Canada, Alberta at 5% level         ('US', 'CA', 0.15), # highlight US California with color for 15% level         ('MX', None, 0.05), # highlight whole Mexico with color for 5% level         ('NZ', None, 0.05), # highlight whole New Zealand with color for 5% level         ('DE', None, 0.01), # highlight whole Germany with color for 1% level         # ... usually few hundreds of such specifications         ]      Imagine that this is the only requirement there is for the maps, so I would not like to install a lot of dependencies for only this task. I would appreciate list of packages I need as well as the smallest data in order to make this work. * * * _Edit-1:_ In order to really narrow down the scope to the minimum: consider this to be an automated way of performning the following manual task: > Open an empty map image (like in this image) in a `Paint`, and for every > data use "Fill with Color" tool to color areas accordingly to the data > values.


Python3 libraries for GIS analysis

I am trying to do some spatial data analysis in Python 3 now, and I was wondering what the best libraries are for GIS in Python 3. The Fiona library was developed with Python 3.3 but the Python 3+ version of Shapely is still in development. Other libraries have been ported over to Python3 but I am not sure if they have been thoroughly tested and validated compared to their Python2 counterparts. I was hoping that other developers out there could provide some suggestions. The basic functionalities that I need are as follows.   * Mapping Shapefiles or GeoJSON   * Associating polygons with points and lines   * Mapping vector data as well as raster data So nothing too exotic. Any help would be appreciated.


ARCGIS Python script (needs a raw input line to let the user choose the output file name)

I just need a quick fix here.   Having trouble finding a way to allow the user to input the new output file name (the one that is created after the buffer function completes). Any ideas?   I am also not sure if I created the variables correctly. Using the latest (ArcGIS 10.1 SP1) python script to change point layer to a 50 m buffer layer change the output file name to whatever so it can make a new one               import arcpy       arcpy.env.workspace="C:/Users/oxana/Desktop/David"       arcpy.Buffer_analysis("builteditptFT", "c:/users/oxana/desktop/david/PDpyth1","50   meters","FULL","ROUND","LIST")      note that this only buffers from the point layer that was converted from the centroids of the newly created buildings   it is prone to generalizing around each point because the buffer radiates 50 meteres from the center of each building (not the outsides)


"Map 1" vs "Map 0"

When I compose a map to print, and get everything all set; then close it for some reason; then reopen it--the scale bar appears messed up. I then go to the Scale Bar settings, and find that under "Map" there is more than one option. Changing that fixes the scale bar. How can I fix this permanently? I don't even know what this other "map" refers to. I did not intentionally add it to my map. How can I delete the "Map" that is messing things up? THANKS -TVZ (Using QGIS 1.8.0, Windows 7)


Plot specific area of map

Is there a way to zoom in a map on an R spatial polygons object. i.e. if I have polygons for a whole country but I want to zoom in on a specific city and just plot that area how do I do that ? At the moment I am just plotting a chloropleth of the whole country like this:                   colours<-(brewer.pal(9, "YlGn"))         brks<-classIntervals(poly$var1, n=9, style="quantile")         brks<-brks$brks         plot(ED, col=colours[findInterval(poly$var1, brks,all.inside=TRUE)], axes=F, border=FALSE)      I wanna zoom in and plot a rectangular area around a particular city but can't figure out how.


How to create images from symbol in ArcObjects?

I'm trying to create images from symbols so that I can create a custom symbol dialog. Found this sample code on the ESRI forums and it workes great on almost all symbols. I'm having problems with a IMultiLayerFillSymbol that has a IPictureFillSymbol on one of it's layers. It's throwing an exception on symbol.Draw(geometry); Does anyone know if there is any additional steps that has to taken to get it to work with a IMultiLayerFillSymbol?               static public Image ImageFromSymbol(ISymbol symbol, Int32 width, Int32 height)     {         try         {             Bitmap bitmap = new Bitmap(width, height);             Graphics g = Graphics.FromImage(bitmap);             IntPtr hdc = g.GetHdc();                  IGeometry geometry = null;             if (symbol is IMarkerSymbol)             {                 IPoint point = new PointClass();                 point.X = (width / 2);                 point.Y = (height / 2);                 geometry = point;             }             else if (symbol is ILineSymbol)             {                 IPolyline line = new PolylineClass();                 IPoint ptFrom = new PointClass();                 IPoint ptTo = new PointClass();                 ptFrom.X = 3; ptFrom.Y = (height / 2);                 ptTo.X = (width - 3); ptTo.Y = ptFrom.Y;                 line.FromPoint = ptFrom;                 line.ToPoint = ptTo;                 geometry = line;             }             else             {                 IEnvelope bounds = new EnvelopeClass();                 bounds.XMin = 1; bounds.XMax = width - 1;                 bounds.YMin = 1; bounds.YMax = height - 1;                 geometry = bounds;             }                  symbol.SetupDC(hdc.ToInt32(), null);             symbol.Draw(geometry);             symbol.ResetDC();             g.ReleaseHdc(hdc);             g.Dispose();                  return bitmap as Image;         }         catch (Exception ex)         {             Debug.Print(ex.Message);             return null;         }     }      Update, Added the following code to make it work: Updated the ImageFromSymbol metod with:               ITransformation pTransformation = CreateTransFromDC(                            hDC, size.Width, size.Height);     pSymbol.SetupDC(hDC.ToInt32(), pTransformation);      Added:                private const int LOGPIXELSY = 90;          [DllImport("GDI32.dll", EntryPoint = "GetDeviceCaps",                ExactSpelling = true, SetLastError = true)]            private static extern int GetDeviceCaps(IntPtr hdc, int nIndex);               private static ITransformation CreateTransFromDC(IntPtr hDC,                int lWidth, int lHeight) {                // Calculate the parameters for the new transformation,                // based on the dimensions passed to this function.                try {                    IEnvelope pBoundsEnvelope = new EnvelopeClass();                    pBoundsEnvelope.PutCoords(0.0, 0.0, (double)lWidth,                        (double)lHeight);                         tagRECT deviceRect;                    deviceRect.left = 0;                    deviceRect.top = 0;                    deviceRect.right = lWidth;                    deviceRect.bottom = lHeight;                         int dpi = GetDeviceCaps(hDC, LOGPIXELSY);                    if (dpi == 0) {                        throw new Exception(                          "Could not retrieve Resolution from device context.");                    }                         // Create a new display transformation and set its properties                    IDisplayTransformation newTrans =                        new DisplayTransformationClass();                    newTrans.VisibleBounds = pBoundsEnvelope;                    newTrans.Bounds = pBoundsEnvelope;                    newTrans.set_DeviceFrame(ref deviceRect);                    newTrans.Resolution = dpi;                         return newTrans;                }                catch {                    return null;                }            }


How to get polygons from PostGIS table to PHP array best?

Currently I fetch polygons from a database table using `ST_AsText()`. The results are in text format, like the following.               POLYGON((84 104,212 48,321 122,243 179,275 238,131 240,84 104))      I am not quite sure if `ST_AsText()` and string manipulation is the way to go for performant polygon reading. Is there another way to fetch the points? How can I get all points of a polygon as an array in PHP?


Will PostgreSQL 9 work with ArcGIS 10?

I realize it is not listed as a supported Database. I am interested to see if anyone has tried to install and use it.


GeoServer GetLegendGraphic request ServiceException

I have a layer in GeoServer which is stored in a PostGIS database. The geometrytype of the layer is GEOMETRY and now the table only contains LINESTRING. But I what to add points and polygons to the table, this is the reason I choose GEOMETRY as the geometrytype. My layer is working fine in the Layer Preview... But when I do a GetLegendGraphic request I get the following ServiceException: **_java.lang.IllegalArgumentException: class com.vividsolutions.jts.geom.Geometry is not supported by this method class com.vividsolutions.jts.geom.Geometry is not supported by this method_** My SLD is quite simple:               <?xml version="1.0" encoding="ISO-8859-1"?>     <StyledLayerDescriptor version="1.0.0" 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"       xsi:schemaLocation="http://www.opengis.net/sld http://schemas.opengis.net/sld/1.0.0/StyledLayerDescriptor.xsd">       <NamedLayer>         <Name></Name>         <UserStyle>           <Abstract></Abstract>           <FeatureTypeStyle>       <Rule>       <Title>Test</Title>       <ogc:Filter>        <ogc:PropertyIsEqualTo>         <ogc:PropertyName>feature_type</ogc:PropertyName>         <ogc:Literal>Test</ogc:Literal>        </ogc:PropertyIsEqualTo>       </ogc:Filter>               <LineSymbolizer>           <Stroke>          <CssParameter name="stroke">#FFA618</CssParameter>          <CssParameter name="stroke-width">2</CssParameter>          <CssParameter name="stroke-dasharray">5 2</CssParameter>           </Stroke>         </LineSymbolizer>       </Rule>           </FeatureTypeStyle>         </UserStyle>       </NamedLayer>     </StyledLayerDescriptor>      My GetLegendGraphic request is:   http://localhost:8080/geoserver/wms?REQUEST=GetLegendGraphic&VERSION=1.0.0&FORMAT=image/png&WIDTH=16&HEIGHT=16&LAYER=namespace:layer I have tried the request on a GeoServer 2.0.1 and a GeoServer 2.0.2 and in both cases I get ServiceException as above mentioned. What can I do to get the request working?


Whether to use File Geodatabase (*.gdb), Personal Geodatabase (*.mdb) or shapefiles?

I am a bit confused on the topic of geodatabases. Previously, I have used geodatabases as an easy way to keep all files together rather than having dozens of shapefiles all over the place. And it's easier when finally exporting all the relevant stuff to the client. But earlier it got me thinking; what is really the difference between gdb's and mdb's, from a user's perspective? Which one should I stick to? Why should I use a geodatabase (should I?) instead of shapefiles when analysing and running functions on my data? Pros and cons in general? My primary focus is within ArcGIS, I don't usually edit stuff outside of that environment. This question is somewhat related to an earlier one, but it was asked quite a while ago and references mainly ArcMap 9.3. COMMENT: Thank you all for contributing. To sum this up: due to performance and limitations, FGDB's seems to be what I should use provided no editing outside of ArcMap is necessary. Shapefiles are still good because of their acceptance over platforms and applications.


Quick way to make a convex hull

I am using Python and ArcGIS 10.1. Is there a function or very quick way to extract points out of point feature class into an array **without** looping through a search cursor (as I have millions of points to process). I eventually want to transform this into a multipoint feature and make a convex hull from it. * * * **UPDATE:** Although, initially I thought the 11 million point convex hull optimisation was good, I think it is really too long, in order to create a random convex hull and then do an intersect and switch selection, it took 4 hours and 30 mins. That's NOT INCLUDING the actual processing of the sub-selection of points to (some 100K) create the final convex hull result. So I think it may be back to the drawing board. The real success here for actually completing long operations and memory intensive ones, has been to upgrade ArcGIS Python geoprocessing to 64bit. Before I did this, it would just hang, now it has all the memory it needs and can really slog through some long iterative processes.


How to create a layer switcher in leaflet

I had used the following layers in leaflet,now i want to switch between the 4 layers ,please tell the code to create a layer switcher control in leaflet               // night layer     //var osmUrl = 'http://{s}.tile.cloudmade.com/a7db5562d0cd4fe0a2d8983fbb3b5519/999/256/{z}/{x}/{y}.png',               //golden layer     var osmUrl = 'http://{s}.tile.cloudmade.com/a7db5562d0cd4fe0a2d8983fbb3b5519/997/256/{z}/{x}/{y}.png'          //black layer     //var osmUrl = 'http://{s}.tile.stamen.com/toner/{z}/{x}/{y}.png',          osmAttribution = 'Map data © 2012 OpenStreetMap contributors',     osm = new L.TileLayer(osmUrl, {maxZoom: 18, attribution: osmAttribution});


What is the most appropriate technology for publishing Census (and other) data for small areas?

We have a web application which shows the (United Kingdom) Census data on thematic maps using client-side code (based on OpenLayers) and KML. // This section can be skipped // begin At the moment it works in the following way   * We store area shapes as KML in database (just shapes, no styles, colours or other data)   * We also store the Census data against each area (e.g. area="Output area E01000123", time period=2011, variable= "% of single people", metric value = 34.12) in the same DB. Of course, data structure is normalized.   * A user selects a region (e.g. London) and a variable   * We get all the small areas within the selected region, generate the KML and apply colors based on metric values.   * The web page displays KML on a map (using OpenLayers-based Javascript code). //end This approach works well when the number of small areas in a region is not too large. In some cases the number of small areas might reach 5000. In this case client-side code is too slow in IE8 (which we have to support). So I guess we have to generate an image from the shapes server-side. What is the most appropriate technology? In the past we used GeoServer with PostGis backend. It was not too fast and quite unstable. I can explain how we used it. UPDATE I decided to use GeoServer + PostGis. GeoServer serves images. When a user selects an area, the client-side code creates a custom SLD and passes it to GeoServer. When no or few areas are selected, it is very fast (~1sec), even with 30000 polygons in the layer. If a user selects e.g. 2000 areas, it takes several seconds (I guess because GeoServer has to process a huge SLD file).


Estimation of values for unsampled points

I want to measure biomass on a coastal mudflat. I can only access points inside the polygon. Are there any methods available that would allow me to estimate values of the points outside of the polygon, based on the values of points inside the polygon?               set.seed(5)     x <- rnorm(50, -1.841, 0.01)     set.seed(50)     y <- rnorm(50, 55.663, 0.01)     xy <- data.frame(x,y, values=rnorm(50))     coordinates(xy) <- c("x", "y")     proj4string(xy) <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84")     plot(xy)          makePolygons <- function(coordsx, coordsy){            coords <- matrix(c(c(coordsx, coordsy)), ncol=2)       p <- Polygon(coords)       p <- Polygons(list(p), ID = "p")       myPoly <- SpatialPolygons(list(p))       spdf = SpatialPolygonsDataFrame(myPoly, data.frame(variable1 = c(2),                                                          variable2 = c(3), row.names = c("p")))       proj4string(spdf) <- CRS("+proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0")            print("polygon is in longlat!!!")            spdf          }          myPoly <- makePolygons(coordsx=c(-1.841960, -1.843464, -1.888623, -1.841960),                            coordsy=c(55.633696, 55.68178, 55.63841, 55.633696))          plot(myPoly, add=T)      ![enter image description here](http://i.stack.imgur.com/yF44s.jpg)


WFS - From PostGIS to Shapefile with GeoServer

From Server side, I'm using GeoServer and from client side I'm using OpenLayer to allow the users see the layers(PostGIS) on the map. The new requirement is that the users may download layers in shapefile format. What is the best way to do this? PS: I would prefer add the option to download them in the same map view.


Mapserver - What are some options for a GUI for mapfile creation?

What are some possible options for a GUI to manage style and mapfile for Mapserver. I am interested in all solution, QGIS plugin, ArcGIS mxd or MapInfo workspace converter or else. Thank you for all your suggestions


DataB and DataB-copy worked on simultaniously - how merge DataB-copy into DataB

ArcGis 10 - EDITED (sorry for not being more specific, quite new to the arcgis world) We need to enter A LOT of Map data (Sewage system) into a Database in a relatively short time. It involves an area of about 14km^2 with pipes, manholes (couple hundred) and relvant buildings. All of which have numerous attributes. These are cataloged in 30 some feature classes. Data entered via scans that are georeferenced and then copied by hand via ArcMap. I assume that they're ESRI-feature classes [ArcCatalog says "Personal-Geodatabase Feature- Class"]. The database (.mdb file) is about 130mb. Following idea: Database "A" is main Database. A copy is made "A-copy". With two liscenses two people can work on A simultaneously (different areas). After a certain time, we want to import all the entered data from "A-copy" back into "A". How would one best attempt to do this this without overwriting the changes made in "A" during this time? (the entered data is as good as idependant from the data already entered; It would be possible to wipe all tables from "A-copy" before starting) Both users would ideally be accesing and entering data into the same tables. [Maybe it would be possible to limit that to a certain extent, not quite sure how topo-lines owuld react if e.g. one user placed all the manholes in one area and then on the next day a second person would connect them in the copy- db and then transfer the pipe feature-class back into the original.] From what I can tell sde is not in use; we have no possibility of multi-user simultaneous access to this particular database. Is there a way to qickly check this? I didn't set up the system and am not versed in doing that either. This wouldn't have to happen on a daily basis, so a multistep process would be still be feasable. Thanks a million


Why the layer can be loaded to Google Map not OpenStreet Map using MapServer?

Updates: Problem Solved !!!!!!!!!!!!!!!!!!!!! Replace "wms_srs" "EPSG:3857" to "wms_srs" "EPSG:900913" in Mapfile !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! I am using OpenLayer WMS to loaded a layer from the PostGIS Raster Table, the JavaScript on openlayer are as following:                       function init(){             var map = new OpenLayers.Map('map');                  var imagery = new OpenLayers.Layer.WMS(                 "Global Imagery",                 "http://maps.opengeo.org/geowebcache/service/wms",                 {layers: "bluemarble",                 }             );              // var imagery = new OpenLayers.Layer.OSM( "Simple OSM Map");     /*It doesn't work with this OSM basemap, showing error:           WMS Server Error: Invalid SRS is given */             map.addLayer(imagery);                  localWMS = new OpenLayers.Layer.WMS( "Local WMS layers",                       "http://127.0.0.1/cgi-bin/mapserv.exe?MAP=C:/Users/gyu1/Desktop/Newfolder/mapserver.map",                       {layers: 'polygon',transparent: true}, {isBaseLayer: false} );             map.addLayer(localWMS);                  map.setCenter(                 new OpenLayers.LonLat(1.147, 52.472).transform(                     new OpenLayers.Projection("EPSG:4326"),                     map.getProjectionObject()                 ), 5             );                  map.addControl( new OpenLayers.Control.LayerSwitcher());         }      The map file linking the Data and styling are as following:               MAP     NAME "polygon"     WEB              METADATA                "wms_title"          "MapServer"                "wms_onlineresource" "http://127.0.0.1/cgi-bin/mapserv.exe?MAP=C:/Users/gyu1/Desktop/Newfolder/mapserver.map"                "wms_srs"            "EPSG:4326 EPSG:3857"             END           END          PROJECTION         "init=epsg:3857"     END                LAYER             NAME "postgis"             TYPE raster             STATUS ON             DATA "PG:host=localhost port=5432 dbname='fusiongeo' user='postgres' password='**' schema='public' table='london' mode='2'"              PROCESSING "NODATA=-9999"             PROCESSING "SCALE=AUTO"                  CLASS                 NAME "boring"                 EXPRESSION ([pixel] < 50)                 COLOR 250 100 0             END             CLASS                 NAME "mildly interesting"                 EXPRESSION ([pixel] > 50 AND [pixel] < 100)                 COLOR 100 0 0             END             CLASS                 NAME "very interesting"                 EXPRESSION ([pixel] >= 100)                 COLOR 0 50 0             END         END          END # end map      It works for the Google Base Map not the OpenStreetMap. I am quite new with the MapServer, openlayer and I don't know how to change the SRS to make it work. Thanks. The data in PostGIS are in EPSG:3857.


How to select landlocked countries in Africa using Select By Location?

I wish to select the landlocked countries in Africa using just Select by Location. I have already dissolved the boarders to get one big polygon and I thought I might combine these two telling ArcGIS somehow that I just want the countries to be selected that are not touching the outside borders, but how can I do the next step?


Field calculatur in ArcMap 10: How to specify field row?

I'm new to arcmap, so I'm not entirely sure about terminology which is probably why I can't seem to find the answer to this seemingly trivial question. In the attribute table I want to be able to create a new field which takes the value from another field and adds the value from the previous row to it. My idea was simply to make a some kind of a `for` loop but to do that it seems I need to be able to specify a row number of a field, and I don't know how to do that. Any help would be appreciated, ranging from a better solution to the answer to this specific question :).


Which interpolation method is the most suitable for DEMs?

Which interpolation method is the best suitable for hill/valley areas and normal ground with undulations and rivers? When working with Aster 30mX30m DEMs, which method will give the most realistic elevations? I am trying with different methods and getting different contour distribution and shapes but don't know which will give the most accurate and realistic values.


Is there software that will do reverse geocoding to convert photographs' GPS locations to addresses?

I am using a standalone GPS tracker that produces GPX files and then clock- sync photos with it using Microsoft Pro Photo Tools. This works pretty well. But I was also using a feature of Pro Photo Tools called "Get location text" which uses the latitude/longitude information from EXIF and translates it using Windows Live services to actuall addresses. At least it used to, but stopped recently, probalby because of some change of the online API or whatever, and a lack of update for the software. So, the questions is: do you know of any other software that can bulk translate location info into addresses and write that back into the photo's metadata? I expect the solution to use Google Maps or Bing Maps APIs - both provide (reverse) geocoding world-wide and publicly, and work quite well through their websites. I find it hard to believe that no one has created an app yet that would do this in bulk for photos.


Color a layer of the map based on the location of dots on another layer

Thanks for helping me in GIS 101 Problem #2, now I have followup questions: I created a map as below: ![enter image description here](http://i.stack.imgur.com/712oY.png) The map is created using the following script:               ct <- readShapeSpatial("housect_37800_0000_2010_s100_census_1_shp/wgs84/housect_37800_0000_2010_s100_census_1_shp_wgs84.shp")          ct_mod <- fortify(ct,region="SLDLST10")     chart <- ggplot(data=ct_mod,aes(long,lat))     chart <- chart + scale_x_continuous(limits=c(-73.8,-71.7),breaks=seq(-74,-71,0.1))     chart <- chart + scale_y_continuous(limits=c(40.9,42.1),breaks=seq(40,43,0.1))     chart <- chart + geom_polygon(fill="grey80",aes(group=group))     chart <- chart + geom_path(color="white",aes(group=group))       # coord_equal make x-y same scale     chart <- chart + coord_equal()          chart1 <- chart + geom_point(data=hosp.list,aes(x=coord_x,y=coord_y))        This is map of the Connecticut Census State Legislative District Map with several hospital located, you can see a few regions dot and a few does not. I want to fill the regions with dot using one colour, and regions without dot using another colour. Can anyone give me a clue? Thanks.


Using openlayers to render UK National Grid maps

Is there a way I can use Openlayers with a 1 pixel = 1m/2m/4m OS Grid map at high resolutions (i.e. OS Street View), but switching to a Mercator projection when you zoom out (i.e. for Google Maps or OSM), but have markers and/or vector overlays reprojected with the change of zoom? As OS grid and OS Street view is based on powers of powers of 10, I've split it into 250 pixel tiles, rather than the 256 pixels used for most web maps. I have a folder structure like: `xy/xy/xy/xy-z.png` where each folder is a power of 10 grid square, i.e. each xy represents hundreds of km, tens of km, km, and position within the 1km square respectively. There are 4 zoom levels in the bottom folder, 1 being the whole 1km square, 2 is 4 500m squares and 4 is 16 250m squares.


Could not install Applications Manager (ArcGIS Server 10)

I try to install ArcGIS Server version 10. Everything worked fine with installing the GIS Services. The Server is running. I can create services and see them in ArcCatalog. But when I try to install the Web Applications I couldn't select the "Applications Manager". I try to install it on Windows Server 2008 (64bit) with IIS7. I double checked the settings and required extensions for IIS and it seems to be correct. I even changed the IIS seting from 64bit to 32bit. But the "Application Manager" in the Select Feature Dialog is still not shown. I appreciate any help with this. Kai


How do I resolve line outlines jutting into other lines in Mapnik/Cascadenik

We are designing a basemap, and I am trying to resolve a problem with outlines that comes up in roads. The idea is to have lighter colors roads, but they don't get pushed out far enough with our light background. I've tried various combinations of outline-cap and outline-join, but still run into tiny outlines sticking into the joints of lines/roads. is there something I am missing to resolve this, or should i just choose a color that stands out above the land?


A machine to measure latitude to 6 decimal places of a degree

I am doing a PhD research project to do with medieval and ancient archaeoastronomy which requires me to be able to measure latitudes up to an accuracy of 6 decimal places of a degree of latitude. Could you point me in a direction of where I might be able to hire or buy a machine that does this? What is the general price range for something like this? Thank you.


Trigger on link length and filling attribute table

I create linear geometry in QGIS 1.8.0 (geometry stored in PostGIS 2.0). Help create trigger functions:   1. in database should not be a lines less than 2 meters (linear geometry with length <2 meters should be corrected in the process of creating the geometry) - if I created linear feature with length <2 meters I receive message "feature "gid" have wrong length"    2. auto fill Length column in attribute table, when new linear objects added For calculating length I use `ST_Length_Spheroid(r.the_geom,'SPHEROID["WGS 84",6378137,298.257223563]')` function. So, example trigger function for intersecting polygons below:               CREATE OR REPLACE FUNCTION check_geom_overlap() RETURNS trigger AS $$     DECLARE     gid integer;     BEGIN     SELECT INTO gid c.gid     FROM town c WHERE c.gid <> NEW.gid AND ST_Overlaps(c.the_geom, NEW.the_geom);          IF gid IS NOT NULL THEN     RAISE EXCEPTION 'Modified feature % overlaps with feature %!', NEW.gid, gid;     END IF;     RETURN NEW;     END; $$ LANGUAGE 'plpgsql';          CREATE TRIGGER check_geom_overlap_trigger BEFORE INSERT OR UPDATE     ON town FOR EACH ROW EXECUTE PROCEDURE check_geom_overlap();      Thanks!


Publish a 'flattened' WMS in ArcGIS Server

we have a WMS which contains several layers, published using ArcGIS Server/WIN 10.2. I am attempting to determine how we could 'flatten' the service to only contain 1 or 2 layers instead of several - which makes it easier for implementation into openlayers, etc - but still maintain all data. identify capability does not need to be available, nor does the option to turn on and off layers. Example:   ![Current WMS config](http://i.stack.imgur.com/T5pJ0.png)   * 'flatten' the WMS to look like: ![new config - one layer only](http://i.stack.imgur.com/HKbq0.png) -or- ![new config, multiple layers](http://i.stack.imgur.com/4LVDZ.png) Thanks Mucho


Why does GPS altitude jump sometimes?

I am currently developing a little application in Android which uses GPS. When i click a button i start to record a "track". If i walk around i save my position every second into a file. This works perfectly except for one thing: I do not only save longitude and latitude but also the current altitude. Looking through the data I recognized that the altitude jumps some times. For example the altitude values look like this: 271.4073616670083; 271.5165709069828; **345.6;** 271.58538936056465; The values around 270 metres are quite good. Jump like these happen irregularily but quite often (every 15 seconds there ist at least 1 jump). And the values always jump up to something around 340 metres. While recording the track with my application I have also been walking around with an app called NMEA Recorder. This app recorded the altitude values as well but they don't jump like the values I receive. I know that the altitude values of GPS (calculated by the WGS84-Ellipsoid) aren't very accurate but that still doesn't explain why the values jumpt sometimes. Would be nice if somebody could tell me why jumps like these happen.


Display a color ramp legend (of a vector layer) in the qgis map composer

I would like to have a color ramp instead of 10 categories of color when displaying my vector layer in the map composer of QGIS. Is it possible ? Thank you in advance


How to use ESRI Javascript Build Roll-Ups?

On their developer guide website (LINK) ESRI mentions you can use build roll- ups to more efficiently load modules, and since my define array is starting to look like the opening text in Star Wars I thought I'd try it out, but so far I haven't made it work. The guide says including "esri/dijit/editing/Editor-all" will load a bunch of other modules (like TemplatePicker, FeatureLayer, edit,...) but I'm not sure I'm doing this correctly because after adding Editor-all my script that creates a new Editor object fails. So after adding the Editor build roll-up this fails:               editorWidget = new Editor(params);      But using the full api path still works:               editorWidget = new esri.dijit.editing.Editor(params);      So what's the point of using roll-ups if I have to go back to a pre AMD way of referencing the modules? I can't find any script examples in AMD style that use build roll-ups on ESRI's dev site, so can someone please suggest what I may be missing or not doing correctly? Thanks in advance!


Is it possible to store a default symbology layer in an arcpy add-in?

Since arcpy is limited in how symbology can be manipulated, I was curious to know if a symbology layer can be stored in an arcpy add-in. Does anyone know if it is possible to store a .lyr file, or anything else for that matter, in a .esriaddin so that it can be used once the addin is installed?


Access to WMS from ArcMap with extra parameters

I'm trying to access AIRS WMS from ArcMap with additional time parameter (&TIME=2012-07-01) to see this summer's CO concentration. But the server returns a response with standard parameters (last month). How to pass this parameters in the ArcMap WMS-connection correctly?


ArcGIS 10.0 Running Python function CADToGeodatabase() crashes ArcMap/ArcCatalog

My code:                   import arcpy, csv         from arcpy import env              dgnfile = str(arcpy.GetParameterAsText(0))         outFeatureClass = str(arcpy.GetParameterAsText(1))  #the fc in my fgdb template to be  populated by data read from the dgn         i = 0         while outFeatureClass[i:i+4] <> ".gdb": i+=1           inFeatureClass = outFeatureClass[:i+4]           CADdataset = "GD_Temp"           refscale = "1000"           arcpy.CADToGeodatabase_conversion(dgnfile,inFeatureClass,CADdataset,refscale)      Apologies to @PolyGeo, @Arabella - This is my code which crashes out. I am able to read geometry out of a CAD dataset which has been created in ArcCatalog manually, and populate an existing FC in my FGDB. But the dgn-to- CADdataset part crashed. When I placed my first question I ASSUMED the CADToGeodatabase routine would work. Is it a simple Python code error? ESRI are their usual helpful self. By the way, I've omitted the code which subsequently reads the output from CADToGeodatabase(). * * * @PolyGeo Ok, following is what I am feeding into CADToGeodatabase routine, as quoted from PythonWin debug session (...=intevening folders): argument 1: 'C:\GeneralWorkData\\...\9202.dgn'   argument 2: 'C:\GeneralWorkData\\...\GD_9202_ED5031N_FGDB_Pop.gdb'   argument 3: 'GD_Temp'   argument 4: '1000' It is supposed to create a new feature dataset in the .gdb, which is what I am trying to achieve here. The Sample Code you mentioned is the method I am using for my work, but it crashes ArcMap/ArcCatalog/PythonWin Debug, with no (helpful) error message. These are the right argument formats according to Resource Centre Help for Arc 10.0! Don't want to look at "Interoperability"(?) right now, as I have code which works to read out of the CADToGeodatabase results so it would mean starting all over again. In other words, I can import dgn data by hand in Catalog then read that data in Python and put it where I like, so the second half of the Python project is done, now I need the first. I posted originally prematurely, didn't know this would crash, apologies again. Is there any chance someone could try this routine with real dgn data, cos I am thinking we are looking at an ESRI bug...? Any help appreciated!


Importing grib, grib2, netcdf to PostGIS

According to grib2 to raster? (grib2 -> raster2pgsql -> PostGIS 2.0 -> GeoServer), it seems like that grib2 can be treated as raster directly. Is it the so called "standard" way? How about grib and netcdf?


Shoreline identification using ERDAS 2011

How can we separate shoreline from shallow water sea and land area in a LISS 3 image using ERDAS 2011 software? thanks.


Why does the added raster layer have strange edges?

I'm using Geoserver in conjunction with the google maps API to add a raster layer onto a map. It loads fine, however there's some strange problem with opacity at the edge of the layer where you can see a boundary. Have you any idea how to fix this? ![Screen shot](http://i.stack.imgur.com/qg948.png) (The original is at obstest.heliohost.org/map2.html.) I'm totally new to gis and to Geoserver so sorry if this is an obvious question! Thanks.


"Schema does not exist" error trying to load raster using raster2pgsql

I receive a "schema does not exist" error when trying to load a raster to a PostGIS database (see full error at the bottom). The schema does exist in the database. This is the syntax I'm using:               raster2pgsql -s 2048 -I -C -Y C:\Users\My_Book\Desktop\LOAD2POSTGIS\Aspect.tif environmental.aspect > Aspect.sql     psql -h localhost -U postgres -d Howbill -f Aspect.sql      I am able to load the raster into the default public schema if I leave out environmental.aspect. However, that leads to other problems when moving the raster table from the public schema to another schema. More details regarding the moving from public schema to other schemas can be viewed here: Can't view PostGIS rasters after changing schemas               **Full error:**          C:\Program Files (x86)\PostgreSQL\9.2\bin>raster2pgsql -s 2048 -I -C -Y C:\Users\My_Book\Desktop\LOAD2POSTGIS\Aspect.tif environmental.aspect > Aspect.sql     psql -h localhost -U postgres -d Howbill -f Aspect.sql     Processing 1/1: C:\Users\My_Book\LOAD2POSTGIS\Aspect.tif          C:\Program Files (x86)\PostgreSQL\9.2\bin>psql -h localhost -U postgres -d Howbi     ll -f Aspect.sql     BEGIN     psql:Aspect.sql:2: ERROR:  schema "environmental" does not exist     psql:Aspect.sql:3: ERROR:  current transaction is aborted, commands ignored unti     l end of transaction block     psql:Aspect.sql:5: invalid command \.     psql:Aspect.sql:6: ERROR:  syntax error at or near "010000010013"     LINE 1: 010000010013B444D4496B144013B444D4496B14C082F38736BB1CD94011...             ^     psql:Aspect.sql:7: ERROR:  current transaction is aborted, commands ignored unti     l end of transaction block     ROLLBACK          C:\Program Files (x86)\PostgreSQL\9.2\bin>


How to stop refresh Google maps?

I am using Subgurim.Controles google maps in my project,i was updating every 40 with new marker to the Gmap. I just want to update position of marker every 40 seconds, not refresh entire Gmap. Here i am using timer for every 40 sec to get new latlng from database table; aspx code:               <asp:Timer ID="Timer1" runat="server" OnTick="Timer1_Tick">                  </asp:Timer>                  <div id="map">                  cc1:GMap ID="GMap1" runat="server" Height="100%" Width="100%" />                  </div>       C# code:                 protected void Timer1_Tick(object sender, EventArgs e)       {  GMap1.resetMarkers(); Subgurim.Controles.GLatLng gLatLng = new         > Subgurim.Controles.GLatLng(Convert.ToDouble(status.Rows[0][0]),         > Convert.ToDouble(status.Rows[0][1])); GMap1.setCenter(gLatLng, 16,         > Subgurim.Controles.GMapType.GTypes.Normal); Subgurim.Controles.GMarker         > oMarker = new Subgurim.Controles.GMarker(gLatLng);         > GMap1.addGMarker(oMarker); }


Automatic selection of legend object in QGIS

I'am using the Corine Land Cover data for a projet. There is around 40 types of land use in this database. In my print composer with my scale (1:20 000) I neve have more than 10 type of land use. I would like to show in the legend, only those 10 type of land use. Is that possible? (This option is available for the table) I'am using QGIS master 1.9. Thank you for your help. Cheers Arthur


Progress bar for Leaflet map

Is there any way to show progress bar while waiting that all layers are rendered on Leaflet map? In this specific case, while geojson layers are rendered. Thanks.


AddJoin failing in Python ArcGIS

My AddJoin operation is failing:               arcpy.CopyRows_management(lup_file,out_dir+'lup')     try:         # Perform copy rows on merged csv file         print 'start'         aa = arcpy.Raster(SEIMF)         arcpy.AddJoin_management(aa, "VALUE", out_dir+'lup', "VALUE", "KEEP_ALL")         print 'end'         #arcpy.AddJoin_management(gis_dir+SIF, "VALUE", out_dir+'lup', "VALUE", "KEEP_ALL")     except:         print arcpy.GetMessages()      error messages: > Executing: AddJoin C:\Users\Documents\Projects\AB\Inputs\gis\sif VALUE > C:\Users\Documents\Projects\AB\Outputs\C\output\lup VALUE KEEP_ALL   >  Start Time: Sat Feb 22 02:27:50 2014   >  Failed to execute. Parameters are not valid.   >  ERROR 000840: The value is not a Table View.   >  ERROR 000825: The value is not a layer or table view   >  ERROR 000840: The value is not a Raster Catalog Layer.   >  ERROR 000840: The value is not a Mosaic Layer.   >  WARNING 000970: The join field VALUE in the join table sif is not indexed. > To improve performance, we recommend that an index be created for the join > field in the join table.   >  Failed to execute (AddJoin). Does anyone know what is happening? I get this error only in the python code and can execute AddJoin just fine if I use the ArcGIS manually.


Match Lines with Data Driven Pages

I have a map book created via the Data Driven Pages tool in ArcGIS. Specifically, I used the Strip Map option when I created the book. And I have employed the 'Angle' attribute of the index layer in order to rotate each page so that the route is (generally) shown west-to-east, with the north arrow turning dynamically with each page. Now, I understand the labeling options for the Data Driven pages. But I don't want to merely label the pages - I want to create vertical match lines that appear on each page. Can this be done, even as each record of the index is turned 'upright'?


Point Within, Query in a buffer polygon in WMS like in ArcGIS ESRI

I need to do a query to filter the result inside a polygon or geometry, like this: http://help.arcgis.com/en/webapi/javascript/arcgis/jssamples/query_buffer.html But, i need to do this in a opengeo WMS server. I have the following query: suite.opengeo.org/geoserver/usa/wms?SERVICE=WMS&VERSION=1.1.1&REQUEST=GetFeatureInfo&LAYERS=usa%3Astates&QUERY_LAYERS=usa%3Astates&STYLES=&BBOX=-79.89257812499999%2C41.68111756290652%2C-64.88525390624999%2C46.41513877649202&FEATURE_COUNT=5&HEIGHT=600&WIDTH=1366&FORMAT=image%2Fpng&INFO_FORMAT=text%2Fhtml&SRS=EPSG%3A4326&X=681&Y=298 And i need to filter the results inside of a polygon. The parsed query:                 SERVICE: "WMS",       VERSION: "1.1.1",       REQUEST: "GetFeatureInfo",       LAYERS: "usa:states",       QUERY_LAYERS: "usa:states",       STYLES: "",       BBOX: BBOX,       FEATURE_COUNT: "5",       HEIGHT: HEIGHT,       WIDTH: WIDTH,       FORMAT: "image/png",       INFO_FORMAT: "text/html",       SRS: "EPSG:4326",       X: X,       Y: Y      I find this XML code, but i dont know if is this way:               FILTER: '<Filter xmlns:gml="http://www.opengis.net/gml"><Intersects><PropertyName>the_geom</PropertyName><gml:Polygon srsName="http://www.opengis.net/gml/srs/epsg.xml#2400"><gml:outerBoundaryIs><gml:LinearRing><gml:coordinates>1 1</gml:coordinates></gml:LinearRing></gml:outerBoundaryIs></gml:Polygon></Intersects></Filter>'      **EDIT:** In ESRI Wms, we can use this:               **point_within**     Description: An SLD style that uses a "within" spatial filter to render only features within a polygon geometry     URL of SLD file: http://sampleserver1.arcgisonline.com/arcgis/wms/slds/point_within.xml     GetMap request: http://sampleserver1.arcgisonline.com/ArcGIS/services/Specialty/ESRI_StatesCitiesRivers_USA/MapServer/WMSServer?VERSION=1.3.0&REQUEST=GetMap&CRS=CRS:84&BBOX=-178.217598,18.924782,-66.969271,71.406235&WIDTH=765&HEIGHT=360&LAYERS=0,1,2&STYLES=,,within&EXCEPTIONS=application/vnd.ogc.se_xml&FORMAT=image/png&BGCOLOR=0xFFFFFF&TRANSPARENT=TRUE&SLD=http://sampleserver1.arcgisonline.com/arcgis/wms/slds/point_within.xml      But, how I can do the same with others WMS? By example the http://suite.opengeo.org/geoserver/usa/wms ? I tried using a FILTER with GML, but I dont know how to do this. Someone help me?


Create Text files from a loop with Python?

I need to create some files to hyperlink to from Arcmap. I have a folder of jpeg files (100.jpg, 101.jpg etc...). I need to loop through them and create a text file of the same name except replace the *.txt at the end with *.fsv (100.fsv, 101.fsv etc...). This text file of the same name (with the txt replaced with fsv) is a set of instructions on how to open that jpeg in a piece of software so I need to write some instruction in it for example 100.fsv would contain: open=100.jpg set to 100 everything in the txt file is static except for the name of the jpg it opens. Thanks for any help.


How to make the Overlays in OpenLayers exclusives using events at checkbox level?

On my map I have 3 overlay layers, all of which need to be mutually exclusives, ie. can't be on (checked) at the same time. I thought of using events to trigger the overlays as the following in my index.html main file:               function init(){       var map, kml_1Layer, kml_2Layer, kml_3Layer;      ...     kml_1Layer.events.register("visibilitychanged", this, layerChanged);     kml_2Layer.events.register("visibilitychanged", this, layerChanged);     kml_3Layer.events.register("visibilitychanged", this, layerChanged);      ...     }          function layerChanged(layer) {          if (layer.object.visibility) {          switch (layer.object.name) {             case "kml_1Layer":                   kml_2Layer.setVisibility(false);                   kml_3Layer.setVisibility(false);                   break;             case "kml_2Layer":                   kml_1Layer.setVisibility(false);                   kml_3Layer.setVisibility(false);                   break;             case "kml_3Layer":                   kml_1Layer.setVisibility(false);                   kml_2Layer.setVisibility(false);                   break;          }       }     }      You can see the web site here: http://genepoulin.net/geo/index.html Obviously the events are not triggered. What am i missing here? Thanks. UPDATED and Resolved... Corrected code (i was wrong naming the switch cases):               function init(){       var map, kml_1Layer, kml_2Layer, kml_3Layer;      ...     kml_1Layer.events.register("visibilitychanged", kml_1Layer, layerChanged);     kml_2Layer.events.register("visibilitychanged", kml_2Layer, layerChanged);     kml_3Layer.events.register("visibilitychanged", kml_3Layer, layerChanged);      ...     }          function layerChanged(layer) {          if (layer.object.visibility) {          switch (layer.object.name) {             case "*given label of layer1*":                   kml_2Layer.setVisibility(false);                   kml_3Layer.setVisibility(false);                   break;             case "*given label of layer2*":                   kml_1Layer.setVisibility(false);                   kml_3Layer.setVisibility(false);                   break;             case "*given label of layer3*":                   kml_1Layer.setVisibility(false);                   kml_2Layer.setVisibility(false);                   break;          }       }     }       You can see the result on my page. Thanks all.


Why do I need to restart tomcat when using WFS from GeoServer across different apps and workspaces?

I've built a couple of apps that consume WFS 1.1.0 from GeoServer 2.2.2. I organise my layers in different workspaces, and they have a range of sources including shapefile, oracle, postgres etc. I've noticed that wfs services only work in one application at a time. If i bounce tomcat then whichever app I go to first will be able to consume wfs. If i hit a second app, no luck. I can see the request for the second app being made in the logs:               14 May 16:06:24 INFO [geoserver.wfs] -     Request: getFeature     service = WFS     version = 1.1.0     baseUrl = http://localhost:8080/geoserver/     query[0]:         srsName = EPSG:3857         typeName[0] = {http://www.openplans.org/usaid}jamoats_geom     outputFormat = text/xml; subtype=gml/3.1.1     resultType = results      Puzzled. Am i missing some config or security setting? Thanks! <\---Update after iant's request for more info--> I inherited this stack from another developer... I'm running nginx/0.7.65 & Tomcat 6 on Ubuntu 10.04.4 LTS. geoserver.war is deployed at /var/lib/tomcat6/webapps. I build a simple web-app (js/css/html) and package as a .war. I psftp this over to my server next to geoserver. I add sites to nginx for each web-app by editing /etc/nginx/sites-available/default. I interrogated the XHR responses for a successful and an unsuccessful WFS request. I'm using two projects here to debug; 1.climvar and 2.usaid. If I load up climvar after restarting Tomcat my wfs features display. If I load up climvar after restarting Tomcat, BUT visit usaid first, no wfs features display. I looked at the responses in more detail. In both responses 82 features came back: In the successful response:               xmlns:climvar=”http://www.fasga.com/climvar"...      In the unsuccessful response:               xmlns:usaid=”http://www.fasga.com/usaid”.       So it seems the wrong namespace is being 'remembered' across my WFS. There is no difference between the request payloads. In the unsuccessful response, all the prefixes are null. Instead of:               <climvar:id>0</climvar:id>       I get:               <null:id>0</null:id>      Here's my js layer definition for climvar:                   var layerOpg = new OpenLayers.Layer.Vector(nsClimVar.Constants.OPG_LAYER_NAME, {             //id : "layerOpg",             styleMap : nsClimVar.Styles.StyleMaps.styleMapOpg,             strategies : [new OpenLayers.Strategy.Fixed()],             protocol : new OpenLayers.Protocol.WFS({                 version : "1.1.0",                 url : "http://" + nsClimVar.Constants.SERVER_IP + "/geoserver/climvar/wfs", // geoserver wfs endpoint                 featurePrefix : "climvar", // geoserver worspace name                 featureType : "opg_new", // geoserver Layer Name                 featureNS : "http://www.fasga.com/climvar",                 srsName : "EPSG:900913" // default             })         });      Layers in usaid are declared in the same way. In GeoServer, both workspaces have the namespaces set correctly, both are enabled, both have WFS enabled. I've tried removing the featureNS and featurePrefix params but this makes no difference. The only workaround I can think of is having all my layers in one workspace which I don't want to do. Any ideas?


GeoJson attributes encoding limitation

I have problem when trying to encode attribute values in latin1 using ogr GeoJson driver. Accent char are not displayed correctly. Is the geoJson encoding is limited to UTF-8? I cant find any info about this in specs. Can someone try to encode a geojson in latin1 to see if it works?


Finding Code in QGIS Select Features Toolset

Can anyone tell me in what directory the pre-installed Python scripts for the QGIS selection tools (i.e., Select by Point, Select by Polygon) are found? I am using on a Windows 7 installation of QGIS 2.0.1.


Does converting polygon to raster generate more accurate results from the Zonal Statistics in Arcgis 10.2

I have heard the following for Arcgis 9.1 and wonder whether it is still an issue in Arcgis 10.2. If so, should I convert the polygon into raster before dissolve or convert the dissolved polygon directly? http://support.esri.com/en/knowledgebase/techarticles/detail/19184


QGIS Routing between two Points without Roads

I'm using QGIS 1.8 and i need to analyze possible routes for roman roads. Is there a Possibility to show routes which comply special criteria, e.g. straightest way between two points, or path which is not higher than 50m... Since i search for possible **roman** roads, there are no streets with which QGIS can work. My idea is, it should analyze the map by my own criteria to find the routes. Is there any Plugin or function to do that?


How to convert lines from dxf to polygons?

I have several dxf files. Each file, contains the 0 layer. The 0 layer contains lines (which make a closed polygon). In each dxf the lines have their color and line width. I do not have autocad or else like. I have Qgis 1.7.4. I can put these dxf files into vector layers (as line vector) into Qgis. I do not know the way (using Qgis 1.7.4 ): 1\. To turn each vector layer (with lines) to vector layer with polygon. To “copy” some of these polygons (vector layers) into a new vector layer. 2\. Can each polygon into the new vector layer to have different border color? thank you


osm2pgrouting stuck (showing classes table created) when importing an osm file

I am trying to import an osm file using osm2pgrouting. while following the pgrouting workshop. (file is 700mb~) im running ubuntu 12.04 on m1.large on amazon. the same file was successfully upload to postgres with osm2pgsql. i am using the following command:               osm2pgrouting -file "/home/ubuntu/israelmap/israel_and_palestine.osm"  -conf "/usr/share/osm2pgrouting/mapconfig.xml" -dbname routing -user postgres -clean      the proccess starts ok, but the way table will never have any data in it:               We need a way of type highway and class service     We need a way of type highway and class service     We DON'T need a way of type  and class     We DON'T need a way of type  and class     We DON'T need a way of type  and class     Split ways     NOTICE:  CREATE TABLE / PRIMARY KEY will create implicit index "nodes_pkey" for table           "nodes"         Nodes table created         NOTICE:  CREATE TABLE / PRIMARY KEY will create implicit index "ways_pkey" for table "ways"         Ways table created         Types table created         Classes table created     create topology     NOTICE:  CREATE TABLE will create implicit sequence "vertices_tmp_id_seq" for serial column "vertices_tmp.id"     CONTEXT:  SQL statement "CREATE TABLE vertices_tmp (id serial)"     PL/pgSQL function "assign_vertex_id" line 15 at EXECUTE statement     #########################     size of streets: 159747     size of splitted ways : 365597     finished      after a very long while (2 hrs) it sais finished but way table is empty! how do i fix this?   * i was able to load a routing table using osm2po. still no success with osm2pgrouting  Thanks for your help!


Can you create feature classes with the same name in different feature dataset

I've searched but can't seem to find any info on this which is surprising as this should be an easy one to answer. I'm using ArcGIS 9.3 and trying to create a feature class called 'cadastre' in multiple feature datasets within the same file geodatabase. ArcCatalog is giving me a 'The specified name already exists. Please enter a different name' error. It's been a while since I used ESRI software but I would have imagined this was possible since they are in seperate feature datasets. Is this expected funcitonality of the geodatabase?


How to get minx, miny, maxx and maxy of a geometry?

I'm new to Oracle Spatial. In MySQL or PostGIS, it's easy to get minx miny maxx maxy of a geometry, any and/or all of them. But how to do that in Oracle Spatial? Do I need to chain up several functions and properties to get any of them? Well this doesn't make sense to me. Thanks in advance.


I am unable to draw a shape layer

When I add a new shape layer or want to edit an existing one, nothing is happening. I even don't see the dots where I'm clicking. Maybe the problem is I deleted some layers (by opening the .qgs as a txt-file) as my project was collapsing by opening. I deleted three layers (maplayer and legendlayer, from start to end). Could there be some problem. The file is at https://dl.dropbox.com/u/83565282/AA_SR_Mendae.qgs When I open a new project everything is fine. Thanks a lot for all your help!! I'm using QGIS 1.8.0 Lisboa


Can I Dynamically Update Symbol Ranges?

**Background:** I'm currently trying to dynamically change the ranges for a layer with new style graduated symbology. **Problem:** I thought it may be as simple as updating the ranges for the existing symbology using the setLowerValue() and setUpperValue() methods on the QgsRendererRangeV2 class as below:               layer = qgis.util.iface.activeLayer()     rangePoint = 0     rangeIncrement = 100     for range in layer.rendererV2().ranges():             range.setLowerValue(rangePoint)             rangePoint += rangeIncrement             range.setUpperValue(rangePoint)      By my reconning the API documentation suggests that this should be possible. However, I am getting the error:               AttributeError: 'QgsRendererRangeV2' object has no attribute 'setLowerValue'      I've tried also using:               range.setLowerValue(float(rangePoint))      with no luck. What am I doing wrong?


What's the best method to get a weighted average for drive times?

I have a closest facility network analyses from the centroid of county to the closest medical facility. What's the best method to get the weighted averages of these drive times, using population as the weight. I am snapping all facilities and centroids to the road network, so now I am getting ridiculous averages like 6 hours. It is correct, but not feasible.


Why are poles missing from GIS topographic data?

Why are poles data missing from most topographic resources ? Data from 60⁰ to 90⁰ lat frequently missing from topographic GIS data. Example with SRTM4.1: ![enter image description here](http://i.stack.imgur.com/YQxmC.png) Common sense or deduction lead me to think that either:   1. **Orbits:** the orbit does not go to poles. Which is absurd for a topographic mission.   2. **Data :** the data on poles is messed up. But come on, in 2010's, space agencies computer scientists are not able to clean this up ? I asked myself this for years. Yet, I have no real, solid answer up to now. Just for fun, if the there are relevant GIS resources for these 2 areas, please share.


a hosted service to scrape and parse xml content and map it?

I am looking for a solution to map (on the web) some public data without scraping it first, esp. repeatedly scraping it to check for updates. Is this easy to do or only the data source itself can build the mapping into their web server themselves? Or to be clear: I understand what it would take to scrape the site myself and create a static map out of it by reuploading the parsed, structured data. I understand it is unlikely that someone is offering such a service ready-made. But still, perhaps worth asking. The concrete example is about Hungarian public procurement contracts (time constraint: election in 13 days), with contracts like this: http://kozbeszerzes.ceu.hu/tender/2013/00/2465.xml (labels are in Hungarian, but the XML elements and attributes are human-readable) The contractor and the bidders are linked, and they have addresses in attributes. (view the source for the entity as the first winner in the above example: http://kozbeszerzes.ceu.hu/entity/t/10011922.xml) It would be great to map the winning bids by contractor or bidder, maybe with some extra filters or something. Thanks for any thoughts!


QGIS 2.2.0.1 Valmiera, Windows 7 64 bit - copying out of GDALInfo-Window in QGIS not possible

When trying to copy information given out by GDALInfo with different kinds of raster data (DEM, normal georeferenced images etc.) through the interface of QGIS 2.2 by right-clicking and choosing copy (I want it in the clipboard) a Python error message keeps appearing. The steps to get this error are as follows: I go over Raster - Miscellaneous - Info in the menu bar. Then I right click the information that's shown in the resulting window and choose "copy all" or "copy". It's basically the equivalent of the QGIS GUI to entering gdalinfo and the location of the raster in the OSGEO4W shell. Only solution is to either enter all data manually in the program i want it in, or to use the OSGEO4W shell and copy the info about the raster from there (since it's not possible to do that out of properties of raster data in QGIS). This is the error message that shows up every time I try, was like this in QGIS 2.0, and it hasn't changed in 2.2:               Fehler bei der Ausführung von Python-Code:          Traceback (most recent call last):     File "C:/PROGRA~2/QGISVA~1/apps/qgis/./python/plugins\GdalTools\tools\doInfo.py",        line 73, in doCopyLine     output.append( r.text() + "\n" )     AttributeError: 'str' object has no attribute 'append'                         Python-Version:     2.7.4 (default, Apr  6 2013, 19:54:46) [MSC v.1500 32 bit (Intel)]               QGIS-Version:     2.2.0-Valmiera Valmiera, c3a2817          Python-Pfad: ['C:/PROGRA~2/QGISVA~1/apps/qgis/./python/plugins\\processing',        'C:/PROGRA~2/QGISVA~1/apps/qgis/./python', u'C:/Users/Lars/.qgis2/python', u'C:/Users      /Lars/.qgis2/python/plugins', 'C:/PROGRA~2/QGISVA~1/apps/qgis/./python/plugins',      'C:\\PROGRA~2\\QGISVA~1\\bin\\python27.zip', 'C:\\PROGRA~2\\QGISVA~1\\apps\\Python27       \\DLLs', 'C:\\PROGRA~2\\QGISVA~1\\apps\\Python27\\lib', 'C:\\PROGRA~2\\QGISVA~1\\apps     \\Python27\\lib\\plat-win', 'C:\\PROGRA~2\\QGISVA~1\\apps\\Python27\\lib\\lib-tk',      'C:\\PROGRA~2\\QGISVA~1\\bin', 'C:\\PROGRA~2\\QGISVA~1\\apps\\Python27', 'C:\\PROGRA~2     \\QGISVA~1\\apps\\Python27\\lib\\site-packages', 'C:\\PROGRA~2\\QGISVA~1\\apps     \\Python27\\lib\\site-packages\\PIL', 'C:\\PROGRA~2\\QGISVA~1\\apps\\Python27     \\lib\\site-packages\\win32', 'C:\\PROGRA~2\\QGISVA~1\\apps\\Python27\\lib\\site-     packages\\win32\\lib', 'C:\\PROGRA~2\\QGISVA~1\\apps\\Python27\\lib\\site-packages     \\Pythonwin', 'C:\\PROGRA~2\\QGISVA~1\\apps\\Python27\\lib\\site-packages\\Shapely-     1.2.18-py2.7-win32.egg', 'C:\\PROGRA~2\\QGISVA~1\\apps\\Python27\\lib\\site-packages     \\wx-2.8-msw-unicode', 'C:\\PROGRA~2\\QGISVA~1\\apps\\qgis\\python\\plugins\\fTools     \\tools']


Illegal argument exception while reading tiff image using GeoTools

I am working on a Java swings based GIS application to project my tiff image as maps . But each time I run my application I get `Illegal argument: "semiMajorAxis=�"` . The tiff images I have are working fine with QuantumGIS. Below is the complete exception log               WARNING: Illegal argument: "semiMajorAxis=�".     org.geotools.data.DataSourceException: Illegal argument: "semiMajorAxis=�".         at org.geotools.gce.geotiff.GeoTiffReader.<init>(GeoTiffReader.java:220)         at org.geotools.gce.geotiff.GeoTiffFormat.getReader(GeoTiffFormat.java:269)         at org.geotools.gce.geotiff.GeoTiffFormat.getReader(GeoTiffFormat.java:236)         at org.geotools.gce.geotiff.GeoTiffFormat.getReader(GeoTiffFormat.java:78)         at org.geotools.tutorial.quickstart.ImageLab.displayLayers(ImageLab.java:94)         at org.geotools.tutorial.quickstart.ImageLab.getLayersAndDisplay(ImageLab.java:82)         at org.geotools.tutorial.quickstart.ImageLab.main(ImageLab.java:59)     Caused by: org.geotools.data.DataSourceException: Illegal argument: "semiMajorAxis=�".         at org.geotools.gce.geotiff.GeoTiffReader.getHRInfo(GeoTiffReader.java:343)         at org.geotools.gce.geotiff.GeoTiffReader.<init>(GeoTiffReader.java:207)         ... 6 more      Below is the code.It works fine with bluemarble.tif image:               AbstractGridFormat format = GridFormatFinder.findFormat( rasterFile );             reader = format.getReader(rasterFile); //exception is thrown here      Below is the gdalinfo of the tiff image:               Driver: GTiff/GeoTIFF     Files: 46N10copy.TIF     Size is 4186, 4186     Coordinate System is:     GEOGCS["Geographic (Longitude, Latitude); user-defined (non-standard); user-defined (non-standard)",         DATUM["unknown",             SPHEROID["unretrievable - using WGS84",6378137,298.257223563]],         PRIMEM["Greenwich",0],         UNIT["degree",0.0174532925199433]]     Origin = (75.500000000000000,22.750000000000000)     Pixel Size = (0.000059722885810,-0.000059722885807)     Metadata:       AREA_OR_POINT=Area       TIFFTAG_MAXSAMPLEVALUE=255       TIFFTAG_MINSAMPLEVALUE=0       TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)       TIFFTAG_XRESOLUTION=200       TIFFTAG_YRESOLUTION=200     Image Structure Metadata:       INTERLEAVE=BAND     Corner Coordinates:     Upper Left  (  75.5000000,  22.7500000) ( 75d30' 0.00"E, 22d45' 0.00"N)     Lower Left  (  75.5000000,  22.5000000) ( 75d30' 0.00"E, 22d30' 0.00"N)     Upper Right (  75.7500000,  22.7500000) ( 75d45' 0.00"E, 22d45' 0.00"N)     Lower Right (  75.7500000,  22.5000000) ( 75d45' 0.00"E, 22d30' 0.00"N)     Center      (  75.6250000,  22.6250000) ( 75d37'30.00"E, 22d37'30.00"N)     Band 1 Block=4186x1 Type=Byte, ColorInterp=Palette       Color Table (RGB with 256 entries)         0: 255,255,255,255         1: 132,115,115,255         2: 189,165,165,255         3: 57,41,41,255         4: 181,90,90,255         5: 115,49,49,255         6: 255,198,189,255         7: 239,222,214,255         8: 173,156,148,255         9: 99,82,74,255        10: 74,57,49,255        11: 189,173,165,255        12: 156,140,132,255        13: 90,74,66,255        14: 247,231,222,255        15: 181,165,156,255        16: 148,132,123,255        17: 49,33,24,255        18: 173,148,132,255        19: 239,222,206,255        20: 206,189,173,255        21: 132,115,99,255        22: 214,206,198,255        23: 156,148,140,255        24: 115,107,99,255        25: 222,206,189,255        26: 148,132,115,255        27: 123,107,90,255        28: 140,115,90,255        29: 140,132,123,255        30: 74,66,57,255        31: 231,222,206,255        32: 132,123,107,255        33: 222,214,198,255        34: 239,231,214,255        35: 181,173,156,255        36: 107,99,82,255        37: 181,173,148,255        38: 74,66,41,255        39: 107,107,99,255        40: 173,173,148,255        41: 173,173,140,255        42: 165,173,123,255        43: 165,173,148,255        44: 123,140,132,255        45: 123,181,181,255        46: 140,198,206,255        47: 140,181,189,255        48: 198,189,206,255        49: 90,66,107,255        50: 132,107,140,255        51: 107,82,115,255        52: 189,156,181,255        53: 206,189,198,255        54: 239,90,115,255        55: 206,90,107,255        56: 239,57,82,255        57: 198,115,123,255        58: 239,82,99,255        59: 206,66,82,255        60: 255,156,165,255        61: 247,123,132,255        62: 189,49,57,255        63: 0,0,0,255        64: 0,0,0,255        65: 0,0,0,255        66: 0,0,0,255        67: 0,0,0,255        68: 0,0,0,255        69: 0,0,0,255        70: 0,0,0,255        71: 0,0,0,255        72: 0,0,0,255        73: 0,0,0,255        74: 0,0,0,255        75: 0,0,0,255        76: 0,0,0,255        77: 0,0,0,255        78: 0,0,0,255        79: 0,0,0,255        80: 0,0,0,255        81: 0,0,0,255        82: 0,0,0,255        83: 0,0,0,255        84: 0,0,0,255        85: 0,0,0,255        86: 0,0,0,255        87: 0,0,0,255        88: 0,0,0,255        89: 0,0,0,255        90: 0,0,0,255        91: 0,0,0,255        92: 0,0,0,255        93: 0,0,0,255        94: 0,0,0,255        95: 0,0,0,255        96: 0,0,0,255        97: 0,0,0,255        98: 0,0,0,255        99: 0,0,0,255       100: 0,0,0,255       101: 0,0,0,255       102: 0,0,0,255       103: 0,0,0,255       104: 0,0,0,255       105: 0,0,0,255       106: 0,0,0,255       107: 0,0,0,255       108: 0,0,0,255       109: 0,0,0,255       110: 0,0,0,255       111: 0,0,0,255       112: 0,0,0,255       113: 0,0,0,255       114: 0,0,0,255       115: 0,0,0,255       116: 0,0,0,255       117: 0,0,0,255       118: 0,0,0,255       119: 0,0,0,255       120: 0,0,0,255       121: 0,0,0,255       122: 0,0,0,255       123: 0,0,0,255       124: 0,0,0,255       125: 0,0,0,255       126: 0,0,0,255       127: 0,0,0,255       128: 0,0,0,255       129: 0,0,0,255       130: 0,0,0,255       131: 0,0,0,255       132: 0,0,0,255       133: 0,0,0,255       134: 0,0,0,255       135: 0,0,0,255       136: 0,0,0,255       137: 0,0,0,255       138: 0,0,0,255       139: 0,0,0,255       140: 0,0,0,255       141: 0,0,0,255       142: 0,0,0,255       143: 0,0,0,255       144: 0,0,0,255       145: 0,0,0,255       146: 0,0,0,255       147: 0,0,0,255       148: 0,0,0,255       149: 0,0,0,255       150: 0,0,0,255       151: 0,0,0,255       152: 0,0,0,255       153: 0,0,0,255       154: 0,0,0,255       155: 0,0,0,255       156: 0,0,0,255       157: 0,0,0,255       158: 0,0,0,255       159: 0,0,0,255       160: 0,0,0,255       161: 0,0,0,255       162: 0,0,0,255       163: 0,0,0,255       164: 0,0,0,255       165: 0,0,0,255       166: 0,0,0,255       167: 0,0,0,255       168: 0,0,0,255       169: 0,0,0,255       170: 0,0,0,255       171: 0,0,0,255       172: 0,0,0,255       173: 0,0,0,255       174: 0,0,0,255       175: 0,0,0,255       176: 0,0,0,255       177: 0,0,0,255       178: 0,0,0,255       179: 0,0,0,255       180: 0,0,0,255       181: 0,0,0,255       182: 0,0,0,255       183: 0,0,0,255       184: 0,0,0,255       185: 0,0,0,255       186: 0,0,0,255       187: 0,0,0,255       188: 0,0,0,255       189: 0,0,0,255       190: 0,0,0,255       191: 0,0,0,255       192: 0,0,0,255       193: 0,0,0,255       194: 0,0,0,255       195: 0,0,0,255       196: 0,0,0,255       197: 0,0,0,255       198: 0,0,0,255       199: 0,0,0,255       200: 0,0,0,255       201: 0,0,0,255       202: 0,0,0,255       203: 0,0,0,255       204: 0,0,0,255       205: 0,0,0,255       206: 0,0,0,255       207: 0,0,0,255       208: 0,0,0,255       209: 0,0,0,255       210: 0,0,0,255       211: 0,0,0,255       212: 0,0,0,255       213: 0,0,0,255       214: 0,0,0,255       215: 0,0,0,255       216: 0,0,0,255       217: 0,0,0,255       218: 0,0,0,255       219: 0,0,0,255       220: 0,0,0,255       221: 0,0,0,255       222: 0,0,0,255       223: 0,0,0,255       224: 0,0,0,255       225: 0,0,0,255       226: 0,0,0,255       227: 0,0,0,255       228: 0,0,0,255       229: 0,0,0,255       230: 0,0,0,255       231: 0,0,0,255       232: 0,0,0,255       233: 0,0,0,255       234: 0,0,0,255       235: 0,0,0,255       236: 0,0,0,255       237: 0,0,0,255       238: 0,0,0,255       239: 0,0,0,255       240: 0,0,0,255       241: 0,0,0,255       242: 0,0,0,255       243: 0,0,0,255       244: 0,0,0,255       245: 0,0,0,255       246: 0,0,0,255       247: 0,0,0,255       248: 0,0,0,255       249: 0,0,0,255       250: 0,0,0,255       251: 0,0,0,255       252: 0,0,0,255       253: 0,0,0,255       254: 0,0,0,255       255: 0,0,0,255      All suggestions & solutions are welcome.


Alternatives to ArcGIS Image Server Extension to serve Mosaic Dataset?

I am wondering if anyone has any resources for alternatives to ESRI's ArcGIS Image Extension for Server? We are looking to serve JP2000 images that have been consolidated using the Mosaic Dataset Method in ArcGIS. any thoughts?


What does Image width is not in range Valid range is from 1 to 2048 mean?

We are facing an issue on ESRI Silverlight (v 3.1) and ArcGIS Server 10 SP4 based web application. Few of the users complaining that map is not being displayed in their machine. When I investigated fiddler log, map service request throws 404 error . Below is the detailed > _Error Specified format is not supported Code: 400 Image width is not in > range. Valid range is from 1 to 2048._ Users are using dual monitors and different browser zoom levels, so few people are OK and few others facing trouble in viewing map layers, how to fix this issue for different screen resolutions and browser zoom levels? I am looking for comprehensive solution for this issue.


leaflet show miles/kilometers

Googlemaps shows you how large an area is on a map by showing a small bar with the number of miles/kilometers at the bottem left. Does leaflet also support this somehow? I do not see it listed in the docs at leafletjs but it seems like something that the library should support.


How to configure mapproxy and mapnik for getfeatureinfo requests

I'm wondering if its possible to configure `mapproxy` and `mapnik` for `getfeatureinfo` requests. I've made this combination work for `getmap` requests, but when I request `getfeatureinfo` I receive `layer xxx is not queryable`. I've tried setting `featureinfo:true`, but no luck there. Any help appreciated.


Has anybody used arcpy.mapping.ExportReport successfully on 64bit machine?

The **online help for arcpy.mapping.ExportReport** says that it has some restrictions on use but does not list **64bit** machines amongst them. I am running the following test from **IDLE** under **ArcGIS 10.1 SP1** (with **64bit Background Geoprocessing** installed) on **Windows 7 SP1**. I've create a simple test.mxd (one layer called Suburbs in data frame called Layers) and a simple test.rlf (included Suburbs, said Finish to take all defaults, saw report and saved it as *.rlf) However, this code (more or less straight from the Help example):               import arcpy     mxd = arcpy.mapping.MapDocument(r"C:\temp\test.mxd")     df = arcpy.mapping.ListDataFrames(mxd, "Layers")[0]     lyr = arcpy.mapping.ListLayers(mxd, "Suburbs", df)[0]     arcpy.mapping.ExportReport(lyr,                                r"C:\temp\test.rlf",                                r"C:\temp\test.pdf")     del mxd      throws this error:               Python 2.7.2 (default, Jun 12 2011, 14:24:46) [MSC v.1500 64 bit (AMD64)] on win32     Type "copyright", "credits" or "license()" for more information.     >>> ================================ RESTART ================================     >>>           Traceback (most recent call last):       File "C:\temp\test.py", line 7, in <module>         r"C:\temp\test.pdf")       File "C:\Program Files (x86)\ArcGIS\Desktop10.1\arcpy\arcpy\utils.py", line 181, in fn_         return fn(*args, **kw)       File "C:\Program Files (x86)\ArcGIS\Desktop10.1\arcpy\arcpy\mapping.py", line 515, in ExportReport         return report_source._arc_object.ExportReport(*gp_fixargs((report_layout_file, output_file, dataset_option, report_title, starting_page_number, page_range, report_definition_query, extent, field_map), True))     RuntimeError: Error in generating report     >>>       Any thoughts on cause and/or possible workaround are welcome!   (My suspicion is that ExportReport may not be supported on 64bit machines.)


Why can't a user with db_owner role version data?

I've made an enterprise geodatabase and added a user with db_owner permissions. When I add a database connection in ArcCatalog and connect to this geodatabase using the user I created, I'm not able to version the data. When I right click on a table or dataset in the enterprise geodatabase in ArcCatalog the menu is all grayed out except for the "Locks..." selection. I've added the user to other enterprise geodatabases with the same permissions and everything worked fine. The user was able to version data. What could be wrong here? If I sign into the geodatabase that isn't working with the SDE login, I am able to version the data.


QGIS GUI: How to get independant windows?

I am working with the georeferencer plugin, and I would like to have it working in fully independent windows. Is it possible ? Thanks for your help. _Context:_ I am working with qgis 2.4 under Ubuntu 12.04. I spot a dock/undock button on the top right corner of the plugin window. However even when undock is pressed, the windows remains linked to the qgis window and cannot be moved freely or displayed in full screen. * * * Here is as overview of the windows layout ![Windows Overview](http://i.stack.imgur.com/R0m5T.png)


Editing Feature Class using ArcObjects

I have a problem with editing feature class using arc Objects.Indeed I am working on it for a while and cant find any solution into doing it..I asked a question before about how to connect to database and I thought the problem is with it but I find another problem..Here is my code:                          // getLayerFromGeoDatabase();             // Create the workspace factory.             Type factoryType = Type.GetTypeFromProgID("esriDataSourcesGDB.SdeWorkspaceFactory");             IWorkspaceFactory workspaceFactory = (IWorkspaceFactory)Activator.CreateInstance(factoryType);             IWorkspaceFactory2 workspaceFactory2 = workspaceFactory as IWorkspaceFactory2;             // Build a connection string.             String[] connectionProps =      {         "dbclient=SQLServer", "serverinstance=CRAZYVIRUS-HP",              "database=test", "authentication_mode=OSA"      };             String connString = String.Join(";", connectionProps);             IPropertySet propertySet = new PropertySetClass();             propertySet.SetProperty("SERVER", "CRAZYVIRUS-HP");             //propertySet.SetProperty("INSTANCE", "CRAZYVIRUS-HP");             propertySet.SetProperty("DATABASE", "land_gis");             propertySet.SetProperty("authentication_mode", "OSA");             //propertySet.SetProperty("PASSWORD", "123456");             //propertySet.SetProperty("VERSION", "sde.DEFAULT");                  //IWorkspace workspace = workspaceFactory2.OpenFromString(propertySet, 0);             IWorkspace workspace = workspaceFactory2.Open(propertySet, 0);             ESRI.ArcGIS.Geodatabase.IFeatureWorkspace featureWorkspace = (ESRI.ArcGIS.Geodatabase.IFeatureWorkspace)workspace; // Explict Cast             ESRI.ArcGIS.Geodatabase.IFeatureClass featureClass = featureWorkspace.OpenFeatureClass("test.DBO.test");                  ESRI.ArcGIS.Carto.IFeatureLayer featureLayer = new ESRI.ArcGIS.Carto.FeatureLayerClass();                       featureLayer.FeatureClass = featureClass;                            featureLayer.Name = featureClass.AliasName;             featureLayer.Visible = true;                  IWorkspaceEdit workspaceEdit = (IWorkspaceEdit)workspace;                  // Start an edit session. An undo/redo stack isn't necessary in this case.             workspaceEdit.StartEditing(false);                  // Start an edit operation.             workspaceEdit.StartEditOperation();                  // Create a row. The row's attribute values should be set here and if             // a feature is being created, the shape should be set as well.             IFeature feature = featureClass.CreateFeature();             feature.Store();                  // Save the edit operation. To cancel an edit operation, the AbortEditOperation             // method can be used.             workspaceEdit.StopEditOperation();                  // Stop the edit session. The saveEdits parameter indicates the edit session             // will be committed.             workspaceEdit.StopEditing(true);      I can connect to database but cant add features in it.I got a error in this line               IFeature feature = featureClass.CreateFeature();      it says               enter code hereUnderlying DBMS error[[Microsoft][SQL Server Native Client 11.0][SQL Server]Violation of PRIMARY KEY constraint 'R7_pk'. Cannot insert duplicate key in object 'dbo.TEST'. The duplicate key value is (0).]      I can connect to database using arc catalog and also I have enabaled versioning for featureclass in ArcCatalog..I am really confused with this problem..I have already used Create Enterprise Geodatabase for this database and also my feature class is empty some times I got this error on the same line I think there is sth wrong with code..Some how this code creates a new feature but is says this new one is duplicated:( **Update** I managed to make this code work               Type factoryType = Type.GetTypeFromProgID("esriDataSourcesGDB.SdeWorkspaceFactory");             IWorkspaceFactory workspaceFactory = (IWorkspaceFactory)Activator.CreateInstance(factoryType);             IWorkspaceFactory2 workspaceFactory2 = workspaceFactory as IWorkspaceFactory2;                  IPropertySet propertySet = new PropertySetClass();             propertySet.SetProperty("dbclient", "SQLServer");                  propertySet.SetProperty("DATABASE", "land_gis");             propertySet.SetProperty("authentication_mode", "OSA");                  IWorkspace workspace = workspaceFactory2.Open(propertySet, 0);             ESRI.ArcGIS.Geodatabase.IFeatureWorkspace featureWorkspace = (ESRI.ArcGIS.Geodatabase.IFeatureWorkspace)workspace; // Explict Cast             ESRI.ArcGIS.Geodatabase.IFeatureClass featureClass = featureWorkspace.OpenFeatureClass("test.DBO.test");                  ESRI.ArcGIS.Carto.IFeatureLayer featureLayer = new ESRI.ArcGIS.Carto.FeatureLayerClass();                       featureLayer.FeatureClass = featureClass;                            featureLayer.Name = featureClass.AliasName;             featureLayer.Visible = true;                  IWorkspaceEdit workspaceEdit = (IWorkspaceEdit)workspace;                  workspaceEdit.StartEditing(false);                  workspaceEdit.StartEditOperation();                  IFeature feature = featureClass.CreateFeature();                  feature.Shape = polygon;                  ISubtypes subtypes = (ISubtypes)featureClass;             IRowSubtypes rowSubtypes = (IRowSubtypes)feature;                  rowSubtypes.InitDefaultValues();                  int EnergyFieldIndex = featureClass.FindField("Energy");             feature.set_Value(EnergyFieldIndex, property_info[1]);                  int corp_typeFieldIndex = featureClass.FindField("corp_type");             feature.set_Value(corp_typeFieldIndex, property_info[2]);                  int PoolFieldIndex = featureClass.FindField("Pool");             feature.set_Value(PoolFieldIndex, property_info[3]);                  int irrigationFieldIndex = featureClass.FindField("irrigation");             feature.set_Value(irrigationFieldIndex, property_info[4]);                  int varietyFieldIndex = featureClass.FindField("variety");             feature.set_Value(varietyFieldIndex, property_info[5]);                  int treesFieldIndex = featureClass.FindField("trees");             feature.set_Value(treesFieldIndex, Convert.ToInt16(property_info[6]));                       int utilizatFieldIndex = featureClass.FindField("utilizat");             feature.set_Value(utilizatFieldIndex, property_info[7]);                       int addressFieldIndex = featureClass.FindField("address");             feature.set_Value(addressFieldIndex, property_info[9]);                  int water_hourFieldIndex = featureClass.FindField("water_hour");             feature.set_Value(water_hourFieldIndex, Convert.ToDecimal(property_info[10]));                  int owner_IDFieldIndex = featureClass.FindField("owner_ID");             feature.set_Value(owner_IDFieldIndex, Convert.ToInt16(property_info[11]));                  int w_sourceFieldIndex = featureClass.FindField("w_source");             feature.set_Value(w_sourceFieldIndex, property_info[12]);                  int w_incheFieldIndex = featureClass.FindField("w_inche");             feature.set_Value(w_incheFieldIndex, Convert.ToDecimal(property_info[13]));                       feature.Store();             workspaceEdit.StopEditOperation();             workspaceEdit.StopEditing(true);      I set feature.Shape = polygon; using this piece of code                   ESRI.ArcGIS.Geodatabase.IWorkspaceFactory workspaceFactory = new ESRI.ArcGIS.DataSourcesFile.ShapefileWorkspaceFactoryClass();     ESRI.ArcGIS.Geodatabase.IWorkspace workspace = workspaceFactory.OpenFromFile(string_ShapefileDirectory, 0);     ESRI.ArcGIS.Geodatabase.IFeatureWorkspace featureWorkspace = (ESRI.ArcGIS.Geodatabase.IFeatureWorkspace)workspace; // Explict Cast     ESRI.ArcGIS.Geodatabase.IFeatureClass featureClass = featureWorkspace.OpenFeatureClass(string_ShapefileName);     Polygon p2 = featureClass.GetFeature(0).Shape as Polygon;     SelectedFeature = p2 as IPolygon;      I get an error that says `Arithmetic overflow error converting expression to numeric` I think the problem is with one of myfields..How can I find which one makes this problem..I know that this error says that a field in the database is attempted to insert or update larger data of the same type or other data type but cant manage to find which one is it


Moving object of fixed dimension without distortion

I am trying to move a polygon object in a shapefile to a new location. I will perform an offset so that the process can be automated by shifting the object to a new geographic coordinate or moved by some defined distance from the current position. If the methods are different, let's focus on the shifting my rigid box to a new coordinate. The object has defined dimensions (say 10m x 10m) but the reference coordinate system is a geographic coordinate system. As such the dimensions change when I move it far from its current location. Facts to consider in the solution:   * ArcPy solutions will work   * Using ArcGIS Desktop 10.1   * A movement of the centroid would be acceptable, because it is the offset from the initial position that the object will be moved; alternatively it could be moved from any arbitrary point as long as the point is known   * The coordinate system is also arbitrary. The solutions can specify or be general in that regard.   * The process needs to be repeated 100s of times


Mapnik - Weird tiles rendering when used with spatialite

I use Mapnik in an GeoDjango App with Spatialite 4.0 as backend. There is something weird with the tiles that I suspect might be related with a bad handling of mbr with spatialite. If you look at those two images, the first is the tiles rendered in OpenLayers with a map.zoomToExtent on the features. If you look at the triangle in atlantic ocean you will see that the feature is not rendered entirely by Mapnik. The second image show the entire feature not tiled. ![enter image description here](http://i.stack.imgur.com/SpgkP.png)![enter image description here](http://i.stack.imgur.com/xHc0o.png) If I zoom in, the triangle continue to not be tiled correctly. If I zoom out a lot it finally appears. I use this query in Mapnik:               query = '(select ' + geometryField + ', "id" from "shapefile_feature" where' + ' shapefile_id=' + str(shapefile.id) + ') as geom'     datasource = mapnik.SQLite(file='C:\mygeosite\sqlite3\sql3.db',                                     table=query,                                     srid=4326,                                     geometry_field=geometryField)          featureLayer = mapnik.Layer("featureLayer")     featureLayer.datasource = datasource


Identify parcels that have two types of land use

I have two polygon layers in a geodatabase. One layer is land parcels. The other layer is land use. I want to identify the six largest parcels that contain both "upland forests" and "wetlands" (i.e. WHERE landuse = 4 AND 6 in a single parcels). How do I proceed???


How to move a value from field1 to field2 using Field Calculator?

Hi I would like to move value in field1 to field2. Purpose is to after that delete the column 'field1' to get a better column naming standard /Jan


Import garmin Topo map to QGIS

I would like to import my Garmin Topo Maps to QGIS. I want to use the Topo Map as a Layer to create Shapefiles. Can anybody tell me if it is possible to import the Garmin maps including waypoints as a vector layer in QGIS?


How to copy data from one table to another using cursors?

I'm going to be as descriptive as possible. I have a table (TA) with values which are null for some rows, filled out for others. I have a table (TB) with the null values filled out. Dec in TA is the shared value in TB (called Nulled in that table) I can't do a table join because this would duplicate the amount of nulled columns in the first table. The gist of what I want it to do is this: if CURR_TYPE (one of the columns) is null in TA, lookup its Dec value, go to TB, look for the same value under Nulled, copy TB's CURR_TYPE, and paste it onto TA's CURR_TYPE. This is what I have so far, and I've tried variations on things I've found here, but I get lost with how to tell python to LOOKUP in another field, copy and               # Import arcpy module     import arcpy               # Local variables:     TA = "F:...\\TA"     TB = "F:...\\TB"     TLDisNull = arcpy.AddFieldDelimiters(TA, "SponsorOrganization") + " IS NULL"          #Update rows     with arcpy.da.UpdateCursor(TA, "CURR_TYPE", TLDisNull) as u_cursor:             with arcpy.da.SearchCursor(TA, "Dec") as s_cursor1:                     for s_row in s_cursor1:                             with arcpy.da.SearchCursor(TB, "Nulled") as s_cursor2:                                     for s_row in s_cursor2


WFS Filter Buffer around a point

I'm trying to make a buffer around a geometry (point/s,linestring,etc) using **WFS** queries without luck. This is what I've been doing:                  <ogc:Filter>          ...                  <ogc:Function name="buffer">            <gml:Point>                 <gml:coordinates cs="," ts=""> x,y </gml:coordinates>                    </gml:Point>             <ogc:Literal>300</ogc:Literal>          </ogc:Function>                  ...           </ogc:Filter>      I read I could do it using Dwithin, but is it possible with the buffer? I'm pretty new with all this. Appreciate any help.


Is there a way to write single bands raster from multi-band layer faster in r?

I have 1000 multi-layer rasters and have to export each separate layer to a tif. Some friends recommended R to do the job. I used the following code:               library(raster)          namesVector<- c("layer1", "layer2", etc) #contain layer names          stack <- stack("multiLayerRaster.tif") #20 layers, total 2 GB          s <- unstack(stack)          for(i in seq_along(s)){writeRaster(s[[i]], file=namesVector[i])}      The code works fine, but it is very slow, it takes about 20 min to export layers from one multiLayerRaster. Is there anyway to make the process faster? Any advice is highly welcomed!


kml to kml conversion with XSL

I have a simple KML file witch I get with ogr2ogr where all placemark are 2D polygon like this one :               <Placemark>     <Style>     <LineStyle><color>ff0000ff</color></LineStyle>  <PolyStyle><fill>0</fill></PolyStyle>     </Style>     <ExtendedData><SchemaData schemaUrl="#BATI">     <SimpleData name="HAUTEUR">5</SimpleData>     </SchemaData></ExtendedData>     <Polygon>     <altitudeMode>relativeToGround</altitudeMode>     <outerBoundaryIs>     <LinearRing>     <altitudeMode>relativeToGround</altitudeMode>     <coordinates>     5.210274804171726,45.193994541274471,0 5.210273193448314,45.193992448578285,0 5.210244307501982,45.194002928913058,0 5.210153551080012,45.193930862320109,0 5.209984563035429,45.194038696922476,0 5.210076592373608,45.194110738456011,0 5.210146416870383,45.194072424624984,0 5.210176223631215,45.194092352007594,0 5.210274804171726,45.193994541274471,0     </coordinates>     </LinearRing>     </outerBoundaryIs>     </Polygon>     </Placemark>      I would like to convert all of them to polygon with extrusion. The value of z coordinates is define by the value of "HAUTEUR" extended data. Heres it's 5 meters :               <Placemark>     <Style>     <LineStyle><color>ff0000ff</color></LineStyle>  <PolyStyle><fill>0</fill></PolyStyle>     </Style>     <ExtendedData>     <SchemaData schemaUrl="#BATI">     <SimpleData name="HAUTEUR">5</SimpleData>     </SchemaData>     </ExtendedData>     <Polygon>     <extrude>1</extrude>     <altitudeMode>relativeToGround</altitudeMode>     <outerBoundaryIs>     <LinearRing>     <altitudeMode>relativeToGround</altitudeMode>     <coordinates>     5.210274804171726,45.19399454127447,5 5.210273193448314,45.19399244857829,5 5.210244307501982,45.19400292891305,5 5.210153551080012,45.19393086232011,5 5.209984563035429,45.19403869692248,55.210076592373608,45.194110738456,5 5.210146416870383,45.19407242462499,5 5.210176223631215,45.19409235200759,5 5.210274804171726,45.19399454127447,5      </coordinates>     </LinearRing>     </outerBoundaryIs>     </Polygon>     </Placemark>      I tought that a XSL Stylesheet is the best and the easiest way to process but I don't figure how to start. Does someone as ever manage it?