Are there GUI resources for gdal and ogr?

For those who are new to the programming side of GIS, I feel that a lot of great tools and resources are inacessible. I am interested in GUI based wizards for handling geoprocessing with libraries like gdal and ogr. Any tips?


Saving Interpolation data in a shapefile as an alternate to Raster Output

Is there a possibility that we store our interpolation data in a shapefile instead of a raster image, so that the shapefile can be used to display the Interpolation as a replacement to the RASTER OUTPUT. Also I intend to use such a shapefile in a HTML webpage with hover effects, so that if a hovers the mouse over some points, they can see the value at the location in a pop-up.


OpenLayer Maps not zooming

Initially the map zoom is set to 8, and that's the only place in the code i can change the zoom, any other location doesn't work. I have a function that draws a marker and a route on the map (works just fine) but i can't change the zoom on that part. When i try to change the zoom using `map.zoomTo(n)` it just zooms to maximum and `map.setCenter(new OpenLayers.LonLat(lon, lat), zoom)` also doesn't work. I am trying to zoom the map after i have redrawn poly layer and marker layer. Does anyone have any idea what could be done here?


Trial Version ArcEditor 10.0 Maplex question

I have the Trial Version for ArcEditor 10.0 with the extensions. Using Maplex/Labeling Manager, after clicking Property, I go to "FITTING STRATEGY". all the lables (stack, overrun, reduce font size, etc..) are all ghosted. I am not given an option to check/uncheck or make changes to anything on that page, nor LABEL POSITION / CONFLICT RESOLUTION pages. While in Trial, are you restricted what you can do? I thought you were able to use the full Maplex extension, or am I mistaken? What are my options, if any to make these areas workable? Thanks, Z ![enter image description here](http://i.stack.imgur.com/Jb3af.jpg) ![enter image description here](http://i.stack.imgur.com/aJEGY.jpg)


how to start qgis with --configpath option on windows

how can I start QGIS on Windows 7 with the --configpath option? I'm used to work with Linux where this procedure is very very simple (just write in a terminal **qgis --configpath** ) with all the options. Someone knows how to do the same thing on windows?


Path Distance calculation limited to downhill only

I am trying to calculate a path between two points with a relative elevation difference of about 100m and 3km separation. The points represent the start and end of a gravity fed waste water pipe. Essentially the path should follow the contours and have a consistent gradual fall. I have run the Esri Path Distance tool using Slope and DEM as cost and input surface (respectively) with the DEM as the vertical input as well using default values. The resultant path minimises slope but I also want it to only travel downhill. Can someone explain how to control the vertical factor parameters to only calculate a downhill path. The Esri help has not really described this.


Statistical Analysis of Shapefile Proximity within Raster Image - ArcMap10.1

I have a raster image that has been reclassified into 15 equally incremented values/levels. This raster represents "sinkhole risk" in a certain county. So level 1 is the lowest risk region and level 15 is the highest risk region. I also have the shapefile data for the actual sinkholes existing in that county. My hope is to come up with some method of statistically analyzing the data by answering the following question: **"What percentage of the sinkholes are in the risk regions 1-3?"** or **"What percentage of the sinkholes are in the risk regions 4-6?"** etc... Any suggestions on a tool to use or a general approach? I just need some guidance in the general right direction! Thank you. PS: I posted a photo below for more clarity if it helps! The photo has been simplified into 3 values/levels, but my real map has 15 values. ![Photo Explanation](http://i.stack.imgur.com/beoOr.jpg)


Inline variable substitution in ModelBuilder iteration to Make NetCDF Raster Layers?

I'm trying to convert lots of annual NetCDF files into rasters using ModelBuilder. I'm iterating over several years so the filename changes every year which works fine (using the inline substitution like %year%). Unfortunately, the variable name within the NetCDF files changes every year as well. So when I'm changing the variable name inside the "Make NetCDF Raster Layer" tool to "something_%year%" the tool won't accept it. Is there an easy workaround?


Create Fantasy Map

I'm currently looking to create a Fantasy map and share it with my friends. We have the map designed already in flash to certain degrees of detail and able to export as a vector. We'd like it to be usable as something like Google Maps, able to zoom in to show greater levels of details. Searching for locations...etc would be an added bonus. We've been looking into Leaflet and Open Street Maps, but seem to be a little confused and overwhelmed. Can anyone suggest some simple to use software where we may be able to import our vectors and build upon them? Online or offline.


Analysis of positioning polygon on polygon using ArcGIS

![enter image description here](http://i.stack.imgur.com/CmJWI.jpg)This might be a really complex analysis and I'm not entirely sure how possible it is but it is worth a shot. Is it possible to do analysis of positioning between either a "polygon and polygon" or "point and polygon"? Basically I have a polygon that delineates a boundary and I have another polygon that shows a building footprint and I want to know whether or not I can fit another building within this boundary. Of course the position of the existing structure plays a critical role in determining the possibility of adding a second structure. I have my doubts about doing this analysis as there could be any number of complex variables in determining whether or not a second structure would fit (ie. shape and size of the first building, shape and size of the second building, shape and size of the boundary etc...) Complete shot in the dark but wondering if anyone has ever done anything like this?


Are there better climate data than Worldclim for Europe?

I am looking for the best available climate data in Europe. Worldclim seems to be the best source of information for global data, but I am wondering if there is a higher resolution (better than 1km cell size and/or 1 month temporal resolution) product in Europe. (I've seen on this site that PRISM is great for the USA, but I couldn't find anything similar for Europe). EDIT : my purpose is to use the dataset for ecological modeling. EDIT: if you use another dataset, I am very interested in knowing why it is used instead of worldclim.


How to get Network Analyst data from script file back into Model Builder?

I am running a script which takes some input parameters and runs some Network Analyst processes. The script is imported into a toolbox and is setup in Model Builder with the necessary input and output parameters. My problem is that I am struggling to understand how to get the resulting Solve operation back out of the script and into Model Builder for subsequent work. Here is my script:               # Import arcpy module     import arcpy          # Check out any necessary licenses     arcpy.CheckOutExtension("Network")          # Script arguments     Network_Dataset = arcpy.GetParameterAsText(0) # Network dataset layer          Input_Locations = arcpy.GetParameterAsText(1) # Layer file          Output_Path = arcpy.GetParameterAsText(2) # Output path          # Local variables:     #Solve_Succeeded = Network_Analyst_Layer__3_     #Output_derived_feature_class = Network_Dataset          arcpy.AddMessage("Creating the Origin Destination layer...")          # Process: Make OD Cost Matrix Layer     arcpy.MakeODCostMatrixLayer_na(Network_Dataset, "odMatrix", "Length", "100", "", "", "ALLOW_UTURNS", "", "", "", "NO_LINES", "")          arcpy.AddMessage("Calculating Locations. This may take a while...")          # Process: Calculate Locations     arcpy.CalculateLocations_na(Input_Locations, Network_Dataset, "500 Meters", [["Streets","SHAPE"],["Streets_ND_Junctions","NONE"]] , "MATCH_TO_CLOSEST", "", "", "", "", "", "", "", "", "", "INCLUDE", "")          arcpy.AddMessage("Adding locations to Origin Destination layer...")          # Process: Add Locations to Origins Sub-Layer     arcpy.AddLocations_na("odMatrix", "Origins", Input_Locations, "", "500 Meters", "", "", "MATCH_TO_CLOSEST", "APPEND", "", "", "INCLUDE", "")          # Process: Add Locations to Destinations Sub-Layer     arcpy.AddLocations_na("odMatrix", "Destinations", Input_Locations, "", "500 Meters", "", "", "MATCH_TO_CLOSEST", "APPEND", "", "", "INCLUDE", "")          arcpy.AddMessage("Starting Network Analysis solve...")          # Process: Solve     arcpy.Solve_na("odMatrix", "SKIP", "TERMINATE", "")          arcpy.AddMessage("Network analysis completed.")          # How to get this out of script???          # Return as model output parameter     arcpy.SetParameterAsText(3, outputTable)


cql_filter for fid

when i run this code:               http://****/geoserver/Kaski/ows?service=WFS&version=1.0.0&request=GetFeature&typeName=Kaski:PIN&maxFeatures=50&outputFormat=json&cql_filter=(fid=97)      i got this error :               <ServiceExceptionReport xmlns="http://www.opengis.net/ogc" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" version="1.2.0" xsi:schemaLocation="http://www.opengis.net/ogc http://schemas.opengis.net/wfs/1.0.0/OGC-exception.xsd">     <ServiceException code="InvalidParameterValue">Illegal property name: fid</ServiceException>     </ServiceExceptionReport>      but i need to filter with fid. How can i do that ?


NetCDF file size calculation

I'm planning to use the NetCDF file format for storing my time-series observation data on a server. I would like to use the "NetCDF-CF" conventions. Before I start writing any scripts for importing data from my SQL Server database to NetCDF, I'd like to make an estimate of the disk space requirements for the NetCDF file. My data looks like this   * 1000 stations (each station has ID, Latitude, Longitude, Elevation)   * 10 quantities (for example, air temperature, air pressure,   precipitation..)   * 10 years of data with hourly time resolution, that is 10 years*365 days*24 hours = 87600 hours per station and quantity   * the observation data value fits in 4-byte integer (if it's a decimal number then I can scale the value to integer.) So I need 10 years*365 days*24 hours*4 bytes = 350400 bytes per station and quantity.   * I presume that I also need to store the "time" information somewhere in a suitable format. Usually 8 bytes is enough to store a time value, so I need 10 years*365 days*24 hours*8 bytes = 700800 bytes for the time column. Now calculating the NetCDF file size: ( [1000 stations] * [10 quantities] * [87600 hours] * [4 bytes per value]) + [700800 bytes for time] = 3504700800 bytes = 3.264 Gigabytes total Is this roughly correct? Or does NetCDF have any significant "storage overhead" like most relational databases do (for example in Postgres there is a 23 byte overhead per data row or in MSSQL it's 7 extra bytes per data row)?


FME - Can I write back into my original file geodatabase

I have a gas pipeline network which I have stored in an File Geodatabase as polylines. I've created an FME workspace to split the pipelines at fittings (valves, reducers, tees, etc) which are stored as point features. This effectively creates smaller line segments when a fitting overlaps a pipe. My question is: can I have FME remove my original pipes and replace them with my split lines? Or do I have to write to another geodatabase and then cut/paste back into the original file? Please note, I am using a tester to query out only a certain area, so I don't want to replace all of the pipes, just the area I'm working on. I'm relatively new to FME and appreciate any input. Thanks, Rob


How do I save an ArcMap map document using VB.Net?

I am trying to save the map document that is currently open. How do I write the code? I have tried many suggestions that I found using Google searches, but nothing works. Currently I have: Dim clickedDoc As IMapDocument = My.ArcMap.Application.Document   clickedDoc.Save(True, False) which causes ArcMap to crash. I did get it to work with some other code, but it wouldn't actually save. Also, this should be a background job. So no save dialog should come up.


Interactive geographic region definitions

## Background An organization has regions within a province, but has no list of coordinates that define region boundaries. I would like to draw custom polygon areas on a map using a web-based and user-friendly interface (such as Google My Maps or Scribble Maps). The regions are for data aggregation. ## Technical Maptitude and OpenGeo are being considered for the final GIS application. The database will likely be PostgreSQL running PostGIS. ## Question What open-source tool would you use to define the coordinates by drawing polygons on an interactive map? That is:               Browser -> Map -> Polygons (shared boundaries) -> Export      If the export step (to PostGIS) was performed automatically, so much the better. A cross-browser solution (IE8, FF, Chrome, Safari) with a working demo would be most helpful. ## Related Links   * http://stackoverflow.com/questions/2325670/draw-polygon-using-mouse-on-google-maps Thank you!


Issues importing feature class to SDE SQL Server 2012

I have received an error indicating there could be issues with my geometry when attempting to import a feature class into SDE. I ran the check geometry tool. No issues.   I have removed all fields except for my lat/lng, which is a double type, my geometry type, which is point, and objectID.   When doing this, the feature class imports without error. What could be the issue? My first guess is having to do with Nulls? When an objectID doesn't allow null values, does this mean that the feature should be free from nulls in other fields also? I am using ArcGIS 10.2 and sql server 2012 for my SDE connection.


Why is ArcMap leaving connections "<IDLE> in transaction"?

The TL;DR version: When you use ArcMap Query Layers to access a PostgreSQL table, ArcMap's database connections are almost always sitting idle in unclosed transactions. Why would this be? I have PostgreSQL 9.1.x (with PostGIS 2.0.x) with some existing tables and data. I successfully use ArcMap (10.1 SP1) to connect to one of those tables (using a Query Layer) and render dots on my map. The ArcGIS enterprise geodatabase tables are not installed. One day I brought up pgAdmin's server status window while ArcMap was running (and nothing else was accessing the database). As expected, there are two connections whose application name is "`ArcSDE(sv=nnnn;sn={xxx})`". But I was surprised to see that both of those connections are listed as "<IDLE> in transaction" in pgAdmin's "Query" column. That seemed fishy. As a rule of thumb, it's generally bad to leave lingering open transactions. And as a programmer, I've witnessed first-hand the dangers of improper failure handling that leaves unclosed transactions. So I closed everything down, restarted PostgreSQL, then started ArcMap again. Same result: two connections, both of which are "<IDLE> in transaction". I cranked up the PostgreSQL statement logging and watched as ArcMap did its magic. Sure enough, the very last statement in the log (once the ArcMap window finishes rendering) is `BEGIN`, with no corresponding `COMMIT` or `ROLLBACK`. And then every time you pan or zoom, the very first thing it does is `COMMIT`. Then it gets all the data and ends its rendering with an un-closed `BEGIN` again. At first glance, it appears benign. Those two connections hold no locks other than the locks on their own virtual transaction IDs. But it sure seems odd. Any ideas what legitimate reason ArcMap might have for behaving this way?


Shapefile to Network xml

I have a shapefile representing the US freight highway network available from Oak Ridge National Lab. I'm trying to prepare this as a `.xml` network for processing with MATSim. The output structure is very simple, something like               <?xml version="1.0" encoding="UTF-8"?>     <!DOCTYPE network SYSTEM "http://www.matsim.org/files/dtd/network_v1.dtd">     <network name="VISUM export national network 2007-11-28">          <!-- ====================================================================== -->              <nodes>             <node id="1000" x="730065.3125" y="220415.9531" type="2" origid="1000" />             <node id="1001" x="731010.5" y="220146.2969" type="2" origid="1001" />             .....         </nodes>          <!-- ====================================================================== -->              <links capperiod="01:00:00" effectivecellsize="7.5" effectivelanewidth="3.75">             <link id="100365" from="226" to="227" length="921.0" freespeed="33.3333333333333" capacity="5600.0" permlanes="2.0" oneway="1" modes="car" origid="183" type="10" />         ...         </links>     <!-- ====================================================================== -->          </network>      Most users of MATSim pull data from OpenStreetMap using `osmosis`, but I have to use _this_ network for a variety of reasons. I have seen lots of suggestions in different places as to what I might try. Some of these include:   1. `v.net` in GRASS   2. `shp2pgsql` suggested by pgRouting tutorial   3. manually with pgRouting, though this advice may be outdated?   4. writing code through the `sp` and `rgdal` libraries in R to build a node and link structure.   5. find someone with a TransCAD license to export the network ORNL provides along with the shapefile, though I don't know TransCAD exports to `.xml`. It seems that with options 1 - 3, I'll then have to figure out a way to write the `.xml` from the database. Option 4 seems difficult, but mainly in terms of my time. I'm happy to learn something new, but I don't want to invest a day towards option 1 only to find that option 3 would have been better. Or is there a simpler way? * * * Post-script: Someone asked the MATSim developers if they had a tool to do the `shp -> xml` conversion, but the answer was that shapefiles are too wild to build one.


SearchCursor (arcpy.da), accessing data on joined tables

Can you access joined tables through a da.SearchCursor? I've found documentation to do this with SearchCursor (arcpy), but nothing for the 10.1 data access module.


Geoserver, openlayers and WMS

I have millions of seismic lines to publish on a map (3,000,000+) so WFS is not really an option at a large scale. For this reason I will have to use WMS. My question is, is it possible using Geoserver as the engine and OpenLayers (latest versions) to present a map and allow the end user to click on an item and highlight the item? Basic map functionality I would have thought. There is a demo on the OpenLayers website but it appears this is the only example that does not work :( http://openlayers.org/dev/examples/getfeatureinfo-control.html Thanks, Chris


Creating perimeter around multiple polygons using ArcGIS for Desktop?

This may have a simple answer that I'm just missing. I have multiple polygons in one shapefile I would like to make a larger Polygon around all of them. Is there a tool that is able to do this? Or even a Polyline would work if that is not possible, I can always do the Feature to Polygon tool after. The images below better explain what I am looking to do (second is just a quick polygon I drew in a few seconds). It is in order to show different sales areas for my company. ![Polys](http://i.stack.imgur.com/PwhC9.jpg) ![enter image description here](http://i.stack.imgur.com/vsZTu.jpg)


How do I perform a 3 parameter geotransformation and projection in Sql Server 2012?

I have a table with latitude longitude (NAD27) columns in it. I compute two other columns, X and Y, representing Web Mercator (WGS84) location. Currently I'm using a Arcmap to do this, by applying the recommended geotransformation for the study area - the 3 parameter (geocentric) geotransformation \- to go from NAD27 to WGS84. I would like to do this entirely within Sql Server 2012. **From what I can tell, Sql Server does not support datum transformations out of the box.** Does anyone know of a Sql library that supports this geotransformation? I would like to simply use the same coefficients in Sql that I'm currently using in Arcmap. I also need to project from WGS84 lat/long into web mercator. I see this formula implemented in javascript, but if someone has a Sql stored procedure that does this, it would be great.


varying extent of each data driven page

I would like to export data driven pages, each with a different extent. I understand that within the Setup for the data driven pages, there is an option 'Extent' and one can choose 'Data Driven Scale'. I'd expect that if I select this and choose a column from the table of the data driven page layer, the extent would change according to the value. Apparantly this doesn't happen. How does this option work exactly? Or is there an other way for varying the extent of each data driven page?


OWSLib and python3

Trying to install OWSlib-master (and 0.7.1) from https://github.com/geopython/OWSLib yields problems that ultimately seem to be related to python3               >c:\Python32\python.exe setup.py install     running install     running bdist_egg     running egg_info     writing requirements to OWSLib.egg-info\requires.txt     writing OWSLib.egg-info\PKG-INFO     Traceback (most recent call last):       File "setup.py", line 34, in <module>         packages      = find_packages(),       File "c:\Python32\lib\distutils\core.py", line 148, in setup         dist.run_commands()       File "c:\Python32\lib\distutils\dist.py", line 917, in run_commands         self.run_command(cmd)       [...]         lines = header.split('\n')     TypeError: Type str doesn't support the buffer API      Are those errors errors of my (ActivePython 3.2) installation, or are these errors of OWSLib being incopatible with python 3? Has anyone succeded in running OWSlib on python 3?


How do I get a decent easy to use watermark from an image in Quantum?

Our users need to be able to put a graphic watermark onto maps they create in the print composer. Currently we get them to georeference the image as a raster and then create a transparency. This is quite complex for some of my beginners. Anybody got any suggestions?


MS Access and geocoding in ArcGIS

We are a small non-profit. We use a fairly complex MS Access application to track everything. We would very much like to geocode addresses in our application and display the results on a map. In other words, we would like to input an address and have it display on a map. We would need to use the GIS data provided by our municipality. They use ArcGIS and the MXD files contain a fair amount of ancillary data for each point. We have access to ArcGIS 9.x. But I am unable to get the map control working in our Access application, and apparently ESRI is eliminating that control anyway. I tried the following already:   * Google Maps static API. Works pretty well, but I have no way of geocoding against the local municipal data.    * MapWinGIS makes great maps based on ArcGIS file, but I have no way to geocode against the local data, as far as I can tell.   * Export MXD to KML. Works OK, but no labels/other ancillary data.   * Export points in ArcGIS and import to Access table. Sloppy, doesn't bode well for future updates to the GIS data provided by the municipality, and have no way of dealing with intersections. I would be willing to import the ArcGIS data into a different GIS product, if it would help (we could presumably just import new versions of the data in the same way.) Any suggestions greatly appreciated.


Strange Behavior for TOCControl PopupMenu

The sample "TOCControlContextMenu" demonstrates how to display a context menu from the AxTOCControl. When I run this sample and right-click on the TOC control the context menu pops up and remains displayed until I select a menu item. In my application the right-click displays the context menu but as soon as I release the right-mouse-button the first item in the menu is selected. I can see no substantial differences in the code. Has anyone observered (and possibility overcome) this behavior before?


Draw elevation from mouse move

I recently tried this example but unfortunately It's not working, I don't know why (I've tried as same as possible) My output is I want to use WMS data from geoserver and applied it in CesiumJs or Openwebglobe (3D platform) Does anyone know how to make a similiar example like this? Thanks Best Regards


network based interpolation methods, which one could I use?

I need to interpolate traffic data consisting in raw number of vehicles passing across certain points, where counters are. It is evident I haven't a uniform distribution of points and values about traffic. I would like to ask which methods I could use to interpolate this network deducing traffic for each location of the overall network. I tried to look on the most used desktop GIS (QGIS, ArcGIS) but it seems there is not a'ready made' solution at the moment. Is there an algorithm that I can implement in Python to achieve this result? Thanks!


Customized PDF Report of Each Manhole from a Database

I need a way to generate individual PDF reports for all the data collected for each manhole in a database. I would end up with hundreds of reports, one for each manhole in the feature class. The report needs to selected attribute data such as manhole type in a chart along with 1-3 pictures and possibly a small map. It doesn't look like the reporting tool built into ArcMap can make something this complicated. Any help would be greatly appreciated!


Best place to download or purchase shapefiles for Louisiana Parish?

I am looking for shapefiles to either download or purchase for the BEAUREGARD Parish in Louisiana. I have been able to find files for the Townships and Sections but I am in desperate need of a shapefile with TRACTS. If anyone know where I could obtain something like this it would be greatly appreciated. Thanks in advance for any help that you can provide.


Failed to convert polygon to raster

I have two steps here. First select layer by attributes. Second convert polygon to raster. It failed when it goes to second step and kept giving me error 999999: Error executing function. I pasted my codes below. Is there anyone who has the same problem. I am now using Arcgis 10.1 and has SP1 installed.               exec('''     import arcpy     from arcpy import env     arcpy.env.overwriteOutput = True     arcpy.CheckOutExtension("spatial")     import string     env.workspace = "G:\\countoverlap"               # Local variables:     FRAP_Fire111 = "G:\\countoverlap\\FRAP_Fire111.shp"          fldName = 'Year'     fcName = 'FRAP_Fire111.shp'     #set creates a unique value iterator from the value field     yearList = list(set([row.getValue(fldName) for row in arcpy.SearchCursor(fcName)]))     print(yearList)           ''')               env.workspace = "G:\\countoverlap"     env.scratchWorkspace = "G:\\countoverlap\\overlap.gdb"     arcpy.MakeFeatureLayer_management(FRAP_Fire111, "lyr")          for year in yearList[1:len(yearList)]:              print(str(year))              # Process: Select Layer By Attribute         arcpy.SelectLayerByAttribute_management("lyr","NEW_SELECTION", '"Year" = year' )         print('Layer_Seclected')              # Process: Polygon to Raster         # Set local variables         inFeatures = "lyr"         valField = "Year"         outRaster = "str(env.workspace)+'\\'+str(year)+'.tif'"         assignmentType = "MAXIMUM_AREA"         priorityField = "NONE"         cellSize = 1080         # Execute PolygonToRaster         arcpy.PolygonToRaster_conversion(inFeatures, valField, outRaster, assignmentType, priorityField, cellSize)         print('Raster_Created')


Georeference to UTM projection using openlayers plugin in QGIS

I have some aerial imagery that I need to georeference using QGIS into a UTM projection. In order to do this I need a basemap layer for the reference points. The openlayers plugin could work great for this but it's in the web Mercator projection. I attempted to georeference to the Web Mercator projection and then use the warp tool to transform to NAD 83 UTM but it doesn't seem to transform correctly. Should my process here work and I'm just missing a setting or step, or should I be going about this a different route all together?


Transforming WFS Projection

I am trying to overlay features from a WFS using OpenLayers, my geometries in database are in EPSG:27700 and my map's projection is EPSG:3857. My mapfile returns data in EPSG:27700. Is there any way I can transform points into the map projection? Here is my code which is not working.               var wfs = new OpenLayers.Layer.Vector(                 'layername',                 {                     strategies : [new OpenLayers.Strategy.BBOX()],                     protocol : new OpenLayers.Protocol.WFS({                         url : 'proxy.php?url=xxx',                         featureType : 'typename',                         "geometryName": "msGeometry",                         "featurePrefix": "ms",                         "srsName": 'EPSG:3857',                         "maxFeatures": 1000,                         "version": "1.0.0"                     })                                     }         );


How to create a schematic view of a line feature?

DEAR PEOPLE FROM THE FUTURE: Here's what we've figured out so far… I have a line feature (polyline with curves) in ArcGIS Desktop, and I want to create a linear representation (straigth line) of the feature, furthermore i want to be able to locate events in the linear representation. this is somethig like i want ![enter image description here](http://i.stack.imgur.com/Gtw4A.jpg) Making some research i found something similar called "column view" unluckly my pc reboot an does not saved the history of the browser....so I lost the document. (I have been searching for it again, but i have not been able to find it)


Duplication of WMS Imagery

I am trying to get a state served Imagery WMS layer on a leaflet example and I am getting duplication of the image itself. Any ideas would be great. Below is the URL of the imagery: http://demac.udel.edu/data/aerial- photography#aerial-mapservices


how to create a shapefile containing attributes from a list?

I have a list of features and each feature is a dict in that list. following you can see the structure: [ {"FID": 0, "Shape": polyline, "Nodes": [0,1]}, {"FID": 1, "Shape": polyline, "Nodes": [0,2]},...] Now, I want to create an output shapefile containing all the attributes here. I should mention that I can create the shapefile using the below code but it does not add the attributes, which means that afterward, I need to iterate over the features and update them.               geom_list=[]     for feature in list_of_dicts:         geom_list.append(feature[u"Shape"])          outfile = "D:/address/to/the/result_features.shp"     result = arcpy.CopyFeatures_management(geom_list, outfile)      this is a hassle and in huge dataset can be really slow. is there anyway to create a shapefile with a list of attributes?


Conflict between ArcPy and GDAL/OGR?

# TL;DR I can't get an installation free version of GDAL/OGR to work from within ArcMap 10.0. I get an `ImportError: DLL load failed` on `_mod = imp.load_module('_gdal', fp, pathname, description)`. Below, more details on the context and a procedure on how to reproduce the error. * * * Someone in the comments specified that arcpy and GDAL work fine together on his machine. What could be the differences between our configurations (software installed, environment variables, ...)? # Description of my issue ## What I am trying to do... I want to deploy an ArcMap 10.0 ToolBox with scripts that use external libraries. I would like to avoid having any manipulation done on the user side: ideally, no additional installation/manipulation should be required for the user. For that reason, I embed the libraries files inside the ToolBox. ## What my problem is... When I try to use gdal/ogr from osgeo along with arcpy, I get the following error:               Traceback (most recent call last):       File "path_to_script\debugDepl.py", line 32, in <module>         from osgeo import ogr       File "path_to_script\osgeo\__init__.py", line 30, in <module>         _gdal = swig_import_helper()       File "path_to_script\osgeo\__init__.py", line 26, in swig_import_helper         _mod = imp.load_module('_gdal', fp, pathname, description)     ImportError: DLL load failed: The specified procedure could not be found.      There are some binding errors that I don't understand. It seems that there is a conflict between arcpy and gdal/ogr (more information below) Related:   * GDAL bindings problem with ArcMap 10.1 Python tool   * gdal ImportError in python on Windows   * Conflict arcpy (Arc 10) and gdal?   * Running GDAL and arcpy together   * How to run OGR in ArcGIS Scripts? * * * # How to reproduce the error If you can help me achieve want I want, it might be useful if I explain every step that lead to this error. I tried to detail as much as I could so that you could reproduce it. **1/** On a Windows 7 64bits machine, Install ArcGIS Desktop 10.0 (32 bits) **2/** Download GDAL‑1.11.0.win32‑py2.6.exe Windows binary **3/** Install that library with your Python installation **4/** Move the associated library files (`path_to_python_install\Python26\ArcGIS10.0\Lib\site-packages`) to another new folder (osgeo folder + several Python files) **5/** Create a Python script in that folder and edit it in your favorite Python editor and check that you can import gdal/ogr               from osgeo import gdal     from osgeo import ogr      **6/** Try importing arcpy from the script before and after importing gdal/ogr               from osgeo import ogr     from osgeo import gdal      import arcpy      That should work, now try the other way around               import arcpy     # Maybe some environment variable settings     from osgeo import ogr     from osgeo import gdal      And that version should create the error. **7/** Create a toolbox (or call the script from within ArcMap) with the version that worked standalone. And try to run it within ArcMap (should fail). ## Notes on the steps **1/** I've found some official Python bindings but only for ArcMap 10.1 **2/** This distribution includes a complete GDAL installation. It is a 32bit version for my 32bit version of Python26. **4/** As specified in these recommendations for scripts deployment, I should copy the library files directly. I want to avoid the user to have to install things as much as possible **5/** This may (or may not) require some environment variable settings:               import os     environList = os.environ['PATH'].split(';')     environList.insert(0, r'path_to_folder\osgeo')     os.environ['PATH'] = ';'.join(environList)     os.environ['GDAL_DATA'] = r'path_to_folder\osgeo\data\gdal'      GDAL Documentation specifies that we should:   * Add the installation directory bin folder to your system PATH, remember to put a semicolon in front of it before you add to the existing path.   * Create a new user or system variable with the data folder from your installation. Developers on this blog page seemed to have that error when the environment variables were incorrect. **6/** It seems that for some reason if you import arcpy first and ogr/gdal second, gdal is no longer able to find the dll's that it requires. This could have been due to some environment variable settings but the modification I proposed at step 5) didn’t change anything. **7/** But then my guess is that arcpy is automatically imported (first) so we’re back to the 2nd version of 6) ## Some more information about my attempt to understand what is going on From PyScripter:               >>> from osgeo import gdal     >>> gdal.__file__     'C:\\Python26\\ArcGIS10.0\\lib\\site-packages\\osgeo\\gdal.pyc'     >>> import arcpy     >>> arcpy.__file__     'C:\\Program Files (x86)\\ArcGIS\\Desktop10.0\\arcpy\\arcpy\\__init__.pyc'      Alright, the library files were found where they are supposed to be. But then…               >>> import arcpy     >>> from osgeo import gdal     Traceback (most recent call last):       File "<interactive input>", line 1, in <module>       File "C:\Python26\ArcGIS10.0\lib\site-packages\osgeo\__init__.py", line 29, in <module>         _gdal = swig_import_helper()       File "C:\Python26\ArcGIS10.0\lib\site-packages\osgeo\__init__.py", line 25, in swig_import_helper         _mod = imp.load_module('_gdal', fp, pathname, description)     ImportError: DLL load failed: The specified procedure could not be found.      Maybe something happened to the environment variables or to the python path when importing Arcpy? Well, it doesn't seem like it...               >>> import os     >>> import sys     >>> envBefore = os.environ     >>> pyPathBefore = sys.path     >>> import arcpy     >>> envAfter = os.environ     >>> pyPathAfter = sys.path     >>> envBefore == envAfter     True     >>> pyPathBefore == pyPathAfter     True      Thank you for reaching the end of that long post.


Is it possible to use comments in QGIS Query Builder?

I am trying different queries in the Query Builder. Is there a way to comment the query? I tried with `/* */` `#` `--` without success : > An error occurred when executing the query. The data provider said: OGR[3] > error 1: SQL Expression Parsing Error: syntax error **EDIT** : Thanks to the answer and comments of @Joseph,   I have made a feature request of this question: http://hub.qgis.org/issues/10914


ArcGIS REST vs. SOAP API

When is it appropriate to use the ArcGIS Server REST API vs. the SOAP API and vice-versa? What do you see as the advantages of one over the other? For example, the SOAP service can be consumed as a Service Reference to integrate into Visual Studio. Is there anything available that would give you that level of integration with REST? More info: ArcGIS GIS Services


Pulling SDE feature classes from current map to fieldmapping as DETable

I'm working with SDE layers in arcmap. trying to avoid pointing to the sde file itself. I've tried workspace.Pathname, it was empty for sde remote workspace. I've also tried IName from workspace, empty again. Any ideas on what I should be doing differently? BTW, it works just fine for GDBs.                       Dim gputilities As ESRI.ArcGIS.Geoprocessing.IGPUtilities = New ESRI.ArcGIS.Geoprocessing.GPUtilitiesClass()             'Dim inputTableA As IDETable = CType(gputilities.MakeDataElement(pWorkspace.PathName + "\Weld", Nothing, Nothing), IDETable)             Dim inputtables As IArray = New ArrayClass()             Dim gpValueTableObj As ESRI.ArcGIS.Geoprocessing.GpValueTableObject = New ESRI.ArcGIS.Geoprocessing.GpValueTableObject             gpValueTableObj.SetColumns(1)             Dim pdet As IDETable = New DETable             For Each pLayer In checkedLayers                 Try                     gpValueTableObj.AddRow(pLayer)                     Dim workspaceFactory As SdeWorkspaceFactory = New SdeWorkspaceFactoryClass()                     Dim workspaceName As IWorkspaceName = workspaceFactory.Create(Nothing, pLayer, pWorkspace.ConnectionProperties, 0)                     Dim pName As IName = TryCast(workspaceName, IName)                     inputtables.Add(CType(gputilities.MakeDataElementFromNameObject(pName), IDETable))                      Catch ex As Exception                 End Try                 'inputtables.Add(CType(gputilities.MakeDataElement(weldLayer.Name, Nothing, Nothing), IDETable))             Next             fieldmapping = New ESRI.ArcGIS.Geoprocessing.GPFieldMappingClass()             fieldmapping.Initialize(inputtables, Nothing)


ArcGIS 10 Linux experience

We have ArcGIS 10 installed on a nice clean Linux box. We have developed scripts to run on the Windows desktop development machines, and modified them to run on Linux. Now, if we run the scripts on the windows boxes, they 'perform' within ioperational tolerances and perform correctlty. However, once we run them on the Linux envronment, they seem to 'cock about'. I am led to believe that ArcPy is synchronous and that calls to various objects have to finish computing, before the next call toprocessing can take place. It seems not. A call to aggregate points to a polygon, then to buffer the polygon, see's an empty buffer returned, when the polygon exists. There is some latency between the Dev box and the ArcSDE instance, but surely if this is truly synchronous, it would only move on if the work had completed? Does anyone have any experience of this?


Assigning points to their respective Reaches using Stream Gradient?

This is a continuation of the question I initially posed here. Also, I've attached a shapefile on my Drop Box account here. I'm trying to break a stream network up into 100m reaches, assigning elevation values pulled from a DEM, to the start and end of each segment, the goal being to determine gradient for each segment. Process-wise, I'm not approaching this correctly. So far, I haven't been able to those z-values to the start and end nodes. What I've been doing is to segment the stream network, and using XTools, convert Features to Points, again, at 100m intervals. They do line up with the segments, as you can see below. From there, I can easily use XTools to pull raster values, and assign them to the points. What I'm not able to do, is to assign those points back into the reaches. I can get one of them, sure, but not both. Thus, I can't calculate gradient. Again, any help would be MUCH appreciated, particularly any sort of working example. ![enter image description here](http://i.stack.imgur.com/irCsx.png)


Is the csv handling in QGIS bugged?

I am using QGIS 2.4.0 and sadly it seems that QGIS imported csv files are somewhat bugged. I cannot use imported csv files the same way as shapefiles. It seems that sometimes it is easier to import the file and then just save the csv file inside QGIS as a shapefile and import it again to make it work. The reasons I am saying this are the following:   * I cannot save the style files associates with csv files.   * In properties tab type and typename are always double and the lenght and precision are always 0   * Calculations are not possible since some columns are not seen as containing Real or Natural numbers. Once saved as shapefile I can suddenly do everything I need. So am I missing something or is there no other way to not to have this intermediary saving step?


Image classification using google satellite images

Currently our expert user classifies locations by viewing the Google earth satellite image of a location. They consider circular areas of a 5km radius at a set zoom level. They evaluate the image for density of buildings or trees in within the circle and assign the location a classification. We get data in the form:               location         | Classification     (lat-long)       |     ----------------------------------     123.123, 12.34   | B     122.123, 12.34   | A     121.123, 12.34   | C     120.123, 12.34   | A     119.123, 12.34   | A      I want to use this as training data to build a system that can automatically classify a location based on the satellite image. To train the system, I would need to extract relevant parameters from the satellite images. I guess the image parameters would need to be rotation independent. Possible parameters:   * total intensity for RGB values,    * histogram   * count of edges   * ?? I plan on using C# and available AI & image processing libraries (e.g. opencv, encog, AForge) to do implement the training and classification stages. Some of the image processing tutorials concentrate on object recognition but that's not really what I'm after (I think). It's more about grouping the images into different categories based on more general properties of the image. Can anyone advise on a) the feasibility of this plan and b) what parameters I should be extracting from the images? Thanks


How to display Google Maps as basemap in a Silverlight application?

I am writing a web mapping application and I use ArcGIS Server services. I want the application to use Google Maps as a base map. I am new to this job. I have searched a lot on the Internet. But I could not find anything useful. Is there anyone who could help me in this matter. Sample site are shown below. http://esup.qom.ir


How can I zoom to extenct according for dynamic WMS (coming from DB)?

I have a WMS layer that is coming from SQL View. I want to be able to zoom to extent of layer when clicking a button (view changes because data keep coming to the DB). How can I get the rectangle that currently contains the points and zoom to it ? (or an alternative is to make a polyline from the points and zoomtoexctent of it.. but its kinda overkill I think)


Tiles cover Polygon on zoom and move

i'm using OpenLayers.Geometry.Polygon.createRegularPolygon to generate some nice circles. After i draw them to the Map i'm Moving them and after some while they where hiding by some tiles. The circle zIndex is higher then the tiles zIndex so there is no problem (hopefully). I'm using Openlayers Version: Release 2.13 Here is some code i'm using:               var circle = OpenLayers.Geometry.Polygon.createRegularPolygon(         point,         rad,         200,         0     );     var style = {         fillColor: '#ea2828',         strokeColor: '#ea2828',         strokeOpacity: 1,         strokeWidth: 4,         graphicZIndex: 1099,         fillOpacity: 0     };     var featureInner = new OpenLayers.Feature.Vector(circle, null, style);     layer.addFeatures(featureInner);      when i want to set the new point im using:               layer.removeAllFeatures();      to remove the old circle and then i put it with the code from above again on the layer. Every time i zoom into the map where i can't see the circle and move to the circle postion i can't see the circle. The Magic at this point is that i can use the controlls hide layer and activate it again and wow the circle is back. I hope you have a idea for a fix. Best regards


How add vector style rule filter with different value in table

I have a simple wfs layer from geoserver, here is code:                 //wfs-t editable overlay       var wfs_layer = new OpenLayers.Layer.Vector("Vector sample", {       strategies: [new OpenLayers.Strategy.Fixed()],       //projection: new OpenLayers.Projection("EPSG:4326"),       protocol: new OpenLayers.Protocol.WFS({       version: "1.1.0",       // loading data through localhost url path       url: ".....",       featureNS :  ".....",       extractAttributes: true,       // layer name       featureType: "....",       geometryName: "geom",       schema: "localhost:8080/geoserver......"       })       });           //vector style:     var vector_style = new OpenLayers.Style({      'fillColor': '#669933',      'fillOpacity': .8,      'strokeColor': '#aaee77',      'strokeWidth': 3,      'pointRadius': 8      });          var vector_style_select = new OpenLayers.Style({      'fillColor': '#FF0000',      'fillOpacity': .10,     'graphicName': 'circle',      //'label': 'klik',     'pointRadius': 8,      'strokeColor': '#f86633',      'strokeDashstyle': 'solid',      'strokeWidth': 4      });          var vector_style_iz_table = new OpenLayers.Style({     'cursor': 'pointer',      'fillColor': '#669933',      'fillOpacity': .8,      'fontColor': '#343434',      'label': '${table1}',      'pointRadius': '8',      'strokeColor': '#aaee77',     'strokeWidth': 3,      'strokeDashstyle': 'solid',      'strokeWidth': '${4}'      });     var vector_style_map = new OpenLayers.StyleMap({      'default': vector_style_iz_table,     'select': vector_style_select      });     wfs_layer.styleMap = vector_style_map;      I have table from "table1" database, where is three different value: "yes", "no" and "null". Everything works fine, in this style I can view vector layer with green fillcolor with label "yes" or "no" from my "table1". I want 3 different style, for "yes"-> green fillcolor, for "no"-> red fillcolor and "null"-> black fillcolor. Is that posible having 3 styles for one vector layer? I think that I must write rule filter for all three different value? Thanks!


How to rotate backgroundgraphic only in openlayer

I know how to rotate image but if I want to keep the background image fixed what I should I do , can someone please let me know how to do that? Many thanks any comments is highly appreciated.


Removing Places from OSM

Doe is exists a way to remove all the places from the OSM map? Also like shops, bars, restorations, hotels etc. I wish to use a maps with a less info. Is it possible to do this from the original OSM server? Maybe like a option in URL or something else? I use Leaflet library for my project, maybe some option in it?


QGIS split shp into different layers

I'm building my first map, so please forgive me beginners mistakes.   I've downloaded QGIS 2.0.1 and following tutorial I was able to download data from http://sdi.geoportal.gov.pl/wfs_prg/wfservice.aspx?request=getfeature&version=1.1.0&service=wfs&typename=WOJEWODZTWA and present is as one layer in application. My SHP file has 16 objects (Polish provinces).   ![enter image description here](http://i.stack.imgur.com/ibXDV.png) My problem is that right now I have one layer (with one layer style) and I can show/hide all provinces at once, what I want to do is to have all provinces as layers (or groups) so I could show only some provinces and change their color. I know that I can select shape using select tool, then export selection into SHP file and import 16 layers into project.   This is fine when layer has 16 objects, but my second map has over 3000 districts and I must be able to hide/show/change color of them individually. Can this be done better/faster? **EDIT**   I managed to split my SHP into layers, this way I have 16 layers and I can add style to every layer individually. I can also hide layers using layers panel. ![enter image description here](http://i.stack.imgur.com/4Hngy.png) Can I add similar select to layer? So I can choose which object to show/hide?


Lambert conformal conic projection - Europa

I am trying to calculate map of Europa in Lambert projection. I am using formula from wikipedia (Lambert wiki), but calculation is verřy slow with comparison to Mercator. I know, there are sin, cos, cotg and pow to calculate (mercator has only one tan and ln), but the difference is almost 30 minutes (Mercator is done in about 3 minutes). Is there some kind of speed up, or similar projection, that can be used ? I also need to put result into image with width and hight. For Mercator, I am using this: http://stackoverflow.com/questions/5983099/converting-longitude- latitude-to-x-y-coordinate. Can it be altered to Lambert (or Lambert similar) projection ? Because if I use same approach only witch changed formulas for projection, Lambert is moved to upper part of image and I can see only part of Europe) **EDIT: Added code**               double longitude = lon * (MPI / 180.0);     double latitude = lat * (MPI / 180.0);          double lambertPhi = this->lambertF * pow(cot(0.25 * MPI + 0.5 * latitude), this->lambertN);          lonOut = lambertPhi * sinf(this->lambertN * longitude);     latOut = this->lambertPhi0 - lambertPhi * cosf(this->lambertN * longitude);      All **this->** variables are precomputed only once


How to keep track of edited location?

With large datasets, how do you guys keep track of where you have edited across the map? For instance, I am editing a city's waterline and water valve data. I estimate using a 1:400 average zoom level. This is a large city, and I can quickly lose track of where I have edited, especially after days of working on it. I'm using ArcMap 10.2 Basic with a file geodatabase. Is there a function to track the location and extent of where my window has been? Or is there a better way to track changes/edited location? If not, I might write script to make an initial polygon, or draw a feature, based on the current window extents, and then merge new extents to it as I pan across the map. Perhaps stop it from drawing if zoom is beyond 1:1000. What do you guys say?


How can I apply the same change to multiple map documents?

I have several map documents (around 30) that reference all the same files and display them all the same way. I created them by setting one map how I wanted it, then using "save as" to create a new map then adjusting the data frame extents in the layout. Of course now that I have them all created, I need to make the exact same change to all of them. Is there a more efficient process than just manually repeating the same steps in each document? If not, is there a better way I could have set up the project to facilitate changes like this? I'm using ArcMap 10.0 Thanks


imposm2 to imposm3

I am using imposm2 to import osm data to postgresql database but recently moved to imposm3 which is still in development phase I want to know is there an easy way to convert imposm2 mapping file (python) to imposm3 mapping file(json) ?


Run a ArcMap (10.0) Addin Button with VBA from a MS-Access form

I need a way to execute an ArcGis-Addin (vb.net) from a MS-Access database form (VBA). The problem is that there is no ESRI.ArcGIS.Desktop.Addin reference in VBA and the "pItem.Execute" method with the "pCmdBars.Find(pUID)" does not work with my Addin-Button. I want to open the ArcMap.mxd file and perform a selection that shows the spatial data of my MS-Access database query. The rest is not a problem (maybe!). Thanks, in advance....


How to display a single in a group of nodes with same attributes?

In QGIS, if I have some nodes with a name attribute is there a way to display only one node that would be in the middle of this group of nodes having the same name? For example, let's say i have 3 nodes with an attribute "ABC" I would like to display one node what would be in the middle of this group.


session id changing ; tomcat and apache web serevr

I'm using a WFS protocol on Vector layer. also I'm working on 2 servers, tomcat (to deploy GeoServer) and Apache web server(because i use php) I have an HTML page works on Tomcat and login.php and another php files on Apache web server the problem is the session id changes every request.. so i got the "undefined user id" error but when i put all my files on Apache.. i got the same session id


What to do with old GIS books/manuals/etc?

We're having an office clear-out and have come across a number of old GIS paper materials. Random examples (to give a taste - there are about 20-30 in all) include:   * Postcodes: The New Geography (Book - 1992)   * Intermap (Product Handbook and Quick Start Guide (2004))   * ARC Macro Language (Self-study Workbook 1993)   * Understanding GIS; The ARC/INFO Method (ESRI Book, 1990)   * Institute of British Geographers; Transactions New Series; Volume 16 Number 1 1991   * arc/info; Computer mapping and geographica information management system (Folder, probably training materials; Nov 1985!) So I'm wondering - what should we do with them? Is there anywhere that would take "historic" GIS materials (museum, library etc)? They're all a bit too specialised for a regular charity shop but it seems a shame to just dump them. (Note: I'm in the UK if that makes any difference)


How to get Quantum GIS display a line or polygon feature created by Python package pyshp

I am not sure if I have discovered a regression, or if I have not enough sleep. While I manage to create shapefiles with point geometries with pyshp, **line and polygon features don't show up in QGIS**. In ArcMap, however, they show up as expected. **Steps to reproduce**   1. Get pyshp (version 1.1.7)   2. Execute one or both of these example scripts, taken directly from the pyshp documentation: _Polygon_              import shapefile     w = shapefile.Writer(shapefile.POLYGON)     w.poly(parts=[[[1,5],[5,5],[5,1],[3,3],[1,1]]])     w.field('FIRST_FLD','C','40')     w.field('SECOND_FLD','C','40')     w.record('First','Polygon')     w.save('test-polygon')      _Polyline_              import shapefile     w = shapefile.Writer(shapefile.POLYLINE)     w.line(parts=[[[1,5],[5,5],[5,1],[3,3],[1,1]]])     w.field('FIRST_FLD','C','40')     w.field('SECOND_FLD','C','40')     w.record(FIRST_FLD='First', SECOND_FLD='Line')     w.save('test-line')        3. Open the resulting shapefile in an empty QGIS session (tested in 1.7.3, 1.8.0 and 1.9.0 Alpha) **Expected result**   A polygon or a line in the map window (at least after a "zoom to layer extent"). **What I get**   An empty map. Did I miss something? **Update** After contacting the package author, version 1.1.9 was mentioned on the homepage. After an upgrade, the issue vanished. _That's_ good customer support!


Convert .osm to .map for mapsforge

I have download latest version of osmosis from http://wiki.openstreetmap.org/wiki/Osmosis#Downloading and mapsforge-map- writer-0.3.0-jar-with-dependencies.jar from https://code.google.com/p/mapsforge/downloads/list. I copied mapsforge-map-writer-0.3.0-jar-with-dependencies.jar into 1- osmosis lib directory: osmosis-latest/lib/default 2- C:\Users\username\AppData\Roaming\Openstreetmap\Osmosis\Plugins 3- C:\Users\username\AppData\Local.openstreetmap\osmosis\plugins I downloaded iran.osm file from geofabrik.de and then I created a batch file as below:               osmosis --rx file=iran.osm --mw file=iran.map bbox=35.574683,51.086426,35.833402,51.664581 type=hd      and I got this error:               SEVERE: Execution aborted.     org.openstreetmap.osmosis.core.OsmosisRuntimeException: Cannot load JPF-plugin 'org.mapsforge.MapFileWriter' for extensionpoint 'Task' at org.openstreetmap.osmosis.core.TaskRegistrar.loadJPFPlugins(TaskRegistrar.java:186)         at org.openstreetmap.osmosis.core.TaskRegistrar.initialize(TaskRegistrar.java:88)         at org.openstreetmap.osmosis.core.Osmosis.run(Osmosis.java:81)         at org.openstreetmap.osmosis.core.Osmosis.main(Osmosis.java:37)         at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method)         at sun.reflect.NativeMethodAccessorImpl.invoke(Unknown Source)         at sun.reflect.DelegatingMethodAccessorImpl.invoke(Unknown Source)         at java.lang.reflect.Method.invoke(Unknown Source)         at org.codehaus.plexus.classworlds.launcher.Launcher.launchStandard(Launcher.java:329)         at org.codehaus.plexus.classworlds.launcher.Launcher.launch(Launcher.java:239)         at org.codehaus.plexus.classworlds.launcher.Launcher.mainWithExitCode(Launcher.java:409)         at org.codehaus.plexus.classworlds.launcher.Launcher.main(Launcher.java:352)         at org.codehaus.classworlds.Launcher.main(Launcher.java:47)     Caused by: org.java.plugin.PluginLifecycleException: plug-in org.mapsforge.map.mapfile-writer requires plug-in org.openstreetmap.osmosis.core.plugin.Core which is unknown or has incompatible version at org.java.plugin.standard.StandardPluginManager.checkPrerequisites(StandardPluginManager.java:460) at org.java.plugin.standard.StandardPluginManager.activatePlugin(StandardPluginManager.java:397)         at org.java.plugin.standard.StandardPluginManager.activatePlugin(StandardPluginManager.java:245)         at org.openstreetmap.osmosis.core.TaskRegistrar.loadJPFPlugins(TaskRegistrar.java:182)         ... 12 more      PS: I switched back to osmosis 0.36 and I got this error:               SEVERE: Thread for task 1-rx failed     org.openstreetmap.osmosis.core.OsmosisRuntimeException: Cannot represent 2147487728 as an integer.         at org.openstreetmap.osmosis.core.util.LongAsInt.longToInt(LongAsInt.java:33)         at org.openstreetmap.osmosis.core.domain.v0_6.CommonEntityData.init(CommonEntityData.java:142)         at org.openstreetmap.osmosis.core.domain.v0_6.CommonEntityData.<init>(CommonEntityData.java:74)         at org.openstreetmap.osmosis.core.domain.v0_6.Entity.<init>(Entity.java:61)         at org.openstreetmap.osmosis.core.domain.v0_6.Node.<init>(Node.java:71)         at org.openstreetmap.osmosis.xml.v0_6.impl.NodeElementProcessor.begin(NodeElementProcessor.java:87)         at org.openstreetmap.osmosis.xml.v0_6.impl.OsmHandler.startElement(OsmHandler.java:90)         at org.apache.xerces.parsers.AbstractSAXParser.startElement(Unknown Source)         at org.apache.xerces.parsers.AbstractXMLDocumentParser.emptyElement(Unknown Source)         at org.apache.xerces.impl.XMLDocumentFragmentScannerImpl.scanStartElement(Unknown Source)         at org.apache.xerces.impl.XMLDocumentFragmentScannerImpl$FragmentContentDispatcher.dispatch(Unknown Source)         at org.apache.xerces.impl.XMLDocumentFragmentScannerImpl.scanDocument(Unknown Source)         at org.apache.xerces.parsers.XML11Configuration.parse(Unknown Source)         at org.apache.xerces.parsers.XML11Configuration.parse(Unknown Source)         at org.apache.xerces.parsers.XMLParser.parse(Unknown Source)         at org.apache.xerces.parsers.AbstractSAXParser.parse(Unknown Source)         at org.apache.xerces.jaxp.SAXParserImpl$JAXPSAXParser.parse(Unknown Source)         at org.apache.xerces.jaxp.SAXParserImpl.parse(Unknown Source)         at javax.xml.parsers.SAXParser.parse(Unknown Source)         at org.openstreetmap.osmosis.xml.v0_6.XmlReader.run(XmlReader.java:108)         at java.lang.Thread.run(Unknown Source)


How to use SEXTANTE in a script outside QGIS?

I´m trying to do some geoprocessing using the QGIS Sextante plugin outside of the QGIS application (and integrate it with different other libraries like arcpy, etc.).               import sys, os     from PyQt4.QtGui import QApplication     QApplication(sys.argv)     from qgis.core import *     import qgis.utils          sys.path.append("C:/Users/.../.qgis/python/plugins")          import sextante     sextante.core.Sextante.Sextante.initialize()     input = "D:/Python_Test/a.shp"     output = "D:/Python_Test/b.shp"          sextante.runalg("qgis:convexhull",input,None,None,output)     print "Finished!"      Unfortunately the process fails with the following error:               Traceback (most recent call last):       File "<interactive input>", line 1, in <module>       File "C:/Users/.../.qgis/python/plugins\sextante\core\Sextante.py", line 376, in runalg         alg = Sextante.runAlgorithm(algOrName, None, *args)       File "C:/Users/.../.qgis/python/plugins\sextante\core\Sextante.py", line 276, in runAlgorithm         if not param.setValue(args[i]):       File "C:/Users/.../.qgis/python/plugins\sextante\parameters\ParameterVector.py", line 59, in setValue         layers = QGisLayers.getVectorLayers(self.shapetype)       File "C:/Users/.../.qgis/python/plugins\sextante\core\QGisLayers.py", line 79, in getVectorLayers         layers = QGisLayers.iface.legendInterface().layers()     AttributeError: 'NoneType' object has no attribute 'legendInterface'      What am I doing wrong? Sextante seems to be initialized correctly because               sextante.alghelp("qgis:convexhull")      gives me               ALGORITHM: Convex hull         INPUT <ParameterVector>         FIELD <ParameterTableField from INPUT>         METHOD <ParameterSelection>         OUTPUT <OutputVector>      I´m using Windows 7 Professional 64Bit and the OsGeo4W installation with QGIS 1.8.0 and python 2.7.4.


How to avoid the CRS Dialog when creating a vector layer memory?

I am trying to create a vector memory in a plugin. Basing on some examples,I was able to create a memory vector layer and to add it in QGis with : QgsMapLayer.instance().addMapLayer(myLayer,True) But when the layer is added , QGis always displays the crs dialog to specify the crs of the layers :![crs dialog](http://i.stack.imgur.com/qWy7I.png) I tried to specify myself the crs with : layer.setCrs(QgsCoordinateReference(4326)) but it still shows me the crs dialog. Do you have any ideas to how avoid this dialog ? Best regards


where I can find the shapefiles for the ecoregion map of the US?

I want to find a shapefile for the map at the bottom of http://www.epa.gov/wed/pages/ecoregions/level_iii_iv.htm I would want to create a map showing some data for the 10 forest ecoregions. Thanks.


Duplicate source information from esri.tasks.PrintTask

We are trying out the PrintTask in ArcGIS Server 10.1 pre-release. When an Esri base map is enabled, we are getting source information in two places in the resulting .pdf. I can confirm that this was not the case with 10.1 beta, and we are using the same .mxd's for the service. Is there any way to remove one of these source references, or change the placement? I didn't see a corresponding parameter in the javascript api. ![Screenshot of duplicate source info](http://i.stack.imgur.com/EecR3.jpg)


Representing a TIN in GeoJSON

Regarding the GeoJSON representation of a TIN, are there   1. any specifications or    2. existing implementations I may not be aware of?  You could probably hack one together using the multipolygon feature type, but this seems inefficient. Generally speaking, I don't see any surface types (regular or irregular) represented in the GeoJSON world, is there any existing or upcoming work in this area?


Live KML Feeds - From the Ground Up, How Would You Build It and Why?

I found the "Live KML from PostGIS to Google Earth" post very informative but a bit daunting for someone new to the open source world of products like PostGIS and GeoServer. If I am building the product from scratch, is the PostgreSQL / PostGIS - GeoServer combination the accepted best path to take? What are some alternatives or strengths / weaknesses?


arcpy.ListFeatureClasses(wildcard, type, dataset) - Help with grabbing from 2 different feature datasets?

I'm trying to join multiple line fcs to their matching named point. (point and lines in same gdb, but separate feature datasets. stand alone python script/arcgis 10.1) My initial line of psuedocode code reasoning was:               point_list = arcpy.ListFeatureClasses()     for point_fc in point_list:         get point_fc name         line_list = arcpy.ListFeatureClasses(point_fc name*, "Polyline", line feature dataset)         for line_fc in line_list:             arcpy.JoinField_management(....)         How can I get two different feature class lists? Even if I declare the line feature dataset as separate from the environment workspace/point feature dataset the resulting list is 0 length whereas, if, in a new script, I set the environment workspace to the line feature dataset, then run the ListFeatureClasses, the correct list is returned. Do I reset the workspace before running the 2nd list? A virtual 2ndcupofcoffeethankyou for your help. @artwork21 had a solution to a previous post solution incorporating dataset list that worked here: ArcGIS 10.1 Python script to snap points to lines


How do I copy all feature attributes to a feature in a new layer?

I'm trying to create a new layer in QGIS based on an existing layer. The new features should have all the attributes of the corresponding features in the original layer. I've been trying to use `outFeat.setAttributeMap(inFeat.attributeMap() )` similar to what Victor Olaya is doing at http://qgissextante.blogspot.de/2012/11/spliting-vector- layer.html, but all the attributes come out as NULL. Could you tell me where the error is in my code?               def test(outputFilename):      layer = qgis.utils.iface.activeLayer()      provider = layer.dataProvider()      fields = provider.fields()      writer = QgsVectorFileWriter(outputFilename, "CP1250", fields, provider.geometryType(), provider.crs(), "ESRI Shapefile")      inFeat = QgsFeature()      outFeat = QgsFeature()      inGeom = QgsGeometry()      while provider.nextFeature(inFeat):       point = inFeat.geometry().asPoint()       inGeom = inFeat.geometry()       outFeat.setGeometry(inFeat.geometry() )       outFeat.setAttributeMap(inFeat.attributeMap() )       writer.addFeature( outFeat )      del writer      newlayer = QgsVectorLayer(outputFilename, "Output", "ogr")      QgsMapLayerRegistry.instance().addMapLayer(newlayer)


How to zoom to wms layer after cql filter on Open layers

I created o layer contain points and I added it to a map as WMS layer on OpenLayers I did a CQL filter an it's work but I'd like to zoom to the selected point when I worked with getextent or getdataextent it doesn(t zoom to the right place it missed with a few Km there is my the code of filter function filterNom(){                           var filterType ="cql";                 var filter = document.getElementById('batima').value;                 //var filter = "cité aïn-chock";                 // by default, reset all filters                 var filterParams = {cql_filter: null};                 if (OpenLayers.String.trim(filter) != "") {                     if (filterType == "cql")                          filterParams["cql_filter"] = "batiment='"+filter+"'";                      }                 // merge the new filter definitions                 //mergeNewParams(filterParams);                 wfs_layer.mergeNewParams(filterParams);                 var bounds =  new OpenLayers.Bounds                 bounds = wfs_layer.getExtent()                 //map.zoomToExtent (bounds);                 //var center = new OpenLayers.LonLat( wfs_layer.getExtent() );                 //center.transform(new OpenLayers.Projection("EPSG:4326"), new OpenLayers.Projection("EPSG:900913"));                 //alert(center+','+wfs_layer.getExtent());                 map.setCenter(wfs_layer.getExtent(),18);                 var center = map.getCenter();                 center.transform(new OpenLayers.Projection("EPSG:4326"), new OpenLayers.Projection("EPSG:900913"));                 alert(center.toString());      } and ther is the code to add the WMS layer                 wfs_layer = new OpenLayers.Layer.WMS("Patrimoines",             "http:// localhost/geoserver/cite/wms",              {layers: 'cite:patrimoine' , transparent: true },             {isBaseLayer: false, opacity: 1, singleTile: true, visibility: true,              projection: 'EPSG:900913'});      Thank you for helping me


Esri ArcGIS 10 - using ArcPy to select buffer line by length

I created a buffer, then converted polygon to line (around a city block). I want to select the internal line. I'm thinking since there is no uniform way that ArcGIS labels the interior or exterior buffer line, I can use length to select the line I want by using an if/then statement. I have not been able to find the script command that simply selects a feature/line. Anyone know?


Query SDE to display records created for 'X' date/time

I am looking to query our SDE to display feature classes modified beginning from January 1, 2013 to today. I am using Microsoft SQL Server Management Studio, and as it stands searching through dbo tables but can't find any information that would display feature classes modified. The purpose is to build a business case by identifying a volume of edits/modifications to our data over time. Any help is truly appreciated. ![Below is an image of my work environment](http://i.stack.imgur.com/VUo3k.jpg)


How do I connect to geoserver on my system using java code?

I am writing the following code... But i am getting an error as... java.io.FileNotFoundException: http:// localhost:8080/geoserver/ at sun.net.www.protocol.http.HttpURLConnection.getInputStream(HttpURLConnection.java:1613) at Try.main(Try.java:25) Am i doing it right? I am new at this geoserver as well as client server coding...               import java.io.*;     import java.net.*;          public class Try{              public static void main(String args[]) throws Exception{                  URL hp= new URL ("http://localhost:8080/geoserver");                       try {                      URL myURL = new URL("http://localhost:8080/geoserver/");                      URLConnection myURLConnection = myURL.openConnection();                      myURLConnection.connect();                      BufferedReader in = new BufferedReader(new InputStreamReader (myURLConnection.getInputStream()));                      String inputLine;                      while((inputLine = in.readLine()) != null )                 {                     System.out.println(inputLine);                      }                      in.close();                  }             catch (Exception e){                 System.out.println("Into CATCH");                 e.printStackTrace();             }                    }     }


Why a severe performance difference using deleteSearchedRows

The issue I am having involves two versioned tables. A function deletes 48 rows from Table A and and 48 rows from Table B using deleteSearchedRows. However, the delete from Table A takes roughly 4 seconds while the delete from table B takes roughly 67 seconds. The query filter for both deletes is a single number field that has been indexed in the table, objectID indexed as well. Table A has roughly 45,000 records and Table B 60,000 records. Both tables also have cascade deletes set to No, so that the only data being deleted would be that table's. Table structures are very similar. Table A has OID field + 3 Long integer fields, 2 short integer fields and a date field. Table B has OID, 3 long integer, 1 Text, and a date field. Running the delete query from SQL Plus is very quick, so it's probably not a drive efficiency issue. I broke apart the select and delete sections to see if the search was slow. I used the Search(queryFilter, false) and it returns instantly, but looping and deleting every IRow.Delete takes slightly over a second per row for Table B. I am at a loss on what the difference could be to cause such a severe performance difference in two tables.   Any ideas/debug methods are appreciated. Using ArcGIS 10.0 with an Oracle database


Run function if item in QComboBox is selected

I got another issue writing a plugin to convert three shape files into an .obs input file. I created a UI (Qt) with a QComboBox to select a layer (that works fine). As soon as one of the layers in the QComboBox is selected, a second ComboBox should be filled with the fields of the layer to select from. To do that I'm trying to call a function, that is to be called, as soon, as a layer is selected in the first ComboBox. Unfortuntelly that does not work. Here is what I tried: First I wrote a function to fill the second ComboBox with the field names:                 def build_Layer_Index_Changed():         QMessageBox.information(None, "DEBUG", "Build Layer routine called. :)")         canvas = self.iface.mapCanvas()         dlg2 = ShowDialog2()          build = canvas.layer(int(dlg2.Sel_Build_Layer.currentIndex()))         if int(dlg2.Sel_Build_Layer.currentIndex()) >= 0:           if not build.geomTypeDict == 2:   # Types: 0:"point",1:"line",2:"polygon", Buildings layer has to be of type 2             QMessageBox.critical(self.iface.mainWindow(), ":-/ Selected buillding layer is not of type vector(polygon)! :-/" )            else:             build_fields = self.getFields(build)             dlg2.Sel_buildhight_field.clear()             dlg2.Sel_buildhight_field.addItems(build_fields)      "dlg2.Sel_buildhight_field" is the second QComboBox. I tried to connect it to the first ComboBox's (dlg2.Sel_Build_Layer) Signal by               QtCore.QObject.connect(dlg2.Sel_Build_Layer, QtCore.SIGNAL("currentIndexChanged"), self.build_Layer_Index_Changed())      but I get an error: "TypeError: build_Layer_Index_Changed() takes no arguments (1 given)". If I remove the brackets after self.build_Layer_Index_Changed, I get no eror, but the function is not called (as the MessageBox doesn't appear. I thought I did it exactly as given in the Example at http://techbase.kde.org/Development/Tutorials/Python_introduction_to_signals_and_slots, but obviously I'm doing something wrong…


How to speed up a query for nearest linestring to point data

I have processed a Strahler stream order layer with a lot of river linestrings (type: Multilinestring) and their corresponding strahler number (see information here). Besides I have a point layer. I want to know which linestring is the nearest to the distinct points from the point layer to extract the strahler number for every point. ![enter image description here](http://i.stack.imgur.com/codah.png) My query here is very slow. For around 500 points and 2000 linestrings it takes over 6 minutes:               SELECT DISTINCT ON         (lname)         lnumber,         lname,         strahler,         min(ST_Length(geom)) AS distance     FROM (         SELECT              ST_MakeLine(ST_ClosestPoint(strahler_streams.geom, landslide.geom), landslide.geom) AS geom,             strahler,             lname,             lnumber     FROM landslide, strahler_streams) AS foo         GROUP BY strahler, lname, lnumber, geom          ORDER BY lname, distance;      I know there is the Indexed Nearest Neighbour Search in PostGIS. But when I try this...               SELECT         landslide.lnumber,         strahler_streams.strahler     FROM landslide a, strahler_streams b         WHERE lnumber=2114             ORDER BY b.geom <#> a.geom LIMIT 1      ... it only works when I select one single point (`WHERE lnumber=2114`) How can I speed up the query for the nearest neighbor search (linestring, point)?


From text data create a FeatureClass with table with python

I'm searching how the output of one program (text file with some writings, but also with data as numbers, dates, times, coordinates, all separated with commas) can be transformed to a Feature Dataset. **1) .txt -> .csv -> gp.TabletodBase** and then it's nearly done.. But get problems writing csv with Excel dialect - all row is writen in one cell. So, I create the delimiter ';' to get the data into cell by cell.               f= open(rawtable, 'r')     f2 = open(outtable, 'wb')     csv.register_dialect('dotcoma', delimiter=';')     writer=csv.writer(f2, dialect='dotcoma')          for line in itertools.islice(f, 6, 7): #reading labels that are in 6'th row         labels= line.strip('\n').split(",")         writer.writerow(labels)      But well, then Arcgis does not recognize the table.. I have been trying to cast the data, like integers, floats, strings.. but that did not help.. > **2) .txt -> .shp...** \----> THAT WORKS!!! The shapefile does not have a type of time.. So was skiping it. _EDITED_ However tried the Pyshp library. But after creating fields blindly I record the tuple of all data in a row and I'm getting an error: `IndexError: tuple index out of range`.. I simply write fields from the list of data, depending to which group they go:               for fieldname in labelkeys:         if fieldname in t_int:             w.field(fieldname, 'N', 8, 0) #1probl was here as I made '8' vs. 8 ! stpd         ...         else:             w.field (fieldname,'C', 25)      Another problem with coordinates. I was writing floats... and it did not like it.. So I had to write integers. Plus one cannot pass a tuple to w.record(). Thanks to StackExchange I've learned the benefits of asterix :) and one simply has to write w.record(*mylist). **3) .txt -> feature class...** Here I have problems writing the entire row at once..(like the list of values) into a NewRow. Since field names are made automatically I can not for each make like `row.fieldname=somevariable`. I have tried while looping to SetValues `row.setValue(labelkeys[unit], data[unit])` and then at the end of a list `rows.InsertRow(row)` , but it seemed that it sets just one value. Though I tried by hand and got an error in data casting.. The error comes that I read data like string '34', cast it to int(34) and try to write into field with type 'DOUBLE' (since I grouped some data and wanted to make more generic types).               if labelkeys[unit] in t_int:         if data[unit]:             row.setValue(labelkeys[unit], int(data[unit]))      So for now.. nothing is really working.. Maybe there is some other way I have not though about? If not, could anyone help to get the things work?...


Combine overlapping rasters with priority to values of specific rasters?

I have 5 rasters that cover a large area. I want to combine all 5 into 1 raster. Each raster has data for portions of the total area that none of the other rasters cover. Raster C overlaps A,B and D. While E overlaps D. No other rasters overlap each other. On the portions that are overlapping, only certain Rasters contain the data that is needed to be displayed. The following shows which Raster are needed to be displayed. Raster C > A Raster C > B Raster C > D Raster E > D I've tried `Mosaic to New Raster`, `Mosaic Dataset` and `Mosaic (Data Management)`, all unsuccessful. The result is only the first Raster. The Blank area has a Value of 0, which I believe is overlapping the other rasters value when I combine them. How do I display the Values > 0 when I combine the Rasters giving priority to the values of Raster C when combined with A,B and D and Raster E when combined with D? The edges of Raster D that overlap C and E have missing data as seen below. ![enter image description here](http://i.stack.imgur.com/IlaVG.png) ArcGIS 10.1, advanced arcinfo license. Edit: I've tried all the suggestions below and here are results: ![](http://i.stack.imgur.com/7bYtS.jpg) Above: This is Radouxju's suggestion using raster calculator with IsNULL, the resulting image is Raster C and I tried to overlap only D, D's area is not present. It is an exact image of Raster C. I also tried Raster to New Mosaic with FIRST and LAST, the resulting image was the same as above, but with colour. ![enter image description here](http://i.stack.imgur.com/Ot9GH.jpg) Above: This is Cody's suggestion using Mosaic Dataset and setting priorities, as you can seen the 0 Values overlap the other Rasters depending on their priority. This image is an example/ test to visualize what happens, the priorities here are B > CE > AD.


How to compile the source code of module MOD14 of NASA

I'm working in a project which uses an Active Fire Detection base on MODIS (MOD14) of NASA. The source code and compiled-file can be downloaded at NASA Now I want to compile the source code, although the runable module is already provided. Have anyone ever working on stuff like this?


Inserting strings in language other than English into MariaDB with GeoNetwork?

I am new to **GeoNetwork** and I testing how to insert a text string into **MariaDB** using a language other than English. The database is showing the inserted text string as `??????????`. I have set XML encoding to `UTF-8` AND MariaDB to `utf8_general_ci`. Can you help me to fix this problem? ![enter image description here](http://i.stack.imgur.com/wbDEp.png) ![enter image description here](http://i.stack.imgur.com/JGT9t.png)


how can i convert an excel table to a shapefile with correct scale in arcmap?

i know that there is already one similar question to this one but my problem is different and none of the provided answers did resolve my problem. i've an excel table with some rows and in each row there is the identification and the xy coordinates of an especific point. the xy coordinates were aquired with a gps and are in wgs84. the data layer is in wgs84 too. my problem is that every time and every way i try to import this data to a shapefile, the resulting points are perfect in form but at a ridiculous scale (1:0.89)... exporting the table to dbf and using it gave the same result... how can i make the data import in a way that the scale is correct? thanks in advance =^) \-----EDIT----- if i use the "go to xy" tool to enter the coordinates of two points manually, using the values of x and y in the excel/dbf file, the two points appear in the right places and at the right distance (which implies a correct scale)... \-----EDIT----- excel file: http://dl.dropbox.com/u/3639616/xls%20to%20shp.xls the farthest points are codigo=394 and codigo=297, with about 26.3 km between them.


Any difference between EPSG:3003 and EPSG:102091?

I am working on QGIS 2.2.0. The PRS I am working on is Monte Mario / Italy zone 1 (EPSG:3003), with the following parameters:               PROJCS["Monte Mario / Italy zone 1",     GEOGCS["Monte Mario",         DATUM["Monte_Mario",             SPHEROID["International 1924",6378388,297,                 AUTHORITY["EPSG","7022"]],             AUTHORITY["EPSG","6265"]],         PRIMEM["Greenwich",0,             AUTHORITY["EPSG","8901"]],         UNIT["degree",0.01745329251994328,             AUTHORITY["EPSG","9122"]],         AUTHORITY["EPSG","4265"]],     UNIT["metre",1,         AUTHORITY["EPSG","9001"]],     PROJECTION["Transverse_Mercator"],     PARAMETER["latitude_of_origin",0],     PARAMETER["central_meridian",9],     PARAMETER["scale_factor",0.9996],     PARAMETER["false_easting",1500000],     PARAMETER["false_northing",0],     AUTHORITY["EPSG","3003"],     AXIS["X",EAST],     AXIS["Y",NORTH]]      A dataset I need to use on the above PRS is in Monte_Mario_Italy_1 (EPSG:102091), with the following parameters:               PROJCS["Monte_Mario_Italy_1",     GEOGCS["GCS_Monte_Mario",         DATUM["Monte_Mario",             SPHEROID["International_1924",6378388,297]],         PRIMEM["Greenwich",0],         UNIT["Degree",0.017453292519943295]],     PROJECTION["Transverse_Mercator"],     PARAMETER["False_Easting",1500000],     PARAMETER["False_Northing",0],     PARAMETER["Central_Meridian",9],     PARAMETER["Scale_Factor",0.9996],     PARAMETER["Latitude_Of_Origin",0],     UNIT["Meter",1],     AUTHORITY["EPSG","102091"]]      Is there any difference between the two PRS? (I see a difference in decimals in parameter UNIT["Degree", ...]). In other words, do I need a re-projection or datum transformation to analyse the data together (e.g. buffering, calculating distances etc) ? Thanks in advance. Alessandro


The maximum record length exceeded?

I downloaded some census tract data which are stored in .gdb. I am able to open them in ArcGIS as shpfiles. I selected some tracts and tried to export them into another folder. But I kept getting this error: "Error exporting data. The maximum record length has been exceeded". Is this because the attribute table has too many columns? What is the solution?


Distance between latlongs varies in different projection systems?

I have shapefiles with correct projections. The SRS is EPSG 32735 (WGS 1984 UTM 35S) Now, I have these two coordinates from this projection system: (31.3751,-23.4388) and (31.3779,-23.4355). Obviously the coordinates (format and value) will vary on different projection systems, and there in lies my confusion. I used the following link to calculate the distance between the two points: http://www.movable-type.co.uk/scripts/latlong.html I copy pasted the lat/longs exactly as shown above and the distance it gives me is 0.313 meters. But I'm not sure if this is correct at all, since the calculator didn't ask me for any projection system? So:   1. Does this mean the projection doesn't matter to calculate the distance?   2. Does it differ between different projection systems?   3. Is the tool I'm using incapable? Thanks, Kaushik


How to create generic building footprints based on attribute information?

I'm trying to create a tool that generically creates building footprints on a given parcel according to certain attribute information of the parcel like building density, max. height, utilisation, etc. Is there a way to generically create polygons in model builder with specified dimensions? Hints and suggestions are much appreciated.


Help with Turn-feature-class

I aleady have a street network dataset. I have create a turn-feature-class in a feature dataset. I want to add the turn-feature-class to the network dataset. How do I do that in ArcCatalog? Many Thanks


raster statistics - regressions between raster data

I have the following problem: I have two set of rasters, with data that should correlate...one are measured data, one are modelled data. The correlation is weak, however, when looked at with a normal regression analysis. I suspect that different parts of the base terrain need different modelling input, or that in some parts there is a shift between the correlating cells. So now I try to visualize where these residuals are highest between my two rasters, and where the real outliers are and whether they are in similar terrain. Is there any tool that would allow me to directly compare the rasters in GIS? I run ArcGIS 9.3 with Spatial Analyst and Statistics, and also have QGIS and GRASS on my Linux-based machine. What would be easiest? I am not trying to do a full-on analysis, I am just trying to see where the odd data are and if there is a systematic spacial pattern behind them. Cheers, Peter


Debugging SQL in arcpy.da.SearchCursor RuntimeError: An invalid SQL statement was used?

I am going through a python script in arcpy 10.1 to Select a field value (i.e. NEAR_DIST) based on three condition as                       1) NEAR_DIST<=radius(i.e.500 or anything from user input)              2) NEAR_DIST<>0               3) NEAR_DIST<> -1      I am using 500 which is to be variable asked from user and to be saved in "radius" variable so, i am only using "radius" variable in the SQL statement in following code-                   import arcpy,os,sys         feature =     "C:\Users\Bakul\Desktop\Od\example_shapefile\Polyline_segments_example.shp"         radius = 500         # create a dictionary to hold the list of near ids for each input id         nearest_dict = dict()         arra = []         whereClause = '''"NEAR_DIST" <= '%{x}%' AND "NEAR_DIST" <>'%{y}%' AND "NEAR_DIST"   <>'%{z}%\''''.format(x=radius,y=-1,z=0)         print radius         with arcpy.da.SearchCursor(feature, "NEAR_DIST", whereClause) as rows: ##NEAR_DIST  is a field, I need sort all NEAR_DIST <=500 (i.e radius) and not equal 0 and -1             for row in rows:## debugger does not pass here,                  nearest_id = row[0]                 arra.append(nearest_id)             del rows             del row         print arra      Debugger Error is                   500         Traceback (most recent call last):           File "C:\Users\Bakul\Desktop\Od\Python\scratch.py", line 89, in <module>             for row in rows:         RuntimeError: An invalid SQL statement was used.      Can you tell where is the problem (i think in SQL whereClause)?


Can QGIS open a SDF file created by HEC RAS?

I used the "export GIS data" option in HEC RAS and got a SDF file. Is it possible to open it in QGIS?


Overview effect on postgis polygon table

Is There a Way, with postgis and qgis, to reproduce the overview effect of tiff Raster on a polygon table in order to increase speed? What if we could create simplified versions of this table,gathered under the same view. Until now, i only play with different scale visibility on different simplified version of the table... Best regards


Intersecting polygons in Oracle

Have a polygon in table 2 with its id, no and geometry. In table 1 I have polygons with the same fields. The polygon in table in intersects with a few polygons from table 1.What I am trying to do is clip of the polygons which overlap polygon from table 2 and insert the same fields and geometry in table 3               insert into table 3       select a.store_id,b.store_id,a.store_number,a.client_id,sdo_geom.sdo_intersection(b.geometry,a.geo metry,0.005)       from table_1 a, table_2 b       where b.store_id=34746       and sdo_anyinteract(b.geometry,a.geometry)='True';        When I do run this I dont get any output. Where would I be going wrong


Insert Vertex into Polyline based on Point layer using arcpy

This is essentially the same question asked here that was not answered... Using ArcGIS 10.1 Advanced I have a point layer with points that represent the intersection of a driveway with a street centerline. There is not always a vertex on the street centerline where the driveway intersects, so I am looking for a solution (preferably using arcpy) that loops through each point and inserts a vertex on the line that intersects with that point. I do not want to use Split Line at Point method because I don't want to split the line. Looking for any ideas on how to accomplish this...


OpenLayers: draw path with only two points

I want to draw a path with no more then 2 points. On the second point I want the path to be finished. I'm using the OpenLayers DrawFeature, I set a callbacks for "point" where I count the created points and if the points are two I'm finishing the sketch with `finishSketch()` method. The problem is that when finishSketch() method calls the `finishGeometry()` method which removes the last added point from the path which cause the following error to show `"Uncaught TypeError: Cannot read property 'geometry' of null "` My question is the any legit way for finishing a path programmaticaly or this is a bug in the DrawFeature control.


Partially fill a sink - ArcGIS

I am trying to fill sinks but only to a certain level. For example - if I want to fill them by 5 meters - this means that a sink that is 15 meters deep - will fill only to 10 meters and a sink that is 3 meters deep would fill completely. Using ArcGIS 10.0 - service pack 4. Any help of this would be appreciated. I thought I created a script to do what I wanted, but my client clarified something different than what I though. Attached is an image of what I have and what I want for a custom fill tool.![enter image description here](http://i.stack.imgur.com/MjEVl.png)


Split polylines at polygon boundaries in QGIS

I have a set of contour lines and a polygon set containing glaciers. Because I want the glaciers to have blue contours, I need to split the contour lines at the glacier boundaries. I can achieve this using the following method: 1) Intersect the contours with the glacier polygon to create the glacier contours. 2) Create a spatial difference between the glacier polygon and an extent polygon. 3) Intersect the contours with the spatial difference from (2) to create the blue conoturs. 4) Merge the two results together if I want it in a single data set. But the following method would be more convenient, and usable for other types of analyses as well. 1) Combine a polyline dataset (or point dataset for that matter) with a polygon dataset. The lines with be split at the polygon boundaries and each output line (or point) will take the attributes from the polygon it is inside. Such a function would make my problem so much simpler to solve, as I can do it in one single operation. Is there a function for this in QGIS?


How to add feature class to MXD with ArcPy (Python)?

I am using ArcGIS 10.1 and I am trying to add a feature class that I create to an MXD, so that I can join to it. This all needs to be done using Python. Is this possible?


Japan Data: Postcode/Census

Hi I'm primarily looking for a GIS data source for postal codes in Japan as well as census data for different scales of geography. Has anyone found a good source of these data sets that is in English? Thank you for your help.


Retrieving the proportion of specific pixel class from CHM rasters

CHM = canopy height model I am very new to R (specially for GIS purposes) and I am struggling a bit in getting the following pixel statistics:   1. Proportion of zero value pixels = [nº of zero value pixels/(total nº of pixels - nº of NA pixels)]   2. Proportion of pixels smaller than 3rd quantile = [nº of <3rd quantile pixels /(total nº of pixels - nº of NA pixels)] Here is a reproducible example:               library(raster)       f <- system.file("external/test.grd", package="raster")       r <- raster(f) #r is the object       plot(r)       summary(r)


How to loop through a list of shapefiles and rename them?

I am trying to upgrade a simple script to list all shapefiles in a folder and if necessary in any folders within and then loop through the list of shapefiles and rename them.               import arcpy          # Set work environment     arcpy.env.workspace = "C:\\Users\\Rachael\\Documents\\Export"          # List shapefiles in folders     shp = arcpy.ListDatasets()          try:       for shp in arcpy.ListDatasets():         if arcpy.Exists("Histori2.shp"):           arcpy.Rename_management("Histori2.shp", "H_BEV_CANS.shp")