title: Get mouse position in OpenGeo

text: I'm using GeoExt as part of OpenGeo Suite using code generated by the Open Geo suite SDK. I would like to display the current mouse coordinate and don’t really know how to do it. As far as I know GeoExt has no plugin for this and don’t know how to add the OpenLayers.Control.MousePosition() to mi code. I found a simillar question here but in my case I dont have a GeoExt.MapPanel, instead I have a gxp.Viewer. Here is my code               /**      * Add all your dependencies here.      *      * @require RowExpander.js      * @require widgets/Viewer.js      * @require plugins/LayerManager.js      * @require plugins/AddLayers.js      * @require plugins/RemoveLayer.js      * @require plugins/FeatureManager.js      * @require plugins/FeatureGrid.js      * @require plugins/OLSource.js      * @require plugins/OSMSource.js      * @require plugins/WMSCSource.js      * @require plugins/WMSGetFeatureInfo.js      * @require plugins/Zoom.js      * @require plugins/ZoomToExtent.js      * @require plugins/ZoomToLayerExtent.js      * @require plugins/Legend.js      * @require plugins/GoogleGeocoder.js      * @require plugins/GoogleSource.js      * @require widgets/CrumbPanel.js      * @require plugins/QueryForm.js      * @require plugins/ZoomToSelectedFeatures.js      * @require plugins/LoadingIndicator.js           */          var app = new gxp.Viewer({         portalConfig: {             layout: "border",             region: "center",                  // by configuring items here, we don't need to configure portalItems             // and save a wrapping container             items: [{                 id: "centerpanel",                 xtype: "panel",                 layout: "fit",                 region: "center",                 border: false,                 items: ["mymap"]             }, {                     // container for the layer manager etc.                     id: "westpanel",                     region: "west",                     xtype: "gxp_crumbpanel",                     collapsible: true,                     //collapseMode: "mini",                     //hideCollapseTool: true,                     split: true,                     border: true,                     width: 200                 }, {                     // container for the FeatureGrid                     id: "south",                     region: "south",                     xtype: "gxp_crumbpanel",                     collapsible: true,                     collapseMode: "mini",             collapsed: true,                     hideCollapseTool: true,                     split: true,                     border: true,                     height:200             }],             bbar: {id: "mybbar"}         },              // configuration of all tool plugins for this application         tools: [{         ptype: "gxp_layermanager",         outputConfig: {             id: "tree",             title: "Capas",             tbar: [] // we will add buttons to "tree.bbar" later             },             outputTarget: "westpanel"         }, {             ptype: "gxp_addlayers",             actionTarget: "tree.tbar"         }, {             ptype: "gxp_removelayer",             actionTarget: ["tree.tbar", "tree.contextMenu"]         }, {             ptype: "gxp_zoomtolayerextent",         actionTarget: ["tree.tbar", "tree.contextMenu"]         }, {             ptype: "gxp_zoomtoextent",             actionTarget: "map.tbar"         }, {             ptype: "gxp_legend",             actionTarget: "tree.tbar"         }, {             ptype: "gxp_zoom",             actionTarget: "map.tbar"         },{             ptype: "gxp_wmsgetfeatureinfo",             outputConfig: {                     width: 400                     }         }, {             ptype: "gxp_googlegeocoder",             outputTarget: "map.tbar",             outputConfig: {                     emptyText: "Buscar(Google)..."                 }         }, {             ptype: "gxp_featuremanager",             id: "featuremanager",             paging:false,             autoSetLayer: true,         //autoLoadFeatures:true         }, {         ptype: "gxp_featuregrid",            featureManager: "featuremanager",            outputConfig: {                id: "featuregrid"             },                outputTarget: "south"         }, {         ptype: "gxp_queryform",            featureManager: "featuremanager",            outputConfig: {             title: "Busqueda",             width: 400             },            actionTarget: ["tree.contextMenu", "map.tbar"],//"featuregrid.bbar"            appendActions: false,            autoExpand: "south",         }, {         ptype:"gxp_zoomtoselectedfeatures",            featureManager: "featuremanager",            outputTarget: "featuregrid.bbar"         }, {         ptype: "gxp_loadingindicator",         loadingMapMessage:"Cargando mapa..."         }],              // layer sources         defaultSourceType: "gxp_wmscsource",         sources: {             local: {                     ptype: "gxp_wmscsource",                     url: "/geoserver/wms",                     version: "1.1.1"                 },                  osm: {                  ptype: "gxp_osmsource"             },                 google: {                 ptype: "gxp_googlesource"             }         },              // map and layers         map: {             id: "mymap", // id needed to reference map in portalConfig above             title: "Mapa",             projection: "EPSG:900913",             center: [-10764594.758211, 4523072.3184791],             zoom: 4,             //controls: [pos],             layers: [{                 source: "google",                 name: "SATELLITE",                group: "background"         }, {                 source: "google",                 name: "ROADMAP",                group: "background"         }, {                 source: "osm",                 name: "mapnik",                 group: "background",             selected: true             }, {                 source: "local",                 name: "CONAGUA:estados",                 selected: true             }, {                 source: "local",                 name: "CONAGUA:Sitios",                 selected: true             }],             items: [{                 xtype: "gx_zoomslider",                 vertical: true,                 height: 100             }]         }          });


title: Repeat the same edit on several layers

text: For a project we have a master polygon layer which is the basis for a few "produced" layers. These produced layers may contain some of the original features and fields, and also additional features/fields. Occasionally there is a need to modify a few polygons in the master layer (geometry/attributes), and the same edits need to be performed for the "produced" layers as well to maintain consistency. Is there a way to automate this by a customized tool in ArcGIS, so I don't need to repeat the edits several times manually? It seems to me that if all modifications are done via geoprocessing, then it's pretty straightforward. Is it possible if the modifications are done as manual edits in an editing session? Will this entail editor extension? Edit*: I'm looking for solutions with ArcObjects or Python. I have ArcGIS for Desktop 10.0.


title: How do I create a polygon layer from multiple touching line segments in QGIS?

text: I'm trying to create a base map for a Maine Coastal area, and wanted to create a polygon layer of the land cover from the coastline in QGIS. Whether I created the coastline vector layer from the USGS national map or downloaded it from NOAA, both options leave me with a multipart line. Ignoring islands the USGS data has ~100 _almost_ connected line segments and the NOAA layer (below) is about 800 connected line segments. I had thought I could disolve these and run linetopoly, but that seems to have some more processing needed that I don't understand, does anyone see what I should be doing to create or get my polygon layer? The raw data looks like this, maybe I could make the coast into one segment?               856      -69.340984,44.013133     ...      -69.345531,44.014013     END     857      -69.356459,44.008586     ...      -69.356459,44.008586     END     858      -69.340544,44.012986      -69.341057,44.013719     etc...      ![failed linestopoly](http://i.stack.imgur.com/SeQGa.png)


title: overpass-api: data export into csv-formate

text: **update;** perhaps we can do this like so: to use Perl for text-mangling - so we can use the XML::Simple module. here's an example of a little script to parse your XML: see the Code:               #!/usr/bin/perl     use strict;     use warnings;     use XML::Simple;     use Data::Dumper;          my $xmlfile = shift || die "Usage: $0 <XML_FILE>\n";          my $ref;     eval {       $ref = XMLin($xmlfile,         ForceArray    => 0,         KeyAttr       => [ ],         SuppressEmpty => '',       ) or die "Can't read XML from $xmlfile: $!\n";     };     die $@ if($@);     print Dumper $ref;      **Explantion:** iterating thru the array/hash it creates the file and helps carving up the data into comma separated lines of data than can be redirected to file. see the basic question: i run the following code in opverpass-api - see here http://overpass-turbo.eu/ i have the options to export of the data to the following formats to GeoJSON to GPX to KML and to get the data from Overpass API loat them to JOSM laden (only for requests, that give back valid OSM-XML with Metadata) GeoJSON to save it as gist note - i did not install the overpass-api on my opensuse 13.1 yet. but i am willing to do so. as for now - running the above mentioned code in the oerpass-api - here. how to treat it to get it exported as csv-formated hope i was able to provide all the necessary things for a clear and concise question. See the output that needs to be transformed:               <?xml version="1.0" encoding="UTF-8"?>     <osm version="0.6" generator="Overpass API">     <note>The data included in this document is from www.openstreetmap.org. The data is made available under ODbL.</note>     <meta osm_base="2014-04-27T13:49:02Z"/>            <node id="297489767" lat="49.4085014" lon="8.6941465">         <tag k="addr:city" v="Heidelberg"/>         <tag k="addr:housenumber" v="23"/>         <tag k="addr:postcode" v="69115"/>         <tag k="addr:street" v="Sofienstraße"/>         <tag k="name" v="ARLT"/>         <tag k="phone" v="+49 6221 20229"/>         <tag k="shop" v="computer"/>         <tag k="source" v="survey"/>         <tag k="website" v="http://www.arlt.com"/>         <tag k="wheelchair" v="yes"/>       </node>       <node id="305144906" lat="49.4060012" lon="8.6929652">         <tag k="addr:city" v="Heidelberg"/>         <tag k="addr:country" v="DE"/>         <tag k="addr:housenumber" v="13-15"/>         <tag k="addr:postcode" v="69115"/>         <tag k="addr:state" v="Baden-Württemberg"/>         <tag k="addr:street" v="Rohrbacher Straße"/>         <tag k="name" v="Heidel-bike"/>         <tag k="opening_hours" v="Tu-Fr 10:00-18:30; Sa 10:00-14:00"/>         <tag k="shop" v="bicycle"/>         <tag k="website" v="http://www.heidelbike.de/"/>         <tag k="wheelchair" v="yes"/>       </node>       <node id="305963167" lat="49.4139877" lon="8.6924247">         <tag k="addr:city" v="Heidelberg"/>         <tag k="addr:country" v="DE"/>         <tag k="addr:housenumber" v="4"/>         <tag k="addr:postcode" v="69120"/>         <tag k="addr:street" v="Brückenstraße"/>         <tag k="name" v="Buchhandlung Schmitt & Hahn"/>         <tag k="shop" v="books"/>         <tag k="wheelchair" v="no"/>


title: Is it possible to preview a raster catalogs raster field in a Windows form?

text: **What I want to achieve:** The attribute table for an ArcGIS Raster Catalog contains a raster field. When a raster field cell is clicked with the mouse, a preview window displaying the raster associated with the row is displayed. I want to achieve the same or equivalent functionality, but within a windows form. **What I've tried:** I have managed to display the full attribute table of a RasterCatalog layer within a .Net DataGrid using this code  'Bind a geodatabase table to a .NET control'. However, whilst the raster column is populated with text 'raster', it is not possible to get a preview of the raster. One solution could be to add the image for each row to the datagrid as detailed here, but I'd prefer if I could keep the ESRI raster field functionality if possible. Any ideas for solutions / alternative approaches gladly welcomed. Cheers.


title: What is an "XYZM" measurement model?

text: SpatiaLite supports a remarkable number of different dimension models. Interestingly, some of their values support a "measure" coordinate (XYM, XYZM). How would this measure-coordinate be used, and what is the correct term for searching to figure out how it is used?


title: KML files to IPad securely? Altitude and shape problems in Google Earth

text: Does anyone know any options for maintaining privacy of the files and getting them onto iPad using Google Earth? Using Google account to link to Google Earth, privacy settings in Google my maps are Public and Unlisted (still public if you know the URL), no security of files here. Using MapBox, this requires use of TileMill and the privacy settings are the same public and unlisted (which are public if you guess the URL), support is very limited for KML/KMZ this way in any case. Using DropBox a cloud based option as is Sugar Sync, this is password protected and files open in Google Earth, how ever there is a problem with polygon/shape data when clicked on does not display the info bubble. Does anyone know how I would get them onto the iPad securely? And does anyone know whether there is a work around for the shape/polygon data displaying in google earth on iOS? what would happen if I changed the Altitude mode??


title: How do I load OSM vector data in QGIS 2.x?

text: In QGIS 1.8 I used the openstreetmap plugin to request data from the server, or open a .osm XML file, which produced 3 layers (polygons, polylines and points) in QGIS. This was great (give or take the 64bit bug which highly affects the freshly mapped area I work on), but how do I do that in QGIS 2.x? I can't find any updated documentation about it. In the vector menu, I've tried importing from osm servers, I get a .osm file, but then I'm stuck. There's a menu option for loading an XML file, but it produces a spatialite db, which apparently only has non-geometry tables when I try to add it as a spatialite layer. If I open the attribute table, I can see node ids, ways and so on, but I can't figure out how to actually make use of them. I've also tried adding a vector layer from a .osm XML file. I get a gpsbabel window that shows nothing, then finally an invalid data source message in QGIS. What am I doing wrong? Any pointers to updated documentation? thanks! PS: I know there are a few questions around this subject, but none of them actually seems to answer the question.


title: How to make a new raster using con() script in Python?

text: I'm new to Python and have been trying to code a raster calculator style script. Basically, I want to convert a raster that has values 0 to -10,000 into a raster that is populated by 1s between certain values. The reason I want to do it in Python is because I eventually want to automate this process so that I can pick out different values (ie 0 to -10, -5 to -25 etc) and end up with a raster of 1s for each. My aim then is to convert these to polygons so that I can then do a number of vector clip operations across about 800 shapefiles. I've written a script which produces a new raster at the end but instead of this grid being full of 1s it has values 0-256, so somewhere in my code there's a bug. I'd be grateful if anyone can spot it! If anyone has a way to speed up the code or a better idea for how I can do all this also then that would also be interesting to know. Here's my code:               import arcpy     from arcpy import env     from arcpy.sa import *          arcpy.env.overwriteOutput = 1 #this also doesn't work by the way, I have also tried                 #TRUE, no errors come up but it still doesn't let me overwrite the output file          #Set environment settings     env.workspace = "C:/folderlocation where inRaster is stored"          # Set local variables     inRaster = Raster("inputraster")          # Check out Spatial Analyst extension license     arcpy.CheckOutExtension("Spatial")          # depths wanted from inputraster raster     var = [-2500,-3000] #min depth, max depth          #Execute and save Con     OutRaster = Con((inRaster <= var[0])&(inRaster <= var[1]),1)          OutRaster.save("C:/drivelocation/outputfilename")      If anyone has any ideas they'd be gratefully received! Many thanks.


title: How to reproject all geometries in a PostGIS table?

text: I created a spatial table with SRID:4326. Now I want to change total projection to SRID:32644 into a new table. The old table should remain unchanged.


title: Is there a climatic geo portal for every point in the world?

text: I'm looking for a geo portal or climate database that provides me with the daily climatic information like temperature, precipitation and humidity. Forecast data is not of my interest. I need information for every point with spatial resolution of not more that 3km. Is there a service to provide me with this information? The most important thing that I'm looking for is the continuity of the database. The data in discrete stations is not enough for what I'm doing. All I want is that when I query a point with latitude and longitude, the database or the geoportal gives me the temperature and rainfall rate of that point (no matter how accurate it is)


title: OpenLayers - TileStache Overlay: tiles loaded but not visible

text: I have been trying this since yesterday. If I load only the tilestache layer, it is visible. If I load both bing and tilestache layer, the tilestache layer is not visible. When I look at the firebug console, the tiles are being loaded.               _map = new OpenLayers.Map({div: 'map'});     bing = new OpenLayers.Layer.Bing({         key: 'AqR9gcXm14-PKVQsWrL8H9A_RnusHZCzQKX_tbJegzU1EOZGafcjzCYbcpkwukQh',         type: 'Aerial'     });          tilestache = new OpenLayers.Layer.XYZ(             'hul',             'http://www.hyderabadurbanlab.com/tiles/tiles.py/hul2/${z}/${x}/${y}.png',             {               tileSize: new OpenLayers.Size(256, 256),               tileOrigin: get_point(17.4099, 78.4761),               sphericalMercator: true,               isBaseLayer: false,               transparent: true,               opacity: 0.6,               transitionEffect: 'resizing',               buffer: 1,               wrapDateLine: true             }         );     _map.addLayers([bing, tilestache]);      Attached is the screen shot of how it looks.![Bing and Tilestache](http://i.stack.imgur.com/6ywvt.jpg) When I load just the tilestache layer, it looks like this:![Just tilestache](http://i.stack.imgur.com/4bOsJ.png) I want the Bing maps as the base layer and the tilestache layer on top of it as an overlay. Thank you.


title: How to add KML data but from variable - not from url?

text: I have a working map that retrieves KML from a url, but now I need to get the KML data from a local variable instead. Is there support for this, and if so what is the syntax. Searched the api and various sites but no luck finding anything about this. Thanks! [Addendum]               map = new OpenLayers.Map({             div: "map",             layers: [                 new OpenLayers.Layer.WMS(                     "WMS", "http://vmap0.tiles.osgeo.org/wms/vmap0",                     {layers: "basic"}                 ),                 new OpenLayers.Layer.Vector("KML", {                     strategies: [new OpenLayers.Strategy.Fixed()],                     protocol: new OpenLayers.Protocol.HTTP({                         url: "kml_large.kml",                         format: new OpenLayers.Format.KML({                             extractStyles: true,                              extractAttributes: true,                             maxDepth: 2                         })                     })                 })             ],             center: new OpenLayers.LonLat(-81, 28),             zoom: 7         });


title: How to increase the display performance of Bing basemaps in ArcMap?

text: I'd like to know if there are ways (within ArcGIS itself) to increase the display performance of basemaps within ArcMap. The internet connection was recently upgraded, so it is much faster than before, and by that I mean instead of taking 1 minute to draw the Bing Imagery basemap, it now takes about 30 seconds. I normally pan and zoom with the basemap off, and once I think I'm in the location I want, I turn the basemap on and take my hands away from the mouse (trying to do anything else guarantees me a spinning circle, and most likely a crash). From this post on the ArcGIS forums, I've changed the location of my display cache from My Documents to my C drive, and I could see an improvement. I then changed the location to our server, and it now takes about 5 seconds to draw when I pan and zoom. While I am quite happy with this, does the fact that I am using 10.1 make a difference? I changed the cache locations for my colleagues who are using 10.0, and their performance did not improve at all. I don't know what else to do, as I can't do anything about the antivirus/firewall part (that's under IT's control). It's a bit frustrating, particularly when the basemaps in ESRI Maps For Office and ArcGIS for AutoCAD draw almost instantly, and I can zoom in without worrying my computer will implode. My current settings in ArcMap options are hardware acceleration enabled, and cache location changed. I've looked at the help on managing your local display cache, but my settings appear to be optimised already, unless I am missing something.


title: Touristic poi-data (api or dump) with popularity measure

text: I'm looking for freely available touristic (hotel, restaurants, bars, attractions, etc.) poi-data that have some measure of popularity. With popularity I mean, explicit popularity like on Tripadvisor (i.e: #10 out of 1500 attractions in Paris), but no api/ dump is available for TA. Or otherwise more soft-measures like the nr of checkins (Foursquare for example) Anyone knows of any data-set (dump or api) which makes this type of info available? First use would be to display markers on a map, where only the most popular are shown. (as showing all 5000 available bars in say Paris is not very practical) FYI, see list below of some sources I'm checking out   * Factual Global PLaces (dataset from former SimpleGeo): likely no popularity, not tried.    * openstreetmap (not tried yet, perhaps my next check)    * foursquare (has an api with implicit popularity (nr of checkins)) but hard/impossible to query the api to get complete coverage (besides it's against their ToS   * Facebook Graph api looking for places (not tried) ANy help appreciated, Thanks.


title: What does SetHistogramData do?

text: I run the following code:               IClassify classify = new NaturalBreaksClass();     ITableHistogram tableHistogram = new BasicTableHistogramClass();     IBasicHistogram basicHistogram = tableHistogram as IBasicHistogram;     tableHistogram.Table = table;     tableHistogram.Field = fieldName;     if ( !string.IsNullOrEmpty( normName ) ) tableHistogram.NormField = normName;     object dataFrequency;     object dataValues;     basicHistogram.GetHistogram( out dataValues, out dataFrequency );     double[] data = dataValues as double[];     int[] freq = dataFrequency as int[];     classify.SetHistogramData( data, freq );      but when I examine the classify object SetHistogramData() appears to have done nothing despite the fact that I can clearly see over 2000 items in both the data and freq arrays. I can read the values from the array directly and set them in the ClassBreaksRenderer using set_Break() but if that is the correct way then why should I call SetHistogramData() at all? What does SetHistogramData do?


title: GeoPDF Library for .NET or JavaScript

text: Anyone know of any open source .NET or Javascript libraries for creating GeoPDFs? I know of QuickPDFLibrary, but i'm preferably looking for something open source. I want to programmatically generate some GeoPDFs.


title: Map generation failing with Kartograph.js

text: I am trying to generate a very simple map using kartograph.js. But some how my code is not working.                   <html>     <head>          <script src="jquery-1.8.3.min.js"></script>          <script src="raphael-min.js"></script>          <script src="kartograph.min.js"></script>               <!---starting karograph.js-->                 <script language="JavaScript">     function loadmap(){     var map = Kartograph.map('#map', 900, 0);          map.loadmap('map.svg', function() {         // add layers.             map.addLayer('countries');             map.addLayer('depth');             map.addLayer('trees');             map.addLayer('crops');             map.addLayer('rivers');             map.addLayer('canals');             map.addLayer('city');             map.addLayer('states');             map.addLayer('coast');             map.addLayer('coast');     });     }     </script>     </head>     <body onLoad = "loadmap()">     <div id = "map">     </div>     </body>     </html>      It basically shows a blank page. I am not sure If this has to be hosted on a local server? however I setup a local server using Python but it still shows a blank page. I checked all the source file addresses and they are all right. Not sure whats happening here?


title: Trouble adding tabs in identity popup JS API

text:
**May 28 Edit** I have edited my code a bit and have come up with a fix to the tab container, however I am now having troubles getting my attributes to show when I click on one of my service layers. This is what I have so far                   var infoWindow = new InfoWindow(null, domConstruct.create("div"));         infoWindow.startup();         map.infoWindow.resize(300,300);          //Create info template to hold the content inside     var vliTemplate = new esri.InfoTemplate();     vliTemplate.setTitle("VLI");     vliTemplate.setContent(getWindowContent);          function getWindowContent(graphic) {     //VQO     //make a tab container     var tc = new dijit.layout.TabContainer({     style:"width:100%;height:100%;",     useMenu: false,     useSlider:false,     }, dojo.create("div"));          //add content pane 1 and 2 - these are the tabs     var cp1 = new dijit.layout.ContentPane({     title: "Current",     content: "Visual:
VLI Number: " + graphic.attributes.Tbl_ResultSummary_VQO_CFLBarea + "
	"     });     var cp2 = new dijit.layout.ContentPane({     title: "Approved",     content: "

VQO: ${Tbl_ResultSummary_VQO_EVQO_CD}
	"     });               tc.addChild(cp1);     tc.addChild(cp2);     return tc.domNode;     }      As you see my cp1 and cp2 have different ways of showing my attributes (in order to try and find something that works) cp1 "Current" tab when clicked on reads "Undefined" as the attribute information, and cp2 "Approved" tab shows text written as **VQO:** Tbl_ResultSummary_VQO_EVQO_CD Now i'm trying to figure out how to properly get my attributes to display in the appropriate tab pane! **May 22** I have a ArcGIS JS API map created with an InfoTemplate popup to display attributes but I need to have a tab in my popup window so that I can display different attributes from the same service. For example on one tab I will have "Current" attributes, and in the other tab I will have "Approved" attributes. This is my code for my InfoTemplate, but I don't know how to wrap the tab function around what I currently have displayed.               var vliTemplate = new InfoTemplate(     "Visual", "
VLI Number: ${VLI_VLI_NO}
	"+     "

VQO: ${Tbl_ResultSummary_VQO_EVQO_CD}
	"+     "

VAC: ${Tbl_ResultSummary_VQO_VAC_FIN_CD}
	"+     "

CFLB Area: ${Tbl_ResultSummary_VQO_CFLBarea:NumberFormat(places:0)} Ha
	"+     "

THLB Area: ${Tbl_ResultSummary_VQO_THLBarea:NumberFormat(places:0)} Ha
	"+     "

Total Area: ${Tbl_ResultSummary_VQO_HA:NumberFormat(places:0)} Ha
	"+     "

Greened Up Height: ${Tbl_ResultSummary_VQO_Height:NumberFormat(places:0)}m
	"+     "

Not Greened Up CFLB: ${Tbl_ResultSummary_VQO_Not_Greened_Up_Ha:NumberFormat(places:0)} Ha
	"+     "

Greened Up CFLB: ${greened:compare6}Ha
	"+     "

Maximum Target Percent: ${Tbl_ResultSummary_VQO_Pct_tgt:NumberFormat(places:0)}%
	"+     "

Maximum Target: ${Tbl_ResultSummary_VQO_Target_Ha:NumberFormat(places:0)} Ha
	"+     "

Surplus/Deficit: ${surplusDef1:changeColor6}
	")          compare6 = function (value,key,data) {      var result = "", diff;           switch (key) {         case "greened":           diff = data.Tbl_ResultSummary_VQO_HA - data.Tbl_ResultSummary_VQO_Not_Greened_Up_Ha;           break;          }             result = diff.toFixed();        return result;          };      I have been looking at the sample on https://developers.arcgis.com/javascript/jssamples/find_drilldown.html but I'm stumped....


title: Identify Toggle Layers with Attachments

text: I was looking at this fiddle http://jsfiddle.net/neWnE/1/ and trying to see get attachments to work. Since it uses the popuptemplate I should be able to add the showAttachments: true to it to get attachments to show in the popup window correct?               var template = new PopupTemplate({             title : "Name: {NAME}",             fieldInfos : [{               fieldName : "ADDRESS",               visible : true,               label : "Address:"             }],             showAttachments: true           });      I replaced the map services in the fiddle with my own that I enabled attachments on but they aren't appearing. Is there something else that I would need to do?


title: Layers with same coordinate system should align/overlap but do not?

text: I put 4 different layers into my ArcMap document, but only one of them will show up on the map! I should be seeing the other 3 layers as well. When I do a _Zoom to layer_ , I can see one layer, but I cannot figure out how to make the other 3 show up at once on top of each other. Before I opened ArcMap, I used ArcCatalog where I right-clicked 3/4 of the shapefiles and clicked _Properties_ -> _Coordinate system_ -> _Import_ and then clicked on my Trees layer, which has the coordinate system I want the other 3 to have. But when I opened ArcMap the 4 layers did _not_ overlap, nor could I move them above or below each other.


title: open shapefiles interpreted in Qgis, in ArcGIS

text: I set up a project with a Landsat TM image, and interpreted three shapefile layers with different line interpretations. I now need to open the layers in ArcGIS. The Landsat will open, but the shapefiles won't overlay. They are not "holding" the projection of the image file when I try to open them. I need help!


title: How to project US States in Transverse Mercator

text: > **Possible Duplicate:**   >  Layers having the same coordinate system should overlap, but they don’t. > What’s wrong? I created two different MXDs, one for City of Tucson and another for US States. If you click on the ArcGIS.com Map link above (in the services link page), you will see that City of Tucson Layer is overlaid on the basemap but the US States is projected somewhere in the ocean. I know this is a projection issue but if you look at the metadata in the services link, both are showing Spatial Reference: 102100 (3857). Now the projection of the City of Tucson is NAD_1983_HARN_StatePlane_Arizona_Central_FIPS_0202_Feet_Intl in the ArcMap and the US States is WGS_1984_Web_Mercator_Auxiliary_Sphere in ArcMap. I tried to change the projection to WGS84 but the layer is still not projected properly. How can I correct it? By "correct it" I mean that the projection should be like City of Tucson where it is projected properly over the basemap in ArcGIS.com. **UPDATE: Projection Tool used** I used the Project tool in the Projections and Transformations Toolbox, which creates a new Feature Class as per the Co-ordinate system we provide. So the original data had co-ordinate system of WGS_1984_Web_Mercator_Auxiliary_Sphere and the new one I assigned was GCS_WGS_1984. The new output was saved in a new File Geodatabase. The Spatial Reference that is now reflected in the Meta Data is 4326, which is not what I wanted. I wanted it to be same as City of Tucson which was 102100 (3857).


title: Projection Issues

text: I am working with NC shapefiles, I downloaded a national cities shapefile and I can't seem to line them up, actually you can't even see the layers in the same view. I have tried many times to re project in Arc Catalog, I even took the layers .prj file and copied and pasted it into the cities .prj file and it still did not work. I made the layer undefined and brought it into a projected data frame thinking it would project on the fly and that did not work. I projected it in WGS1984 thinking arc will know how to make it line up with NCs' state plane and that didn't work either. The only difference I have noticed comparing projection data is that the source extents are different. I can't see anyway to change them.   an example of the difference when I am in ArcMap and pull up it's source properties   the file where the projection is correct Left: 363485.468756 m the file where the projection is wrong Left:-178.877400 m any help would be greatly appreciated! Niki


title: Defining projection, reprojecting, and wrong area values

text: Apologies for duplication--this was mentioned in a previous post, but was suggested that this specific problem should be posted separately. My data: http://water.usgs.gov/GIS/dsdl/aquifrp025.tar.gz (USGS circa 2003 data on aquifers in the contiguous U.S.) Program: ArcGIS 10.1 The shapefile has no coordinate system or spatial reference. The metadata specifies NAD 1983 and in decimal degrees. My goal is to render this in such a way as to obtain 'normal' area/length/etc. values in meters for the purposes of calculating overlays with other spatial data. I cannot project the data because it has no coordinate system (error: no input coordinate system specified; paraphrasing)...batch/define doesn't work correctly either (the former same error, the latter (technically for use with rasters?) defines the projection, but reprojecting using the project tool doesn't do anything substantive (i.e. areas for huge swaths of the U.S. "46.1" etc. etc.). Running out of patience and options. Any help=appreciated.


title: Layers with the same Geographic Coordinate system don't overlap

text: I am a GIS newbie. I'm having a problem that, from the looks of things on this website, is a fairly common problem. I have 3 layers, all in the same GCS that do not line up. I have read all the the strings that others have posted on this website regarding this kind of problem, but I still haven't been able to figure out how to fix my particular issue. I am just a social studies student, trying to do a project. I don't have GIS training, so I feel very lost at sea, with all the complex technical terminology. Here are the details: I am working with data for the City of Chicago, USA. I have a street map for Chicago (chicagostreets1999, derived, I think from the US Census), a map of census tracts for the year 2000 (tracts2000, from the US census), and a zoning map (Zoning_nov2012, provided online by the Chicagy City Hall). Census tract boundaries run along streets. The tract boundaries in tracts2000 line up nicely with the street map from chicagostreets1999. Zoning boundaries also follow streets, but the zoning map does not line up properly with these two other files. When I compare the zoning map to the street file, its lines are about 175 feet (53 meters) south, and 148 feet (45 meters) east of the relevant lines on the street map. I looked at the "source" information on each file in arcGIS. Here is the info on their Geographic Coordination Systems. _**_**_Tracts 2000_ ** _*_ ** _*_ *** Data Type: Shapefile Feature Class Geometry Type: Polygon Geographic Coordinate System: GCS_North_American_1983 Datum: D_North_American_1983 Prime Meridian: Greenwich Angular Unit: Degree ![Misalignment, zoom in][1\ _**_*** _Chicago streets 1999_ ** _*_ ***** Data Type: Shapefile Feature Class Geometry Type: Line Geographic Coordinate System: GCS_North_American_1983 Datum: D_North_American_1983 Prime Meridian: Greenwich Angular Unit: Degree _**_**Zoning_nov2012* _*_ ** _*_ ** _*_ ** _*_ ** _*_ ** _*_ **** Data Type: Shapefile Feature Class Geometry Type: Polygon Projected Coordinate System: NAD_1983_StatePlane_Illinois_East_FIPS_1201_Feet Projection: Transverse_Mercator False_Easting: 984250.00000000 False_Northing: 0.00000000 Central_Meridian: -88.33333333 Scale_Factor: 0.99997500 Latitude_Of_Origin: 36.66666667 Linear Unit: Foot_US Geographic Coordinate System: GCS_North_American_1983 Datum: D_North_American_1983 Prime Meridian: Greenwich Angular Unit: Degree I tried setting the coordinate system for all three files by double clicking the "Layer" meta-category in the data view window, and then selecting "Coordinate System" within the data frame properties window. The coordinate system was already defined as: _**_*** _Layers_ ** _*_ ***** NAD_1983_StatePlane_Illinois_East_FIPS_1201_Feet Projection: Transverse_Mercator False_Easting: 984250.000000 False_Northing: 0.000000 Central_Meridian: -88.333333 Scale_Factor: 0.999975 Latitude_Of_Origin: 36.666667 Linear Unit: Foot_US GCS_North_American_1983 Datum: D_North_American_1983 This seemed to be the same as all three layers, so I did not change it. I asked the GIS librarian at my university what the problem might be and what, if anything could be done, to fix the misalignment. He said that he did not know how to fix it, and that the problem was probably due to the City Hall GIS people scanning in their zoning map "at a different scale" than the Census people did. I don't know whether that is true, but I really would like to get these files to align properly with each other. It will make the process of analysis a lot easier (ie not horribly visually confusing), and also more accurate. I am happy to share the files in question with anyone who might be able to help. The Zoning data, for instance, is freely available here on the City of Chicago's City Hall government website: https://data.cityofchicago.org/Community-Economic-Development/Boundaries- Zoning-Districts/p8va-airx I can email the street file to anyone who is willing to help. Thanks so much for your advice!


title: Spatial data not aligning after Define Projection?

text: I have three shapefiles and two buffers made from two separate shapefiles. I was trying to do a spatial join and it kept coming up with an empty file even though there was clearly overlap. I realized that some of the shapefiles had different projections. I changed them all to NAD 1984 for zone 18 and I also tried WGS, and then GCS WGS, but two of the shapefiles and a buffer file are now unaligned with a shapefile and its buffer. They are all the same coordinate system and they aren't aligning when they were before I changed the coordinate system. Is there a way to change them back? The undo button just says "remove layer." Or suggestions to fix this? I don't want to start from scratch with the original data.


title: Why did my shapefile distance get smaller after changing the projection?

text: > **Possible Duplicate:**   >  Layers having the same coordinate system should overlap, but they don't. > What's wrong? I'm working with a shapefile that has a geographic coordinate system and I need to change it to a projected coordinate system to use spatial analysis tools. Steps: 1) I defined the coordinate system using the define projection tool to State Plane (US Feet) California Zone 6. 2) Then I used the Project tool to make define the coordinate system as the state plane system using a different name. I made sure not to use the same name and followed the guidelines listed here: http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//00130000004v000000.htm When I measured the county with the measurement tool, it said it was less than a mile across (it measured about 34 miles across when the coordinate system was defined using the GCS). I'm not sure what I'm doing wrong. Any help is greatly appreciated. Thanks!


title: Why is shapefile not drawing in expected location after (or before) using Define Projection tool?

text: I have **downloaded this zipped shapefile** of bus routes from the City of Houston data portal. I then performed these steps in **ArcGIS Desktop 9.3.1** :   1. Extracted the shapefile from the *.zip   2. Added this shapefile as a layer to a data frame which has data in a coordinate system of `NAD_1983_UTM_Zone_15N` and noticed that it did not draw where I expected it to   3. Used **ArcToolbox>Data Management Tools>Projections and Transformations>Define Project** to change its coordinate system to `NAD_1983_UTM_Zone_15N` and observed that it still did not draw where I expected it to Can you explain why this shapefile may not be drawing in the expected location? As background the properties of a shapefile that is drawing in the expected location is:               Data Type:  Shapefile Feature Class      Shapefile:  D:\Cloud Apps\Google Drive\...\tgr48201trt00.shp     Geometry Type:  Polygon          Projected Coordinate System:    NAD_1983_UTM_Zone_15N     Projection: Transverse_Mercator     False_Easting:  500000.00000000     False_Northing: 0.00000000     Central_Meridian:   -93.00000000     Scale_Factor:   0.99960000     Latitude_Of_Origin: 0.00000000     Linear Unit:    Meter          Geographic Coordinate System:   GCS_North_American_1983     Datum:  D_North_American_1983     Prime Meridian:     Greenwich     Angular Unit:   Degree


title: convert decimal degrees to meters

text: I have a shape file. The coordinates (Lat/Long) are projected as decimal degrees. But, I want to re-project it to meters. How can I do it with ArcMap ? Please Thank you.


title: NAD conversion in Mexixo

text: I am looking for the file "Mexico.gsb" which used with ArcGIS when datum shifting from NAD27 to NAD83 in Mexico.


title: When labelling polygons, not all labels are displayed - QGIS 1.7.0

text: I am trying to create a simple map of Mexico showing the state polygons and their label. The state shapefile is from here, and the specific file is the first one from the top (Áreas Geoestadísticas Estatales (7.43 Mb)). I edited the attribute table to create a new column with shorter state names called `NOM_SHORT`. To add labels I used the "ABC" dialog as below: ![enter image description here](http://i.stack.imgur.com/tzlxZ.png) However not all states were labelled, including Oaxaca and San Luis Potosí, as shown below (painted in light grey). Why? ![enter image description here](http://i.stack.imgur.com/fttEl.png) Here is attribute table: ![enter image description here](http://i.stack.imgur.com/lKriN.png)


title: arcgis server 10, vb.net, IMathOp problem

text: i am using arcgis server 10 and visual studio with vb.net to build web mapping application. inside my code i have the following line:               Dim RasterCal As ESRI.ArcGIS.SpatialAnalyst.IMathOp = _        New ESRI.ArcGIS.SpatialAnalyst.RasterMathOpsClass      or even if i do it this way               Dim RasterCal As ESRI.ArcGIS.SpatialAnalyst.IMathOp = _        New ESRI.ArcGIS.SpatialAnalyst.RasterMapAlgebraOp      in both cases i get initialization error which says:               The runtime application type must be specified before license      initialization.      plz can you tell me what is wrong with the above examples of lines? thanks!


title: Open Street Map: converting to country -> administrative unit -> city database

text: To moderators: If you think this question is too broad, please point me into right direction before closing, I'm **desperately** looking for help. I'm stuck on my task and I don't know what to do. The job is:   * extract data from Open Street Map   * create database model containing following schema and data Schema looks like this:               +---------+     +------------+      +----------+     | country |     | adm_unit   |      | city     |     +---------+     +------------+      +----------+     | id      |--\  | id         |--\   | id       |     | name    |   \-| country_id |   \--| adm_id   |     +---------+     | name       |      | name     |                     +------------+      | coord_x  |  //coordinates - longitude and latitude                                         | coord_y  |                                         +----------+      So far I managed to:   * download map itself (only for Poland to test it)   * import it to database using `osm2pgsql` tool   * perform some selects and realize, I don't know how to do it. The main problem is: I can't see any relations between different items of map. Taking Poland as example: I can find Warsaw (capitol city), as well as "masovian vovoidenship" performing few selects, but there is no relation between these two. There are lot of empty columns in my db instead. Also none of column names suggests that it should hold any relation to "parent". I don't know how to find those relations. Is it something with imports that i'm missing? Should I change anything in my `default.style` (I haven't changed anything inside it)? Or maybe there's some kind of other approach to solve my problem? Maybe there are free databases that contain data I need and are less painful to work with? `Geonames` is no-go, their database is full of errors. I'm really green when it comes to `gis`, I'm PHP programmer, that's completely different world for me. Unfortunately, even if I want to explore it, my deadlines won't allow that.


title: load table with data in ArcPad 10

text: I'm using ArcPad 10 and arcgis 10.1 I need to load a table with this information so that when lifting the geographic coordinate this information related to this coordinate taken with trimble gps ..... I am using ArcPad 10 and ArcGIS 10.1 I have a table "customers" dbf format. with information (name, address and customer code) I need to georeference the location of these customers and for that, create an electronic ticket for use in ArcPad 10 on my gps trimble My question is, if you can load the table "customers" in ArcPad then in the field when saving the coordinate I can associate with the information in the table


title: Esri GIS Bookmark but with out a dijit

text:
I need to create a "Bookmark" like drop down on my site. Inspiration comes from samples:1) http://developers.arcgis.com/javascript/samples/widget_bookmarks_dropdownbutton/ and 2) https://developers.arcgis.com/javascript/jssamples/widget_bookmarks_dropdownbutton.html and 3) http://developers.arcgis.com/javascript/sandbox/sandbox.html?sample=widget_bookmarks_editable (-minus the editing part). I'm using bootstrap's drop down menu for the book marks so I can't use the bookmark digit directly BUT the sample I provide above does exactly what I need it to do. Basically I want to click my drop down menu, and select from a predetermined list (each entry contains a xmin,xmax,ymin and ymax) and have the map pan and zoom to fit the coords on the screen. I cannot find any examples on how to do this without using a dijit. I can't really post code because I don't know where to start. In my mind, I picture something along these lines: `function ZoomToRegion(RegionID){ // console.log('zoom to xy coords'); var wkid = 102100; if (RegionID == 1){ var xmin = '-8849573.386742238'; var ymin = '4366083.055648041'; var xmax = '-8448431.862301651'; var ymax = '4573380.276357428'; var regionname = 'Central '; } else if (RegionID == 2){ var xmin = '-9323482.962110244'; var ymin = '4350184.153764758'; var xmax = '-8922341.437669657'; var ymax = '4557481.374474147'; var regionname = 'Far Sourthwest'; } ... not not sure what to put here for the panning and zooming to fit }` and then my html code is: `

Default View

Central

Eastern

Far Southwest

Near Southwest

Northern

Northwest

` This seems like it should be pretty straight forward but wihtout a non-dijit example, I'm at a loss. * * * **Here is what I ultimately went with:**

Default View

Central

Eastern

Far Southwest

Near Southwest

Northern

Northwest

               function ZoomToRegion(RegionID){     //  console.log('zoom to xy coords');            var wkid = 102100;           if (RegionID == 1){             var defaultExtent = new esri.geometry.Extent({                 "xmin": -8849573.386742238,                   "ymin": 4366083.055648041,                 "xmax": -8448431.862301651,                 "ymax": 4573380.276357428,                 "spatialReference": {                     "wkid": 102100                 }             });         }         else if (RegionID == 2){             //  Far Southwest             var defaultExtent = new esri.geometry.Extent({                 "xmin": -9323482.962110244,                   "ymin": 4350184.153764758,                 "xmax": -8922341.437669657,                 "ymax": 4557481.374474147,                 "spatialReference": {                     "wkid": 102100                 }             });              }         else if (RegionID == 3){             //  Near Sourthwest             var defaultExtent = new esri.geometry.Extent({                 "xmin": -8998166.969728548,                   "ymin": 4356910.61225385,                 "xmax": -8597025.445287962,                 "ymax": 4564207.832963239,                 "spatialReference": {                     "wkid": 102100                 }             });         }         else if (RegionID == 4){             //  Northern             var defaultExtent = new esri.geometry.Extent({                 "xmin": -8803099.673544865,                   "ymin": 4561150.351831801,                 "xmax": -8401958.149104279,                 "ymax": 4768447.572541188,                 "spatialReference": {                     "wkid": 102100                 }             });              }         else if (RegionID == 5){             //  Northwestern             var defaultExtent = new esri.geometry.Extent({                 "xmin": -8988994.526334316,                   "ymin": 4450163.786761777,                 "xmax": -8186711.4774533175,                 "ymax": 4864758.228180463,                 "spatialReference": {                     "wkid": 102100                 }             });              }         else if (RegionID == 6){             //  Eastern             var defaultExtent = new esri.geometry.Extent({                 "xmin": -8700368.307529585,                   "ymin": 4374644.002815981,                 "xmax": -8299226.783088999,                 "ymax": 4581941.223525369,                 "spatialReference": {                     "wkid": 102100                 }             });              }         else if (RegionID == 7) {             var defaultExtent = new esri.geometry.Extent({                 "xmin": -9274930.161743507,                   "ymin": 4318875.546979173,                 "xmax": -8303140.358937347,                 "ymax": 4793763.5163091915,                 "spatialReference": {                     "wkid": 102100                 }             });         }         cppr.App.map.setExtent(defaultExtent.getExtent(),true);     }


title: IQueryDescription SQL Query Returning Limited Records

text: I created a SQL query that is working in SQL Management Studio, but as soon as I transfer the query into ArcMap's IQueryDescription to try and run it with a User Form, the results that I'm getting back are very limited and there is no clear pattern to the results it is returning. (409 records vs. 15 records) I am even copying and pasting the query that works into the ArcObjects code, but I get a limited number of records back and no errors are thrown. Has anyone run into this? What direction should I look for a solution? I included the query below, although the query itself works just fine in SQL Management Studio.               SELECT CADDATA.CALLINDEX.NCALLHIST2, CADDATA.CALLINDEX.NHISTSEQUENCE, CADDATA.CALLINDEX.SZCALL, CADDATA.CALLINDEX.SZCALLTYPE, CADDATA.CALLINDEX.SZCALLDESC, CADDATA.CALLINDEX.SZGEOGROUP, CADDATA.CALLINDEX.SZPRIORITY, CADDATA.CALLINDEX.SZDISPOSITION, CADDATA.CALLINDEX.LTCREATED, CADDATA.CALLINDEX.LTENTERED, CADDATA.CALLINDEX.LTDISPATCHED, CADDATA.CALLINDEX.LTENROUTE, CADDATA.CALLINDEX.LTONSCENE, CASE WHEN RIGHT(SZLOCATION, CHARINDEX(',', SZLOCATION)) = '' THEN SZLOCATION ELSE LEFT(SZLOCATION, CHARINDEX(',', SZLOCATION) - 1) END AS Location, cast(LTCREATED as date)  as gedatedidid FROM CADDATA.CALLINDEX INNER JOIN CADDATA.CALLINDEX_MXSEQUNCE ON CADDATA.CALLINDEX.NCALLHIST1 = CADDATA.CALLINDEX_MXSEQUNCE.NCALLHIST1 AND CADDATA.CALLINDEX.NCALLHIST2 = CADDATA.CALLINDEX_MXSEQUNCE.NCALLHIST2 AND CADDATA.CALLINDEX.NHISTSEQUENCE = CADDATA.CALLINDEX_MXSEQUNCE.NHISTSEQUENCE WHERE (CADDATA.CALLINDEX.SZCALLTYPE IN ('1818', '1825', 'AH', 'BC', 'BS', 'BUS', 'CH', 'F', 'FP', 'PAT', 'PC', 'PW', 'SUB', 'TL', 'TS', 'TST', 'VHC', 'VP')) AND (CADDATA.CALLINDEX.SZGEOGROUP IN ('D1', 'D2', 'D3')) AND (cast(LTCREATED as date) >= '2013-6-4') AND (cast(LTCREATED as date) <= '2013-6-4')


title: Add Numbers to United States Map Help R, R-Studio?

text: Can someone help me with creating a map in `R`? I need to make a map of the number of retail stores a company has in each state. All I need is the store count displayed in each state. For example: if there are 150 stores in California and 250 in Texas, I need the map to display 150 within the California borders and 250 within the Texas borders. I created a map with:               m = map("states", region = c('Alabama', etc....))      I tried adding the numbers as a label but that didn't work.               m = map("states", region = c('Alabama', etc....), labels = c('24','56',etc....))      I must be missing something. I have a .csv file with the state name and the number of stores in each state. Is there an easy way to translate this into a graph? I'd like to avoid doing it by hand in Excel. * * * Thanks! That worked. Does anyone know how to deal with the labels that do not fit within a states boarders? Also, does anyone know how to color code states based on the number value? The image of the map and code I am using is below: ![US Map With Counts](http://i.stack.imgur.com/Nt2mo.png) `library("maps")` `map.text("state", regions=c("alabama", "arizona", "arkansas", "california", "colorado", "connecticut", "delaware", "district of columbia", "florida", "georgia", "idaho", "illinois", "indiana", "iowa", "kansas", "kentucky", "louisiana", "maine", "maryland", "massachusetts:main", "michigan:north", "minnesota", "mississippi", "missouri", "montana", "nebraska", "nevada", "new hampshire", "new jersey", "new mexico", "new york:main", "north carolina:main", "north dakota", "ohio", "oklahoma", "oregon", "pennsylvania", "rhode island", "south carolina", "south dakota", "tennessee", "texas", "utah", "vermont", "virginia:main", "washington:main", "west virginia", "wisconsin", "wyoming"), labels=as.character(c(29, 38, 13, 173, 21, 12, 5, 1, 108, 59, 7, 40, 46, 3, 9, 40, 23, 8, 24, 20, 49, 9, 14, 27, 2, 4, 12, 7, 21, 12, 60, 73, 2, 108, 18, 14, 71, 1, 35, 1, 49, 117, 10, 4, 39, 27, 18, 10, 2)))`


title: How to collect tweets outside a given bounding box?

text: I checked the twitter API here. I know how to get tweets in United States within a bounding box such as:               [-124.47,24.0,-66.56,49.3843].      How can I easily get all the tweets NOT in this bounding box?


title: Arcobjects Multithreading in C++ - 'Aborted' when modifying Geometry

text: My project is required to be cross-platform, so I am using the C++ API for ArcObjects. My application creates a new thread which listens to a network socket broadcasting location information. I am attempting to read in the information and add points to my DynamicLayer. My problem comes when I try to use a new IPoint in the second thread. I can make the call:               IPointPtr newpoint(CLSID_Point);      But on the next line when I execute:               newpoint->PutCoords(x,y);      The application fails with a nonspecific error 'aborted'. Relavant information: I am creating a new thread by calling `pthread_create` on a static method. I have read the ESRI document on multithreading (which like the rest of the documentation is woefully absent of c++ specific information), but didn't see anything I thought could be the problem. Any thoughts?


title: What each stack element does

text: Im using the opengeo stack, and want to double check my understanding of each element.   * PostGIS - The database for storing my data   * GeoServer - The link between the online part / the database Here's where I am not fully sure:   * Geoexplorer, viewing, editing data visually in the POST database   * Client SDK - The same as geoexploer but I have to add the functionality? If I want to add Graphs to the map, i.e. the user clicks on a point, how would I go about that? using SDK and the GeoEXT?


title: Explicitly locking a file geodatabase

text: I need to lock a file geodatabase (or only a single Mosaic dataset) at the beginning of a python script (and release it in the end).. I looked and looked and couldn't find any explicit arcpy command to do this.. Is it even possible to tell the gdb to lock itself? Or is the locking mechanism works autonomously?


title: GeoServer Layer from PostGIS not updating

text: I created a WMS layer from a PostGIS table with a fraction of total polygons in GeoServer. Everything displayed fine. Then, I added the rest of the polygons to the table in postgis. For some reason its not updating the layer. I tried restarting geoserver and its still not updating. How do I get geoserver to update the layer? I use Open layers WMS layer and OpenLayers.Protocol.WFS.fromWMSLayer(polygons) to get polygon data. -N UPDATE: I tried to click Compute from data/Compute from bounding box links to get the new MBR. Unfortunately, it did not update, so here's as much info as I can think to provide. I have a set of polygons scattered around the US. They are from a pair of SHP files in WGS_1984. I started testing the system with a small set of 100 polygons and now that the system works, I wanted to add the the other polygons to the data. I used shp2sql to create a new table in postgis. Then I used a INSERT INTO polygons (SELECT from ...) to add them to the table that geoserver uses. Now there is about 50k polygons in the table. I tried updating the MBR as answers below suggest. Here are the current settings of the layer in geoserver: Native SRS: EPSG:4326 Declared SRS: EPSG:900913 When I hit Compute from data for Native bounding box it gives me: (Min x,Min y, Max x, Max y) -10Mil, 5Mil, -10mil, 5Mil and for the Declared: -95.xxx,40.xxx,-95.xxx,43.xxx  For whatever reason, they are now updating correctly. The MBR should cover most of the US now and not just that little portion in the middle.


title: How to use ICalculator with a pre-expression in C#

text: Another wonderful predicament. I am trying to use an instance of `ICalculator` to update a field based on another field. Textbook stuff. Unfortunately, the example in the SDK for C# is a bit _lacking_. Here is what I have puzzled out so far from the C# GeoDatabaseUI 10.x and VB.NET ICalculator 9.3 examples/documentation ( **NOTE** : I am doing this in C#, in case I was not clear).   1. If you want to use a string in the Expression property, do **_not_** use triple quotes (must be legacy information), but escaped quotes for normal string-within-string notation.    2. You can use VB.NET code in the pre-expression when working in VB.NET (this confuses me). With these assumptions in mind, I have been attempting to use python or C# in the pre-expression property, but I am having no luck. `snip....` Has anyone tried this yet? If so, what is the proper way to use a python function and/or a C# method when working with the `ICalculcator`? **UPDATE:** Okay, trying the method kenbuja suggests below, I changed it to an `IGeoProcessor` method with the following code:               public void GPCalculate(ITable table)         {             IGeoProcessor gp = new GeoProcessorClass();             IVariantArray varray = new VarArrayClass();             string input = workspacepath+"\\"+targetfile;                  string codeblock = "states = {'2':'AK', '1':'AL', '5':'AR', '60':'AS', '4':'AZ', '6':'CA', '8':'CO', '9':'CT', '11':'DC', '10':'DE', '12':'FL',"+                 "'13':'GA', '66':'GU', '15':'HI', '19':'IA', '16':'ID', '17':'IL', '18':'IN', '20':'KS', '21':'KY', '22':'LA', '25':'MA', '24':'MD', '23':'ME',"+                 "'26':'MI', '27':'MN', '29':'MO', '28':'MS', '30':'MT', '37':'NC', '38':'ND', '31':'NE', '33':'NH', '34':'NJ', '35':'NM', '32':'NV', '36':'NY',"+                 "'39':'OH', '40':'OK', '41':'OR', '42':'PA', '72':'PR', '44':'RI', '45':'SC', '46':'SD', '47':'TN', '48':'TX', '49':'UT', '51':'VA', '78':'VI',"+                 "'50':'VT', '53':'WA', '55':'WI', '54':'WV', '56':'WY'}\n" +                 "def getstate(id):\n" +                 "  id = str(id)\n"+                 "  if (len(id) == 14):\n" +                 "    id = id[0]\n" +                 "  if (len(id) == 15):\n" +                 "    tid = id[0:2]\n" +                 "  return states[id]";             // First Parameter in_table             varray.Add(input);             // Second Parameter calculate field             varray.Add("PSTATE");             // Third Parameter expression             varray.Add("getstate(!ID!)");             // Fourth Parameter expression_type (optional)             varray.Add("PYTHON_9.3");             // Fifth Parameter Codeblock (optional)             varray.Add(codeblock);                  IGeoProcessorResult result = new GeoProcessorResultClass();             settesttext("Starting Calculation");             result = gp.Execute("CalculateField_management", varray, null);         }      The above code works perfectly now (thanks to an edit @blah238).


title: How to search an ITable with conditions

text: At the moment I am searching for the most effective method of searching an ITable with one or more conditions (like where clauses). I can perform it manually by iterating through each row and proof manually if the field "xy" has the value "ab". However, because the table has tens of thousands of row this will not be really effective. Does anyone has a effective and elegant solution? Thanks for your help!


title: How to calculate statistics for several subcatchments (polygons) in the same time?

text: I have a shapefile of groundwater management zone which is composed of 40 subcatchments (40 polygons). I'd like to calculate the max., min. and mean elevation for each subcatchment. Would you please let me know if there is any way to calculate all of these parameters (for all of the subcatchments) using the shapefile and the DEM? Thanks in advance!


title: Dealing with bridges and tunnels when modeling raster least cost paths?

text: I'm currently using ArcGIS 10 Spatial Analyst's Cost Path tool to model pedestrian routes. I have sidewalks, fences, buildings and a ton of pedestrian scale features, all combined into a cost surface in which sidewalks are low cost, parks are slightly higher, crosswalks are still higher, and buildings, fences, and the like are barriers. The model works great to define optimal routes, compare routes to each other, assess the walkability of areas, and build walkingsheds. But I have one major problem: bridges and tunnels. Because the cost surface forces all the features into a single 2D plane, I'm unable to model a bridge with a sidewalk that crosses over a street with a sidewalk. I must choose which one should get cut off. An alternative is to find all the bridges and tunnels in the City and model all combinations of under/over, but the model takes forever to run as it is! Is there a third party tool that can perform in least cost path in 3D, or am I missing a workaround of some kind?


title: How is Canada Vector Data (CanVec) hydrology data collected and updated?

text: Canada is vast and collecting any geographic data must be a monumental task. Natural Resources Canada provides a free access to a comprehensive Canada-wide topographic data now refereed to as CanVec (Canada Vector) data. I often download and use CanVec datasets which are updated on yearly basis (current version 12 soon to be version 13) but I am curious how current thi data really is, namely hydrology. How is this data captured and collected and how frequently are features such as waterbodies, wetlands updated? I suspect image analysis using multispectral imagery but clearly some manual work is needed to enforce hydrological correctness. Is is even possible such data can be up to date?


title: geotools / jts crs transformation result not at correct position

text: I copied the example from geotools, which shows how to transform a shapefile into a different projection. (http://docs.geotools.org/latest/userguide/tutorial/geometry/geometrycrs.html) But when transforming from Gauss Krueger 4 into Pseudo Mercator (3857) the shapes are not at the correct position, they are displaced several meters. What can I do to fix this? The sample code contains :                CoordinateReferenceSystem dataCRS = schema.getCoordinateReferenceSystem();                 CoordinateReferenceSystem worldCRS = map.getCoordinateReferenceSystem();                 boolean lenient = true; // allow for some error due to different datums                 MathTransform transform = CRS.findMathTransform(dataCRS, worldCRS, lenient);      ... and then ...               Geometry geometry2 = JTS.transform(geometry, transform);      additional information: The sample data of the source as the following CRS               PROJCS["PD_83_3_degree_Gauss_Kruger_zone_4",        GEOGCS["GCS_PD/83",          DATUM["D_Potsdam_Datum_83",            SPHEROID["Bessel_1841", 6377397.155, 299.1528128]],          PRIMEM["Greenwich", 0.0],          UNIT["degree", 0.017453292519943295],          AXIS["Longitude", EAST],          AXIS["Latitude", NORTH]],        PROJECTION["Transverse_Mercator"],        PARAMETER["central_meridian", 12.0],        PARAMETER["latitude_of_origin", 0.0],        PARAMETER["scale_factor", 1.0],        PARAMETER["false_easting", 4500000.0],        PARAMETER["false_northing", 0.0],        UNIT["m", 1.0],        AXIS["x", EAST],        AXIS["y", NORTH]]      and the target has the following CRS:               PROJCS["WGS 84 / Pseudo-Mercator",        GEOGCS["WGS 84",          DATUM["World Geodetic System 1984",            SPHEROID["WGS 84", 6378137.0, 298.257223563, AUTHORITY["EPSG","7030"]],            AUTHORITY["EPSG","6326"]],          PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]],          UNIT["degree", 0.017453292519943295],          AXIS["Geodetic latitude", NORTH],          AXIS["Geodetic longitude", EAST],          AUTHORITY["EPSG","4326"]],        PROJECTION["Popular Visualisation Pseudo Mercator"],        PARAMETER["semi_minor", 6378137.0],        PARAMETER["latitude_of_origin", 0.0],        PARAMETER["central_meridian", 0.0],        PARAMETER["scale_factor", 1.0],        PARAMETER["false_easting", 0.0],        PARAMETER["false_northing", 0.0],        UNIT["m", 1.0],        AXIS["Easting", EAST],        AXIS["Northing", NORTH],        AUTHORITY["EPSG","3857"]]      Any idea? D3


title: Geoserver Tile caching with multiple instances

text: I work with `GIS` Project and I used `Geoserver`. I want to use caching capabilities but caching working only one instance. when log out and log in again cached tiles are empty. Are there any suggestion?


title: Get shape filename in attribute field using field calculator

text: I apologize for crossposting against the ESRI geoprocessing forum, but it seems like this one gets answers. I'd like to have a modified append tool (ArcToolbox) that takes multiple shapefiles with point attribute data, and copies the inputted shapefile's name(s) to a "filename" attribute column, and then appends it to a target dataset. That way I can track where the data came from. The only way I can think to do this is something along the lines of the following:               for each shp in featureclasses:       outFeatureClass = os.path.join("C:\GIS\Scripts\temp\temptrks", shp.strip(".shp"))       arcpy.CopyFeatures_management(shp, outFeatureClass)       arcpy.AddField_management(outFeatureClass, "filename", "TEXT", "","", "254")       arcpy.CalculateField_management(outFeatureClass, "filename", str(shp), "PYTHON")      Then perform the append for each file in the new workspace. What I'd like to do is be able to have a user drag and drop multiple shapefiles from multiple directories into a Multiple Value input box (like that in the Append tool). I have two questions:   1. Is there an easier/better way to perform the task?   2. If not, what should I set the parameter type to (Multiple Value is not a selection). Any advice is appreciated. Thank you all


title: Viewing / Adding Attachments in Features Attributes - Possible or Not in MOSS?

text: Is it possible to add / view media attachments (photos) in an SDE Feature Class? I'm using the Editing Widget Tool which works. The problem is that I can't add attachments or view existing attachments. I just want to know if attachments are a possibility in Sharepoint Feature Classes or Shapefiles. As a bonus, what could be my issue or how do you enable attachments? My attachments work in ArcGIS.com editing / viewing.


title: ArcPy in_memory\ vanishing from ToC

text: ArcGIS 10.1 - attempting to add a derived raster layer to the ToC. This works in 10.2. In 10.1 I get a screen flash and nothing added.               rasterized = r'in_memory\rasterized'     arcpy.PolygonToRaster_conversion(poly_shp, field, rasterized, None, cellsize)     res = arcpy.MakeRasterLayer_management(rasterized, 'tmplyr')     rasterlyr = res.getOutput(0)          mxd = arcpy.mapping.MapDocument("Current")     df = arcpy.mapping.ListDataFrames(mxd, '*')[0]     arcpy.mapping.AddLayer(df, rasterlyr, "TOP")      Where: poly_shp = unicode vector name field = unicode field name from the poly_shp rasterized = str cellsize = float Adding:               arcpy.RefreshTOC()     arcpy.RefreshActiveView()      has no effect. Edit: * * * I know that the in_memory workspace is functioning properly as the following try is hit on subsequent runs.               try:             arcpy.Delete_management(rasterized)             arcpy.AddMessage('Overwriting in_memory raster')     except:             pass


title: How do you Find Linear Mileage?

text: I have never doen this before, but how do you find the linear mileage of microseismic lines? I tried researching it online in differant forums/tutorials with no luck. I have ArcView 9.3, no other extensions. I went into the Toolbox Linear Referencing Box and from there, I tried a few of the tools, but no luck. This is for PostPlot source points (Microseismic survey). Simply put, as above, I need to find the linear mileage with the postplot points that goes on/through the surface owners. How do I do this? Be advised, I am not a programmer nor a techie. if at all possible :) keep any answers as simple as possible unless it is unavoidable. If you need more information, please let me know as I did my best to provide pertinent information. ## I have included a pdf of the data if this helps any. A cheap and easy way to do this is use the Measure Tool. But, I figure that HAS TO BE A better way to do it then individually for each parcel. Is there? thank you, Z ![enter image description here](http://i.stack.imgur.com/O5FJf.png) ![enter image description here](http://i.stack.imgur.com/M4NiM.png)


title: Is the Intersection Operation Commutative?

text: In map overlay, is intersection of map layers a commutative operation ie layer_A intersection layer_B = layer_B intersection layer_A? If layer A is a square and layer B is a triangle. In this case, (A intersectn B) == (B intersectn A) ?


title: Building orientation and simplification using Opencarto Java library

text: I need to find out orientation of buildings and perform simplification of sides of buildings with a tolerance. I came to know from StackExchange that OpenCarto has useful libraries for geometrical functions.Basically i have stored data in PostGIS. Is there any possibility to use Opencarto Java Library to perform this on buildings? I thinks it is possible to connect Java with PostGIS through JDBC driver. Cheers, Sanjeewa.


title: Coords not aligning

text: The two shp files are not aligning like they should (note the img). The layer and the two shps all have the same coords - NAD_1983_StatePlane_North_Dakota_North_FIPS_3301_Feet The surface shp came with the defined coordinates so I imported those coords into the post plot which had "undefined" coords. They did not align. I checked the Layer coords and they were the same as the surface shp. I deleted the coords from the post plot shp file which originaly had "undefined coords" and started from square one with the post plot shp coords. I re-imported the coords and I get the same results. I used ArCatalog to put in the coordinates. I have Arcview 9.3 What am I missing or not doing? I've always been able to fix this issue except on this occasion. Thanks, Z ![enter image description here](http://i.stack.imgur.com/mte7V.png)


title: Updating Apps on OpenGEO suite Community Edition

text: I would like to update PostgreSQL and PostGIS on Opengeo CE to a newer version, maybe the other apps too (geoserver, openlayers, etc). Can it be done? Can someone point me the right direction? Thanks,


title: Predict user's future location - location awareness mobile computing

text: This question is a part of an online course (Mobility Data Management) that I am currently auditing. A part of the project is to implement a system that can predict user's future location. This kind of application has recently got some interests from database and mobile application developers. This is the problem statement: Assume that we are given a road network G (V,E), the recent motion history (trajectory) T, and current instant velocity v of a user, a temporal duration dt, and (optional) a historical trajectory database D of other users’ movements. Our goal is to estimate user’s location after dt. I found that the following paper has a good description about what a road network and the trajectory data is: http://www.sis.pitt.edu/~xliu/papers/ceus.pdf Does anyone has an idea on how to formulate a simple solution of this mobility awareness problem (predicting user's future location), and how to start the implementation (how you formulate the problem in an OOP way, such as: what classes are needed to represent the problem correctly).


title: How do I convert GML to ESRI featureclasses?

text: I have a series of xml files containing gml tags such as :               <gml:posList srsDimension="2">454780.3125 525867.9375 454762.0625 525994.5625      and I need to convert these into esri featureclasses. How would I start to import them? Thanks


title: Avoid backtracking while creating loop in Google Map API 3

text: I am trying to build a solution where i would create loop by google map. My first and only business rule is there should not be any repeated route. But with my demo http://www.mindcrew.in/road%20runner/example.html it does show. Please help me on this .


title: Geocoding through a language barrier

text: This is related to a previous question of mine. It has come to my attention that the streets file I use as my address locator in Canada uses French (or simply accented) characters for the names. These show up in similar manners to the following:               QuÃ©bec     De la RiviÃ re-du-Lac     FranÃºois-de Villars      The Address product, however, uses standard (Oops, my Ameri-centrism is showing?) English characters, causing the previous examples to display like so:               Quebec     De la Riviere-du-Lac     Francois-de Villars      I believe the geocoding would fair much better if the street names actually matched via the same character set. How would you recommend I go about converting from one to the other? I have considered regex'ing it, but I thought I would ask around before I start. P.S. It would be nice if I didn't get garbage text from French characters, regardless. All the info is stored in Geodatabases, by the way, if that matters at all. I'm on Arc10, as well. P.P.S. It would appear that I can display accented characters just fine, but in the Shapefiles and GDB Feature Classes, they become garbage text. Is there some setting I am missing here? **CURIOUS:** It would seem that in one street product the names come up garbage, but in the other they display just fine.


title: How can I specify the ranges in the gdaldem's color configuration file?

text: I want to render an IMG file according to a specific renderer, and export to a jpeg File. I have asked a question about this, and got a satisfactory answer. The issue now is that I need to write a color configuration file which matches the symbology that is required. My Symbology is supposed to be like this: ![enter image description here](http://i.stack.imgur.com/IevlT.png) Based on this answer, I have found that you can write the colors config file somewhat like this:               40  140 0   3     35  204 9   9     30  255 67  11     25  255 165 0     20  255 246 0     15  54  254 0     10  0   255 209     5   1   127 255     0   4   0   213     -4  133 1   243      My problem is that I don't understand if the given values for the pixels are the lower bound, or the upper bound. I don't know how to specify the ranges according to the symbology that I need. Where can I Find documentation for Gdaldem's color configuration file? Or alternatively, how can I specify the ranges in the configuration file?


title: configure mozswing in Java

text: I'm using OpenLayers - 2.11,GeoServer 2.13, ExtJS 3.4, and Java for developing GIS based web applications. My web map application is loading through "mozSwimg" in Java. While running application i'm getting Unresponsive Script warning in javascript file. Please guide me how to configure mozswing to solve this issue. Thanks..


title: Work flow for Lidar and OCAD

text: I am updating our maps for an orienteering club in Ilinois. I was very excited to find out that we have access to .las lidar data in 2000x2000ft files. I checked out OCAD11 and it can handle .las files so I thought I was golden. Except the .las files are in feet and OCAD11 does not handle feet and the conversion to utm is not right. I have contacted OCAD and they are unable to solve this issue. So I am looking for a possible work around or alternate work flow for develop georeference contours in utm. I have looked at Fusion and also the article "Using Lidar point clouds for Orienteering base map generation" by Terje Mathisen. I do not see a clean path for dealing with multiple .las files covering a park and also dealing with the Illinois state plane to utm. My current plan is to create a dem in Ocad11 export to a dem file. Import the File to qgis, do the conversation to utm, export to xyz file, import back to OCAD, create contours. Can someone suggest a cleaner simple work flow. I would like to avoid buying OCAD11. We currently use OCAD9. I would like to be able to make the contours externally and import them into OCAD9. Thanks in advance


title: Where to store geodata to be layered over an OSM tile map?

text: I'm working on a project for which most data can come straight from OSM. But some data is not relevant to OSM and will need to be stored separately then layered over the top. Probably the tiles will be generated in TileMill, and served with Leaflet. Either I'll do the layering within TileMill, or on the fly in Leaflet. The question is - where to store the extra geodata, and how to manage it? The data will mostly be POIs, some of which will have lengthy explanatory text and even photos associated with it. As a ballpark, let's say there will be 1000 POIs in total. (It's not necessary that the photos be stored in the same system.) This data will need to be maintained by a loose team of volunteers, so I'm looking for some kind of (free) platform that supports online editing. It would be great if it were possible to automatically layer a whole set of related geo objects, rather than having to specify them all individually. I'd also like to avoid having to run a server of my own just for this purpose. One possible solution I thought of was using Google Maps's "my places" (or whatever they call it now) to create the objects, then download them say to Dropbox as .kml's. That seems a bit clunky. What would be some better ways?


title: creating a st_contains query with geographies not geometries

text: I'm trying to find all the geography points that are contained in a geography polygon, let say all the traffic lights within a given city or neighborhoods, how can i convert this working st_intersect query to work with contained and remain with geographies not geometries               select count(*) from table1     where table1.UniqueID in (     SELECT table1.UniqueID     FROM table1 INNER JOIN table2 ON ST_Intersects(table1.point, table2.polygon)     where table2.UniqueID = 2);


title: arcpy.CopyFeatures between file GDBs?

text: I'm trying to figure out why I can't copy a featureclass from one file GDB to another. I can copy the featureclass to within the source GDB. There's no name conflict in the target GDB. I don't think I've set any GDB constraints that would stop me. I'm running 10.0 ArcGIS. Any thoughts? Here's what I've done to explore a bit and the error that I get. Thanks for any ideas!               [Dbg]>>> arcpy.env.workspace     u'C:\\Workspace\\wk_NHDPlus/NHDPlus03_Final.gdb/PRMS'     [Dbg]>>> arcpy.env.overwriteOutput     True     [Dbg]>>> arcpy.Exists('nsegment_flowlinesFL')     True     [Dbg]>>> arcpy.GetCount_management('nsegment_flowlinesFL')     <Result '5530'>     [Dbg]>>> desc = arcpy.Describe('nsegment_flowlinesFL')     [Dbg]>>> desc.dataElement.catalogPath     u'C:\\Workspace\\wk_NHDPlus\\NHDPlus03_wk.gdb\\Output\\nsegment_flowlines'     [Dbg]>>> arcpy.CopyFeatures_management('nsegment_flowlinesFL', 'C:\\Workspace\\wk_NHDPlus\\NHDPlus03_wk.gdb\\Output\\junk')     <Result 'C:\\Workspace\\wk_NHDPlus\\NHDPlus03_wk.gdb\\Output\\junk'>     [Dbg]>>> arcpy.CopyFeatures_management('nsegment_flowlinesFL', 'C:\\Workspace\\wk_NHDPlus\\NHDPlus03_Final.gdb\\PRMS\\junk')     Traceback (most recent call last):       File "<interactive input>", line 1, in <module>       File "C:\ArcGIS\Desktop10.0\arcpy\arcpy\management.py", line 1943, in CopyFeatures         raise e     ExecuteError: ERROR 000210: Cannot create output C:\Workspace\wk_NHDPlus\NHDPlus03_Final.gdb\PRMS\junk     Failed to execute (CopyFeatures).


title: How do I select by location based on selected features with Python?

text: Is there a way with Python to do Select by Location with the Use selected feature option? I'm trying to automate a process and am running into barriers. I need to be able to select the 2 lines within the same feature class that intersect a point, but I need to do this 1 at a time so I can capture the field attributes for each line. The 2 lines will then be dissolved and eventually appended back into the original line feature. Basically, I'm trying to merge the lines and keep the attributes. The exact tool I need is the Merge tool in the Editor drop down, but apparently that isn't accessible with Python for some dumb reason. There's also an Unsplit tool but of course it's only for Advanced, which I don't have. What I'm trying to do: From the Points feature class, get a list of the OBJECTID values and then use Select by attributes to select each point one at a time. But then I need to be able to do Select by Location on the Line using only the Selected point feature. Then once the 2 lines are selected, I'll copy their field attributes and dissolve them and then append them back into the original line and add the values back. ![enter image description here](http://i.stack.imgur.com/lLMcO.png) I think I got it to work finally. The only thing I haven't figured out yet is how to delete the original lines that are being replaced when I append the newly dissolved lines. Because I think the only way I can use the Select by Location is when a layer file and not a feature class, so I'm not sure how to go about it yet.


title: Openlayers MultiPolygon holes

text: I have trawled over many pages about the openlayers handlerOptions for drawing polygon holes. I cant create a polygon hole in existing polygons using the following:               var draw = new OpenLayers.Control.DrawFeature(     farmpaddocks,     OpenLayers.Handler.Polygon,     {     handlerOptions: {holeModifier: "altKey", multi: true, drawingHole:true,   },     title: "Draw Feature",     displayClass: "olControlDrawFeaturePolygon",          featureAdded: setAttributes          }          );      All I get are tiny dots on the vertices....no polygons lines are drawn when pressing "alt" key. Any ideas??


title: Is there a way to control transparency for the basemap layer in ArcGIS for Android?

text: The options seem sparse, but maybe I'm missing something (I'm new to both ArcGIS for Android and Android Honeycomb). Is there a way to make the basemap transparent? Update: this question is on the Esri application available in the android market. ![Screenshot from ArcGIS Android 1.0.1, build87](http://i.stack.imgur.com/nZkEh.jpg)


title: How to update raster or feature layers being used by ArcGIS Server by a Model

text: I have a mxd file that contains a few raster and feature layers that was published by ArcGIS Server as a map service. Meanwhile, a model was created to update the contents of the raster and feature layers in the mxd file. The intention is that all the users across my organisation will be able to see changes every time the model runs. The problem that I have is the geoprocessing by the model failed due to the error that is the data "May be locked by another application." every time I ran it. Could anyone please help with offering an idea how to get around this issue? p.s. I also tried to use File Geodatabase to store the layers but ended up with the same error. Thanks! Cheers, Alex


title: How to use javascript in fusion tables to query

text:
I'm having a quite bit of difficulty creating some functions for my fusion tables. I have two tables here: Points Parcels I'm trying to create a number of functions: 1) I want to be able to select the points layer using a drop down menu to select by its ID number. 2) I want to create a function where when I pick my point, I can tell fusion tables to find me the parcels that are within 1km of the point. It would preferable if I can change the distances, so I can do searches for 5km for example. Here's my code:               <!DOCTYPE html>

           Zone: <input type="text" id="search-string-l0">           <input type="button" onclick="changeMapl0()" value="Search">

           SPP_ID:            <select id="search-string-l1" onchange="changeMapl1(this.value);">             <option value="">--Select--</option>             <option value="                    <META http-equiv="Content-Type" content="text/html">          <meta http-equiv="content-type" content="text/html; charset=UTF-8">                    <body style=      Thanks for your help!


title: How do I visualize a time and space sequence of events in ArcGIS?

text: I would like to visualize this scenario: I have 5 weather stations. I would like to present a sequence of rainfall events for them. For example, rain was starts at station 2 at 6:09 and then went to station 5 after 8 minutes and ... Please guide me on how I can visualize the aforementioned scenario as a snap shot via GIS. Suggestions with tools description are highly appreciated.


title: Accumulate value loop

text: I am new to Phyton and currently trying to get my first script to work. Could someone tell me what is wrong with this code:                    total = 0     def accumulate(increment):      global total      if total:       if (total >= 1103):         total = (total - 1103)        else:         total += increment      else:         total = increment      return total           This one (the example from ESRI) works:                    total = 0     def accumulate(increment):      global total      if total:       total += increment      else:       total = increment      return total           I am trying to run this code in Arcgis 10.1! Thanks so much!


title: QGIS & GeoMedia .mdb .gws compatibility?

text: i cannot open Geomedia file directly from QGIS 1.8.0! Gdal driver >= 1.9.0 reference is made to below link qgis-connecting-geomedia-mdb-warehouse anyone have any experience in using GeoMedia file in QGIS???


title: segment lines according to point using pgRouting?

text: I'm working on postgresql and I am trying to create table to generate graphs. I have a table that contains all the objects of a city (points) and another contains the city street. My goal now is to segment the street in small segments that begins and ends with two objects (points). I try with different PostgreSQL query but I can not create a table of links that restrained the source and destination of each seqments ST_pointN actually gives the first and last point in a linestring against me by I want to have all the object on a street with the order to which I may have small segments. I want to see is what pgRouting can help me segment the streets and if so how? thank you


title: What are the advantages either way between GML and ShapeFiles?

text: After a good start many months ago of getting into GIS through QGIS, I had to drop what I was doing and did not touch it for quite some time. I am now able to get back into it and have a need to load GML vector maps. It strikes me that they have the same properties as ESRI Shapefiles, but there must be a difference. Can someone please give a simple explanation of the difference and the advantages either way? I see that one can be converted into the other, so if there must be a reason to do so.


title: Use GeoExt to edit Features attributes on a Map

text: I'm new to geoext and the documentation doesn't really help me. I've got a map with points on it. These points have attributes like 'name'. I coded a map with geoExt that has a vector layer which is poplulated by features(points) from a server. I can add features to the vector layer (not to the server) and select them. Now I'd like to select the features and edit their attributes using a window showing up. Saving them and in the next step saving them to the server. I've been trying to do this for a few hours now, but I can't get used to ext js/geoExt "philosophy" and I can't find examples that help me do what I wan't. Could somebody give me advice how to proceed? Edit: i worked on the code and i'm now able to read the attributes of a feature an show it in a popup. i still like to know how to edit the values and save it to the server Here's my (new) existing code:               var app, items = []          var osmLayer = new OpenLayers.Layer.OSM();     var vectorLayer = new OpenLayers.Layer.Vector("Vector Layer");     var map = new OpenLayers.Map({         projection: new OpenLayers.Projection("EPSG:900913"),         displayProjection: new OpenLayers.Projection("EPSG:4326"),         // todo: benutzerdefinierbar         units: 'm',         controls: [             new OpenLayers.Control.ZoomPanel()         ]     });          Ext.onReady(function() {         app = new Ext.Viewport({             layout: "border",             items: items         });     });          // registering a event that calls popup creation     vectorLayer.events.on({         featureselected: function(e) {             createPopup(e.feature);         }     });          // create a popup showing the properties of a feature     var popup;     function createPopup(feature) {         // so kann man properties schreiben         // feature.attributes.placename = "testtestPlacename";         popup = new GeoExt.Popup({             title: "Features Properties",             location: feature,             width:200,             // steht unter dem grid             //html: "hier steht html text",             maximizable: false,             collapsible: false,             unpinnable: false,             items: //new Ext.grid.PropertyGrid({                 new Ext.grid.EditorGridPanel({                 xtype: "editorgrid",                 //title: "testTitle",                 source: feature.attributes,                 autoHeight: true             })         });         popup.show();     }          // populate vectorLayer with features from remote source     var store = new GeoExt.data.FeatureStore({         layer: vectorLayer,         proxy: new GeoExt.data.ProtocolProxy({             protocol: new OpenLayers.Protocol.HTTP({                 url: '...',                 format: new OpenLayers.Format.GeoJSON({})             })         }),         fields: [             {name: 'title', type: 'string'},             {name: 'description', type: 'string'}         ],         autoLoad: true     });          items.push(new GeoExt.MapPanel({         region: "center",         map: map,         title: 'xyz',         layers: [             osmLayer,             vectorLayer         ],         tbar: buildMapToolbar(map, vectorLayer)     }     ));


title: What Reliable Reverse Geocoding API Can I Use?

text: I am working in a company where our team requirement is to do reverse geocoding lookups. Below is our requirement. To give you more background, the 3 geolocation-related needs I have are as follows:   1. Lat/Long lookup: returning country, city, state, region, postal code   2. Address lookup: returning lat/long   3. Postal Code lookup: returning country, city, state, region, lat/long (I assume the lat/long of the postal code’s centroid is what is usually provided) Which provider can I look into which will fulfill my above needs? I am looking for something which is used worldwide and are more reliable in terms of accuracy and coverage.


title: Selected Literature on coordinate transformations?

text: I am looking to hone my skills in the **theory** and the mathematicas involving the projections systems, projection transformations ( from xy to uv, φλ and vice versa). Can you share the books/pdfs/online-references that you have read and found them useful?


title: Morphometric and hydrological analysis with Python and ArcGIS

text: How can I write a program using python, that calculates some morphometric and drainage indexes using numeric attributes in the attribute table with ArcGIS? The data I have in my table are area and perimeter of the basin, and channel order of the watercourses (number and km of each order). For example, I need to calculate the Relationship of bifurcation with this formula: RB=Nu/Nu+1 (Nu - considered channel order; Nu+1 - channel order immediately ) above the considered one). So, I want to create a program that would be able to get the right value from the attribute table (in ArcGIS) and run with the formula. How can I do that? Thanks


title: Help opening GIS files

text: I was given 8 files that supposedly contain GIS data. They appear to be a file split into parts and are named: "stuff.zip [1_8].dat" , "stuff.zip [2_8].dat" , etc.. I tried stripping the ".dat" but that didn't help to unzip the files. I know this is probably a zip problem rather than GIS problem but I was hoping someone might recognize how the heck someone would send GIS data in this form. Opening with a text editor shows the following type of data, any idea if this is an old ESRI file type that was given the wrong extension or something? Thanks for any ideas. > UEsDBBQAAgAIAJlIQTmhyPOiMJdSALDsUgAXAAAAU2V3ZXJfSW52ZW50b3J5X01hcC5wZGYAM4DM > fyVQREYtMS41DQoxIDAgb2JqDQo8PA0KL0xlbmd0aCAyIDAgUg0KL0ZpbHRlciAgL0ZsYXRlRGVj > b2RlDQo+Pg0Kc3RyZWFtDQp4nOy8WbOqSNcuej8j/A8LUBBQ6QUVe7ABGxCxp7FB7BGbc2LXxfnt * * * Edit: it looks like this is a problem with how the data was encoded for the email message? After searching around a bit, I'm not sure how to deal with this. Obviously this is in the completely wrong forum at this point but thanks for any help. Here is a header before the alphanumeric string begins: > Content-Type: application/octet-stream; name="stuff.zip" Content-Transfer- > Encoding: base64 Content-Disposition: attachment; filename="stuff.zip"


title: How to delete features/records in ArcGIS Server using Python

text: How can I use python to delete records from a table, or features from a featureclass, when they are stored in ArcGIS Server? using REST "delete features" in arcgis api for javascript shows how to do this in JavaScript - what is the Python syntax?


title: ArcGIS10 Geocode Parser Pretype Problems

text: If I leave any important information out, let me know and I will append what I can. I am running into this problem when using the geocoder to locate various addresses in Quebec with ArcMap 10 SP1. The type is **US Dual-Ranges** , and I am using the FSA (no LDU) for Zone/Zip. My Address Locator appears to recognize all of the French road types. Diacritics were removed to normalize the street names between locator and addresses, and apartment numbers were removed. The issue I am currently wrestling with is that many addresses that fail with absolutely no matches quite suddenly find several as soon as the prefix type is removed. The one we shall focus on is Ch or "Chemin." After receiving the first address match results, I opened the review window to attempt a few manual matches. Here is what I found regarding the pretypes. With the CH prefix present, there are no matches: ![enter image description here](http://i.imgur.com/8QTup.png) http://i.imgur.com/8QTup.png When I remove the CH prefix, however, the story changes: ![enter image description here](http://i.stack.imgur.com/djqZ1.png) http://i.imgur.com/DsZCu.png From what I understood of the Geocoding process, the address given should be parsed for its component parts and matched against the locators. How is it, then, that removing one of the most integral parts of the match (the prefix) allows it to match? Even worse, why does leaving it in prevent any matches? The prefix type is even present in the locator, as shown in the second image. What is going on?


title: Need to remove small areas, but only if they do not have adjoining boundary

text: I'm generalizing some lat/lon data in grass64 and as a precursor to the generalization, I am removing small areas from the data using 'v.clean'. The following command is working well:               $ v.clean input=foo output=foo_rmarea_100km tool=rmarea thresh=100000000 type=area --overwrite      However, the one problem I have is that it is removing a small area that has a shared border with the rest of the larger data set. Is it possible to specify that only areas which DO NOT have a shared border with another area AND which are smaller than the defined threshhold be removed? Thanks