title: Why does the 3D Analyst 2D export renderer create anomalies?

text: A lot have been improved in the 3D Analyst 10 including the ability to export large scenes as large 2D images. However, some problems still exist one of which is: When I export a scene that contains a very large TIN the resulting image appears to have some triangles shaded darker which causes a "spotty" appearance. It's almost as if some "faces - triangles" were shaded incorrectly, or reversed?   My current graphics card is NVIDIA Quadro 4000. (Used to have the NVIDIA Quadro FX 1700 which caused similar issues with the 9.3 version of 3D Analyst) I should also mention that it does not matter what type of image i export, compressed or not the result it more less the same... Does anyone have any ideas as to how to rectify this? EDIT/ UPDATE: It seems that this is only happening when the TIN is exaggerated, in this case 5x EDIT/UPDATE 2: If i drape imagery over the terrain and turn the IN off, there are no issues. This problem pertains only to the display of TIN Arc Scene screen capture: ![alt text](http://i.stack.imgur.com/o9G1Q.jpg) Exported Image (jpeg) ![alt text](http://i.stack.imgur.com/97ddg.jpg)


title: Iterate through TOC elements in ArcGis Server 9.3

text: I'm developing website with a custom tool in ArcGis Server 9.3 using Visual Studio 2008 and the ADF components, the thing is that when I click on my custom tool it calls a class, this class needs to check each layer available to see if it is enabled and visible (this properties changes when you zoom out or zoom in), does anyone knows how to accomplish this? Thanks.


title: FME - Is it possible to add custom command line input arguments to a workspace?

text: I have a workspace whose command line call is like so:               fme.exe DwgImporter.fmw               --SourceDataset_ACAD Z:\FME\Data\DWGs\*.dwg               --DestDataset_MSSQL_ADO IHS_TPW               --SourceDataset_GEODATABASE_SDE sde      Is it possible to somehow add my own custom input arguments that can then be read from within the workflow? Something like:               fme.exe DwgImporter.fmw               --SourceDataset_ACAD Z:\FME\Data\DWGs\*.dwg               --DestDataset_MSSQL_ADO IHS_TPW               --SourceDataset_GEODATABASE_SDE sde               --CustomArg1 arg1               --CustomArg2 arg2      Or does anyone know of another way to pass in variables at runtime?


title: Can I incorporate an existing DBF file with coordinates for Tiff tiles into geodatabase like spatialite?

text: I have 285 Tiff files that make up the 1:10000 OS map of Norhtern Ireland. I am having problems running QGIS when I try to load all 285 tiff files in for base maps (I think the memory on computer is problem here) My thinking was if I create a geodatabase this will speed up accesssing these base maps. I also have with each tiff file a corresponding AUX , RRD, TFW file. There is a historical DBF file that appears to index each Tiff file with max and min x.y coordinates. I beleive this was used for ARC 3?? Is there any easy way to get this DBF file into a spatial database to allow me to access the Tiff files quicker? I tried pyrimids on the Tiffs but they are not 'pyrimid-able'. I also tried to create a geodatabase using spatialite but becasue the Tiffs are not geo-referenced (this info sits outside the Tiff file on the TFW file) it would not allow me to load them into my DB. Any assistance would be greatly appreciated


title: Show loading icon with custom tool

text: I have a custom tool in vb.net for a web application, i'm using Asp.net, the class is in vb.net and the gis service is in Arcgis Server 9.3 sp1, the tool makes a selection of layers and takes a lot of time, so i want to show a loading icon on the map while it's getting the info from the server. Any idea? I already tryed to modify the progress bar from the map but it works only within the map events, i want to show the progress bar from the custom class, from vb.net Thx


title: Strange line appearing

text: I have this line in my project's layer which i cant make disappear: ![strange line](http://i.stack.imgur.com/c5lRZ.png) I hope you can see it. In QGIS is much more visible. An it also appears on the print composer. Any suggestions would be very usefull. Thanks in advance. Tiago Brandão.


title: How do I add (and map) a non-spatial table to QGIS from SQLite?

text: I have an sqlite file from my Wigle Android app (wiglewifi.sqlite). It's not a spatial file, just tables with lat/lon and other data. I can add spatialite layers just fine, but cannot seem to find a way of adding this table in the same way I'd add a delimited text layer - how to go about doing this? Clarification: It's possible to add as a table (as per @Giovanni-Manghi's answer below), OR it's possible to export the table as a CSV and then add as a delimited text layer, but ideally I'd like a **point layer based directly on the SQLite table** (in the same way as in ArcGIS I would add a table then right-click and "Display XY data...")


title: How to join an external table to a set of polygons, and symbolize on the joined values, in a web map?

text: Say I have a table showing a number of statistics - for example, various gun crime by US State measures. I also have a dataset of state boundaries. How can I join the table to the features in a web map, and symbolize the polygons using attributes from the table? The joined table must be specified within the web application (rather than hard-coded via ArcMap) since the statistics will be derived on-the-fly and will be specific to each user. I need the ability to change the symbology on-the-fly, for example by choosing a different crime statistic to map, in which case the map should update immediately. **ArcGIS Server** Using ArcGIS Server 10.1 the REST API supports table joins. Are there any sample scripts which step through how to set this up? **Other environments?** Can the same thing be achieved using FOSS software? Thanks. (This is a modified version of my earlier question since that hasn't yet received any responses)


title: How to change the QGIS help dialog?

text: I am working with python plugins for QGIS. I am done with my one of my .py file using python 2.5 and pyqt4. It is just help file for QGIS, which i changed. I wanted to incorporate the same in QGIS. I installed QGIS 1.7.1. I know it to be done through python plugins. But how to go about it? where to copy my file, so that it will get open from QGIS? QGIS already has that page (help file ), I need to change it. Help me out.


title: How to create a efficent linestring out of geography points using postgis

text: I'm trying to create a linestring out of geography points. At the moment I'm doing it like this:               select ST_GeographyFromText(ST_Makeline(a.record_point)) from (     select (record_point::geometry)     from records order by record_timestamp) a        Problem is that most points have unnecessary decimal places (6 is sufficient) and some points are exactly the same. I can set the decimal place precision when converting to KML but is there another way to do this and how to get rid of overlapping points? **My own solution:** It is PostgreSQL specific because of the window function. Basically I just leave out spatially equal sequential geometry points and reduce the precision with geojson converter (kml converter has the same option).               select st_asgeojson(ST_MakeLine(points.record_point),6) from (     select (record_point)::geometry as record_point     lag( (record_point)::geometry, 1) over w as prev_point,     from records      window w as (order by record_id)     order by record_timestamp     ) points where not st_equals(record_point, prev_point)


title: IGPutilities in IGPFunctions

text: Quick question, in Esri's area calculation sample of IGPFunction implementation they use IGPUtilities but they don't release the instance. On the singletons page in SDK help they list GPUtilities as a singleton object. So shouldn't our instances of IGPUtilities be released? And also should we not use the Activator class to create the instance? In a IGPFunction implemented class do people have any preference on creating a class instance of IGPUtilities (like in the Esri Area Calculator example) or create them in the methods when needed? Thanks Rob


title: What is the longitude and latitude of Tokyo, Japan

text: Please help me I need to know for a project


title: "Missing dependency.This algorithm cannot be run :-(" QGIS + GRASS Windows 7

text: Windows 7, OsGeo4W pack installed including QGis 2.4.0 and GRASS. I Qgis when i try to execute a GRASS geoprocessing algorithm i get this popup: "Missing dependency.This algorithm cannot be run :-( This algorithm requires GRASS GIS 7 to be run. Unfortunately, it seems that GRASS GIS 7 is not installed in your system, or it is not correctly configured to be used from QGIS Click here to know more about how to install and configure GRASS GIS 7 to be used with QGIS" The link doesnt help me. I get the same message at my computer at work with Qgis 2.2.0 The GRASS plugins are activated under geoprocessing settings Need help! Thank you Do i need to configure something?: ![enter image description here](http://i.stack.imgur.com/ccyHF.jpg) When i run MSYS i get this: ![enter image description here](http://i.stack.imgur.com/6KtPQ.jpg)


title: Alternative solution to ESRI's Geometric network for wet Infrastructure

text: I'll be glad for some help with "market research" on GIS for water pipes networks   what are the differences between these solutions in compere to ESRIs geometric network GHydraulics WaterCAD V8i WaterGEMS V8i are they based on graph and connectivity table ?, do they allows flow simulation and finding the nearest valve (isolation trace) ? is there any more solutions ?


title: How to make a new raster out of two other rasters

text: I have two raster images, one is the Corine Land Cover for the year 2000 and another one which is the anthropogenic effect on land. The first image reclassified in 9 categories:   * Artificial surfaces   * Agricultural areas   * Forests   * Scrub and/or herbaceous vegetation associations   * Open spaces with little or no vegetation   * Inland wetlands   * Maritime wetlands   * Inland waters   * Marine waters And the second image has seven categories:   * Closed τree or shrub vegetation, canopy cover > 70%   * Open tree or shrub vegetation, canopy cover > 40-70%   * Scattered tree or shrub vegetation, canopy cover <40%   * Bare of any vegetation area   * Αbandoned fields, grasslands and scrubland zones of deciduous oaks and fir   * Βrushwood   * Cultivated areas What i want to do is to built a new raster with the information from these two rasters. For example i want to keep the information from "waters" from corine, and combine the information for Forests from corine and from the other raster the "Closed tree or shrub vegetation". I need you help.


title: TileCache and Openlayers Problem - Projections?

text: I have a very known problem, but I can't solve it. I have set the same maxExtents in the config as in the Openlayers init JS. I have done the same thing with the projection (EPSG:900913) and done all the transformation on OpenLayers FW. But there is still an error, that the tiles can't be shown because they don't start at the right point. It seems that they are not constraint. Why? Any idea? Where is the keypoint in the solution! Thanks. Martin Edit: * * * I don't see any pics as tiles. I just see red-loadfail boxes. If I click with the right mouse button on the tile and let me show the picture I get this message: "An error occurred: Current x value 967644.673573 is too far from tile corner x 939258.202600" * * * My tilecache.cfg config is: * * * [osm] type=MapnikLayer mapfile= /var/www/tilecache/osm.xml #layers=osm spherical_mercator=true #tms_type=google srs=EPSG:900913 bbox=654558.60577333,5745163.8837113,1173307.4327978,6075303.6111284 #bbox=5.88,45.78,47.81,10.54 #maxReolution=19567.879237500 maxResolution=156543.03390000001 resolutions=156543.03390000001, 78271.516950000005, 39135.758475000002, 19567.879237500001, 9783.9396187500006, 4891.9698093750003, 2445.9849046875001, 1222.9924523437501, 611.49622617187504, 305.74811308593752, 152.87405654296876, 76.43702827148438, 38.21851413574219, 19.109257067871095, 9.5546285339355475, 4.7773142669677737, 2.3886571334838869, 1.1943285667419434, 0.59716428337097172, 0.29858214168548586 maxResolution=78271.51695 extent_type=loose levels=20 * * * OSM CODE * * *                 function init() {         var tilecacheBbox = new OpenLayers.Bounds(5.88, 45.78, 10.54, 47.81).transform(             new OpenLayers.Projection("EPSG:4326"), // transform from WGS 1984             new OpenLayers.Projection("EPSG:900913") // to Spherical Mercator Projection           );         alert(tilecacheBbox);              var map_options = {             maxExtent: tilecacheBbox,             restrictedExtent: tilecacheBbox,             maxResolution: 156543.03390000001,             resolutions: [156543.03390000001, 78271.516950000005, 39135.758475000002, 19567.879237500001, 9783.9396187500006, 4891.9698093750003, 2445.9849046875001, 1222.9924523437501, 611.49622617187504, 305.74811308593752, 152.87405654296876, 76.43702827148438, 38.21851413574219, 19.109257067871095, 9.5546285339355475, 4.7773142669677737, 2.3886571334838869, 1.1943285667419434, 0.59716428337097172, 0.29858214168548586],             //maxResolution: 78271.51695,             projection: new OpenLayers.Projection("EPSG:900913"),             //units: "meters",             numZoomLevels: 20         };              map = new OpenLayers.Map("basicMap", map_options);         map.addControl(new OpenLayers.Control.PanZoomBar());         map.addControl(new OpenLayers.Control.LayerSwitcher());         var mapnik = new OpenLayers.Layer.OSM();         // var tilescache = new OpenLayers.Layer.TileCache();         tiles = new OpenLayers.Layer.WMS("OSM@Mapnik", "tilecache/tilecache.cgi?",{             layers: 'osm', format: 'image/png', srs: 'EPSG:900913'         });              map.addLayer(mapnik);         map.addLayer(tiles);         map.setCenter(new OpenLayers.LonLat(8, 46.5)  // Center of the map           .transform(             new OpenLayers.Projection("EPSG:4326"), // transform from WGS 1984             new OpenLayers.Projection("EPSG:900913") // to Spherical Mercator Projection           ), 8 // Zoom level         );       }


title: How to get the data in dbf file while creating shapefile using arcGis

text: I am using arcObject dll to create shapefile. I am successfully able to create the shape file but I am having some issues with populating the shapefile data into .dbf file. At present the .dbf file just contains the column names but the dataset is empty. ![enter image description here](http://i.stack.imgur.com/DyJZJ.png) I adding the data in db file as follows                public IFields CreateFieldsCollectionForFeatureClass(ISpatialReference spatialReference, string ViewType, DataSet ds)         {             //IFeatureLayer featureLayer = (IFeatureLayer)annotationLayer;             IFeatureClassDescription fcDesc = new FeatureClassDescriptionClass();             IObjectClassDescription ocDesc = (IObjectClassDescription)fcDesc;                  // Create the fields using the required fields method.             IFields fields = ocDesc.RequiredFields;                  // Locate the shape field with the name from the feature class description.             int shapeFieldIndex = fields.FindField(fcDesc.ShapeFieldName);             IField shapeField = fields.get_Field(shapeFieldIndex);                  // Modify the GeometryDef object before using the fields collection to create a new feature class.             IGeometryDef geometryDef = shapeField.GeometryDef;             IGeometryDefEdit geometryDefEdit = (IGeometryDefEdit)geometryDef;                  // Alter the feature class geometry type to lines (default is polygons).             if (ViewType == ViewsType.Powerplant.ToString())             {                 geometryDefEdit.GeometryType_2 = esriGeometryType.esriGeometryMultipoint;             }             else if (ViewType == ViewsType.Pipeline.ToString())             {                 geometryDefEdit.GeometryType_2 = esriGeometryType.esriGeometryPolyline;             }             else if (ViewType == ViewsType.Facility.ToString())             {                 geometryDefEdit.GeometryType_2 = esriGeometryType.esriGeometryMultipoint;             }             geometryDefEdit.HasM_2 = true;             geometryDefEdit.GridCount_2 = 1;                  // Set the first grid size to zero and allow ArcGIS to determine a valid grid size.             geometryDefEdit.set_GridSize(0, 0);             geometryDefEdit.SpatialReference_2 = spatialReference;                  // Because the fields collection already exists, the AddField method on the IFieldsEdit interface             // will be used to add a field that is not required to the fields collection.             IFieldsEdit fieldsEdit = (IFieldsEdit)fields;             foreach (DataColumn col in ds.Tables[0].Columns)             {                 IField filed = new FieldClass();                 IFieldEdit editFiled = (IFieldEdit)filed;                 editFiled.Type_2 = esriFieldType.esriFieldTypeString;                 editFiled.Name_2 = col.ColumnName;                 editFiled.DefaultValue_2 = "test";//ds.Tables[0].Rows;                 //foreach (DataRow d in ds.Tables[0].Rows)                 //{                 //    editFiled.DefaultValue_2 = d["ID"];                 //}                      //editFiled.DefaultValue_2 = ds.Tables[0].Rows[0].ItemArray[0].ToString();                 editFiled.Length_2 = 50;                 fieldsEdit.AddField(editFiled);             }             //fields = (ESRI.ArcGIS.Geodatabase.IFields)fieldsEdit;                  return fields;         }      How can i get the data of that field in .dbf file. Thanks in advance


title: GEOS and Shapely installation on Windows

text: When trying to install Shapely on my Windows 64bit computer, I cannot get the GEOS library to work. So far, I have run the OSGeo4W installer from which I installed GDAL (I believe the geos library is included in that package). After that, I checked and I have geos_c.dll on my `C:\OSGeo4W\bin` directory, but either I have missed some configuration steps or the library does not work. I need Shapely to work, so I also ran `pip install shapely` after installing GDAL, and it apparently worked (although it could not find the C library for GEOS). In my code, I can import Shapely, but when I try to use it, I get an error telling me `"geos.dll"` is not found. Any help with this will be very appreciated. Thanks!


How can I interact with ArcObjects from a windows service? It appears that I can't work with an AxMapControl object from a console or service, because it complains about the activeX state. How does one interact with ArcObjects from a service?


title: TIN interpolation crash in Qgis

text: I have problems using the function TIN(linear) interpolation in Qgis. The program crashes when larger areas is interpolated. Map properties: vector map, Geometry typ = Point, features in this layer: 523851, EPSG:32633 - WGS 84 / UTM zone 33N, files = .shp When i divide the map into more layer(like 80 000 points each) and then interpolate them together(by adding more layers in the interpolation meny), it works, except some areas still makes the program crash. Any idea how to solve this? Or have someone had/have the same problem? Maybe the problem is: can´t handle large files or the algorithm have problems triangulate between some points patterns.


title: Connecting to ArcGIS Server with IIS 7.5 (ArcObjects/COM)

text: This code works fine on a development box, but fails in a server environment:               GISServerConnection gisServerConnection = new GISServerConnection();     gisServerConnection.Connect("localhost");      The exception is this:               Access is denied. (Exception from HRESULT: 0x80070005 (E_ACCESSDENIED))     at ESRI.ArcGIS.Server.IGISServerConnection.Connect(String machineName)      My best guess is that I need to add a user to the `agsusers` group. Our server (Server 2008 R2) uses IIS 7.5, so I tried adding the application pool identity `ASP.NET v4.0 Classic` to the group. This failed, so I tried adding `IUSR` to the group, which also failed. What could I be doing wrong?


title: How to reclassify only a part of a raster but keep the rest unchanged (ArcGIS or R)

text: I'm trying to reclassify only a certain part of a raster...not the entire raster. Here's the problem. I have a LULC raster of the entire US. There are some issues with the raster that I want to fix, namely some cells that should be water are classified as wetlands. These are mainly off the coast, not inland. However, some cells that are classified as wetlands should really be wetlands (mainly inland) and I don't want these to change. I plan on using a ZIP Code boundary file to determine which cells are considered "off the coast" and which ones are not. I have access to ArcGIS 10.1 and R with the appropriate raster packages I need so a solution in either software would be fine. Can anyone help me figure this out?


title: Disable sounds/beeps in QGIS?

text: Is there any way to disable, mute, or otherwise change the incessant message 'beep' in QGIS? I need volume on the PC for other unrelated notiifications yet am tired of apologising for the echoing cacophony to workmates whenever I'm using QGIS. Thanks


title: Insert a loading message into DataGrid

text: I'd like to show a loading message (Loading...) that shows while the datagrid is being populated. It would show when the search button is clicked and disappear when the grid results have populated. I've seen a few possible options including dojox.widget.Standby however, it looks like there is something that is already embedded in the dojox.grid.DataGrid I can't find any documentation on how to access it or show it in my app. Anybody had any luck inserting a loading message in their app?


title: How to upload lines created by joining longitude, latitude pairs?

text: The data in a Mapinfo table contain longitude and latitude pairs. I have joined the point pairs using :   1. Making the table mappable   2. Setting the projection to WGS84 and    3. update TABLE set obj = createline(x1,y1,x2,y2)    4. The lines are created and   5. Map Catalog is created using Easyloader I wanted to upload it using Easyloader but it shows an error message that it could not open the table mytable.tab (mytable is a tab file). What is the process of uploading any table to the database? The database is MSSQL. Edit --------------------------------------------------- After Peter Horsbøll Møller's suggestion I saved a copy and the error was > A .Net Framework error occoured during execution of user defined routine or > aggregate "geometry": System.FormatException:24117:The lineString input is > not valid because it does not have enough distinct po Data Values:18418


title: geoserver installation instructions for osx

text: I recently setup Geoserver on my MacBook (OS X Mavericks) in anticipation of loading it onto an eventual OS X Mavericks Server machine. While installing on my macbook I found a nice blog with step by step instructions that I followed. Now that my server is here, I can't find that nice blog anymore. Simply downloading the geoserver.dmg file to my mavericks server box and dragging the geoserver.app to my Applicaiotns folder is not cutting the mustard. Can anyone point me to some documentation for setting up geoserver on a production mac server?


title: how to identify more than 500 features?

text: I am new to arcgis. and use Arcgis 10 I have done the closest facility analyst. And there are many routes laying on my network. I want to know how many routes use the road section(one feature in network layer), and so on for all of the network. When I use identify, it can just return at most 500 features. Question1 : **How can I get more?**   I've google for this, I know I have to edit the "maxcount". but I can't find that! I have no idea. Please answer specific, thanks a lot! Question2 :Is there any easier way to identify the whole network? like model builder or other tool?


title: Deleting large numbers of fields in ArcGIS 10

text: It there a good strategy for deleting a few thousand fields from the attribute table of a shapefile in ArcGIS using either ArcMap or ArcCatalog? I know how to manually right-click and choose 'delete field;; I'm wondering if there is any other strategy. If I could export the shapefile with columns selected, that would also be fine, but I don't think that can be done.


title: How can I set fill colors using hex value in arcgis.com webmap?

text: Has anybody been able to use hex values (e.g. #FF0000) for setting polygon fill symbols in an arcgis.com 'webmap'? I see predefined colors to select from, but no hex input box. Thanks in advance for any help.


title: Cannot uninstall OpenGeo Suite 3.0.2

text: I cannot uninstall opengeo suite 3.0.2. I tried uninstalling and it would timeout on me. So I deleted the directory and now it won't let me uninstall or reinstall. My goal is to get OpenGeo 3.0.2 working with postgreSQL 9.1 and postGIS. I had it working on Windows 8 (32 bit) and am now trying to move my data to Windows 8.1 (64 bit).


title: Can Spatial Analyst extension be used with ArcGIS Server Standard?

text: I have a model written in ModelBuilder that uses a number of tools in Spatial Analyst extension, e.g. Raster Calculator and Focal Statistics, etc. It runs perfectly with ArcMap. Can I publish it on ArcGIS Server Standard as Geoprocessing service and trigger it from there? I am with 10.1. Having read the matrix, I found it stated ArcGIS Server Standard only supports Schematic and Network Analysts. Do I have to run my Model off ArcGIS Server Advanced?


title: Whats the best way to create bathymetric contour lines?

text: I have a raster dataset with both Topographic and Bathymetric stored in it. What is the best way to create exact contours of the lake-bed? This becomes tricky when I want to make the contour for depth = "zero", when in fact it is not as easy as it seems, since I get contour out of the lake "boundary". Then the question becomes how we really define a lake boundary. Thanks P.S. I'm using ArcGIS-10.1 and model builder.


title: Cannot start OpenGeo Suite 3.1

text: I try to start the suite and it says: "The Suite executable (C:\Program Files\OpenGeo\OpenGeoSuite\opengeo- suite.bat) cannot be found. The 'suite_exe' value must be set in your config.ini to a valid executable." OpenGeo is actually installed in C:\Program Files (x86)\ Does this make a difference? C:\Program Files (x86)\OpenGeo\OpenGeoSuite\opengeo-suite.bat exists and is correctly identified in the first config.ini file. There are two config.ini files. One is in dashboard\Resources and the other is in webapps\dashboard. Both have not be modified. The first config.ini file seems to correctly specify the link to opengeo- suit.bat:               suite_host=localhost     suite_port=8080     suite_stop_port=8079     suite_exe=C:\Program Files (x86)\OpenGeo\OpenGeo Suite\opengeo-suite.bat     suite_dir=C:\Program Files (x86)\OpenGeo\OpenGeo Suite          geoserver_data_dir=C:\Users\Administrator\.opengeo\data_dir     geoserver_username=admin     geoserver_password=geoserver          pgsql_port=54321     pgadmin_path=C:\Program Files (x86)\OpenGeo\OpenGeo Suite\pgsql\9.1\bin\pgadmin3.exe     pgshapeloader_path=C:\Program Files (x86)\OpenGeo\OpenGeo Suite\pgsql\9.1\bin\shp2pgsql-gui.exe      The second config.ini file:               suite_host=localhost     suite_port=8080     suite_stop_port=8079     suite_exe=@SUITE_EXE@     suite_dir=@SUITE_DIR@          geoserver_data_dir=@GEOSERVER_DATA_DIR@     geoserver_username=admin     geoserver_password=geoserver          pgsql_port=@PGSQL_PORT@     pgadmin_path=@PGADMIN_PATH@     pgshapeloader_path=@PGSHAPELOADER_PATH@


title: QGIS Custom CRS has some glitches

text: I am georectifying scanned versions of Survey of India maps, which are in Everest spheroid and india_nepal datum. The exact parameters are: +proj=longlat +a=6377301.243 +b=6356100.2284 +towgs84=295,736,257,0,0,0,0 +no_defs I encountered two issues, one minor and one more significant.   1. I define a custom CRS using these parameters. Give it a name "Everest_India_Nepal_custom". Save the CRS. Use it when georectifying a toposheet. Save the toposheet with this new CRS. (gtif file). But when I reopen that gtif file in a new project, its properties show the CRS parameters, but do not show the custom name I had given for ease of recognition. When saving it said :User 100001: Everest... But when opeining it says "User 100009: *Generated CRS....". So it does not remember the name I gave. Although the name still shows up in CRS options.   2. bigger problem: if I try to validate this custom CRS by using the "calculate" button at the bottom of the custom CRS window, i get rubbish. But the CRS definition is correct, and it works--when my rectified topo is placed on google earth, it sits perfectly (error <10m). Any help greatly appreciated.


title: Polygon Neighbors not found in Proximity Toolset of ArcGIS Desktop 10.0?

text: I cannot find the Polygon Neighbors tool. I have searched for it as I usually do and I have also browsed the proximity toolset but it is not there. The other Tools, e.g. Buffer and Near are there and work fine. Here is a screenshot: ![enter image description here](http://i.stack.imgur.com/kNnSc.png)


title: How to Comment Out Definition Queries

text: I use Definition Queries a lot in my daily GIS life, but one of the many things I've found annoying about ArcMap is the inability to comment out Definition Queries. I'd like to have this option as sometimes I want to turn on/off the Definition Query or I want to have multiple Queries on a single Shapefile, but only one active at a time. So far to get around this, I usually just have a notepad document in the background and copy and paste as necessary. But the ideal solution would be to use the definition query tab in the layer properties to replace this document. I might have multiple different queries on a layer, but all commented out. As I understand, Definition Queries are written in SQL. To comment in SQL I believe you either start a line with: **--** or encompass a statement in brackets: **{ }**. I've tried both of these:               Original: "DATA_AQUISITION" LIKE 'MISSING XSECTION'     Dashes:   --"DATA_AQUISITION" LIKE 'MISSING XSECTION'     Brackets: {"DATA_AQUISITION" LIKE 'MISSING XSECTION'}      The last two statements return a SQL expression error and won't draw in ArcMap. I'm not looking for a programming solution, as that would just be too much work. Any help would be appreciated, thanks!


title: Converting polygon values to a vector grid using PostGIS

text: What is the best way of converting a value from irregularly shaped poylgons to a polygon grid using PostGIS? I have a set of polygons that look like this: ![enter image description here](http://i.stack.imgur.com/6o1eY.png) and I would like to overlay a grid on these polygons, and assign values from the polygon to the grid value, proportional to the area of the polygon that intersects the grid. For example, where a grid intersects more than one polygon, I want to use the fraction of each polygon area that intersects with the gridcell to proportionally calculate the value. ![enter image description here](http://i.stack.imgur.com/J43LV.png) My approach for doing this is as follows. Assuming each polygon has a value X_i and an area A_i, I calculate the density using D_i = X_i / A_i. Then, for each gridcell, I multiply the fraction of the polygon area that intersects that gridcell by the appropriate density value. I think this approach is satisfactory but I have not figured out a way of abstracting it so that it is easier for me to do this regularly. I am curious to learn if there is a better way to do this, if there is a term for this process and if there is a function that exists to perform a calculation of this type.


title: Why do I get 'KeyError: 0' when using Line Intersection tool?

text: I'm using Quantum GIS version `1.8.0-Lisboa` and trying to create intersection point using Intersection tool. (Vector -> Analysis tools -> Line Intersections ) But when I do that for my shapefile, It pops up the following error in the dialog box:               An error has occured while executing Python code:          Traceback (most recent call last):       File "/usr/share/qgis/python/plugins/fTools/tools/doIntersectLines.py", line 96, in accept         self.compute(line1, line2, field1, field2, outPath, self.progressBar)       File "/usr/share/qgis/python/plugins/fTools/tools/doIntersectLines.py", line 174, in compute         outFeat.addAttribute(0, atMap1[index1])     KeyError: 0          Python version:     2.7.3 (default, Aug  1 2012, 05:27:35)      [GCC 4.6.3]               QGIS version:     1.8.0-Lisboa Lisboa, exported          Python path: ['/usr/share/qgis/python', '/home/sibi/.qgis//python', '/home/sibi/.qgis//python/plugins', '/usr/share/qgis/python/plugins', '/usr/local/lib/python2.7/dist-packages/identi-0.1.3-py2.7.egg', '/usr/local/lib/python2.7/dist-packages/pyipmsg-0.0.1.0-py2.7.egg', '/usr/local/lib/python2.7/dist-packages', '/usr/local/lib/python2.7/dist-packages/rope-0.9.4-py2.7.egg', '/usr/local/lib/python2.7/dist-packages/ropemode-0.2-py2.7.egg', '/usr/local/lib/python2.7/dist-packages/ropemacs-0.7-py2.7.egg', '/usr/lib/python2.7', '/usr/lib/python2.7/plat-linux2', '/usr/lib/python2.7/lib-tk', '/usr/lib/python2.7/lib-old', '/usr/lib/python2.7/lib-dynload', '/usr/local/lib/python2.7/dist-packages', '/usr/lib/python2.7/dist-packages', '/usr/lib/python2.7/dist-packages/PIL', '/usr/lib/python2.7/dist-packages/gst-0.10', '/usr/lib/python2.7/dist-packages/gtk-2.0', '/usr/lib/pymodules/python2.7', '/usr/lib/python2.7/dist-packages/ubuntu-sso-client', '/usr/lib/python2.7/dist-packages/ubuntuone-client', '/usr/lib/python2.7/dist-packages/ubuntuone-control-panel', '/usr/lib/python2.7/dist-packages/ubuntuone-couch', '/usr/lib/python2.7/dist-packages/ubuntuone-installer', '/usr/lib/python2.7/dist-packages/ubuntuone-storage-protocol', '/usr/lib/python2.7/dist-packages/wx-2.8-gtk2-unicode', '/usr/share/qgis/python/plugins/fTools/tools']      Is it some bug in the Quantum GIS program or I'm missing something here?


title: Lower the level of a label in QGIS

text: I am finding when I use label vector lines using the newer labeling tool that the labels always appear on top of the lines of that layer. And they appear on top of all other layers in my project. Is there a way to have these labels lower down (i.e. not on top of other layers). Thanks, David


title: How to Force the Showing of Labels of a Vector File in QGIS

text: I imported some OSM data. I then label the lines (roads) with their "name" field using the new label layering settings dialog box. When I zoom to different scale levels I find that while all the lines remain, some of the labels appear / disappear. In any particular zoom level if I adjust the size of the text of the labels, more or less labels show (i.e. when I make the text larger less labels show and when I make the text smaller more labels show). While in general this is good because it creates less clutter and keeps labels from overlapping, etc. there are times when I want to see some of these labels and I can't unless I adjust the text size or zoom in further. I also can't find any control over this either. So, I guess I am wondering if there is a way to force all the labels to show no matter what my text size or scale is. Another option would be a way to turn on a particular label that I wanted. In fact the latter option would be best. Thanks in advance for any help with this. David


title: can't geocode with PostGIS 2.1.1

text: I don't know where's the mistake... I am using the process from the PostGIS documentation (http://postgis.net/docs/manual-2.1/postgis_installation.html#loading_extras_tiger_geocoder) . I created a database, installed the extensions, installed wget and 7zip to match the locations outlined in the doumentation, created profile in loader_platform (based on 'windows' profile, I am on windows 8), however when I try to run the Loader_Generate_Nation_Script generated by PGAdminIII I get lots of errors. The files are downloaded but I don't see any indication that any unzipping is done. If I ignore the errors and continue the process as explained in the documnetation, of course goeocoding doesn't work. I have the Loader_Generate_Nation_Script pasted here, and also the output from the Command Prompt Window below. Can anybody help ? LoaderGenerateNationScript:               set TMPDIR=\gisdata\temp\     set UNZIPTOOL="C:\Program Files\7-Zip\7z.exe"     set WGETTOOL="C:\wget\wget.exe"     set PGBIN=C:\Program Files\PostgreSQL\9.3\bin\     set PGPORT=5432     set PGHOST=localhost     set PGUSER=postgres     set PGPASSWORD=password1     set PGDATABASE=my6thGIS     set PSQL="%PGBIN%psql"     set SHP2PGSQL="%PGBIN%shp2pgsql"     cd \gisdata          cd \gisdata     %WGETTOOL% ftp://ftp2.census.gov/geo/tiger/TIGER2013/STATE/ --no-parent --relative --recursive --level=1 --accept=zip --mirror --reject=html      cd \gisdata/ftp2.census.gov/geo/tiger/TIGER2013/STATE     del %TMPDIR%\*.* /Q     %PSQL% -c "DROP SCHEMA IF EXISTS tiger_staging CASCADE;"     %PSQL% -c "CREATE SCHEMA tiger_staging;"     %PSQL% -c "DO language 'plpgsql' $$ BEGIN IF NOT EXISTS (SELECT * FROM information_schema.schemata WHERE schema_name = 'tiger_data' ) THEN CREATE SCHEMA tiger_data; END IF;  END $$"     for /r %%z in (tl_*state.zip ) do %UNZIPTOOL% e %%z  -o%TMPDIR%      cd %TMPDIR%     %PSQL% -c "CREATE TABLE tiger_data.state_all(CONSTRAINT pk_state_all PRIMARY KEY (statefp),CONSTRAINT uidx_state_all_stusps  UNIQUE (stusps), CONSTRAINT uidx_state_all_gid UNIQUE (gid) ) INHERITS(state); "     %SHP2PGSQL% -c -s 4269 -g the_geom   -W "latin1" tl_2013_us_state.dbf tiger_staging.state | %PSQL%     %PSQL% -c "SELECT loader_load_staged_data(lower('state'), lower('state_all')); "         %PSQL% -c "CREATE INDEX tiger_data_state_all_the_geom_gist ON tiger_data.state_all USING gist(the_geom);"         %PSQL% -c "VACUUM ANALYZE tiger_data.state_all"     cd \gisdata     %WGETTOOL% ftp://ftp2.census.gov/geo/tiger/TIGER2013/COUNTY/ --no-parent --relative --recursive --level=1 --accept=zip --mirror --reject=html      cd \gisdata/ftp2.census.gov/geo/tiger/TIGER2013/COUNTY     del %TMPDIR%\*.* /Q     %PSQL% -c "DROP SCHEMA IF EXISTS tiger_staging CASCADE;"     %PSQL% -c "CREATE SCHEMA tiger_staging;"     %PSQL% -c "DO language 'plpgsql' $$ BEGIN IF NOT EXISTS (SELECT * FROM information_schema.schemata WHERE schema_name = 'tiger_data' ) THEN CREATE SCHEMA tiger_data; END IF;  END $$"     for /r %%z in (tl_*county.zip ) do %UNZIPTOOL% e %%z  -o%TMPDIR%      cd %TMPDIR%     %PSQL% -c "CREATE TABLE tiger_data.county_all(CONSTRAINT pk_tiger_data_county_all PRIMARY KEY (cntyidfp),CONSTRAINT uidx_tiger_data_county_all_gid UNIQUE (gid)  ) INHERITS(county); "      %SHP2PGSQL% -c -s 4269 -g the_geom   -W "latin1" tl_2013_us_county.dbf tiger_staging.county | %PSQL%     %PSQL% -c "ALTER TABLE tiger_staging.county RENAME geoid TO cntyidfp;  SELECT loader_load_staged_data(lower('county'), lower('county_all'));"         %PSQL% -c "CREATE INDEX tiger_data_county_the_geom_gist ON tiger_data.county_all USING gist(the_geom);"         %PSQL% -c "CREATE UNIQUE INDEX uidx_tiger_data_county_all_statefp_countyfp ON tiger_data.county_all USING btree(statefp,countyfp);"         %PSQL% -c "CREATE TABLE tiger_data.county_all_lookup ( CONSTRAINT pk_county_all_lookup PRIMARY KEY (st_code, co_code)) INHERITS (county_lookup);"         %PSQL% -c "VACUUM ANALYZE tiger_data.county_all;"         %PSQL% -c "INSERT INTO tiger_data.county_all_lookup(st_code, state, co_code, name) SELECT CAST(s.statefp as integer), s.abbrev, CAST(c.countyfp as integer), c.name FROM tiger_data.county_all As c INNER JOIN state_lookup As s ON s.statefp = c.statefp;"         %PSQL% -c "VACUUM ANALYZE tiger_data.county_all_lookup;"       Windows8 Command Prompt output:               Microsoft Windows [Version 6.3.9600]     (c) 2013 Microsoft Corporation. All rights reserved.          C:\Users\user1>set TMPDIR=\gisdata\temp\          C:\Users\user1>set UNZIPTOOL="C:\Program Files\7-Zip\7z.exe"          C:\Users\user1>set WGETTOOL="C:\wget\wget.exe"          C:\Users\user1>set PGBIN=C:\Program Files\PostgreSQL\9.3\bin\          C:\Users\user1>set PGPORT=5432          C:\Users\user1>set PGHOST=localhost          C:\Users\user1>set PGUSER=postgres          C:\Users\user1>set PGPASSWORD=password1          C:\Users\user1>set PGDATABASE=my6thGIS          C:\Users\user1>set PSQL="%PGBIN%psql"          C:\Users\user1>set SHP2PGSQL="%PGBIN%shp2pgsql"          C:\Users\user1>cd \gisdata          C:\gisdata>     C:\gisdata>cd \gisdata          C:\gisdata>%WGETTOOL% ftp://ftp2.census.gov/geo/tiger/TIGER2013/STATE/ --no-pare     nt --relative --recursive --level=1 --accept=zip --mirror --reject=html     --2014-03-30 11:08:52--  ftp://ftp2.census.gov/geo/tiger/TIGER2013/STATE/                => `ftp2.census.gov/geo/tiger/TIGER2013/STATE/.listing'     Resolving ftp2.census.gov... 148.129.75.35     Connecting to ftp2.census.gov|148.129.75.35|:21... connected.     Logging in as anonymous ... Logged in!     ==> SYST ... done.    ==> PWD ... done.     ==> TYPE I ... done.  ==> CWD /geo/tiger/TIGER2013/STATE ... done.     ==> PASV ... done.    ==> LIST ... done.              [                                    ] 197         --.-K/s   in 0s          2014-03-30 11:08:53 (4.83 MB/s) - `ftp2.census.gov/geo/tiger/TIGER2013/STATE/.li     sting' saved [197]          --2014-03-30 11:08:53--  ftp://ftp2.census.gov/geo/tiger/TIGER2013/STATE/tl_2013     _us_state.zip                => `ftp2.census.gov/geo/tiger/TIGER2013/STATE/tl_2013_us_state.zip'     ==> CWD not required.     ==> PASV ... done.    ==> RETR tl_2013_us_state.zip ... done.     Length: 8599274 (8.2M)          100%[======================================>] 8,599,274    251K/s   in 24s          2014-03-30 11:09:17 (344 KB/s) - `ftp2.census.gov/geo/tiger/TIGER2013/STATE/tl_2     013_us_state.zip' saved [8599274]          FINISHED --2014-03-30 11:09:17--     Downloaded: 2 files, 8.2M in 24s (344 KB/s)          C:\gisdata>cd \gisdata/ftp2.census.gov/geo/tiger/TIGER2013/STATE          C:\gisdata\ftp2.census.gov\geo\tiger\TIGER2013\STATE>del %TMPDIR%\*.* /Q          C:\gisdata\ftp2.census.gov\geo\tiger\TIGER2013\STATE>%PSQL% -c "DROP SCHEMA IF E     XISTS tiger_staging CASCADE;"     NOTICE:  schema "tiger_staging" does not exist, skipping     DROP SCHEMA          C:\gisdata\ftp2.census.gov\geo\tiger\TIGER2013\STATE>%PSQL% -c "CREATE SCHEMA ti     ger_staging;"     CREATE SCHEMA          C:\gisdata\ftp2.census.gov\geo\tiger\TIGER2013\STATE>%PSQL% -c "DO language 'plp     gsql' $$ BEGIN IF NOT EXISTS (SELECT * FROM information_schema.schemata WHERE sc     hema_name = 'tiger_data' ) THEN CREATE SCHEMA tiger_data; END IF;  END $$"     DO          C:\gisdata\ftp2.census.gov\geo\tiger\TIGER2013\STATE>for /r %%z in (tl_*state.zi     p ) do %UNZIPTOOL% e %%z  -o%TMPDIR%     %%z was unexpected at this time.          C:\gisdata\ftp2.census.gov\geo\tiger\TIGER2013\STATE>cd %TMPDIR%          C:\gisdata\temp>%PSQL% -c "CREATE TABLE tiger_data.state_all(CONSTRAINT pk_state     _all PRIMARY KEY (statefp),CONSTRAINT uidx_state_all_stusps  UNIQUE (stusps), CO     NSTRAINT uidx_state_all_gid UNIQUE (gid) ) INHERITS(state); "     CREATE TABLE          C:\gisdata\temp>%SHP2PGSQL% -c -s 4269 -g the_geom   -W "latin1" tl_2013_us_stat     e.dbf tiger_staging.state | %PSQL%     Unable to open tl_2013_us_state.shp or tl_2013_us_state.SHP.     tl_2013_us_state.dbf: dbf file (.dbf) can not be opened.          C:\gisdata\temp>%PSQL% -c "SELECT loader_load_staged_data(lower('state'), lower(     'state_all')); "     NOTICE:  INSERT INTO tiger_data.state_all(region,division,statefp,statens,stusps     ,name,lsad,mtfcc,funcstat,aland,awater,intptlat,intptlon,the_geom) SELECT  FROM     tiger_staging.state;     CONTEXT:  SQL function "loader_load_staged_data" statement 1     ERROR:  syntax error at or near "FROM"     LINE 1: ...,aland,awater,intptlat,intptlon,the_geom) SELECT  FROM tiger...                                                                  ^     QUERY:  INSERT INTO tiger_data.state_all(region,division,statefp,statens,stusps,     name,lsad,mtfcc,funcstat,aland,awater,intptlat,intptlon,the_geom) SELECT  FROM t     iger_staging.state;     CONTEXT:  PL/pgSQL function loader_load_staged_data(text,text,text[]) line 24 at      EXECUTE statement     SQL function "loader_load_staged_data" statement 1          C:\gisdata\temp>%PSQL% -c "CREATE INDEX tiger_data_state_all_the_geom_gist ON ti     ger_data.state_all USING gist(the_geom);"     CREATE INDEX          C:\gisdata\temp>%PSQL% -c "VACUUM ANALYZE tiger_data.state_all"     VACUUM          C:\gisdata\temp>cd \gisdata          C:\gisdata>%WGETTOOL% ftp://ftp2.census.gov/geo/tiger/TIGER2013/COUNTY/ --no-par     ent --relative --recursive --level=1 --accept=zip --mirror --reject=html     --2014-03-30 11:09:18--  ftp://ftp2.census.gov/geo/tiger/TIGER2013/COUNTY/                => `ftp2.census.gov/geo/tiger/TIGER2013/COUNTY/.listing'     Resolving ftp2.census.gov... 148.129.75.35     Connecting to ftp2.census.gov|148.129.75.35|:21... connected.     Logging in as anonymous ... Logged in!     ==> SYST ... done.    ==> PWD ... done.     ==> TYPE I ... done.  ==> CWD /geo/tiger/TIGER2013/COUNTY ... done.     ==> PASV ... done.    ==> LIST ... done.              [                                    ] 198         --.-K/s   in 0s          2014-03-30 11:09:19 (12.6 MB/s) - `ftp2.census.gov/geo/tiger/TIGER2013/COUNTY/.l     isting' saved [198]          --2014-03-30 11:09:19--  ftp://ftp2.census.gov/geo/tiger/TIGER2013/COUNTY/tl_201     3_us_county.zip                => `ftp2.census.gov/geo/tiger/TIGER2013/COUNTY/tl_2013_us_county.zip'          ==> CWD not required.     ==> PASV ... done.    ==> RETR tl_2013_us_county.zip ... done.     Length: 73500446 (70M)          100%[======================================>] 73,500,446   186K/s   in 3m 59s          2014-03-30 11:13:18 (301 KB/s) - `ftp2.census.gov/geo/tiger/TIGER2013/COUNTY/tl_     2013_us_county.zip' saved [73500446]          FINISHED --2014-03-30 11:13:18--     Downloaded: 2 files, 70M in 3m 59s (301 KB/s)          C:\gisdata>cd \gisdata/ftp2.census.gov/geo/tiger/TIGER2013/COUNTY          C:\gisdata\ftp2.census.gov\geo\tiger\TIGER2013\COUNTY>del %TMPDIR%\*.* /Q          C:\gisdata\ftp2.census.gov\geo\tiger\TIGER2013\COUNTY>%PSQL% -c "DROP SCHEMA IF     EXISTS tiger_staging CASCADE;"     DROP SCHEMA          C:\gisdata\ftp2.census.gov\geo\tiger\TIGER2013\COUNTY>%PSQL% -c "CREATE SCHEMA t     iger_staging;"     CREATE SCHEMA          C:\gisdata\ftp2.census.gov\geo\tiger\TIGER2013\COUNTY>%PSQL% -c "DO language 'pl     pgsql' $$ BEGIN IF NOT EXISTS (SELECT * FROM information_schema.schemata WHERE s     chema_name = 'tiger_data' ) THEN CREATE SCHEMA tiger_data; END IF;  END $$"     DO          C:\gisdata\ftp2.census.gov\geo\tiger\TIGER2013\COUNTY>for /r %%z in (tl_*county.     zip ) do %UNZIPTOOL% e %%z  -o%TMPDIR%     %%z was unexpected at this time.          C:\gisdata\ftp2.census.gov\geo\tiger\TIGER2013\COUNTY>cd %TMPDIR%          C:\gisdata\temp>%PSQL% -c "CREATE TABLE tiger_data.county_all(CONSTRAINT pk_tige     r_data_county_all PRIMARY KEY (cntyidfp),CONSTRAINT uidx_tiger_data_county_all_g     id UNIQUE (gid)  ) INHERITS(county); "     CREATE TABLE          C:\gisdata\temp>%SHP2PGSQL% -c -s 4269 -g the_geom   -W "latin1" tl_2013_us_coun     ty.dbf tiger_staging.county | %PSQL%     Unable to open tl_2013_us_county.shp or tl_2013_us_county.SHP.     tl_2013_us_county.dbf: dbf file (.dbf) can not be opened.          C:\gisdata\temp>%PSQL% -c "ALTER TABLE tiger_staging.county RENAME geoid TO cnty     idfp;  SELECT loader_load_staged_data(lower('county'), lower('county_all'));"     ERROR:  relation "tiger_staging.county" does not exist          C:\gisdata\temp>%PSQL% -c "CREATE INDEX tiger_data_county_the_geom_gist ON tiger     _data.county_all USING gist(the_geom);"     CREATE INDEX          C:\gisdata\temp>%PSQL% -c "CREATE UNIQUE INDEX uidx_tiger_data_county_all_statef     p_countyfp ON tiger_data.county_all USING btree(statefp,countyfp);"     CREATE INDEX          C:\gisdata\temp>%PSQL% -c "CREATE TABLE tiger_data.county_all_lookup ( CONSTRAIN     T pk_county_all_lookup PRIMARY KEY (st_code, co_code)) INHERITS (county_lookup);     "     CREATE TABLE          C:\gisdata\temp>%PSQL% -c "VACUUM ANALYZE tiger_data.county_all;"     VACUUM          C:\gisdata\temp>%PSQL% -c "INSERT INTO tiger_data.county_all_lookup(st_code, sta     te, co_code, name) SELECT CAST(s.statefp as integer), s.abbrev, CAST(c.countyfp     as integer), c.name FROM tiger_data.county_all As c INNER JOIN state_lookup As s      ON s.statefp = c.statefp;"     INSERT 0 0          C:\gisdata\temp>%PSQL% -c "VACUUM ANALYZE tiger_data.county_all_lookup;"     VACUUM          C:\gisdata\temp>


title: python.exe has stopped working

text: A python script was written about 18 months ago by a person who has now left. It produced the required outputs then. I've been asked to run it again but with different (finer resolution) data inputs. The input dataset has been split into 20 sub-sets of approx 2,700 data points each. However, the script crashes ("python.exe has stopped working") after approx 300 data points have been processed (range 295 to 306 and does NOT always fail on the same record). As its old(ish), the script was written using arcgisscripting and not arcpy. Broadly it does the following using cursors: a) For a given point, calculate the cost distance (using gp.CostDistance_sa) with a cutoff of 60 minutes travel time. b)Calls gp.ExtractValuesToPoints_sa to extract all the individual values at each data point and outputs a feature class to a file geodatabase. c) Reads the feature class created in b) above and writes the values to a CSV file (omitting any points with "No Data" (value -9999)). Repeats a), b) and c) for all remaiming data points in the input file. Processing time is approx. 1 minute per data point on average. The PC has a quad core Intel i7-2720QM CPU running at 2.20GHz with 8GB RAM running Windows 7 (64 bit). Python version is 2.6.6 (shell also states "[MSC v,1500 32 bit (Intel)] on win32). Not sure if relevant but ARCMAP 10.0 (SP4) is also installed. I've tried running it on a different PC (so far without crashing). Currently the job is running successfully (but more slowly) on an older PC and has reached 419 records without crashing. This has an Intel Core 2 DUO E7500 processor running at 2.93GHz with 4 GB RAM and 64bit Windows 7. Python version 2.5.1 (shell also states "[MSC v,1310 32 bit (Intel)] on win32). ARCMAP 9.3 is installed (no mention of any Service Packs). Please can someone offer some advice about why the script seems to work for a while then crash and how to resolve it? The fact that a different PC appears (so far) to handle the script suggests something "environmental". Thanks in anticipation. * * * Thanks to everyone for your suggestions. I'll try these. As an update, the PC running ARCGIS 9.3 is still successfully processing the data and has reached 1,300 data points processed (and still counting). A colleague also ran the data on their PC running ARCGIS 10.1 - it crashed after 267 records on two seperate occasions. Although not conclusive, the common thread seems to be that ARC 9.3 will process the data but ARC 10.x will not. Any further contributions gratefully received.


title: Licensing of QGIS libraries

text: what type of license are subject to the following QGIS libraries? libqgis_analysis.so libqgis_core.so libqgis_gui.so I would like to use the following libraries released in QGIS package in a "close" (not open source) application. Are there any derived-license isusses ??


title: How to read/write ENVI metadata using GDAL?

text: I am trying to use GDAL to read and write images in ENVI format (for example a .bsq file with an associated .hdr file). The GDAL Formats Info page states "Starting with GDAL 1.10, all ENVI header fields will be stored in the ENVI metadata domain.". If that is the case, how should I access them? I am using the Python bindings with GDAL 1.9.2 (the latest at this time of writing), and have tried getting all of the metadata:               image.GetMetadata()      which gives me a list of the band 'names' as displayed in ENVI, but that's about it. I've also tried specifically asking for the ENVI metadata domain, as instructed on the GDAL information page linked above:               image.GetMetadata('ENVI')      but that doesn't return anything. Is there a way to get other metadata from the ENVI file, or will I have to parse the header file myself? Furthermore, even writing the band names (the ones that seem to be returned when I run `GetMetadata()`) using `SetMetadata()` doesn't seem to put them in the .hdr file. Can anyone shed any light on this?


title: How to handle/ avoid intersect errors when clipping the difference between two layers?

text: Using QGIS + Grass, I've added a shapefile of coastline and islands to my map. I then buffer this shapefile to create a coastal zone of a known width. At this stage the buffer consists of the offshore (buffered) area + the area occupied by the land/ islands. As I am only interested in the offshore part, I then "clip" the land shapefile from the buffer. I have tried this using both Vector>Geoprocessing Tools>Difference and GRASS v.overlay.xor. Either approach results in geometry errors (lots of them). e.g.,               segments 0 and 24 of line 43 intersect at xx.xxx, yy.yyyy      This appears to be due to narrow connections between land areas (yellow in the picture below) resulting in intersecting nodes after the ocean (blue) is clipped. ![enter image description here](http://i.stack.imgur.com/7LgeN.png) My land data comes from an external source, and I don't have a lot of scope to simplify or adjust this. What is the best way to handle these errors, or to avoid them in the first place? I'm not a GIS expert, so there's a good chance that I'm approaching this problem from the wrong end. Would appreciate any tips or advice to avoid these errors.


title: How to use Openweathermap's wms in ArcGIS Online?

text: I would like to add this wms service to ArcGIS online. http://wms.openweathermap.org/service Unlike in QGIS, when I added it to the web services, there is no way for me to choose which layer to add in. Is there a way to specify it in ArcIS online?


title: Route analysis in R based on the OSM transport network: possible?

text: I'm analysing the likely paths taken by commuters on their way to work, based on Census flow data, job density distributions and the transport network. At the moment I'm using pgRouting and the pgRouting add-on for QGIS (described here) for my analysis: ![routing in qgis](http://i.stack.imgur.com/p7Ojc.jpg) However, I'd like all the steps I take to be full replicable (see excellent article on this here), then I'd prefer to do my analysis from the command-line and save the script files for future reference. R-geo packages fulfil this requirement well, as they are free, lightweight, and stable: download a tiny script file and you should be able to re-do all the things I have done (see here for an unrelated example of this on Rpubs). The R packages I've identified for doing this are:   * gdistance published December 2012, looks promising, but no sign of compatibility with osm transport networks   * e1071 seems rather abstract and possibly overkill for what I want to do   * pairdist.lpp function in spatstat package looks promising, but has little documentation Before diving in to these, thought it worth asking people who may have faced similar dilemmas. So: is it possible to do my routing analysis in R, or would it be best to learn PostGIS properly and do it all in pgRouting?


title: Is it possible to do scale dependent rendering in Google maps with KML?

text: I have some data that I would like to display in greater detail as the user zooms into the interactive Gooogle map. Does google support dynamically turning layers on and off depending on the scale?


title: WMS CQL Filter in Geoserver not working on adding to the layer

text: I draw in the init() function a WMS layer an add it to the map:               var lyr_osm = {};     var map = {};          $(document).ready(function () {         map = new OpenLayers.Map('map',                         {                         numZoomLevels: 18,                         projection: googleMercator                         }            );         });              lyr_osm =new OpenLayers.Layer.OSM("OSM");                 lyr_state = new OpenLayers.Layer.WMS("State",                "http://54.228.241.187:80/geoserver/us_admin/wms",                {                'VERSION' : '1.1.0',                layers: "us_admin:states",     //           style: style_state_polyg,                format: "image/png",                transparent: true                },                    {isBaseLayer: false,                     visibility: false,                     sphericalMercator: true,                                             opacity: 0.5}                );              map.addControl(new OpenLayers.Control.LayerSwitcher());         map.addControl(new OpenLayers.Control.MousePosition());         map.addLayers([lyr_osm, lyr_state]);         var mapCenter = new OpenLayers.LonLat(-97, 40);         mapCenter.transform(wgs84, googleMercator);         map.setCenter(mapCenter,4);     }      Then onclick of an input button I pass a location from a google Autocompleter, transform it to EPSG900913, and use it to filter the WMS layer:                $(document).ready(function () {     $('#search1').click(function()                                                {            sdata = $('#searchTextField').val();         if( sdata == undefined || sdata == 'Enter an address...' || sdata < 4 ){         alert('Please insert or choose an address to proceed.');         } else {                 latlong = String(ads.d);            longlat = latlong.split(",");            long = String(longlat[1]);            slong = long.replace(")", "");            lat = String(longlat[0]);            slat = lat.replace("(", "");            position = new OpenLayers.LonLat(slong, slat);            position.transform(wgs84, googleMercator);            var mlat = String(position);            mlatt = mlat.replace("lon=","")            mlattt = mlatt.replace("lat=","");            mlatttt = mlattt.replace(","," ");            lyr_state.mergeNewParams({cql_filter: "contains(geom,POINT(-10625033.067 3488845.830))"});       lyr_state.mergeNewParams({visibility: true});            stbound=lyr_state.getDataExtent();       map.zoomToExtent(stbound);     });     });      The problem is that onclick returns me an aocean area, and not what I want. The new parameters are updated in the DOM but I do not see the result. Any hint? Thanks! Geoserver is 2.2 with no error showing. Java = Oracle Corporation: 1.7.0_11 (Java HotSpot(TM) 64-Bit Server VM) Any hint? Thanks!


title: Show representation properties in object inspector in ArcMap 10

text: When you inspected a feature in the object inspector in ArcMap 9.x (opened by clicking the "Attributes" button in the editor toolbar), you were able to see the properties of ALL the features representations, as individual tabs in the dialog (behind the attribute tab). In ArcGIS 10.0, I have not been able to do this. Only, if the representation is active, can I see its properties in the inspector. This has ENORMOUS (!!!) costs to my work. Until now, my workaround has been to add the layer several times with the different representations, but I am kind of hoping, that there is some way to make. However doing this on a country-wide road network, makes it painstakingly slow to work in ArcMap, so i really hope that one of you has some way of making my problem go away.


title: Generic comparison for kriging in geoR and ArcGIS

text: I am going to do a project based on kriging. The targeted area would be as big as United States and my grid would be 10km*10km. I am wondering the differences of efficiencies for the methods implemented in geoR package and ArcGIS. Have anyone done any research on the performance yet. I would prefer geoR if it works well.


title: How to add a float or integer number (i.e 255) to a raster layer?

text: Is there any code to add a float or integer value with a raster layer? I can easily do mathematical operations between two raster data. But I am unable to do the same between float or integer number and raster layer. For example if want to execute (256-band3), how can I execute?


title: How does versioned views work with multiple editors?

text: I understand that in the database environment, by executing `version_util.set_current_version('version1')` Editor1 is able to write SQL statements against [table_name]_vw views to edit the data at version1. But if Editor2 (from another workstation) executes `version_util.set_current_version('version2')` the [table_name]_vw views should be populated with the data at version2, right? In this case, if Editor1 run a query against the [table_name]_vw, will he still get the data at version1? Or is he still editing the right version?


title: stratigraphic ranking with field calculator?

text: I would like to populate a feature class float attribute with a numerical ranking of a text attribute of a geologic map. I have a field named [sheet_unit] with strings such as "H10, Q" and "J10, Pl" with I would like to relabel in another attribute [Stratigraphic_Rank] with a numerical value such as 0 and 10, etc. I was hoping to use an IF THEN statement but I cannot for the life of figure out how to use the field calculator to do this. Anyone willing to help me start? This and my several variations of it don't work: If [sheet_unit] = "H10, Q" THEN [Stratigraphic_Rank] = 0 ElseIf [sheet_unit] = "J10, Pl" [Stratigraphic_Rank] = 10 EndIf


title: Slippy map tile question and good practice

text: I have created a small slippy map (Leaflet) and have given it a bounding box. If the user tries to pan/slide outside the given bounds then the map bounces back. Also, the slippy map only has tiles within the bounds so if the user pans out it will all be grey area only.All this is what I want but my question is, in the browser console I get " GET: "file location" net::ERR_FILE_NOT_FOUND" which is expected as I don't have those tiles outside of the map bounds. So my questions are: What are some slippy map good practices? Should I include tiles outside of the bounding box even though they will not be used? Are these GET errors bad for performance? Is there a way to not get these GET errors and if so how? Is this ok? Using Leaflet API Any advice/help is greatly appreciated. Thank You


title: Repetitive automated mapping: which platform?

text: I could give a lot of examples, but imagine you have do eg a study in 10 different locations and have to generate a lot of maps with the same layout, legend. Maps can contain vector and raster data. It should be scriptable, because a large number of maps have to be made/remade in the future. Which platform would you use?


title: Why do area calculations return negative values?

text: I have created a shapefile with some polygons and am now trying to calculate the area of these polygons. The original CRS was WGS84, I have reprojected it to WGS 84 / World Mercator so unit is metres. On-the-fly projection is on. If I'm now using the field calculator to calculate the area of the polygons using $area, I first get some really weird values in my output column (far too high) and if I click 'save' or end the toggle mode, the values change to negative values which are the same for all polygons. Switching to another CRS and reprojecting the file again e.g. to GDA94/MGA55 EPSG:28355 and calculating the area creates again funny values, some of them negative. I can't figure out where the problem is and am grateful for any help! Thanks a lot!


Which one is best open source GIS software where we can use Google/bing/yahoo maps as a basemap? i want to update my shp files using any of these basemap.


title: How to Publish a Time Aware Mosaic dataset as an ImageService with Temporal component?

text: I have a Time Aware Mosaic dataset , which contains 5 rasters. I have entered the dates for the given rasters, by adding a Text Field with the dates. (The dates of the images are from 3/1/2013 to 3/5/2013) The Mosaic dataset works great with the Time Slider window in ArcMap. But when I publish it to ArcGIS Server as an Image Service (through ArcCatlog), I do not see the correct Time Extent on the REST end point. The REST end point shows the Time Info as: > Time Info: > > Time Extent: [2013/01/01 00:00:00 UTC, 2013/01/01 00:00:00 UTC]   >  Time Reference: N/A What steps have I missed? Where do you set the Time Extent/Time Reference while publishing the image service? My Configuration: ArcGIS Dekstop 10.0 (ArcInfo)+ ArcGIS Server Advanced+ Image Extension


title: GDAL SQL syntax to add field an put values

text: I try to make a new field in my attribute table and put in values from another field after converting from character to integer but can't get the right syntax. I run this as shell-script in MSYS.               cd D:/GIS_DataBase/CorineLC/shps_app_and_extr/     myfile=extr_and_app.shp          name=${myfile%.shp}          ogrinfo $myfile -sql "ALTER TABLE $name ADD COLUMN code_num int(3)"     ogrinfo $myfile -sql "UPDATE TABLE $name SET code_num = CONVERT(code_06 As int(3))"      Error Message (MSYS):               Kay@KAY-PC /c/users/kay/desktop/bash     $ sh calc_field_shp.sh     Warning 6: Unsupported column type 'int'. Defaulting to VARCHAR     INFO: Open of `extr_and_app.shp'           using driver `ESRI Shapefile' successful.     ERROR 1: SQL Expression Parsing Error: syntax error     INFO: Open of `extr_and_app.shp'           using driver `ESRI Shapefile' successful.      # Edit - Another go with SQLite:               cd D:/GIS_DataBase/CorineLC/shps_app_and_extr/     myfile=extr_and_app.dbf          name=${myfile%.dbf}          ogrinfo $myfile -sql "ALTER TABLE $name DROP COLUMN code_num"          ogrinfo $myfile -sql "ALTER TABLE $name ADD COLUMN code_num integer(3)"     ogrinfo $myfile -dialect SQLite -sql "UPDATE $name SET code_num = CAST(code_06 As integer(3))"      Error message:               Kay@KAY-PC /c/users/kay/desktop/bash     $ sh calc_field_shp.sh     INFO: Open of `extr_and_app.dbf'           using driver `ESRI Shapefile' successful.     INFO: Open of `extr_and_app.dbf'           using driver `ESRI Shapefile' successful.     ERROR 1: SQL Expression Parsing Error: syntax error     INFO: Open of `extr_and_app.dbf'           using driver `ESRI Shapefile' successful.


title: Aligning text/features on a map

text: Sorry for the simple question .... but something odd is happening in ArcMap & I don't know why. When I want to align text or features relative to each other, they are aligning relative to the layout screen instead. For example, if I wanted my north arrow to be perfectly centered in a text box, they are perfectly centered, however, ArcMap is moving them to the center of paper space instead of moving them relatively positioned in relation to the "last" element I selected (the way the command is supposed to work). Am I overlooking something simple here? It's very annoying, as I have to manually move everything back where it belongs every time.


title: Ways of expressing geographic questions in machine readable form

text: A core concept of GIS is answering questions about datasets. From the point of view of a database; SQL with spatial extensions is a way of asking such questions. What other ways can questions be expressed in a machine readable text based form? What are the benefits of different approaches?


title: Does the 'FLEXLM_TIMEOUT' environment variable work with ArcGIS 10 Desktop?

text: In earlier versions of ArcGIS, you could create/set the FLEXLM_TIMEOUT to a higher value for slower connections. I'm currently trying to troubleshoot a system (one where I have minimal control) with slow connections to the license manager. I was able to create a FLEXLM_TIMEOUT user variable and fired up ArcMap several times with no timeout error. But maybe I'm lucky and the network I'm on is having a good afternoon. Does anyone know if this environment variable is actually used by ArcGIS(FLEXnet Publisher)?


title: How to do K means clustering of anistropic diffused image?

text: I am trying to detect edges and i have done a anistropic diffusion of the image to get the diffused image. How do i do the K means clustering now of this diffused image in CIE labspace ? (Not orignal RGB image)


title: How can I add/define a legend for RGB images in QGIS?

text: I have a RGB image in QGIS that is the product of some raster analysis (e.g. QGIS's relief operation, results from sound/traffic models) to which I want to add a legend. I'd like to denote several RGB combinations with a individual labels. How should I go about it? I can't find a way to add legends to multiband color images in QGIS. Ultimately I want to publish the image as a WMS service through GeoServer. Defining a legend on the server (through SLD?) is therefore also an option.


title: Conversion of raster to polygon

text: Using C# and ArcObjects, what is the recommended way to convert a raster file into a polygon. I need this for calling the Near tool so that I can move a station which may be inland to the nearest water feature. I have googled and really not found anything useful yet. There is a lot for taking a polygon to a Raster but not the other way around.


title: QgsMapCanvas displays layer only after i resize the window

text: Im medlling with a standalone PyQGis application (based on 2.0.1). I have no problems with the setup, however when i set up my layers as stated in the docs:               layer = QgsVectorLayer(path, name, provider)     if not layer.isValid():       raise IOError, "Failed to open the layer"          # add layer to the registry     QgsMapLayerRegistry.instance().addMapLayer(layer)          # set extent to the extent of our layer     canvas.setExtent(layer.extent())          # set the map canvas layer set     canvas.setLayerSet( [ QgsMapCanvasLayer(layer) ] )      Everything works as it should, with this one annoying detail. The map canvas widget is not actually rendered until i resize the window containing the widget. When setting up, calls to refresh() method doesn't make a difference. Is this some kind of a bug on 2.0.1?


title: How to examine length of vector lines in qgis/grass?

text: I have vector layer of lines. How can I examine length of lines? ![enter image description here](http://i.stack.imgur.com/2xQeI.png)


title: How to disable and enable drawing in ArcMap using Python?

text: The script that I'm writing repositions two data frames and sets their extents. When doing this, it redraws the entire Active View 4 times, slowing the script substantially. Pressing F9 or clicking the 'Pause Drawing' button before running the script disables the drawing and makes the script run much faster, but I want the script to do this automatically. I want to disable ArcMap 10 drawing at the beginning of the script, and enable it at the end. How do I do this?


title: How to develop WMS services in java?

text: Can anyone point me to any resources regarding developing a custom WMS service in java? I am seeking for a tutorial for an experienced programmer, but including the details of setting up the server and the like as well as some code samples.


title: "Item not found in this collection" when trying to set the Shape on IFeatureBuffer

text: I'm working on code to open a 3rd party data file and convert it into a collection of features in a custom layer based on FeatureLayer. I have had no issue importing the data itself (nodal data for each 'cell' as a feature, such as elevation). Because the data files are on average rather large (10k+ cells), I'm using a cursor implementation with an IFeatureBuffer that is filled with each cell one at a time. However, when I try to set the geometry of the cell via IFeatureBuffer like so, IGeometry cellGeometry = getCellGeometry(cell); featureBuffer.Shape = cellGeometry; I subsequently get the error "Item not found in this collection" called on the second line. My first thought was that I was creating the cell geometry incorrectly in some way, so I tested it by creating a simple point as the geometry. This failed with the same error. I next thought it might be an issue with ZAware. I tried setting it to both ((IZAware)cellGeometry).ZAware = true; and ((IZAware)cellGeometry).ZAware = false; with the same issue both times. I'm thinking it might have to do with how the fields in the FeatureClass were setup, but I'm quite stuck at the moment. UPDATE: I've tried updating how the fields are created to fix any issues that may exist therein, and the issue remains. The fields are generated like so,                   {             if (this.ImportedField == null)             {                 return null;             }                  IObjectClassDescription fcDesc;             if (this.DatasetType == esriDatasetType.esriDTTable)                 fcDesc = new ObjectClassDescriptionClass();             else                 fcDesc = new FeatureClassDescriptionClass();                  IFields fields = fcDesc.RequiredFields;             IFieldsEdit fieldsEdit = (IFieldsEdit)fields;                  //Setup shape field             IField shapeField = fields.Field[fields.FindField("Shape")];             IFieldEdit shapeFieldEdit = (IFieldEdit)shapeField;             shapeFieldEdit.Name_2 = "SHAPE";             shapeFieldEdit.AliasName_2 = "SHAPE";             shapeFieldEdit.Type_2 = esriFieldType.esriFieldTypeGeometry;             shapeFieldEdit.Length_2 = 0;             shapeFieldEdit.Precision_2 = 0;             shapeFieldEdit.Scale_2 = 0;             shapeFieldEdit.IsNullable_2 = true;             shapeFieldEdit.Required_2 = true;             shapeFieldEdit.Editable_2 = true;             shapeFieldEdit.Domain_2 = null;             shapeFieldEdit.DomainFixed_2 = true;             shapeFieldEdit.GeometryDef_2 = BuildGeometryDefinition(BuildGenericSpatialReference(true), esriGeometryType.esriGeometryPolygon, true, false);                  var fieldNames = this.ImportedField.CoordinateData.ScalarData.Select(nodeData => nodeData.Name).ToList();             fieldNames.AddRange(this.ImportedField.CellData.ScalarData.Select(cellData => cellData.Name));                  //Add fields from ImportedField             foreach (var fieldName in fieldNames)             {                 IFieldEdit fieldEdit = new FieldClass();                 fieldEdit.Length_2 = 1;                 fieldEdit.Name_2 = fieldName;                 fieldEdit.Type_2 = esriFieldType.esriFieldTypeSingle;                 fieldsEdit.AddField(fieldEdit);             }                  return fields as IFields;         }


title: ArcGIS for Desktop on Linux?

text: Will ArcGIS 10 work properly on Linux using **wine or crossweaver**?


title: WMS with Access Control (ACL) for different users? Or solution for thousands of points on web map?

text: I have a set of Sites in the thousands. Not all users see same the same Sites on the map - thus ACL. Data is stored in a non-GIS database (mssql, ora, etc) and right now is retrieved with a standard SQL query and returned to client (OpenLayers) with a JSP. The current code is already optimized to do client-side clustering (so less of a mess on the client, but all the data is there), bounding-box retrieval (so a sub-set of data but have to go to server for every pan/zoom for new data) and even server-side clustering (to reduce data shipped to client). What I would like to instead is generate images/tiles for all my data, similar to how Google Maps shows tons of little red dots for all the results - http://bit.ly/d73qrw [google maps search for "coffee"] and display on the client. And when clicked-on do a quick ajax call to a WMS service to get info. That's the idea at least. But here's the problem - I can't just setup a standard WMS service in front of my data because not all users see same data. Is there a way to generate these tiles on the fly or make WMS ACL-aware? Thank you, -Vadim edit - 9/22/2010 - So I found out how google generates the tiles, or rather the technology behind it. They are using Google Fusion Tables. Store N rows in their tables and then if the data is lat/lon aware, the Fusion Table product can generate the tiles on the fly (!). This is the kind of thing I'm looking for - performant map with tons of data. But of course I need it to be ACL controlled. Does it make sense to write a custom light-weight implementation of the WMS spec or modify an existing product? Though geoserver seems like a lot to "just" modify to support ACL. edit - 9/27/2010 - Some more info since adding bounty. My data is in Oracle. Ora spatial is not enabled. Right now the data is extracted at the business level and converted to data, sent to client where the client puts the "dots" on the map. ACL is done at the business logic level, not DB or ActivDir or anything like that. Authentication is simple, but Authorization is not and so had to be captured in code. Would like to know how to best create a WMS service to serve up 1000's of "dots" on the map where each user will see a different subset of dots. Is the answer a CQL_FILTER? But then how are the params set? One idea that I have now is to do a 2-step process. First run in- house query to get list of ID's that user allowed to see, then build a WMS request string with those ID's in CQL_FILTER param. Is there anything that simplifies this process? And if I go with this, how can I add this layer as a "WMS" layer to an Open Layers client since to OL the end-point is my code to get IDs from DB not the actual WMS service on e.g. GeoServer?


title: Analysing a network distance to two or more points using GRASS or QGIS

text: I am doing a piece of analysis work looking at distances along a network in GRASS. Most of the analysis involves finding which parts of the network are within 600m of a series of points. GRASS does this very well using the v.net functions. I've got stuck on trying to identify which parts of a network are within 600m of more than two points in the series. Does anyone have any ideas how I could do this within GRASS? I'm using the GRASS plugin in QGIS.


title: ArcMap 10 - Extract raster cells from imagery based on RGB value

text: I've got some aerial photography (format: IMAGINE Image) which is detailed enough to see to fire hydrants and a shapefile of fire hydrants which was digitized based on some CAD data. My project is to QC the fire hydrant shapefile locations by spot checking with the imagery. Is it possible to extract RGB values from the imagery that match the fire hydrants RGB values? Fire hydrant (red): ![Fire hydrant](http://i.stack.imgur.com/d9Uz1.png) I'm looking to select cell values (which i have written down) that are similar to the cell values of the fire hydrant and then extract those cells (either a raster file of those cells or a shapefile which plots location of similar cells). I'm assuming that there are very few features in the imagery that have similar reddish RGB Values. After looking at Arron's Answer and playing around with both the supervised and unsupervised approach (see my comment below), i was unable to get the tool to perform exactly what i wanted until I started too look at the confidence raster output from the Maximum Likelihood Classification tool. I am not to sure what exactly the raster output is supposed to symbolize but from just looking at the level 14 cell value, but it captured the all the fire hydrants. raster output: ![raster output](http://i.stack.imgur.com/5Q6B5.png) location of fire hydrants: ![location of fire hydrants](http://i.stack.imgur.com/uDhHD.jpg) My next step is to use the tool Raster to Polygon convert the fire hydrant raster footprint to a polygon. I am running the tool right now and in the past hour it has only made it to 11%. The area I am working in is large, at 1x1 mile so I understand that it could take some time if there is a lot of small raster cells to convert to a polygon. Are there any suggestions on running some tools to clean up some of the data so it speeds up the process of converting the raster to a polygon? I may not have a need to run this tool this since the raster footprint did such a good job at capturing the fire hydrants but I am interested in speeding up the process since I foresee this being used in other applications for the future.


title: Get "approximate" lat/long coordinate for UTM grid zone

text: I have data assigned with a UTM grid zone designators, such as `10S` for example. Now I need an approximate latitude/longitude pair, preferably the center point of the given grid zone. Can you give me some advice, on how to map those values?


title: Openlayes - Loading 3 different KMLs via proxy; one appears and the others don't

text: I'm loading 3 KMLs using an ashx handler taken from http://code.google.com/p/iisproxy. The 3 KMLs I'm trying to load are Openlayers sundials ( http://openlayers.org/dev/examples/kml/sundials.kml ), Marine Traffic KML ( http://www.marinetraffic.com/ais/ge_marinetraffic.kml ), and flightwise ( http://mw1.google.com/mw-earth- vectordb/flightwise/flightwise_root.kml ). The Sundial is working for me, but the other 2 don't. I'm afraid the problem is because the other 2 KMLs have file links inside while sundials KML has data only. Is there a way I can make the other two work on my IIS? This is the code I used:               <!-- language: lang-js -->          function init() {              var map = new OpenLayers.Map("map");         OpenLayers.ProxyHost = "/proxy.ashx?url=";              var osm = new OpenLayers.Layer.OSM();              var flight = new OpenLayers.Layer.Vector("Flightwise", {             strategies: [new OpenLayers.Strategy.Fixed()],             protocol: new OpenLayers.Protocol.HTTP({                 url: "http://mw1.google.com/mw-earth-vectordb/flightwise/flightwise_root.kml",                 format: new OpenLayers.Format.KML({                     extractStyles: true,                     extractAttributes: true                 })             })         })              var marine = new OpenLayers.Layer.Vector("Marine", {             strategies: [new OpenLayers.Strategy.Fixed()],             protocol: new OpenLayers.Protocol.HTTP({                 url: "http://www.marinetraffic.com/ais/ge_marinetraffic.kml",                 format: new OpenLayers.Format.KML({                     extractStyles: true,                     extractAttributes: true                 })             })         })              var sundials = new OpenLayers.Layer.Vector("Sundials", {             strategies: [new OpenLayers.Strategy.Fixed()],             protocol: new OpenLayers.Protocol.HTTP({                 url: "http://openlayers.org/dev/examples/kml/sundials.kml",                 format: new OpenLayers.Format.KML({                     extractStyles: true,                     extractAttributes: true                 })             })         })              map.addLayers([osm, flight, marine, sundials]);              var controls = [                         new OpenLayers.Control.Navigation(),                                                 new OpenLayers.Control.LayerSwitcher(),                                                 new OpenLayers.Control.ScaleLine(),                                                 new OpenLayers.Control.MousePosition(),                         /*new OpenLayers.Control.OverviewMap(),*/                         new OpenLayers.Control.KeyboardDefaults()                     ];              map.addControls(controls);                   map.zoomToMaxExtent();     }


title: ArcGIS how can I limit study area?

text: I'm ArcGIS beginner. At the moment I have 2 layers:   * roads for whole country   * zone of the city I would like to delete the roads that are outside the city. Any idea how can I do it?


title: Creating rasters?

text: I am doing my masters in archaeology, but have used GIS for some time now, however this is my first attempt at generating and using rasters. I appear to have been successful for some, but not so well for others. All the books and tutorials I have looked at don't really explain what the different values are when creating a raster and I'm unsure if I should be changing the default values to something else to produce the image I want and if I am to change them what to? I could spend a life time changing the values, but still wouldn't be sure I understood what they did. Can anyone recommend a book, blog, tutorial that would help me with this? I am using GRASS, QGIS and ArcMap/Scene, but if there is another program out there that could be useful I'm willing to try. Thanks Christina


title: Generate raster from isolines

text: I have a vector dataset of mean annual precipitation polylines, and I would like to generate a raster from them. I'm not sure the best way to go about it. ArcGIS has the **Polyline to Raster** tool but from reading its help file I get the impression it's not really designed to do what I need it to do. Obviously some kind of interpolation is required. I am hoping that SAGA-GIS might have a useful module. Either that or I guess I could convert the lines to points and then krige them. I would appreciate learning from others' experiences on this one.


title: Call GeoExt action by click ExtJs button

text: I have a problem.   I have a ExtJs button and want to call anyone GeoExt action. For example this               {xtype:'tbbutton',                             cls: 'x-btn-icon',                             id: 'testbutton',                             tooltip: "Рисовать",                             icon: url_servlet+'externals/gxp/src/theme/img/silk/delete.png',                             handler: function(){                                 new GeoExt.Action({                                     text: "max extent",                                     control: new OpenLayers.Control.ZoomToMaxExtent(),                                     map: app.mapPanel.map                                 })                             }                                }      But its not work and map not zoom to max extent. So where i can mistake. In ideal i want to call action from actionColumn's buttons/ I will glad if u help/


title: Surface covered by buffers in QGIS

text: In QGIS, I have plotted buffers around points. I would like to calculate the area which falls under buffers without over counting due to buffer overlaps. Is there a quick and easy way to do this?


title: file is not a supported raster data source

text: I'm using Ubuntu 12.04 and QGIS 2.01. I'm receiving the error message "filexx is not a supported raster data source". filexx is a .RST Idrisi format and has 1.3 Gb in my notebook disk (Ext4 format) and copied on a 4Tb external disk (NTFS format). Is there a limitation to manage the size of the file? I was able to read many .RST smaller files. Is there a problem with the disk format?


title: ArcGIS Viewer For flex add remove layer programitically

text: I want to add remove operational layers to ArcGIS Flex Viewer Problematically with some security verification. if user is Manager all layers will be added and visible if user from public then due to security reseason i need to remove some layers. how i can handle it in arcgis flex viewer. if some relevant tutorials exits then please refer me. or guide me .


title: convert wgs84 coordinate to Google Map coordinate system(EPSG 3785) with geotools and java

text: Could you please point me out where can I make such a transformation: convert wgs84(EPSG 4326) coordinate to Google Map coordinate system(EPSG 3785), and vise versa. Thank you.


title: How to split polygons without ArcView?

text: I'm trying to split polygons into equal areas. While fishnet sort of works; the areas I'm trying to split don't touch; they're national parks. Fishnet spans the whole area, not just the polygons in question. Ultimately, I'm searching for a method to split polygons into a user-given number of equal area polygons. I've read everything I can about the Huber methods to do this in ArcView, but without access to that (or experience with it) this isn't the answer for me. Thanks very much.


title: Append to text file using ArcGIS and Python

text: I have wrote a small routine using Python for ArcGIS (9.3) that will get the count of 2 data layers, calc a percentage, and then write the text back to a text file (CSV) to create a sort of "report" that can be opened in excel and formatted. I want to be able to run this routine on many different dataset pairs and append the data to the end of my CSV report. The problem is that I am having is that when I attach this Python script to an ArcGIS Script Tool, it automatically erases the file. I am testing that with               if os.path.exists(file):      right off the bat and this response is always "file does not exist". I have tried setting               gp.OverwriteOutput = 0      with no success, and am appending to the file using               fwrite = open(file, 'a')      Does anyone have a workaround that will allow me to choose an output file but NOT overwrite the output?


title: Belgium Flood And Precipitation Mapping

text: Similar to this question, does any one know a source for free downloadable flood and precipitation maps for Belgium? Thanks all in advance. C


title: how to make openlayers example work in jsFiddle

text: Hy, I'm working on a jsFiddle with a openlayers example in it. http://dev.openlayers.org/releases/OpenLayers-2.11/examples/snapping.html At the moment it's not working because it's not getting a response for the http request to get data. How do I fix that? The jsFiddle is here: http://jsfiddle.net/TcuxA/6/ Go to the line "// create three vector layers" in the script. There are 3 requests for data. If you type the URLs in your browser you get the JSON, but my firebug gives 3 errors when I run the jsFiddle. I tried fixing with jsFiddle echo, but that didn't work. I don't know how to change the script to load the data otherwise. Why can't I get the json by these URLs? What are good sollutions?


title: GDAL and ArcObjects in The Same Project - A Recipe For Disaster?

text: I am working on a project where I am attempting to use GDAL and ArcObjects together in the same project. Although the project itself is nothing more than a toolbar in ArcMap, I would like to use some GDAL based raster libraries developed for other non-ArcMap projects that offer a SIGNIFICANT performance boost over the equivalent ArcObject based methods. I've been getting all types of crazy errors when mixing the two. However, when I execute the same chunks of code in localized (either GDAL alone or ArcObjects alone) projects, everything works as expected. Can ArcObjects and GDAL co-exist in the same project? Any tips to avoid complications? FYI - I am using the C# GDAL bindings distributed with FWTools 2.4.7.


title: Concatenate Fields into One Field in QGIS

text: I have seen some posts on this, but only for ArcMap. I am looking to take a range of values from two fields and concatenate them into one field to make a unique value. For example:               XMAX     YMAX       GRIDID     332500   5366000    33255366     333000   5366000    33305366      In ArcMap you can do something like `(Mid([XMAX]),1,4)+(Mid([YMAX]),1,4)`. I know this can be done in a PostGIS database, but for the moment it needs to be done on a shapefile and I am looking to use QGIS. Thanks


title: ArcGIS JavaScript API: Override infowindow 'Zoom To' function or fire a second event

text:
Those of you familiar with the Popup version of InfoWindow in ESRI's JavaScript API know that the Popup has a 'Zoom To' link on it that zooms you 4 levels closer to the linked feature. I'd like to fire another event after the default zoom to function executes. Here is the link in the HTML.

         Zoom to

      I'm able to add a function to the Zoom To click with jQuery like this:               $('.action.zoomTo').click(function() {alert('Clicked!');});      But the the function above occurs _before_ the default function of 'Zoom To'. I wouldn't mind their click function executing, so long as it fires _after_ mine does. Any ideas?


title: How to clean up contours?

text: I have .dxf file with contour lines representing elevation. But this contours are too "ugly" and elevation model looks just bad. As you can see on image, contours are not continuous... ![Contours and DEM created from them](http://i.stack.imgur.com/ORdLb.png) Is there any way to clean up my contour lines so generated DEM will looks as it should? So I tried terrain over TIN method and here is result. Unfortunately that pesky spots are still there. This is created from whole contour file, while first image is just cropped part of it ![This is created from whole contour file, while first image is just part of it](http://i.stack.imgur.com/cFOoE.png)


title: Get ESRI JSON response locally from Arc document?

text: I'm trying to use the ArcGIS Flex SDK for a mobile app that will sometimes be used in environments with no connectivity. I've figured out that layers can be brought in locally by saving the JSON response given by the server and simply using the local file as the layer URL. Is there a toolbox, plugin, or .py script that will take existing features or shapefiles and convert them to the ESRI JSON format that the server uses? I've tried esri2open, but that only does GeoJSON and a generic JSON object, and I'd prefer to have something that automatically gives the exact same JSON structure as the REST services would when loading an online map. Ideally, it would be something that could run from the Arc desktop without the additional steps of publishing a service and downloading the response. Does something like this exist? I've noticed that 10.2 may have that capability, but for now I'm stuck with 10.1.


title: grass - add on problem

text: I can't install the addon in grass 6.4.3 to my linux mint. i've already installed the SVN and than it says: GRASS 6.4.3 (fabian_idrisi):~ > g.extension r.fuzzy Fetching from GRASS-Addons SVN (be patient)... A r.fuzzy/local_proto.h A r.fuzzy/shape.png A r.fuzzy/main.c A r.fuzzy/boundary.png A r.fuzzy/description.html A r.fuzzy/fuzzy.c A r.fuzzy/Makefile A r.fuzzy/set.png U r.fuzzy Checked out revision 60053. Compiling ... Makefile:13: /usr/lib/grass64/include/Make/Module.make: No such file or directory make: _*_ No rule to make target `/usr/lib/grass64/include/Make/Module.make'. Stop. ERROR: Compilation failed, sorry. Please check above error messages. I try to locate that directory but my usr/lib/grass64/ folder only contain scripts, bin, driver,lib licale ect. What am i doing wrong? Thanks


title: Shapely WKT precision/rounding

text: I'm using Shapely's wkt.dumps function to get WKT from geometry objects. I can't seem to find a way to limit the precision of the coordinates. Any help is appreciated.


title: PostGIS automatically creating geometries after importing using QGIS DB Manager

text: I have a shapefile consisting of linestrings, in which some of the features have null geometry. Yet, after importing it to PostGIS, using QGIS' DB Manager, a geometry is created for those features (what it does is create a small segment on top of another existing linestring). Has anyone ever seen this happen? This problem doesn't happen when I import the shapefile using OpenGeo's pgShapeLoader. Is it a bug of DB Manager, or am I not configuring it correctly for importing? I'm using QGIS 2.2.0 for mac, and PostGIS 10.1.0


title: Will ArcGIS 10 support Python 2.7 or Python 3.1?

text: I'm curious to know whether ArcGIS 10 will support newer versions of Python than 2.6. And for that matter newer versions of NumPy (ie, 1.4, 1.5 etc). I'm sure that this documentation is on the ESRI site somewhere, but I can't find it. I have found that ArcGIS 10 will respect an existing Python 2.6 install as long as you have Numpy 1.3.0 installed (ref), and I'm guessing, at least until Numpy 1.5.0 is out and supporting Python 3, then no Python 3 for ArcGIS, but I'd still be interested to know if upgrading is possible when it's available.


title: How to open an output set that is a list of double values in a geoprocessing service

text: we are trying to reduce the amount of data we are returning from the server by eliminating returning a polyline and instead returning a collection of values of type GPDouble. When the transaction is kicked off inside our GeoProcessing Service via Execute, the first thing we do is create the parameters we are either going to process or return. for example to get the target list we execute this code:                       IFeatureClass targetClass;             filter = null;             _gputilities.DecodeFeatureLayer(GPUtilities.GetParameterValueByName(parameterValues, "TargetList"), out targetClass, out filter);             if (targetClass == null || targetClass.FeatureCount(null) == 0)             {                 message.AddError(2, "Could not open Target input dataset.");                 return;             }             message.AddMessage("Target List has been defined.");      I changed my OutList parameter definition as follows (old definition commented out)                           //IGPParameterEdit3 outputPoints =                 //    GPUtilities.CreatePolylineRecordsetParameterEdit("OutputList", "Output List", "MsramSDE.DBO.Standard_Polyline",                 //        esriGPParameterDirection.esriGPParameterDirectionOutput,                 //        esriGPParameterType.esriGPParameterTypeDerived, new GPFeatureRecordSetLayerTypeClass(), true);                 IGPParameterEdit3 outputPoints = GPUtilities.CreateParameterEdit("OutputList", "Output List",                      esriGPParameterDirection.esriGPParameterDirectionOutput,                     esriGPParameterType.esriGPParameterTypeDerived, new GPDoubleTypeClass(), true, "0.00");                 parameters.Add(outputPoints);      None of our other GeoProcessing Services return anything but a FeatureLayer. So I lack any example and could not find one via Google for opening up an Output Item that gives me a way of writing out the double values. Anyone ever done this and how do I do it? I've looked through the GPUtilities calls to see if there was an obvious DecodeBaseType or some such call but could not find one. Should this be defined as a record set instead?