title: Problem loading OS OpenData multipoint shapefile into PostGIS

text: I am trying to load the OS OpenData Strategi shapefile dataset into a PostGIS database. I have had no problems with the polyline and polygon shapefiles, but I am unable to get the point shapefiles to load. I am using the following shp2pgsql command line to create the table and load the data:               shp2pgsql -c -I -s 27700 admin_font_point strategi_point | psql -d opendata      The error I get is:               Shapefile type: MultiPoint     Postgis type: MULTIPOINT[2]     SET     SET     BEGIN     NOTICE:  CREATE TABLE will create implicit sequence "strategi_point_gid_seq" for serial column "strategi_point.gid"     NOTICE:  CREATE TABLE / PRIMARY KEY will create implicit index "strategi_point_pkey" for table "strategi_point"     CREATE TABLE                              addgeometrycolumn                              -------------------------------------------------------------------     public.strategi_point.the_geom SRID:27700 TYPE:MULTIPOINT DIMS:2      (1 row)          ERROR:  new row for relation "strategi_point" violates check constraint "enforce_geotype_the_geom"     ERROR:  current transaction is aborted, commands ignored until end of transaction block      The last error just repeats for each insert. The generated SQL looks like:               SET CLIENT_ENCODING TO UTF8;     SET STANDARD_CONFORMING_STRINGS TO ON;     BEGIN;     CREATE TABLE "strategi_point" (gid serial PRIMARY KEY,     "code" int4,     "legend" varchar(42),     "file_name" varchar(16),     "number" numeric(11,0),     "name" varchar(180),     "number0" varchar(64),     "admin_name" varchar(50),     "type" varchar(40),     "ferry_from" varchar(50),     "ferry_to" varchar(50),     "ferry_time" varchar(10),     "ferry_type" varchar(20),     "restrictio" varchar(20),     "access" varchar(30),     "amended" date,     "usage" varchar(64),     "location" varchar(30),     "gis" varchar(80),     "owner" varchar(60),     "north" varchar(60),     "south" varchar(60),     "east" varchar(60),     "west" varchar(60),     "clockwise" varchar(60),     "anticlockw" varchar(60),     "imperial" int4,     "metric" int4);     SELECT AddGeometryColumn('','strategi_point','the_geom','27700','MULTIPOINT',2);     INSERT INTO "strategi_point" ("code","legend","file_name","number","name","number0","admin_name","type","ferry_from","ferry_to","ferry_time","ferry_type","restrictio","access","amended","usage","location","gis","owner","north","south","east","west","clockwise","anticlockw","imperial","metric",the_geom) VALUES ('5734','Administrative Detached Attribute Point','gb_north','5307',NULL,NULL,'ISLE OF MAN',NULL,NULL,NULL,NULL,NULL,NULL,NULL,'20000413',NULL,NULL,NULL,NULL,NULL,NULL,NULL,NULL,NULL,NULL,'-9999','-9999','0101000020346C000000000000302A0A4100000000304F1C41');      and then lots more inserts all looking much the same. I have tried all sorts of variations; with and without SRID; with and without .prj file; created a new test database with no other data; using simple geometries; but all to no effect. I have made zip file of shapefile available if anyone has the time to try this themselves: http://www.passback.org.uk/tmp/admin_font_point.zip I am using a fully up to date Fedora 14, Postgres 8.4.6, PostGIS 1.5.1. Some updated notes based on the suggestion to use SPIT: If I use shp2pgsql in "create table only" mode I get a table with the following constraints:               Check constraints:         "enforce_dims_the_geom" CHECK (st_ndims(the_geom) = 2)         "enforce_geotype_the_geom" CHECK (geometrytype(the_geom) = 'MULTIPOINT'::text OR the_geom IS NULL)         "enforce_srid_the_geom" CHECK (st_srid(the_geom) = 27700)      Looking at the table created by SPIT I get slightly different constraints:               Check constraints:         "enforce_dims_the_geom" CHECK (st_ndims(the_geom) = 2)         "enforce_srid_the_geom" CHECK (st_srid(the_geom) = 27700)      I guess I could use shp2pgsql -p to create the table and then use psql to remove the enforce_geotype_the_geom constraint and my bulk import from the command line using shp2pgsql -a should work. I would still like to understand why that constraint is causing a problem? Thanks, Keith.


title: Efficient way to provide cached tiles?

text: I am creating an app in Android that requires cached areas. I allow the user to cache an area of reasonable size. The server provides these tiles in a zip in the format /z/x/y.png. I took a (very) naive approach and simply wrote a python script to grab the images from the localhost using apache2 with mod_tile (127.0.0.1/osm/z/x/y.png), save to a directory, then zip. Under a load test, this failed terribly. Is there a more efficient way to create the png files locally? (Mapnik api and python?) Or is a queue for users an absolute necessity to prevent server overload?


I believe I know the answer, but I wanted to ask. Is there a way to create a Personal Geodatabase in QGIS? I know QGIS can view Personal Geodatabases, but can one be created? I have a project where the client is requiring a small spatial database (they are very low on the tech level, so I do not want to jump into PostGIS etc etc). They also only use ArcGIS. I would like to do the project in QGIS, hense my question. Thanks in advance....


title: multiple export in tiff using "export map..."

text: Is there a way to export multiple maps in *.tiff using "Export Map" and setting Data Driven Pages? Because Data Drive Page work in "Export map" only with PDF.


title: Logging geoserver's request logs into postgres

text: I'm trying to force geoserver to write request.log(s) into a postgis database. Unfortunately, geoserver's documentation remains silent about logging into database at all. Searching the web I found one answer on mailing list for similar question. I didn't try this, but that's because I've got feeling that's not thing I'm looking for, because this solution seems to allow to save to database logs of geoserver's activity not incoming request what is my point... I also found geoserver's plugin called monitor, which could be thing I'm looking for, (again) unfortunately in v 2.3 adding it brakes answers for Request=GetFeatureInfo, so it's useless for me. Finally, there's always a solution to use PostgreSQL COPY to import data from logfiles into database, but I think there have to be more elegant way. The point is to log requests into database, so after little post-processing on request strings, I'll be able to write simple reports of IPs using my server or most commonly downloaded areas. More specified question is: where can I find settings for request logs and how to reconfigure it to put data into database? Can anybody help?


title: Points to Cells

text: ![Cells with assigned values to it](http://i.stack.imgur.com/xXlWz.jpg) How do I create such map from point data? I am using Arc GIS 9.3 and 10


title: Triggers for delta tables

text: I have been facing an issue to capture the delta updates on SDE. We have written triggers on A and D tables to capture the delta updates(to capture the version name based on state id's), however I couldn't able to get the very first delete of any feature from a new version that has been created on Default. After debugging, I have observed that, after the very first delete of a feature, the State Id in versions table is not being updated however the State Id in SDE_States is updated. So, the state id in versions table and the state id in SDE_States table is different. I would like to say that the same query is working fine from second deletes and so on.... Having issue with getting the version name for only the first delete. For your reference, please do find the query below.               Select @ObjectId=SDE_DELETES_ROW_ID, @SdeStateId=DELETED_AT, @SdeCurrentStateId = SDE_STATE_ID FROM inserted;          SELECT @Version= Name FROM sde.SDE_versions v,sde.SDE_states s INNER LOOP JOIN          sde.SDE_state_lineages l ON l.lineage_name = s.lineage_name WHERE s.state_id = v.state_id and lineage_id = @SdeStateId and s.state_id = @SdeStateId          IF @version IS NULL     SELECT @Version= Name FROM sde.SDE_versions v,sde.SDE_states s INNER LOOP JOIN          sde.SDE_state_lineages l ON l.lineage_name = s.lineage_name WHERE s.state_id = v.state_id and lineage_id = @SdeCurrentStateId


title: Working with GB SRTM Digital Elevation Model (DEM) 90m

text: I have downloaded the **freely available DEM file** from sharegeo open. I would like to know for a set of points at what elevation they lie. I have imported the DEM file into ArcGIS but it seems strange; the highest point is 63. I was wondering if anyone else has experience with this dataset or similar and could shed some light.


title: How to save KML-layers style (Q-GIS)?

text: Is it possible to save a KML-layers style? My layer is set to 50 % opacity and I want to save that setting into the KML-file. The plugin qgistokml can´t yet do it, and only work with points.


title: How to eliminate small features in a grid?

text: (Update) My final goal is to take some logged point data (1 point about every 20 ft or so there are over 15,000 points) and interpolate that point data. Then I would like to convert the interpolated raster into distinct polygon zones based on the number of classes the data is in using the Natural Breaks classifications. So if I want it in 5 classes, there will be several polygons that make up these 5 classes, but within the tabular data I can see which of the 5 classes each polygon belongs to. In addition to this, I would like to eliminate the smaller polygons created by the interpolation, like say anything under 1 acre in size. Maybe there is a GP process I can run to do this? Let me know if I need to explain this better, and thanks for the help so far. * * * I was thinking the "output cell size" setting would be the trick here but it defaults to a very small number (0.00003) and if I make it any bigger like even say "0.003" it results in the interpolated raster having virtually no data where both high and low values is a really high number and very low number. If I use the default settings I get a decent kriged raster but I would like to adjust to exclude smaller areas, like say less than 1 acre in size. Basically I want the kriging interpolation to ignore areas that are smaller than an acre or even a 1/2 acre. Currently it is showing areas that are easily less then 0.001 acres in size of which I would like it to ignore those. Also if I zoom into it looks like the cell size grids are 8ft by 8ft. How do I get that larger to say 100 ft by 100 ft? Like I said I have tired putting larger number in the output cell size but get back basically nothing but extreme values. I must be doing something wrong. Any help would be appreciated.


title: Projection parameters from WGS84 UTM37 to Ain El Abd UTM37 and the reverse

text: I am developing application using JavaScript apis, I want the correct projection parameters to convert from WGS84 UTM37 to Ain El Abd UTM37 and the reverse. When I click in the map I get Ain El Abd Coords from the map onClick event. I need to convert these to WGS84 for display in the Web App.


title: sub-dividing KML files

text: I do have a kmz file which marks the boundaries of a larger district and it's sub-districts. Still new to GIS and in the process of getting familiar with QGIS, I wonder can anyone advise me how can I "extract" the boundaries of each of these sub-districts from that kmz? E.g. are there any QGIS plugins which could be used for this task? Thanks !


title: How to set file path to the modelbuilder tool?

text: I have got the following problem: I want to dissolve all the shapefiles in my directory and save them under a different name (consisting of the same name and an addition "_diss"). I have managed to create the full path to the new file, but I dont know how to put this as an argument to the dissolve tool. I have tried "%Full path%" but it says that it is not able to save it. ![enter image description here](http://i.stack.imgur.com/UzDXk.png) ![enter image description here](http://i.stack.imgur.com/tSrqI.png) ![enter image description here](http://i.stack.imgur.com/lNYxh.png) **How can I save the shapefiles with the names I have specified?**


title: Arcpy: replaceDataSource, cad_workspace

text: I am having trouble replacing the data source on CAD files using arcpy. Everything I have tried runs into this error:               Runtime error type 'exceptions.ValueError' : Layer: Unexpected error      Here is a list of arguments I have tried:               layer.replaceDataSource('C:/CAD_DIR/CAD_File_1.dxf', 'CAD_WORKSPACE', 'Polyline')          layer.replaceDataSource('C:/CAD_DIR/CAD_File_1', 'CAD_WORKSPACE', 'Polyline')          layer.replaceDataSource('C:/CAD_DIR/', 'CAD_WORKSPACE', 'CAD_File_1/Polyline')          layer.replaceDataSource('C:/CAD_DIR/', 'CAD_WORKSPACE', 'CAD_File_1.dxf/Polyline')          layer.replaceDataSource('C:/CAD_DIR/', 'CAD_WORKSPACE', 'CAD_File_1.dxf.Polyline')          layer.replaceDataSource('C:/CAD_DIR/', 'CAD_WORKSPACE', 'CAD_File_1.Polyline')      It appears to take more time to think about going wrong when I put the 'CAD_File_1*' part in the 'Dataset_Name' section, that is to the right of 'CAD_WORKSPACE'. But it still fails. I have noticed that it treats the CAD file as a Feature Dataset and the particular feature type, in this case Polyline, as the Feature Class which may be the issue. I have also tried this when pointing at a 'Polyline' layer:               layer.replaceDataSource('C:/CAD_DIR/', 'CAD_WORKSPACE', 'CAD_File_1.dxf)      It doesn't produce the same error or in fact any error but the data does not update.


title: How to solve a Run-time error 91?

text: I have a function which a part shown below. The problem is that I get a Run- time error 91 (object variable or with variable not set) and the debug highlights the last line of the code ( Set pFeatureClass = pFeatureLayer.FeatureClass) shown below while the layer file exist on the table of contents. The strange thing is that when I use this function by giving directly the arguments as shown in the example with asterisks below it works fine but when the function is used within a loop in a procedure it presents that error 91. Could you help me please?               Public Function GetAveX(BlockID As String, Solindex1 As String, Solindex2 As String _     , MPindex1 As String, MPindex2 As String, ParcelID As Long) As Double          'Get the current map document     Dim pMxDoc As IMxDocument     Set pMxDoc = ThisDocument          'Get the active map (data frame)     Dim pMap As IMap     Set pMap = pMxDoc.FocusMap          ' Get the proper layer      Dim pLayers As IEnumLayer     Set pLayers = pMap.Layers               Dim pLayer As ILayer     Set pLayer = pLayers.Next     Do Until pLayer Is Nothing     If pLayer.Name = "VectorB" & BlockID & "MP" & "_" & Solindex1 & "_" & MPindex1 Then         Exit Do     End If     Set pLayer = pLayers.Next     Loop          'Get the feature class     Dim pFeatureLayer As IFeatureLayer     Set pFeatureLayer = pLayer          Dim pFeatureClass As IFeatureClass     Set pFeatureClass = pFeatureLayer.FeatureClass          'rest of the code     ....               '**********An example that works fine     Dim i As Double     i = GetAveX(25, 27, 36, 19, 20, 45)     MsgBox i


title: How to improve accuracy or extrapolate missing or bad GPS points

text: Is there an easy to use freeware tool to postprocess my location data, in order to correct 1) location errors (or improve accuracy not based on differential GPS correction) and 2) deal with missing locations, based on ecological information? e.g. spikes in the location data are indicative of location errors, as are large angles between 3 fixes at small step sizes Brown hare location data is based on small GPS collar (U-blox LEA4S), and has been measured every 6 minutes in an open , relatively flat enviroment, including acceleration data every minute for 35 days, 24hr a day


title: Do ArcGIS SpatialReference object factory codes correspond with EPSG numbers?

text: I've got a project where I'm reading in geometry objects from text files to put them in ArcSDE. These geometry objects have an associated EPSG code, and I'm looking for a way to ensure they're in the right projection when I've finished with them. Having taken a look at the factory codes that I'm interested in, (GCS GDA94; EPSG:4283, and VicGrid 94; EPSG:3111) for a project I'm working on it appears that the EPSG codes correspond directly to the ESRI factory codes for the spatial reference (here and here). This is great, as it means I don't have to worry about navigating to projection files to create SpatialReference objects, for example in the code below:               import arcpy          # Create a spatial reference object using a factory code     #     sr = arcpy.SpatialReference()     sr.factoryCode = int(epsg_code.lstrip("EPSG:"))     sr.create()      Unfortunately I can't find any reference in the ESRI documentation to say whether this holds true all the time or not - I'm guess it does, but I'd rather know for sure. Does anyone know where I'd be able to find out?


title: No PopUp with WMSGetFeatureInfo and Mapserver

text: I've searched for answer for the last couple of days with no success. When I click on an object on the screen I get this reply from the server               http://localhost/cgi-bin/mapserv.exe?LAYERS=testpoint&QUERY_LAYERS=testpoint&STYLES=&SERVICE=WMS&VERSION=1.1.1&REQUEST=GetFeatureInfo&BBOX=517948%2C157505%2C535448%2C170005&FEATURE_COUNT=1&HEIGHT=500&WIDTH=700&FORMAT=image%2Fpng&INFO_FORMAT=text%2Fhtml&SRS=EPSG%3A27700&X=357&Y=241&MAP=%2FWorkTemp%2Ftest.map      following this link returns the value for the object I just picked. I get no errors in the mapserver or openlayers log file. OpenLayers code               layer = new OpenLayers.Layer.WMS( "WMS data","http://localhost/cgi-bin/mapserv.exe?", {layers: 'testpoint' , map: '/WorkTemp/test.map', transparent: true},{isBaseLayer: false, opacity: 1, singleTile: true, visibility: true});          osMap.addLayer(layer);          info = new OpenLayers.Control.WMSGetFeatureInfo({            url: 'http://localhost/cgi-bin/mapserv.exe?',            title: 'Identify features by clicking',            queryVisible: true,            maxFeatures: 1,            vendorParams: {                 map: '/WorkTemp/test.map'            } ,            info_format: 'text/html',            layers: [layer],            eventListeners: {            getfeatureinfo: function(event) {                  if(event.features.length > 0) {                          osMap.addPopup(new OpenLayers.Popup.FramedCloud(                          "chicken",                           osMap.getLonLatFromPixel(event.xy),                          null,                          event.text,                          null,                          true                          ));                  }            }            }     });          osMap.addControl(info);     info.activate();      test.map layer               LAYER       NAME "testpoint"       TYPE POINT       STATUS ON       CONNECTIONTYPE OGR       CONNECTION "MYSQL:test,host=localhost,user=root,password=xxxxxx,port=3306"       DATA "SELECT SHAPE,test from testpoint"       DUMP TRUE              TEMPLATE "c:\WorkTemp\test_template.html"                 LABELITEM "test"        PROJECTION             "init=epsg:27700"         END       CLASS         NAME "test"         STYLE           SYMBOL  'circle'           SIZE 10           COLOR 255 215 0         END         LABEL          COLOR  0 0 0          FONT arial          TYPE truetype          SIZE 8          POSITION AUTO          PARTIALS FALSE          OUTLINECOLOR 255 255 255         END       END     END # layer      This is the same problem that I'm having :- Use OpenLayers WMSGetFeatureInfo on Mapserver Layer?. But I have already set my WMS layer to be OpenLayers.Layer.WMS I downloaded and installed OSGEO4W as my starting point. I'm calling http://openspace.ordnancesurvey.co.uk/osmapapi/openspace.js for my openlayer commands. any suggestions would be gratefully appreciated. Char


title: How to release locks to temp files after ArcPy script crashes in ArcMap?

text: Sometimes when I'm testing an ArcPy script and it ends with exception, ArcMap is still holding locks to some temp files (I run my scripts through ArcMap). I cannot run the script any more, because it cannot delete/overwrite the temp files and in the end I have to restart ArcMap to release the locks. Is there a way to do in a more clever fashion?


title: Create point in intersections of two layers

text: I have two layers. One of them is 'Road_line' and another one is 'forest_shape'. I want to create a point layer that has points in intersection of that layers. I couldn't find any tools for it. Could you please help me Sincerely yours


title: Define gdalwarp tile size

text: gdalwarp has `-co TILES=YES` (256x256 as default) Can I set different tile size? Thanks


title: How to convert from GeoTIFF to PNG/JPG and Web Mercator in Java?

text: I'd like to convert raster image from (geo)tiff to an png/jpg and afterwards if it is necessary ... I'd like to convert the coordinates with help of EPSG- Codes to the OSM/Google Projection(EPSG: 3857/ 900913). All this i'm doing with java programming. Anyone a good hint... like an java libary?? Thanks a lot.


title: Need help with Merge GeoProcess task in ArcEngine (C#)

text: Help! Having the hardest time getting the correct syntax when adding layers as gp parameters to call the Merge_management GeoProcessing Task. Below is just one of MANY, MANY ways I have tried with no luck. Can anyone help? I've successfully called many GeoProcessing tasks via code in ArcEngine, but this is the first one to take an array ( ie. [layer1,layer2,etc]) as one of the parameters. Intial Attemp:                   public void RunMergeGeoProcess(ILayer SourceLayer1, ILayer SourceLayer2, string ResultFeatureClassName, string GeoDBAbsFileName)         {                  IFeatureLayer featureLayerSource1 = (IFeatureLayer)SourceLayer1;             IFeatureLayer featureLayerSource2 = (IFeatureLayer)SourceLayer2;                       _geoProcessor.SetEnvironmentValue("Workspace", GeoDBAbsFileName);             //_gpParams.Add(featureLayerSource1);             ////_gpParams.Add(featureLayerSource2);                  IWorkspace refWS = ((IDataset)featureLayerSource1.FeatureClass).Workspace;             string dbPath = refWS.PathName;                  IDataset myDS1 = featureLayerSource1.FeatureClass as IDataset;             string mysDS1 = dbPath + "\\" + myDS1.Name;             IDataset myDS2 = featureLayerSource2.FeatureClass as IDataset;             string mysDS2 = dbPath + "\\"+ myDS2.Name;                  string myResult = dbPath + "\\" + ResultFeatureClassName;                  string param1 = "[" + mysDS1 + ", " + mysDS2 + "]";                  _gpParams.Add(param1);             _gpParams.Add(myResult);                       try             {                 _geoProcessor.OverwriteOutput = true;                 _geoProcessor.Execute("Merge_management", _gpParams, null);                 //returnMessages(_geoProcessor);             }             catch (Exception ex)             {                 MessageBox.Show("There was a GeoProcessing Error." + ex.ToString());                 returnMessages(_geoProcessor);             }                  _gpParams.RemoveAll();              }      Working Version:                   public void RunMergeGeoProcess(ILayer SourceLayer1, ILayer SourceLayer2, string ResultFeatureClassName, string GeoDBAbsFileName)         {                  IFeatureLayer featureLayerSource1 = (IFeatureLayer)SourceLayer1;             IFeatureLayer featureLayerSource2 = (IFeatureLayer)SourceLayer2;                  _geoProcessor.SetEnvironmentValue("Workspace", GeoDBAbsFileName);                  IWorkspace refWS = ((IDataset)featureLayerSource1.FeatureClass).Workspace;             string dbPath = refWS.PathName;                  IDataset myDS1 = featureLayerSource1.FeatureClass as IDataset;             string mysDS1 = dbPath + "\\" + myDS1.Name;             IDataset myDS2 = featureLayerSource2.FeatureClass as IDataset;             string mysDS2 = dbPath + "\\" + myDS2.Name;                  string myResult = dbPath + "\\" + ResultFeatureClassName;             string param1 =  mysDS2 + ";" + mysDS1  ;             _gpParams.Add(param1);             _gpParams.Add(myResult);                  try             {                 _geoProcessor.OverwriteOutput = true;                 _geoProcessor.Execute("Merge_management", _gpParams, null);                 //returnMessages(_geoProcessor);             }             catch (Exception ex)             {                 MessageBox.Show("There was a GeoProcessing Error during the Merge Operation." + ex.ToString());                 //returnMessages(_geoProcessor);             }                  _gpParams.RemoveAll();              }


title: Exporting table to XYZ ASCII file via ArcPy?

text: I'm looking for a way to export an ArcGIS table (created with the **Sample tool** ) to a text file via ArcPy. I can do this in ArcGIS via the context menu by right-clicking the table, but have not found a way to script this.


title: ArcObjects 10 Java SDK

text: I am being tasked with converting a AGS georpocessing toolkit currently written in AO 10 C# SDK to AO 10 Java SDK. Are there any resources out there which show an end to end example of creating an AGS 10 Java deployable toolkit? Thanks in advance


title: Projection for polygon intersection confined into one hemisphere?

text: I have polygons where each vertex is a (lon,lat) geographic coordinate. I want to find out if they intersect or contain each other. The polygons are guaranteed to never exceed a hemisphere (so they'll never 'wrap around' beyond a plane that cuts the Earth in half), but they can vary dramatically in scale. I think one way this might be done is by transforming the polygons using some specific projection, doing the math for intersection in cartesian coordinates and then 'unprojecting' back to geographic coordinates. Is there a projection I can generate that is fast to apply, even if its at the cost of somewhat not-so-accurate (~50m-100m) results? I'm not using any database tools or large libs here, just C++ as I need a small, portable set of code. It would be handy if I could come up with the projection on the fly based on the geographic coordinate polygons rather than having to refer to a list of 'proper' SRSs.


title: ArcGIS 10 add-in: ComboBox add-in component focus issue

text: **Update 9/9/11:** Here is a minimal example project that will reproduce the issue for me (VS 2008 solution): http://wfurl.com/09061bf I have an ArcMap 10 add-in with a `ComboBox` class that I place on a toolbar. When I attempt to display a `MessageBox` after the user submits a value by typing in the ComboBox and pressing `Enter`, the ComboBox steals the keyboard focus from the MessageBox, so when the user tries to dismiss the MessageBox by pressing Enter again, the ComboBox's `OnEnter` event handler runs again instead of the MessageBox being closed. Without any checks against this, you might end up with numerous MessageBox instances on top of one another with no way to dismiss them except by clicking OK on each one. Attempts to work around this behavior have failed and include:   1. Disabling the ComboBox before showing the MessageBox and re-enabling it afterwards. The enabled state doesn't change until after the MessageBox is closed and the focus issue still occurs. Employing `Application.DoEvents` has no effect.   2. Setting a class-level boolean that tells the OnEnter handler to simply return when true (Set to True just before showing the MessageBox and False afterward.) The focus issue remains, meaning it takes 3 Enter keystrokes to dismiss the MessageBox, but at least the OnEnter code doesn't run again.   3. Setting the MessageBox's parent window to the ArcMap application using the `IWin32Window` overload (actually this is necessary to get the MessageBox to show in the middle of the ArcMap window instead of on top of the ComboBox).   4. Explicitly setting the MessageBox's MessageBoxOptions property to 0x1000 (`MB_SYSTEMMODAL`) or 0x2000 (`MB_TASKMODAL`)   5. Using `IApplication.StatusBar.set_Message()` instead of MessageBox. The ComboBox's HintText message, even if not specified, or set to an empty string, instantly overrides anything I write in the set_Message main pane (0), so this method is essentially useless with Addin ComboBoxes. Anyone have any other ideas as to what else I can try, or must I once again look at a WinAPI solution to deal with the ArcGIS framework's shortcomings?


title: Unable to get Mapbox leaflet-omnivore to load markers from .csv

text: I'm using the Mapbox.js API to build a really simple web map using a .csv file to place markers with tooltips. I'm trying to use the suggested leaflet- omnivore library to easily parse the csv and place the markers. Using this example guide I have been unable to get omnivore to place any markers at all. I've also noticed that just importing the library adds considerable load time to the map with the script:               <script src='//api.tiles.mapbox.com/mapbox.js/plugins/leaflet-omnivore/v0.0.1/leaflet-omnivore.min.js'></script>      Here is (some) of the .csv I'm trying to place:               Name,latitude,longitude,Sport,Contact,Website,Location     Oregon Avenue Octopi,39.915048,-75.170185,Softball,,www.octopisoftball.org,"Marconi Plaza Philadelphia PA 19148"     West Philly Waste,39.943065,-75.215689,Softball,,,"Kingsessing Rec Center Philadelphia PA 19143"     Pittsburgh Pounders,40.466298,-79.961581,Baseball,,,"Arsenal Park Pittsburgh PA 15201"      I've tried variations on the snippets:               omnivore.csv('/data/teams.csv').addTo(map);      and               omnivore.csv('/data/teams.csv', latfield:'lat', lonfield:'lng', delimiter:',').addTo(map);      My csv is in a proper location and should be able to be read. The obvious question is whether this is a problem with my code or with omnivore. I have had success inputting the markers individually but would much rather use a csv for easy adding and editing by multiple people. I'm also open to suggestions outside of omnivore if anyone has any. Thanks!


title: How do I copy polygon features with a WHERE clause using ArcPy.da

text: This question is similar to Copy a polygon geometry in Arcpy, except within the context of using Arcpy.da instead. I want to copy the geometry from one feature class to another using ArcPy.da (I'm not concerned about attributes of the feature class right now). I have a Python script that reads features from one feature class based upon a given WHERE clause using a search cursor and then inserts the features from the search cursor into another feature class:               import arceditor     import arcpy               inputfeatureclass = arcpy.GetParameterAsText(0)     insert_field = arcpy.GetParameterAsText(1)     where_field = arcpy.GetParameterAsText(2)     expression = arcpy.GetParameterAsText(3)     insertfeatureclass = arcpy.GetParameterAsText(4)               def get_specified_rows(inputfeatureclass, insert_field, where_field, expression):              searchcursor = arcpy.da.SearchCursor(inputfeatureclass, (insert_field, where_field), where_clause=expression)               return searchcursor          def insert_specified_rows(searchcursor, insertfeatureclass):              insertcursor = arcpy.da.InsertCursor(insertfeatureclass, ("SHAPE@"))         for row in searchcursor:             arcpy.AddMessage(row[0])             insertcursor.insertRow(row[0])          searchcursor = get_specified_rows(inputfeatureclass, insert_field, where_field, expression)     insert_specified_rows(searchcursor, insertfeatureclass)      The polygon **Shape** field is the value passed for my insert_field parameter. The program fails on insertcursor.insertRow(row[0]) with the following error: sequence size must match size of the row. I'm getting the following message on arcpy.AddMessage(row[0]): **[3088680.146493972, 10159411.210196389]**. This does not look like a polygon array, perhaps I am reading the polygon incorrectly using my search cursor?


title: Import Point Data to Oracle through FME with considering few other Point file attribute

text: I have three Point file (*.mif); Pole, Switch and Transformer Pole is the main feature which will be my Output, but I need additional information within pole layer attribute using information from Switch and Transformer layer.   * Pole layer have 3 attribute field: PoleID, Switch and Transformer   * The attribute I will use from Switch layer is SwitchType – (i.e. LV/HV etc.)   * From transformer, I want only the status of availability within pole (Yes / No )  ## Condition: a. I want to search any Switch which is available within 2mtrs buffer to a Pole b. If switch identified, then I will read the attribute field of the switch; SwitchType – (i.e. LV/HV etc.) and update in the particular Pole attribute in the Switch field. c. For Transformer, I want to search any Transformer which is available within 2mtrs buffer to a Pole d. If Transformer Identified, then the Pole attribute field (Transformer) will updated with Yes otherwise No e. Finally all the poles with updated attributes will import to Oracle Kindly help me to update these attribute updation using different point layer


title: Empty raster returned after GRASS i.atcorr with Landsat8 data

text: I've searched at length on previous forums but no solution seems to work for me. I'm trying to run the i.atcorr module on GRASS (version 6.4.3 on a Macbook Pro, Mavericks OSX) with Landsat8 imagery. Parameters didn't seem to be loading from the provided metadata file (MTL_txt), so I went back to the GRASS manuals and tried to input values by hand, here is what I used: 8 = geometrical conditions=Landsat ETM+ 05 04 03.81 98.117 23.131 = month day hh.ddd longitude latitude ("hh.ddd" is in GMT decimal hours) 1 = atmospheric mode=tropical 0 = aerosols model=none (ASSUMING NO AEROSOLS) 0 = visibility [km] 0.732 = aerosol optical depth 550 nm -1.408 = mean target elevation above sea level [km] -1000 = sensor on board a satellite 0 = THIS IS WHAT I'M WORRIED ABOUT: There is no pre-defined coding for Landsat 8 satellites and I am unsure how to define my own spectral conditions?? I checked I was in the correct region (i.e. used g.region) and I've tried using both radiance and reflectance at sensor values (computed using i.landsat.toar, which seems to have worked nicely) but the rasters returned are empty :( **Here is what I'm running (in verbose)** i.atcorr -a --overwrite --verbose iimg=Part1toar_B1@Charlie icnd=/Users/charlottewhitham/Documents/PhD/My research/my GIS/grassdata/NGH/Charlie/.tmp/Charlie.local/819.2 oimg=P1toar_atcorrB1 * **_*_ ** _*_ ** _*_ ** _*_ ** _*_ ***** 6s version 4.2b * * *   * geometrical conditions identity *   * \------------------------------- *   * etm+ observation * * *   * month: 5 day: 4 *   * solar zenith angle: 23.55 deg solar azimuthal angle: 103.24 deg *   * view zenith angle: 0.00 deg view azimuthal angle: 0.00 deg *   * scattering angle: 156.45 deg azimuthal angle difference: 103.24 deg * * *   * atmospheric model description *   * \----------------------------- *   * atmospheric model identity : *   * tropical (uh2o=4.12g/cm2,uo3 =.247cm-atm) * * *   * aerosols type identity : *   * no aerosols computed * * *   * spectral condition *   * \------------------ *   * constant *   * value of filter function : *   * wl inf=27640754752056197120.000 mic wl sup= 0.000 mic * * *   * target type *   * \----------- *   * homogeneous ground *   * constant reflectance over the spectra 0.000 * * *   * target elevation description *   * \---------------------------- *   * ground pressure [mb] 862.22 *   * ground altitude [km] 1.408 *   * gaseous content at target level: *   * uh2o= 2.056 g/cm2 uo3= 0.244 cm-atm * * *   * atmospheric correction activated *   * \-------------------------------- * Atmospheric correction... Atmospheric correction complete. (Fri Dec 6 12:16:10 2013) Command finished (22 sec)  **And here is a description of the output:** | Layer: P1toar_atcorrB1@Charlie Date: Fri Dec 6 12:12:58 2013 | | Mapset: Charlie Login of Creator: charlottewhitha | | Location: NGH | | DataBase: /Users/charlottewhitham/Documents/PhD/My research/my GIS/grassda | | Title: ( P1toar_atcorrB1 ) | | Timestamp: none | |----------------------------------------------------------------------------| | | | Type of Map: raster Number of Categories: 0 | | Data Type: CELL | | Rows: 7381 | | Columns: 7551 | | Total Cells: 55733931 | | Projection: UTM (zone 47) | | N: 2666715 S: 2445285 Res: 30 | | E: 522315 W: 295785 Res: 30 | | Range of data: min = NULL max = NULL | | | | Data Description: | | generated by i.atcorr | | | | Comments: | | i.atcorr -a iimg="Part1toar_B1@Charlie" iscl=0,255 icnd="/Users/char\ | | lottewhitham/Documents/PhD/My research/my GIS/grassdata/NGH/Charlie/\ | | icnd_lsat1.rtf" oimg="P1toar_atcorrB1" oscl=0,255 As you can see the range of data is NULL to NULL... not good. I would be very grateful if somebody could please help me with this - I'm really unsure what to try next! Thank you very much in advance and please let me know if I need to provide more info, Cheers, Charlie


title: Creating outline shp - georeferencing existing map/pictures

text: So I need to create shapefiles (that need to be adapted and transfered for other programs) with the district outline of a couple of Caribbean countries. Is georeferencing (using the _Georeferencing toolbar_ ) of pictures/maps online a good solution? Is there any other method I can approach, possibly more professional? What other measure should I take in order to make the shapefiles created can be used after (without basemap)? (Only thing off is Coordinate system, which in my case is WGS 1984) Realise this is a newbie question again, sorry, but help will be very much appreciated. Thanks


title: How to handle the New ESRI Ocean Basemap

text: Some of you may have seen this recent (4/28/14) ESRI blogpost - The New Ocean Basemap. ESRI, in their infinite wisdom, has split up their Ocean Basemap into 2 separate services. One with just the map tiles. And one with just the labels. And they will 'retire' the old basemap in 2015. My first thought: Cool, now I'll be able to drop the labels from their basemap if I want. My second thought: But what about those times when I don't want to drop the labels...? With basemaps being mutually exclusive, am I really going to need to build special circumstances into my mapping applications (to say whenever the ESRI Ocean Basemap is active, make sure this other ESRI Ocean Label overlay/layer is also active...etc). What a PITA. Are they really making it this difficult to display their Ocean Basemap **with** the Labels??? Have I mis-read their blogpost?   Am I missing the better/obvious solution? BTW I've created an openlayers example in a jsfiddle so you can see the 3 ESRI offerings in place (Original Ocean Basemap, New Ocean Basemap, and the New Ocean Labels).


title: Add x y data from existing attribute table

text: In Quantum GIS, I already imported an excel or dbf table. Now I need to tell Quantum GIS which columns are Latitude and Longitude. How do I do that? I know there is Add Delimited Text Layer which allows that, but it seems Unreasonable to me to convert XLS or DBF file to text file before importing! There must be some direct way of telling QGIS that _these_ parameters of the imported attribute table are actualy x y data.


title: How to restrict level crossings in route?

text: I want to evaluate if it is useful to bypass a level crossing through a nearby tunnel or overpass when the LC is closed or if it is worth waiting. I want to use the OSM data for realistic scenarios and calculate the most time efficient route for closed and opened LC. Is it possible to solve this problem with pgRouting? Or can you give me some hints where to start?


title: map stops at map.graphics.add(graphic);

text: function getLocRes(results) {     //map.graphics.clear();     var symbol = new esri.symbol.SimpleFillSymbol(esri.symbol.SimpleFillSymbol.STYLE_SOLID, new esri.symbol.SimpleLineSymbol(esri.symbol.SimpleLineSymbol.STYLE_SOLID, new dojo.Color([80, 0, 0]), 2), new dojo.Color([80, 0, 0, 0.5]));     var combinedExtent = new esri.geometry.Extent();     console.log("er",results.features.length);     for (var i = 0;i < results.features.length; i++) {console.log("i=",i);         var polyGraphic = new esri.Graphic(results.features[i], symbol);         console.log("polyGraphic",polyGraphic);         polyGraphic.setSymbol(symbol);         //map.graphics.add(polyGraphic);     }     combinedExtent = esri.graphicsExtent(results.features);     map.setExtent(combinedExtent, true);      } I made a function that shows multiple results in polygon. when I test it with out any map.whatever sort of working however, neither of map.graphics.clear() or map.grahpics.add() make the map work am I missing something here?


title: How to publish from ArcMap to GeoServer (not using GeoCat) - GeoCat Alternative?

text: How can I publish a map from ArcMap 10 to OpenGeo's GeoServer without using GeoCat Bridge? Are there any alternative solutions to accomplishing this? Thanks!


title: what criteria should I consider in buying a pocket GPS

text: It's my first time making this kind of choices, and I really don't know how to go with a choice or another. I basically know that it should have good **_durability_** , a good **_battery life_** and good **_coverage_**. What are other **_practical_** considerations I should be aware of ? And what about **_accuracy_** ? I actually need to have a pocked GPS to pick points from the electrical network (the nodes are the electrical posts), all to create a cartography of 1/25000 scale. I should cover 500 points or so. Thanks. P.S : this question is better being a wiki, I think.


title: OpenGeo Suite 3.1 doesn't work anymore when setting up Client SDK 3.0.2

text: We've installed OpenGeo Suite 3.1 on Ubuntu 12.04 LTS and everything works fine. When trying to set up the OpenGeo Client SDK 3.0.2 (following http://suite.opengeo.org/opengeo-docs/apps/clientsdk.html), however, the OpenGeo Suite doesn't work anymore. After inserting the SDK in PATH and adding ANT_HOME and JAVA_HOME in /etc/environment, tomcat6 starts up, but doesn't serve anything in the browser anymore. The severe error in the log looks as follows:               24 Jul 15:35:24 WARN [config.CustomEditorConfigurer] - Passing PropertyEditor instances into CustomEditorConfigurer is deprecated: use PropertyEditorRegistrars or PropertyEditor class names instead. Offending key [org.geotools.util.Version; offending editor instance: org.geoserver.platform.util.VersionPropertyEditor@3f705316          24.07.2013 15:35:24 org.geoserver.platform.GeoServerExtensions checkContext     SCHWERWIEGEND: Extension lookup occured, but ApplicationContext is unset.          24 Jul 15:35:24 ERROR [geoserver.global] -     - GEOSERVER_DATA_DIR: /usr/share/opengeo-suite-data/geoserver_data      Any help would be greatly appreciated. Are there conflicts between the versions? A new SDK doesn't seem to exist yet. EDIT: The error from the log doesn't seem to have to do anything with our problem, since it is also logged in another working geoserver environment.


title: display a choropleth map using D3 and Leaftlet

text: I am trying to get a leaflet map to display a choropleth map using d3. I have written this code using samples I have found however the map just displays as the black fill, not the color classes I defined. This is the code I am using. Anyone have an idea about what I'm doing wrong?               <!DOCTYPE html>               <meta charset="utf-8" />          <meta name="viewport" content="width=device-width, initial-scale=1.0">          <link rel="stylesheet" href="http://cdn.leafletjs.com/leaflet-0.5/leaflet.css" />                    <div id="map" style="width: 900px; height: 500px"


title: Find locations of raster maximum within polygons: Basin Delineation

text: I'm trying to identify the Location (lat/lon) and upstream area of every raster point that flows into the Arctic Ocean (there are a lot so it has to be fully automated). I'm starting with the following rasters:   * Flow Direction   * Accumlated Flow Area My approach so far has been to use the above rasters to create a delineated basin shapefile. From there I was hoping I could use the shapefile to look at the Flow Accumulation Raster and identify the location of the maximum value in each basin. In terms of output, I just need a table that shows [Basin#, Longitude, Latitude, UpstreamArea] for every basin. Any ideas or tips would be apreciated. [Note that I've already looked here but I don't know how to restrict the raster to only the area bounded by the polygons in the shapefile] * _Using: ArcGIS10, Python, etc._ Thanks, Joe


title: Is there are QGIS plugin to allow the 3d visualisation of geological borehole data similar to the functionality of Target for ArcGIS?

text: I'm very new to QGIS, I have started using it as my new company is not willing to pay for ArcGIS and Target at this stage so I'm hoping QGIS can substitute not just in the interim but for good. What I am looking for is a way to load in my borehole database and create 3d visualizations of the data, make cross sections and if possible model seams. Is there any capability within the QGIS world to undertake this kind of work?


title: Displaying information on a map using an InfoWindow

text: I am wanting to display information in an InfoWindow when the user clicks on a point on the map. I have implemented the InfoWindow in XAML, but I do not understand the code behind it. Any information will be appreciated. Here is my (what I feel is relevant) XAML:                <esri:InfoWindow x:Name="MyInfoWindow"                                  Padding="2"                                  CornerRadius="20"                                  Background="{StaticResource PanelGradient}"                                  Map="{Binding MyMap}"                                  ContentTemplate="{StaticResource MyFeatureLayerInfoWindowTemplate}"                              MouseLeftButtonUp="MyInfoWindow_MouseLeftButtonUp" />      If you need anymore information, let me know. I have the information being displayed correctly at the moment in a MapTip. Here is the XAML for that:               <esri:FeatureLayer.MapTip>                             <Border CornerRadius="10" BorderBrush="Black" BorderThickness="2" Margin="0,0,15,15" Background="{StaticResource PanelGradient}">                                 <StackPanel Margin="7">                                     <TextBlock Text="{Binding [WVAL_ID], StringFormat='ID: \{0\}'}" FontWeight="Bold" Foreground="White"  />                                     <TextBlock Text="{Binding [SIZE], StringFormat='Size: \{0\}'}" FontWeight="Bold" Foreground="White" />                                     <TextBlock Text="{Binding [WVAL_USE], StringFormat='Use: \{0\}'}" FontWeight="Bold" Foreground="White" />                                     <TextBlock Text="{Binding [WVAL_TYPE], StringFormat='Type: \{0\}'}" FontWeight="Bold" Foreground="White"  />                                     <TextBlock Text="{Binding [OPENS], StringFormat='Opens: \{0\}'}" FontWeight="Bold" Foreground="White"  />                                     <TextBlock Text="{Binding [TURNS], StringFormat='Turns: \{0\}'}" FontWeight="Bold" Foreground="White"  />                                     <StackPanel Orientation="Horizontal" >                                          </StackPanel>                                 </StackPanel>                             </Border>                         </esri:FeatureLayer.MapTip>


title: Using ArcGIS ModelBuilder with a "While" Iterator

text: I am struggling to get the ArcGIS ModelBuilder "While" Iterator to work as I would expect it. I have a hunch it's user error, but even reading the help hasn't seemed to get it working for me. I have a simple model that calls a Python script that simple does a `time.sleep(15)` and then sends back the conditional parameter of "TRUE" or "FALSE" using `arcpy.SetParameterAsText(1, "FALSE")` or similar for True statements in an if/else statement. I've also tried every combination of TRUE, True, and true, as outputs from my Python script without success. ![Model with WHILE Iterator](http://i.stack.imgur.com/G4R7u.jpg) I can see in the results that my conditions are being set correctly:               Executing (While): IterateCondition FALSE TRUE     Start Time: Tue Nov 20 16:10:06 2012     Succeeded at Tue Nov 20 16:10:06 2012 (Elapsed Time: 0.00 seconds)      and               Executing (While): IterateCondition TRUE TRUE     Start Time: Tue Nov 20 16:26:41 2012     Succeeded at Tue Nov 20 16:26:41 2012 (Elapsed Time: 0.00 seconds)      but the WHILE Iterator continues to fire the sleep script, even if I change the output parameter of my script. Any ideas on why I cannot seem to get the sleep script to stop firing? What should I be passing from my SetParametersAsText() in order to get the While Iterator to stop firing the Sleep script?


title: Is it possible to export cross sections from QGIS to HEC-RAS

text: Is it possible to create cross sections along a river in QGIS and then export the data to HEC-RAS? Thanks


title: Create Array From Selected Features Based on the Order in Which They Were Selected

text: **I need to create an array of features from a selection based on the order in which they were selected.** As in: I select feature 5,2,3,4. I need to have an array so I can apply attributes (with python) to to these features based on order in which they were selected. This is too apply part numbers to pre-mapped survey plans. To add a part to each is time consuming when their are 50 parts. This would allow me to select them in the order of part number, then apply the correct part with a button I have already created with python. I've been trying for a long time to figure this out and I have run out of solutions. I am pretty well versed in Python and with that ArcPy. Seem's like this should be easier then it is. Thanks for any help.


title: Link attribute value to Fanout in FME?

text: I would like to use attribute value from Workbench as an attribute that says Yes or No to the Fanout options in Workbench settings. So apply value from WB to the advanced fanout settings of WB. Any idea? Someone, please?


title: If statement in ModelBuilder to check if field exists

text: I'm trying to create the following in ModelBuilder:               Set input feature class       if field1 exists         execute the processes in branch1       else         add field1         execute the processes in branch2         merge branch      While I have no problem with the processes in the two branches, the `IF` check at the start doesn't seem to be working the way I want it to work. I am using this code in the `Calculate Value` code block               def ifFieldExists(Input_FC, SUMMARY):       import arcpy       fields = arcpy.ListFields(Input_FC)       for f in fields:         if f ==" SUMMARY":           return "true"         else:           return "false"      and expression               ifFieldExists(r"%Input_FC%", "%SUMMARY%")      This is the precondition for `branch1`, so if it returns `true`, the field exists so `branch1` executes. The other branch's precondition returns the reverse. When I run it on an input that does contain the field, the precondition returns `false` instead of `true`, so the wrong branch executes.


title: How to convert fake donut holes in polygons to real ones?

text: I have this problem where I have to import shapefiles into our system which must be geometrically valid, but the software that everyone uses to produce them seems to output files where every polygon that should contain a hole is instead made up from a single outer ring which goes around the boundary and then goes inward along some arbitrary path, traces around what should be the inner ring, and then doubles back on itself back to the outside. What I seek is a way to resolve these by somehow dissolving the vertices along the arbitrary path in and back out of the hole, to make a seperate inner and outer ring which will form a valid polygon. I have tried the trial version of ArcGIS and it seems to do this automagically when you snap the vertices on the double-back part together, but my business doesn't have a copy of ArcGIS and I might have a hard time convincing my boss to buy it, so I was hoping QGIS would allow me to do such a thing? Screenshot: ![description](http://i.imgur.com/zGwtG.png) I couldn't figure out how to get any WKT out of this sorry, but from playing around with the vertices I'm quite sure that it's one single ring that doubles back on itself for a bit.


title: How to compute a map-to-screen transformation?

text: Every projection calculation I have come across dealing with the Lambert Conic projection converts the given Lat & Long to meters. So now the magical question: **how on earth do you take the meters and convert to X, Y in pixels?** I have a map of the US. The following is the projection info:               ullat=53.3638     ullon=-130.964     lrlat=18.2977     lrlon=-69.1058          sp1=0     sp2=39     rlon=-95.5     rlat=36.5     width=640     height=480      I converted this http://svn.onelab.eu/deprecated.PLCWWW/imports/import-001-head-04-24/plot- latlong/plot-latlong over to php, but get returned meters. I am working on an open source project. Most of the users will not be able to install Proj4, so that is not an option (even there the return is in meters). Anybody?


title: Understanding null shape type

text: I have some shapefiles that I am reading and then loading into tables in the database. The shapefiles with no zoom (eg. "au_postcode") have no shapes with type null. When I tried to load a shapefile with zoom (eg. "au_postcode_100"), I saw some null shapes. It would be great if someone can please explain to me how the null shapes work and why would a shape file be null. Its kinda confusing because to me a shape is at least a point and if there is a physical location, there would be a shape for it. I am new to GIS and if something in this question does not make sense or I've used the wrong terminology, please let me know. Thanks


title: How to create a polygon from other polygon boundaries?

text: I want to write an algorithm that creates a polygon inside a hole made by a set of other polygons. It's something like a flood fill and it's like the "boundary" tool in Autocad. The image shows an example: the algoritm should fill the enclosed white space in the center (after a click on it). Any theory, pseudocode or examples would help a lot. ![enter image description here](http://i.stack.imgur.com/IIc4P.jpg)


title: Looking for a sample dataset of international addresses

text: I need a test dataset of 100 addresses in each country's local langauge. This is approximately 10,000 addresses. The don't need to be geocoded or resolved to address components. For example, 1 microsoft way seattle wash xxxxx, or 5814 pontiac dr san jose, ca 95123 Почтовый адрес: 125009, Москва, Театральная пл., 1. are comlete addresses with enough data. Similar ones for addresses in each world langauge would be fine. Can I get such a non geocoded, random data set for low or no cost? This is to test a google geocoded address cache (as per the API recommendations)


title: Link non-spatial SQL Server query to a Geodatabase?

text: Is it possible to build an object (of whatever kind) in ArcEngine code that links the geodatabase data with a non-spatial SQL Server query? How would you go about doing this?


title: DIY easy and ~accurate small-scale landscape survey

text: ## Summary Is there any equipment that I can rent 'inexpensively' (less than $400) that would allow 1-2 people to accurately (+/- 6") capture the topography of hilly terrain of perhaps 250' on a side in an afternoon? ## Details My wife and I are about to embark upon a redesign of the landscape around our home in the mountains. Our lot started out on a generally 1 in 3 slope, and this has (of course) been changed in many interesting ways during building. There are small hills, patios, semi-flattened areas followed by large falloffs. As part of our designing, we want to capture exactly what the landscape looks like now, digitally. (Then we'll sculpt in the computer what we want, taking into account existing boundaries and surrounding contours, and see roughly how much dirt we need or will need to remove.) We performed a ghetto survey this weekend with strings and poles, roughly sampling a grid of elevations every 10'. The results we got produced a reasonable-looking representation of the terrain. However, the 10' sampling was a bit too coarse for my tastes and—more importantly—we were only able to sample about 1/4 of the area we need before our measuring poles weren't tall enough or we ran out of reference objects to attach to. So: what hardware exists that I can rent to do this task for me? Ideally I'd like to take a pole, walk all over the property and every few feet stick it on the ground, press a button, hear it go "beep!" and have it record an XYZ sample that is accurate within +/- 6". After doing this for a couple hours I'd plug it into the computer and download some DEM, CSV, or other accessible data of all the points in a non-proprietary format. It looks like Trimble offers some neat gadgets in this arena, and I see that some companies rent this sort of gear out, but not being a landscape survey engineer professional I cannot figure out exactly which, if any, products meet my needs. Can someone who knows this product space please help suggest what I should be looking for? Thanks! _Note: I originally asked this question on Home Improvement, butas suggested there and approved here on meta I'm cross-posting it to this site._


title: Using Python to combine middle parts of two fields into third field

text: Below is how I used the Update Cursor in my final script. As you can see getValue was used a few times, perhaps there is a cleaner way to do this? But my concatenated field wouldn't populate properly otherwise.                   sourceRows = arcpy.UpdateCursor(inFeatures)         source = sourceRows.next()               while source:             fieldString1 = source.getValue(fieldName1)             fieldString2 = source.getValue(fieldName2)             concatString= (str(source.getValue(fieldName2))[1:7] + str(source.getValue(fieldName1))[1:7])             source.setValue(fieldName3, concatString)             sourceRows.updateRow(source)             source = sourceRows.next()         del source, sourceRows      I have a script that creates three fields in a feature class. The script adds the X and Y coordinates to the first two fields. I would like data found in the middle parts of the first two fields to populate the third field. Does anyone know how to find the middle part of a field using Python? I see where VBScript seems to have a MID function but I haven't found anything for Python yet. Below is what the third field (fieldName3) should look like.               fieldName1     fieldName2   fieldName3     12345.321       36789.321   67892345     24680.5425      98765.432   87654680      Do you see how I want data to be combined in fieldName3? Below is the script as I have it now.               import arcpy     from arcpy import env     arcpy.env.overwriteOutput = True          env.workspace = "D:/gisfiles/test.gdb/"          inFeatures = "fc"     fieldName1 = "xCoord"     fieldName2 = "yCoord"     fieldName3 = "Combined"     fieldPrecision = 18     fieldScale = 11          expression1 = "float(!SHAPE.CENTROID!.split()[0])"     expression2 = "float(!SHAPE.CENTROID!.split()[1])"          arcpy.AddField_management(inFeatures, fieldName1, "DOUBLE", fieldPrecision, fieldScale)     arcpy.AddField_management(inFeatures, fieldName2, "DOUBLE", fieldPrecision, fieldScale)     arcpy.AddField_management(inFeatures, fieldName3, "DOUBLE", fieldPrecision, fieldScale)          arcpy.CalculateField_management(inFeatures, fieldName1, expression1,                                 "PYTHON")     arcpy.CalculateField_management(inFeatures, fieldName2, expression2,                                 "PYTHON")          sourceRows = arcpy.UpdateCursor(inFeatures)     source = sourceRows.next()          while source:             **#Instead of 12345 I need an expression to combine fieldName1 and fieldName2**             source.setValue(fieldName3, "12345")             sourceRows.updateRow(source)             source = sourceRows.next()     del source, sourceRows


title: How to make beautiful maps in R?

text: There are quite some nice possibilities to analyse spatial data in R and in the context of my current project I would like to use R more frequently to do that. Until now I plot my maps with ggplot2 package that brings a lot of practical tools to plot and explore data. Still my maps are not nearly as good looking as the ones I used to plot in ArcGIS. So for publishing purposes I wonder if there are any good tutorials, books, practical tips, packages etc. that I could use to experiment a little bit and make my maps more good looking. Any advice would be warmly welcome :)


title: How to avoid corrupt geometries when translating between TopoJSON and GeoJSON?

text: I'm trying to display some geojson vector features on a leaflet map. The size of my GeoJSON is too big (25.4mb of data, 8.47mb gzipped), so I'm trying to use a TopoJSON format to reduce the size without loosing detail. I converted my geojson data with topojson tool and just basic params: topojson --id-property okato_code -p name -o adm.topojson adm.json That significantly reduced the file size (down to 1.54mb of data, 201.3kb gzipped). But when I translate it back ( using topojson.object() ) to GeoJSON on the client, it looks overgeneralized and like it has some incorrect topologies when zooming in. I thought there was no generalization happening during translation to/from topoJSON, then why are the geometries different? How to make them identical? My goal is just to reduce the size of the file without any simplification of geometries. Here's the map - http://basemap.ru/json/ You have to wait for the big GeoJSON to download, it's included for comparison, but after it loads, you can see a distinct difference on high zoom levels. GeoJSON layer is orange and TopoJSON layer is white. I've included the layer switcher to turn topojson and geojson on and off. The data itself is here: http://basemap.ru/json/json/adm.geojson http://basemap.ru/json/json/adm.topojson The script with translation is http://basemap.ru/json/mapscript.js If it makes sense, my map is in Asia Lambert Conformal Conic projection (http://spatialreference.org/ref/esri/102012/)


title: How do I fix this QGIS sextante error? NameError: name 'sextante' is not defined

text: A colleague of mine posted this the other day. We haven't got any nibbles yet and I wanted to elaborate on the issue. We are having trouble running the r.mapcalculator function and receive the error below ![enter image description here](http://i.stack.imgur.com/GQ0v9.png) I have gone into the file "HistoryDialog.py" and changed a reference to Sextante on line 101 to 'sextante' (original text shown below). I think I am clutching at straws here, because after this change I receive another error altogether. ![enter image description here](http://i.stack.imgur.com/LaPnY.png) A good googling returned nothing related to the 'sextante is not defined error'. Surely someone else has come across this before. Also for us it is happening on Ubuntu and Windows. We installed the latest dev version through Osgeo4w installer but no luck!!


title: How to calculate altitude with given lat,long position?

text: I'm new to GIS. I want to calculate the elevation of any point on earth giving lat, Long (for example like google elevation Api, or the earth tools: http://www.earthtools.org/webservices.htm # cheigit) based the ASTER GDEM MODEL V2 I went to download data from there, and gives me geotif files. Then what should I do? Thanks


title: What is difference between Geo BI (Geo-business Intelligence) and GIS (Geographic Information System)?

text: I tried finding information with difference between Geo Business Intelligence and Geo IS but unfortunately no clearly answer being found. What is different between Geo BI and GIS and its advantage/disadvantage?


title: LinearRing getVertices implementation

text: I'm fairly new to OpenLayers so please excuse me if this is an ignorant question. I'm curious why the LinearRing getVertices implementation               getVertices: function(nodes) {         return (nodes === true) ? [] : this.components.slice(0, this.components.length-1);      },      does not return the last point since I was under the impression a LinearRing was supposed to represent a closed ring. I was using getVertices to serialize to another format which requires closure and found this interesting. Can someone enlighten me?


title: How to create shaded relief in ArcGIS?

text: How can I create shaded relief from contour lines (.shp format)? I see Hillshade in ArcToolbox, but it asks for a raster file. I have ArcGIS 10.0.


title: GIS research areas in Telecommunications?

text: I am doing MS in GIS and Remote Sensing; i have a telecommunications Engineering Bachelors degree. Can any one please guide me what sort of GIS research can i do for telecommuncations industry???


title: Modelling two-way single-track/one-lane width road in ArcGIS Network Analyst?

text: I don't know how make one-lane roads(only room for one car at a time) for one car waiting for another car to pass by in ArcGIS. Yesterday, I asked about the wait time when the two cars meet in the road, one car should wait for the other car to pass that as How to calculate wait time on a one-lane road with two-way travel? But I realized that I can't even make a one-lane road... How can I make this road in ArcGIS Network analyst (especially, VRP (Vehicle Routing Problem))?


title: Why does QGIS not detect the CSV/DBF column headers?

text: I am trying to import a csv or dbf into qGIS so I can perform a join on a shapefile. However, qGIS is not detecting the first row as header names. I know I've done similar joins before and there's never been a problem detecting column headers. Any help would be much appreciated. Thanks.


title: Displaying 3D subsoil models in ArcScene

text: I was wondering if anybody had any ideas on how to display a 3D subsoil model in arcscene? The data we will have will be exported from rockworks in point form - and then we need to convert it to - TINs (? i assume). I have never done any extensive work in arcscene before and i couldn't find anything too helpful on the internet!


title: How to determine which EditEvent is fired BEFORE the feature is created?

text: I have written an AddIn-Tool and want to handle the OnSketchFinished-Event. The Way               ((IEditEvents_event)_myEditSketch).OnSketchFinished += OnFinishedSketch;      doesn't really work for me, cause I need the get the event before the feature is created. I rather need something like the OnFinishSketch that was fired by the IEditTask. Does anyone have an idea how to solve this problem? I'm very thankful for any suggestions.


title: Can a raster map in GRASS be saved in JPEG/BMP format?

text: I'm trying to save a raster map in GRASS in jpeg/bmp format without using the GUI for d.out.file/d.save. Is there a script to do it? I tried r.out.gdal but it produced a map that is black and white.


title: How can I import ALOHA .pas plume file into QGIS?

text: I am trying to find a way to import ALOHA's plume model (.PAS) into QGIS? I have not been able to find anything showing how to accomplish this outside of using ESRI-based tools.


title: Why is st_intersects giving bad results?

text: I am using PostgreSQL version 9.2.2 with PostGIS 2.0.2 (or just 2.0, not sure), and need to do simple queries to find points that intersect with a polygon. The following test case returns only partial results for an unknown reason. If the minimum longitude value of the polygon is changed from 0 to 1 to 2 to 5 to 20 (leaving the latitudes limits as they are), different numbers of the points are returned, even though all of the points fall into the polygon for all these values.               CREATE TABLE testing.tsttab     (       lat real,       lon real,       tstpoint geography(Point,4326)     );          INSERT INTO testing.tsttab (lat, lon, tstpoint) values (36.1, 22.1, ST_GeographyFromText('POINT(22.1 36.1)'));     INSERT INTO testing.tsttab (lat, lon, tstpoint) values (36.2, 22.2, ST_GeographyFromText('POINT(22.2 36.2)'));     INSERT INTO testing.tsttab (lat, lon, tstpoint) values (36.3, 22.3, ST_GeographyFromText('POINT(22.3 36.3)'));     INSERT INTO testing.tsttab (lat, lon, tstpoint) values (36.4, 22.4, ST_GeographyFromText('POINT(22.4 36.4)'));     INSERT INTO testing.tsttab (lat, lon, tstpoint) values (36.5, 22.5, ST_GeographyFromText('POINT(22.5 36.5)'));     INSERT INTO testing.tsttab (lat, lon, tstpoint) values (36.6, 23.1, ST_GeographyFromText('POINT(23.1 36.6)'));     INSERT INTO testing.tsttab (lat, lon, tstpoint) values (36.7, 23.2, ST_GeographyFromText('POINT(23.2 36.7)'));          select * from testing.tsttab     where st_intersects(     st_geogfromtext('POLYGON((2 35, 2 37, 40 37, 40 35, 2 35))'),tstpoint);      If I convert to GEOMETRY rather than GEOGRAPHY, it performs correctly. Am I mis-using something?


title: ArcObjects - Create and Draw a Transparent Polygon

text: I'm trying to create a transparent Polygon in ArcMap; however, it's never transparent. Here's the code I have.               // Create polygon     IPolygon polygon = new PolygonClass();     // create & add points,......          // Create lineSymbol     ISimpleLineSymbol lineSymbol = new SimpleLineSymbolClass();     // set color & style,....          // Create fillSymbol     ISimpleFillSymbol fillSymbol = new SimpleFillSymbolClass();     fillSymbol.Outline = lineSymbol;     fillSymbol.Color = //set to color red//;     fillSymbol.Color.Transparent = 0;     fillSymbol.Style = esriSimpleFillStyle.esriSFSSolid;          // Draw Shape     myMap.DrawShape(polygon, fillSymbol);      RESULT ==> A solid red polygon is drawn. The Transparent setting is ignored. Thoughts?


title: NASA World Wind Pick From A Group Of Icons

text: I am using an IconLayer in NASA World Wind (Java SDK) to show positioned items on a map which works fine except for when several are in a group and start to overlay one another. The display of the icons in a group like that is not so much of a problem as I need to allow the end user to easily select the correct icon in that "pile" easily. A workaround would be to allow for a right-click on an icon in that group would generate a context menu that shows a list of the items in that "pile" so they could select the one they want. Creating a context menu on right-click is easy enough but I do not see how I could create a buffer from that point to do a simple geometry relationship like "contains" in NASA World Wind. I have searched the API a bit and have not seen anything that looks like it has potential so I was hoping that someone here might now of either how to do this in NASA World Wind or a simple way of including or achieving this functionality.


title: Double tap on the map not work in the chrome of Nexus 7

text:
I am stuck with a problem which is the double tap not working in the chrome of Nexus7. Say you have code like below. please help to review it .               <!DOCTYPE html>                   <title></title>         <meta http-equiv="X-UA-Compatible" content="IE=7, IE=9,IE=10" />         <meta name="viewport" content="width=device-width, initial-scale=1.0, maximum-scale=1.0, user-scalable=0">         <link href="http://js.arcgis.com/3.9/js/dojo/dojo/resources/dojo.css" rel="stylesheet">         <link href="http://js.arcgis.com/3.9/js/dojo/dijit/themes/claro/claro.css" rel="stylesheet">         <link href="http://js.arcgis.com/3.9/js/esri/css/esri.css" rel="stylesheet">

                     Could someone please help me out ? **Updated** I test above script in below devices. Please review it.   1. Surfari in IPad2 (Double tap works fine.).   2. IE11 in Surface (Double tap not works) (and Why the feature esri-touch test is false? and the feature esri-pointer test is true.).   3. Chrome in Window7 computer(Mouse double click works).   4. Chrome in Nexus7 (Double tap not works).


Where can I download the last default.style for osm2pgsql? I get `http://svn.openstreetmap.org/applications/utils/export/osm2pgsql/default.style - this directory does not exist` when I try to download it.


title: problems of display modis ndvi with TIMESAT

text: I use MODIS MRT tool mosicia tow tiles MODIS NDVI data, then use TIMESAT image view to display the mosaic data, encounting file size waring : ![enter image description here](http://i.stack.imgur.com/Fp4W3.png) and after drawing data, the sequences of two tiles is wrong, original data display sequence is right to left, but display result is left to right. ![enter image description here](http://i.stack.imgur.com/jlMsH.png) but when i use ARCGIS or envi to display the mosaic data, the display sequence is ok. ![enter image description here](http://i.stack.imgur.com/s8FiS.png) so i want to know what is wrong with my mosaic data?


title: Leaflet: Is there a way to use click event just like mouseout?

text: I made a map with leaflet by following the official Leaflet tutorial on interactive choropleth maps. I would like to make the map compatible with smartphones, so I added a click event, which acts like hovering on a feature.               function onEachFeature(feature, layer) {     layer.on({         mouseover: highlightFeature,         mouseout: resetHighlight,         click: highlightFeature     });     }      The consequence: When a user is clicking, the `hightlightFeatures` remain highlighted. Is there a way to prevent this? In other words: Is there an event on the lines of `mouseout`, in order to activate `resetHighlight` by clicking another feature?


title: How to find all intersections and count lines arriving at each intersection?

text: I'm having a hard time getting off the ground with this, and am trying to move away from ArcGIS / arcpy. I have a shapefile with a network of polylines. I want to find/count all 2-way intersections, 3-way intersections, 4-way, etc. This would involve, I think:   * splitting all lines at intersections   * deleting identical segments (with some spatial tolerance)   * create points at the end of each new line segment   * buffer the line ends, and count the number of points near/at each line end I would like to do this in a python script. Is this overall plan easily implementable in pyQGIS? An outline of the calls I would be using would be very helpful.


title: How to count the number of fields in a table using arcpy?

text: arcpy contains a ListFields method, which "Lists the fields in a feature class, shapefile, or table in a specified dataset". The resulting object has a `count` method, but I can't see how to return the field count from this (shouldn't this return an integer?). Aside from iterating through the fields and incrementing a counter, how can I return the field count?               import arcpy     fields = arcpy.ListFields(<table>)     count = fields.count          [Dbg]>>> count     <built-in method count of list object at 0x0EDEABE8>


title: Images exported from ArcGIS Online as File Geodatabase inaccessible to ArcGIS for Desktop?

text: I am a new user of ArcGIS Online. I have a system that includes features and images of the objects and I use Esri ArcGIS for Desktop to perform data collecting. Before exporting the whole database as FGDB I use attachment management to download the images as FGDB zip file. How can I reach the images in the exported FGDB file using ArcGIS for Desktop?


title: node.js GeoLib fails on polygon that straddles the anti-meridian

text: geolib.isPointInside does not appear to work correctly when the polygon longitudes straddle the anti-meridian. I maintain a system with lots of polygons that do straddle the antimeridian (I'm in Australia :-) so it's important to me. see the test code below that illustrates the problem when it's run in node:               var geolib = require('geolib');          // -------------     // Test 1 This succeeds     // -------------     var polygon = [         {latitude: -46, longitude: 1},         {latitude: -46, longitude: 2},         {latitude: -47, longitude: 2},         {latitude: -47, longitude: 1}     ];     console.log('Test 1: Polygon:',polygon);     console.log('should be true', geolib.isPointInside({ latitude:-46.5, longitude:1.5}, polygon))     console.log('should be false', geolib.isPointInside({ latitude:-47.9, longitude:1.2}, polygon))     console.log('should be true', geolib.isPointInside({ latitude:-46.8, longitude:1.1}, polygon))     console.log('------------');          // -------------     // Test 2 This succeeds     // -------------     var polygon = [         {latitude: 10, longitude: -10},         {latitude: 10, longitude: 10},         {latitude: 20, longitude: 10},         {latitude: 20, longitude: -10}     ];     console.log('Test 2: Polygon:',polygon);     console.log('should be true', geolib.isPointInside({ latitude:15, longitude:-5}, polygon))     console.log('should be false', geolib.isPointInside({ latitude:15, longitude:-15}, polygon))     console.log('should be true', geolib.isPointInside({ latitude:15, longitude:-5}, polygon))          // ----------     // Test 3 : AntiMeridian 1: This fails     // ----------     var polygon = [         {latitude: 10, longitude: 175},         {latitude: 10, longitude: -175},         {latitude: 20, longitude: -175},         {latitude: 20, longitude: 175}     ];     console.log('Test 3: Antimeridian 1: Polygon:',polygon);     console.log('should be true', geolib.isPointInside({ latitude:15, longitude:177}, polygon))     console.log('should be false', geolib.isPointInside({ latitude:15, longitude:-170}, polygon))     console.log('should be true', geolib.isPointInside({ latitude:15, longitude:-176}, polygon))     console.log('------------');          // ----------     // Test 4: Antimeridian 2: This fails     // ----------     var polygon = [         {latitude: 10, longitude: -175},         {latitude: 10, longitude: 175},         {latitude: 20, longitude: 175},         {latitude: 20, longitude: -175}     ];     console.log('Test 4: Antimeridian 2: Polygon:',polygon);     console.log('should be true', geolib.isPointInside({ latitude:15, longitude:177}, polygon))     console.log('should be false', geolib.isPointInside({ latitude:15, longitude:-170}, polygon))     console.log('should be true', geolib.isPointInside({ latitude:15, longitude:-176}, polygon))     console.log('------------');


title: Hotkey in ArcGIS for Desktop to delete layer from Table of Contents?

text: I'm using ArcGIS 10.2 for Desktop and wish to remove/delete a layer from the Table of Contents. I cannot find such hotkey. I have tried the `DELETE` button but that does not do it. Also such hotkey is missing in the ArcGIS Help.


title: Is it possible to import an autocad feature to a feature layer without exporting to .shp/.dwg

text: Does anyone know if it is possible to import a feature (point, line or polygon) drawn in Autocad directly to a feature layer in ArcMap, so in essence a "feature sucker", rather then exporting the feature to a shp or dwg then importing into ArcMap then copying the feature into the correct feature class. EDIT: The datatype that I want to target is a feature class on ArcSDE (9.2). Versions are ArcGIS 9.3.1 and AutoCAD 2011. Thanks.


title: Google Maps/Google Earth - KML Overlays and Custom Orientation

text: I am am trying to get the floor plans for our various buildings imported into 4 separate KML files, one for each floor level so that I can toggle them with javascript. I'm able to overlay our floor plan images and properly scale them in Google earth, but once exported the KML doesn't display properly. Please take note that whatever the orientation of the Google map is, I can't get my overlays to display properly. On the actual API based map, I have it set to be oriented with West at the top instead of North, which seemingly throws us another curve ball, as when I attempt to load the KML per the API's docs and refresh the API based map, I end up with nothing but "There are no images for this zoom level." My KML:               <?xml version="1.0" encoding="UTF-8"?>     <kml xmlns="http://www.opengis.net/kml/2.2" xmlns:gx="http://www.google.com/kml/ext/2.2"     xmlns:kml="http://www.opengis.net/kml/2.2" xmlns:atom="http://www.w3.org/2005/Atom">     <GroundOverlay>     <name>Untitled Image Overlay</name>     <Icon>         <href>http://img694.imageshack.us/img694/2600/walther.png</href>         <viewBoundScale>0.75</viewBoundScale>     </Icon>     <LatLonBox>         <north>37.21924494476235</north>         <south>37.21876115002708</south>         <east>-93.2637540015874</east>         <west>-93.26445522810076</west>         <rotation>87.56327787656277</rotation>     </LatLonBox>      My API v3 Javascript:               <script type="text/javascript">         function initialize() {             var options = {                 center: new google.maps.LatLng(37.22090,-93.26371),                 zoom: 18,                 minZoom: 18,                 maxZoom: 21,                 heading:-90,                 mapTypeId: google.maps.MapTypeId.SATELLITE             };             var map = new google.maps.Map(document.getElementById("map_canvas"), options);                  //here's where it breaks!!!             //it just won't load the kml...or it is and just breaks beyond reason!             var ctaLayer = new google.maps.KmlLayer('test.kml');             ctaLayer.setMap(map);         }       </script>      Here's what it looks like in Google earth (correct): http://img811.imageshack.us/img811/6348/gearth.jpg Here's what it looks like in Google maps (for testing only, using it to load the KML...until I also figure out how to get my API based map to even display the layer without breaking): http://img838.imageshack.us/img838/2237/gmapsr.jpg Does anyone here have any idea what I need to do to get the display working properly, with regular google maps and more importantly my API based google map?


title: Plotting a circle on a Web map to view the resulting coverage over a central location

text: Not exactly a GIS question but I hope no one will mind. Does anyone know of a website (something equivalent to Google Map) where you can indicate a location, give a radius and you can view the resulting coverage over that central location ? Thanks.


title: 3D Analyst missing from Tool Boxes in ArcGIS 9.3

text: I have a missing 3D Analyst from ArcGIS 9.3 (ArcMap - ArcInfo). I need it for a task to be done. How can i install it? Does anybody know it? Should i donwload the extension or is there any other solution?


title: Looking for an Android Object-relational mapping programming layer, PostGIS oriented

text: I'm working on a Android app which should manage data on a remote PostGIS database. Before diving into actual programming, I'm searching for ORMs for Android. I've found a couple, namely androrm and ormlite. However, these don't seem to give me an interface to query objects by the geometry column. For example, find all objects within a radius from the current location... I'd need to write SQL for this... What are my best options here? Is there a way to completely avoid PostGIS SQL? Or not really? Thanks in advance. (EDIT, 19/08/11: changed the title)


title: How to get an array of simple lines from one multipoint linestring in PostGIS?

text: How can I extract some LINESTRING consisting of 3 or more POINTs to several LINESTRINGs each of 2 POINTs in PostGIS?


title: Explanations for Various GIS Job Titles

text: I've been out of school and working in the GIS field (mostly the field of a desk) for over five years now. The one consistency I've found when working for various companies and applying for various jobs is that there is no consistency in the way our job titles correlate to the actual tasks we accomplish on a daily basis. Does anyone have a good list of the general tasks / job description for some of the GIS job titles:   * GIS Technican   * GIS Analyst   * GIS Administrator   * GIS Developer   * GIS Applications Specialist   * GIS Professional   * GIS Manager   * GIS / IT Technican   * GIS Technical Architect I know each title varies by region and company, but a general list would be a great resource. If I have left any out, please feel free to add.


title: Easiest way to create historical street maps using QGIS

text: I want to make some historical street maps using QGIS. I understand how to upload a historical OS map from DigiMaps to QGIS, but I want to actually 'draw' the streets myself in order to have some detailed maps. The area I am dealing with was subject to a slum clearance in the 1960s and I am looking to map it how it was in the late-nineteenth and early-twentieth centuries. Could anyone offer advice? Thanks in advance!


title: Error: already exists.(stg_e_filealreadyexists) when using IMapDocument.Open method?

text: I got the errror: already exists.(stg_e_filealreadyexists) when using the method IMapDocument.Open() as below codes:               if (File.Exists (strDirectory + "\\temp.mxd"))         File.Delete(strDirectory + "\\temp.mxd");     System.IO.FileStream oFileStream = new System.IO.FileStream(strDirectory + "\\temp.mxd", FileMode.Create);     oFileStream.Write(MyImage, 0, MyImage.Length);     IMapDocument mapDoc = new MapDocumentClass();     mapDoc.Open(strDirectory + "\\temp.mxd"); // Throw exception      Do you know why and how to fix this error ? Thanks and regards, Tai


title: how to transfer from traditional GIS software developing to internet area

text: I am a GIS developer, formerly doing some secondary development with GIS software. I now want to obtain some internet-based GIS area in my career. Somebody advised me to first read some networking/linux books. When prepared, then tranfer to some company to provide experience and so on. My question is. Is there any opensource project belonging to this type: server side provide some service, for example, provide some navigation service(or other GIS area service). I want to learn something about these software platforms and how to realize the functionality provided by server to client, Scheme about the whole thing. Ideall6 such type of thing is using in c++/c. Easier for me to learn. My former area and expertise has been using desktop software. Does anyone have any advices? thx.


title: Is it possible to draw a point interactively and execute an intersection operation (drill down) to identify spatial relationship with other layers?

text: Is it possible to draw a point interactively and execute an intersection operation (drill down) to identify spatial relationship with other layers using GeoServer? I am a newbie with GeoServer…


title: Citing background layer from ArcGIS Online that is used in ArcMap

text: How do you cite a background map layer that is used in ArcMap that is retrieved from ArcGIS Online? For example the background layer is suppose to be credited from the U.S. Geological Survey and in the description it was digitized from the U.S. Geological Survey Professional Paper 1183.


title: Looking for a list of tools/add-ins for ArcGIS Explorer

text: I'm hoping to compile a list of add-ins & resources for ArcGIS Explorer. I intend to install it at home, and apparently, there are a number of server feeds out there that can significantly augment the capabilities of this **free** ESRI software (IMS/WMS servers, toolbox items, etc, etc). I thought many of you might be interested in having access to such a compilation as well, especially since it's a free software. I'm quite anxious to see what's out there!! Try your best to give at least a brief description of the add-in/server feed you post in your answers. Thanks everyone


title: How to sample a 3D raster along a 3D path (or 4D raster along 4D path)?

text: I have a 3d raster dataset of atmospheric variables (lat x lon x altitude). I would like to sample this raster dataset along a 3d path. Any idea how I can achieve this, preferable with qgis or postgis or rasdaman? **Edit:** Actually, my data is 4d (time x lat x lon x altitude), but having a 3d solution would already be great.


title: QGIS MergeShapes Error

text: I am using the MergeShapes plugin with QGIS 1.8.0 on a Windows 7 computer. The issue I am having is when I perform a merge from a directory based on one geometry type, and then switch the geometry type to process a second time (ie: process all polygons in c:\Temp and then process lines in the same directory). This is a limitation when you need to preform multiple merges based on geometry. Any help would be appreciated. Below is the error from MergeShapes. An error has occured while executing Python code: Traceback (most recent call last): File "C:/Users/ryangarnett/.qgis//python/plugins\mergeshapes\mergeshapesdialog.py", line 141, in accept self.progressFiles.setRange( 0, self.inputFiles.count() ) AttributeError: 'NoneType' object has no attribute 'count' Python version: 2.7.2 (default, Jun 12 2011, 15:08:59) [MSC v.1500 32 bit (Intel)] QGIS version: 1.8.0-Lisboa Lisboa, 6416f38 Python path: ['C:/PROGRA~2/Quantum GIS Lisboa/apps/qgis/./python', 'C:/Users/ryangarnett/.qgis//python', 'C:/Users/ryangarnett/.qgis//python/plugins', 'C:/PROGRA~2/Quantum GIS Lisboa/apps/qgis/./python/plugins', 'C:\PROGRA~2\Quantum GIS Lisboa\bin\python27.zip', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\DLLs', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\plat-win', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\lib-tk', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\qgis\bin', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages\PIL', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages\win32', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages\win32\lib', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages\Pythonwin', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages\wx-2.8-msw- unicode', 'C:\Users\ryangarnett\\.qgis\python\plugins\qgis2kml', 'C:\Users\ryangarnett\\.qgis\python\plugins\qgis2kml', '/usr/share/qgis/python', '~/.qgis/python', '.', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\qgis\python\plugins\fTools\tools', '~/.qgis/python', '.']


title: Proxy Issues for WFS and GeoServer - Running on Local Machine

text: I see that a number of similar questions have been asked, but I cannot find the answer to mine so I will be as specific as I can. I am trying to create an OpenLayers map that can run without an internet connection. I have created my own OpenStreetMap tiles and am storing them locally on my hard drive. I have also downloaded local copies of the OpenLayers.js and OpenStreetMap.js files, etc. The basic map functionality is all working. I am now trying to plot a number of latitude/longitude points on the map. These points are stored in a PostGIS database, and I have loaded them into a layer in GeoServer. I am attempting to create a WFS layer in my OpenLayers map that will plot the data from the GeoServer layer, but as of right now the points are just not showing up on the map. I understand from previous questions that I probably have to deal with Access-Control-Allow-Origin. I have installed Apache Server, downloaded and modified the proxy.cgi file to allow 'localhost' and 'localhost:8080', and placed the proxy.cgi file in the cgi-bin directory. My question is: Do I need to use Apache Server if I am trying to run this all locally? I have not yet been able to get this to work and I imagine it is because I don't know what I am doing with Apache. I appreciate any suggestions you may have. Thanks, -Rob _**_**EDIT ** _*_** Thanks for the suggestions! iant, I have successfully moved everything into the GeoServer www directory, and I can open my map by double-clicking on the .html file. Strangely enough, I get a 404 NOT_FOUND error when I try to use a url like the one you posted. Unfortunately, my WFS layer still won't show up. In GeoServer, I have the following settings: \- Workspace Name: MyWorkspace \- Workspace URI: myworkspaceuri \- Layer Name: testpoints \- Native SRS: EPSG:4326 I am including my code below:               // Start position for the map (hardcoded here for simplicity)         var lat=42.3583;         var lon=-71.0603;         var zoom=10;              var map; //complex object of type OpenLayers.Map              //Initialise the 'map' object         function init() {                  // allow testing of specific renderers via "?renderer=Canvas", etc             var renderer = OpenLayers.Util.getParameters(window.location.href).renderer;             renderer = (renderer) ? [renderer] : OpenLayers.Layer.Vector.prototype.renderers;                  map = new OpenLayers.Map ("map", {                 controls:[                     new OpenLayers.Control.Navigation(),                     new OpenLayers.Control.PanZoomBar(),                     new OpenLayers.Control.Permalink(),                     new OpenLayers.Control.ScaleLine({geodesic: true}),                     new OpenLayers.Control.Permalink('permalink'),                     new OpenLayers.Control.MousePosition(),                                         new OpenLayers.Control.Attribution()],                      maxExtent: new OpenLayers.Bounds(-8182919, 5035362.5, -7784059.5, 5304535),                 maxResolution: 700,                 numZoomLevels: 18,                 units: 'm',                 projection: new OpenLayers.Projection("EPSG:900913"),                 displayProjection: new OpenLayers.Projection("EPSG:4326"),                      layers: [new OpenLayers.Layer.OSM(                             "Local Tiles",                              "openlayers/tiles/${z}/${x}/${y}.png",                              {numZoomLevels: 18,                              alpha: true,                              isBaseLayer: true,                             //tileOptions: {crossOriginKeyword: null}                             }                         ),                              new OpenLayers.Layer.Vector("Points", {                             strategies: [new OpenLayers.Strategy.BBOX()],                             protocol: new OpenLayers.Protocol.WFS({                                 version: "1.1.0",                                 url: "localhost:8080/geoserver/wfs",                                 srsName: "EPSG:4326",                                 featureNS: "localhost:8080/workspaceuri",                                 featurePrefix: "MyWorkspace",                                 featureType: "testpoints"                             }),                             renderers: renderer                          })                          ]                       } );                           var switcherControl = new OpenLayers.Control.LayerSwitcher();                 map.addControl(switcherControl);                 switcherControl.maximizeControl();                  if( ! map.getCenter() ){                 var lonLat = new OpenLayers.LonLat(lon, lat).transform(new OpenLayers.Projection("EPSG:4326"), map.getProjectionObject());                 map.setCenter (lonLat, zoom);             }         }


title: How to enter metes and bounds descriptions in ArcMap with only a Basic license?

text: I am using ArcGIS for Desktop 10.1 Basic without extensions. How can I convert a metes and bounds description such as below? For example, 47 degrees, is there a way you can convert the degrees minutes seconds while using CONTROL G (Direction / Length)? I switched in the Editing Options the direction type to Quadrant Bearing and Direction Units to DMS.               SW, SOUTH 0 DEGREES 02' WEST, 110.51 FEET;      SOUTH 47 DEGREES 52' WEST, 246 FEET;      SOUTH 74 DEGREES 41' WEST, 167.25 FEET;      NORTH 75 DEGREES 19' WEST, 274.85 FEET;      NORTH 41 DEGREES 49' WEST, 339.45 FEET;      NORTH 89 DEGREES 14' WEST, 42.60 FEET;      AND EAST 876.75 FEET TO BEGINNING, CONTAINING APPROXIMATELY 4 ACRES,LESS 4 ACRES IN THE NE CORNER OF SAID SW/4


title: How to draw a boundary feature of an area?

text: Which shape file is more appropriate for drawing a boundary of an area. I mean would I take polyline feature or polygon feature for drawing a boundary? I want to digitize a map of metropolitan area. First I have made a manual map and now I am trying to digitize.In this map there are moujas, thanas, rivers, roads, different important location such a market place, police station. So in this case, there is overlapping of borders. Because a thana border can be a district border at the same time. So In this case which feature I chose.


title: What is a simple way to create a popup table on map

text: Using arcgis javascript, is there a way to create a popup table after i perform some processing? Basically my user will select an extend and all address that fall within the extend will be shown.


title: Select empty Postgis (Postgresql) geometry column

text: I try to select rows in Postgis that have an empty geometry column. This works but gives me the opposite:               SELECT * FROM foo where geom <> '';      And this fails:               SELECT * FROM foo where geom = '';          >>>ERROR: parse error - invalid geometry     >>>SQL state: XX000     >>>Hint: You must specify a valid OGC WKT geometry type such as POINT, LINESTRING or POLYGON      And this returns nothing:               SELECT * FROM foo where ST_IsEmpty(geom);      Does anyone know the correct syntax?


title: How to remove 'Edit Feature' option in 'View Map' task of ArcGIS Mobile?

text: We have a Windows Mobile application which is based on ArcGIS mobile's out of the box tasks (version 10.0). We are using 'View Map' task directly. These are the options available to 10.0 ![enter image description here](http://i.stack.imgur.com/9XSfh.png) While upgrading to version 10.1 we have this additional 'Edit Features' which is not desirable for our workflow. ![enter image description here](http://i.stack.imgur.com/nQQya.png) So I want to get rid of this Edit Feature option. Is there a configuration file where I can get rid of it? so technically I am asking to control the behavior of this : {"assemblyQualifiedName":"ESRI.ArcGIS.Mobile.Client.Tasks.ViewMap.ViewMapTask, ESRI.ArcGIS.Mobile.Client","name":"View Map","description":"Browse the map display"}, \- thanks


title: In cartodb, how can I insert a new row in a table via API?

text: We've discovered cartodb and we want to use it in our project, but in the API documentation we've only seen the sql.execute() method (docs) in order to insert a now row in a table. Is this the way to insert points dynamically in our map?