title: Question about geotiffs

text: My colleague and I are georeferencing a bunch of .tif images. After each image, we'll click "update georeferencing". The resultant .tif files are different between our computers and I'm trying to learn why. If I use my colleague's computer to georeference the file, I can keep the georeferenced .tif and delete the .tfwx, tif.aux.xml, and .tif.ovr files, and the .tif will still appear in the correct location. ArcGIS recognizes the correct projection. But if I georeference the file on my computer and delete the .tfwx etc. files, the georeferencing information from my file is completely lost and it does not appear in the correct location. The coordinate system information is lost on my .tif file. I'm flabbergasted about why our computers are performing this task differently. I've followed the exact same steps on the exact same files on both computers. Can you offer any advice? We're both using ArcGIS 10. I've run gdalinfo on the various geotiffs, and the results look identical between the two resultant .tif files (right two panels). See gdalinfo results here - https://dl.dropbox.com/u/17669584/error%20messages/stav_vs_matt_computer.JPG.


title: Trouble Displaying GeoJSON file in OpenLayers

text: I am noobie to OpenLayers and can not for the life of me figure out how to display a geoJson file. This seems like it should be easy but I have been banging my head for a day. I'm sure I am missing something obvious. I have my geoJSON file (countries.json) on a server and basically followed the instructions in "OpenLayers:beginners guide" word for word so I dont know why the thing wont display. Here is my code:               function init(){         map = new OpenLayers.Map('map_element', {});              var vector_strategies = [new OpenLayers.Strategy.Fixed()];         var vector_format = new OpenLayers.Format.GeoJSON({             externalProjection: new OpenLayers.Projection("EPSG:4326"),             internalProject: new OpenLayers.Projection("EPSG:900913")              });         var vector_protocol = new OpenLayers.Protocol.HTTP({                  url: 'ne_50m_admin_WM.json',             format: vector_format         });              vector_layer = new OpenLayers.Layer.Vector('Countries',{             protocol:  vector_protocol,             strategies: vector_strategies         });              map.addLayer(vector_layer);         if(!map.getCenter()){         map.zoomToMaxExtent();      I keep getting thrown this error: _mapProjection is null_ Does anybody have any ideas what may be going astray?


title: Calculating distance from points to nearest "eastern" polygon

text: I am using ArcGIS 10, and trying to do the following: I have a layer file of points (industrial facilities), and I would like to calculate the distance from each point to the nearest "eastern" polygon (U.S. state). I know there are tools for calculating distances from points to polygons, but I don't know how (or, if you can) specify direction of the polygon. Any advice? Thanks in advance.


title: Map Service Update - clearing browser temp files

text: I'm noticing whenever I make a change to an ArcGIS map service (add layer, remove layer...etc) I always have to clear server cache and my internet browser temp files in order for my flex viewer app to not give an error. This can be problematic for any user viewing the application before and after an update has been made to the service. What is the strategy to prevent these errors from occurring on users accessing an app that may need periodic updates? Changing the application url (I'm guessing there must be something better than this)?


title: problem with search (formPanel)

text: am new to GeoEXT and am trying to do the search with form panel. and i want to search features from 2 postgis table, published in geoserver, this features correspond to the field on the formpanel and show them on two separated grid panels, when i hit the search button i get nothing on my gridpanel :s here's my java script code:                     // define the data source                  var protocol= OpenLayers.Protocol.WFS({             url: "http://localhost:8080/geoserver/wfs",             featureType: "titres_salemed_valid",             featureNS: "http://geoserver.org/ancfcc"         });              var protocol2= OpenLayers.Protocol.WFS({             url: "http://localhost:8080/geoserver/wfs",             featureType: "bornes_salemed_valid",             featureNS: "http://geoserver.org/ancfcc"         });                  formPanel = new Ext.form.FormPanel({             title: "Recherche T/R",             height: 200,             region: "north",             items: [{                 xtype: "textfield",                 name: "nature__eq",                 value: "T",                 fieldLabel: "Nature"             },{                xtype: "textfield",                name: "indice__eq",                value: "20",                fieldLabel: "Indice"             },{                xtype: "textfield",                name: "num__eq",                value: "",                fieldLabel: "Num:"             },{                xtype: "textfield",                name: "complement__eq",                value: "",                fieldLabel: "Complement:"               }],              // here action on search button               buttons: [{text: 'Localiser',             handler: function(){        // here to populate the store 2 with feature from table 2 defined by protocol 2                              formPanel.getForm().doAction(searchAction2, {                         callback: function(response) {                                              features2 = searchAction2.response.features;                                         store2.loadData(features2);                               };                               };     // here to populate the store with feature from table  defined by protocol                                formPanel.getForm().doAction(searchAction2, {                         callback: function(response) {                                              features2 = searchAction2.response.features;                                         store2.loadData(features2);          };          };                  }         }]          });                var searchAction = new GeoExt.form.SearchAction(formPanel.getForm(), {         protocol: protocol,         abortPrevious: true     });                  var searchAction2 = new GeoExt.form.SearchAction(formPanel.getForm(), {         protocol: protocol2,         abortPrevious: true     });                var cols = [ {name: 'titref', type: 'string'},                  {name: 'mappe', type: 'string'},                  {name: 'type', type: 'string'}                ];              var reader = new GeoExt.data.FeatureReader({},cols);              var store = new GeoExt.data.FeatureStore({             reader: reader,             fields: cols,             autoLoad: false         });                   // pour les titres               var cols2 = [ {name: 'name', type: 'string'},                  {name: 'x', type: 'string'},                  {name: 'y', type: 'string'}                ];              var reader2 = new GeoExt.data.FeatureReader({},cols2);              var store2 = new GeoExt.data.FeatureStore({             reader: reader2,             fields: cols2,             autoLoad: false         });       // create grid panel configured with feature store         gridPanel = new Ext.grid.GridPanel({             title: "results",             region: "center",             store: store2,             height: 100,             columns: [{                 header: "titref",                 width: 200,                 sortable: true,                 dataIndex: "titref"             }, {                 header: "mappe",                 width: 60,                 sortable: true,                  dataIndex: "mappe"             }, {                 header: "type",                 width: 60,                 align: 'right',                 sortable: true,                 dataIndex: "type"             }],             sm: select         });              // create grid panel configured with feature store for bornes layer         gridPanel2 = new Ext.grid.GridPanel({             title: "tableau des bornes",             region: "south",             store: store2,             height: 250,             columns: [{                 header: "name",                 width: 100,                 sortable: true,                 dataIndex: "name"             }, {                 header: "x",                 width: 120,                 sortable: true,                 dataIndex: "x"             }, {                 header: "y",                 width: 120,                 align: 'right',                 sortable: true,                 dataIndex: "y"             }],             sm: select2         });


title: Dynamic labeling looks awful in flex map

text: I published a service to ArcGIS Server 10.1 from an .mxd with no halos around the labels and created a flex application using the flex viewer. When I browse to the map here and zoom in on an individual trail segment, the labels become very grainy and appear to have a white halo around them, making them very difficult to read. It becomes exceptionally noticeable when imagery is used as a basemap. Is there any thing I can do to stop the labels from getting all messy looking or is this just how they look when using flex?


title: Pictometry Accuracy and/or Alternative

text: I'm looking for some insights into the accuracy (x,y,z) of oblique imagery from Pictometry in terms of providing accurate measurements for property evaluation. Is the high cost of this product justified by its potential to provide not only high visual/pixel resolution, but also highly accurate measurements? Essentially, I'm curious about alternatives that are more cost effective and provide equal or better accuracy in the 2-D plane. I find a lot of sales/marketing info from Pictometry, but not much quantitative metadata in terms of consistency for what DTMs they use to rectify the images. There seems to be a lot of momentum behind this company and I want to know if it's justified.


title: request two layers in search action of one PanelForm

text: here's my case: i have a panelForm with 4 fields, and a search button, that search for the features that correspond to the 4 fields given by the user. and i have two layers that both contains the 4 fields. after the search i have two gridpanel, each one hold data(attributes) from each layer corresponding to the found feature, now i only get the features of only one gridPanel, because the panel form accept only one protocol. Plz how can i get the other attribute corresponding to the founded feature from the other layer. thnks in advance


title: Tabulate Area crashes Python and ArcGIS

text: I have a script that runs `Tabulate Area` on a single feature class multiple times with different rasters. In total `Tabulate Area` runs 5 times. 2 of the input rasters are provided as arguments to the script and `Tabulate Area` works as expected. The other 3 rasters are results of raster functions (i.e. `Reclassify, SetNull, EucDistance, etc.`). Every time I run the script, Python crashes when trying to run `Tabulate Area` on one of those preprocessed rasters. It happens on a different raster every time, seemingly randomly. There is no error message, Python simply crashes. When I hit `Get Details` in the generic Windows error dialog, it seems like there is a different Fault Module every time as well. When I try to run `Tabulate Area` manually in ArcCatalog, ArcCatalog crashes and asks if I want to send a report to ESRI. Again, no error message. **UPDATE:** Here is some of the code involved:               def GetBuffer(inRaster, lcClass, bufDist, ringBuffer=True):              setnull = SetNull(inRaster, 1, "Value <> " + str(lcClass))              dist = EucDistance(setnull, bufDist)              buf = Reclassify(dist, "Value", RemapRange([[0, 100, 1], ["NODATA", 0]]))              if ringBuffer:                  reclass = Reclassify(setnull, "Value", RemapValue([[0, 0], [1, 1], ["NODATA", 0]]))                  ring = arcpy.gp.Minus_sa(buf, reclass)                  return ring              else: return buf          def SubtractRaster(inRaster1, inRaster2):              sub = arcpy.gp.Minus_sa(inRaster1, inRaster2)              setnull = SetNull(sub, 1, "Value <> 1")              return setnull          lcRaster = arcpy.Raster(landcover)          treecanopy = GetBuffer(lcRaster, TreeCanopyClass, TreeCanopyBuffer, False)          parkinglots = SubtractRaster(GetBuffer(lcRaster, ParkingLotClass, ParkingLotBuffer), treecanopy)          buildings = SubtractRaster(GetBuffer(lcRaster, BuildingClass, BuildingBuffer), treecanopy)          impervious = SubtractRaster(GetBuffer(lcRaster, ImperviousClass, ImperviousBuffer), treecanopy)          ...          TabulateArea(fc, "GID", parkinglots, "Value", plTabArea , 1)          TabulateArea(fc, "GID", buildings, "Value", buildTabArea , 1)          TabulateArea(fc, "GID", impervious, "Value", impTabArea , 1)      **UPDATE 2:** I installed ArcGIS 10.1 Service Pack 1 after seeing bugs related to `Tabulate Area` in 10.1. When I first ran the script with SP1, it worked. When I ran it again later, I had the same issue where Python simply crashed without providing an error message.


title: Simple article to explain how GIS work?

text: I'm just a user of GPS and GoogleMaps/OpenStreetMap. I don't have the maths to understand how they work, but am curious to learn. Is there a good, newbie-accessible article that explains in simple terms how GIS work? Thank you.


title: How to change the language setting/interface of Quantum GIS (to English) using Windows XP in French?

text: I want QGIS interface in English in order to share the information with the GIS community. Using other language than English when chating with the GIS community is very difficulte. When I installed QGIS, it automatically set the language to French. Very frustrating to not have the opportunity to select my own language. I am using Windows XP in French but want Quantum GIS in English. Thanks in advance! Doum


title: how to correct the coordinate of points

text: I have some a layer which contains some points. The coordinate of the points should be GPS points. However, the coordinate of these points are not correct. What should I do? I am using QGIS. If I use Coordinate Reference System WGS84/EPSG4326, the coordinate of these points are around (178.xxx, 0.3xxxx). However, the coordinate of these points should be around (103.xxxx, 1.3xxxx). The files are .shp files. The files should use WGS84 coordinate reference system.


title: CQL Filter on WMS

text: I have CQL_Filter which filters WMS layer and displays only the points I asked for, The problem is that all features even the hidden, keep showing its popup when I click on their spot. So, as I understand, I should apply filter also to my protocol or the control: **1- for the protocol:** Coz I use Protocol for selecting like so:               //select is empty vector layer      var selectControl = new OpenLayers.Control.GetFeature({              protocol: OpenLayers.Protocol.WFS.fromWMSLayer(wms, {                  url: "http://IP/geoserver/wfs",                  featureType: "wms"              }),              hover: false          });         selectControl.events.register("featureselected", this, function(e) {               select.addFeatures([e.feature]);                 createPopup(e.feature);                   else              e.feature.destroy();               });       **2-Or on the control:**               Var info=new OpenLayers.Control.WMSGetFeatureInfo({           layers: [wms],      autoActivate: true,           maxFeatures: 3,          infoFormat: "application/vnd.ogc.gml",       eventListeners: {              "getfeatureinfo": function(e) {      var items = [];                  Ext.each(e.features, function(feature) {                     items.push({      source: feature.attributes                      });                      createPoup(feature);…       I have no idea how to accomplish that. Any help? Thanks in advanced. **EDIT**               wms.mergeNewParams({'CQL_FILTER': "type LIKE '%bird%'"});      I already have this line,it filtered the view only and returned all features.


title: OpenLayers DrawFeature override double click

text: Can anyone explain how I can ovveride the dblclick event when I'm using an OpenLayers.Control.DrawFeature? I see there is parameter in the constructor for handlerOptions - Used to set non-default properties on the control’s handler. My handler is just an OpenLayers.Handler.Path. I want to run some custom validation when a user double clicks so that they can continue the sketch should they fail. Double click --> "you've drawn something stupid are you sure you want to finish?" --> No ---> carry on sketching Thanks, James


title: KML parsing error: parse error line 1, column 0. Not well-formed, invalid token

text: A week ago I opened an old .kmz file in google earth to add a few pin-drops. Everything saved correctly and updated for a good 3 days, then I started to get a parsing error in google earth upon opening the file. No idea what this is or how it developed its self, but I know absolutely nothing about how to fix these, I know nothing about programming or coding, and I have tried previous recommendations on the google groups help page with no progress made (https://groups.google.com/forum/#!forum/kml-support-advanced). I was wondering if there was a detailed list of instructions on how to repair a file?


title: Installing PostGis / GeoServer in Mac OSX

text: On a Mac Server I am working to setup the ports. My understanding is that there is a system instance of PostgreSQL running and then a geoserver instance. Version is 10.7.5. I need help approaching how to configure ports in this server instance. This link was used as a reference - http://suite.opengeo.org/opengeo- docs/installation/mac/misc.html Numerous config files are referenced, but even logged in as root, the file paths don't seem correct. I am working to allow remote connections to the server from the local area network, so am deducing I need to change port numbers. I want to make sure that I don't inadvertantly change system Postgresql ports, but rather the postgresql instance for the geoserver. Hoped an easier install.


title: Tilemill 2 (not v1): How to set map extents or bounds for export?

text: I'm using Tilemill 2 and trying to set bounds or extents in order to reduce file size but I can't figure where to do that. I'm also getting a 422 error when trying to upload to my mapbox account and am guessing that is a file size error but I don't actually know.


title: DD83 to State Plane conversion

text: I am having an elementary issue with a conversion I have completed successfully in the past in ArcMap10. I am trying to convert DD83 to State Plane 83 using Data Management Tools > Feature > Project. I am trying to convert coordinates in a shapefile I created from an Excel file I added and projected. Is there something I am missing with the formatting of my original Excel file. When I complete the State Plane conversion, the new shapefile only shows the original DD83. Any help would be appreciated. Thanks.


title: Aerial Photograph Footprint Size Calculation

text: I was wondering what an easy formula to calculate the dimensions of the ground covered by an aerial photograph. I know the following: camera sensor size (mm x mm) camera focal length (mm) height above ground (m) How can I combine these variables to calculate metres by metres covered by each individual image captured? Thanks, Jarrett


title: Extract Raster Values by Polyline End Points

text: I am looking for a way to extract elevations of end points for line data in ArcMap 10.2. I have a road dataset that I want to use as a Network Dataset in Network Analyst. I want to be able to set the elevation change from the beginning of the line to the end of the line. So I can implement the elevation similar to what they have here in step 20. The data available: Roads Dataset, 5m DEM Thanks!


title: How to create linestrings with a definite angle and length, that are fixed to a point

text: I want to create linestrings with a definite angle (e.g. 160°) and length (e.g. 2m) that are fixed to a number of points of another linestring. So, I want to use the ST_DumpPoints function to find the points and bind the created linestrings to them. Is there a way to declare an angle (α) during the linestring creation? Here is an example picture: ![enter image description here](http://i.stack.imgur.com/7krOR.png) I want to create the blue linestrings. **EDIT** The angle (α2) in the picture isn't really exemplary. But rather an azimuth of 160° (like α1). **UPDATE** The answer from Evil Genius helped me calculating the maximum width of a polygon with a given aspect.


title: How to derive acquisition date from short unsigned file using ENVI

text: I need to derive a processing date from a file which is unsigned short file using ENVI software. The files are in ENVI standard format (image and header file format). I have searched this on Google as well but so far unable to figure it out. Any help would be greatly Appreciated, Cheers! Jibran


title: Select by Attributes - If string has a letter

text: I'm running ArcMap 10.1. I'm using a File Geodatabase. I have a string coordinate field that I want to clean and some records have letters in it (mostly 'N' or 'E' or 'n' or 'e' e.g. '6213210N' or 'N6584312' but also some refer to map a grid reference e.g 'm29'). I want to select the records that have a letter regardless of it's position or what letter. Why wouldn't this work?               UPPER("CoordinateField") LIKE '%[A-Z]%'      As far as I can see it works in for similar solutions but not in ArcGIS. As seen **here**. Am I using incorrect notation for something stored in a File Geodatabase? Hoping someone can shed some light on this. Much appreciated.


title: Publish ESRI Feature Services from spatial view on ArcGIS Server 10.1

text: I have a feature class layer and a handful of non spatial plain tables in a 10.1 SDE geodatabase on MS SQL Server. What I want to achieve is to publish a feature service from ArcGIS Server so as for users to edit both features and attributes from the non spatial tables. I tried to make a spatial view by joining the feature table with other non- spatial tables by referencing foreign keys. I was able to publish the spatial view as feature service but not to make any edits from the REST service (Add, Update, Delete). Just wondering why this was not successful and if a spatial view cannot be edited through feature service? Any workaround to achieve my task? Thanks in advance. Cheers, Alex


title: Optimize ST_CoveredBy performance

text: I have a PostGIS 2.0.3 table with 113,862 records with multipolygons represented both as geometry and geography. There is a gist index on each of them. I'm running the following queries:               SELECT *, ST_AsText(geog) As wkt     FROM ds_pg_data_noise_poly WHERE      ST_CoveredBy(ST_PointFromText('POINT(12.432404 55.697712)', 4326), geog)      LIMIT 100;      The above query use Geography in a WSG84 coordinate system and returns the correct 8 rows in 33872 ms.               SELECT *, ST_AsText(geom) As wkt      FROM ds_pg_data_noise_poly WHERE      ST_CoveredBy(ST_GeomFromText('POINT(715692.377 6177776.241)'), geom)     LIMIT 100;      The above query use Geometry in a UTM32 coordinate system and returns the correct 8 rows in 10494 ms. Is there anyway I can optimize the speed of this query? I believe I read somewhere that I could use bounding boxes to decrease the amount of records to iterate. Is this correct, and how would that look? Any suggestions are welcome. _UPDATE:_ By minimizing the returned columns I managed to reduce the query time to 2168 ms.               SELECT noise_cl, noise_in      FROM ds_pg_data_noise_poly WHERE          ST_CoveredBy(ST_GeomFromText('POINT(715692.377 6177776.241)'), geom);      I still do not find this acceptable, but it is definitely better. As each record is placed within a subset, I can actually reduce the records to iterate even more.               SELECT noise_cl, noise_in      FROM ds_pg_data_noise_poly WHERE     (subset = 'subset1' OR subset = 'subset2') AND      ST_CoveredBy(ST_GeomFromText('POINT(715692.377 6177776.241)'), geom);      The above query returns in 1353 ms. I am still not quite happy with the result. Any suggestions are more than welcome.


title: Openlayers WMS GetFeature Info Pop Up Not Working / Mapserver

text:
I am new to OpenLayers and Mapserver .. I have been trying for days and am not able to make it work I want to make a Pop Up for Get Feature Info but nothing appears when i click any place This is the code for HTML/OL part My OpenLayers Map                                          <link rel="stylesheet" href="style.css" type="text/css">                       <link rel="stylesheet" type="text/css" href="C:/ms4w/Apache/htdocs/GeoExt-1.1/GeoExt/resources/css/popup.css">

                     Mapfile is MAP NAME "QGIS-MAP" SIZE 420 300 UNITS dd IMAGECOLOR 255 255 255 IMAGEQUALITY 95 IMAGETYPE png EXTENT 60.899437 23.702916 77.843079 37.097012 #PROJECTION # "init=epsg:4326" #END PROJECTION "init=epsg:900913" END OUTPUTFORMAT NAME "png8" DRIVER "AGG/PNG8" MIMETYPE "image/png ; mode=8bit" IMAGEMODE RGBA EXTENSION "png" END WEB TEMPLATE getfeatureinfo.html IMAGEPATH 'C:\ms4w\Apache\htdocs\output/' IMAGEURL '/output/'                METADATA       "ows_enable_request" "*"       'ows_title'           'QGIS-MAP'       'ows_onlineresource'  'http://localhost/cgi-bin/mapserv.exe??map=C:/ms4w/Apache/htdocs/getfeatureinfo.map&LAYERS=ALL&MODE=MAP'       'ows_srs'             'EPSG:4326 EPSG:900913'       "wms_include_items" "all"     END      END LAYER NAME 'pak_adm2' TYPE POLYGON DUMP true TEMPLATE getfeatureinfo.html EXTENT 60.899437 23.702916 77.843079 37.097012               CONNECTIONTYPE postgis     CONNECTION "dbname='post_map_open' user='postgres' password='user' host=localhost port=5432 sslmode=disable"     DATA 'geom FROM "pak_adm2" USING UNIQUE gid USING srid=4326'          METADATA       'ows_title' 'pak_adm2'       "wms_enable_request" "*"       "wms_inlude_items" "all"       "ows_featureinfo_textplain_include_items" "*"     END          STATUS ON     TRANSPARENCY 100          PROJECTION        "init=epsg:4326"      END          CLASS        NAME 'pak_adm2'         STYLE          WIDTH 0.91           OUTLINECOLOR 0 0 0          COLOR 0 178 0        END     END        END END


title: How to ensure that Tilemill interprets my postgis query as intended?

text: One of my layers (a postgis query from my local osm2pgsql database) is being translated incorrectly in Tilemill. The query is attempting to select the way, the area of the polygon, and the sport=* tag as follows:               ( SELECT way, sport AS sport, way_area AS area, COALESCE(leisure, "natural", highway, amenity, tourism) AS type       FROM planet_osm_polygon       ORDER BY way_area DESC     ) AS data      This query works as expected in pgadmin and returns the results that I expect. In my project.mml file, this layer is represented as follows: > { "Datasource": { "dbname": "cleveosm2pgsql", "extent": "-20037508.34 > -20037508.34 20037508.34 20037508.34", "geometry_field": "way", "host": > "localhost", "id": "landuse", "key_field": "", "port": "5432", "project": > "osm-bright-imposm", "srs": "+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 > +lon_0=0.0 +x_0=0.0 +y_0=0.0 +k=1.0 +units=m +nadgrids=@null +wktext > +no_defs +over", "table": "( SELECT way, sport, way_area AS area, > COALESCE(leisure, \"natural\", highway, amenity, tourism) AS type\n FROM > planet_osm_polygon\n ORDER BY way_area DESC\n) AS data", "type": "postgis", > "user": "skors" }, "class": "", "geometry": "polygon", "id": "landuse", > "name": "landuse", "srs": "+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 > +lon_0=0.0 +x_0=0.0 +y_0=0.0 +k=1.0 +units=m +nadgrids=@null +wktext > +no_defs +over", "srs-name": "900913", "status": "on", "advanced": {} }, When I open my project in Tilemill, I receive the following error: > Postgis Plugin: ERROR: column "sport" does not exist LINE 1: SELECT > ST_AsBinary("way") AS geom,"sport","type" FROM ( SELE... ^ > > in executeQuery Full sql was: 'SELECT ST_AsBinary("way") AS > geom,"sport","type" FROM ( SELECT way, way_area AS area, COALESCE(landuse, > leisure, "natural", highway, amenity, tourism) AS type FROM > planet_osm_polygon WHERE way_area > 10000 ORDER BY way_area DESC ) AS data > WHERE "way" && ST_SetSRID('BOX3D(-9089279.907446878 > 4999593.14607681,-9030576.269723864 5058296.783799823)'::box3d, 900913)' I have already cleared my cache (within the layer's menu in tilemill) and restarted tilemill. Any Idea why Tilemill is modifying my query to change the sport column from sport to "sport"?


title: converting projections using dotspatial or proj4

text: how can i convert projections using either dotspatial or proj4? i am making a program using sharpmap and dotspatial in VB.net have to convert a EPSG 4030 into a Lambert conformal conic projection and back to WGS84 on the fly also can i convert a system on the fly using a WKT?


title: Automatic matching of two street networks

text: I'm faced with the problem of **matching two street networks**. I have these two differently sourced networks. The links in each network can be segmented, i.e. there is not necessarily a single whole link between two intersections, instead the link can be broken into an arbitrary number of segments with different attributes. **I should transfer the attributes from one network to the features of the other network** and I figured for doing that I will have to match corresponding streets/street segments. My current idea is to densify the links in either network and then do an analysis of the **cumulative distances between vertices in the two networks?** If the nodes of two segments overall have small distances from each other, the two links probably correspond? However, I can see problems if two corresponding road segments are expressed in links of wildly different length in the two networks. That notwithstanding, I'm unsure as to how to best go about implementing a method that would analyse the two networks and maybe also offer a way of visually exploring the quality of the matching. On the plus side I'm familiar (to different degrees) with many platforms ( _ArcGIS, QGIS, SAGA, R, Python, ..._ ) and can install any promising tool. I appreciate any pointers you may be able to give.


title: Recreating A Lat/Long Grid Using a Specified Formula

text: I'm trying to recreate the Gridded Mercury emissions 2000 (0.5 x 0.5 degree grid) designed by the Arctic Monitoring and Assessment Programme. The file is in ASCII format, but without any coordinates. The arithmetic to produce the grid is provided, but I don't quite understand how to create the actual grid with the information provided. For you math junkies this is probably a fairly simple problem. I have access to ArcGIS 10 to complete this work. This is the arithmetic used, and the cellcodes are on the downloadable dataset.               259200 (720 x 360) cells     Z05 Cellcode = (j*1000) + i     j = row number starting at 1 for 90S to 89.5S latitude, to 360 for 89.5N to 90N latitude     i = column number starting at 1 for 180W to 179.5W longitude, to 720 for 179.5E to 180E longitude.     (coordinates represent the center of the gridcell)     The latitude and longitude of the center of a gridcell is given by:     latitude = ((j-181)/2) + 0.25     longitude = ((i-361)/2) + 0.25


title: Error message of "Unable to open input file <D:\waveform\Grass/f_format>: Too many open files"

text: I'm using GRASS 7 under Windows 8.1 x64. Sometimes when using the r.to.rast3 Grass Gives me the error message of "Unable to open input file ... Too many open files". I have reduced the number of opened files to 3, but still receiving the same error message. I think it has something to do with the region set parameters. Did anyone of you has faced such problem before? Thanks.


title: Update column in CartoDB with value from another table

text: I am trying to update the column 'lt_dispatch_level' in my table "inland" with the values from another table "current_obs" that has a corresponding column 'lt_dispatch_level'. However when I run this SQL query: `INSERT INTO inland (lt_dispatch_level, timestamp) SELECT lt_dispatch_level, timestamp FROM current_obs;` I get the values from current_obs added in as a second row, with null values being aded in the first row. I attached an image of what this looks like. ![enter image description here](http://i.stack.imgur.com/Ydgv3.png) What do I need to change in my SQL query to get my values in the first row?


title: Unable to drag and drop current map elements into custom geoprocessing tool?

text: I know there mus be a simple solution to this but for some reason i cant figure it out. I have scripted an arcpy tool and the tool works fine in Arc10.2. The only problem is when I use the tool i cant drag and drop current map elements into the parameters, I have to manually navigate to the file through explorer each time, also there isn't a drop down menu of current map items to select from either? Is it something to do with how i have set the tool parameters or is it a scripting issue?


title: Shorten multiline in QGIS

text: I have a multiline dataset in a spatialite database. How to shorten that multiline based on a specified distance or percentage in QGIS?


title: QGIS 2.2 Openlayers Plugin shift between shapefile and base layer for large format

text: I am attempting to create a 24x36 map of New York City, and overlay it with a shapefile from City Planning. I'm using the Openlayers plugin, and I've tried all of the different types of base street map layers. I have been able to successfully convert the shapefile from a NAD 83 projection (EPSG 2263) to a WGS 84 projection (EPSG 3857), but I'm still suffering from some type of shift whenever I draw the map out on the print composer on a large scale piece of paper. The larger I try to print, the more dramatic the shift the misalignment. ![Aligned base layer and shapefile in browser](http://i.stack.imgur.com/axZtL.jpg) ![Misaligned layers in print composer](http://i.stack.imgur.com/J7fdE.png)


title: Summing the number of intersections with a separate buffer layer

text: I have dabbled with ArcGIS and QGIS a bit, and I consider myself a self-taught beginner with a whole lot to learn. So, I have two layers: 1) census blocks, and 2) point data. My goal is to calculate the sum of the point data within the census blocks to eventually create a choropleth map and to run regression models with the sums of the points as a dependent variable; however, some of the points are directly between blocks and are not truly within them, so the points in polygons tool seems to miss a lot of points. Because I am interested in the effect of these points on the proximal blocks, I had the idea of creating a 15m buffer around the points, and summing up the number of buffers that each census block intersected, but I am at a loss for how to do this in QGIS. So far, I have tried to join attributes by location, intersections, unions, and they all don't seem to get me where I need to be. If anyone has an idea how to do this, I would greatly appreciate. I attached a sample picture of my layout with the intersecting census blocks selected. Please excuse my labeling. ![enter image description here](http://i.stack.imgur.com/jZ9va.png)


title: Finding shortest cumulative route (i.e. best 'hub') between networked points

text: How do I find the shortest cumulative route (i.e. best 'hub') between networked points? I have a dataset where some but not all of the points in the network connect to each other, however it is possible to go from any one point to another via the network. My problem is that these points have different attributes, such that I need to find the shortest cumulative route between the points in each group (i.e. the best 'hub' that overall anyone in those points can get to the quickest/shortest average distance from all the other points in the group), as well as the overall 'hub' for all the points. I'm not sure how to do this in QGIS----much less how to format my dataset for it, since it's not a stereotypical road or transportation network----or whether there is a better alternative with which to accomplish the analysis. The spatial distance between the points is actually irrelevant, because the 'length' component for distance is the number of nodes between points (i.e. going through point B to reach point C from point A is shorter than going through points D and E when traveling between A and B).


title: How to get Z attribute of point using Python and ArcGIS Desktop 9.3?

text: I am making a python code for arcgis 9.3. In this code i create a new point with x,y and z attributes as following :               rowsLine = gp.updatecursor(line_Feature_Class)     rowLine = rowsLine.next()          rowsPoint = gp.InsertCursor(point_Feature_Class)     feat = rowsPoint.NewRow()          while rowLine:              #...              ptObj = gp.CreateObject('Point')         ptObj.x = xPremierPoint         ptObj.y = yPremierPoint         ptObj.z = zPremierPoint         feat.Shape = ptObj         rowsPoint.InsertRow(feat)      I want to add the z value in a field "zvalue" in my "point_Feature_Class". The problem is i can't use an insert cursor (rowsPoint) and an update cursor at the same time. So, at the end of the code, I want to delete the insert cursor, create an upsdate cursor, to loop in my "point_Feature_Class" and get the z attribute of each point and add it in the fiels "zvalue". The code should be something like that :               rows = gp.updatecursor(point_Feature_Class)     row = rows.next()          While row :              Zvalue = #???? row.shape.getz ????         row.SetValue("zvalue", Zvalue)          rows.updateRow(row)     row = rows.Next()      Thanks for your help ;)


title: Inserting newlines in Rectangle Text elements via ArcPy causes overlap?

text: I ran across an issue the other day when I tried to use ArcPython's mapping module to edit rectangle text elements with newlines (\n) in an ArcMap document. Here's what the output looked like: ![enter image description here](http://i.stack.imgur.com/mt2OH.png) Here's the code I used to generate that output. The first column are rectangle text elements Text1, Text2, Text3 going down; the second column are "plain" text elements Text4, Text5, and Text6 going down.               import os     import arcpy          HomeDir = r"C:\Desktop"     arcpy.env.workspace = HomeDir          CurrentMXD = arcpy.mapping.MapDocument(r"C:\Desktop\TextTest.mxd")     OutputFilename = r"C:\Desktop\TextTest.pdf"     if os.path.exists(OutputFilename):         os.remove(OutputFilename)          for TextElement in arcpy.mapping.ListLayoutElements(CurrentMXD, "TEXT_ELEMENT"):         TextElementName = TextElement.name              String1 = "The quick brown fox jumped over the lazy dog.\nShe sells sea shells by the sea shore."         String2 = "The quick brown fox \njumped over the lazy dog.\nShe sells sea shells by the sea shore."         String3 = "The quick brown fox jumped \nover the lazy dog.\nShe sells sea shells by the sea shore."              if TextElementName == "Text1":             TextElement.text = String1         if TextElementName == "Text2":             TextElement.text = String2          if TextElementName == "Text3":             TextElement.text = String3         if TextElementName == "Text4":             TextElement.text = String1         if TextElementName == "Text5":             TextElement.text = String2         if TextElementName == "Text6":             TextElement.text = String3          arcpy.mapping.ExportToPDF(CurrentMXD, OutputFilename)      So far, it looks like the presence of the messed up text depends on whether the line is longer enough to wrap, and whether the line before the newline is longer than the line after the newline. Any ideas about what could be going wrong? Is there a workaround? I could use plain text elements and worry about wrapping lines using Python, but I'm hoping I can figure something out.


title: Using newline/page break tag in ArcGIS Desktop Text_Element in arcpy.mapping?

text: I want to update a Text_Element with a newline or page break using ArcGIS Desktop 10.0. Does anyone know the tag or syntax?


title: Lining up .asc XY points with DXF using QGIS 2.0

text: I've exported XY points (nodes of a polyline with attributes) (.asc) from Star APIC and attempted to import it into QGIS as a delimited text file. I cannot however get it to line up correctly with the corresponding DXF file I exported from the same database. My DXF lines up with aerial photographs, so I am sure it is correctly projected. The CSV should be in the same projection (since it comes from the same database). I am considering using the Affine Transformation plugin to rescale and transform. Is this the best solution?


title: joining polylines on geometry and attribute

text: I have dowloaded an Open street Map shapefile containing streets and their attributes, for my whole country. I have qgis and spatialite_gui and a very basic understanding of SQL The problem I have is that OSM splits roads of the same name according to a variety of factors such as attribute changes , this means our main street is divided into multiple sections. each section has it's own OSM_ID I wish to join these sections back together ( I don't care about the attributes, except the name) based on their proximity and their names. (My software that will consume this data from another api will determine if the intersecting road terminates at the intersection from the linestring, i.e. it's not just for display purposes.) I have tried using the qgis-line dissolve method but dissolving on name alone causes a merge between all streets in the country having the same name. eg there are about 30 "Victoria streets". so I need to merge the streets that have the same name AND have a similar geometry (e.g. within 10-20m ) How/what steps should I take to achieve this, I will then want to convert back to a shapefile which will be uploaded to public data server which others can use. My table is called CentrelinesTBL with fields               osm_id     Highwayname     Speed     Name  ( <--Streetname)     Geometry - looks to be a blob, eg: it appears in the table as "BLOB sz =384 geometry" but if I use blob explore I can see...          {"type":"LineString","coordinates":[[175.4088968683988,     -39.40728494205201,0],[175.4091779512512,-39.407522943693,0],     [175.4107932119636,-39.40868038774494,0]]}


title: Error using 'Join attributes by location' via Python

text: Trying to run a simple script using `qgis:joinbylocation` but getting either errors or it's not processing correctly. I don't want to take a summary of the attributes, just the attributes of the first located feature and I want to keep all records. The _parameter selections_ should be correct but it's the `sum,mean,min,max,median` that's causing problems. I tried removing them but python wasn't too happy. Also (for some reason) python doesn't recognise `mean` or `median` but is fine with the other statistical terms. Any advice please?               ##Test=name     ##File_1=file     ##File_2=file     ##Output=output vector          mean=None     median=None          processing.runalg("qgis:joinbylocation", File_1, File_2, 0, sum,mean,min,max,median, 0, 1, Output)


title: Interpolate Temperature Using DEM

text: I want to "co-interpolate" temperature with DEM. I have SRTM 90m aggregated to 1km and I want to "co-interpolate" my DEM with points with temperature data. How do I do this? I am using ArcGIS 10.


title: How to import code list from .mdb to domain model in Enterprise Architect?

text: I have build an UML model in domain model structure in Enterprise Architect. Now i want to import a couple of code lists to this model. The problem is that the tables are in .mdb ArcGis workspace. I try to export tables from .mdb to xml document and then to import xml to Enterprise Architect- but the particular fields are in esri format type and i need them in enumeration (if the values is fixed) or codeList (if the values may be extended) to fulfill the terms of interoperability (INSPIRE).


title: Why are my labels overlapping in ArcGIS 10?

text: Please excuse what may be a very newbie question. I'm not new to GIS, but just learning ArcGIS. I have been led to believe that the standard labeling engine doesn't allow any overlap: http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#/Positioning_labels_with_the_Standard_Label_Engine/00s800000020000000/ however, I am seeing significant amounts of overlap between labels in just one layer. This is the image I get: ![alt text](http://i.stack.imgur.com/FaRXa.png) This is within a Data Frame with a fixed scale (1:2500) and no reference scale . I wondered whether that had anything to do with it, but playing around with either of these parameters doesn't seem to have any effect. Is there anything I'm missing in the label manager options - I have checked and double checked that the "Place Overlapping Labels" check box is off. When I increase the buffer size, I get fewer labels, but some of them still overlap.


title: WMS + Google Earth + tilted view = alignment failure?

text: I'm working on a custom web-based application that generates spatial data and that includes a web-based map viewer. We have a custom servlet that handles WMS requests (not a complete WMS implementation, just enough for our needs) and uses MapServer/SWIG to generate image data. The web front end uses OpenLayers. So far, so good: everything works and gives us nice-looking interactive maps. We also built a feature that creates a KML file with a `<GroundOverlay>` that refers to our WMS URL -- the user clicks a link and it opens in Google Earth and they can see the data. Again, so far, so good. The problem comes when the user starts to tilt the camera and pan around in Google Earth. Google Earth starts mis-registering our data, at first by just a few meters at a time but eventually by tens of kilometers, if you tilt, zoom in and out, and pan around enough. I did some searches and I found several sites that indicated that Google Earth's support for tilting WMS imagery was poor, but little in the way of specifics. Has anyone else tried this and experienced similar issues? More importantly, how should I approach this problem? Given that our users are very interested in doing exactly this kind of exploration (plot data with topography, tilt the camera down, and fly around in our subject area) and that our users are already pretty committed to using Google Earth as the client, how should we make the data available so that they can view it?


title: Product License not initialized error when trying to open an SDE connection in ArcObjects

text: I am trying to connect to create connection in using ArcObjects using the code below.               IPropertySet propertySet = new PropertySet();     IWorkspaceFactory sdeWsFactory = new SdeWorkspaceFactory();          propertySet.SetProperty("SERVER", server);     propertySet.SetProperty("INSTANCE", instance);     propertySet.SetProperty("DATABASE", database);     propertySet.SetProperty("USER", user);     propertySet.SetProperty("PASSWORD", password);     propertySet.SetProperty("VERSION", version);          IWorkspace WS = sdeWsFactory.Open(propertySet, 0);     IFeatureWorkspace fws = WS as IFeatureWorkspace;      When it tries to execute the Open, I get the error message: The Product License has not been initialized. How do I tell which "product" it is talking about? We are runnuing ArcSDE 9.3.1 SP2.


title: How can I dock ArcToolbox in ArcMap 10 custom template and make it stay?

text: Here is the situation: Because UI Controls cannot be saved to Normal.mxt like they used to I had to create an empty document with all my custom VBA code that is mostly called from the local document-level UI controls. (Until I have time to get proficient with .NET things will have to run this way) All is good although it would have been nice to have Normal.mxt-level UI Controls back... Whenever I need custom functionality i start creating the new document within my Custom Template. One annoying problem I have is that I have to dock ArcToolbox every time i open a document created from the custom template. It will just not stay docked whether I save the custom template with it docked or the Normal.mxt. Normal template documents are fine, ArcToolbox is where i want it, just the custom template is misbehaving. Any ideas? Thanks![alt text](http://i.stack.imgur.com/4di38.jpg)


title: Writing list of feature classes with fields to one csv file per feature class?

text: I have a database with multiple feature classes. Now I extracted a list of feature classes and a list of fields. I figured out how to write the list of feature classes and the corresponding fields to csv. But the csv file gets really long. Instead I would like to write each feature class name with the corresponding field names to a separate csv file. My code looks like this right now:               def csvwriter (csvfile, mylist):         myfile = open(csvfile,'ab')         wr = csv.writer(myfile)         wr.writerow(mylist)          datasets = arcpy.ListDatasets()     for ds in datasets:         dslist = []         dslist.append(ds)         csvwriter (csvfile1,dslist)         for fc in arcpy.ListFeatureClasses("*","",ds):             fclist = []             fclist.append(fc)             csvwriter(csvfile1,fclist)             for field in arcpy.ListFields(fc):                 fieldlist = []                 fieldlist.append(field.name)                 csvwriter(csvfile1,fieldlist)


title: Deleting feature classes in filegeodatabase issue

text: I have filegeodatabase with 20 featureclasses and need to delete few of them, but option DELETE in AcrCatalog if I select this fc, and click right mouse is disabled for the most of fcs in fgdb. Why is the reason for that? Can I fix it? This fgdb is created by xmlscheme of gdb. Best regards,


title: gps point closer to the boundary of the polygon

text: I would like to capture the gps points which are closer the polygon boundary (inside), not point in polygon. I have arcgis arcinfo v9.3. i can do with inner buffer of the polygon and intersection and finding the remaining gps points. Is there any other way ?


title: arcpy join field tool is not licensed

text: I am using ArcGIS and arcpy to process some shapefiles. When I try to use the function arcpy.JoinField_management() the program stops and Python says that the tool is not licensed. This confuses me because the documentation for the tool does not ask for any special licenses. http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//001700000065000000 An extract of my code and the error follow (the code is quite long but I run it ineractively to check it as I go). Any help would be appreciated. Regards, Alex               import arcpy     arcpy.env.workspace = 'C:/Data/'     arcpy.CheckOutExtension('Spatial')     outZSaT = arcpy.sa.ZonalStatisticsAsTable('mb06.shp', 'MB_NUM', 'slope', 'sum_slope.dbf', 'DATA', 'SUM')     arcpy.JoinField_management('mb06.dbf', 'MB_NUM', 'sum_slope.dbf', 'MB_NUM')          Traceback (most recent call last):       File "<pyshell#2>", line 2, in <module>         workspace + '/sum_slope7.dbf', 'MB_NUM')       File "c:\Program Files (x86)\ArcGIS\Desktop10.0\arcpy\arcpy\management.py", line 3667, in JoinField         raise e     ExecuteError: Failed to execute. Parameters are not valid.     ERROR 000824: The tool is not licensed.     Failed to execute (JoinField).


title: How to determine the distribution of pollutants using marine tissue samples?--3D Interpolation

text: I have point data representing tissue sampling performed by NOAA NCCOS for Gulf of Maine waters. I am assuming a temporal resolution of one sampling point per location for one year. I have searched for literature to address how to interpolate this data across water depth, but have yet to come to a conclusion. A simple Kriging shows patterns of overlay, which as a simplified approach, is fine. Are there more refined methods to determine the distribution of pollutants using marine tissue samples? I would prefer methods that have been published or applied. I am using ArcGIS 10.2 and I am open to any open source software. I am also proficient in Python and ArcObjects.


title: Creating a 10-mile radius buffer zone that averages a variable in my shapefile?

text: I´m trying to create a "buffer zone" for every polygon in my shapefile, which will identify a 10 mile radius around the centroid of each polygon, and calculate the AVERAGE "demand for services" (one variable in my shapefile attribute table) for all the surrounding polygons within this radius. All of the relevant polygons are in the same layer/shapefile. Ideally, this average would be calculated in a new field in the attribute table, so that I can symbolize this variable and get a more even "heat map" of demand (which is not as tied to the municipality boundaries.) Can anyone help me out with instructions on how to do this?


title: Feature Distance for each Raster Cell

text: I have a polyline layer of roads and a raster of land ownership. For each raster cell I would like to calculate the distance to the nearest road. I am on ArcMap 10.2 with Spatial Analyst and 3D analyst.


title: adding line breaks to labels in ArcGIS Online

text: I'm using the new labeling functionality in ArcGIS Online. My question is how do I add a line break to the label if I want to use more than one field?


title: How to truncate file geodatabase tables with python?

text: Everything is in the title. I need to truncate file geodatabase tables (1 or all) with python. What is the code for that? Thank you in advance!


title: Deleting features (rows) from shapefile (attribute table) using ModelBuilder?

text: I run into a problem and can't figure out how to solve this one. I use the ModelBuilder and have a shapefile with loads of information. This information needs to be reduced. Therefore I want to select several rows (not all of them) out of several columns by attribute to delete them. So for example, column X contains data with attributes a, b and c. I only want to delete c. Since I use the ModelBuilder, I need a tool from the toolbox which does that (preferably) in one step. Any ideas?


title: How do I add a Natural Earth label layer to a TileMill project?

text: I'm using TileMill from MapBox to generate tiles and export Mapnik xml files. So far so good. I've downloaded the entire Natural Earth dataset, and have so far successfully imported layers for lakes, coastlines, land and rivers. Now I would like to add labels, using the "Physical Labels" dataset. When I attempt to add a layer using the shape file 'ne_10m_geography_regions_points.shp', I get a few dozen little red dots (all of which seem to be on a coastline), but no labels. I've tried combinations of the 'ne_10m_geography_regions_polys.shp' and points, but don't get anything nearly resembling the examples on the Natural Earth download site. I'm a bit new at this, so I'm not sure if the problem is TileMill, the Natural Earth data, or me. I suspect it's me but I don't understand why this layer would be different than the others I've successfully used.


title: Reclassify based in neighbour polygons

text: I would like to do a simple operation in QGIS and I am not sure about which is the best option. I have a polygon shapefile and each one has one class attached. I want to reclassify one of them based in the neighbours. For example, I have a polygon with class street and I want to reclassify to house if one of the polygons around it (or one that contains it) belong to class urbanization. Maybe there is already something similar. Has anybody an idea? Thank you very much


title: How to simplify a polygon down to X amount of points?

text: If I have a polygon that has 100000 points, how can I calculate what the tolerance should be to simplify it down to 1000 points using ST_SimplifyPreserveTopology?


title: Exporting feature class from Geodatabase to shapefile, shortens/abbreviates field values in resultant file?

text: I'm exporting feature class (File Geodatabase) to shapefile, after exporting the field values in resultant shapefile is shortening, like if a field value is Andhra Pradesh then in Feature class changed as a AP. I think that that the field value is 10 character limitation. How we can by- pass this limitation. How I can export feature class with all field values?


title: Starting ArcMap says system clock has been set back more than 24 hours?

text: In order to test my log4net logging daily, I changed my clock to be 2 days later. When I changed back the clock to normal, I cannot start ArcMap, I get a message: > The system clock has been set back more than 24 hours I found a possible solution provided by Esri, but it not work. Is there another solution for it? I am using Wndows 7 and ArcGIS 10.2 for Desktop.


title: Looking for a good solution to have GPS position on KMZ files on Android platform

text: I'm looking for the best solution / app on the android platform to show my GPS position over a KMZ file, possibly offline. Here is what I've tried to do:   1. Google Earth Android app I downloaded KMLZ to Earth app and Google Earth app. Now I can open my KMZ files stored in a Google Drive with Google Earth, but GPS position seems to be not accurate. Also KMZ file doesn't look very well.   1. Maps Engine Android app I manage in converting my KMZ file to Google MapsEngine. You can find how to looking at my answer here. Now the GPS position is ok, vector drawing is ok too, but I need to be always online.   1. Other apps I've also tried many KMZ viewer, but got some problems opening my files on most of them, and none of these apps seems to be very good. Can anybody suggest a better solution/app?


title: Cannot Set TOC View/Tab to Display with VBA?

text: I have some code from 9.3 that would check the CurrentContentsView.ContextItem to ensure that the TOC was in Display (or Drawing Order) view. If not it would set it to Display via:               If Not pDoc.CurrentContentsView.ContextItem = 1 Then         Set pDoc.CurrentContentsView.ContextItem = 1         pDoc.ContentsView(1).Refresh Nothing     End If     pDoc.ContentsView(1).Refresh Nothing      Now that the TOC is a DockableWindow (which I can access) and there are several IContentsViews, this no longer works. It's surprising that this is impossible to find in the forums and the VBA SDK as as brief and cryptic as could be. How do I access the ContextItem (which seems stuck at NULL) to set it? Also, is the fact that "currentContentsView" automatically reverts to a lower-case "c" when the object apparently shouldhave an upper-case "C" significant? Any help would be MUCH appreciated. For some reason, my code runs 3x faster in Display View and I want to automatically set it if the user forgets (or doesn't know). Oh...and yes, I AM aware of VBA's demise and the existence of Python. Noted. Thanks in advance.


title: How to export raster value into excel or ASCII

text: I'd like to get all pixel values of time series NDVI image to find regression coefficient to smooth the abnormal NDVI value through the time. So that firstly, I want to export values of raster into excel or ASCII file. And then using the exported values I want to run polynomial regression as a function of time in one of the statistic software to find the coefficients. Finally, I want to run regression equation with the coefficients on whole raster image again. Generally I'm going to find best fitting equation to the value of raster image. How can I do this task using one of the softwares (ARCGIS or ERDAS)? Thank you very much


title: Error with the Contour (Spatial Analyst) tool in ArcGIS

text: I am having an issue with the Contour (Spatial Analyst) tool in ArcGIS. When I manually run the tool in ArcToolbox, it will sometimes say that the tool has failed, giving the following error message: ExecuteError: ERROR 010157: Unable to open feature class F:\2014\data\083c15_contour.shp Failed to execute (Contour). The strange part is that when I go to the output location of the tool the shapefile is there, and I can open it in ArcMap without an issue. When I use Python to run the tool, I receive the same error message, and I can't figure out why. This is the loop I have that is meant to convert a list of DEM files into contour shapefiles:                   for dem_file in dem_files:             file_split = dem_file.split("_")             file_name = job_folder + "\\" + file_split[0] + "_contour.shp"             arcpy.sa.Contour(dem_file, file_name, contour_interval)      When I run the script, I can go to the output location and bring the shapefile into ArcMap, even though there was an error. So I don't think this is an issue with my code because it also happens when I run the tool manually. Does anyone know how I can work around this error?


title: GeoDjango order_by('distance') not working as expected?

text: Using GeoDjango, I am trying to sort a list of users based on the distance from a point. It does generate a list, but its not ordered by distance, though I specify that in my query. This is the query I am using.               origin = GEOSGeometry('SRID=4326;'+ 'POINT('+ str(buyer.buyer_current_location.x)+' '+str(buyer.buyer_current_location.y)+')')     close_sellers = Seller.objects.exclude(profile__user = user).filter(seller_current_location__distance_lte=(origin, D(mi=distance_mi)),profile__user__is_active = True).distance(origin, field_name='seller_current_location').order_by('distance')[:20]      and these are the models that I have               class Buyer(models.Model):         # This field is required.         profile = models.ForeignKey(UserProfile,unique=True)         # Other fields here         screen_name = models.CharField(max_length = 200,default='buyer')         buyer_current_location = models.PointField(null = True, srid=4326)         objects = models.GeoManager()          class Seller(models.Model):         # This field is required.         profile = models.ForeignKey(UserProfile,unique=True)         # Other fields here         seller_current_location = models.PointField(null = True, srid=4326)         objects = models.GeoManager()      Does anyone have any clue regarding what might be going wrong? # Edit The list seems sorted, but only messes up on for sellers in a certain part of town (its totally crazy!)


title: How to create correct vector grid in meter

text: I want to create a vector grid, with squares of 100x100m but have some serious trouble in solving this apparently easy task. I already read all the posts, but didnt find a good solution. Here is what I tried:   1. Using the research > vector grid tool: I always fail to have the right X and Y Values, it is never 100m, even if I have "m" as my default lenght.I need 100m but it says always something around 59,... if i measure it with the measurement tool. and also the vertical lenght seems to be slightly bigger than the horizontal lenght.   2. Using processing tools: create grid - I have the same problem with lenght.    3. Using the MMQGIS create Grid tool: Same problem as 1 and 2. how can I change the way i calculates distances? I think I have to look somewhere else, but dont know where. FYI: Running system is Win 8, 64 bit, QGIS Version 2.4, Thanks!


title: Growth of a PostGIS database for Social Media - Move to Cassandra?

text: I currently have a large and growing database of geo-social media. I'm using Postgres with postgis, and so far it's been quite good. My use case thus far has been to query the database for all posts within specific areas, such as every tweet within each post/zip-code in the US/Canada (or other boundaries -- the posts are not geo-coded to the post/zip code. selection is through intersection of coordinates). In the past this was relatively quick and easy to do, however, now that my data table has grown to over 350m posts, my queries take forever to complete. Each post code can take up to a day to retrieve the tweets from within its geometry. I am currently debating what to do next. I know of a few solutions, but I'm curious what you think of them and if I'm missing anything.   1. throw more memory at postgres, and reduce disk searching. I'm currently working on doing this, but I wonder if it is a good long term solution.   2. find a better way to partition data into more tables. Currently, all of my posts are in the same table, with each column indexed. Given that I need to search through all of the posts to find the ones in the given areas, is this the wrong way to be doing it?   3. create a hybrid of Cassandra-PostGIS. Bulk store my data into Cassandra, and then pull out roughly the data I need (say, one state/province at a time) to a staging data table in PostGIS, then query individual post/zip codes using PostGIS. Is this overkill/horribly inefficient? I have looked into Postgres-XL, but it has been kiboshed by others in this project for cost issues. All help is appreciated, Thanks!!


title: Join leads to null value outside of edit session

text: I am using Arcmap 10.2 on Windows 7 I joined a shapefile with an Excel file, and the attributes are shown as until I start an edit session. Then, the values are as they should be. They match up by verifying manually the attribute table and the Excel file. I clicked "Save edit" and then ended the edit session, at which point the attributes went back to . Why is it doing that? How can I fix it?


title: R SP library: identicalCRS and parameters in different order

text: The following returns FALSE:               library(sp)     points <- matrix(rnorm(30, 50, 10), ncol=2)     proj1 <- CRS("+proj=utm +zone=10 +datum=WGS84 +ellps=WGS84 +units=m +no_defs +towgs84=0,0,0")      proj2 <- CRS("+proj=utm +zone=10 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0")     sp1 <- SpatialPoints(points, proj1)     sp2 <- SpatialPoints(points, proj2)     identicalCRS(sp1, sp2)      However, the proj4strings have the same options, just in a different order. Are the projections truly not equivalent? That is, I'm thinking of replacing identicalCRS with something like               length(setdiff(unlist(strsplit(sp1@proj4string@projargs, ' ')),                     unlist(strsplit(sp2@proj4string@projargs, ' ')))) == 0


title: Google Earth export from QGIS

text: I'm using Quantum GIS. Is there a way to export a shapefile to Google Earth format (*.kmz) using QGIS? Thanks.


I create Spatial ETL ToolBox in ToolBox in ArcGIS Catalog. When I run tool from catalog I see log massages. But when I run this tool from python script log are not print to python console.               import arcpy     arcpy.ImportToolbox(r"c:\Toolbox123.tbx")     arcpy.SpatialETLTool(param1,param2,...)      How to print log from FME when run from python script?


title: Creating an Address Locator with a Database of Geo Information

text: I am trying to build an Address Locator for ArcGIS 10.1. I would like to use a Microsoft Server MySQL database where location information is stored using an alias, an address, and XY coordinates. What is the best way to create a geocoding service using this setup? An option I have considered is writing a script to dump the SQL database into a format that plays well with ArcGIS on a regular basis but wanted to make sure there wasn't a better way to do this.


title: How to use GDAL to convert Int16 data to Byte

text: I have some data with Int16 bands in HDF format. I want to convert this to Byte bands in GeoTIFF format. The data range for the Int16 is -2000 to 12000, so I want to map this to the Byte range 0 to 255. How do I do this with GDAL tools? (if it is possible) Otherwise how can I write a Python script to do this.


title: min and max values are not correctly calculated for DEM

text: I'm working with DEM in QGIS and need to know the min and max values of my data. In the Properties window, under the Metadata tab, the statistics section shows such values. However, I usually find some pixels with values higher than the max and lower than the min. Is there another way to obtain these values? (I've tried this with several DEMs) Cheers, Andrés


title: Display Selected Features Attributes in User Form

text: I'm trying to convert my ArcGIS 9.3 code to ArcGIS 10 via Visual Basic 2008 Express; however, I'm having a couple issues. The goal is for a user to select a feature in ArcMap, click on this custom button, and a user form appears displaying attribute information of the selected feature. The code worked great in 9.3, but not in 10. Currently, I can get the user form to open by itself...but if a user selects a feature and clicks the button it crashes ArcMap.               **User Form Code:**         Open User Form         Public Class TransLineEdit             Inherits ESRI.ArcGIS.Desktop.AddIns.Button             Dim pForm As New frmEdit                Protected Overrides Sub OnClick()                 Try                     pForm.Show()                 Catch ex As Exception             System.Windows.Forms.MessageBox.Show(ex.ToString, "OnClick")                 End Try             End Sub         End Class                   **User Form Code:**         Imports ESRI.ArcGIS.ArcMapUI         Imports ESRI.ArcGIS.Carto         Imports ESRI.ArcGIS.Geodatabase         Imports ESRI.ArcGIS.Geometry         Imports ESRI.ArcGIS.ArcMap         Imports ESRI.ArcGIS.DataManagementTools         Imports ESRI.ArcGIS.Desktop         Imports ESRI.ArcGIS.Editor         Imports ESRI.ArcGIS.Framework              Public Class frmLineSelectAttribute     'Declare variables that are used all thru the application     Public pMxDoc As IMxDocument     Public pMap As ESRI.ArcGIS.Carto.IMap     Public pEditor As IEditor     Public Application As IAppROT               Private Sub cmdCancel_Click(ByVal sender As System.Object, ByVal e As System.EventArgs) Handles cmdCancel.Click              'Kill form         Me.Close()              'Clear selected features         'Dim pActiveView As ESRI.ArcGIS.Carto.IActiveView         'pActiveView = pMap         'pMap.ClearSelection()         'pActiveView.Refresh()     End Sub     Shared Function AccessLayersData(ByVal pLayer As ESRI.ArcGIS.Carto.ILayer) As IFeatureClass         Dim pFeatureLayer As ESRI.ArcGIS.Carto.IFeatureLayer = pLayer         If Not pFeatureLayer Is Nothing Then             Return pFeatureLayer.FeatureClass         Else             Return Nothing         End If     End Function               Private Sub UserForm_Initialize()              'Get Improvement Projects FeatureLayer         Dim pFeatLyr As ESRI.ArcGIS.Carto.IFeatureLayer         pFeatLyr = GetLayerByTOC("TransportationImprovementProjects_Lines")              'Get selected feature (only one)         Dim pFeat As IFeature         pFeat = GetSelFeat(pFeatLyr)              'Populate form controls with feature's attribute values         txtPrjYear.Text = pFeat.Value(pFeat.Fields.FindField("YEAR"))         txtPrjType.Text = pFeat.Value(pFeat.Fields.FindField("PROJECTYPE"))         txtAgency.Text = pFeat.Value(pFeat.Fields.FindField("AGENCY"))         txtLocation.Text = pFeat.Value(pFeat.Fields.FindField("LOCATION"))         txtDescrip.Text = pFeat.Value(pFeat.Fields.FindField("DESCRIPTION"))         txtCost.Text = pFeat.Value(pFeat.Fields.FindField("COST"))         txtDate.Text = pFeat.Value(pFeat.Fields.FindField("DATEEDITED"))         txtUser.Text = pFeat.Value(pFeat.Fields.FindField("USEREDITED"))              'Dim pDoc As IMxDocument         'Dim pItem As ICommandItem         'pItem = Project.ThisDocument.CommandBars.Find(arcid.Query_ZoomToSelected)         'pItem.Execute()         'pDoc.ActiveView.Refresh()          End Sub     Public Function GetSelFeat(ByVal pFeatLyr As ESRI.ArcGIS.Carto.IFeatureLayer) As IFeature              'Gets selected feature from feature layer and returns feature         'Call this function only when you are sure that only one feature is selected              'Initialize the required variables         Call Initialize()              Dim pFeatSel As ESRI.ArcGIS.Carto.IFeatureSelection         pFeatSel = pFeatLyr              Dim pSelectionSet As ISelectionSet         pSelectionSet = pFeatSel.SelectionSet              pSelectionSet.Search(Nothing, True, pFeatSel)              'Return the selected feature         GetSelFeat = pFeatSel.NextFeature          End Function          Public Sub Initialize()              'Initalize global variables         pMxDoc = Application.Document         pMap = pMxDoc.FocusMap              'Instantiate the editor         pEditor = Application.FindExtensionByName("ESRI Object Editor")          End Sub          Public Function GetLayerByTOC(ByVal lyrName As String) As ESRI.ArcGIS.Carto.IFeatureLayer              'This function finds a feature layer based on its TOC Name         Call Initialize()              Dim lyrCntr As Integer         Dim pFeatLyr As ESRI.ArcGIS.Carto.IFeatureLayer              For lyrCntr = 0 To pMap.LayerCount - 1                  'Ensure that the layer is valid             If pMap.Layer(lyrCntr).Valid = True Then                      'Ensure that the layer is a feature layer                 If TypeOf pMap.Layer(lyrCntr) Is ESRI.ArcGIS.Carto.IFeatureLayer Then                     pFeatLyr = pMap.Layer(lyrCntr)                          If UCase(pMap.Layer(lyrCntr).Name) = UCase(lyrName) Then                         GetLayerByTOC = pFeatLyr                         Exit Function                     End If                 End If             End If         Next lyrCntr              'Layer not found. Show a message         MsgBox("Layer " & lyrName & " not found !", vbCritical, "Error")          End Function              End Class


title: Simple pixel analysis

text: I am very new to QGIS and I simply have two tif images which I need to do pixel analysis on so I can compare the change in vegetation over time. I have no idea how to go about this and if anyone has some advice for me I would be very grateful


title: Can you use active layer as input for dynamic text?

text: I'm wondering if you can lift the name of a layer in the ToC (within a group) and use it in dynamic text. Is there existing syntax to do this?   Perhaps ...               <dyn type="ToC" property="activeLayer"/>      Using ArcMap 10 sp3.


title: What projections should I use to make my own Globe?

text: Searching for an answer to this question, I found instructions posted by Gulf of Maine Research Institute showing how to create a globe. ![enter image description here](http://i.stack.imgur.com/hVBEb.gif) Using manual methods ... ![enter image description here](http://i.stack.imgur.com/oLCn2.gif) What approach would I take to create a globe using GIS? What projection should I use for each individual gore? If I wanted fewer seams near the poles, is there some other projection I could use? Could I do a series of projections to create the gores based on a soccer ball and stitch them together? ![enter image description here](http://i.stack.imgur.com/1M3HN.jpg) How would I determine the point of tangency for each pentagon and hexagon, along with their vertices in latitude/longitude? ![enter image description here](http://i.stack.imgur.com/oZxw8.png) Would some other non-soccer-ball isohedron be more suitable?


title: How to turn WKBtype integer into text for memory layer creation in QGIS

text: I'm creating a memory layer in QGis python and I need to match the geometry of the memory layer that I am building to the geometry of another layer. After some investigation it seems to me that, as per QGIS 2.2 Documentation and the QGIS 2.2 API Documentation, the URI defining the geometry type approximately matches the QGis::WkbType member enumeration text. I know I can easily query the WKBType of my input layer with               myWkbType=mylayer.dataProvider().geometryType()      However, this returns an integer whereas the constructor requires text. Is there a simple built-in pyQGIS way to turn these constants into text or do I need to build a hacky little reverse mapping function for it?


title: Using IBufferConstruction to create single and multiple buffers (VBA and ArcObjects)

text: I need a bit of help with using the iBufferContruction interface, my overall goal is that the user clicks on the map, a graphic is created and then buffered,at multiple distances and then stored in a featureclass So far, i can create a single buffer (i know if it created because the length that gets returned is correct), however when I store it in my FC the shape is empty. Here's the code:               Private Sub Buffer_MouseDown(ByVal button As Long, ByVal shift As Long, ByVal x As Long, ByVal y As Long)          'Get the location of the mouse and set it to the point     Set m_pPoint = m_pMxDoc.ActiveView.ScreenDisplay.DisplayTransformation.ToMapPoint(x, y)          'new marker element on the map     Set m_pElement = New MarkerElement          'set the element property     m_pElement.Geometry = m_pPoint          Set m_pGraphics = m_pMap          m_pGraphics.AddElement m_pElement, 0 'Grpahic for the point click.          Set m_pAC = m_pMxDoc.ActiveView          'refresh the map     m_pAC.PartialRefresh esriViewGraphics, Nothing, Nothing     '     Dim pPoly As IPolygon4     Set pPoly = New Polygon          'Dim pPoly As IPolygonElement     'Set pPoly = New PolygonElement          Dim pBC As IBufferConstruction     Set pBC = New BufferConstruction          'Dim pBCProp As IBufferConstructionProperties     'Set pBCProp = pBC     'pBCProp.ExplodeBuffers = True          Set pPoly = pBC.Buffer(m_pElement.Geometry, 200)     Debug.Print pPoly.Length ' RETURNS 1256.63706143592          'pBC.ConstructBuffers     Dim pFeat As IFeature     Dim pfeatureCursor As IFeatureCursor     Dim pfeatureBuffer As IFeatureBuffer          Set pfeatureBuffer = m_pBufferFClass.CreateFeatureBuffer     Set pfeatureCursor = m_pBufferFClass.Insert(True)     Set pFeat = pfeatureBuffer          pfeatureCursor.InsertFeature pfeatureBuffer     pfeatureCursor.Flush     Set pFeat.Shape = pPoly          m_pAC.PartialRefresh esriViewGeography, Nothing, Nothing      Where am I going wrong? Also has anyone got any code examples with creating multiple buffers at various distances?


title: Parts needed to run Geoserver/PostGIS from a program like Microsoft Access

text: I have a fair bit of experience in GIS and Python, but am brand new to developing on the web. I have rented some space on a server and have a PostgreSQL/PostGIS database as well as Geoserver on it. The other side of the coin is a Microsoft Access database (which will later be a web based application, but that's later) that has a bunch of data related to the Spatial data (polygon layers) up on the server. What components do I need to develop so that the Access database can request specific data from Geoserver and have it display on a webpage (that is seen by the user through the Access program). I also want the Access program to be able to request a KML layer that will display in Google Earth. For the first option I think that I need to use the OpenLayers GetMap examples and build a map from that. The map will need to be able to get all of a particular layer, as well as one polygon from the same layer displayed differently so that it stands out. I remember that there is a way to integrate Google Maps as a base layer. Are there other/better options? For the second option (use the Access database to query Geoserver and display the result in Google Earth) I'm at a loss. It seems like it can be done, but how? I just need to be pointed in the right direction as most of these technologies are new to me and I'm a bit lost. Thanks a lot!


title: How to read files one by one with FME

text: In one folder I can have one or more xml files. I would like (need) to use them one by one as input dataset (reader). Actually each xml file will be used within according fme workbench regardless some values within xml file. This question is linked with my old question Variable XML dataset as reader in FME wher I got answer for one file in folder. If I have two xml files, for instance, fme read them both at the "same time" and use values from both files for one fme wb.


title: how does arcmap calculate scalebar inside a WGS84 layout?

text: ArcMap lets you put a scalebar inside a layout referncing a WGS84 data frame. How is the scale calculated? Is it a "mean" scale? And, in case it's right, what "mean" scale could mean?


title: Is getting wrong attribute values from a feature in qgis

text: I make two polylines features(in my example both consists of two points) with the add feature tool. Then I (with code) add some attribute values on the first feature, then changes the geometry of the second feature before I add values to the second values. When I use layer.featureAtId() to get the second feature and then reads it, I get the values from the first feature. I don't get this problem if I don't change the geometry or if I commit the two new polylines before doing the mentioned operations. In the attribute table it is correct. If I use `layer.nextFeature(feat)` to iterate through all the features I get the right value.               # creates two features with the add feature tool     layer = qgis.utils.iface.activeLayer()          #select the first feature with the select tool      featId1 = layer.selectedFeaturesIds()[0]          # add values  to the first feature     for i in range(5):         layer.changeAttributeValue(featId1, i, (i*10))               #select the second feature with the select tool      featId2 = layer.selectedFeaturesIds()[0]          # changes the geometry of the second feature     newGeom = QgsGeometry.fromPolyline([QgsPoint(1,1), QgsPoint(2,2)])     layer.changeGeometry(featId2, newGeom)     qgis.utils.iface.mapCanvas().refresh ()          # add values  to the second feature     for i in range(5):         layer.changeAttributeValue(featId2, i, (i*10 + 1000))               # get the second feature by using featureAtId...         feat = QgsFeature()     layer.featureAtId(featId2, feat, True, True)          # print values from the second feature     for i in range(5):         print feat[i].toString()          # for some reason the attributeValues to the first feature get printed out.      # in the attribute table it looks correct.      # feat.geometry().asPolyline() seems to return the right geometry     # if I save edits after creating the two new features it works as it should         I used QGIS version 1.8.0-Lisboa and QGIS version 1.9.0-Master


title: How can I change mouse cursor while it hovers over the map in Oracle Mapviewer

text: I'm using Oracle Mapviewer. I want to switch mouse pointer from default to crosshair when user enters "Interactive mode". I'm developing in GWT and tried do do something like this:               public static void setCursorCrosshair() {         removeCursorStyles();         getMapWidget().addStyleName("cursorCrosshair");     }      But it changes mouse pointer only in few parts of map, it looks like it omits map tiles and change cursor to crosshair on FOI's that I've added just before. How can I make it to display crosshair over the whole map?


title: Deleting columns from joined attribute table in QGIS

text: After joining the attribute tables from two layers, I have a lot of duplicate columns on my table, making it hard to read. However, even when in "edit" mode I am unable to select and delete columns. Is there a way to do this?


title: how to see a base map created by map builder in database?

text: I am new in oracle, I tried to create a style, then different themes and then a base map. I saved all my work, but I can't see the base map in the database. Did I miss something? The interesting part is that I have a spatial table, which has a geometry field, but I can't see that in Map Builder. Can anyone help me? I'm really stuck here...


title: How to extract elevation information from contour lines?

text: I'm working with OS windows 7 and using ArcGIS 10.2 and I have a point shape, which I must extract elevation for every point in it but the only useful data is a contour line shape. Not sure if it's useful but I need your help once again to lead me on how do I extract the data I need from it.


title: Mapnik + OpenStreetMap + Shapefile with Raster Images + Mod Tile + PostGIS

text: I am new to GIS. I have a Shapefile (.shp, .shx and .pbf). The .pbf file contains path to PNG files located on my machine only. I was able to save the Shp data to PostGis via Shape2Gis tool. Now that a tables are created with PostGIS. How can I use this table to show data with Mapnik? I am using Mapnik + OpenStreetMap and Mod_tile.I have been previously using osm2pgsql and the maps are being rendered. I have a shapefile, shx file and pbf file. Not sure how to use it to render maps..


title: How do I change polygon opacity using a esri WFS in leaflet?

text: Note...Solved and here for anyone having this problem...damn simple problem, I had not used a higher case character for fillcolor and fillopacity....they should be fillColor and fillOpacity , it works after changing that. mmm I seem to be able to change the opacity of the polygon outline but not the fill itself. Can find what is going wrong with my code. I try to use the fillcolor and fillopacity but nothing changes on my map?               var LMDB1 = new L.FeatureGroup();              L.esri.featureLayer("https://services.arcgis.com/HIDDEN/arcgis/rest/services/HIDDEN/FeatureServer/1", {                  style: function (feature) {             if(feature.properties.TAG === 'Lowest Health'){                 return { color: "#AD0909", weight: 1, opacity: .1, fillcolor: "#000000", fillopacity: 1};                 } else if(feature.properties.TAG === 'Low to Average Health'){                 return { color: "#FF8000", weight: 1, opacity: .1 };                 } else if(feature.properties.TAG === 'Average Health'){                 return { color: "#FFFB0D", weight: 1, opacity: 1 };                 } else if(feature.properties.TAG === 'Average to High Health'){                 return { color: "#C4FFD3", weight: 1, opacity: 1 };                 } else if(feature.properties.TAG === 'Highest Health'){                 return { color: "#0E8F2E", weight: 1, opacity: 1 };                 } else {                 return { color: "#FFFFFF", weight: 1, opacity: 1 };                 }             },             onEachFeature: function (feature, layer) {                 layer.bindPopup(L.Util.template(poly_templ, feature.properties));}                  }).addTo(LMDB1);


title: ESRI ArcMap Add-in with app.settings not recognizing app.config changes

text: I have an ESRI Add-in developed for ArcMap which requires a configuration file. After spending a while trying to read config values from a single App.Config file (and always getting null) I believe that the Add-in cannot read values from here as it is a class library, and looks for the calling application (ArcMap)'s config file when I ask for a key's value (hence the null). To get around this I used an App.Settings file, which the application can read fine. Creating this also introduces an App.Config file into the environment and Visual Studio seems to keep the two files in-sync during development. Now that the Add-in is being deployed I need to be able to change configuration values (e.g. log file location). I have tried opening / extracting the .esriaddin file and updating the App.Config file in there but the Add-in retains the same configuration values it had when compiled. I know that the new App.Config values are being persisted in the .esriaddin file because I can view them again after closing the archive. Does anyone know a reliable way to configure an Add-in and permit this configuration to be updateable once deployed? Any suggestions very welcome as it seems ridiculous that I should need a custom config file for this. App.Settings values are at the application level, and currently both App.Settings and App.Config have build action: none / do not copy.


title: ArcGIS 10 - How can you open geodatabase files using a Browse Dialog Command (VB.NET/C#)?

text: I am currently using the following method to open raster files from a form:               Public Function BrowseRaster(ByVal sTitle As String, ByRef sFolder As String, ByRef sName As String) As Boolean           Dim pGxDialog As IGxDialog = New GxDialog         Dim pRasterFilter As IGxObjectFilter = New GxFilterRasterDatasets()              Dim pFilterCol As IGxObjectFilterCollection         Dim pEnumGx As IEnumGxObject         Dim pGxObject As IGxObject              sFolder = ""         sName = ""              pFilterCol = pGxDialog         pFilterCol.AddFilter(pRasterFilter, True)         'pGxDialog.StartingLocation = ""         pGxDialog.RememberLocation = True         pGxDialog.AllowMultiSelect = False         pGxDialog.Title = sTitle                       If pGxDialog.DoModalOpen(0, pEnumGx) Then                 pGxObject = pEnumGx.Next                 Dim sFile As New FileInfo(pGxObject.FullName)                 'sName = pGxObject.BaseName                 sName = pGxObject.Name                 sFolder = sFile.Directory.FullName                 Return True             End If          End Function      How can I modify this function to also open files from inside of a file geodatabase?


title: Deleted features not synchronizing in one-way (ArcSDE to File Geodatabase) geodatabase replication?

text: I am testing three one way (ArcSDE 10 -> File Geodatabase) replicas. Each replica is derived from a layer created by selection by location from a citywide road dataset, so that each of the three layers is distinct from one another. When I add features (all line layers) to the ArcSDE parent replica, all the features appear in the children after synchronization, as I would expect. The problem occurs when I delete any new test feature that is beyond the extent of the original dataset. If I first add a new feature beyond the dataset extent to the parent, this feature will be added to the child, but then when I delete the same feature in the parent, this feature will remain in the child after synchronization. I do not have this issue when testing replicas that use the entire (non- selected) dataset. I have tried some different combinations of extent settings in the replica creation and in the environment settings but none have worked yet. I suspect this is related to the fact that these test features do not intersect the original extent...but I can add features outside of the extent, and I do not have this problem when testing a replica using the 'full' dataset. Any ideas?


title: Offline Map Service with Mapnik and Leaflet

text: I'm relatively new to GIS, so I hope my asking doesn't appear to be ignorant in any way. I currently work on an offline map service solution for usage in areas without internet connection. For that purpose I want to download the data from a provider (i.e. geofabrik.de), load it into a PostGIS enabled database by using osm2pgsql, process the data with Mapnik and store the rendered tiles on the local machine. The map data doesn't need to be up-to-date but should be able to get updated once or twice a year. From there, Leaflet will be used to display the rendered tiles and also satellite tiles if possible, as a slippy map. The map user should also be able to edit data on the map and to save them as vector layer to the database. ![enter image description here](http://i.stack.imgur.com/OX4fW.png) I already did something similar by using Geoserver, GeoWebCache and OpenLayers but I'd like to have at least two approaches for comparison. My questions are:   * Where is the difference in using Mapnik or Geoserver? This may seem to be a redundant question but I couldn't answer it - so far.   * Where do I get the satellite tiles? I read something about OSM providing servers for them but only give out higher capacity downloads to chosen projects.   * Are there any performance issues I should look out for? The provided systems are old Panasonic Toughbooks with limited CPU and RAM. I already tried the Geoserver way on them. It worked, since there is no real-time rendering. But maybe I'm overlooking something.   * Any flaws with the approach?  Thanks for investing time in this!


title: Projecting gps data to get length of resulting line shape in ArcMap 10.0

text: I have gps data in WGS84 geographic coordinates. I will use the gps data to build a line shape in ArcMap 10.0. The data falls in UTM 15N. To get most accurate length of line, should I project the point data to the UTM and then build the line in the UTM. Or should I build the line in WGS84 then project the line data to the UTM. Or does it not matter? What are the effects? If UTM data doesn't best preserve distance then I will add that UTM must be used.


title: Does arcpy.AggregatePoints_cartography create multiple polygons?

text: Does this function allow you to put in an array of points, which may, or may not, be representative of many polygons, and create many polygons from it, or just the one, linked?


title: Get mouse position in GeoExt

text: The title says it all. I'm using GeoExt as part of OpenGeo Suite. I would like to display the current mouse coordinates. Update: I was getting confused as to what parts of the code are GeoExt and which parts are OpenLayers. Here's some code. mapPanel is a GeoExt object. mapPanel.map is an OpenLayers object.               var mapPanel = new GeoExt.MapPanel({          // etc.     });          mapPanel.map.addControl(new OpenLayers.Control.MousePosition());


title: How to connect PostgreSQL from osm2pgrouting

text: I am using OSGeo Live. It is an Ubuntu system. I installed it in my PC. I was following one pgRouting workshop. http://workshop.pgrouting.org/index.html In one step, I needed to connect osm2pgrouting with postgis to import data. But, the commands are not working. It is showing following connection fail error. I did not change anything. All the passwords are as default.               $ osm2pgrouting -file "data/sampledata.osm" \     >                           -conf "/usr/share/osm2pgrouting/mapconfig.xml" \     >                           -dbname routing \     >                           -user postgres \     >                           -clean     host=127.0.0.1 user=postgres dbname=routing port=5432     connection failed     $      I tried to fix it. I edited **/etc/postgresql/8.4/main/pg_hba.conf** , and change the postgres connection method to **trust** from **indent**. But, it does not work. The result is same, as before.               # Database administrative login by UNIX sockets     # indent     local   all         postgres                          trust      Thanks, in advance.


title: Using Google Earth to query database

text: I've been reading about using Google Earth to query data (this article) and I'm interested in learning how to make one of these. I just have a couple of questions: 1. Is this possible with the free version of Google Earth? 2. Can Google Earth connect to any database application (i.e. Access/SQL Serve)? 3. Is it a good query tool (in comparison to querying layers in ArcGIS or a relational database)? 4. How long does it take to make one of these tools? Thanks for your help everyone!


title: Overlaying WMS from Mapserver in OpenLayers with epsg:900913 from Google/Bing/OpenStreetMap

text: ADDED: Another curious effect is that no matter to what scale I zoom on the baselayer google/bing map, the WMS overlay is always the same-sized shot of MN. That is, as the baselayer map is scaled up or down, the WMS overlay map never scales -- it always shows the same bounds. I am trying to convert an Openlayers page from running overlay maps with OpenLayers.Layer.MapServer to one running WMS with OpenLayers.Layer.WMS. The information I am providing through these interfaces is overlaying on a Google/Bing/Openlayers baselayer. The Mapserver implementaion works fine, but the WMS implementation seems to generate an overlay from my data that is ~ 2x too large for the overlay with the baselayer maps. ** Please excuse the use of links here and for pastebin... I wanted to embed images and make links directly to pastebin, but being a new user, I was given restrictions of no images and only two links. Here is a picture of the working overlay: ![enter image description here](http://i.stack.imgur.com/8bCue.png) Here is a picture with the problematic WMS overlay: ![Problem](http://i.stack.imgur.com/n68ta.png) The mapfiles for mapserver are mostly identical, as I used the working Mapserver mapfile to create the mapfile for WMS, but with some necessary mods to the WMS mapfile to make it WMS enabled. Here is the working Mapserver map file: http://pastebin.com/VQBW5YSx And here is the WMS mapfile: http://pastebin.com/swMagqma Any help is appreciated. Thanks, Matt


title: Where can I download a xyz bathymetry file of the Northwest Atlantic?

text: I'm looking to download the bathymetry of the Northwest Atlantic ocean as a txt file or any other file format that I can easily upload and use in R. I once used the Global topography - Scripps Institution of Oceanography but it looks like I'm running into an error message this time...? I have been looking at other websites, but I'm having a hard time to locate a new site where I could download a xyz file of my area of interest. I would appreciate any recommendations. Thank you very much in advance,