How to move points to another polygon based on attribute?

I have a point layer and a grid cell/fishnet layer. Each point falls within a grid cell. Now I want to move some of the points to another (the nearest) grid cell based on an attribute that exists in the attribute tables of both layers. The attribute contains information on the neighbourhood. Some points fall within a grid cell where the value of this attribute is not equal for both layers. So, I want to move these points to the nearest grid cell that has the same value regarding this attribute. Does anyone have an idea? I use ArcGis 10.2 Advanced. Thanks for your help!


Change timestamps according to transportation used

I have a GPX file with about 500 numbered positions. Unfortunately the time stamps are not correct. I do however know the positions and times the transportation form was changed. For example:               Position 1 to 75:      starts at 07:00 stops at 07:45     transportation by foot at an average speed of 10 km/h          Postition 76 to 250:     starts at 07:45 stops at 10:00     transportation by train at an average speed of 180 km/h      How can I change the timestamps and use different increments for different forms of transportation? What tools could I use?


Can i export BUS routes from OpenStreetMap as shapefiles to work on ArcMap?

The OSM has all the routing for the buses in my city, even showing the ID of the bus if you zoom close enough. I know how to download the shapefiles with the roads, rails, etc., from Geofabrik, but the bus routing isn't included there. If I can't get access to this routes in ArcMap, I'll have to draw them one by one, which doesn't make that much sense since that work has already been done. Please, can someone help me?


Arithmetic on rasters with variable names in Map Algebra/Spatial Analyst?

I'm a Python (and StackExchange) novice trying to figure out how to perform calculations on temporary rasters with variable names. My attempt at a minimum working example is below, but I would eventually like to do this with thousands of larger grids, so I expect that I will need to keep everything temporary until the very end to increase efficiency. The code below fails on the last line. The error message is "RuntimeError: ERROR 000732: Input Raster: Dataset bin2 does not exist or is not supported". If I replace the last line with the permanent filenames:               tmpout = Minus("D:/tmp/tmp" + `year`,"D:/tmp/tmp" + `yearold`)      then it runs fine. I am running IDLE in ArcGIS 10.0 with Python 2.6.5 (r265:79096, Mar 19 2010, 21:48:26) [MSC v.1500 32 bit (Intel)] on win32.               import arcgisscripting     import arcpy     from arcpy import env     from arcpy.sa import *          arcpy.env.extent = arcpy.Extent(0,0,2,2)          arcpy.CreateRandomRaster_management("D:/tmp/","tmp1","","",1)     arcpy.CreateRandomRaster_management("D:/tmp/","tmp2","","",1)     bin1 = arcpy.Raster("D:/tmp/tmp1")     bin2 = arcpy.Raster("D:/tmp/tmp2")          for year in range(2, 3, 1) :         yearold = year - 1         binyear = "bin" + `year`         binyearold = "bin" + `yearold`         print "fails in next line"         tmpout = Minus(binyear,binyearold)      I _think_ the problem is that `bin1` and `bin2` are being interpreted as strings as opposed to pointer to temporary rasters, but I can't figure out how to make my intent understood. Replacing `binyear` with `Raster(binyear)`, for example, did not help. Any help would be appreciated.


Rotation by attribute in Geoserver SLD style

I have a table with field calculated by function:                  NEW."rotationAngle" := degrees(                              ST_Azimuth(                                 ST_StartPoint(NEW.the_geom)                               , ST_EndPoint(NEW.the_geom)                              )                           ) - 90;      Now i publish this table as Geoserver layer and set style with `Roataion` tag:                 <Rule>                 <TextSymbolizer>           <Geometry>             <ogc:Function name="startPoint">               <ogc:PropertyName>the_geom</ogc:PropertyName>             </ogc:Function>           </Geometry>           <Label>             <ogc:PropertyName>annotation</ogc:PropertyName>           </Label>           <LabelPlacement>            <PointPlacement>              <Rotation>               <ogc:PropertyName>rotationAngle</ogc:PropertyName>              </Rotation>            </PointPlacement>           </LabelPlacement>           <Fill>             <CssParameter name="fill">#FFFFFF</CssParameter>           </Fill>           <VendorOption name="spaceAround">10</VendorOption>           <VendorOption name="autoWrap">85</VendorOption>           <VendorOption name="group">yes</VendorOption>         </TextSymbolizer>         <LineSymbolizer>           <Stroke>             <CssParameter name="stroke">#000000</CssParameter>           </Stroke>         </LineSymbolizer>       </Rule>      And result: ![enter image description here](http://i.stack.imgur.com/A7HRm.jpg) How can you see label angle not equals line angle.   When i publish table as a SQL view:               SELECT gid, annotation,"fontSize", the_geom, degrees( ST_Azimuth( ST_StartPoint(the_geom), ST_EndPoint(the_geom) ) )-90 as angle from annotations      With same style i get result: ![enter image description here](http://i.stack.imgur.com/Lk7ZP.jpg) Label angle equals line angle.   Why it's hepping? What i do wrong in first case? **UPDATE** I checked angle values in both cases and they are same.   Table's geometry have a custom projection. When i publish it on Geoserver i reproject it in `EPSG:4326`. **Update** I changed function to:                  NEW."rotationAngle" := degrees(                              ST_Azimuth(                                 ST_Transform(ST_StartPoint(NEW.the_geom),4326)                               , ST_Transform(ST_EndPoint(NEW.the_geom),4326)                              )                           ) - 90;      But not help. **Update** I found interesting behivor. I setted Label text as `rotationAngle` property and open both layers in Layers Preview and results: View Layer: ![enter image description here](http://i.stack.imgur.com/SSeOh.jpg) Table Layer: ![enter image description here](http://i.stack.imgur.com/n4JQy.jpg) Table Layer in application (map projection EPSG:900913): ![enter image description here](http://i.stack.imgur.com/HaLyO.jpg) Just wft going on? **Solved** Of couse it was stupid. My application's map showed in EPSG:900913. So why i calculate rotationAngle in EPSG:4326? I dont know. This is correct function:                      NEW."rotationAngle" := degrees(                              ST_Azimuth(                                 ST_Transform(ST_StartPoint(NEW.the_geom),900913)                               , ST_Transform(ST_EndPoint(NEW.the_geom),900913)                              )                           ) - 90;


Use GDAL/Python to add rasters (tiff) to an ESRI File Geodatabase

I'm using GDAL python to create a bunch of tiff files from a netcdf file.               def modelout_to_geotif(netcdfile, outdir, VarName):         """ Given the path to the netcdf output,         produces geotiffs for each band of the given variable name.              Input::             netcdfile: string path to netcdf file             VarName: string variable name in netcdf file (e.g. 'GridPrec')              Output::             creates geotiffs for each 'band' (e.g. time step). Output is written             to current working directory.                   """         try:             nci = gdal.Open('NETCDF:{0}:{1}'.format(netcdfile, VarName))             ncd = Dataset(netcdfile)         except:             print "Could not open input file: {0}".format(netcdfile)             return              try:             geotransform = nci.GetGeoTransform()         except:             print "Could not get geotransform.\n"             return         try:             projection = nci.GetProjection()         except:             print "Could not get projection.\n"             return              numbands = nci.RasterCount         times = [o2d(t) for t in ncd.variables['time']]         print "Found {0} bands to process for {1}".format(numbands, VarName)         for i in range(numbands):             time = times[i]             band = nci.GetRasterBand(i+1)             raster = band.ReadAsArray()             y,x = raster.shape                  output_file = "{0}_{1}.tif".format(VarName, time.strftime('%Y%m%d_%H'))             subdir = os.path.join(outdir, VarName)             if not os.path.exists(subdir):                 os.mkdir(subdir)             output_file = os.path.join(subdir, output_file)             if os.path.exists(output_file):                 continue             # Create gtif             driver = gdal.GetDriverByName("GTiff")             #TODO: Make this more generic             dst_ds = driver.Create(output_file, x, y, 1, gdal.GDT_CFloat32 )                  # top left x, w-e pixel resolution, rotation, top left y, rotation, n-s pixel resolution             dst_ds.SetGeoTransform( geotransform )                  # set the reference info              dst_ds.SetProjection( projection )                  # write the band             dst_ds.GetRasterBand(1).WriteArray(raster)      Now I would like to add those to a mosaic catalog in an ESRI File Geodatabase. I've managed to get it working using arcpy, but it's really slow, and I'd prefer to use GDAL. Is there a way? Below is my arcpy script:               import os     try:         import arcpy         from arcpy import env     except:         raise ImportError("Unable to import arcpy - are you using ESRI Python?")          env = 'P:/Enki'               def add_to_mosaic(filegdb, raster_dir):                  print env         #Add Raster Dataset type Raster to FGDB Mosaic Dataset         #Calculate Cell Size Ranges and Build Boundary         #Build Overviews for Mosaic Dataset upon the 3rd level Raster Dataset pyramid         #Apply TIFF file filter         #Build Pyramids for the source datasets         varnames = ['GridPrec']#,                     #'GridSnowOut',                     #'GridRad',                     #'GridTemp',                     #'GridSoilOut',                     #'GridRunoff']              for VarName in varnames:             gdbname = filegdb              mdname = os.path.join(gdbname, VarName)             messages = 'Failed'             sr =  arcpy.SpatialReference()         #        sr.name = 'WGS_1984_UTM_Zone_32N'             sr.factoryCode = 32632             sr.create()             print sr              try:                 noband = 1                 res = None                 #TODO: get from input data                 print 'trying to create %s' % VarName                 res = arcpy.CreateMosaicDataset_management(gdbname, mdname, sr, noband)                  except:                 print res.getMessages()                                 rastype = "Raster Dataset"             inpath = os.path.join(raster_dir, VarName)             updatecs = "UPDATE_CELL_SIZES"             updatebnd = "UPDATE_BOUNDARY"             updateovr = "UPDATE_OVERVIEWS"             maxlevel = "2"             maxcs = "#"             maxdim = "#"             spatialref = "#"             inputdatafilter = "*200701*.tif"             subfolder = "NO_SUBFOLDERS"             duplicate = "EXCLUDE_DUPLICATES"             buildpy = "BUILD_PYRAMIDS"             calcstats = "CALCULATE_STATISTICS"             buildthumb = "NO_THUMBNAILS"             comments = "Add Raster Datasets"             forcesr = "#"                  res = arcpy.AddRastersToMosaicDataset_management(mdname,  rastype, inpath, updatecs, \                                                        updatebnd, updateovr, maxlevel, maxcs, \                                                        maxdim, spatialref, inputdatafilter, \                                                        subfolder, duplicate, buildpy, calcstats, \                                                        buildthumb, comments, forcesr)             print res.getMessages()


Algorithm behind custom tool of ArcGIS for Desktop called xyz2grid21?

I am a researcher from Romania and I have used `xyz2grid21 ArcGIS tool` (by Min-Lang Huang) in order to convert xyz files (LIDAR data) to ASCII and ESRI grid in order to obtain a DTM. I need to know what algorithm does the XYZ2GRID tool use in order to rasterize the data? Is it using a mean or minimum value or does it interpolate somehow?


Creating weight raster for flow accumulation in ArcGIS for Desktop?

I am studying flood flow using ArcGIS for Desktop. When creating flow accumulation in hydrology, there is a weight raster. However, İf you do not have it, the programme sets its value to 1 which is not a suitable approximation for the simulations. How can I create a weight raster in ArcGIS for Desktop?


geometryType property of Vector Layers in OpenLayers

## **ORIGINAL:** I have a couple questions that are focused around the geometryType property with vector layers. When you create a vector layer you can set the geometryType property to limit the geometry types the layer supports... > {String} geometryType allows you to limit the types of geometries this layer > supports. This should be set to something like “OpenLayers.Geometry.Point” > to limit types. -- OpenLayers API Docs but in the Vector.js file the code is looking for an object when checking the geometry                if (this.geometryType &&               !(feature.geometry instanceof this.geometryType)) {                 var throwStr = OpenLayers.i18n('componentShouldBe',                           {'geomType':this.geometryType.prototype.CLASS_NAME});                 throw throwStr;               }      _Is this a typo in the API Docs or am I doing something wrong here?_ On a related note if I set the geometryType of a VectorLayer to be **OpenLayers.Geometry.Polygon** and then try to edit geometry of features in that layer it doesn't want to work because it tries to insert **OpenLayers.Geometry.Point** geometries into the layer to represent the draggable vertices. _Is there a something built in to get around this or do I need to write something to handle this?_ Thank you for any help. ## **UPDATE:** This seems to work but I'll be running more tests. _Vector.js_                   // Override for geometryType.         var geometryType = options && options.geometryType !== undefined ? options.geometryType : this.geometryType;      _ModifyFeature.js_               // Where "addFeatures" exist I've tacked on the geometryType     this.layer.addFeatures([this.radiusHandle], {silent: true, geometryType: OpenLayers.Geometry.Point});      Thanks for the help!!


ogr SetAttributeFilter uppercase?

How do I use the ogr `SetAttributeFilter` with uppercases? Code example:               lyr = shp_obj.GetLayer()     sel_str = "%s = '%s'" % (field, value)     lyr.SetAttributeFilter(sel_str)      Can I just do:                sel_str = "UPPER(%s) = '%s'" % (field, value)      ?


cannot download 3rd party plugins in qgis

1. I want to download SDA4PP plug-in in QGIS. Unfortunately I cannot see the button to download 3rd party plugins, as indicated in the install instructions at that link. I am using this program on a company network. I have put in some proxy settings which mean that I can access the global/official repository. Before I did this I could not even access this. When I go to the plugins tab though there are only a few listed and when I put something in the search box nothing comes up, i.e. it only searches currently downloaded plugins (those that came with the install). Anything I'm missing?


All LGAs for Australia in KML or KMZ format

Can somebody please show me where to find a KML or KMZ file that includes the digital boundaries for all Local Government Areas in Australia? Thanks!


Openlayers Plugin, Project CRS and Grids in Composer

Apologies in advance for the potential for this to be a duplicate question, I've had a good look, and I can't see a solution thatI understand well enough to help me. I have my project CRS set to a local UTM projection, however when I add a map from Open Layers the Project CRS changes to EPSG:3857. When I add a grid in composer, I need to display the longs and lats of the local UTM, not the Open Layer EPSG. I understand that I need to have the Project CRS the same as the Open Layers CRS, but is there a way I can 1) save and reproject the section of the Open Layer raster into my local UTM, so it still lines up with my vector layers or 2) display the appropriate coordinates on the grid, despite the Project CRS being in EPSG:3857. This is an issue that I have with many of my maps, where I use aerial imagery for my base layers, but need to display local coordinates in my final map. Will there be an easy solution to this in QGIS 2.0?


GeoExplorer: "Could not create layer record for layer group"

I installed GeoServer via OpenGeo Suite, and still working on it. I created a "Layer Group" and add layers from local PostGIS. I tried to preview it from "Layer Preview" and show perfectly through the "OpenLayers" option, but unable to show in "GeoExplorer". I got the javascript warning message says: `"Could not create layer record for layer 'testgroup'. Check if the layer is found in the WMS GetCapabilities response."` on GeoExplorer.js:1441 To view the source code shows like:               1441: {name:"fixed",type:"boolean"},{name:"selected",type:"boolean"},{name:"restUrl",type:"string"},{name:"infoFormat",type:"string"},{name:"getFeatureInfo"}];c.fields.each(function(a){i.push(a)});b=new (GeoExt.data.LayerRecord.create(i))(h,b.id);b.json=a}else window.console&&0<this.store.getCount()&&void 0!==a.name&&console.warn("Could not create layer record for layer '"+a.name+"'. Check if the layer is found in the WMS GetCapabilities response.");return b},getProjection:function(a){var b=this.getMapProjection(),         "Could not create layer record for layer 'testgroup'. Check if the layer is found in the WMS GetCapabilities response."     1442: c=b,a=a.get("srs");if(!a[b.getCode()]){var c=null,d,e;for(e in a)if((d=new OpenLayers.Projection(e)).equals(b)){c=d;break}}return c},initDescribeLayerStore:function(){var a=this.store.reader.raw;this.lazy&&(a={capability:{request:{describelayer:{href:this.url}}},version:this.version||"1.1.1"});var b=a.capability.request.describelayer;if(b)a=a.version,1.1<parseFloat(a)&&(a="1.1.1"),a={SERVICE:"WMS",VERSION:a,REQUEST:"DescribeLayer"},this.describeLayerStore=new GeoExt.data.WMSDescribeLayerStore({url:this.trimUrl(b.href,      I am sure that I can preview them separately in GeoExplore from "Layer Preview" page. It seems everything goes well when I checked it from `http://localhost:8080/geoserver/wms?SERVICE=WMS&REQUEST=GetCapabilities`, without any error massages or exceptions shown, and the layer group description show at the bottom of the page. Does anyone run into the same problem? **//UPDATE** I found that I missed a javascript error massage when I view the `REQUEST=GetCapabilities` page. It says: `"Uncaught TypeError: Cannot set property 'boxSizing' of null"` on ContentVeil.js:11 To view the source code shows like:               11: veil.style.boxSizing = "border-box";         "Uncaught TypeError: Cannot set property 'boxSizing' of null"     12: veil.style.borderStyle = "solid";      I am wondering where the error comes from. Furthermore, I cannot load the **"field list"** when editting limit condition in style's rules in GeoExplorer.


QGIS 2.4.0 How to display only selected features on map?

Not sure if this option ever existed, but the new QGIS version hides a lot of things somewhere else and this seems so simple that I'm kind of convinced it was possible to do this before. I have quite a large data set with ship tracks. I now want to analyse individual ship tracks (made from point layer with points2one) - however, it is so crowded that it is hard to look at individual tracks. I don't want to save each single track as a new layer, so I wonder if it is possible to select one track in the attribute table and display only that one on the map/hide the others? Thanks a lot!


Combining multiple year yield data in ArcGIS for Desktop?

I am using ArcGIS 10.1 for Desktop and trying to find the simplest way to combine multiple year of crops yield data into one file and get the average yield. The files are point files from a yield monitor.


How to split polygons at intersections in ArcGIS 9.3?

I have two polygon layers. I'd like to split the features in one based on where they intersect the other layer. I've tried using the Clip tool, but there is no option to simply split the features -- you're either left with all of the overlapping parts or non-overlapping parts. I've found several Cut Features tools in some obscure extensions in 10.1, but unfortunately I don't have access to those. Can anyone advise on how to do this in out-of-the-box 9.3.1? **UPDATE** I'm using an ArcInfo license. Here is a visual of what I'm trying to achieve: ![enter image description here](http://i.stack.imgur.com/X2vEe.png)


Import error:No module named _gdal

I have installed python 2.5.1 and FWTools2.4.7 in my laptop(Windows 7) when i run my pythonscript in FWTools i get following error Import error:No module named _gdal can any one suggest how to solve this problem?


ogr2ogr does NOT save "styles", say "color", information; how do I get around it?

I'm using `ogr2ogr` to convert KML to GeoJSON and I'd like to extract style information along with the data. The bash script I am using is:               filekml="/tmp/poly.kml"     filegeojson="/tmp/poly.geojson"          ogr2ogrProg="/usr/local/fwtools/FWTools-linux-i686-3.1.0/bin_safe/ogr2ogr"          # randomNum=14089918888     randomNum=`date '+%s'`          # Grab the National Weather Service Polygon (KML) data file, using SH script in Linux:     getFromURL="http://www.srh.noaa.gov/data/radar/poly.kml?${randomNum}"     wget ${getFromURL}     echo y | mv poly.kml\?${randomNum}  ${filekml}          #-- use OGR program above to convert the downloaded KML file to geoJSON format     rm -f ${filegeojson}     $ogr2ogrProg -f "GeoJSON"  ${filegeojson}  ${filekml} -nln sr_polygon          # ---- examine the KML file and its' sister geoJSON file     vi /tmp/poly.kml /tmp/poly.geojson      You will notice KML contains lines like these: ==> color ffa9a9a9   ==> scale 0.5   ==> href http://www.srh.noaa.gov/ridge/graphics/FFW.gif   BUT--- geo JSON file has absolutely NO information on STYLE (that was embedded in KML).


searching data from wms geoserver in jquery mobile

hye all, need someone help, i use template from openlayer using jquery mobile, im really new in openlayer, i want to search data in wms, the coding is below $('#searchpage').live('pageshow',function(event, ui){ $('#query').bind('change', function(e){ $('#search_results').empty(); if ($('#query')[0].value === '') { return; } $.mobile.showPageLoadingMsg();                   // Prevent form send         e.preventDefault();              var searchUrl = 'http://ws.geonames.org/searchJSON?featureClass=P&maxRows=10';         searchUrl += '&name_startsWith=' + $('#query')[0].value;         $.getJSON(searchUrl, function(data) {             $.each(data.geonames, function() {                 var place = this;                 $('<li>')                     .hide()                     .append($('<h2 />', {                         text: place.name                     }))                     .append($('<p />', {                         html: '<b>' + place.countryName + '</b> ' + place.fcodeName                     }))                     .appendTo('#search_results')                     .click(function() {                         $.mobile.changePage('#mappage');                         var lonlat = new OpenLayers.LonLat(place.lng, place.lat);                         map.setCenter(lonlat.transform(gg, sm), 10);                     })                     .show();             });             $('#search_results').listview('refresh');             $.mobile.hidePageLoadingMsg();         });     });     // only listen to the first event triggered     $('#searchpage').die('pageshow', arguments.callee);      }); really need anyone help Thanks


What is the process to calculate Helmert parameters?

I have some double points (the same points in 2 different datums). I'd like to calculate the Helmert parameters but I don't understand how. NOTE: I don't need software. I need an algorithm in order to write a small Delphi application.


How do I export symbology to ArcGIS Online?

We're planning to put up a bunch of our maps into ArcGIS Online and we'd like to create a streamlined process to do so, with as much automation (Python) as possible. This is our first test with a map of a city park with lots of unique symbols. We figure we'd get a map up there, then figure out what can be scripted. Anyhow, so, we have an .mxd of a park where pretty much each layer has fancy symbols for things like parkbenches, bathrooms and the like. Upon fiddling around with ArcGIS Online, it looks like it only takes .txt, .csv, and .gpx files. I managed to get layers into "My Content" but it doesn't let me add it into a map. When I go "Add>Search for Layers" and tell it to search in My Content, nothing shows up. I then shared the layer in My Content so it was public on ArcGIS Online, and tried seeing if the search func. would see the layer in ArcGIS Online. No luck. How do I get it to take attribute data, but also the symbology data from my layers? I hope that made sense. Will definitely clarify as much as needed.


How to use ST_Intersection?

Here's a quick summary about what I'm trying to do: I have 3 tables in Postgres, 'a' and 'b', each have a Polygon column, and 'c' has a Point column. What I'm trying to do here is to get the geometries intersections between 'a', 'b' and 'c', and to display such geometries on an OpenLayers vector layer. I already know how to display any kind of geometry from a String in OpenLayers, but I'm having troubles with the PostGIS' ST_Intersection function, I'm doing this:               SELECT ST_Intersection(a.geom, b.geom) as inter from a, b;      where a.geom and b.geom are both the geometry columns, and I get this error message:               NOTICE:  TopologyException: found non-noded intersection between 515172 2.14408e+06, 497067 2.13373e+06 and 501321 2.13546e+06, 471202 2.14843e+06 500621 2.13576e+06      ERROR:  GEOS Intersection() threw an error!      Also I tried to express the resultant geometry as text using ST_AsText like this:               SELECT ST_AsText(ST_Intersection(a.geom, b.geom)) as inter from a, b;      but it send me this error message:               HINT: No function matches the given name and argument types. You might need to add explicit type casts.      I don't know what I'm doing wrong, I just want to get the Polygons' WKT to display it on OpenLayers, here's how I display a geometry from a WKT:                                   var in_options = {                             'internalProjection': new OpenLayers.Projection("EPSG:4326"),                             'externalProjection': new OpenLayers.Projection("EPSG:4326")                         };                               var fea= new OpenLayers.Format.WKT(in_options).read(data); //data is the string with the WKT                         vectorLayer.addFeatures([fea]); //this piece of code works great                         map.zoomToExtent(bounds);      UPDATE: I tried the next:               SELECT ST_Intersection(a.geom, b.geom) as intersect_ab FROM a INNER JOIN b ON      ST_Intersection(a,b) WHERE ST_Overlaps(a.geom, b.geom)      AND ST_isvalid(a.geom)='t' AND ST_isvalid(b.geom)='t';      but I get the next error message:               ERROR: Function st_intersection(a,b) does not exist.     HINT: No function matches the given name and argument types. You might need to add explicit type casts.      I added the isvalid to verify only valid polygons are being evaluated, but it's telling the error is in the ST_Intersection(a,b), both a, b and c have the same SRID so I'm really confused, sorry if I'm asking too much, but I'm quite new with PostGIS so I hope I'm not bothering you a lot. Thanks.


Access Metadata from ArcSDE-tables using python

Is it possible to access the metadata of a feature class in a ArcSDE-data, using arcpy / python (on Arcgis 10.0). I suppose I could export the metadata to xml and parse the xml, but is there a more direct approach? I need to make a summary table of all the data in our database, with several elements from the metadata (author, last update etc. ). We always enter metadata for each new dataset in our DB.


Map tiles and max/min latitude and longitude

I'm trying to build a map based on a center point lat/long, the map will always be a fixed size say 1024 by 1024. The user can chance the center point and zoom. If I populate the map with map tiles from map quest for example, is there someway I can know the max/min latitudes of the map currently being displayed? I want to know this so I can accurately project things onto the map. Thanks


How to extend Geoserver raster data sources

I am new to web GIS world. I have installed Opengeo Suite on my PC running with Windows 7 operating system. So, geoserver is locally installed (http:localhost:8080/geoserver). I have tried to import shapefiles into geoserver. Now I want to import .ECW and .JPEG image into Geoserver. When I click _stores_ and _add stores_ , I can see the following list of rasters,![enter image description here](http://i.stack.imgur.com/Mm70M.png) I looked into the following path and found that image extraction files are installed there, C:\Program File\OpenGeo\OpenGeo Suite\webapps\geoserver\WEB-INF\lib ![enter image description here](http://i.stack.imgur.com/lniGx.png) From internet search, I found that some people have advised to install gdal- data in the above mentioned path. I extracted gdal-data into the _lib_ directory, ![enter image description here](http://i.stack.imgur.com/oUVHc.png) Now I restart OpeGeo dashboaard and geoserver, still I cant see .ECW or :JPEG foramt in the raster data sources. Can anyone help me out?


How to convert shape files (.shp) to Cad (.dwg ) format using Quantum Gis 1.7.0?

How to convert shape files to .dwg format using Quantum GIS 1.7.0 ? Thank you


Request to run mapbasic 10.5 code in mapbasic 9.1 version

I have created mapbasic codes using 10.5 version. Now I have to run in a lower version (9.1) Please anyone help me how to do this.. How to copy paste the program codes from 10.5 version to 9.1 version? Thanks, Jagan


Set feature attribute by name via QGIS python api

I try to set feature attributes by attribute name with this code               pr = layer.dataProvider()     pr.addAttributes([QgsField("test", QVariant.Int)])     layer.updateFields()     for feature in layer.getFeatures():      attrName = 'test'      feature[attrName] = 1      but in result I have NULL in all fields. If I use               feature.setAttributes([1])      it works properly. What I do wrong? Why feature[attrName] = 1 doesn't work? **Update** : Find this solution               pr = layer.dataProvider()     pr.addAttributes([QgsField("test", QVariant.Int)])     layer.updateFields()     for feature in layer.getFeatures():      attrName = 'test'      feature[attrName] = 1      pr.changeAttributeValues({feature.id() : {pr.fieldNameMap()[attrName] : 1}})


Stray horizontal lines after converting from GeoJSON to topojson

After converting a GeoJSON file to topojson, and loading the topojson into Leaflet, stray lines appear in the resulting map. Like this: ![enter image description here](http://i.stack.imgur.com/QBjcy.jpg) What is causing this? If I view the GeoJSON, these stray lines do not appear. I've played with the quantization behavior, and this does not seem to help. thanks. (I am using Leaflet-Omnivore to digest the topojson.)


Showing data in form view (instead of grid view) using ArcGIS API for JavaScript?

Can anyone suggest how to bind data with the help of query in a Form view of the ArcGIS API for JavaScript? I have done it in Grid view but as the requirement of a client is that I have to bind it to a Form view.


How to create animated timeline in ArcGIS for Desktop?

I can't find any options in ArcGIS timeslider to create a timeline like presented in this video. The video itself is a guide for creating time animation, but no information is given how to create such timeline like in the end. I have tried to search some information about that and can't find anything usefull. Does ArcGIS support such functionality? **UPDATE 1:** I need an autorefreshed timeline (which is presented not as usual date text). It should be animated line where each date is a vertical line on a timeline. Also I need to label each date on it as a text. See picture below (frame by frame animation). ![enter image description here](http://i.stack.imgur.com/qEMpX.gif) Added a post at ArcGIS Idea, to add such functionality. Upvotes are welcome.


QGIS: drawing order of lines and polygons (line between polygon)

I have 2 layers. One with forests and lakes (polygon) and one with rivers (line). Now I want the features to draw in that order: The river have to be UNDER the LAKES and the forest have to be UNDER the river and OVER the lakes. Tricky? Because I have to export as vector-pdf, I can not use the feature blending modes. (or is there a hidden way to export blending modes to vector? I need the borderlines of the lake and the river. I can not physically change the layers. ![enter image description here](http://i.stack.imgur.com/UA8Aa.png) ...or is there a other way to draw the river flow into the lake?


Inserting full bleed images in browser background instead of maps using CartoDB?

I am building a product that will utilize high fidelity photography. I would like to show my photos in a slideshow format that takes up the entire area of the browser. CartoDB and the Odysssey.js would work well. Can I substitute the maps for photography in the background of my browser?


Grass7 - r.regression.multi and new raster

I am on the learning curve with Grass7. I managed to make r.regression.multi work, but there is one thing that would be of high interest for me: it is possible to create a new raster showing the results (residuals) of the analysis in a new raster? This would enable me to nicely locate the areas with high residuals.


How to add the current document path on Print Composer?

I'm trying to demonstrate to my company that we can begin using QGIS in conjunction with (and sometimes even instead of) ArcGIS and MapInfo. I've started setting up some layout templates but have run into a. In Print Composer, is there a way to add the current project path? In ArcGIS, the dynamic text is:               <dyn type="document" property="path"/>      I'd prefer not to have to copy/paste the document path every time I make a map. * * * **EDIT:** _To address answer and comments by artwork21 and Nathan W:_ I'm using QGIS 2.0.1. I would also like this path to appear to anyone who may need to open the project for editing without modifying their `startup.py` file.


Automated quality control and assurance - geodata?

I've recently started looking at the issue of automated quality control and assurance of geographic data. What software exists that can be used for accomplishing this? Esri has got their product 'Data Reviewer' which seems interesting, are there any comparable alternatives to this product? I would like to tag this question with quality-control and quality-assurance, but can't do so.


Don't route through bollards, concrete blocks, etc with osm2po

Is it possible to prevent road segments being used when these are blocked by bollards, concrete blocks or the like? I have one example where there the road Solhøjvej is blocked with a so called concrete block as you can see here http://www.openstreetmap.org/?lat=57.010248&lon=9.936238&zoom=18&layers=M Here is a photo of the "block": https://maps.google.com/?ll=57.010343,9.934331&spn=0.000929,0.00284&t=m&layer=c&cbll=57.010252,9.934346&panoid=SrWxGTFwoDbDZeDgXHl8vg&cbp=12,10.73,,0,10.62&z=19 I’m exporting the Openstreetmap to PostgreSQL using osm2po and performs routing, but this does not know of these concrete blocks or bollards or stuff blocking the road. I’m aware that this is perhaps not a simple task, to detect and split road segments when such obstacles occur, but that would be very useful. Am I just missing an option?


Developing Carto Mobile App

I have the project to develop a mobile application (Android) being able to get offline maps, go on the field, collect datas about the area we're at and then synchronize it with our data server when we're back to the office. This is for my end of study work (sry, don't know a best translation in english), and i looked at apps like gvSIG Mini, qGIS for Android, qmap, .. But none of this really fits my need, as it is hardly adaptable .. Any ideas on api/existing apps/anything ?


How long for surveying quality GPS to get accurate 3D position information

I recently read an Ordnance Survey Blog article that seemed to suggest that it was necessary to leave the GPS device in use for four hours to gather requisite data. I am aware that "domestic" GPS units benefit from being left in one position to get better positional accuracy but had assumed that survey quality equipment was "accurate" immediately. What am I missing please? Nige


Using FME to merge lines in multiple SHP layers

I have multiple layers of complex laine data that I need to flatten int a single layer of lines. The issue is that some of the line data overlaps and intersects and is not end node to end node. How in FME can I take this data anf flatten it into a single set of joined lines?


How to use a database in Qgis

What is the primary use of database in Q-Gis? How to connect to a database in Qgis?


Unable to create a circular polygon with ArcGIS Explorer SDK

I am having trouble creating a circle with a specified radius of 200 feet. I have tried several ways to translate the points along the circular arc. here is my current code:               var display = ESRI.ArcGISExplorer.Application.Application.ActiveMapDisplay;          var center = display.TrackPoint();     var convertedRadius = Unit.Convert(200, Unit.Feet, center.CoordinateSystem.Unit);         var firstPoint = GeometryOperations.Move(center, convertedRadius, 0);          var points = new List<Points>();     points.Add(firstPoint);          for (int i = 1; i < 360; i++)     {         var point = GeometryOperations.Rotate(firstPoint, center, (Math.PI / 180) * i) as Point;         points.Add(point);     }          var poly = new Polygon(points, display.CurrentCoordinateSystem);     var circleGraphic = new Graphic(poly, Symbol.Fill.Solid.Red);     display.Graphics.Add(circleGraphic);      I only seem to have a problem creating a polygon. Using the above loop I can place 360 markers in a circle around the center point. Alternately, If there is a way to create a buffer from a point using programmatic methods that would suffice. I am trying to query all shapes that are within 200 feet of a target point.


Add print/export button to axToolbarControl

I have a map that's being rendered inside a custom application, and now I need to add to the ESRI.ArcGIS.AxToolbarControl an export (to PNG) button and a print button, but honestly don't know how. I've already looked at this and this samples, but couldn't make them work.


Where is DatabaseSupport directory in the ArcGIS client?

May I know where is this folder stated in the guide? DatabaseSupport directory in the ArcGIS client installation location and move it to a location on the Oracle server. http://resources.arcgis.com/en/help/main/10.1/index.html#//019v0000000s000000


convert shapefiles to cad files

hi i'm having a difficulty in converting shapefiles to cad file. i'm currently using the ArcGIS 9.3.1 version and i tried to convert it using the toolbox>conversion tool>to cad>export to cad....the conversion was partly successful but then it doesn't have any attributes when i tried to open it in ArcMap. When i try to open it in Autocad i can't see a thing... how can i resolve this? thank you.


How convert a genuine UK shapefile(other than the one from natural earth) to d3/topojson?

I am following Mike's page (http://bost.ocks.org/mike/map/) to make UK's map. I've converted Mr. Crone's UK postcode shapefile into topojson but failed to display it correctly. http://www.sharegeo.ac.uk/handle/10672/51?show=full               ogr2ogr -f GeoJSON postcodes.json gb_pca_2010.shp     topojson -o uk.json --id-property Label postcodes.json      The resulting json is of 2.2M big. I have used the following code to display it               var height, svg, width;     width = 960;     height = 1160;          svg = d3.select("body").append("svg").attr("width", width).attr("height", height);          d3.json("uk.json", function(error, uk) {       var projection, subunits;       if (error) {         return console.error(error);       }       subunits = topojson.feature(uk, uk.objects.postcodes);       projection = d3.geo.mercator().translate([width / 2, height / 2]);       return svg.append("path").datum(subunits).attr("d", d3.geo.path().projection(projection));     });      The result is that firefox freezes and in the end, I get half black and half white display. :( I have also tried mapshaper to simplify it to 1MB but still not working. What else I could have tried?


10.1 beta - register dll-s

I just tried to register my extension for ArcGIS 10 in the 10.1 beta and it failed when I tried to register the dll: "Could not load file or assembly 'ESRI.ArcGIS.ADF.Local' Version=10.0.0.0." Do I have to re-compile it against 10.1 dll-s in order to make it work? For the 9.x version all I had to do was to compile them against the base level (9.2) and it worked in 9.3 as well.


Examples of geovisualizations of global connectivity

I stumbled upon interesting visualization of global connectivity patterns: > 'Major road and rail networks in Africa, along with transmission line and > underwater cable data.' ![enter image description here](http://i.stack.imgur.com/eQEPh.jpg) (Source, more pictures in the gallery. Original.) Can you point me to other examples of visual representations of connectivity in global scale? Connectivity in this case applies to both transportation and communication. And geovisualizations might include static maps as well as interactive web applications.


ArcObjects: How to Speedup Spatial Query to Feature Layer?

I have ROAD (POLYLINE) Shape File with more than 400,000 Records. I am working on a code which goes like this:        1. Loop through all Feature of the Shape File   2. Take First Feature   3. Find the LAST VERTEX of that feature. I get a POINT   4. Query to FEATURE Layer all the POLYLINES Connected to that POINT and get back result in feature cusros   5. Do some Process on returned feature of RESULT Cursor   6. Go to Step1 until Loop is not completed. When I do STEP-4, it obviously take lot of time because, I create a BUFFER of Say 1 meter on that POINT and Make a spatial FILTER and Run LAYER.SEARCH QUERY on entire featurelayer using that SPATIAL FILTER. Everytime querying 400,000 records making it very slow. Is there a better way to do this thing ?   * I have to use ArcObjects VBA.    * I do not want to Create NETWORK from the layer.


Preparing CSV files for use in ArcGIS for Desktop?

I wonder if there are any recomendations how to properly prepare CSV files for use in ArcGIS for Desktop. I ask because I have some troubles using CSV files because ArcGIS attributes wrong field types to my columns and also misinterprets special characters such as á or ê. I have read in the ESRI forum that there is a so-called schema.ini file that defines somehow the field types e.g "Col22=V002 Text" see here http://forums.esri.com/Thread.asp?c=93&f=1149&t=64464 That's kind of funny because I have often seen these .ini files on my disc but never actually wondered what they are good for. It is kind of weird that Excel stores such metadata in an extra file since other programs like R don't do so. I already tried to manipulate this .ini file with little success since I didn't find out how to apply for example "string" type. There are some information on MS sites, see here: http://msdn.microsoft.com/en- us/library/windows/desktop/ms709353%28v=vs.85%29.aspx but I couldn't find a solution. Also I didn't really like the idea to work with this .ini file because it is quite a bunch of work to define and type all the fieldnames when I have let's say 50 columns. And these .ini files might get lost, etc. **So, is there any better way to properly prepare CSV files for use in ArcGIS for Desktop?**


repeating records in attribute table

I have a point shapefile containing records of bird sightings with location, species and land cover type that has about 30,000 records. I noticed today that one record is repeated about 60 times, with the same gps points, same location name, same everything. But when I highlight all of the repeating records, they are actually all different records. I can't seem to get the attribute table to accurately reflect the records; I've tried copying the file, renaming it, exporting the incorrect records to a new shapefile, restarting Arcmap multiple times, etc. Is the file corrupted? What can I do to fix these records?


Leaflet with EPSG4326 Tile Layer (using maptiler.org)

I have successfully been using Klokan Petr Pridal's MapTiler (http://maptiler.org ) to take GeoTiff files and create TMS tile pyramids for Leaflet (option #1 which projects to EPSG3857, also called Google Maps compatible Spherical Mercator). But my source data is in EPSG 4326, so I am now using option #3 for WGS84 Plate Caree (Geodetic). I believe it is using gdal2tiles under the hood. The issue I have is that I do not see any of these overlay tiles in Leaflet. Even though am being careful to make sure the map is created with the right CRS. I have tried tuning TMS: false, and TMS: true. No luck. One thing I see is that the numbering even for the top levels z/x/y is weird: (e.g. 0/1/2.png ). Should a Krakow, Poland map really be at coordinates: 5/35/24? or 10/1136/795? With the Spherical Mercator I have 10/568/346!!! Also, all the specs for TMS tile numbering that I have seen are only good for EPSG3857 (spherical Mercator) which I believe has bounds of _85/-85 lat. But EPSG4326 would go to 90. So what gives in the numbering?               var map;     var krakowLayer = null;     var mapOptions = {         center: new L.LatLng(0, 0),         zoom: 1,         tms: true,         crs: L.CRS.EPSG4326     };         map = new L.Map('map', mapOptions);         krakowLayer = new L.TileLayer("./tiles4326/{z}/{x}/{y}.png", { tms: false, opacity: 0.8 });          function addLayers() {         var baseLayers = {             "CloudMade": cloudLayer,         };         var overlays = {             "Krakow Layer": krakowLayer         };         L.control.layers(baseLayers, overlays).addTo(map);     }


Using osm2pgsql data with pgrouting 2.0

I have successfully manage to make the workshop from of pgrouting (wrote for v1) with pgrouting v2. Now I want to use pgrouting 2 on my postgis database created with osm2pgsql. The problem is first to create the segments and this blog post was very helpfull for this. However,I get:               gis=# select pgr_createTopology('network', 0.01, 'the_geom', 'gid');     ERROR:  We can not determine the srid of "network.the_geom"!      Is there a way to set the srid properly so that `pgr_createTopology` finds it. Or is there an alternative way? I can't believe osm2pgrouting is the only importer for pgrouting. Thanks


A couple GeoExplorer questions

I just started working with the OpenGeo Suite. I'd really like to use it for a County web map, but have a couple questions... 1) Is there a way to add a lat/lon or other X/Y location to the URL and have the map center and zoom to the point? Something like "?lat=46.265&lon=-90.206&zoom=7" at the end of the URL. 2) Is the Query tool available in a published map? It's not in the Choose Tools window when publishing. Thanks. === UPDATE === I will drop part 2 and stick with part 1. I guess I forgot to mention the crucial bit that I'm working in GeoExplorer. I did have the tag, but don't fault anybody for missing it. The bounding box isn't working in GeoExplorer. I have the map published in an iframe. Here's a screenshot of the test page without the bbox... i.imgur.com/Wc022Ee.jpg. Here's a screenshot with the bbox... i.imgur.com/7MES4cn.jpg. That's just a blank map that hides behind the real map and only shows up when the real map doesn't. I do have a correct bbox... i.imgur.com/MbAOgVD.jpg in the code. (I guess I need a higher rep before I can add more direct links, so you'll need to copy paste the image URLs if you want to see them.) Here's the source behind the page...
Iron County Test Map

This is a test map, version T-001.  A disclaimer might go here.
                  I tried putting the "?bbox=" at the end of the "src=.../1". Thanks for any help. And thanks for the welcome.


How can I convert an Integer field to a Real field in a shapefile?

How can I convert an Integer field to a Real field in a shapefile ? I'm searching a tool, or a way to do that (may be under QGis ?!). Regards


Server side clustering - cluster (feature) geometry

I am trying to implement server side clustering and show it using `OpenLayers`. I am using this good material as an example. As a clustering result, a get some list of features, some of them are clusters. Like this (it is an output of previous example):               Array (          [0] => Array (              [0] => Array ( [id] => marker_3 [lat] => 59.431602 [lon] => 24.757563 )              [1] => Array ( [id] => marker_4 [lat] => 59.437843 [lon] => 24.765759 )              [2] => Array ( [id] => marker_6 [lat] => 59.434776 [lon] => 24.756681 )          )          [1] => Array ( [id] => marker_5 [lat] => 59.439644 [lon] => 24.779041 )          [2] => Array ( [id] => marker_2 [lat] => 59.432365 [lon] => 24.742992 )          [3] => Array ( [id] => marker_1 [lat] => 59.441193 [lon] => 24.729494 )      )      where 0 element is cluster containing 3 features. But how can I compute this clusters geometry? Or how can I visualize it using `OpenLayers`?


getting the specified geographic transformation, or at least the GCS, of a data frame

First, the TLDR digest:   * I'm getting a **bogus value** from `dataFrame.geographicTransformations` on some machines; how do I get the correct value as shown in the properties dialog? and/ **or**   * Given a spatial reference with a PCS, how do I find out (in Python) **the GCS** or datum that it's based on? And now the details... (Reposted from http://forums.arcgis.com/threads/76061-getting-the-specified- geographic-transformation-or-at-least-the-GCS-of-a-data-frame) I'm writing an arcpy script in ArcGIS Desktop 10.0. This script opens a user- specified MXD and processes user-specified layers in it. Part of the processing involves projecting the data from whatever GCS it's in, to whatever PCS the data frame has in its spatial reference. Sometimes, the input GCS and the output PCS are in the same datum (e.g. WGS_1984), and so I don't need to supply a geographic transformation to the `Project_management()` tool. Other times, they are not in the same datum; e.g. the data frame coordinate system might be Europe_Lambert_Conformal_Conic, which uses the datum ED_1950. If I don't supply the right geographic transformation (or lack thereof), the Project tool fails. How do I know what transformation to use? In the case where the datums are different, I can require that the data frame in the MXD file have a transformation method specified in its properties, e.g. ED_1950_To_WGS_1984_10. If that was the end of the story, my job would be simple: if the data frame has a transformation method in its properties, pass that to `Project_management()`, otherwise omit that parameter. The script would just look at the value of `dataFrame.geographicTransformations`. The problem is that on some machines, when I ask for `mydataframe.geographicTransformations`, **I get a spurious value of "NAD_1927_To_NAD_1983_NADCON".** So if the data frame has no transformation methods specified in its properties dialog, `dataFrame.geographicTransformations` returns `["NAD_1927_To_NAD_1983_NADCON"]`. Which causes `Project_management()` to fail, if that transformation doesn't apply to the GCS and PCS actually being used. For example, when the data's GCS is GCS_WGS_1984, and the data frame's PCS is Asia_Lambert_Conformal_Conic, no geo transformation is needed (both have datum D_WGS_1984), and none is specified. So if I pass `dataFrame.geographicTransformations[0]` to `Project_management()` as the geographic transformation parameter, I will cause it to fail because the NAD transformation is inapplicable. If the data frame _does_ have a transformation method specified in its properties dialog, then (on these same machines) `dataFrame.geographicTransformations` returns a list with "NAD_1927_To_NAD_1983_NADCON" as the first element, followed by the one specified in the properties dialogue. Again, this is on some machines. On others, the spurious "NAD_1927_To_NAD_1983_NADCON" doesn't appear, so I can rely on `dataFrame.geographicTransformations` to be set, or not, as appropriate. So my question is, where is this spurious transformation coming from, and how do I get rid of it, so that if I don't need a transformation, I can avoid passing a bad one to Project_management()? There's nothing in the Environment Settings under Output Coordinates / Geographic Transformations. Or, if I can't get rid of it ... maybe I can detect when the data's GCS and the data frame's PCS have the same datum, so I don't need a geographic transformation. In that case I would just not pass any such parameter to `Project_management()`. The trouble with that approach is, I don't know how to get the datum, or the GCS, of a PCS spatial reference. I can look at `dataframe.spatialReference.PCSCode`, but `spatialReference.GCSCode`, `GCSname`, `datumCode` and `datumName` are blank (see http://help.arcgis.com/en/arcgisdesk...00000p6000000/ \- these properties are only available with a GCS). I know that **ArcGIS desktop knows** the GCS being used by the data frame, because when I open data frame properties, the Coordinate System tab tells me not only the PCS but also the GCS (including datum). But the **GCS is not available** thru `dataframe.spatialReference`, that I can tell. So, if someone could tell me   * how to accurately find the specified geographic transformation of the data frame (i.e. get rid of the uninvited NAD geographic transformation) (this would be the preferred solution); and/or   * how to determine the GCS or datum of a PCS (can my script look it up in some standard place, based on the PCS code?) ... I would be very appreciative.


How to import a simulation from ParaView to ArcMap?

I made a simulation and reworked the output with **ParaView**. I simulated an avalanche on a DEM and what I want to do is to compare the area of the simulated avalanche with that of the real avalanche. To do this I would use ArcMap 10 where I already have the real area as a shapefile. How can I bring what I have on ParaView to ArcMap? Thank you very much, Michele


Get road type (clazz) in osm2po

Can I get clazz of a segment in Java? For example,               RoutingResultSegment rrs = graph.lookupSegment(path[j])      How can I get `rrs`'s clazz?


How to extract zipped shapefile and check its shape type using arcpy module without writing on disk

I have a zipped shapefile saved on my disk. I want to extract it in memory and check the shape type. I'm using arcpy module and want to stick to only default modules comes with python installation. I tried something like ` import cStringIO, zipfile, arcpy stringIO = cStringIO.StringIO() file_ = r"C:\Polyline.zip" zipFile = zipfile.ZipFile(file_) for name in zipFile.namelist(): if name.endswith(".shp"): stringIO.writelines(zipFile.read(name)) print arcpy.Describe(stringIO).shapeType ` But it is not working. I got the trace back: ` Traceback (most recent call last): File "", line 8, in File "C:\Program Files (x86)\ArcGIS\Desktop10.1\arcpy\arcpy_ _init_ _.py", line 1190, in Describe return gp.describe(value) File "C:\Program Files (x86)\ArcGIS\Desktop10.1\arcpy\arcpy\geoprocessing_base.py", line 374, in describe self._gp.Describe(*gp_fixargs(args, True))) RuntimeError: Object: Describe input value is not valid type ` How can I unzip a shapefile in memory and finds its shape type?


Does map export using GarminCustomMap plugin for QGIS work on Mac OS?

When I try to export a map using GarminCustomMap (1.4) plugin QGIS 2.01 on Mac OS 10.9 quits with no messages. I suspect this might be a Maverick related issue. Is anyone successfully using the plugin with 10.9? I have reported the crashes on the plugin issues page.


Opening InMemory Table to execute some SQL gives Runtime Error 445 - Object doesn't support this action?

I have a InMemory that gets created, what I then want to do is run INSERT, UPDATE and DELETE queries on the table. The code i have at the moment is below, I create a new InMemory table, add it to the map (which works).                   Dim pInMemTable As ITable         'create Temp DP table         Set pWSF = New InMemoryWorkspaceFactory         Dim PName As esriSystem.IName         Set PName = pWSF.Create("fred", "InMemory", Nothing, 0)         ' Debug.Print ">>" + PName.NameString + "<<"         Set pFWS = PName.Open              Dim pFlds As IFieldsEdit         Set pFlds = New esriGeoDatabase.Fields         pFlds.AddField MakeField("OBJECTID", esriFieldTypeOID)         pFlds.AddField MakeField("DP_REF", esriFieldTypeString, 10)         pFlds.AddField MakeField("WANTED", esriFieldTypeString, 1)              Dim pUID As New UID         pUID.Value = "esriGeodatabase.Object"              Dim ptable As ITable         Set ptable = pFWS.CreateTable("DP_Temp", pFlds, pUID, Nothing, "")        ' Debug.Print ptable.              Dim pSAT As IStandaloneTable         Set pSAT = New StandaloneTable         Set pSAT.Table = ptable         pSAT.Name = "DP_Temp999"              Dim pMxDoc As IMxDocument         Set pMxDoc = ThisDocument              Dim pSATC As IStandaloneTableCollection         Set pSATC = pMxDoc.FocusMap              pSATC.AddStandaloneTable pSAT         pMxDoc.UpdateContents              Set pInMemTable = pSAT              Dim path As String         path = "InMemory/DP_Temp999"         If pWSF Is Nothing Then             MsgBox "pwsf is nothing"         Else            ' MsgBox "pwsf is something"             'MsgBox pWSF.WorkspaceDescription(False)                  'pwsf.Open(             Set pWS = pWSF.OpenFromFile(path, 0) ' **********PWS doesnt get set***********             If pWS Is Nothing Then                 MsgBox "PWS IS NOTHING"              Else                 MsgBox ("Yes!!")             End If              End If      but the code fails where PWS is getting set as it always returns as Nothing. I have tried various ways of setting the path value e.g. "InMemory/DP_Temp999", "In_Memory/DP_Temp" and trying the "\" as opposed to "/" but it still comes back as nothing. Does anyone have any ideas please? or is there a better way of executing SQL on a table that will allow SELECT, DELETE and INSERTS. ******UPDATE**** Right I have been working on this for a few days have have been able to set pWS, by casting by ptable to a idataset object and then setting pWS to the pdataset propertyset properties which works. However I have then tried the following code, to see if SQL can be executed                   If pWS Is Nothing Then             MsgBox "PWS IS NOTHING"         Else             'MsgBox ("Yes!!")             Dim pWorkspaceProperties As IWorkspaceProperties             Set pWorkspaceProperties = pWS                  Dim pWSProperty As IWorkspaceProperty             Set pWSProperty = pWorkspaceProperties.Property(esriWorkspacePropertyGroup, esriWorkspacePropCanExecuteSQL)             If pWSProperty.IsSupported = True Then                 MsgBox "SQL is supported"                 Dim sql As String                 'sql = "select * from 'DP_TEMP';"                 sql = "UPDATE dp_temp set dp_temp.[DP_REF] = 'Test';"                 pWS.ExecuteSQL sql             Else                 MsgBox "SQL is not supported"             End If         End If      When the code is run, the pWSProperty.IsSupported returns as True, however when the .ExecuteSQL method is run, I get the following error Runtime Error 445 - Object doesn't support this action. So I am a bit confused as to whether SQL can be run or not.


QGIS 1.8 compiled against PyQt 4.8+?

I'm currently running QGIS on Win7 64bit. My installation says it was compiled against and running against PyQt 4.7.1, although according to the Osgeo4W installer and PyQt4.QtCore.PYQT_VERSION_STR, I actually have PyQt 4.8.3 installed under the Osgeo4w directory. I'm running into a bug in my PyQt code (related to QNetworkAccessManager) that, according to the bug tracking log, was fixed in PyQt 4.8.0. It seems, then, that the PyQt 4.7.1 module is being used, and the 4.8.3 module is being ignored. I'd like to resolve this issue. I've seen a few posts on various forums that suggest that the generic distribution of QGIS 1.8 is compiled against PyQt 4.8. I've tried reinstalling QGIS using the Osgeo4w installer after removing all references to PyQt from my system path, but I still end up with QGIS compiled and running against PyQt 4.7.1. Could anyone shed some light on this issue for me? Should QGIS 1.8 installed through Osgeo indeed be a distribution compiled against PyQt 4.8+? That is, am I doing something wrong in the installation process? Or will I need to compile QGIS manually if I want it to be compiled against PyQt 4.8? Thanks very much. Alex


Creating a vector based design from a raster database

I have a raster table that I've created using Postgis raster. It has 365 records and the table looks like this:               +-----+------------+---------+------------+     | rid |  rast      |filename |date_of data|     +-----+------------+---------+------------+     | 100 |            | file56  |2002-01-09  |     | 101 |            | file57  |2002-01-10  |     +-----+------------+---------+------------+      The second column _rast_ is raster files in WKT format and data type is _raster_. All my raster files have the same extents and cellsizes. I want to turn this raster design into a vector design like the following:               +-----+---------------+----------+     | ID  |pixel_centroid |array_366 |     +-----+---------------+----------+     | 1200|               |{26,20,..}|     | 1201|               |{28,68...}|        +-----+---------------+----------+      Imagine the raster files have 1600 pixels and accordingly I will have 1600 records in my new table. 'Pixel_centroid' is the point geometry of the pixels and could be calculated using               select ST_centroid((ST_PixelAsPolygons(rast, 1)).geom)  from raster_table      and array_366 is pixel values of all the records in the order of date_of_data:               select (ST_PixelAsPolygons(rast, 1)).val  from raster_table      might be used but I don't know how to put all the values in their place and populate the array.


Google Earth plug-in: how to get name/description of a model, from lat/long or click

In Google Earth, if you move the mouse over a 3D model (building), it turns light blue. If you click on it, you get a description (or sometimes just an approximate address). Google Earth plug-in web apps like Monster Milktruck or First Person Cam are wonderful educational tools for exploring an area, and learning about it spatially. For example you are walking or driving around Rome, and you see buildings like the Coliseum or the Forum, but you don't know what they are; so it would be very useful to be able to click on them as you can in Google Earth, and get a description. Is there any way to do this in the G Earth Plug-in API? I see that a hitTest can tell you whether you clicked on a model, and if so, what the lat/lon of your hit is. But surprisingly, it doesn't tell us _which_ model was hit. Nor can I find any way to look up the model in Google 3D Warehouse using the lat/lon. Can someone tell me a way to do this? Or is it for sure not allowed?


Does Bing maps silverlight API work on Windows phone 7?

Has anyone seen any good resources on using Bing Maps Silverlight APIs on WP7?


Can I get the metadata about the themes in Oracle MapViewer?

I am trying to use the MapViewer xml api to generate a dynamic theme (for generic highlighting support), but don't want to expose table names and such to users, I would like then to communicate in terms of themes. Is there some way that I can obtain the mappings from theme to table so that I can generate dynamic themes based on `jdbc_query` element? The `jdbc_query` needs me to give it some SQL which identifies the instances to be included in the dynamic theme. I want the user to be able to specify a theme name and then be able to construct an sql statement which gets the geometry of the nearest neighbour to a particular point on that theme. To do this I need to know which table the data in the theme came from and which column contains the geometry in that table. Can I find that information using an `info_request` perhaps?


bicycle / bike map road grade markers - semi-automated?

So I'm making a bike commuting map for our area, and one of the features to show (common on a lot of bike maps) is some type of road grade marker. We're planning to show three categories: 'kinda steep, pretty steep, and crazy steep' - we're not planning to use those names - but you get the drift. (Actually, four categories - anything less than the first category just doesn't get a marker.) The markers are one, two, or three arrows - here's a clip: ![enter image description here](http://i.stack.imgur.com/0an1D.png) The problem: since we live in the foothills - not much flat stuff around here - the task of getting all the grade markers right is pretty daunting and tedious. So - it would be nice to automate it. So I've got the DEM for the area from the NED, and got a nice slope raster from it (thanks to folks that helped figure that out). Overlaying that on the bike map (the roads are a vector layer of course), you can see a way to conceptually automate this task. Basically, making a vector layer that is the boolean AND of the roads layer and the slope layer of certain values (whether by categorizing the slope layer, or by making separate slope layers per grade range ahead of time) should get it done. I'd want the output to be a vector layer. But, I'm not quite clear on how to do all that. Suggestions on workflow? I realize this concept doesn't account for direction (which way is uphill) so I'll plan on doing that by hand (flipping the vector selections as needed, or, doing so with an attribute that is used by the style to decide which way to draw the arrows) UPDATE: this will be done in QGIS 2. A thought: make a vector layer (of shapes) out of raster categories, i.e. a vector shape of areas with slope between 5 and 10 degrees (this would cover one of the categories of steepness), then do a boolean AND of that layer with the roads layer; so now you have a bunch of disjoint lines; that's exactly what is needed; just style these lines with the double-arrow (or as appropriate for that category), and flip direction of each segment as needed (which is easy to do from memory or from a map if needed).


Force a collection of unknown number of points to be a polygon

I am working with point clusters and want to create the rough shape of the cluster in the form of a polygon using `ST_ConvexHull()`. I'm having problems with single outliers though, as they don't convert to a polygon but to a `POINT`:               SELECT ST_AsText(ST_ConvexHull(ST_Collect(point))) ...          POLYGON((13.01741532 52.40118295,13.382100866593 52.516280103729, ...))     POINT(12.5 51.5)      Is there a way to force small geometries (single or double points) into a polygon? I want to simplify the subsequent processing by guaranteeing polygons; accuracy is not a high priority in this case, making a minimal triangle out of a single point is fine, even if not exact.


Looking for a web service for change detection

I'm looking for service that provides me a change detection service based on a database of remotely sensed data. Has anybody happen to know such a service? The service I am looking for is a geoportal that alarms the user of a specific change (like oil spill or major coastal change or even fire in forest) in a specific area or globally. It can be either freely available service or on- demand service.


How to determine intersection between geojson feature and given coordinates?

I'm overlaying a leaflet map with features of geojson, fetched with d3 and added as SVG elements to the DOM. My geojson contains too many features to render at once, so my goal is to add to the map only features, that intersect current cursor location, so I can have a smooth hover effect. I can get current map coordinates of the cursor from a mousemove event, but I have to determine if current coordinates intersect a feature in geojson and exactly what feature is intersected. Thinking further, maybe I have to make a small buffer, so it's easier for the user to hit the closest feature. Is there a fast way in leaflet or d3 (or even pure javascript) to get the intersected feature from geojson feature collection?


FEMA LOMA Data Availability?

Does anybody know where there might be some good/current FEMA data? I am interested specifically in the LOMA data for Florida. I have found a WMS Service that is hosted by FEMA but I'm having trouble actually bringing it into ArcMap.


Use Display Expression for KML Popup, SHP to KML ArcGIS

In Arcmap I've created a complex label expression in the Display expression for my layer. I'd like to have that expression displayed when I convert the layer to .kmz file for display in Google Earth. I'm using the ESRI Layer TO KML conversion tool. What happens in the popup right now is the display expression jammed, without the correct formatting, into the feature name box and also the feature description box, followed by the complete list of fields. I've been looking at the HTML popup tab, but don't see any way to use the Display Expression I already have in that tab. Is it possible for this to happen? It looks like otherwise I'll have to write a XSL template. EDIT: Thanks to @KHibma for the comment. What I've ended up doing is   1. removing all of the fields from display in the shapefile properties in ArcMap before the export so they are not included in the kml popup   2. Labeling the features in arcmap with one field so that field value is the "title" of the feature in the popup.   3. adding HTML formatting tags to the Display Expression (which I was writing in Python). this expression appeared, formatted, in the blue box below the feature title. I can't add a screenshot of the result right now, but will be able add it later. This worked fine on a small scale, but when trying to export an entire parcel fabric for a county I would get a C++ runtime error. One interesting thing I found was that, even before the export, if there was an error in the FindLabel function and it was not able to create a label, it would come back as empty. So I added a try/except statement within the function like:               def FindLabel([field1],[field2]):         try:             #code for function         except:             return "there was a problem creating this display"      which worked great. However, when I tried to import sys or traceback in the function and set the error handling message as the return in the except statement, it still returned blank (i.e. didn't return the text of the error handling message). In the end, I fixed all the problems with the expression, but still got the runtime error when making the big export.


Is there any accuracy difference between standard DTM editing tools like DT Master and automatic removal of artefacts using ERDAS?

I would like to know whether there is any accuracy difference while editing DTM points using standard photogrammetry software like Inpho DT Master and automatics removal of artefacts in ERDAS imagine.


When I used geowebcache, my map overviewmap can not use it?

I know: gwc / service / wms geowebcache use data, overviewmap is use wms original data, I tried modifying openlayers.js,some difficulties. You encountered the same problem and can help me look at it, thank you！


QGIS conversion of Microstation .DGN to .SHP

I'm attempting to convert a .dgn file to .shp. I'm reading that support is only provided for v7 microstation files, is this still the case? I've been using ogr2ogr -f "ESRI Shapefile" output.shp input.dgn and it throws the following error. Unable to open datasource `input.dgn' with the following drivers. -> ESRI Shapefile -> MapInfo File -> UK .NTF ... The list goes on. ANy help is appreciated.


What are the requirements for using Google's GEvent?

For a current project, we need to know when a Google map is done loading all of its tiles. I found the following question, which was for the same basic problem: register event "loadend" on layer So now, we're using the code from that answer:               GEvent.addListener(mGoogleHybridLayer.mapObject, "tilesloaded", function() {      console.log("Tiles loaded");     });      However, this gives a reference error for `GEvent`. Are there some extra requirements for using `GEvent`? Do we need to import it from somewhere? I can't post code we're using, because there isn't any yet - even in the browser's debug console, `GEvent` is undefined. Using Openlayers 2.11. EDIT: maybe I should restate the question into one that can actually be answered: what are all the requirements for being able to use GEvent?


Can I overlay layers in Dealul Piscului 1970 and WGS84 in QGIS?

I have a 10 Gb ECW file (mosaic of a whole county) that is in Oblique Stereographic projection with Datum: Dealul Piscului 1970. The problem was that no matter the Default CRS for new project was, when I wanted to add this raster the coordinates were all wrong, although the input projection for the raster was correct. This problem was solved by deleting the ECW's ".xml" file that was accompanying the ECW file. Once deleted, the problem was solved and the QGIS was showing the right coordinates. Once the above problem was solved. I wanted to add a vector layer of 30 mb that represents the whole waterway of the Danube Delta plus inland towns. This layer is in WGS84. Now when I try to import this layer does not overlay with the ECW file. The ecw file have been defined with the right projection also the vector layer have been defined with the right projection. In settings/options/CRS/ the "on the fly" is checked and the project projection is WGS84 and also tried the vice-verse with Oblique Sterographic Dealul Piscului. With all the permutations available I still can not overlay the layers. Do you have any other suggestions? Thank you


Convert a number in attribute table to binary

I would like to convert whole integers from decimal to binary in ArcGIS using the field calculator. I am looking for the equivalent of DEC2BIN excel command. All values are 255 or below so 8 bit is fine. It must have 0 padding to the left. So for example 4 must be 00000100. I can do it by exporting the attribute table to Excel and using the Dec2Bin command in Excel and rejoin the data but I would like it in a model so want to use the field calculator. Dec2Bin (Value, 8) works in Excel. Attached is how I would like it to look when complete. ![enter image description here](http://i.stack.imgur.com/DiMme.png) So basically change the value to and 8 digit binary value stored as a text field with leading 0s so all field have 8 values. I found this code in arcscripts but cannot get it to function. Dim valI as integer Dim valB as string Dim x as integer valI = [fieldName] for x = 7 to 0 step -1 if valI >= 2^x then valB = valB & 1 valI = valI - 2^x else valB = valB & 0 end if Next valB = valB & "00000000" fieldName = valB


Can I choose specific SDE version to create replica geodatabase from?

I'm trying to automate the process of exporting a replica geodatabase using the geoprocessing tools and eventually a Python script. I can't figure out how to create the replica from a specific version. Is this possible?


mapCanvas.saveAsImage not functioning as expected in QGIS 1.8

I have a plugin for exporting maps as PNG images. I believe this code was working correctly within 1.6 (but I cannot guarantee I haven't made another change that would have broken it there, too). Basically, I load a set of base map features and then loop through a set of point layers that get applied and the whole thing is exported as an image. I'm not getting the full canvas at all, but only a single attribute loaded from the base map and nothing from the point layers. (bad result, expected result) In the process of debugging this, I discovered a curious fact. If I include a QMessageBox.information right before I do the export, it works correctly! Doesn't work:                   self.mapCanvas.resize(QSize(394,350))         self.mapCanvas.zoomToFullExtent()         currentImagePath = theExportPath + folder         currentImageName = attribute.replace(' ', '_') + ".png"         self.mapCanvas.saveAsImage(currentImagePath + "/" + currentImageName)      Does work:                   self.mapCanvas.resize(QSize(394,350))         self.mapCanvas.zoomToFullExtent()         QMessageBox.information(self.iface.mainWindow(),"Debug","Break")         currentImagePath = theExportPath + folder         currentImageName = attribute.replace(' ', '_') + ".png"         self.mapCanvas.saveAsImage(currentImagePath + "/" + currentImageName)       I tried to add a simple self.iface.mainWindow() in place of the message box but that didn't do the trick. I'm thinking there is some kind of refresh I could do, too, but I tried mapCanvas.refresh() and mapCanvas.updateMap() and these in combination with mainWindow() but could get nothing to work. Did this change in 1.8? Any suggestions to make it work? Including the QMessageBox call will get me going but I have to press "OK" for each map to be exported. With all the export combinations, I'd be doing this more than 1000 times! Thanks! B


How to create a multipart polygon?

I have two polygons that consist of two separate polygons/lines in QGIS). How can I make this into a multipart polygon, which consists of 1 line in QGIS?


Cutting features in ArcGIS 10

What´s the easiest/smartest way of getting the opposite output of the clipping tool? I just need the area of feature _A1_ which **is not** covered by feature _A2_ in a new Layer/File. Now, I get the result by union´ing both features and deleting the "wrong" parts by hand.


how to make transparent layers with Google earth layer?

with your help, i solved most of my problems, but here is another one - i guess, a smaller issue... I managed to get Google earth layer, but when I add my GPS gpx files, are all hidden under the Google map... How to put them onto Google Earth layer? thank you in advance (and sorry for being a dumbass) Krepki


Why does intersect not calculate polygon areas?

this is strange: I intersected 2 shape-files and the resulting shape does _not_ contain the Shape_Area and Shape_Length of the intersected polygons. Instead the field Shape_Area contains the polygon-area of one of the 2 original Layers. Recalc of the geometries with vector->geometry->add/export column seems to add the areas. That is quite a drag and very slow for my 350.000 polygons. So I am hoping for a systematic approach that adds the areas while generating the intersect. Is intersect done by QGIS (qgis.bin) or by a gdal-tool? Either way: is there a "calculate areas"-option hidden somewhere? THX, stn


Extract polygon from several FC in one process

I've built a polygon and used Clip (Analysis Tools) in order to extract information from several FC that intersect the polygon. Of course, I used the clip tool for each feature class. My boss wants a single tool that clips all these feature clases at once and saves the result in different FC. Is there a tool in ArcGIS that does that? Select by Graphics is cool, but then I have to export the data "manually". I've created this model, but in the case of the model I don't have the liberty of selecting the area I want to clip everytime. Thank you


Good place to get Country, State, City, and Zipcode information?

I'd like to populate a database with common location information. I'm not entirely certain that addresses across the globe have similar or even the same entities as the Country, State, City, and Zipcode that we are so familiar with here in the United States. Does anyone know of a web service or other method that would be useful for obtaining such data and populating a local database with its contents? (or conversely just contacting the web service to obtain this information as needed)


Housenumber(point) angle from line (street)

using arcgis 10 with Point data (Housenumber without angle attribute) and line data (street without angle attribute) I want to create housenumber labels from point with correct angle from nearest streetline. is there any toolbox or script available for this, and if there is a script can you guide me step by step beacuse I dont know where to put script, and python would be the best I believe. my thinking of this is calculate AVERAGE streetline angle and then somehow put this values to nearest point? any help is good help :)


Problem in filling table of ArcSDE feature class in SQL Server 2008?

I have an ArcSDE point feature class and I add it to ArcMap 10. After that I add some points to it and fill every field for all points in the attribute table and Save Editing. However, the SQL Server table that related to the feature class is empty while the attribute table in ArcMap and ArcCatalog is full (as expected). I checked the connection to the spatial database and that was correct. How do I fill the table in SQL Server when the attribute table in ArcGIS is full?


How can I visualize a route in QGIS?

I have a route in QGIS that visits every road/edge at least once. By using a marker line I can see the direction of the route, but it doesn't show the order of the segments (see the screenshot below). Is there any style/visualization that will allow me to see the order of the roads in the route? Currently the route is a single polyline, but I can split it into segments if that will make it easier (the route is generated with Python code). ![Route example](http://i.stack.imgur.com/reec7.png)


How can I programatically create and add features to a memory layer in QGIS 1.9?

I had a working plugin in QGIS 1.8 which read data from an MSAccess database and added it to a series of memory layers. There is some processing involved in the interim so I don't think simply using QGIS to read directly from the database is an option. I would like to move from QGIS 1.8 to 1.9 (primarily due to the improved quality of the print composer). The plugin doesn't work with the new API. I've tried a variety of methods that came up in google searches. One, modifying the below code - from http://www.qgis.org/pyqgis- cookbook/vector.html#memory-provider, i.e. adding geometry and attributes to the dataprovider and then updating the layer - to suit the new API worked a little but the attributes were not visible until I manully entered editing mode (similar to http://hub.qgis.org/issues/3713). An alternative approach, detailed in reply #1 of the above link, added the layer and attributes correctly but I was unable to add features to the layer. Given that it should be a pretty simple task I'm hopeful that someone on here can offer a working example of how this should be done. (P.S. I am not a professional programmer and most of my coding is fairly crude - I welcome any guidance but ask that you excuse some ignorance on my part)               # Receivers = a list of lists returned from a database query          # create layer     vl = QgsVectorLayer("Point", item, "memory")     pr = vl.dataProvider()          # add fields     pr.addAttributes( [ QgsField("Rec_No", QVariant.Int), QgsField("Include",  QVariant.String), QgsField("Label",  QVariant.String), QgsField("X", QVariant.Double),                         QgsField("Y", QVariant.Double), QgsField("Z", QVariant.Double), QgsField("Height", QVariant.Double),                         QgsField("Project_Re", QVariant.String), QgsField("NCA", QVariant.String),                         QgsField("DayCrit", QVariant.Int), QgsField("EveCrit", QVariant.Int), QgsField("NightCrit", QVariant.Int) ] )          for i in range(len(Receivers)):                 # add a feature       fet = QgsFeature()       X = Receivers[i][3]       Y = Receivers[i][4]       fet.setGeometry( QgsGeometry.fromPoint(QgsPoint(X,Y)) )            # Details = a list of results returned from a database query specific to each result in 'Receivers'            if Receivers[i][3] != 0:         Include = 'Yes'       else:         Include = 'No'            fet.setAttributeMap( { 0 : QVariant(Receivers[i][0]), 1 : QVariant(Include), 2 : QVariant(Receivers[i][2]),                              3 : QVariant(Receivers[i][3]), 4 : QVariant(Receivers[i][4]), 5 : QVariant(Receivers[i][5]), 6 : QVariant(Receivers[i][6]),                              7 : QVariant(Details[0]), 8 : QVariant(Details[1]), 9 : QVariant(Details[2]), 10 : QVariant(Details[3]), 11 : QVariant(Details[4]) } )       pr.addFeatures( [ fet ] )          # add a style     vl.loadNamedStyle('C:/OSGeo4W/apps/qgis/python/plugins/Gopher2QGIS/styles/Receiver_Style.qml')          # update layer's extent when new features have been added     # because change of extent in provider is not propagated to the layer     vl.commitChanges()     vl.updateExtents()     vl.updateFieldMap()          QgsMapLayerRegistry.instance().addMapLayer(vl)


Web mapping application with time slider and raster data

I need to create a web mapping application that can show raster layer of historical maps with a time slider. Basically, something like this: http://goo.gl/GYHOX I have different maps of different years of the same are, that can be showed on a time series and moving on with a time slider. The interface can be simple, the important thing is the time slider. I have programming experience only in HTML/CSS, little of PHP. Any suggestions of which tools to use? I had a look at TileMill but the problem is the time-slider that I have no clue on how to build it. Thanks in advance for who will help me.


How specify a corridor minwidth?

I’m trying to make a corridor between two substations with ArcGIS, a macro corridor and a micro corridor. I have two questions in this respect and would be grateful for your suggestions. a) How can I specify, that my corridor must have a minwidth of 100m? b) How can I make a linestraighting of my costpath within microcorridor? Thank you very much.


How to output visually correlated clusters using Morans I Spatial Autocorrelation?

![this is a screenshot of the geoprocessing output message](http://i.stack.imgur.com/Lqt7v.png) ![this is a screenshot of the attribute table result](http://i.stack.imgur.com/aDTkC.png) How do I use **Morans** I to output visually correlated clusters? I can seemingly generate statistical data, but not one that I know how to visually output on a map. I'm trying to run a morans I analysis on a collection of municipal level boundaries. As I understand the Morans I process, the 'input' feature class is whatever layer is going to be auto correlated, and the input field is whatever variable is going to be measured for the auto-correlation. In the case of the attached photos, I tried to auto-correlate on the basis of the variable R_vicMargNum. Approximately 35 out of the 80 rows in my attribute table do not have any values for that variable; I'm still hunting for the data which can properly fill the fields. 1) I don't understand why the fields which had a value for that variable are 'null' for columsn f7-f9. 2) I don't understand the value of the statistical output.


Change number of decimal digits for features in ArcGIS 10?

I am editing vertices in polygons and points in ArcMap 10, and realised that they only hold 3 decimal digits for the X, Y coordinates, but I would like at least 6 decimal digits. Even if I type in more digits, it gets rounded up to 3. How do I change this? Thanks.


How to create a buffer from a click on map on ARCMap

I'm searching a way to creating a buffer of a specified radius starting froma a click on the map.


Densify a JTS LINESTRING (add points)

I have a LineString I would like to "soften" in a way that if two coordinates (vertices) in the LineString are very far apart, I want to **automatically** add some coordinates in between. The line string should stay the same but with more coordinates describing it. Do you know of a way? ![enter image description here](http://i.stack.imgur.com/1FLCj.png) I have a LineString like the one on the left - and I want to add some coordinates in the between the existing coordinates (like on the right) - hopefully by JTS and hopefully automatically Thanks!!


Flex Identify tool version 3 1085 error

I am building a flex web map app. I am using the Identify tool version 3, here. Below is the syntax I am using, based on the PDF help file. It, however, gives a 1085 error. Does anyone know if I am using the correct syntax?               <layers onlythese="true">         <layer>             <name>Australia offices</name>             <fields>                 <field name="Address" />                                 <field name="City" />                                 <field name="State" />                                 <field name="Postcode" />                         </fields>             <links>                 <link includeinresults="false" field="" alias="" tooltip=""/>                      <linkprefix></linkprefix>                     <linksuffix></linksuffix>                     <iconfield></iconfield>                     <iconprefix></iconprefix>                     <iconsuffix></iconsuffix>                 </link>             </links>             <zoomscale>15000</zoomscale>             <forcescale>true</forcescale>                  </layer>                  <layer>                         <name>QLD rail</name>             <fields>                 <field name="FNAME" alias="Name"/>                                 <field name="ELECTRIFIE" alias="Electrified"/>                         </fields>             <links>                 <link includeinresults="false" field="" alias="" tooltip=""/>                      <linkprefix></linkprefix>                     <linksuffix></linksuffix>                     <iconfield></iconfield>                     <iconprefix></iconprefix>                     <iconsuffix></iconsuffix>                 </link>             </links>             <zoomscale>15000</zoomscale>             <forcescale>true</forcescale>         </layer>     </layers>


Openlayers map changelayer event firing twice for params

I have registered **changelayer** event of my map.                   var registerEvents = function () {             map.events.register("changelayer", null, onMapLayerChanged, true);         };              var onMapLayerChanged = function (event) {             console.log(event);         };      when I changed a layer's params as following, console writes two times.               layer.mergeNewParams({'CQL_FILTER': "name=wate"});


Problem with OpenLayers Tile Stitching Printing - All Black Image

I am trying to get an image of my OpenLayer application. I am using this http://www.mapsportal.org/olprint/ and http://trac.osgeo.org/openlayers/wiki/TileStitchingPrinting as a guide. When I click my print button I get a print failed message and when I go to the TMP folder (c:\wamp\www\tmp) there are IMG files and a JPEG but the JPEG is all black.


leaflet and geojson tiles

New to mapping. Trying to create a heat map with data that updates every hour for US Counties with leaflet.js similar to http://leaflet.cloudmade.com/examples/choropleth- example.html. But, adding all the US Counties in 1 shot makes map jerky in IE. So found this example of using GeoJSON tiles with leaflet.js, http://bl.ocks.org/3725681 Only problem is that many tiles are not rendered completely. I think some of the z/x/y.json files are missing coordinates, but not sure. Can someone point me to the right direction of the problem here.


Count the number of symbols in a 'new symbology' renderer from pyqgis

Background : I create a lot of symbology using 'saturation' ramps - for want of a better term, e.g. where the hue and value of each symbol stays the same, but the saturation changes. I couldn't figure out how to do this easily in qgis (EDIT: see comment - its easy when you know how), and each time I classify items, the colours reset to a colour ramp rather than a saturation ramp. I figured this should be relatively easy to do using pyqgis, which it turns out has been, using this code (built for 5 symbols):               from PyQt4.QtGui import QColor              # import this namespace as I want to refer to QColor later     layer = qgis.utils.iface.activeLayer()     i = 0                                       # used to increment the saturation value     for sym in layer.rendererV2.symbols():      # loop through each symbol in the layer using 'new symbology'         i += 50                                          sym.setColor(QColor.fromHsv(40,i,255))      Problem: . The above code works fine, but is not dynamic, in that I've set it up to increase the saturation by 50 each time (giving 5 graduations between 0 and 255). I want to set the saturation to graduate dynamically based on the total number of symbols. The mathematics of this are pretty straight-forward, but I'm not able to get the count from the symbols() method. I've tried:               >>> rendererV2.symbols().count()     Traceback (most recent call last):       File "<input>", line 1, in <module>     TypeError: count() takes exactly one argument (0 given)      ... which as you can see returns an error (although the QT documentation seems to list count() as a method of a QList<> object). I've tried:               >>> rendererV2.symbols().count(QgsSymbol)     0      ... which as you can see returns 0, which isn't true. There are 5 items in the list. I'm sure I must be doing something daft, but being new to python, pyqgis and qt, it is beyond my ability to see what at present. Help?


How to save QGIS composer template with image - e.g. company logo?

I can't figure out how to save QGIS Print composer template with image. In our map compositions we always have company logo in one place. It is easy to add image (logo) to composer and place it to spot I need. But when I save template, it saves without image (logo). Is there some setting for it? Am I misssing something or there is no way so far to save template including image?


How to Project shapefile from NAD 1983 UTM Zone 12 to WGS 1984 UTM Zone 13N?

I have a question and I have a shapefile that is in NAD UTM Zone 12 and Datum is North American Datum of 1983. What I would liket to reprojection it to NAD 1983 UTM Zone 13N and change the Datum to WGS 84 ? I have tried different methods but cane out the same. it seems like the Datum doesn't change, is it ? Should it be just one step to use the project to reproject to NAD83 UTM Zone 13 ? or is there more than just one step to get it there ? Currently I have 9.3.1 .


How can I conduct Geographically Weighted Analysis using an exponential decay function?

I have some points (representing 30 study sites) and I want to calculate the weighted mean of several variables (landscape layers) using a negative exponential decay (weights function) to give more weight to the landscape features closer to the points (study sites) than the features (grid cells) further away. (The series of underlying grid layers represent landscape features such as roads; habitat etc.) Also, how can I use spatial statistics to find the optimal distance (bandwidth) for my distance-decay function? In the literature everyone resorts to creating buffers (vector) at different distances (say 1km, 3km & 5km) to quantify the variables and then uses statistics to determine which buffer distance is significant. Other ecologists (Rhodes et al. 2006) have used the “negative exponential distance weighted density” of each variable and provided a scale parameter (L) which controls how rapidly the influence (ie weighting) of the variable declines with distance, eg L = EXP(-0.002*Distance) I’m using ArcGIS spatial statistics; grid (raster) spatial analysis tools; and investigating GWmodel for R-stats from an earlier post.


mapping and update data in database

We have two feature layers:   * layer A - drainage segment layer (along road outline)    * layer B - road centre line layer. In my SQL Server database I have a table storing the mapping information of segment ID (from layer A) and Road name (from layer B). ![enter image description here](http://i.stack.imgur.com/3sn5F.jpg) The problem is that when I create a new polyline feature (a segment) in layer A, I need to update the database table also. I am stuck here. I do not know how to update the database table at the same time as when a user creates a new feature in layer A. Should I create a new geoprocess (using python script) to do so? Is there any other solution for it? I am using vb.net arcmap addin (vs 2012), arcmap 10.2 and sql server 2008 r2.


Error in Conditional Tool of ArcGIS for Desktop?

I am trying to extract stream raster from a flow accumulation raster using Conditional Tool in Spatial Analyst. Here's my input:   1. Input Raster: FLOWACC.tif   2. Expression: Value >= 3.5   3. True Raster: 1   4. Output raster: Stream.tif Error gives me invalid expression. Basically my objective is to isolate all values >=3.5 and consider them as stream raster. I am using ArcGIS 10.2 for Desktop.


Qgis doufur, query problem

i've a question about a query I need to query results that are multiple of 50 like 100,150,200,300 How can i perform the right query?


How to create single contour shapefile from a contour shapefile?

Is there a way to create single contour shapefiles from a main contour shapefile. The main shapefile has quite a few contours and I need to isolate them and create individual files by elevation. Selecting contour by contour and creating a shapefile would take a while. Anybody knows of a faster way? Thanks!


Interacting between Openlayers and Python

I have a openlayers website like this one that gives the user the ability to draw polygons and then creates a GeoJSON file. I want to process this GeoJSON with a python script. The python script returns in the end a string that should be returned to the website and displayed as a info box. What kind of frameworks/apps/steps should I look into to accomplish this project?