Is there a way to increase the search radius of nearest neighbors for the Nibble tool?

I am working with a raster map vegetation types. I have about 20 types with unique raster code values that can be sorted into three classes: 1 - vegetation types I want to expand using the nibble function (e.g. forest) 2 - vegetation types I consider static that I do not want to either nibble or be nibbled (e.g. water, developed areas) 3 - vegetation types that I want to eliminate through the nibble function (early seral stages, such as grass and chaparral) The first time I ran this I was under a deadline and didn't have time to figure out a better solution, and wasn't able to find one googling/searching. I'm hoping to find a way to make this a shorter geoprocessing slog. As inputs to the Nibble, I reclassified my veg layer so that the type 2 (static) values were nodata, since I didn't want them to nibble. That was my input layer. For the input MASK, I reclassified the original raster but made the type 3 values nodata, since they are the ones I want to eliminate through nibbling. The catch is that I tend to have large amounts of type 3 clustered together, and on top of that there are several places where the type 2 and type 3 values are clustered near each other. The upshot is that the nibble gave me an incomplete output. It ran correctly, as far as I could tell, but it would stop before nibbling away all of the masked values. I think the problem is that the nibble function has some predetermined search radius for the nearest neighbor, and if a cell is surrounded by type 2 and 3 cells then it can be so far away from a type 1 cell that the nibble doesn't reach it. However, looking into the tools, the whole Nibble tool is just one function that is called with those two inputs, so I can't modify the python code to change the search radius. I am not skilled enough in python to write a loop that would run the nibble over and over again. So, is there a way to get the Nibble to fill in more of the masked values at a time? Is there a place I can get the Nibble code to try and modify it? I would imagine at some point a distance threshold was specified. Thanks in advance! I hope my explanation was clear! Edit: More information on the data I am using               +-------------+-------------------+--+-------------+--------+-------------+--------+-------------+---------+     |    Name     | Integer in Raster |  |   Static    | Value  |  Dissolve   | Value  |  Nibblers   |  Value  |     +-------------+-------------------+--+-------------+--------+-------------+--------+-------------+---------+     | SMCM        |                 0 |  | URB         |      4 | MCH         |      3 | SMCM        |       0 |     | OCFW        |                 1 |  | BAR         |      8 | GRASS       |      7 | OCFW        |       1 |     | SMCX        |                 2 |  | WAT         |      9 | MCP         |     10 | SMCX        |       2 |     | MCH         |                 3 |  | AGR         |     11 | ADS         |     19 | MHW         |       5 |     | URB         |                 4 |  | MRIP        |     12 | ASPE        |     23 | BLU         |       6 |     | MHW         |                 5 |  | MED         |     13 | ASPW        |     24 | DFTO        |      14 |     | BLU         |                 6 |  | Pixel Count | 322591 | Pixel Count | 333625 | RF          |      15 |     | GRASS       |                 7 |  |             |        |             |        | LPN         |      16 |     | BAR         |                 8 |  |             |        |             |        | SCN         |      17 |     | WAT         |                 9 |  |             |        |             |        | WWP         |      18 |     | MCP         |                10 |  |             |        |             |        | YPN         |      20 |     | AGR         |                11 |  |             |        |             |        | SAGE        |      21 |     | MRIP        |                12 |  |             |        |             |        | CMM         |      22 |     | MED         |                13 |  |             |        |             |        | LSG         |      25 |     | DFTO        |                14 |  |             |        |             |        | Pixel Count | 3744816 |     | RF          |                15 |  |             |        |             |        |             |         |     | LPN         |                16 |  |             |        |             |        |             |         |     | SCN         |                17 |  |             |        |             |        |             |         |     | WWP         |                18 |  |             |        |             |        |             |         |     | ADS         |                19 |  |             |        |             |        |             |         |     | YPN         |                20 |  |             |        |             |        |             |         |     | SAGE        |                21 |  |             |        |             |        |             |         |     | CMM         |                22 |  |             |        |             |        |             |         |     | ASPE        |                23 |  |             |        |             |        |             |         |     | ASPW        |                24 |  |             |        |             |        |             |         |     | LSG         |                25 |  |             |        |             |        |             |         |     | Pixel Count |           4401032 |  |             |        |             |        |             |         |     +-------------+-------------------+--+-------------+--------+-------------+--------+-------------+---------+      I tried the first answer below and it works, but doesn't necessarily solve the problem. Some of my patches of veg that need to be dissolved are pretty large. Therefore, I need an enormous neighborhood to eliminate the whole thing. Given this, it would seem I'm looking at using the raster calculator in an iterative way as well. So, now my question becomes, which is better, and if both, in what order? For example, I could Nibble for a few iterations, then switch to the raster calculator and iterate that over some reasonable-sized neighborhood (not sure how big that would be - 10x10? 50x50) until I eliminated all the type 3 pixels. Or I could do things in the opposite way. Raster Calculator is a little faster, but does it do a better job with respect to landcover assessment if there is a large neighborhood?


How to insert a quotation mark through the Field Calculator of attribute table of a Shapefile?

I have a shapefile of polylines with information of diameter. But the problem is that i didn't enter the inch sign(") with the values of diameter. What would be the solution that qutation mark would place along with diameter values?


How do I force the label to display above the lines

How do I force the label to stand out above the red lines to make the label more readable ? Here is the sample of the image ![enter image description here](http://i.stack.imgur.com/C8ujG.jpg)


Which CRS fits "Lambert Conformal Conic Projection, Standard Parallels 46.40 and 49.20"?

I am trying to georeference a VFR chart in qgis and I can't figure out which CRS to use. Map projection - Lambert Conformal Conic Projection, Standard Parallels 46.40 and 49.20. Any ideas which one on the list will work?


Does anyone use both ArcGIS and Python(x,y)?

I'm trying to get ArcGIS 10.1 and **Python(x,y)** to peacefully co-exist with a unified installation of Python. Surely someone else is using both these tools... right? Pointers would be much appreciated. Here's where I've gotten to:   1. VERSION: Since ArcGIS 10.1 uses Python 2.7.2 (which contains NumPy 1.6.1), I went back to the pythonxy install mirrors and found version 2.7.2.0 of Python(x,y) which includes Python 2.7.2 and NumPy 1.6.1.   2. INSTALL: I crossed my fingers and just installed Python(x,y) after copying my original "C:\Python27\ArcGIS10.1\" folder to a safe location. It installs most of the files at "C:\Program Files (x86)\pythonxy\", and the Python files in "C:\Python27\". This is the same folder where my Arc install of Python lives, except Arc placed its version in the sub-folder "\ArcGIS10.1\". So to recap, I now have a "python.exe" and associated files in "C:\Python27\" as well as "C:\Python27\ArcGIS10.1\".   3. NEXT? Both Spyder and ArcGIS 10.1 seem to work so far but I'm guessing this won't last long. Two separate installs on the same machine seems like a bad move in principle. I wonder: Is it possible to change Arc's preferences to point to the version in the parent folder, and condense versions? I've gone through folder by folder and in almost every case the Python(x,y) install includes all the files that ArcGIS 10.1 Python does (and then some). The exception is that the Arc install sometimes includes '.pyc' and '.pyo' versions of some '.py' files. But I'm guessing this isn't a critical difference since both directories contain the actual '.py' file.)   4. Forgot to mention: I added "C:\Python27\ArcGIS10.0\" to the list of system paths, per the instructions in Method 2 at http://forums.arcgis.com/threads/34418-is-it-possible-to-install-other-python-distributions-on-top-of-the-arcgis. I really don't know what this does, just experimenting here... Please advise! If anyone can offer further guidance on this, it would be much appreciated!! (PS - I may post the same question at forums.arcgis.com, so if interested in the answer check there too) (PPS - Tried to add 'python(x,y)' as keyword but was unable to create a new keyword)


How to add QGIS Python API Intellisense support in Aptana Studio3/Pydev

I'm developing QGIS plugins using PyDev in Aptana Studio 3 and I would like to know if it's possible to make it aware of the QGIS API and add some "intellisense" functionnalities when programming. For example, I would like to see the available methods of an object when typing the dot "." after it... UPDATE: I also see this message: > module not in pythonpath When Hovering any QGIS object in code... Is it supposed to show me something more usefull? ![enter image description here](http://i.stack.imgur.com/SJRTz.png) Thanks a lot Etienne


Diffusion/Kernel Interpolation with barriers?

I need to interpolate water quality data in a rather small area with peninsulas, islands and developed coastline that will act as barriers between my testing/sampling points. I understand I have to use non euclidean method to do so, but does anyone know the difference between "Diffusion Interpolation with Barriers" and "Kernel interpolation with barriers"?


Using threading module with Arcpy: Manager/Worker architecture

I asked this question on Stack Overflow on how to setup a manager/worker scenario. But I am getting some error when trying to use threading on my Python install with ArcGIS 10 (just the python install from ArcGIS 10.0). Is anyone using threading successfully in arcpy. My end goal is to be able to poll an arcgisjobs directory on ArcGIS Server 10.0 and call worker processes for each new folder. Each new process should be called and set to go and return control back to the manager to call more processes before previous ones are finished.


What is the easy way for extracting a dted raster file?

i have some srtm data which was converted to dted file with globalmapper or arcmap ascii to raster. i want to extract raster from vector with easy way.. i have tried "extract by mask" method but it depends raster files for i/o. do you know any easy way extracting raster? A section is raster file B section is vector file... **B section is what i want to extract from A Raster File.** **Raster File:** ![enter image description here](http://i.stack.imgur.com/qbqPK.jpg) **Raster and Vector File:** ![enter image description here](http://i.stack.imgur.com/QM01z.jpg)


Promotion radius in miles around venues using PostGIS

Bear with me as I'm new to PostGIS and GIS in general. I also don't know the domain well enough to create a good title for this post! Here is what I'm trying to accomplish: I have Venues in my system. Each Venue has a promotion radius in miles. Given an arbitrary point (lat/lng), I want to find all the Venues which encompass the arbitrary point in their promotion radius. Here's what venues tables look like:               venues     id | name | latitude | longitude | promo_radius      Is this possible with PostGIS? Are there any resources that would help me further understand the query? Thanks!


Explode MULTILINESTRING into individual segments in PostGIS 1.5

I have a `MULTILINESTRING` in PostGIS consisting of several `LINESTRING`'s. Consider each `LINESTRING` to have a set of `POINT`'s > P1, P2, ... PN forming the line segments > L1, L2, ... LN-1 where LX is the line segment connecting PX to PX+1. **Using PostGIS 1.5 how can I explode the`MULTILINESTRING` into individual line segments of each `LINESTRING`?**


Geo and Population Coverage Calculation

We are calculating Population and geocoverage in mapinfo. We have calculated population coverage as an area overlap. **How do we calculate population coverage if we have population defined for each polygon?** ![coverage](http://i.stack.imgur.com/I4ju6.png) **1.** **Are there any particular method to calculate population coverage for each boundary?** **2\. Any documentation and detail for analysis (not only steps but reason) would be helpful.**


How to dynamically add all feature layers using ArcGIS JavaScript API?

I have a mapserver > http:// mymapserver/services/rest/xxx/mapserver/ which contains 7 layers with /mapserver /1 /through mapserver/7 If a define an ArcGISDynamicMapServiceLayer and 7 FeatureLayers in OnDemand Mode, everything gets displayed and I can query every layer based on the infoTemplate. There is however one small glitch: the ArcGISDynamicMapServiceLayer has 1-x default visible layers which are not queryable and have a slight offset. Is there a way, to define the ArcGISDynamicMapServiceLayer, retrieve all included FeatureLayers and add them on the fly to a) the basemap and b) the Layerlist, while retrieving c) the Legend? Bonuspoints for keeping the hints borderline simple; I am a javascript novice and I hate dojo. Kind regards this is the actual code:               var visible = [];     var ml = [];          require(["dijit/layout/BorderContainer", "dijit/layout/ContentPane", "dijit/layout/AccordionContainer", "esri/dijit/Legend", "dojo/_base/array", "dijit/form/Button"]);     require(["dgrid/OnDemandGrid", "dgrid/Selection", "dojo/store/Memory", "esri/map", "esri/layers/FeatureLayer", "dojo/_base/declare", "dojo/number", "dojo/parser", "dojo/domReady!", "dojo/on", "dojo/dom", ],     function (Grid, Selection, Memory, Map, FeatureLayer, declare, dojoNum, parser, array, on, dom) {              parser.parse();         window.infoTemplate = new esri.InfoTemplate("Attributes", "${*}");         window.map = new esri.Map("map", {             basemap: "topo",             center: [12.8, 51],             zoom: 9         });              on(dom.byId("btnStreets"), "click", function () {             window.map.setBasemap("streets");         });         on(dom.byId("btnSatellite"), "click", function () {             window.map.setBasemap("satellite");         });         on(dom.byId("btnHybrid"), "click", function () {             window.map.setBasemap("hybrid");         });         on(dom.byId("btnTopo"), "click", function () {             window.map.setBasemap("topo");         });         on(dom.byId("btnGray"), "click", function () {             window.map.setBasemap("gray");         });         on(dom.byId("btnNatGeo"), "click", function () {             window.map.setBasemap("national-geographic");         });              window.MapServiceURL = "http://theserver/ArcGIS/rest/services/xxx/MapServer/";         window.HWSUrl = "http://theserver/ArcGIS/rest/services/xxx/MapServer/1";         window.outFields = ["OBJECTID", "HOCHWERT", "SchadenNr", "STATUS"];              var fl = new FeatureLayer(window.HWSUrl, {             id: "schaden",             mode: esri.layers.FeatureLayer.MODE_ONDEMAND,             outFields: window.outFields,             infoTemplate: infoTemplate         });              /** excluded FeatureLayers 2-7 **/              ml = new esri.layers.ArcGISDynamicMapServiceLayer(window.MapServiceURL);              dojo.connect(fl, "onLoad", function (fl) {             fl.maxScale = 0;             fl.setSelectionSymbol(new esri.symbol.SimpleFillSymbol().setOutline(null).setColor("#AEC7E3"));         });              dojo.connect(map, "onLoad", function (map) {             // show the border container now that the dijits              // are rendered and the map has loaded             dojo.style(dijit.byId("container").domNode, "visibility", "visible");         });              //add the legend         dojo.connect(map, 'onLayersAddResult', function (results) {             var layerInfo = dojo.map(results, function (layer, index) {                 return { layer: layer.layer, title: layer.layer.name };             });             if (layerInfo.length > 0) {                 var legendDijit = new esri.dijit.Legend({                     map: map,                     layerInfos: layerInfo                 }, "legendDiv");                 legendDijit.startup();             }         });              if (ml.loaded) {             buildLayerList(ml);         } else {             dojo.connect(ml, "onLoad", buildLayerList);         }         map.addLayers([fl]);                    });          /*building LayList function*/     function buildLayerList(layer) {         var items = dojo.map(layer.layerInfos, function (info, index) {             if (info.defaultVisibility) {                 visible.push(info.id);             }             return "<input type='checkbox' class='list_item' " + (info.defaultVisibility ? "checked=checked" : "") + "' id='" + info.id + "' onclick='updateLayerVisibility();' />" + info.name + "
";         });         dojo.byId("layer_list").innerHTML = items.join(' ');              layer.setVisibleLayers(visible);         map.addLayer(layer);     }          function updateLayerVisibility() {         var inputs = dojo.query(".list_item"), input;              visible = [];              dojo.forEach(inputs, function (input) {             if (input.checked) {                 visible.push(input.id);             }         });              if (visible.length === 0) {             visible.push(-1);         }         ml.setVisibleLayers(visible);     }


How to animate raster time series using open-source GIS?

I have a bunch of rasters taken about 10 minutes apart, and I have a nice QGIS project that produces pretty enough maps of them. I'm wondering if there's a way to animate this raster series, with a color map and overlaid with (static) vector data, using open-source GIS. A QGIS solution would be ideal since I already have the project set up there, but other suggestions are welcome as well.


GeoEyes stereo Image

We have 4 stereo GeoEyes images. We work with ErdasImagine 2014 when we generate a blk file after reading the images we see this error: "Unsupported Compression type". Please help us. Best Regards


Points appear at different location when using google map overlay

I have used vector point from database and also polygon country vectors from database.Points appear with in this country vectors but as I switch to google earth map as base layer the points at smaller zoom level would appear outside boundary where as at other zoom level it would appear inside. I have tried to add projection too but in vain.Do I need to reproject google earth's base layer and How?               map = new OpenLayers.Map('map', {                         controls: [                             new OpenLayers.Control.Navigation(),                             new OpenLayers.Control.PanZoomBar(),                             new OpenLayers.Control.LayerSwitcher({'ascending':false}),                                                    new OpenLayers.Control.ScaleLine(),                             new OpenLayers.Control.MousePosition(),                             new OpenLayers.Control.OverviewMap(),                             new OpenLayers.Control.KeyboardDefaults()                         ],                                   projection: new OpenLayers.Projection("EPSG:900913"),                  displayProjection: new OpenLayers.Projection("EPSG:4326"),                  numZoomLevels: 21,                  maxExtent: new OpenLayers.Bounds(-20037508, -20037508,20037508, 20037508.34)                          });                                     var border = new OpenLayers.Layer.WMS( "Border",                     url,                      {layers: 'GisData:BORDER'}                 );                          var dsBand = new OpenLayers.Layer.WMS( "District Band",                     url,                      {layers: 'GisData:BND'}                 );                                                 var Site1 = new OpenLayers.Layer.WMS( "Site 1",                     url,                      {layers: 'GisData:sites', transparent: "true",format: "image/png"}                              );                     dwnSite3G.setIsBaseLayer(false);                          var Site2 = new OpenLayers.Layer.WMS( "Site 2",                     url,                      {layers: 'GisData:sites2', transparent: "true",format: "image/png"}                              );                               var g_pmap = new OpenLayers.Layer.Google(                     "Google Physical",                     {type: G_PHYSICAL_MAP}                     );                          var gmap_gs = new OpenLayers.Layer.Google(                         "Google Streets" // the default                          );                          var gmap_ghy  = new OpenLayers.Layer.Google(                         "Google Hybrid",                         {type: G_HYBRID_MAP}                     );                          var gmap_gsat = new OpenLayers.Layer.Google(                         "Google Satellite",                         {type: G_SATELLITE_MAP}                     )


QGIS - Export PDF map with georeferenced data

I am used to working with ArcMap but have been playing around with QGIS. When I export a PDF map with ArcMap, the PDF is spatially referenced. This is a handy tool, I can put the PDF map on my iPhone and it tracks my location on the map (even when I'm out in the woods away from cell coverage). It's very handy for forestry work. Are the PDF maps exported with QGIS also spatially referenced? Is there a plugin that allows for this feature?


Using SearchCursor to access & export values in a raster table

I am desperately trying to export the values in a raster table to a .txt file. I've seen some other folk post about this, but for some reason can't get the code to work in my case. I have a raster, with table already built, called "shed". I simply want to export two fields, "Value" and "Count", to a .txt or .csv file (called "try" in this case). So far, I haven't even been able to get a script that even looks at the raster's table. So far I have:               outfile=open("E:/Weathering_GIS/try.txt",'w')     rows=arcpy.SearchCursor('shed',"","","Value;Count","")     for row in rows:         val=row.getValue('Value')         count=row.getValue('Count')          outfile.write(val)      But "val" and "count" always fail to become anything- they just remain empty values! Am I missing something? Maybe some sort of import I have overlooked?


SVG iconed features not rendered correctly during zoom

Im using Geoserver with PostGIS database to put points to map using GeoExplorer. The points are styled to use SVG icon with fixed size (uom metre, size) so they scale together with the map as user zooms in and out. For now, there are no filters or any other issues. Style code :               <sld:PointSymbolizer uom="http://www.opengeospatial.org/se/units/metre">             <sld:Graphic>               <sld:ExternalGraphic>                 <sld:OnlineResource xmlns:xlink="http://www.w3.org/1999/xlink" xlink:type="simple" xlink:href="load_icon.php?id=${param}"/>                 <sld:Format>image/svg+xml</sld:Format>               </sld:ExternalGraphic>               <sld:Size>80000</sld:Size>             </sld:Graphic>           </sld:PointSymbolizer>           <sld:TextSymbolizer/>      load_icon.php is returning the SVG code for given param :               <svg xmlns="http://www.w3.org/2000/svg" version="1.1" width="220" height="110">     <g transform="scale(0.5)">     <g transform="translate(90,50)">                     <rect x="1.5" y="1.5" height="103" width="167" fill="#fff" stroke="red" stroke-width="3"/>     <path d="M3,3 l164,100 m0,-100 l-164,100" stroke="red" stroke-width="3" fill="none"/>     </g>     <text x="173" y="49" text-anchor="middle" font-family="Verdana" font-size="25" fill="red">XXXXX< /text>     <text x="270" y="153" font-family="Verdana" font-size="25" fill="red">G< /text>     </g>     </svg>      I use the plain GeoExplorer to view the icon on the map, but same issue is when connecting the layer as WMS to Leaflet and other libraries. The icon is showing ok in some resolution : ![icon ok](http://i.stack.imgur.com/zsgpN.png) but gets cropped as user is zooming in : ![icon croped](http://i.stack.imgur.com/MooMk.png) and even more cropped with further zoom : ![icon more cropped](http://i.stack.imgur.com/zupKx.png) with tiles that do not contain the part of the icon. There is caching turned off for this layer, no help. What I am doing wrong here ?


accessing mapinfo data without using mapinfo api

I need to be able to access various mapinfo layers but without using the mapinfo api. I have found that I can convert the mapinfo layers to an access database to get the attribution and I can use sharpMap to load and view the layer. The problem is now I'm not sure how to get a location shown in sharpmap and then query the objects to get to the attributes, or if its even possible. It seems that mapinfo might keep this kind of link internaly perhaps


How to install fTools and Table manager

I've just make clean install of Qgis Lisboa under Ubuntu Precise. Latest Qgis version under Windows XP is very unstable (crashes every few minuter while editing shp file). After new installation there are some problems for me:   1. I can't see or add independent repositories. So I can't add them under "Repositories" tab. Can I add them anyway? What I need is Table Manager at least. Possibility of install some other plugins could be usefull as well.   2. I can't find "Join attributes" under fTools menu (it was the same problem after upgrading to Lisboa under Win XP, but deleting fTools directory and reinstalling this plugin solve problem but only under Win XP). So, I can't work under Windows due to crashes nor under Ubuntu :( Could you help me to make Qgis full functional under Ubuntu? TIA Best regards, Cyfranek


Python: shp file to GML using OGR

I am trying to convert a shapefile to GML 3.1.1 . I am using gdal 1.9.1 however the GML conversion it generates for a geometry feature is GML 2 and it also swaps the long/latitude coordinates. Are there simple workarounds to keep the coordinate positions the same, while changing the output to GML 3.1.1? Thanks.


How to create TINs in ArcGIS from triangular elements

I am working with FEM models with triangular elements, and I'd like to view the results as TIN surfaces in ArcGIS (and in ArcScene). In my external source, I have **nodes** (with X,Y,Z coords) and **elements** (forming a triangle from three nodes), which are the fundamentals of a TIN. How can I bring these together in ArcGIS to convert my triangular elements into a new TIN surface? This should ideally be a 1:1 conversion, but the documentation isn't helping me. Here is an illustrated example of the data I have: ![Simple TIN](http://i.stack.imgur.com/TaS5A.png)               Nodes: n, X, Y, Z     1, 51.06, 84.98, -7.46     2, 30.45, 80.67, -7.47     3, 44.83, 72.70, 19.46     4, 64.99, 69.91, 61.48     5, 30.33, 63.40, 71.31     6, 46.95, 61.96, 48.44     7, 62.04, 51.97, 82.51     Elements: n1, n2, n3     6, 7, 4     4, 3, 6     3, 5, 6     5, 3, 2     1, 3, 4     2, 3, 1      I'm currently using ArcGIS 9.3 with the 3D Analyst Extension (among others), and would like some hints for a solution targeted for either Python or VBA. * * * _**Alternative approach_** : The raw data -> Esri TIN conversion doesn't look very easy, and I can't make heads or tails of the API for 3D Analyst, so I'll hold the Python/VBA solution off. I have the same data represented in a Shapefile where all features are triangular PolygonZM shapes and the vertices have the height-field. (If you ask, this Shapefile was made using OGR via Python from the original dataset structured similarly as above). Viewing this file in ArcScene is really really slow, as I have ~100k triangular polygons in the Shapefile. This PolygonZM Shapefile should ideally convert triangle-by-triangle to a TIN dataset. I have the data, so interpretation or interpolation is not wanted! So, how can I convert this PolygonZM Shapefile (pictured below) into a TIN? In ArcScene, in the "3D Analyst" menu, there is a "Create TIN From Features" tool, but I don't understand the language "soft clip" etc. The dialog provides no help. ![ArcScene: Create TIN From Features](http://i.stack.imgur.com/qaCeA.png) * * * # LandXML It turns out the best format is LandXML, thanks to the related answers below. For the example above, the file looks like this (note: the coordinate order is Y X Z):               <?xml version="1.0" ?>     <LandXML version="1.2" xmlns="http://www.landxml.org/schema/LandXML-1.2" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://www.landxml.org/schema/LandXML-1.2 http://www.landxml.org/schema/LandXML-1.2/LandXML-1.2.xsd">       <Surfaces>         <Surface name="My TIN">           <Definition surfType="TIN">             <Pnts>
84.98 51.06 -7.46

80.67 30.45 -7.47

72.70 44.83 19.46

69.91 64.99 61.48

63.40 30.33 71.31

61.96 46.95 48.44

51.97 62.04 82.51
             </Pnts>             <Faces>               <F id="1">6 7 4</F>               <F id="2">4 3 6</F>               <F id="3">3 5 6</F>               <F id="4">5 3 2</F>               <F id="5">1 3 4</F>               <F id="6">2 3 1</F>             </Faces>           </Definition>         </Surface>       </Surfaces>     </LandXML>      Then use the LandXMLToTin_3d tool to process the result to an Esri TIN file: ![Resulting TIN](http://i.stack.imgur.com/SmMke.png)


How can you make RasterToASCII_conversion observe environment setting output cell size?

I have two rasters (an ESRI grid and an IMAGINE image) that differ in projection, extent and cell size. I want to convert the IMAGINE to an ASCII, and want it to have the extent, cell size and projection of the ESRI grid. I've set global environment variables (extent, output CRS, cell size, snap raster) to match the values in the ESRI grid, and the ASCII that results from RasterToASCII_conversion has the correct projection and extent, but ignores the cell size setting. However, if I use the Raster and ApplyEnvironment tools (see code below) from the Python console, the output is as expected (correct cell size, extent, etc).               import arcpy     arcpy.env.extent = 'c:/path/to/grid'     arcpy.env.cellSize = 'c:/path/to/grid'     arcpy.env.outputCoordinateSystem = 'c:/path/to/grid'     arcpy.env.snapRaster = 'c:/path/to/grid'     r = arcpy.Raster('c:/path/to/imagine')     arcpy.CheckoutExtension('Spatial')     rclip = arcpy.sa.ApplyEnvironment(r)     rclip.save('c:/path/to/outfile')      While I could do the above and then convert the resulting grid to ASCII, I'd prefer to do the whole procedure without using Spatial Analyst. Is this possible? I'm using ArcGIS 10.


How to "lower the resolution" of a lat long csv

I have a csv with 12.5 million date/lat/long/value points, it's a point every 800 meters ( in a grid ) of the united states what's the best way to convert the points to every 1600 meters? or 3200? etc... Really need to reduce the file size


Spatially enabled Postgres database

I am trying to create a spatially enabled PostGIS database. I am following the PostGIS documentation, http://postgis.net/docs/manual-1.5/ch02.html#id2648455. In the short version, i executed,               createdb yourdatabase     createlang plpgsql yourdatabase      But I did not execute the last three commands,               psql -d yourdatabase -f postgis.sql     psql -d yourdatabase -f postgis_comments.sql     psql -d yourdatabase -f spatial_ref_sys.sql      Is my database still spatially enabled? If not, then can you tell me where these .sql files are located so I can execute those last commands, since I am getting the following error,               postgis.sql: No such file or directory      I am using Ubuntu 12.04


Qgis 1.8 Lisboa does not have Join attributes (only Join attributes by location)

I installed 1.8 Lisboa. I tried to join to the attribute table more data (more columns ), but Vector/Datamanagement Tools/ Join attributes by location there were also not "Join attributes" alternative. b) I installed 1.4 and there was Vector/Datamanagement Tools/ Join attributes 2) what to do, when attribute table has only Nulls, how to edit, the edit button does not work nor 1.4 and 1.8 operation system: Windows XP IBM ThinkPad T60


How to rotate a feature in openLayers

I have an svg image and I want to rotate it around itself so I do the following:               var origin = new OpenLayers.Geometry.Point(vertices[0].x,vertices[0].y);                 feature.geometry.rotate(90 ,origin);     feature.layer.drawFeature(feature);      but nothing happens. Does anyone know why? Thanks in advance.


Is there a workaround for 80+ layers in OpenLayers?

My Map web app reads client list from SQL server and for each client it draw points on a different layer so that the user could turn on and off that layer and see his entire customers. Problem is that once I have around 80+ layer than popups are appearing below the layers. Is there a workaround how I can show different features and turn them on / off without screwing with the popups ?


How to cache PostGIS views?

Does anyone knows if there exists a caching system for views or queries in postgis? I have a system that gives information through a webpage using a postgis view as backend. The problem is that, sometimes the view is big enough that takes 5 seconds per user to serve the info. I know there's this system named geocache that works very good for raster data using wms protocol. I was wondering if there's a similar tool for vector / table data.


How to convert polygon layer to txt files?

Is there any way to convert polygons that do not have coordinates in the attribute table to a txt file with coordinates. I will be using the txt file in a different software (not ArcGIS), so I need the txt file with polygon coordinates. Hope someone can help me. Thank you in advance. Lina


Retrieving routedescription geometry as GeoJSON

I'm building a small web-map and would wich to display a polyline representing a specific route via roads. With leafletjs it is easy to display the geometry on/over the map. So that's not the problem. I'm looking for a webservice where I can 'download' the route as a polyline/Multipoint. Does anyone know if this exists?


Can esri ipad app consume web mapping service?

Can the ESRI app for the iPad consume regular web mapping services, or do I need to publish services for this app as Mobile project services?


When a class in an ArcGIS extension should be an ArcGIS class

Visual Studio templates provided by the ArcGIS SDK offer a "ArcGIS Class" which is a simple C# class that has the attributes Guid, ClassInterface, and ProgId before the class declaration.   While adding a new class to an ArcGIS extension project, when should I use an ArcGIS class and when should I use a regular C# class?


Traceback (most recent call last):       File "/home/kasun/.qgis2/python/plugins/floogin/floogin.py", line 246, in findRoutes         path=self.astar.astarSearch(graph, start, end)       File "/home/kasun/.qgis2/python/plugins/floogin/astar.py", line 110, in astarSearch         for tile in graph[current]:     TypeError: 'QgsGraph' object does not support indexing      Above Error occured when I run the following code. I'm new to qgis and I don't know why the error occurs.I wanted to find path between two selected points.I would really appreciate if you could please help me out to find solution to the problem. Thank You.               def astarSearch(self, graph, current, end):         openSet = set()         openHeap = []         closedSet = set()              def retracePath(c):             path = [c]             while c.parent is not None:                     c = c.parent                     path.append(c)             path.reverse()             return path              openSet.add(current)         openHeap.append((0,current))         while openSet:             current = heapq.heappop(openHeap)[1]             if current == end:                     return retracePath(current)             openSet.remove(current)             closedSet.add(current)             for tile in graph[current]:                     if tile not in closedSet:                         tile.H = (abs(end.x-tile.x)+abs(end.y-tile.y))*10                          if tile not in openSet:                                 openSet.add(tile)                                 heapq.heappush(openHeap, (tile.H,tile))                         tile.parent = current         return []               def findRoutes(self):     vl=self.canvas.currentLayer()     director = QgsLineVectorLayerDirector( vl, -1, '', '', '', 3 )     properter = QgsDistanceArcProperter()     director.addProperter( properter )     crs = self.canvas.mapRenderer().destinationCrs()     builder = QgsGraphBuilder( crs )          global x1     global y1     global x2     global y2     pStart = QgsPoint( x1, y1 )     pStop = QgsPoint( x2, y2 )          tiedPoints = director.makeGraph( builder, [ pStart, pStop ] )     graph = builder.graph()     tStart = tiedPoints[ 0 ]     tStop = tiedPoints[ 1 ]     idStart = graph.findVertex( tStart )     idStop = graph.findVertex( tStop )          start, end = idStart, idStop     path=self.astarSearch(graph, start, end)     if path is None:             print "No path found"     else:         print "Path found:", path


QGIS is not working with Postgres Plus Advanced Server

We have a problem to integrate GIS software with our GIS database which is populated on Postgres Plus Advanced Server 9.2 from EnterpriseDB. We are using QGIS to maintain and analyze our GIS repository, the software is working just fine when the database in the postgres community version, but when the database in ppas; QGIS return the following error while trying to fetch the list of GIS layers: > Erroneous query: >   >   >      SELECT > l.f_table_name,l.f_table_schema,l.f_geometry_column,upper(l.type),l.srid,c.relkind >      FROM geometry_columns l,pg_class c,pg_namespace n >      WHERE c.relname=l.f_table_name >        AND l.f_table_schema=n.nspname >        AND n.oid=c.relnamespace >        AND has_schema_privilege(n.nspname,'usage') >        AND > has_table_privilege('"'||n.nspname||'"."'||c.relname||'"','select') >      ORDER BY n.nspname,c.relname,l.f_geometry_column >   > > returned 7 [ERROR: schema "dbms_profiler" does not exist] And the following is the answer from EnterpriseDB: > We looked at query and it seems, issue is the way that query built. By > default, PPAS stores package names in namespace table and so it is > retrieving the names of packages as name of schemas hence query is failing. > It works in PostgreSQL. Though, you are in postgres compatibility, we still > retains some packages. > > So we would recommend you to contact QGIS vendor on the same as it is not > tested with PPAS. Your help will be highly appreciated.


Flannery compensation in Leaflet JS

I'm creating a bubble map using Leaflet JS. Everything is working smoothly (LeafletJS is awesome!) but now I'd like to try flannery compensation on the bubbles. I can't find the math anywhere though. Best explainations I've been able to find say you have to multiply radius by either 1.4 or 1.8, which just seems a bit useless. What is the math used for Flannery compensation?


Snag loading Gdal 1.11 for Qgis 2.2, on a Mac

Can't load Q2.2 on a Mac without Gdal 1.11. but when I try to load 1.11 'Complete' version, it fails and gives me an error message to contact the developer. Anyone else had this kind of trouble? David A UK


csv to INSPIRE gml problem in FME Desktop 2014

I have .csv file with id, name, lat, long and other. I just cannot make gml writer to write geometry tags. I mean something like . How can I read this things from my .cvs? What am I missing? fme is new for me... Thank you.


How to extract vectors from map using QGIS?

I am trying to extract building footprints from a raster map (as shown in the example below). Please note that the maps are out of copyright, have been scanned and georectified by a library and have been provided to me for exactly this purpose. ![Scaled example - original geotif available on request](http://i.stack.imgur.com/fQotP.jpg) So far I have been using the NYPL Map-Vectorizer which works by first adjusting the brightness/contrast of the image, then setting a colour threshold to extract a black and white image, which in turn is processed using gdal_polygonize.py (followed by some simplification and removal of roads). The issue is that this gives polygons of the _inside_ of the building. Terrace buildings end up with a gap between them. **How can I post-process to remove the gap between the terrace buildings? Is there an altogether better method for extracting building footprints from this map?** Note: Solutions need to be able to run as a script. I'm able to test proposed techniques in QGIS.


Automating Intersect between two polygon layers

I have a spreadsheet that has county-level occurrences of around 1000 species within the US in the column headings (Species State County). Thus, for some counties there may be multiple species. I also have a county shapefile with the County field linking the shapefile to the spreadsheet. Finally I have a shapefile with multiple, separate polygons and I want to know which rows from the spreadsheet intersect or contain a polygon from the second polygon shapefile. I am looking for the best way to automate this process so that I create a new column in the spreadsheet that has a 1 if the county occurrence of the species intesects with the second polygon shapefile. For example, Species State County Intersect sp. A NC Durham 1 sp. A NC Orange 0 where species A county-level occurence intersects with the polygon in Durham County, but not in Orange County. I have not specified a platform, but I have worked with R, ArcGIS, and Python. Currently I am most comfortable working in R (but not for spatial data, there I am a newbie), but not sure how to automate this process. thanks for any information.


Weighting a data frame within another frame for map creation

I've currently got an inset data frame sitting in a corner of a larger data frame. Some labels are being half covered by the smaller data frame. I was wondering if there are any tips on how to weight it (either in Maplex or Standard) so that the labels that are half covered are either removed, or moved so that they aren't behind the smaller data frame? Thanks


Remove overlap in raster using ERASE command line?

I am trying to erase a raster and is there a erase command for the Raster side? I am using ArcGIS 9.3.1. I have 2 rasters that are in the Landcover and there is overlap in between and would like to remove the overlap and then merge them later. Thanks, Robert


Waiting for Vector & WFS loading

_user borrowed from friend_ since i use vector layer with wfs protocol, there is a little while before features loaded on the map. So i need to tell user to wait a moment until they load. what best way to do that ? i found this http://osgeo- org.1560.n6.nabble.com/Patch-quot-waiting-quot-mouse-cursor-on-WFS-loading- td3961983.html but not sure if it will work. any body tried it ?


Trying to make the private land white overlay hillshade but

I'd like to have anyone's feedback on this draft map I created. Tried to make the private ownership white but since it's RGB is not near to 255. I've seen someone created this map similiar to it and leaves it as light grey color. Thanks. http://www.cartotalk.com/index.php?showtopic=7125 ![Gray Image](http://i.stack.imgur.com/RqiR0.jpg) I am using ArcMap 9.3.1. I have Photoshop 7 on my computer but no training background on it. Have use couple of times in the past not related to GIS.


Is it possible to use a Multipatch to implement IFeatureRenderer.Draw in a custom renderer

I have a custom renderer which need to render modified geometry. The input features are all coming in via the `IFeatureCursor` as polygon elements. (This is similar to this example, though I'm not rendering points.) I can easily render Polygons, Lines, and Points via either passing the new geometry into `IFeatureDraw.Draw` or calling `IDisplay.DrawPolygon`/`IDisplay.DrawPoint`/etc. However, `IDisplay` does not provide a mechanism for drawing a multipatch, and passing this in as the geometry to `IFeatureDraw.Draw` results in nothing drawn to the screen. Is there a way to directly "draw" a multipatch element (preferrably with texturing), within this method?


Does a geotagged image contain information about its accuracy?

When i view my current location on Google maps on my android it shows me a circle of uncertainty. If i take a geotagged picture it will off course contain information about its Lat or Lon. Does it also contain information concerning accuracy of this geotagging ?


Toggle editing option in plugins

I am working on python plugins for QGIS.I developed the plugins same as RT SQL LAYER plugins. When we draw the layers using query builder of " **Add Postgis Layer** " option, **toggle editing** option is available for the same. But when i try to draw layers using **RT SQL LAYER** plugin ,it is not available for editing. Same happens with my plugins too. How do I make **toggle editing** option available for my plugins.?? My plugins are in **python**.


extract vertices from dbgeometry value

Who knows the elegant way to extract coordinates of polygon vertices from well known text representation of DbGeometry value? Here's my polygon for example:                {POLYGON ((52.808019 -1.345367, 53.095124 -1.345367, 53.095124 -0.864716, 52.808019 -0.864716, 52.808019 -1.345367))}      Thanks in advance


New Install Quantum GIS - .ind and .idm files not valid

I downloaded and installed the Windows version of Quantum GIS to view Pinellas County data. I can add the .shp, .shx, and .dbf layers, however do not see any data showing as far as road names, etc. I tried to add the .idm and .ind layers that come with the Pinellas County downloads and they show as not valid. Do I need those to see the data on the map? If so, how do I add them.


django + apache mod_tile

I'm not sure if I'm doing this right and I need your opinion on this current situation. I'm creating a Web Interface using django and postgres as the database. With regards to the slippy map, I went with using apache and mod_tile. The tiles where first designed in tilemill and then used the generated xml file with the generate_tiles.py. Am I doing this right or are there other alternatives? Just in case, the web interface will be used to draw lines on the map as routes and store them for querying.


How to install SciPy for QGIS Standalone

A number of QGIS plugins require the Python package SciPy to be installed. This package is available using the OSGeo4W package based installer or precompiled for Windows, but is not provided as part of the Windows standalone installer. The SciPy windows installer does not recognize the Python instance installed with the QGIS standalone install, and so will not install to it. Can anyone provide a simple solution to installing the SciPy package into the QGIS standalone python instance that is suitable for an inexperienced python user? Any help would be really appreciated.


How to generate StreetFrom and StreetTo data from street segments?

## **Background:** I have a polyline feature class of about 4000 street centerline segments for which I would like to add data on what streets the segments intersect with. I would like to add some new fields to my dataset that capture information on what street the polyline segment begins at and what street the polyline segment ends at (essentially the intersections) ## **What I have so far:** I wrote a quick Python script to calculate these values by doing Select By Location -> Boundary Touches:               # Import arcpy and set path to data     import arcpy     arcpy.env.workspace = "c:\Users\me\Documents\ArcGIS\Default.gdb"          #create a master search layer to encompass all streets.      searchLayer = arcpy.MakeFeatureLayer_management('StreetsCP', 'searchLayer')      #set the fields that we are going to be updating     fields = ('FULL_NAME', 'OBJECTID', 'StreetFrom', 'StreetTo')     #create a cursor to encompass all rows in streets     with arcpy.da.UpdateCursor('StreetsCP', fields) as uCursor:             #loop over every row in streets             for uRow in uCursor:                 #set where clause to create a feature layer from the row                   whereClause = '"OBJECTID" = ' + str(uRow[1])                 #create feature layer consisting of just the row we are on                 inputLayer = arcpy.MakeFeatureLayer_management('StreetsCP', 'inputLayer', whereClause)                       #select the streets that touch the boundary of the street we are on                 arcpy.SelectLayerByLocation_management(searchLayer, "BOUNDARY_TOUCHES", inputLayer, "", "NEW_SELECTION")                      touchingStreets = []                 uniqueStreets = []                      #populate the unique streets list                 with arcpy.da.SearchCursor(searchLayer, 'FULL_NAME') as sCursor:                     for sRow in sCursor:                         if sRow[0] != uRow[0]:                             touchingStreets.append(sRow[0])                     #eliminate duplicate street values                     uniqueStreets = list(set(touchingStreets))                 if len(uniqueStreets) <= 2:                       i = 2                     for street in uniqueStreets:                             uRow[i] = street                          i += 1                       uCursor.updateRow(uRow)                      arcpy.Delete_management("inputLayer")      ## **The Problem:** The issue I am running into is the directionality of the street. In my street data I have address directionality using Left From/Left To/Right From/Right To, like in the below table: ![Example Table](http://i.stack.imgur.com/qpsg7.jpg) The Street From and Street To should reflect this same directionality. **What are some ideas for capturing the directionality of the street?** ## **Some ideas I've been exploring:**   * Is there something exposed through ArcObjects that stores the endpoint and beginning point of a polyline using these metrics that I can use to select by location?   * Can I somehow determine the directionality just in ArcPy?   * As a last resort, Network Analyst might be able to offer me this functionality, but due to not having used Network Analyst before the ramp up time may be too much for my purposes.   * If there is functionality that may help me in GRASS/QGIS I'd be willing to explore that as well.   * I also have an "intersections" point file that could be used.


Rendering a QGIS WMS layer in Firefox hogs the machine

I am less than a novice beginner on GIS Systems and am groping for answers. I followed instructions for WIN 7 Install at both http://underdark.wordpress.com/2012/04/06/qgis-server-on-windows7-step-by- step/ and http://hub.qgis.org/projects/quantum- gis/wiki/QGIS_Server_Tutorial#Windows All appears to work well. I can retrieve a capabilities-XML. Next, I add an WMS Layer to a shapefile in QGIS. Again, no issues. Then I attempt to display via Firefox 18.0.2 with address call > > [http://localhost]/qgis/qgis_mapserv.fcgi.exe?SERVICE=WMS&VERSION=1.3.0&SRS=EPSG:4326&REQUEST=GetMap&map=C:/OSGeo4W/apps/qgis/bin/VA_VOTES_BUSINESS.qgs&BBOX=-36.215710,6.817432,-9.700450,58.551008&WIDTH=500&HEIGHT=1000&LAYERS=VA_VOTES_BUSINESS&FORMAT=image/png The entire operating system becomes unstable at the page call. Firefox displays `Connected to LocalHost`; there is 100% CPU usage (Pentium Dual Core 2.2 GHz); and 80% memory usage (4 GB). I have to power off to regain control. I see nothing but errors in the Apache error log but of course I'm not familiar with the references. Spent 10 hours+ trying to troubleshoot and this is my last resort. Any direction would be most helpful and appreciated. QGIS 1.8.0 Lisboa Desktop


Unable to publish geocoding service [ERROR: 001369]

I manage to create a geo locator on arcmap and tested it. but when I publish to arcgis server using arcmap or arcgis server manager. I encounter the error 1) Fail to publish service. The error is very generic on logs. may I know is there anyway to to troubeshoot this issue?


What might slow Inkscape down when working with QGIS SVG or PDF output?

When I import an SVG or PDF exported from QGIS print composer (2.2.0) into Inkscape 0.48 the vector graphic has large black lines and squares. These change size when I zoom in and out, and the graphic is very difficult to manage in Inkscape (slow reaction time when moving around the page and during difficult to group and ungroup vector layers). The map is not large (800 KB) and because it works perfectly other graphic programs (e.g., CorelDraw X5 and LibreOffice Draw 4.1) I believe it is a direct problem with Inkscape. Any suggestions?


Is there a QGIS bezier curve tool?

I'm a recent convert to QGIS but have been confused about the lack of support for drawing curved shapefiles, especially when Arc has almost too many. A lot of the work that our consultancy do is based around natural features such as lakes/rivers so there aren't many straight lines and being able to accurately represent or work with these features requires smooth curves as opposed to the tools i've been able to find on QGIS so far. Other threads i've found relating to this have stated that being able to draw curves would require a fundamental redesign of QGIS, how is this though? Other threads also recommend to use a rough shapfile first then simplify/smooth over to produce your curve but i would be worried about using this simply because you can't see the finished lay of the line when you draw it, which could be frustrating if your boundary moves after smoothing. I really love QGis and dont want to have to go back to Arc, but being able to draw curves is a big sticking point sadly. Any input/idea's would be greatly appreciated. Thanks, Albert


rasterlite build issues

i am able to build rasterlite library but when it is building .exe for tools(rasterlite_load) it gives linking error for external symbole _SPLite3_ like this rasterlite_load.obj : error LNK2019: unresolved external symbol _SPLite3_last_in sert_rowid referenced in function _export_tile there are many like this but all are due to splite I don't know what is missing. Searching google gives sqlite reference. Development Environment - windows and vc libraster version -1.1c Thanks


ArcGIS 10 -- Polygon Centroid Options and Appropriateness (FeatureToPoint, Labels, etc)

I've found a couple of other Stack Exchange questions related to this, but just like the ArcGIS online help, they were not satisfying. For instance, How does ArcGIS calculate the "inside" point for "Feature to Point" Doesn't tell me how the "Inside" point is determined. While the ArcObjects may refer to centroids and label points, I have at least 4 outcomes in my situation that vary differently, wildly in some cases. I have unidentifiable points inside zip codes, so I'm assigning each one to the centroid of that zip code. Then I will calculate straight line distances from these centroids to various locations of interest. Therefore, my choice of centroid is important. Some zip codes are also of odd shapes and may possibly be discontinuous. At least that is what I assume is generating some of the odd behaviors. Now, I have found four different ways to identify possible centrality of these polygons.   1. The polygon geometry contains a centroid information that can be directly extracted (and supposedly is used for labeling?).   2. The labeling location of the polygon (which I found does not always match where the above point is located--hence, the confusion).   3. Feature to point method with 'centroid' option.   4. Feature to point method with 'inside' option. Some of these or most of these may overlap. In some zip codes, they are all completely different. I would have expected 1. and 3. to be similar and maybe even 1. and 2. to be identical, depending on how labeling is set up (I'm using default). My question is not to ask "so how are these calculated?" That's neither here nor there since their algorithms are proprietary, and there's no point in trying to figure it out. The real question is practical: **how should one go about determining the appropriate centroid to use?** My conclusion is the feature to point method with the 'inside' option will be the most appropriate because it is guaranteed to reside within the polygon boundary. I also assume it will be relatively central. But the descriptions about this are nebulous. Hell, consider a crescent shaped polygon. The centroid would be outside the crescent but in the middle of the encapsulating envelope. Would an "inside" option literally just "push" it over into the crescent? We may not know the algorithm, but does the "inside" option actually identify centrality? Is it better than 1.? Maybe I'm missing something in my approach to this problem? Can the projection impact these comparisons? Is there a post-processing review of the polygons I should consider? Note, I'm completely open to using other tools (Python, R, QGIS, etc.). Your answers and advice are appreciated.


PostGIS query vs QGIS filter

I have a PostGIS dataset with around 3 million shapes. I only need to view about 1,500 of these shapes at any given time, but the shapes I need to see depends on which feature is selected in a different dataset. My question is: what is the speed difference between filtering my large dataset within PostGIS and then serving it up to QGIS, vs loading the entire dataset into QGIS and then using the baked-in layer filtering? I'm running PostGIS 2 and QGIS 1.8. I assume PostGIS is faster, but how much faster? I'd like a second opinion before I spend my time coding for one solution over the other. Thanks!


How to integrate CAD and GIS layers in ArcGIS?

I am supposed to integrate CAD files and different GIS layers in a project (as a job interview) in ArcGis. I am writing to know whether there will be any special consideration regarding to this procedure which should be taken into an account that effect the accuracy.


What happened to Live GPS Tracking in 1.8.0-Lisboa?

Dang, I was just getting used to Live GPS Tracking in 1.7.4-Wroclaw. I downloaded the new version 1.8.0-Lisboa on another computer. The Live GPS Tracking thing was in the View pull-down menu in 1.7.4. But, in 1.8.0-Lisboa, it's not there!!! Was Live GPS Tracking taken out of 1.8.0-Lisboa? or maybe Live GPS Tracking is hidden somewhere?


How to rotate the screen in ArcGIS-Online?

I need to print out my map in ArcGISOnline, however I'm using a grid system and need to rotate the screen so that I can print each map individually and have them line up. Any thoughts on how to do this? I will attach a screenshot to help better understand what I mean. ![enter image description here](http://i.stack.imgur.com/o0Hjb.png)


Can ArcPy 10.1 be used to generate true curve elliptical polygons in a file geodatabase?

As background, this question arose from trying to answer a question on **How to generate overlapping polygons from lines output from Table To Ellipse tool?** Using the ArcMap GUI it is very easy to digitize true curve polygons into a file geodatabase feature class using the **Ellipse** Construction Tool but ... Is it possible to write true curve elliptical polygons while reading rows (arcpy.da.SearchCursor) from a table containing a centre point, a major axis, a minor axis and an azimuth for each? I had hoped that there might be a token available with **arcpy.da.InsertCursor** to do this, but `SHAPE@` seems to be limited by what the **Geometry** object supports, and that does not appear to include true curves.


How do I edit a shapefile DBF to repair an extra row?

Something got screwed up and now the rows of my DBF beyond #1210 are offset by one, meaning that each shape has the attributes that should be assigned to the next shape. If I were in a spreadsheet program I would be able to select rows 1210+, copy, then select 1211+ and paste. But I can't do this in the QGIS attribute table editor, and I can't find an editor that will save shapefile DBFs without damage. * * * Apparently I have too little mojo to answer my own question this early. * * * Never mind - I fixed it in Excel. I had previously tried swapping columns in Excel and it didn't work, but I guess rows can be edited. The only minor issue is that the column titles are now all caps. But all the expressions that call them still work, so I guess that's not a problem.


gdalcompare.py not found in GDAL 1.10.1

I cannot find the script gdalcompare.py within the osGeo distribution of GDAL by http://www.osgeo.org/. It is the 1.10.1 though. Is there any way of downloading that script separately from any other source?


Solid Laptop Specs for running ArcGIS

I'm getting a Lenovo thinkpad w520 with 8gb of ram/2000m quadro GPU/2760qm Corei7 to use as my secondary work laptop. I mostly deal with large amounts of ortho photos and some GPS software. Do you guys think this is overkill or any of these specs useless? I'm just hoping this was a good purchase, I didn't think a mid consumer laptop would be able to run GIS very well (without crashing)....given that I do not do much 3D work, just alot of different GIS tasks, is the Quadro 2000m going to benefit me more than say a Nvidia 540m mid level card? Thanks for your help!!!


Alternatives to traditonal map server tile caching

Are there any alternatives to map server tile caching? It seems like there have been many improvements to computers since tile caching has become main stream. With SSD storage and large amounts of memory available on servers, is anyone forgoing all hassle of pre-rending tiles (imagery in particular) and does anyone have any benchmarks comparing dynamic map rendering vs cached?


Reprojection problem with gdalwarp

I have a GeoTiff file with WGS84 datum projection and I'm trying to reproject it using gdalwarp. Using this command               gdalwarp -t_srs '+proj=tpers +azi=19 +tilt=30 +lon_0=-109.5 +lat_0=23 +h=1274199.4 +wktext' atnight.tif atnighttilted.tif      returns this               ERROR 1: Too many points (10000 out of 10000) failed to transform,     unable to compute output bounds.      I'm using GDAL version 1.10.0. This is the contents of the .tfw file:               0.0066667901     0.0000000000     0.0000000000     -0.0066669136     -32.8672753199     72.8560317049      I've just started using gdalwarp. Thanks in advance, Tom The output of `atnight.tif` using `gdalinfo` is:               Driver: GTiff/GeoTIFF     Files: atnight.tif            atnight.tfw     Size is 9123, 4957     Coordinate System is:     GEOGCS["WGS 84",         DATUM["WGS_1984",             SPHEROID["WGS 84",6378137,298.257223563,                 AUTHORITY["EPSG","7030"]],             AUTHORITY["EPSG","6326"]],         PRIMEM["Greenwich",0],         UNIT["degree",0.0174532925199433],         AUTHORITY["EPSG","4326"]]     Origin = (-32.870608714976186,72.859365161672642)     Pixel Size = (0.006666790125743,-0.006666913589395)     Metadata:       AREA_OR_POINT=Area       TIFFTAG_DATETIME=2012:11:26 15:13:57       TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)       TIFFTAG_SOFTWARE=Adobe Photoshop CS6 (Macintosh)       TIFFTAG_XRESOLUTION=72       TIFFTAG_YRESOLUTION=72     Image Structure Metadata:       COMPRESSION=LZW       INTERLEAVE=PIXEL     Corner Coordinates:     Upper Left  ( -32.8706087,  72.8593652) ( 32d52'14.19"W, 72d51'33.71"N)     Lower Left  ( -32.8706087,  39.8114745) ( 32d52'14.19"W, 39d48'41.31"N)     Upper Right (  27.9505176,  72.8593652) ( 27d57' 1.86"E, 72d51'33.71"N)     Lower Right (  27.9505176,  39.8114745) ( 27d57' 1.86"E, 39d48'41.31"N)     Center      (  -2.4600456,  56.3354198) (  2d27'36.16"W, 56d20' 7.51"N)     Band 1 Block=9123x115 Type=Byte, ColorInterp=Red     Band 2 Block=9123x115 Type=Byte, ColorInterp=Green     Band 3 Block=9123x115 Type=Byte, ColorInterp=Blue


Spatial Join Get ZERO Value in ArcMap

The dataset contains about 70,000 point features. And I want to use spatial join tool in arctool box, to join a polygon feature. But it went that the values in the joined column are all zero, values in the 'join_count' are also zero . And the result was right when the point dataset is not that large(actually 40,000 is OK). What's the problem? Do I have any solutions? p.s. I don't want to use the join function in the right-click join. Although it went right, but I want to do spatial join using Geoprocessing.


How to export multiple Shapefiles to one KML file in QGIS?

As the title says, is it possible to export multiple Shapefiles to one single KML in Quantum GIS (Similar to "Map to KML" in ArcMAP)?


How can I reclassify a raster to floating point values?

I have a raster of landuse values that I need to reclassify into pollution concentration values. I have the pollution values in a separate table. They are decimal values (ex. 0.0056), and my landuse values are integers 1 through 6. I have a table mapping the landuse values to the pollution values. I tried using the Reclass By Table tool, setting my pollution values as the "Output value field", but the tool description says that this field must be an integer field. **How can I reclassify the values into decimal values?** I feel like I must be missing a much easier, alternate solution. My goal is to have a pollution concentration raster that I can use later in further calculations.


How to query nearest date in ArcGIS Server REST?

How to query nearest date (today) in ArcGIS Server REST? Table:               ID    NAME            DATA     1     test           2013-06-01 13:00:00     2     test1          2013-06-01 16:00:00     3     test2          2013-07-01 18:00:00      Expected result : query ID = 3


GIS Layer Versioning In GeoServer

i have developed the application for editing the map features online but i want the features to be saved as version of the layer not in original layer.this work is done using GeoExt +GeoServer+OpenLayers.


Is it possible to use commands from ArcMap's main menu in Model Builder

Is it possible to use the functions from ArcMap's main menu (for example Add Data, Start Editing) in my model? I can not find a way to add any other functions except Toolboxes.


Why is there no Map.SnapToLevels property in the ESRI JS API?

The ESRI Silverlight and WPF API's both have a **Map.SnapToLevels** property that can be set on or off. It "gets or sets whether zooming will snap to scale levels if the map contains cached layers". But with the JS API, you always snap to levels and there appears to be no way to turn that behavior off. Is there a way to disable the snap to levels behavior in the JS API? Or can someone provide the rationale behind ESRI's decision to not include that property in the JS API? (I'd like to be able to explain _why_ I cannot provide the same type of zooming in JS API as with WPF API.)


Checking "connectedness'' of line shapefile in ArcMap

I inherited a line shapefile in ArcMap that is used to represent a road network. The problem is that the road network is pretty large and it is impossible for me to: a. Tell if the end of each line (road segments) is "connected" to either an edge, vertex, or the end of another line (road segment); and b. Tell if the entire road network is "connected" - i.e., can I trace a route from any point on the road network to any other point on the road network? Is there a way to accomplish this in ArcMap? That is, is there a function that will allow me to highlight line segments that do no connect to other line segments, or something similar that would allow me to answer the two questions above?


How to connect lines in a network?

I have a shapefile of lines and another of polygons. When two or more lines intersect a polygon they should be connected (snapped), but sometimes during manual digitization this doesn't occur. I need to find and correct all these errors. I used `v.clean` from GRASS GIS and it automatically snapped a lot of vertices, which was good, but some of those disconnected lines are further from the tolerance I gave and if I increase the tolerance too much it begins to generate other types of errors. ![enter image description here](http://i.stack.imgur.com/cUNFp.png) Any ideas?


What arguments does MapBasic's ColumnInfo() take?

MapBasic's `tableInfo()` function has a KnowledgeBase docs page here listing the arguments and all the attributes it takes, and a usage example in the reference document. `columnInfo()` has neither of these (that I can find). I can't find any docs for it anywhere, other than a passing mention in the reference docs confirming that it exists: > The functions **TableInfo( )** , **ColumnInfo( )** and **NumTables( )** let > you determine information... > >   * **ColumnInfo( )** returns information about a column in a table, such as > the column’s name, the column’s data type, and whether the column is > indexed. > That's it. What arguments does it take (e.g. is it `columnInfo( nTabNum, nColNum, Attr )`? Or maybe `columnInfo( objColumn, Attr )`?) and what are the attribute names (e.g. are they named COLUMN_NAME, COLUMN_INFO_NAME, TAB_INFO_COL_NAME... etc?)


How can I use a sub-query in where statement of a query in ArcGIS 10.2 Map service?

I'm migrating my **ArcGIS** map services from version 10.0 to 10.2. It uses ArcSDE to manage spatial data in a **SQL Server** database. I was able to create the map service and everything seems OK. However if I add a sub-query in the WHERE statement of a query form in any map service, I receive "Unable to complete operation." error as the following image: ![enter image description here](http://i.stack.imgur.com/qb6vx.png) The same query was working in ArcGIS server 10.0. Also I can run the same query to select data, when I use the map document directly from ArcMap 10.2. The ArcGIS Server 10.2, ArcSDE, and SQL Server are installed on a server which is internally accessible from the web server. The ArcGIS Web Adaptor is installed on the web server. Both servers have windows Server 2008 and IIS. I appreciate your advice. Thank you


ArcObjects - cheapest way to check if license has been checked out

So I noticed a big performance hit when I add the following line to my code:               ESRI.ArcGIS.RuntimeManager.BindLicense(ESRI.ArcGIS.ProductCode.EngineOrDesktop);      I have two separate classes that call the BindLicense method. After the first time it is called, it isn't necessary to call again, but neither class can assume that has happened. What I want to know is if there is a cheaper way of checking if a license exists before attempting to bind? I feel like I've seen code to do this in Python, but don't know if there is the equivalent in ArcObjects...


QGIS API - Print Composer Template produces empty map?

I have a **print composer template** which I'm trying to populate with different subset queries from a CSV-derived vector layer using the **Python** API. However, all the output PDFs contain everything **except the layers on the map**! I'm doing this with a **standalone app** , and it functions fine in that it does produce PDFs and the data is all loaded and interrogated just fine (I can subset query, count features etc.) Here's the method I'm using to render the composition, as you can see I've tried to **set the layers** on the map in different places (actually there are 2 maps - one shows the extent rectangle of the other, but I can't get either to show anything) - this is my first QGIS scripting experience so I'm sure I'm missing something very obvious...               def compose(layers, species_binomial):         # vlayer is the CSV vector layer, outline is the ADM0 country vector.         vlayer, outline = layers              # construct a QStringList for the renderer setLayerSet         layers = QStringList(vlayer.id())         layers.append(outline.id())              # 1 maybe set the registry layers???         QgsMapLayerRegistry.instance().addMapLayer(vlayer)         QgsMapLayerRegistry.instance().addMapLayer(outline)              # 2 maybe set the canvas layer???         canvas = QgsMapCanvas()         canvas.setExtent(outline.extent())         canvas.setLayerSet([QgsMapCanvasLayer(outline), QgsMapCanvasLayer(vlayer)])              # 3 set the renderer layers???         renderer = canvas.mapRenderer()         renderer.setLayerSet(layers)         renderer.setExtent(outline.extent())              composition = QgsComposition(renderer)         with open(os.path.join(ROOT, 'species_composer.qpt'), 'rt') as fp:             composer_template = fp.read()         doc = QDomDocument()         doc.setContent(composer_template)         composition.loadFromTemplate(doc)              # 4 maybe set the map layers of one of the maps???         map1 = composition.composerMapItems()[1]         map1.setLayerSet(layers)              title = composition.getComposerItemById('title')         title.setText(" ".join(species_binomial))  # this works fiiiine              out_path = os.path.join(ROOT, 'out',  "%s_%s.pdf" % species_binomial)              # make a lovely blank PDF with a title and a blue map rect         ...      What am I doing wrong? Here's an example output... ![PDF output from print composer template](http://i.stack.imgur.com/5dSHY.png)


Convert GML to DXF using ogr2ogr

I cannot make ogr2ogr.exe work. I'm using GDAL 1.11. I keep getting the following error message: > ERROR 1: Did not get at least 2 values or invalid number of set of > coordinates 709993.98 6181946.65 32.36 709990.65 6181946.07 32.36 709992.47 > 6181935.8 32.36 709995.78 6181936.38 32.36 709993.98 6181946.65 32.36 ERROR > 1: Invalid exterior ring I run the following command: > C:\Python34\Lib\site-packages\osgeo\ogr2ogr.exe -skipfailures -nlt > MULTILINESTRING25D -dim 3 -f "DXF" -a_srs "EPSG:25832" -spat 711000 6168300 > 713000 6168600 -geomfield posList bygning.DXF Bygning.gml I'm working on a relatively big file 4.2GB. **UPDATE: Example of data to import**               <gml:surfaceProperty>         <gml:Surface srsName="EPSG:25832" srsDimension="3">             <gml:patches>                 <gml:PolygonPatch>                     <gml:exterior>                         <gml:LinearRing>                             <gml:posList>                                 710820.22 6180637.52 36 710818.14 6180633.83 36 710825.47 6180629.71 36 710827.55 6180633.4 38.25 710828.6 6180637.3 38.25 710825.18 6180638.36 38.25 710824.22 6180635.27 36 710820.22 6180637.52 36                             </gml:posList>                         </gml:LinearRing>                     </gml:exterior>                 </gml:PolygonPatch>             </gml:patches>         </gml:Surface>     </gml:surfaceProperty>      Does the GML3 driver in org2ogr support surfaces? **UPDATE 2 - making it simple** I tried making a gml consisting of only one feature, same namespace etc. as before. I then output the feature to WKT. Output WKT "POLYGON ((711383.39 6182528.23,36.49 711379.31,6182528.35 36.48,711379.06 6182520.08,36.48 711383.15,6182519.97 36.48,711383.37 6182527.56,36.49 711383.39,6182528.23 36.49))" Input GML http://pastebin.com/76AGLWH0 Cleary the commas is placed in the wrong positions. The z-coordinate of the first point is parsed as the x-coordinate of the second feature. If I add -dim 3 or -nlt GEOMETRY, it adds zeroes. BR Thomas


How to create a time slider using arcgis javascript that only shows a year at a time?

I am looking at this example on how to create a time slider using the arcgis javascript api. Instead of showing the data cumulative from the start date. I want to just show the data for that particular year. I have tried changing the ThumbCount to 1, but all that does is set the start time to the start time. Other than that, I can't seem to find any more information on how to get the behavior that I want. I understand that I can use 2 thumbs on the time slider in order to get desired result. However, I am looking to see if it is possible by using only 1 thumb slider.


Installation issue with OpenGeo-Suite

Ok so last time had a big issue with running WFS. Seems the issue was neither with my OpenLayer code nor Geoserver. The issue was with the servers not able to communicate with each other. The tomcat 6 server on which geoserver was running and apache server which is the default web server for my ubuntu. After some digging around it seems that I need to have tomcat 6 pre-installed on my machine before installing OpenGeo Suite. So having reinstalled the package, at end, the terminal gives me the following message               NOTICE: Please run /usr/share/opengeo-suite/tomcat-admin-setup.sh to complete this      installation.      So i executed the following command (from the required directory),               sh tomcat-admin-setup.sh      And I am getting this error,               tomcat-admin-setup.sh: 143: [: Ubuntu: unexpected operator     Unknown distro. The OpenGeo Suite is supported on Red Hat Enterprise Linux 5 and 6,      CentOS 5 and 6, and Ubuntu Linux 10.04 or later.      Does this mean the OpenGeo Suite is not compatible with Ubuntu 12.04??


ESRI Javascript API - Emulate 'Red Glow' Effect From Flex API

I am trying to emulate the red glow effect so common in the Flex API Examples. Really, any type of blinking/flashing style marker will do. I am assuming I will have to use a picture marker symbol with an animated gif? I have probably infected my computer with who knows what downloading trial software to convert from swf to gif and all types of other non nonsensical conversions. I am hoping someone can point me to an example or give me some insight on how to accomplish this.


How to repeatedly mosaic raster in Erdas?

I have found script to mosaic raster in R. But the output raster size is almost 9 times bigger (92mb) than mosaic it in Erdas (12mb). So I prefer mosaic raster in Erdas. But it wastes so much time to mosaic long term (time series) satellite data one by one. So I am looking for possibility to repeatedly mosaic raster to create multiple outputs automatically by single attempt as like batch in Erdas.


Opening PDF generated by Python AddIn using Report (*.rlf) file automatcally?

As a test I have:   1. Authored a test.mxd with one layer called Localities in a single data frame called Layers   2. Authored a test.rlf using View | Reports > Create Report of ArcGIS 10.2.2 for Desktop   3. Created a Python AddIn which creates the expected PDF file from the Report (*.rlf). This is the code which all works to create the PDF report from the Localities features within the current extent:               import arcpy      import pythonaddins          class ReportForExtent(object):         def __init__(self):             self.enabled = True             self.checked = False              def onClick(self):             mxd = arcpy.mapping.MapDocument("CURRENT")             df = arcpy.mapping.ListDataFrames(mxd,"Layers")[0]             lyr = arcpy.mapping.ListLayers(mxd,"Localities",df)[0]             arcpy.mapping.ExportReport(lyr,r"C:\temp\test.rlf",r"C:\temp\test.pdf","EXTENT",extent=df.extent)             del mxd      Is there a way to have the Python AddIn not just create the PDF file but also open it in Adobe Acrobat Reader?


How to install QGIS 2.2 instead of 2.4 to get a working OpenLayers plugin again?

I use ubuntu . I don't know how to return installed qgis 2.2 I try uninstall qgis 2.4 and try find how to install old version qgis 2.2 Please .. tell me how to install old version qgis 2.2 for ubunut 12.4


Examples of Sanborn map WMS?

I'm wondering if anyone has heard of or come across any examples of scanned and georeferenced Sanborn Fire Insurance maps being served as a WMS. I know there are copyright issues, so presumably only for maps before 1922. I'm curious about the possibility of doing this myself... Thanks!


Trouble using ogr2ogr w/ -s_srs to combine inconsistent EPSG files into my postgis DB

I'm having trouble combining shapefiles from TIGER and different local governments into a single postgis db to query legislative locations. I've imported all data to multipolygon geom, but found some datasets' polygon shapes returning really large floats values lat/lng. Looking more closely, I noticed inconsistent EPSG coordinate systems, based on searching prj2epsg.org with contents of my .prj files:               #Census Data (EPSG:4269 - GCS_North_American_1983)         GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137,298.257222101]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]]          #Texas Data (EPSG:3081 - NAD_1983_Lambert_Conformal_Conic)     PROJCS["NAD_1983_Lambert_Conformal_Conic",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Lambert_Conformal_Conic"],PARAMETER["False_Easting",1000000.0],PARAMETER["False_Northing",1000000.0],PARAMETER["Central_Meridian",-100.0],PARAMETER["Standard_Parallel_1",34.91666666666666],PARAMETER["Standard_Parallel_2",27.41666666666667],PARAMETER["Latitude_Of_Origin",31.16666666666667],UNIT["Meter",1.0]]      I can't seem to get ogr2ogr to explicitly transform the target coordinate systems with `-a_srs EPSG:3081 -s_srs EPSG:3081 -t_srs EPSG:4269`:               ogr2ogr -f PostgreSQL PG:"host=..." file.shp -nlt MULTIPOLYGON -s_srs EPSG:3081 -a_srs EPSG:3081 -s_srs EPSG:3081 -t_srs EPSG:4269 -lco GEOMETRY_NAME=the_geom -nln table_name.      I've tried a number of different variations on the above, but at best I get `ERROR 1: AddGeometryColumn failed for layer table_name, layer creation has failed`. **What am I missing here?** I'm still a bit of a novice on work with shp files. Any help or suggestions would be really appreciated! Example source files:   * From Texas: ftp://ftpgis1.tlc.state.tx.us/DistrictViewer/House/PLANH309.zip   * From Tiger: ftp://ftp2.census.gov/geo/tiger/TIGER2013/CD/tl_2013_us_cd113.zip


Connecting Grass with R

I need help to connect R with Grass. In RStudio when I send the line: loc <\- initGRASS("C:/Program Files/R/R-3.1.0/bin/i386", home=tempdir(), override=TRUE) It says: Error in system(paste("g.dirseps.exe -g", shQuote(Sys.getenv("WINGISRC"))), : 'g.dirseps.exe' not found Guys, Can you help me?


How to translate .bil files in a batch mode with QGIS?

I have some .bil files that I would like to convert into ASCII Gridded XYZ format files. It seems to work well with the Raster>Translate tool of QGIS when I process one file at a time. But I would like to do it by batch. So I checked the batch box but got this error message: Processing of the following files ended with error: > > C:/Users/Fall Line Tower/Desktop/Bil To Be Processed/Only > Bill/PRISM_tmax_stable_4kmD1_19810101_bil.bil ERROR 4: `C:/Users/Fall Line > Tower/Desktop/Bil To Be Processed/Onli > Bill/PRISM_tmax_stable_4kmD1_19810101_bil.bil' not recognised as a supported > file format. GDALOpen failed - 4 I have the same error message for each .bil file within my source directory. Any idea where the problem comes from? Does the batch mode take only one type of format? Not the .bil? I am using QGIS 2.2 on Windows 7. Thank you


Managing large number of catchments/watersheds

Does anyone know of an easy way to manage a large number (10,000's) of catchments/watersheds? I have defined a stream network and for each segment, I have a corresponding catchment (in a unique shapefile). The name of the catchment is the x and y co-ordinate of the catchment outlet which is also assigned to the attribute table of relevant stream segment. The end user will only want to view the stream network and maybe a few of the catchments (certainly not all!) at a time. I need a way to make it easy for the end user to view the catchments of interest (I guess they could search for them). My only idea (which I am a little unsure of how to implement) is to build a tool (in ArcGIS) that loads the catchment when they click on a stream segment (with the tool selected) If anyone has a better idea... For my idea, I think I would have a button similar to the identify tool the end user would use this tool to select a stream segment, the tool would then look in the attribute table to find the file location of the catchment and add this as a shape file. Is such a thing possible to do? I have coded up some python scripts before but have never attempted to integrate such functionality into ArcGIS. Any thoughts or direction would be appreciated. We only have the basic version of ArcGIS, although if this sort of functionality is available in something like ArcHydro I'd be interested to know!


What are the advantages/disadvantages of populating a PostGIS database with Raster information?

...as opposed to keeping the rasters as locally stored GEOTIFFs? I have a large project that will require making available approximetely 50,000, 10mb GEOTIFFs as WMS to an OpenLayers based web map. Only 1 WMS will be displayed at any given time so performance is a non-issue. Considering MapServer and GeoServer for this application. I would like to know what are the advantages and disadvantages of storing this type of data. Will updating a PostGIS database when the raster data changes be more difficult, rather than replacing a single (or series of) raster files. WMS load times drastically improved with PostGIS?


Looping multiple feature classes in python for ArcGIS 10

I currently have a large list of feature datasets and dbf tables in a geodatabase. I would like to join each shapefile with its matching DBF file using the "Make Query Table" tool and then "Copy Features" to produce a new shapefile with specified fields. I am trying to automate this process so I don't have to manually select field names for each set of data (output table would use the same fields). By using the Iterate Datasets in ModelBuilder, I was hoping to do this. However, the Make Query Table will not connect to the geodatabase and still requires me to plug in a specific table. Does anyone know I can use looping to run this process? Are there any similar examples using python? Example of Geodatabase: AAA1_part1.shp AAA1_part1.dbf AAA2_part2.shp AAA2_part2.shp ..... ..... Any information would be greatly appreciated! Thanks


Calculating Field based on name of shapefile

I have multiple shapefiles to which I have added a field called "gname1". For each shapefile, I want "gname1" to be set to the name of that particular shapefile, for every feature within the shapefile. I realize there are multiple ways to do this, but I want to know how to do this within the "Field Calculator" tool. So, in the "gname1=" box under Field Calculator, what is the expression I use to denote that I want the field to be set to the shapefile name? I am very new to python, and not entirely familiar with the procedure for using python scripts to calculate fields through this tool in Arc10. I have seen other questions similar to this one, but I was confused on how to do this _within_ the field calculator tool. Thanks for any help!


Get the different values assigned to a zip code

I have a map to which I added a zip code layer and another layer which contains values assigned to every square kilometer within each zip code. I'm trying to find out how many different values there are in each zip code and what they are. So for example, if a zip code has 10 values assigned to it in total, 5 of them being 1.5, 3 of them being 2.1 and the remaining two being 6.4, this zip code basically has **3 different values** , namely **1.5, 2.1 and 6.4**. That is what I want to know. Is there a way to do this without manually checking all the values in every zip code? * * * This is a screenshot of my data. The value that I'm working with is under the column **BI**. I need a table which gives me a list of the long/lat in every zip code and the corresponding values. I can't do anything on excel yet since I don't know which lat/long goes with which zip code. ![data spreadsheet](http://i.stack.imgur.com/wTcJs.png)


ArcGIS Javascript API - Refresh

The ArcGIS Javascript API has a nice refresh interface for ArcGISDynamicMapServiceLayer:                   var dynamicMapServiceLayer = new esri.layers.ArcGISDynamicMapServiceLayer("http://sampleserver1.arcgisonline.com/ArcGIS/rest/services/Demographics/ESRI_Census_USA/MapServer");         map.addLayer(dynamicMapServiceLayer);         dynamicMapServiceLayer.setDisableClientCaching(true);              if (dynamicMapServiceLayer.disableClientCaching) {             aspect.after(dynamicMapServiceLayer, "onUpdateEnd", function () {                 $.delay(15000, "UpdateDynamicLayer", function () {                     dynamicMapServiceLayer.refresh();                 });             });         }      But this is not generalized across the other layer types. Anyone have a solution for refreshing layers of any type?


Is there a simple way to create a new custom vector marker symbol for ArcMap?

I would like to create a set of custom marker symbols for environmental sampling locations. I need there to be "vector" marker symbols not picture symbols. I've tried picture marker symbols and there are numerous drawing and exporting issues associated with this approach especially when using transparent colour. I believe I used a font creator in the past and drew a custom true type font but there seem to have been some issues with accessing the custom font from ArcMap as well as some font size problems. Ideally, I would like to draw a vector graphic in Illustrator/ CorelDraw then somehow use this graphic as a marker symbol in ArcMap. Can this be done? I would really appreciate if someone could **describe the process of creating a custom "vector" marker symbol** ; tools used to create the symbol as well as accessing the symbol in ArcMap. Vector graphic approach is preferred but a true type font approach is also fine.   Thanks!


leica zeno 5 gps handheld

I recently purchased the Zeno 5 gps handheld device form leica geosystems. I would like to ask if anyone knows what other software I can use except the one they offer that integrates with ArcGIS desktop. They suggested the 'Zeno office' and the 'Zeno field' or this software http://www.microsurvey.com/products/fieldgenius/price2007-canada.htm and the 'Zeno connect'. It wasn't clear on the website that you have to buy this extra software so that the zeno 5 gps device will function properly. From my experience I know that when you buy a gps handheld device you dont need to buy extra software to display your position, the satellites, your accuracy, your lat and lon in a user interface. Is there a less expensive alternative software I can purchase and install in this device? Thanks in advance!


Print preview showing distortion in layers from GeoServer via MapFish and OpenLayers?

I have an application which shows the print preview before printing the map and the problem is that the layer is distorted in the first image but when I zoom into the map the distortion goes away. I am using MapFish as the printing module and I assume that OpenLayers deals with the print preview part. I am new to this so I am not sure if this a code based problem or a setting based problem (either in GeoServer or MapFish) ![enter image description here](http://i.stack.imgur.com/0AIRS.png)


Split polygon with holes to get several hole-less polygons

Unlike this thread I actually need to split the polygons into several hole- less polygons. I'm not sure on how to approach to this problem. I have polygons with holes of many diferent shapes, like buildings, esplanades, rivers, etc. And I need to split them in a coherent way (if the polygon is a river it should be splitted perpendiculary to his axis) and with less splits as possible. I need to do it programatically. I'm using ArcObjects 9.3, but if there is no other way I could try to use a specific algorithm to solve this problem.


Calculate MBR envelope for different geometry shapes in different rows

I have some different kinds of geometry (Polygons and line strings) stored in my mysql database. On different rows. I know I can retrieve the envelope of one row with: `SELECT Envelope(zoneShape) FROM geometryTable WHERE refid=4` If there are three different shapes associated with reference id 4, i'm getting three results, obviously. How can I (preferably with one sql query) get the envelope/MBR of all the geometry with e.g. refid=4?


How to overwrite a shapefile or feature class using arcobjects?

I can copy a shapefile but I want to be able to overwrite the copied shapefile to prevent multiple copies being created each time I run the program. I would also like for the saving and renaming of the copied shapefile to happen behind the scenes i.e. the user doesnt see a gxDialog box. Any suggestions on how to overwrite the copied file, either as a shapefile or feature class would be greatly appreciated.


Use output of Add Field as input for Calculate Field in ModelBuilder

I have created a simple model in ArcGIS 10 which, after some initial processing of a feature class, outputs `outputfc` and adds a field to it. Regardless of the feature class input (which the user specifies), the field `newField` that is added is always the same. I then need to do some calculations in this newly added field, but because I know this field will only have the name that I specified, I select `outputfc` as the input table for Calculate Field, and an input variable field with the name of the field `newField` as the value. This works fine to achieve the outcome I want, but if I allowed the user to specify the name of the added field, how would I connect this to the next Calculate Field process?


Possible ways to display maps offline, inside mobile hybrid app

I want to develop a mobile hybrid application to run on both Android & iPhone devices in which i want to display my own map data. I have the SHP files with the data i need. I have Geoserver and PostGIS for loading and serving the map data online. So far everything is fine. My problem is that i still havent figure out the correct walkthrough for displaying the map data OFFLINE. Although i have read many answers here in this website, none of them seems to solve (as per my understanding my issue). The goal is to display points of interrest on map. offline. Some Points of interrest doesnt have address.Just Lon/Lat. I want to cover these requirements:   1. How to offer the map data & the content offline to users device   2. How to offer to user routing options to navigate him to the destination from current location   3. How to add content data from a WordPress CMS, inside the app? (via JSON?GeoJSON?GeoRSS?)  Have been searching and looking many Tile services and many approaches to achieve the requirements, but still havent figure out the correct combination. For the development of app I've been using PhoneGap XDK, Eclipse ADT. And for displaying the maps on mobile devices I'm using OpenLayers mobile. Looking forward your suggestions or walkthroughs. Thank you in advance, with respect


QGis not upgrading from 1.9.0. Ubuntu, Raring, using nghtly PPA

I'm currently running 1.9.0 on ubuntu raring, through the nightly ppas:               deb     http://qgis.org/debian-nightly raring main     deb-src http://qgis.org/debian-nightly raring main      But it doesn't update with software updater, even though I know there are more recent versions it should be pulling down. Can anyone think off-hand why that might be or what I could do to fix it? I noticed that on the new qgis homepage, the nightly PPAs are not listed for raring... have they been discontinued for some reason, or temporarily broken? Thanks Simon


How to set min and max value for stretchType: esriRasterStretch_PercentMinimumMaximum

How can I set the min and max values for the stretch type "Percent Clip" in arcobjects? It's possible to set the Stretch type:               rasterStretch.StretchType = esriRasterStretchTypesEnum.esriRasterStretch_PercentMinimumMaximum;      But I have no idea how I can set the min and max values (like in the Properties Dialog of ArcMap). ![Layer Properties with Stretch type Percent Clip and min/max percent values](http://i.stack.imgur.com/T9GLr.png)


Is there any "proper" way of setting the viewport for QGis Time Manager?

I'm working on video visualisation for some geo data. Everything works fine for me except the viewport that the time manager plugin is using for image frames. It makes image frames size of actual map window size which is something like 1248x636 in my case now. When I resize the QGIS window the time manager plugin starts exporting image frames of different size. Is there any "proper" way of setting the viewport so I don't have to play around with application window size? Or should I perhaps bulk edit (crop?) all the frame images before merging them with mencoder? I'm using Quantum GIS 1.8.0 with Time Manager 0.7