How can I calculate different polygons by different attributes?

I'm doing my master thesis and want to calculate different polygons (in QGIS 1.8 for Mac) in four different columns which I defined by different attributes. i.e. some of the attributes in the first column are named "1112" and some "1113",.... So how can I calculate the area of the polygons in column "1" with the attributes "1112"? The numbers here are only examples, my thesis is about greenspaces, so it would be to complicate to explain here...all of the numbers have a different meaning and have to be joined. I need to sum these attributes to make a statement for the thesis. Can anyone please help me? Thanks a lot   Markus @M_Gebhart


ArcScene TIN smooth

I created a TIN with ArcScene out of laserscan point data. I want to smooth the spikes (image below), to get a nice crossover on the elevations of the surface. Do you have any advice how to do that? ALready tried to convert it into a raster, but i didn't know which options would be the best. ![TIN](http://i.stack.imgur.com/p89QB.png)


Compile Qgis source : Links errors (Windows)

I am actually trying to compile the source of the qgis git : https://github.com/qgis/Quantum-GIS In order to compile it, I follow the steps of the tutorial here : http://www.qgis.org/api/INSTALL.html (starting on the step 4 "Buildings on Window") I manage to configure the CMake with all the dependencies and to generate the project. To generate the project I use MSCV10 but I also tried with MSCV express 2008 like in the tutorial. In the two cases, when I generate the project, I always have the following errors : > LINK : fatal error LNK1181: impossible to open the file > '....\core\Release\qgis_core.lib' > > LINK : fatal error LNK1181: impossible to open the file > '....\analysis\Release\qgis_analysis.lib' I don't understand why the project tries to open a librairie which is not created yet to compile its code. Do you have an idea how to solve this problem ? EDIT : Apparently, the qgis_core.lib and others must be created during the build (this is why there are theses errors). But the lib is not created because of the errors : > .\flex_qgsrastercalclexer.cpp(43) : fatal error C1083: imposible to open the > hearder file: 'unistd.h' : No such file or directory Do you know how to add this header to the project ? I am not very used to MSCV Thank you very much Regards


Topology - To be, or not to be, that is the question

Topology, I mean "real" topology, is not a new concept and it's available since a long time in GRASS gis. I won't speculate the reasons myself but, after a long period of spaghetti model dissemination, recently some attention is leaning towards Topology implementation. One tool of major importance is Postgis2.0 and one great example of a kind of topology in action is openstreetmap. Qgis developers is also generating some discussion about the subject. My experience with Grass was great, there is nothing more pleasant than two areas sharing the same boundary and "edit time" topology check. This presentation also describes great advantages of using real topology. Now the question:   Considering web mapping and desktop gis and with a nice argumentation based on personal experience, industry requirements, tools available, etc. Do you think that real topology is to replace spaghetti model and be widely adopted in a near future?


Joining attributes of polygons and points by location in QGIS

I have two shapefiles....   * points containing a number field 'response_time'   * a polygon grid of square polygons containing fields including 'PNTCNT' with the number of points per polygon, calculated already using the tool 'points in polygon' There are between 0 and 600 points within the bounds of each polygon in the grid. I want to add fields to the polygon grid for the mean, min and max response times within each polygon, in addition to the PNTCNT field already there. I have tried the 'Join attributes by location' tool, making the polygon layer the target and the points the join layer, and vice versa. I've selected 'take summary of interesting features'. But I can't get the result I want. All I end up with is the polygon grid again but with an extra field 'COUNT' which replicates (or duplicates) the 'PNTCNT' field. No min, max and mean of values from the points layer. What am I doing wrong? Is there a better way to do this? I have tried this using QGIS 1.7 and 1.8.


QGIS Python version

I have brewed python 2.7.3 onto my Mac (Lion), and installed shapely via `pip install shapely`. The Polygonizer plugin, however, does not find the shapely library. I see when QGIS loads a Python console, it loads 2.7.1. Is there a way to change this? Do I need to brew QGIS to get it to play nicely with my other python libraries?


Setting up TileCache to make use of Mapnik and OpenStreetMaps

First of all, I'd say that I'm completely newb at maps related stuff... I am trying to install TileCache and configure it to use Mapnik to render tiles from OpenStreetMaps. I managed to setup the Tilecache basic CGI conf, but now I do not know how can I continue with setting up Mapnik and my own tiles because I am not able to find any tutorial on the net explaining the steps for a complete TileCache setup for OSM and Mapnik (neither in TileCache's web it is explained), and, since it is the first time I deal with WMS and maps rendering and so on, I find difficult to understand the proccess... Could somebody help me, please? BTW, I need this to run on Windows Vista/7 Thanks in advance!


Does Polygonize in QGIS 1.8.0 work correctly?

Does anyone else have issues with the Polygonize function in QGIS 1.8.0. When i use it is to extract details from larger raster images, through raster calculator, then polygonizing the resultant image. It works fine in 1.7.4 but never in 1.8.0. never asked before because i just made do with 1.7.4 Cheers. Maybe a bug issue or something but do not know where to suggest it.


Status of the TMS specifications

I am trying to implement some TMS specs as per http://wiki.osgeo.org/wiki/Tile_Map_Service_Specification. I am particularly interested in the service discovery part ( and ). The page is very clear wrt status of the specs, ie experimental. I am wondering if it is worth implementing it at all in fact. Bottom line is that I am struggling to find a server that would actually comply to the specs above mentioned. Does anyone know of such a server? + state of the specs? should I just stick to the good old url pattern and drop all the server service discovery stuff? Thanks for any pointers, Thomas


How can I export an image from an ArcGIS map service where the bounding box wraps?

I'm trying to export maps from an ArcGIS REST Service using WGS84 and have run into difficulty when the bounding box ends up exceeding latitude 180. I was hoping that the map would wrap, but it doesn't. In some cases the resulting image ends up zoomed very far out and not anywhere I intended it to be. This is an example of where the map gets clipped where I'd prefer it wrapped. This is an example of where I've exceeded 180 and then the image is zoomed out and not on the location I intended. Using the coordinates of a location I buffer it with a rectangle for the bounding box, which means that in some cases the bounding box will overlap -180 or +180. The resulting bounding box then doesn't get my desired result. What is the correct way to get this to work?


Remove layers references from geoprocessing tool

When I run my tool it is calling some geoprocessing tools and performing some operations on some layers but when i delete those layers from arcmap then those layers are removed from arcmap but not from Geoprocessing tool. Please tell me how to remove layers references from geoprocessing tools programmatically.


Zip code highlighted in google map

When we search for a zip code in Google Maps we get the area highlighted in the map with a dotted border. How is this calculated? How do they know how to draw the polygon? I'm looking at implementing something equivalent, an overlay on each zip code of the visible area of the map. Any ideas are highly appreciated ;)


Find Gap between two Continuous One Way street in Road Feature Class

I wish to find gap in the continuity of One Way street.I have One Road(Polyline) Feature Class and one way road are represented by one column called ONE Way.If ONE WAY="Y" that means road is one way so i wish to fetch out those polylines that fall between two continuous oneway road and its ONE WAY cloumn has no value. Let me tell you , what i have done so far.Open the feature class as pfeatureclass and done query on it.Selected those features which are oneway.               prdQf=New Queryfilter          prdQf.whereclause="one way="& "YES"     prfcur=pfeatureclas          s.search(prdQf,false)          proadfeature=prfcur.Nextfeature     do until  proadfeature is nothing     pgeom=proadfeature.shape     pPolyline=pgeom     pPolyline.QueryToPoint(topt)     penvelope=topt.Envelope     penvelope.Expand(0,0,False)     Dim sPf as Ispatialfilter=New spatialfilter     sPf.Geometry=pEnvelope          sPf.GeometryField=pFeatureClass.Shapefieldname     sPf.Spatialrel=esrispatialRelEnum.esrispatialRelEnvelopeIntersects     pfc3=pfeatureclass.search(sPf,false)      Find all the edges that are intersecting it and one way is blank and then again find all intersected edges and their one way attribute. if one way value exist, that means there is gap. But its not providing the correct error.???


Serialize python list containing geometry

In simple case if list doesn't contains geometry all works fine:               result = [dict(f.fields, id=f.id, label=f.label) for f in features]     json.dumps(result)      But if `result` list contains geometry (Shapely) I get error `object at 0xb4a87bcc> is not JSON serializable`:               import geojson     result = [dict(f.fields, id=f.id, label=f.label, geom=f.geom) for f in features]     geojson.dumps(result)      If I try to serialize not list but one dict - all works fine:               geojson.dumps(result[0])      Is it possible to serialize list of dicts?


How to install Postgresql and postgis on Ubuntu 13.10?

Dear all open GIS geeks, I have upgraded my OS from Ubuntu 12.10 to 13.10. I had a clean installation of 13.10 version. Can someone please suggest which version of Postgresql and PostGIS can be installed on Ubuntu 13.10 version and how to instal it? I have looked at the steps given in How to install PostgreSQL 9.2 with PostGIS 2.0 on Ubuntu 11.10 (or higher)?? But the three repos are not adding after update. Hence that answer is not helping me.


Multiple Cursors

Below is the sql statement which I have explained below.                Declare         Cursor c1 is          select store_number ,0 as total_area,0 as BG_ID ,0 as store_geom         from table1  where CLIENT_ID= 1 and ORG_ID = 1 ;              Type C1_TAB_TYPE is table of c1%ROWTYPE;               c1_list c1_TAB_TYPE;                Begin                    For r1 in c1                     Loop                select  BG_ID,store_number,total_area                   BULK COLLECT INTO c1_list from                         (                     SELECT BG_ID,store_number,b.geometry store_geom,                    ( (SDO_GEOM.SDO_AREA(SDO_GEOM.SDO_INTERSECTION(a.geometry,b.geometry, 0.005), 0.005, 'unit=sq_mile') /SDO_GEOM.SDO_AREA((a.geometry), 0.005, 'unit=sq_mile'))*100) Total_area                     FROM table2 a, table1 b                     Where store_number != r1.store_number                     AND sdo_relate(a.geometry, b.geometry, 'mask=anyinteract') = 'TRUE');                                if  total_area=100                         then                                              For i in 1..c1_list.count loop                                     insert into temp_Prop_area_100                                     values c1_list(i);                                End Loop;                           else if                            For i in 1..c1_list.count loop                                     insert into temp_Prop_area_block                                     values c1_list(i);                                End Loop;                        End if                      End Loop;              End;      So I select a record(the biggest geometry of all three tables) from table 1 and overlay over smaller polygons from table 2. I get the list of the overlapping geometries from table 2. This may have a few of them intersecting and a few of them completely under the polygon from table 1. Then I want to put them in to two different tables . The ones that are completely under in temp_prop_area and the ones which intersect in temp_prop_area_block. Now the issue I am having is, when I get the id's for all intersecting polygons(temp_prop_area_block) I want to overlay each of this polygon over the polygons from table 3 which area smaller polygons.And similarly find out which polygons do they intesect and what is the area of intersection. The polygon size are in this order polygons from table 1 as the biggest, the from table 2 or temp_prop_area_block and then the polygons from table 3.


change font size of elements in dropdown menu store

I am creating of list of building names from a query, which then populates a dojo.ItemFIleReadStore. This store then populates a dropdown menu (using dijit.Form.ComboBox), where the user makes a selection. My problem is that I cannot change the font size of the dropdown list. I can change the font size of the placeholder in the ComboBox, and also of the final selected result, but not the font size in the dropdown list itself. Here is the HTML for the combo box:               <div id="list">       <select id="menu" dojotype="dijit.form.ComboBox"       value="autoComplete="true" placeholder="Options:" class="comboBoxClass"        required="false" onchange="executeSelectedQuery(this)">       </select>     </div>      and here is the basics of the JS for the store item:                var dataItems = {         identifier : 'name',         label : 'name',         items : values     };          var menuStore = new dojo.data.ItemFileReadStore({         data : dataItems,         clearOnClose : true     });      I can change the font-size of the placeholder with .dijitInputField and change the font-size of the result string with : input[type="text"] in css, but neither of these will change the font size of the stored drop down list of items. Ideas? (I am using ArcMap 10.0 and ArcGIS Javascript API 3.4)


Is it possible to automate the detection of clearcut areas in optical satellite imagery?

Given optical satellite or aerial imagery of a resolution similar to the image below, is it possible to accurately detect clearcut areas? Is their research or documentation that supports a specific methodology (e.g. object-oriented classification which considers contrast and geometry)? How can this process be automated over areas ~1000km^2? ![Clearcuts](http://i.stack.imgur.com/oI7n1.jpg)


ESRI JSON-thingie from ArcGIS online

I try to download a web map from arcgis online into arcgis desktop. From the site I get a file named item.pkginfo and arcgis desktop opens it. Unfortunately it only gets me the basemap image and not the features. According to the knowledgebase article 3824 from Esri, one can reach the data via the url after the tag but what i get is something looking like a long string of json-like data. (The esri-page doesnt say anything about the format one get.) Haw anyone tried to extract data from arcgis online the same way and can tell the data-format of that string? The datastring is this: http://www.arcgisonline.com/sharing/content/items/4dde65a4ffd846ff862ed6c508853e37/data?f=pjson for reference. I admit its a rather vague question but I can't really formulate it better.


Create perimeter from given area

I need to develop a web-based application to manage areas inside a bigger area. I already have the initial perimeter, and it has irregular borders like that: ![Map1](http://f.cl.ly/items/1f140A040d1K0s0j2V1d/map1.png). I will need to tell the app to create a perimeter based only in a given area.               createPerimeter(1000); // create perimeter with 1000m2      That should create a perimeter with 1000m2 adjacent to the border and to any other existing area. ![Map2](http://f.cl.ly/items/2U2b2R2S3D122v0M1x3J/map2.png). And so on. The above method should create adjacent areas until the end of the container area ![Map3](http://f.cl.ly/items/1p3z2I2T3Z2s1c3V0Z2q/map3.png). I'm looking for a solution regardless the language or platform, directions to how to start it. I believe POSTGIS has many functions that could help me, but I have not found a solution so far. Regards.


ArcGIS - How to programmatically create new dynamic layer for service?

I'm trying to create a number of widgets that will visualise different data by geographic area, e.g., the _Potential Losses for Instrument X in Year Y for the Counties of State Z._ Ideally this would be a map of the state (tick can do that) with the showing the county boarders (use javascript api to display basemap layer ) and then colour each county according to loss. :( Going down the ArcGIS route I am assuming I need to create a dynamic layer here with all of the polygons for each county of each state defined and then somehow colour them according to data pulled out of a DB showing loss for that instrument, for that county for that year. Could someone please   * A) validate that that is the correct way to go and   * B) point me in the right direction to some resources to help me learn how to achieve this.


How to gradient buffer a background layer

I have a map that was created by a contractor working for my employer before I began. I am not sure what GIS software the map was created in but it has an effect I would like to be able to reproduce in arcGIS 10.1. The map seams to have a gradient, buffered layer on the outside of the Australia coastline layer. This layer grades out from dark grey to white. Sitting off the coast in my layout view I have the Australia with the extent area showing. What I cant get is the gradient layer to gradient out behind the Aus extent part. In the maps I make the Aus is in a separate data frame and I cant make the background of that data frame transparent. I can't post the maps as I'm a new user but if someone can assist me I can send you directly the maps to explain better what I mean. Any assistance will be greatly appreciated. Cheers


OpenLayers Rule is not working as expected

I'm trying to add a rule to paint some roads of my GeoJSON file with a different style using OpenLayers. The problem is that even after adding a property "featureType" to the feature in response to the event "featureadded", the rule doesn't seem to be working. Any ideas why this happens? This if a feature example in GeoJSON format after changing the "featureType":               {         "type": "Feature",         "properties": {             "id": "42",             "state": "texto",             "featureType": "segmento_logradouro"         },         "geometry": {             "type": "LineString",             "coordinates": [                 [                     574951.47711036,                     9673700.3858644                 ],                 [                     574954.47530982,                     9673558.2712099                 ],                 [                     574814.15957496,                     9673551.0755312                 ]             ]         },         "crs": {             "type": "name",             "properties": {                 "name": "EPSG:29193"             }         }     }       And this is the code I used to add the rules:               var ruleSegmentoLogradouroStyle = OpenLayers.Util.applyDefaults({         fillColor : "#0000FF",         fillOpacity : 0.9,         strokeColor : "#FF0000",         strokeOpacity : 0.9,         strokeWidth : 5     }, OpenLayers.Feature.Vector.style["default"]);          var ruleSegmentoLogradouro = new OpenLayers.Rule({         filter : new OpenLayers.Filter.Comparison({             type : OpenLayers.Filter.Comparison.EQUAL_TO,             property : "featureType",             value : "segmento_logradouro",         }),         symbolizer : ruleSegmentoLogradouroStyle     });          var ruleDefault = new OpenLayers.Rule({         symbolizer : OpenLayers.Util.applyDefaults({},                 OpenLayers.Feature.Vector.style["default"]),         elseFilter : true     });          defaultStyle = OpenLayers.Util.applyDefaults({         fillColor : "green",         fillOpacity : 0.4     }, OpenLayers.Feature.Vector.style["default"]);          MapControl.editLayer.styleMap = new OpenLayers.StyleMap(defaultStyle);     MapControl.editLayer.styleMap.styles['default'].addRules([ ruleSegmentoLogradouro,     ruleDefault ]);


Plot points map of population

I want to plot kind of isarhytmoc points map like in this tutorial, but with population data. However, unfortunately, I have no enough skill and only came with such code:               require(sp)     require(rgdal)     require(RColorBrewer)     require(ggplot2)     library(plyr)     library(maptools)     library(rgeos)     require(gpclib)     gpclibPermitStatus()      gpclibPermit()               rus<-url("http://www.filefactory.com/file/4h1hb5c1cw7r/n/RUS_adm1_RData")     print(load(rus))     proj4.str <- CRS("+init=epsg:3413 +lon_0=105")     gadm.prj <- spTransform(gadm, proj4.str)               popul <- read.csv2(file="C:\\unempl11.txt", header = TRUE,              sep = ";",quote = "", dec=",", stringsAsFactors=F)          df <- fortify(gadm.prj, region = "ID_1")     df2 <- merge(df, gadm.prj, by.x="id", by.y="ID_1")                                   p <- ggplot(df2, aes(x = long, y = lat, group=group)) +            geom_point(data    = popul,                         mapping = aes(x=lon, y=lat, colour=abs),                       size    = 3,                       alpha   = 0.8) +            scale_colour_gradient2(name  = "Population",                                   low   = "darkred",                                   mid   = "white",                                   high  = "blue",                                   guide = "colorbar") +            ggtitle("Population in Russia)")     p + geom_path(data    = df2,                    mapping = aes(x=long, y=lat, group=group),                   size    = 0.125)      Here is the data I use, if anyone can help me with some ideas i will be grateful.


Multiple text formats in a qgis label

Is there any way to have multiple text formats in a single qgis label? I am making a map that shows mountain summits. Many of these summits have more than one name and I would like to show each name in a label and have each name rendered with different formatting to show what type the name is. For example, summits, in the US, often have an "official" name from the GNIS database, a name as shown on a 1:24000 USGS topo map (often just the elevation of the summit), and a local nickname. I would like to be able to have a label that renders each name in a different type font and/or font size.


How to add a legend for mapfile (Raster Data)

What I would like to do is shown in the figure below. How can I add this color ramp legend in a mapfile which will also be linked to a raster. I don't know how to get it working while drawing a figure based on that mapfile. Thanks ! ![ArcGIS](http://i.stack.imgur.com/8pDXY.png)


qgis field calculator - changing datatype

![enter image description here](http://i.stack.imgur.com/oevrn.jpg)I have seen a solution that I don't get to work - I need to change data type from string to integer and I did what is shown i the image below Can someone point out what I do wrong? thx Jan


What percentage of ArcGIS users are using version 10.1 already?

I'm planning on building an add-in using ArcPy, but some of the functionality I need is only available in 10.1. I want to distribute this add-in and was wondering if making it 10.1 only would limit adoption. So would like to know what percentage of ArcGIS users are using version 10.1 already?


Reduce number of vector field arrows in qgis

OK, now that I have figured out how to get my vector fields in qgis, I can't figure out how to reduce the number of vector arrows on the map. Is there a setting for how many vectors are sampled? Thanks.


How can I view the Certificate of a code-signed ESRI AddIn (*.esriAddIn) without ESRIRegAddIn.exe?

For users with ArcGIS 10 installed, it's easy to view the certificate for a signed AddIn since it is presented to you in the Confirm Add-In file installation window. Is there a way to view this certificate information if you don't have ArcGIS 10 installed? As an alternative, can the ESRIRegAddIn.exe be copied to a computer without ArcGIS 10 and successfully run in order to view the certificate information?


Specifying ArcGIS Server WMS GetFeatureInfo tolerance?

is there a way to specify pixel tolerance when requesting feature info via GetFeatureInfo WMS request for an ArcGIS Server WMS-published map service? I have not been able to find any such information in the documentation. The OGC WMS standard does not specify any request parameters for coordinate tolerance, I just wonder if there's similar setting somewhere which could make retrieving feature info a bit more convenient for the user, especially for point features.


Load DTED data files into ArcGlobe

When I add DTED data files (*.dt1 and *.dt2) to my ArcGlobe application, they are displayed as solid gray boxes. Do I need to do something special to the files and/or the ArcGlobe control to see these properly?


Is it possible to query for values in fields with the column names partially matching another field's values in psql?

I am trying to query the values of a series of columns NAME_# with the # matching the values of another column LEVEL and then populating another column with what I retrieve


Tunnel distance approximation

Having two points (longitude, latitude) on earth, if I am sure their distance is under 100 km; how much error will I have using tunnel distance approximation? (Defined at here) Note: Currently I am using haversine formula, but since my points are normally in a 2 km range (at 99% of time and for other 1% it goes from 2 to 100 km) I think I do not need that complex calculation.


How to create a "Datum-less" projection in QGIS?

I want to use QGIS to do some georeferencing and image registration. On ArcMap10 I can define an "Unknown Projection", basically a default metric grid, and use it as a grid for georeferncing. In QGIS however, when I tried to georeference an image with GDAL, it gave me the following error: > Failed to compute GCP transform: Transform is not solvable I've looked around and figured that the problem lies with the projections. I've imported a 4 line CSV with 4 points demarcating the image extent in meters (BL 0,0, UR 450,450), but the layer properties tell me it's in WGS 84. because I do not know the actual coordinates of the image, I have to use just a basic internally coherent 5*5 metric grid, that I can georeference with corresponding marks in the image. This can only work with a datumless metric projection, which I do not know how to define (and save it as a custom defined projection) I'm new to QGIS, so the problem might be remarkably easy to solve, but after a few hours of searching, I still haven't found it. appreciate any help, H.


Polyline length calculation in ArcGIS 10/VB script

I used to run this length calc on my data in ArcGIS 8 and 9, and now it don't work in ArcGIS10. Can you tell me what I have done wrong?               Dim dbllength as double     Dim pcurve as icurve     Set Pcurve = [shape]     Dbllength = Pcurve.Length     Length = dbllength


How to split a polygon into quadrants?

I have a Alberta province (Canada) section polygon shapefile in GCS_NAD83 with 255,000 polygons. I like to split each section into 4 equal quadrants (NW, NE, SW, SE) ![enter image description here](http://i.stack.imgur.com/cC9KR.jpg) The Alberta province has several UTM zones, so I guess I can convert it to Lambert Conformal Conic or Albers Equal Area CS if the section scan be split only in projected units. There are people selling this information, but I was wondering if this is something I can create using the data I have? Please let me know. Thank you.


Creating Longitude and Latitude Lines in Google Maps

I'm working on an independent high school project where i'm trying to build a simple differential equation model for diffusion. I'm trying to overlay the model on top of a map to visualize spread, so for example, oil in water. Right now I have a rectangular box that have the coordinates of ( 31.0195,-102.8072) NW, (19.2019, -102.8072) SW, (31.0195, -75.3963) NE, and (19.2019, -75.3963) SE. I am trying to build a 2D array of lat/long coordinates spaced evenly. For example, an array of lat/long coordinates that are 500m apart. Is there a formula/algorithm where I can convert these coordinates to a distance, break the distance up into 500m chunks, and then get lat/long coordinates along the way (i.e. so I have lat/long coordinates every 500m)? Thanks


Using Fiona to write a new shapefile from scratch

I am new to Fiona, and trying to wrap my head around reading a csv file and writing out an ESRI shapefile. I have most of it down, but am missing two things; (1) defining the projection/spatial reference and (2) the attributes actually being populated with the values from the csv file. So my questions are: 1. How do I define the spatial reference, either by defining an EPSG or a custom projection 2. Why is my data not writing as attributes? Thanks in advance for any guidance.               import csv     from shapely.geometry import Point, mapping     from fiona import collection     from fiona.crs import from_epsg               schema = { 'geometry': 'Point',       'properties': {'dx': 'float:13.3',                 'dy': 'float:13.3',                 'dline': 'str',                 'dtrace': 'int',                 'ddepth': 'float:9.4',                 'dtrash': 'float:9.4',                                'sx': 'float:13.3',                 'sy': 'float:13.3',                 'sline': 'str',                 'strace': 'int',                 'sdepth': 'float:9.4',                 'strash': 'float:9.4',                                  'bx': 'float:13.3',                 'by': 'float:13.3',                 'bline': 'str',                 'btrace': 'int',                 'bdepth': 'float:9.4',                 'btrash': 'float:9.4',}}          with collection(     shapeout, "w", "ESRI Shapefile", schema) as output:     with open(finalout, 'rb') as f:         reader = csv.DictReader(f)         for row in reader:             point = Point(float(row['dx']), float(row['dy']))             output.write({                 'properties': {'dx': 'float:13.3',                 'dy': 'float:13.3',                 'dline': 'str',                 'dtrace': 'int',                 'ddepth': 'float:9.4',                 'dtrash': 'float:9.4',                                'sx': 'float:13.3',                 'sy': 'float:13.3',                 'sline': 'str',                 'strace': 'int',                 'sdepth': 'float:9.4',                 'strash': 'float:9.4',                                  'bx': 'float:13.3',                 'by': 'float:13.3',                 'bline': 'str',                 'btrace': 'int',                 'bdepth': 'float:9.4',                 'btrash': 'float:9.4',},                 'geometry': mapping(point)             })      With the code updated to the "newer" (aka, not outdated) process:               with fiona.open(shapeout, 'w', crs=from_epsg(3996), driver='ESRI Shapefile',                  schema=myschema) as output:     with open(finalout, 'rb') as f:         reader = csv.DictReader(f)         for row in reader:             #geometry             point = Point(float(row['dx']), float(row['dy']))             #attributes             prop = {'dx': float(row['dx']),                     'dy': float(row['dy']),                     'dline': str(row['dline']),                     'dtrace': float(row['dtrace']),                                         'ddepth': float(row['ddepth']),                                         'dtrash': float(row['dtrash'])},              print{'geometry': mapping(point), 'properties':prop}      Gives me this:               {'geometry': {'type': 'Point', 'coordinates': (-245215.3333, -975226.3333)}, 'properties': ({'dtrace': 24387.0, 'ddepth': 1433.2458, 'dtrash': 1433.25, 'dx': -245215.3333, 'dy': -975226.3333, 'dline': '      lsl0904.201302'},)}     {'geometry': {'type': 'Point', 'coordinates': (-245215.0, -975251.0)}, 'properties': ({'dtrace': 24386.0, 'ddepth': 1432.7268, 'dtrash': 1432.73, 'dx': -245215.0, 'dy': -975251.0, 'dline': '      lsl0904.201302'},)}      But as soon as I change to:               output.write({'geometry': mapping(point), 'properties':prop})      I get the error message: AttributeError: 'tuple' object has no attribute 'keys'


Running algorithm from Processing takes much much longer time than doing it from Qgis GUI

I have a road layer with 32000 features, and I want to find all the intersections. when I use processing->SAGA -> convert polygon/line vertices to points, it is all done in less than 5 seconds. However, when I try to do the same thing in PyQgis:               vlayer = QgsVectorLayer("Point", "temp_points", "memory")     c=processing.runalg("saga:convertpolygonlineverticestopoints",layer,vlayer) # get intersection points      it takes about 15 minutes. I would appreciate any help on this.


Fastest way to merge a feature

I am writing a arc engine application that can take a shapefile, and merge a set of polygons together. Here is the code that I am using to merge the polygons:               private static void mergeFeatures(IFeature merger, IFeature mergee)     {         ITopologicalOperator union = (ITopologicalOperator)mergee.Shape;         mergee.Shape = union.Union(merger.Shape);         merger.Delete();         mergee.Store();                 }      Is this the best solution to merge polygons? In my application this part of my code takes 60% of the execution time. Is it slow because I am using a shapefile? Does a file geodatabase help with this? Or is this even just the totally wrong approach for merging polygons. Just to clarify, I am taking a large set of polygons and am trying to merge ones under a certain number of acreage. From these polygons that are small, I need to merge polygons with similar attributes (for example on a vegetation layer, ones that are similar vegetation.) So this can incur a large number of small polygons getting merged together, and eliminating some of the merges can be difficult. If there is an easier way to bulk merge, then that would be best, but thus far, I have not found one.


Does IDistanceOP.EucDistanceFull require a projected coordinate system?

I have designed an implementation of ArcObjects' IDistanceOP.EucDistanceFull. It takes its input as a feature class and converts it to a GeoDataset with which EucDistanceFull can work with. It works fine with feature classes in projected coordinate systems, but doesn't execute property on feature classes with unknown or geographic coordinate systems (It will output a blank or black raster). What could be the issue here? I don't see anywhere in the library documentation that states that EucDistanceFull requires projected coordinate systems.


Change default leaflet marker color

I'm using esri-leaflet, and leaflet inherently to add a point feature class (geojson) to a map. What is the best way to add and style the layer? I've see pointToLayer on the leafletjs.com, but don't want to use the circleMarker. Below is a code snippet of what I am working with at the moment:               new L.esri.featureLayer(schools, {         style: function (feature) {             return{                 weight: 2,                 opacity: 1,                 'marker-color': '#FFF'             }         },         onEachFeature: function (feature, layer) {             layer.bindPopup(L.Util.template(m_template, feature.properties));         }     }).addTo(map);      I would like to use AwesomeMarkers, but I don't know how to implement that with a featureLayer. From my understanding featureLayer extends L.GeoJson. The 'marker-color' has not worked. A live preview of this is available HERE


How to capture accurate GPS points using mobile devices?

I'm currently working in Vietnam as a climbing guide, but am putting my GIS skills to use by doing a side project for the company involving the mapping of climbing routes. Along with another employee who is an app developer, our goal is produce a digital guidebook for the area. For this project to work for me, I need to be able to take fairly accurate GPS points, as there are features just a few meters ap[art from each other that need to be distinguishable, but obviously survey grade equipment is not and option for me here. I've been using the app GPS Averaging to take several hundred points for each location. Despite taking these averages, I'm having very little luck achieving accuracy in my output points, and I suspect my precision is low as well because the in app reported error gets very low (0.8m) after about 20 collected points, then goes up to about 3m after 200-300 collected points. However when I import these points, the actual errors appear to be well beyond 3 meters. What can I do to improve GPS point accuracy from a mobile device? Any suggestions would be hugely appreciated. I've currently only worked on this with android devices, but have access to ios devices too if they provide superior GPS chips.


Difference between 10.1 and 10.2 Map/Feature Service layer ID's?

I noticed that when I set up a 10.2 ArcGIS server and published a map service, the layer ID for the dynamic map service was zero, but each of the FeatureLayers underneath started with 1 and continued on. I wrote some code which made the assumption that the first feature layer started with a 1, but when I tested my code with a 10.1 server, it didn't get the first feature layer since it started with 0. Is this a dependency that changed from 10.1 to 10.2? Thanks for any help :) For now, I will (and should have before) try to handle both cases.


How to display a layer's label source column?

I have a CSV layer that displays points on the map. I use one of the columns as the Label. In Print Composer I would like to display the name of the column used to generate layer Labels. Is this possible? **Update (additional information)** I have a CSV data source in this format:               Id        Lat        Long        PropA        ProcB        PropC        ...     1         11.1       33          456          45           10.3     2         11.2       33          123          75           13.4     3         11.3       33          345          85           14.5     4         11.4       33          678          95           15.6      At the moment, for each column PropA... _N_ , I do this:   1. Update layer Label to Prop ** _N_**.   2. Update text description label in Print Composer.   3. Save map as a PDF. I will be grateful for help automating these steps!


MapInfo cannot save in MWS format

I'm using MapInfo 11.0.4 Release Build 404 Whenever I want to save a workspace file in the MapInfo workspace format (*.mws) it saves the file as expected and with the right file extension. However, when I look at the content of the file in a text viewer, I can see that it is a .wor file in reality. It's not saving with the XML format that should be in a *.mws file. Since there is no other way to save workspace I cannot create the MWS format at all and I need it to use with the MapInfo Spacial Server product. Any help would be appreciated.


Using PROJ.4 library to transform from local coordinate system coordinates to global coordinate system using ground control points?

I have a point cloud whose coordinates are with respect to a local coordinate system. I also have ground control points with GPS values. Can I convert these local coordinates to a global coordinate system using PROJ.4 or any other library? Any code in Python for the above stated problem would be a great help.


gdalwarp not clipping neatline properly

I have a PDF file (10mb) which contains a NEATLINE polygon in metadata.               $ gdalinfo I5601_GeoPDF.pdf      > Driver: PDF/Geospatial PDF Files: I5601_GeoPDF.pdf Size is 5906, 3543 > Coordinate System is: PROJCS["UTM Zone 56, Southern Hemisphere", > GEOGCS["GDA94", ... GeoTransform = 128363.7573358897, 42.324133621968, > 0.8865642156959999 6471730.899290663, 0.8865642156959999, -42.324133621968 > Metadata: CREATION_DATE=D:20080515105609+10 CREATOR=BAE SYSTEMS MapMaker.exe > 2.9.830.19227 **NEATLINE=POLYGON ((191125.203093351592543 > 6340529.201165121048689,188350.332458852295531 > 6473000.141725543886423,378308.662319099530578 > 6476979.201503321528435,381083.532953590154648 > 6344508.260942890308797,191125.203093351592543 6340529.201165121048689))** > ... Now if I extract the NEATLINE data into a .cvs and use gdalwarp to cut it:               $ for i in *.pdf; do gdalinfo $i|grep NEATLINE | awk -F= '{print "id,WKT\n1,\"" $2 "\""}' > $i.csv; gdalwarp -crop_to_cutline -cutline $i.csv -overwrite $i $i.tif; done      I get a geotiff that isn't clipped correctly: picture of results I've tried this on US 250k maps and it works. I'm wondering whether anybody can tell me if I'm doing something wrong, or is the metadata in the PDF wrong ?


What are mystery green lines in Layout View produced during ArcObjects script execution?

Could any of ArcObjects guru tell me, what are these green lines on my Page Layout and what to do, or not to do, so they do not appear? See image attached. I do used ArcObjects from Python. My script is converting data frame grids and graticules into graphics. As result I receive Group Element from which I need to extract single elements (text, lines etc.) in order to play with them. However, after this operation (which is carried out for 3 separate grids) I receive those mysterious green lines, which are not part of GraphicContainer. Actually, green lines appear only after conversion and extraction from one- graticule (WGS). There is some guidance **here** about how to get rid of them, which works (cut and paste or delete and undelete all page elements). However, I would still like to avoid them in the first place. Also, after some configurations of cutting and pasting I figured that green lines do vanish when those 3 image elements (lower right corner) are deleted and restored. Code suspected to be responsible for this:               pGCL.Reset() #reset graphics containere     pElem = pGCL.Next()     while pElem:      pEP3 = CType(pElem, esriCarto.IElementProperties3)      if pEP3.Name == '' and pEP3.Type == 'Group':       pGE = CType(pElem, esriCarto.IGroupElement3)       pGE_no = pGE.ElementCount       for i in range(pGE_no):        if pGE.ElementCount > 0:         pGEe = pGE.Element[0]         pGCL.MoveElementFromGroup(pGE, pGEe, 0)        else: break      pElem = pGCL.Next()      I know there is some discussion going on, whether GIS.SE is appropriate place for such technical, software specific questions. I apologize for this, however, I did not found satisfactory answer anywhere else. ![Green lines in page Layout](http://i.stack.imgur.com/PjY4B.jpg)


Easy Extraction of features from a Feature dataset ( Not feature class)

My Feature dataset (residing inside a GDB) contains more than 30 feature classes and every month I like to extract features from this GDB based on one attribute, which is repeating in all the feature class : Eg : Job_Month = "May" Is there any tool available in ArcGIS, which will extract all the features in a Geodatabase ( or Feature dataset) at one time using an SQL . The extract data wizard does not have any filter/Sql for the attributes. Expected simple solutions, but seems not an easy one..


Mapserver + tilecache + openlayers attribute search

I'm using Mapserver + tilecache + openlayers for my web map application.   I'm using Mapserver as a WMS Server.   Now, I want to perform an attribute search, i.e. type something and the application searches for the text in the database (which is postgresql + postgis). The result should be shown in a table and also on the map with markers.   I have created a control for that, but I'm not sure how to perform the search and retrieve the results.   The solutions I'm currently thinking on:   1. Use GetFeatureInfo on Tilecache (although I'm not sure whether tilecache supports it)   2. Use GetFeatureInfo on Mapserver WMS   3. Write a Python application that performs the query and passes the answers to OpenLayers and I make the required adjustments to show it how I want.  Any idea will be greatly appreciated.


How to restore default values in Mapserver?

I'm using Mapscript reading a configuration of my map file, so i set some values in my mapObj. But in certain point of my code i need to restore the default values locate in map file. I tried to use the method "applyConfigOptions()", but didn't work. Is there a way to do that?


QGIS very slow startup on Mac

I am experiencing startup times of c5 minutes (!) for QGIS 1.8 on my Mac Mini with 2.4GHz processors and 8G RAM. I have just 4 python plugins and only a handful of shapefiles/data layers so I can not see why every start takes 5 minutes from clicking on file to having something showing on screen. Is this normal for QGIS or is there some mysterious solution? I have searched on Google, various forums (fora) and discussion boards etc but no solutions. Any help gratefully appreciated! Eric


Assigning Probability to cells

I am working with the image below and I am trying to find how many pixels are in the grey corners of the image, and I want to assign them a probability of 5%. The bottom layer is my DEM and the top layer is a raster of distances. I know how many pixels I have on my raster layer, however I am not sure of how I can get the pixels outside of the circles and how to assign those pixels a probability. I would appreciate any help I can get! Image Link


Why can't I access a table in a personal geodatabase (.mdb) with arcpy after installing 64-bit geoprocessing?

I wrote and ran a simple script last year - it cycled through some line features, buffered them and then intersected with some other lines - it all worked fine. I'm now trying to tweak the file to use a larger buffer, and it's not finding the database I am using to manage the loop and store the results. This code:               rows = arcpy.SearchCursor("C:\\CEDAR_GIS\\SPEEDY3_GPS_DATA\\shapefiles\\Analysis\\Intersections.mdb\\final_to_school_routes_LINES")       results in this error:               Traceback (most recent call last):       File "C:\CEDAR_GIS\SPEEDY3_GPS_DATA\shapefiles\Analysis\GISbuff-GPSline_IntersectionTS.py", line 29, in <module>         rows = arcpy.SearchCursor("C:\\CEDAR_GIS\\SPEEDY3_GPS_DATA\\shapefiles\\Analysis\\Intersections.mdb\\final_to_school_routes_LINES") #opens search cursor on Overlays_to_run       File "C:\Program Files (x86)\ArcGIS\Desktop10.1\arcpy\arcpy\__init__.py", line 1133, in SearchCursor         return gp.searchCursor(dataset, where_clause, spatial_reference, fields, sort_fields)       File "C:\Program Files (x86)\ArcGIS\Desktop10.1\arcpy\arcpy\geoprocessing\_base.py", line 359, in searchCursor         self._gp.SearchCursor(*gp_fixargs(args, True)))     IOError: "C:\CEDAR_GIS\SPEEDY3_GPS_DATA\shapefiles\Analysis\Intersections.mdb\final_to_school_routes_LINES" does not exist      But that table is in that database, and when I ran the script last year there was no error. I have since upgraded to ArcGIS 10.1 with SP1, which I think meant an upgrade to python. Has there been a change in the handling of .mdb files, or am I doing something stupid?


How to select features by matching attributes of two SHPs?

I have one line shape file with lot of road sections and a point shapefile, some of these points are outside my area of interest. The files have one attribute in common and I would like to do selection, where I select only those points that have matching attribute value with the line shape file. Is this possible in some meaningful way? * * * I believe that in this particular case I could do line buffers and use those to select the point in buffer area - but I was wondering if its possible to match the values "straight away" with each other.


Concatenate only when all fields to be concatenated have values

I have two fields that I want to concatenate. Field1 contains null values and Field2 has values for every record. I want to concatenate only when both fields have values. I tried using a function with an `if/else` statement but that doesn't work as the resulting concatenation will return values from field2.               def concat_fields(field1, field2):         if field1 == ""             return               else:             return field1 + field2      I know this is wrong. Any suggestions? EDIT: ![enter image description here](http://i.stack.imgur.com/smd57.png)


Which spatial DBMS should I select?

I found there are a number of RDMS have the spatial capability e.g. mysql, postgreSQL I considering which I should go forward! I would like to store the lat lon coord from the tracker to the db and make a query to find the nearest place or compare with the defined location. Could anyone give me an idea on this topic? Thanks a lot.


gdal_merge.py aligns incorrectly at the seam of two GeoTiffs

I am trying to merge four GeoTiffs from the National Elevation Dataset (EPSG: 4269) into a 2x2 DEM using gdal_merge.py. For some reason, when it merges the upper two NED rasters, gdal produces a seam between the two of approx. 1 px that it fills with incorrect values (either 0 or some large negative number). It seems like this is probably an alignment issue. How can I eliminate this artifact? Thanks for the help.


Problems combining (union) Multipolygons in geodjango

I'm using geodjango and postgis (1.x),   What is the best way to combine (union) a list of multipolygons. in what i assume is rather inefficient i'm looping trough like this               combined = multipolygon          for item in items:          ____combined = combined.union(item.geom)  #geom is a multipolygon      Usually this works fine, but often i'm getting the error Error encountered checking Geometry returned from GEOS C function "GEOSUnion_r". Here is the geo json version of the item the error is thrown on if it helps               { "type": "MultiPolygon", "coordinates":      [ [ [ [ -80.077576, 26.572225 ],            [ -80.037729, 26.571180 ],            [ -80.080279, 26.273744 ],            [ -80.147464, 26.310066 ],            [ -80.152851, 26.455851 ],            [ -80.138560, 26.538013 ],            [ -80.077576, 26.572225 ]      ] ] ]      }      does anyone have anyideas? the end goal is to take find all the locations (another table) which fall within this list of n polygons (using coordinates__within=combined_area) thanks in advance for any help.


OpenLayers - setting opacity ti a webservice layer (if statement)

I've got an openlayers application which calls the layers from MSSQL via a webservice. The problem is that I want to set an opacity level for one of the layers (If layer name = X, set opacity level to Y). How do I do that ? Thanks Lauren.


Looking for a simple method to find the shortest path with Shapefiles

I am having shapefiles for a particular map location, the shapefile consists of points and polylines. Now I want to find the shortest path between the points. What is perfect and simple method to find shortest path.


How to visualize pgrouting query results in mapguide open source

I would like to make a spatial query using postgis 2.0 and pgrouting 1.05 then visualize results in Mapguide open source.Any suggestions?


GeoServer Where clause filter

I use postgis, geoserver and openlayers for gis development. I simply reach from goeserver to postgis table or view and show it with openlayers. but I need some filtering operation on my spatial data. For example I have historical spatial data and filter them to start_date and finish_date like               select ....... from ...... where date between :start_date and :finish_date      the problem is start_date and finish_date are dynamic. User change them from web my open layer web project. How can I do it? from geoserver?from postgis or openlayers? thanks


Map algebra from QGIS python console?

What command do I use to do map algebra on rasters from the QGIS python console?


Beginning Arcpy: Asking for a quick review of my script

I cobbled this code together from ESRI documentation and stackexchange answers. I tried to use SearchCursor(points, fields=fieldname) but I get this error: > TypeError: SearchCursor() got an unexpected keyword argument 'fields' I'm > not sure why this doesn't work.                   # Local variables:         points =  arcpy.MakeFeatureLayer_management("C:\\alltogether.shp")         homebase_shp =  arcpy.MakeFeatureLayer_management("C:\\homebase.shp")                   fieldname = "Device"         Device = arcpy.AddFieldDelimiters(points, fieldname)         delimitedfield = arcpy.AddFieldDelimiters(points, fieldname)         devices = set([row.getValue(fieldname) for row in arcpy.SearchCursor(points)])              for i in devices:                 attr_sel = arcpy.SelectLayerByAttribute_management(points,"NEW_SELECTION", delimitedfield + " = " + str(i))                 totalpts = arcpy.GetCount_management(attr_sel)                      # Process: Select Layer By Location                 inside_sel = arcpy.SelectLayerByLocation_management(attr_sel, "WITHIN", homebase_shp, "", "SUBSET_SELECTION")                 insidepts = arcpy.GetCount_management(inside_sel)                      print "device %d has %d total points and %d inside the base" % (int(i),int(totalpts[0]),int(insidepts[0]))      a comment I saw earlier suggested doing it with fields= > Since you are only interested in one field, for better performance I would > specify the >optional fields argument, e.g. myList = > set([row.getValue(fldName) for row in >arcpy.SearchCursor(fcName, > fields=fldName)]) – blah238 Jan 29 '12 at 1:03 Select distinct values from a single column of an attribute table (or layer)


mysql database of bars, pubs and restaurants in the united states and maybe canada

I am mostly seeking Bars and Pubs but Restaurants are a very nice bonus. I need a database with zipcode, city, state, longitude, latitude, establishment name, establishment contact information, maybe a category with sub categories if it's restaurants as well. I found one with Oddity Software. But after closer inspection after buying it, the database was garbage. It had very long obtrusive category names and/or sub category names, that would make indexing it for any need a royal pain in the arse. Then it had so many entries that were not restaurants, pubs, or bars.. like auto repair centers, and a minimum of at least 80 categories/sub categories that made no sense to the type of database it claimed to be. All of which had multiple entires with that category/sub category that didn't apply. So I see all these new sites popping up here and there trying to be like yelp or doing something completely different like I want to do. However I don't know where these people are getting their initial databases from. I thought Oddity Software would be my golden fix until I purchased it, found out what I found out about it, then they refused to discuss things further with me. So I am wondering here today. Who knows of a database similar to that of which I seek?


Pros and Cons for developing GIS functionality in ArcGIS Web API (Silverlight) or ArcGIS for Desktop

I have an application coming together using the ArcGIS Server API for Silverlight. The app, currently, takes in two geometry points from the user and passes the geometry to an ArcGIS Elevation server online to get elevation values. Then the SL client sends the data over to an internal ArcGIS Server to do custom slope calculations and linear referencing on a semi-static polyline with a gp service. What are pro and con arguments to move the UI from silverlight over to desktop application. Perhaps a C# extension/addin that allows the user to select points similar to how they interact in the web map. And documentation or white papers on the subject. This would be strictly an ArcMap centric dev question. Now it is an internal app. But deployment issues would be a good topic. As well as maintenance and complexity of development. Any others points would be very welcome as well. I guess this may not be a strict question to get an answer to maybe it needs to be moved?a Thanks.


Dynamically adding WFS layer features with Geoserver

I am investigating using Geoserver to output new WFS layer features dynamically. We have a backend process that will be generating new WFS features periodically that will be placed in a datastore (most likely PostGIS). I plan on creating a WFS layer for the feature data in PostGIS. The idea is to have the backend process dynamically add new features into PostGIS and my expectation is that clients making WFS features requests using Geoserver will get the new feature data I am adding at runtime. Is this the case? I can't seem to find any documentation discussing what happens when the datastore's data is updated.


Does CopyFeatures from in_memory via Python script tool cause General function failure for you?

The Python code below looks pretty straightforward.               import arcpy          if arcpy.Exists("D:/temp/test.gdb"):         arcpy.Delete_management("D:/temp/test.gdb")          if arcpy.Exists("in_memory"):         arcpy.Delete_management("in_memory")          arcpy.CreateFileGDB_management("D:/temp","test","CURRENT")          arcpy.CreateFishnet_management("D:/temp/test.gdb/test","0 0","0 1","1","1","1","1","#","NO_LABELS","#","POLYGON")          arcpy.CopyFeatures_management("D:/temp/test.gdb/test","in_memory/test2","#","0","0","0")          arcpy.CopyFeatures_management("in_memory/test2","D:/temp/test.gdb/test3","#","0","0","0")      and when I run it on a PC (ArcGIS 10.0 SP2; Windows XP SP3; 32bit) it completes in 4-5 seconds. This is the case when I run from IDLE or as a Python script tool in ArcMap (where I added it to a toolbox with no parameters). However, when I move it to a virtual server (ArcGIS 10.0 SP3; Windows Server 2008 R2 Standard; 64bit) it still runs fine from IDLE, but when run as a Python script tool in ArcMap it just sits at:               Executing: test     Start Time: Thu Mar 15 10:31:50 2012     Running script test...      When I move the tool I find a small Error dialog with a white "i" in a blue filled circled. To dismiss that dialog I click the OK button and/or the "x" at its top right. When it goes, a second Error dialog appears - this time with the text below:               General function failure [test3]      The Python test code provided here has been stripped out from a much larger tool that I need to publish as a Geoprocessing Service, which in turn requires that I save it as a tool layer within a map at 10.0 (10.1 won't be an option for quite a while). The culprit seems to be the CopyFeatures from the in_memory workspace which will take some working around, and if I need to eliminate in_memory from this polygon overlay intensive application, then its performance is going to take a big hit. **Would someone be able to run the same test as I describe above on a virtual server (it may also manifest on a physical server but I don't have one to test on) running ArcGIS 10.0 and Windows Server to see if they can reproduce what I am seeing, please?** If it works, or if it does not, can you specify your ArcGIS and Windows versions and whether you ran it on a physical or virtual server, please?


Exposing ArcInfo level liscenses to ArcView desktops via ArcGIS Server and custom python script tool

I have a script that uses an ArcInfo level tool. The script works great on my machine but uses an ArcInfo level tool which the ArcView users don't have. Is exposing these ArcInfo level tools possible throught ArcGIS server? What I am having trouble conceptualizing is how my python script tool will call that gp service. Please help me think through the liscensing and calling the gp service from a custom python script tool on the toolbox. I have ArcInfo, I added the tool to arcmap's TOC saved the mxd and opened the tool and ran it once. I then published the tool throught ArcGIS Server. I am stumbling on how to proceed. Thanks.


How to calculate line centroids for 460,000 features?

I have a large shapefile of lines, and want to run a mmqgis function (Plugins --> mmqgis --> Modify --> Convert Geometry Type.) The problem is the file is very large and crashes the program. What options do I have? Would running it as background python script work (and if so how do I find the commands for mmqgis function?)? There are 460'000 features and i'm running a 3.40GHz i7 with 12GB of RAM. Is there a way of this not breaking down?


Near / Near Table creates -1 values

I have two feature classes.   1. Point feature of hand-sampled GPS locations   2. Line feature of terrace walls digitized from aerial photograph. ![sample points and terrace walls](http://i.stack.imgur.com/1rhft.png) Because of the less than perfect ways I generated these features, I need to allocate each GPS point to the nearest terrace wall from which I sampled it. I defined the Input as the Point Feature, and the near feture as the terrace lines. ![near tool](http://i.stack.imgur.com/X8Pmq.png) But the process fails to generate results. If it's a Near Table, then it's an empty one, and when I use the Near Tool itself on the point feature, it adds two fields with -1 in them. ![Point feature with -1](http://i.stack.imgur.com/lO7oZ.png) When I run the tool on the features in an opposite way, it does generate results, and for every terrace walls gives me the nearest point - obviously, not what I'm looking for. Any idea why the process fails? HDunn.


How to add Shapefile to map?

How to Add Shapefile to map in ArcGIS Silverlight application runtime ? I published the shapefile feature as map service on arcgis server and consumed it my SL Application through the service Url but the shapefile feature does not added to the map .. what's wrong?


How to open *.aux file in QGIS?

My university desires to switch from using ArcGIS to using QGIS and therefore I have to figure out whether the exercises that are used in lectures can be performed using QGIS. I have an .AUX file that can be used in ArcGIS. When I want to open this file in QGIS, the file shows up in the layers list but no image is drawn on the canvas. Also, there is no attribute table. In ArcGIS, there is an image drawn and there is an attribute table. What could cause these differences and how to overcome this problem? Thanks


ArcGIS - Dimension Feature problems

I would like to use the Dimension Feature tool in ArcGIS - this button ![enter image description here](http://i.stack.imgur.com/Qt5hE.jpg) as I want to have a label on every side of some polygons (they are all reasonably simple polygons). When I click on my polygon to add the dimensions (and don't move or click the mouse again), they come up like this which is how I want them: ![enter image description here](http://i.stack.imgur.com/hdWfQ.jpg) However as soon as I click the left mouse button to place them, they go like this: ![enter image description here](http://i.stack.imgur.com/VpcTD.jpg) Am I missing an option somewhere? it does not allow me to select 'Dimension Feature' in the Construction tools section of the create features dialog. Its only available as a button (that I found by adding the button manually). The closest construction tool I can select is Dimension Edge, but that means I have to add dimensions for each edge of the polygon. Its a slower process.


Installing ArcGIS 10.1 on SQL Server 2008 R2 Cluster

We're running into some performance issues when accessing large enterprise geodatabases on a SQL Server 2008 R2 cluster. The database itself performs well and returns queries fine but presenting the information back to ArcGIS for Desktop 10.1 is prohibitively slow. Connection is direct connect. Data type in SDE is DBO. Machines are VMs with 16Gb RAM. Are there any tips or tricks or known documents for installs of these on SQL Server clusters where a node instance is presented? Thanks Doylersa


Arcgis IPolygon(System.__comObject) and ILine2 (Esri.Arcgis.LineClass) intersection

Currently i am developing a Server Object Extension. On my project i am creating a new line class with LineClass myLine = new LineClass(); and getting a polygon from client as JsonObject than converting this polygon to IPolygon with               JsonObject jsonFieldPolygon;     if (!operationInput.TryGetJsonObject("field", out jsonFieldPolygon))          throw new ArgumentNullException("field");     IPolygon4 fieldGeom = (IPolygon4)Conversion.ToGeometry(jsonFieldPolygon,esriGeometryType.esriGeometryPolygon);      Polygon i have converted from jsonobject is becoming a com object, but line i have created is a class in my own application space. Thus i cant do operations with these two geometries, i guess both must be com object for operations like intersect and etc. I am getting an access violation when i try to intersect these 2 geometries. my code is below               JsonObject jsonFieldPolygon;     if (!operationInput.TryGetJsonObject("field", out jsonFieldPolygon))         throw new ArgumentNullException("field");     IPolygon4 fieldGeom = (IPolygon4)Conversion.ToGeometry(jsonFieldPolygon, esriGeometryType.esriGeometryPolygon);          JsonObject jsonFieldPath;     if (!operationInput.TryGetJsonObject("fieldPath", out jsonFieldPath))         throw new ArgumentNullException("fieldPath");     IPolyline fieldPathGeom = (IPolyline)Conversion.ToGeometry(jsonFieldPath, esriGeometryType.esriGeometryPolyline);          ITopologicalOperator myOp = (ITopologicalOperator)fieldGeom;          IGeometry resultGeom = myOp.Intersect(fieldPathGeom, esriGeometryDimension.esriGeometry1Dimension); // this line works because both geometries are com objects               IPoint startPoint = new PointClass();     IPoint endPoint = new PointClass();     LineClass myLine = new LineClass();          myLine.SpatialReference = fieldGeom.SpatialReference;     endPoint.SetEmpty();     startPoint.SetEmpty();          startPoint.PutCoords(fieldGeom.Envelope.XMin, fieldGeom.Envelope.YMin);     endPoint.PutCoords(fieldGeom.Envelope.XMax, fieldGeom.Envelope.YMax);          myLine.PutCoords(startPoint, endPoint);               IGeometry resultGeom2 = myOp.Intersect(myLine, esriGeometryDimension.esriGeometry1Dimension);// this line throws accessviolation exception because one geometry is com other is a class in my app      I hope there is a way to solve this


creating geotiff in ArcGIS to publish in geoserver

Pls give some answer, **Short question** How create geotiff image (to publish using Geoserver)? **Bit more information** I want to publish geo-reference image in geoserver.Then i chose image format as geotiff and create using arcGIS 9.3.It create tiff,tfw and aux files.But when publish using geoserver it gives this message, "Could not list layers for this store, an error occurred retrieving them: Unable to acquire a reader for this coverage with format: GeoTIFF" What is the error,I think this is not the correct method to create geotiff(to using GS).True? If answer is yes What is the correct method? Thanks for your time


Mass-update the projection in compressed GeoTIFF files

I have a huge set of compressed GeoTIFF files. Someone forgot to set the projection while creating the GeoTIFF files and I would like to add it now. If I use > gdalwarp -t_srs source.tif target.tif the result is an uncompressed file. I would like to avoid the (de)compressing if possible, since I expect it will degrade the final result due to the recompression and because it will require more time to decompress and compress each file again. I assume that the projection information is stored in the header. **Is there an easy way to just update the projection information in a GeoTIFF file?**


Openlayers requests to geoserver don't produce cached tiles

I've created a grouped layer made from 4 MrSID rasters. I have a Openlayers map with basic Google layer. When I preview the grouped layer in geoserver , I see that tiles are being created in the GWC folder (tried both EPSG:4326 and EPSG:900919). When I view the layer in my openlayers application (on same localhost machine) , no cached tiles are being created (I've cleared the browser cache to make sure the tiles are being generated from geoserver). In my Openlayers I'm using GeoExt2 ... but still its regular WMS call. Any idea why no tiles are being generated ?!


How to copy a group of layers between QGIS projects?

In ArcGIS I am used to being able to have a couple of MXD documents open and copy-paste groups of layers from one MXD's table of contents to another. This can save a LOT of time ensuring consistent styles, but also in setting up the drawing order, pulling in the data, defining its CRS etc. for the group I would like to do the same thing in QGIS. I have two QGS documents open but there is no apparent option to save or copy a group of layers and dragging/dropping doesn't work (I get a little 'No Entry' symbol). Is this possible and, if so, how? Many thanks


Can I update my Goole Earth browser in real time?

I am about to start a project that will entail updating my Google Earth browser with data from a MySQL database package (WAMP) like is the case for web pages, but after thinking about it I really did not figure it out, so I was wondering if this is possible in real time/ I will be glad with any constructive contribution Thanks in advance. Thanks for your reply, I really realised how vague it was, thanks for pointing that out, you are awesome. For my data flow: * I will like to insert geo-coordinates on html via php to MySQL(using float INT, is there a better solution?) * Then retrieving the data and feeding it in Google Earth in My Places. So far I have done so but upon the addidtion of new data, the kml file must be loaded again before featuring in Google Earth. **My problem now is to know if an automatic browser refresh is possible, so that I won't have to load every kml file after every new data is fed into the database.** As for my SRS; I (Admin) will be using WAMP on Win 8, Php 5.3, MySQL 5.5. The users will input via HTML (frontend); and Google earth as my backend browser. The end point is I will like to be able to locate those in my network in a given perimeter with their geo-coordinates, and dynamically too. Hope I have given enough useful information, thanks in advance.


Hidden Features of ArcGIS

ArcGIS is one of the most popular software-packages available for GIS. What are the lesser-known but useful features of the ArcGIS? Which little/hidden tricks do you known that make your experience with this software better?


I want to edit corresponding postgresql data of a vector from openlayers- how to do it?

I want to edit the corresponding text/numeric data type columns of a vector through openlayres, when a user click/mouse over a vector. like, http://dev4.mapgears.com/bdga/bdgaWFS-T.html# any tutorial page for how to do it? I gone through it , it was difficult for me to understand from its js codes. -- asked in stackoverflow, but no answers. based on a comment, i posted the same question here.


How do you view multiple tables at the same time in ArcMap 10?

My original question was: > How do you undock the tables in the ArcMap attribute window? because that's how I thought it should work in the new ArcGIS 10 docking environment. Either it isn't intuitive or I am missing something obvious. I have an ArcMap user who went up to ArcGIS 10 yesterday and he wants to view two tables at the same time. If he opens two tables, they are loaded as tabs in the table window. He(we) try to drag one tab elsewhere and nothing happens. The only window control button is an 'X' which closes the table. If I right-click on the tab, **Close** is the only option.


Open Layer and Vector seems not working

I configured geoserver (2.4.4) with a layer in EPSG:4326 projection. In my client side I added the openLayer code to read and show this layer as in the following:                  map = new OpenLayers.Map ("map", {                     controls:[                         new OpenLayers.Control.Navigation(),                         new OpenLayers.Control.PanZoomBar(),                         new OpenLayers.Control.Permalink(),                         new OpenLayers.Control.ScaleLine({geodesic: true}),                         new OpenLayers.Control.Permalink('permalink'),                         new OpenLayers.Control.MousePosition(),                                             new OpenLayers.Control.Attribution()],                     maxExtent: new OpenLayers.Bounds(-20037508.34,-20037508.34,20037508.34,20037508.34),                     maxResolution: 156543.0339,                     numZoomLevels: 19,                     units: 'm',                     projection: new OpenLayers.Projection("EPSG:4326"),                     displayProjection: new OpenLayers.Projection("EPSG:4326")                 } );        .....        var kmlLayer = new OpenLayers.Layer.Vector("KML" ,{                     strategies: [new OpenLayers.Strategy.Fixed()],                     protocol: new OpenLayers.Protocol.HTTP({                         url: "http://localhost:8080/geoserver/wms/kml?layers=avl:mobiles"},                         format: new OpenLayers.Format.KML({                             extractStyles: true,                             extractAttributes: true,                             maxDepth: 2                         })                     })                 });              map.addLayer(kmlLayer);      Unfortunately I'm not able to see anything for the overlay layer. Instead if I add the folowing layer                   var avl = new OpenLayers.Layer.WMS("AVL Mobiles",                                        "http://geoserver-1:8080/geoserver/wms",                                        {                                            layers: "avl:mobiles",                                            transparent: true                                        }, {                                            opacity: 0.7,                                            singleTile: false                                        });         map.addLayer(avl);      all seems to be OK. I need to use a Vector layer insted of a WMS one, as I need to add the refresh strategy and WMS layer doesn't support it. Any idea why nothing is showed in map? Any help will be apreciate.   Thanks in advance.   I discovered that the problem could be on request OpenLayer sent to geoserver. In fact I found this exception in geoserver log:               org.geoserver.platform.ServiceException: GetMap requests must include a BBOX par     ameter. but I was not able to add BBOX parameter to my vector layer.


How to import PostGIS layer with geography type to QGIS?

I tried to import in QGIS a layer from a PostgreSQL table with a geography type column created through the extension PostGIS. I noticed that QGIS doesn't seem to like the geography column... Here the step to create the column... `ALTER TABLE traffico ADD COLUMN geog geography(POINT,4326);` then to update the column... `UPDATE traffico SET geog = ST_GeogFromText('SRID=4326;POINT(' || lon || ' ' || lat || ')');` but then when I connect QGIS to my server I see this: ![Import PostGIS layer to QGIS](http://i.stack.imgur.com/56OSw.png) There's no way I can find to tell QGIS to treat my column as a geography type instead of geometry type. Indeed if I try to upload the layer like this the point coordinates do not correspond to reality.


How to use ArcPy cursor to define each feature sequentially as current selection of layer?

I want to implement a script that runs through each feature of a feature class, zooms to the extent of the feature, sets the scale to 1:5.000 and adjacent takes a picture of the scene. The steps are in general very easy. My idea is to use the following code lines:               import arcpy.mapping     myMXD = arcpy.mapping.MapDocument("CURRENT")     df = arcpy.mapping.ListDataFrames(myMXD)[0]     layer = arcpy.mapping.ListLayers(myMXD, "", df)          # Zoom and Scale     df.zoomToSelectedFeatures()     df.scale = 5000      My Question is: How can I run through the whole list of features so that the script will zoom on each feature, set the scale back to 1:5000 and takes a picture? My first idea was to use the Search Cursor to run through the attribute table, but how can I than declare that the feature that is actually "visited" by the cursour is the selected one? Thanks for any help and suggestions. Regards


netCDF does not display properly in QGIS 1.8

I have a netCDF file of fire emissions on a regular grid of 12 km by 12 km grid cells with 200 east-west grid cells by 161 north-south grid cells centered on the continental United States. The map projection is Lambert Conformal. I am using QGIS v1.8. When I open the raster it displays as a very, very tiny raster located in mid-northern reaches of the State of Kansas. Can someone tell me what blunder I am performing?


DLL problem with PyQGIS

I am having problems with Python Shell importing the module "qgis.core". When I type "import qgis.core" the Shell ("idle.pyw") gives me this error:               Traceback (most recent call last):       File "<pyshell#0>", line 1, in <module>         import qgis.core     ImportError: DLL load failed: Impossibile trovare il modulo specificato (Trad "Impossible to find the specified module").      I have already set the environment to point the right folders following "PyQGIS cookbook" instructions. In my case, the paths are:               PYTHOPATH=C:\"QGIS_path"\apps\qgis\python;     Path=C:\"QGIS_path"\apps\qgis\bin.      Am I missing something? What are supposed to be the right libraries the Shell is unable to find? I am working with windows 7 64 bit, QGIS Lisboa (version 1.8), Python 2.7.2 (the one that is included in the QGIS package).


How to create a routing application with ArcGIS Server?

I have the task to find out what is necessary to build a routing application using ArcGIS Server. I am not a GIS beginner, but have had no experience with ArcGIS Server. As this project is at the "ideas" stage, I would like to start my research here by asking if anyone could give me some guidelines regarding what would be a sensible way to start a routing project using primarily (but not necessarily limited to) ESRI Software. Which extentions? What type/quality of data? Planning tips? Which organisations etc etc?


How to rename .met file with model builder which uses an iterator to rename .sid files?

![iterator for renaming .sid rasters](http://i.stack.imgur.com/lvxen.png) I have a model that iterates through a folder containing .sid raster datasets, amongst other types of rasters (so I'm using a *.sid wildcard). I note, however, that ArcCatalog/Model Builder/Rename tool cannot rename the .met file of this .sid dataset, hence the renamed dataset does not have projection information. Does anyone know how to force model builder to rename the .met component of a .sid file? Possibly I'd have to save it as a Python script and add some code to that, if someone can show me. Wil


Good advice regarding Mobile GIS Hardware and software for outdoor tasks?

Could anyone offer advice regarding outdoor mobile hardware and software for GPS-photos, data capture, and textual input. We don´t want to necessarily write long text, so the ability to write into shapefile format would be optimal. Obviously the price is a big factor and we have a budget of around 2-3,000 Euro/Dollar for everything. The colleague who needs this also has an iphone...so this could also be an option. I know that this has been covered in some extent in other posts, but I also would like to know if anyone has used pure mobile GIS opensource software on Android or Apple based operating systems. The Ipad also has GPS build in and would also be a great platform for GIS mobile tasks.


Spatial analysis and clustering near features

Im working on marking behaviour in group living animals and im interested in how marking behaviour is affected by certain characteristics of the neighbouring groups. I have plotted the territories of each group from 95% density isopleths created using a-LoCoH. The marks of each individual from a group are plotted as a separate points layer. Marking behaviour does not appear to be random, and they do seem to cluster around boundaries with neighbouring groups. What I would like to investigate whether individuals mark more near some groups than near others, for example, do individuals mark less near groups that are closely related? is it possible to show that points are clustering near some features but not near others? Any help that can be provided would be greatly appreciated, Dave Humphries ##### update For demonstration, below is an image of four territories (A,B,C, and D) and the marking locations of two individuals from group A. Groups A and B are closely related. Groups C and D are unrelated to group A. ![enter image description here](http://i.stack.imgur.com/eE7TD.jpg)


calculating terrain elevation line

I found too many hits on the web to even see the right direction to go. I have to provide a small proxy server serving requests for a line of terrain elevation data between two coordinates. It is a Linux system where GDAL is installed already. There are standard raster elevation maps as inputs. The proxy needs to calculate about 1000 points/sec at runtime. I'm not that much familiar with python but what I've read so far is that GDAL/python seems the way to go. Is my feeling right to go this way? thanks for feedback Wolfgang R.


Styling map elements based on their properties

I am currently generating tiles based on OSM data using Mapnik and Python and am quite happy using the XML styling but I am at a loss as to how to vary the styling on an element by element basis. I would like to render OSM map tiles where certain elements are styled according to their attributes, for example, a roads colour is based on its length, or an area is coloured based on its size. Is this possible?


How to find polygon neighbors by rook and queen contiguity in a feature class?

I have a feature class containing several polygons. I would like to find all polygons that have a common border to a polygon of my choice. This can be done with the ISpatialFilter and esriSpatialTouches. This is known as the queen contiguity. However, the rook contiguity does not contain the polygons, which have only a point in common. It seems that this is something I can not figure out with the DE-9IM. My idea is to test all polygons which touches my selected polygon with the Intersect method from the ITopologicalOperator interface and see if the result is from the type Point. If this is the case I have to remove this polygon. Is there a better way to achieve this?


How can I change the size of labels and points in GeoServer?

I am beginner in GeoServer. I was able to do some task by my professor they want me to change different size of points in map of GeoServer example I have my database named "PLACES" then these are the data inside the database.               Name          Place     Addr:city        Baluarte                Vigan     Ayusan Norte            Vigan     Barangay 8    Village   Vigan     Chinos                  Los Baños     Casubha                 Los Baños      Then I want to enlarge the labels and points of addr:city before the small cities in the field name how can I do that is that possible?


How to find out the Basin Length?

I am doing Project on Drainage basin, to find out the Elongation ratio I need to find out Basin length, Can any one of you tell me how to find out BASIN LENGTH using ArcGIS or any other way. My data base is like below Image ![enter image description here](http://i.stack.imgur.com/8zXQa.jpg) (i just uploaded this image from internet)


How to create a heatmap using values?

How do I create a 'heatmap' in QGIS based on the values of the different points? I did do this once before with rainfall, but can't remember how!


Impossible to open the .qgs file I've just saved with QGIS 1.8

I'm using QGIS 1.8.0 on Mac OS X 10.6.8 (last stable version of QGIS I downloaded yesterday). After working yesterday and this morning on it, I'm trying to open it again, and QGIS asks me the username and password to access the PostGIS database and crashes after 10 seconds. I also have QGIS 1.7 and 1.9 installed on my computer. Here is the crash message, if this can help:               Process:         QGIS [342]     Path:            /Applications/QGIS.app/Contents/MacOS/QGIS     Identifier:      org.qgis.qgis     Version:         1.8.0 (1.8.0 ())     Code Type:       X86-64 (Native)     Parent Process:  launchd [124]          Date/Time:       2012-07-13 14:54:37.581 +0200     OS Version:      Mac OS X 10.6.8 (10K549)     Report Version:  6          Interval Since Last Report:          253852 sec     Crashes Since Last Report:           18     Per-App Interval Since Last Report:  15946 sec     Per-App Crashes Since Last Report:   7     Anonymous UUID:                      4D1290E1-470F-4017-A8E7-1EA41258D661          Exception Type:  EXC_BAD_ACCESS (SIGSEGV)     Exception Codes: KERN_INVALID_ADDRESS at 0x0000000000610163     Crashed Thread:  0  Dispatch queue: com.apple.main-thread          Thread 0 Crashed:  Dispatch queue: com.apple.main-thread     0   libdelimitedtextprovider.so     0x000000011a2edeaf QgsDelimitedTextProvider::~QgsDelimitedTextProvider() + 35     1   libdelimitedtextprovider.so     0x000000011a2ede62 QgsDelimitedTextProvider::~QgsDelimitedTextProvider() + 18     2   org.qgis.qgis_core              0x0000000103e2d5b4 QgsProviderRegistry::provider(QString const&, QString const&) + 192     3   org.qgis.qgis_core              0x0000000103e50e2b QgsVectorLayer::setDataProvider(QString const&) + 89     4   org.qgis.qgis_core              0x0000000103e5b825 QgsVectorLayer::readXml(QDomNode const&) + 463     5   org.qgis.qgis_core              0x0000000103dfc763 QgsMapLayer::readXML(QDomNode const&) + 3999     6   org.qgis.qgis_core              0x0000000103e1e815 QgsProject::addLayer(QDomElement const&, QList<QDomNode>&, QList<QPair<QgsVectorLayer*, QDomElement> >&) + 1291     7   org.qgis.qgis_core              0x0000000103e1c910 QgsProject::_getMapLayers(QDomDocument const&) + 1378     8   org.qgis.qgis_core              0x0000000103e1f908 QgsProject::read() + 3760     9   org.qgis.qgis_core              0x0000000103e1ea23 QgsProject::read(QFileInfo const&) + 77     10  org.qgis.qgis                   0x00000001000240f4 QgisApp::addProject(QString) + 140     11  org.qgis.qgis                   0x0000000100026525 QgisApp::openProject(QString const&) + 55     12  org.qgis.qgis                   0x000000010001622e QgisApp::openFile(QString const&) + 108     13  org.qgis.qgis                   0x00000001000169ac QgisApp::event(QEvent*) + 40     14  QtGui                           0x0000000104d0193d QApplicationPrivate::notify_helper(QObject*, QEvent*) + 189     15  QtGui                           0x0000000104d07dc4 QApplication::notify(QObject*, QEvent*) + 1060     16  org.qgis.qgis_core              0x0000000103da0762 QgsApplication::notify(QObject*, QEvent*) + 100     17  QtCore                          0x0000000104a4a17c QCoreApplication::notifyInternal(QObject*, QEvent*) + 124     18  org.qgis.qgis_core              0x0000000103da0ada QgsApplication::setFileOpenEventReceiver(QObject*) + 158     19  org.qgis.qgis                   0x000000010000bf52 QgisApp::QgisApp(QSplashScreen*, bool, QWidget*, QFlags<Qt::WindowType>) + 5508     20  org.qgis.qgis                   0x0000000100008f52 main + 7938     21  org.qgis.qgis                   0x0000000100006c04 start + 52          Thread 1:  Dispatch queue: com.apple.libdispatch-manager     0   libSystem.B.dylib               0x00007fff895ffc0a kevent + 10     1   libSystem.B.dylib               0x00007fff89601add _dispatch_mgr_invoke + 154     2   libSystem.B.dylib               0x00007fff896017b4 _dispatch_queue_invoke + 185     3   libSystem.B.dylib               0x00007fff896012de _dispatch_worker_thread2 + 252     4   libSystem.B.dylib               0x00007fff89600c08 _pthread_wqthread + 353     5   libSystem.B.dylib               0x00007fff89600aa5 start_wqthread + 13          Thread 0 crashed with X86 Thread State (64-bit):       rax: 0x000000011a2f34a0  rbx: 0x000000011fbe61b0  rcx: 0x00007fff710b7650  rdx: 0x00000000005e8300       rdi: 0x0000000000610163  rsi: 0x000000011fb00000  rbp: 0x00007fff5fbfdee0  rsp: 0x00007fff5fbfdec0        r8: 0x0000000000000001   r9: 0x000000011fbe67b0  r10: 0x0000000000000000  r11: 0x0000000000000015       r12: 0x000000011fbe1fd8  r13: 0x000000011fbe1f10  r14: 0x000000011fbe5180  r15: 0x000000011fbe1f50       rip: 0x000000011a2edeaf  rfl: 0x0000000000010206  cr2: 0x0000000000610163


Can I collect Seismic data with an array of cell phones?

Suppose I have an array of cell phones, each having:   * an accelerometer    * a clock   * a GPS receiver   * Windows Phone 7 (WP7) How hard would it be to write a program so that synchronizes the phones to a common clock (preferably a GPS clock)? How hard would it be to write a program that starts recording accelerometer readings when it receives text message? I would want it to collect for a duration specified in the text message, and then forward the collected data (including location) back to an address specified in the text message. How hard would it be to write a program (on a PC or a cell phone) that recieves the data, processes it to compute a "fix" (ellipse) on the source signal location, then sends that fix back to each cell phone as a text message? Alternatively, could I record audio spectrograms instead of accelerometer readings?


Most timely source of California Fire Map data?

Fire season always has me scouring the web for satellite fire map information. Unfortunately many of the locations I find are outdated or not maintained well. What would be the best site to go to for the most timely fire data?


Is there any problem with Service Pack 4 for FME?

Let me try to explain a problem. I am testing exactly the same workbench on the to computers. One has installed FME 2012 without SP, another has the same but with SP4. Tested workbenches behaviour are different. First, one without SP4 process data min. 10x faster. Other thing is difference in information about processing in log window. For example, for FME 2012 without SP:               Reading table/feature class 'name1'      Reading table/feature class 'name2'      Processed CANDIDATE 7070      Reading table/feature class 'name33'      Processed CANDIDATE 19460      Processed CANDIDATE 27470      And for FME 2012 SP4 at the same part of log:               Sorting 566 features...     Calculating inside points...     Determining uppermost containment...     Constructing donuts from 566 features...      And that takes long time. Input data are in File geodatabase API format. Any clues what should I do? Why constructing donuts and holes? I was quite satisfied with the way and processing time before. Thank you.


FME - process only matching name datasources

I'm trying to find out how to deal with the below scenario with FME 2012 on board. I have a bunch of geotiff's tiles (A) and a bunch of others geotiff's tiles (B) that I would like to mosaick together. The name of each tile from (A) match eaxactly one tile in (B). I would like the process to be run tile by tile due to the amount of the tile I have to process (roughly 13k). What whould be the scenario with FME to: - take a pair of tiles with exact names from A and B, - run a workbench, - output the result, - take another pair of tiles with exact names from A and B, ... and so on. The problem is that I actually don't know how to force FME to read a tile A then find a coresponding tile in B (not reading all tiles from B) and run a workbench for the matching pair. Thanks in advance for help with that. G.


How to select all parcels of a street? [ArcGIS Server]

I have an Arcgis Server Map Service with a `Street` layer and a `Parcel` layer, I need to get the info of each parcel of a street that is selected by the user, does anyone knows how to do that? I'm a little bit lost. I'm using Arcgis Server 9.3, Visual Studio 2008 with Web ADF for .NET, thanks.


Importing arcgis.rand error from python shell

Not sure what's going on here with this module. Maybe a pathing issue?   1. If I import the module from the python window in ArcMap it works.   2. If I import it in a script, and use that script in a scriping tool the tool works.   3. If I import it in a script and run the script, it fails. Error message says no module named arcgis.   4. If I call the script tool from a Python add-in tool it fails and gives the same error message. This worked on my computer at home but does not work on my computer at work. Any ideas? I'm using ArcGIS 10.1 and Python 2.7.2.


How can I get the map to repeat to fill the screen?

I'm programming my first ArcGIS application. I've created a simple map. When I zoom out to a world view, the map doesn't fill the width of the screen. Ideally I'd like it to repeat until the available space is filled. How can I do this? Here's the code I have so far:               var map = new esri.Map("map", {          basemap: "topo",     center: [0, 0],     zoom: 2     });      And here's the result: ![enter image description here](http://i.stack.imgur.com/exJB8.png) How can I get the map to repeat to fill the white space on either side? (Crossposted: http://forums.arcgis.com/threads/79757-How-can-I-get-the-map-to- repeat-to-fill-the-screen) **Update** : I've set the spatial reference to (I think) Web Mercator, and wrapAround180 to true, but this has no effect. New constructor:               var map = new esri.Map("map", {                 basemap: "topo",                 center: [0, 0],                 zoom: 2,                 sliderPosition: "bottom-left",                 spatialRefence: new esri.SpatialReference({ "wkid": 3857 }),                 wrapAround180: true             });      What am I doing wrong?


Bounded box on postgis

I have a few MULTIPOLYGON geometries that are very complex to work with, and I need to find POINT geometries from a ~40M rows that are inside these MULTIPOLYGON geometries.   Right now I'm working with ST_Within but things are going on extremely slowly. I don't need all of the POINT geometries, just a few of them, but they all have to be inside of the MULTIPOLYGON (I need a few points per geometry). I thought that if there was something reverse to an envelope, that is to find the largest box that is within a geometry. I couldn't find it in the docs or Google or here. Is it possible to create such a box? If so, how? I do not want to use simplify since it is:   1. Not accurate enough.   2. Still slow.


How to make polygon from cells of a raster?

I have 100 rasters, each raster has the same cell size and extent and i want to make a polygon of each cell. Can anyone help me doing this?


MapFish Print: does it support ArcGIS Server?

we are looking for a map print solution, MapFish printing module quite fits. but as far as the searches I've done, it seems MapFish doesn't support ArcGIS Server. One thing I am not sure about is, whether it supports ArcGIS Server WMS or not? wish someone in this site who knows this can give me some advice. thanks!


Random value calculation in QGIS

Is there a way to do a random value calculation within a range in QGIS? I have been able to do this within ArcMap Field Calculator with the following code:               Dim max, min     max=20     min=5     x=(Int((max-min+1)*Rnd+min))          __esri_field_calculator_splitter__     x      Thanks


How to store and edit data with a defined topology using open software and formats?

I'm looking for a way to store and edit spatial data that preserves topological information. In my concrete application, I need to represent triangular meshes, which consist of unique nodes and triangular polygons. So far, I have been working with two separate shapefiles, one for the nodes and one for the triangles. This has the drawback that the two aren't topologically connected: If I move a node, the affected polygons in the other shapefile can't be updated accordingly. Is there a standard open format for storing and editing spatial data of different types (points and polygons) that are connected topologically, similar to what I understand ESRI Coverages to be?


How to migrate a Geoserver setup to Tomcat?

I have some problems with tomcat/geoserver. I was create project in stand- alone geoserver and all works fine. This is my code.               <%@ page language="java" contentType="text/html; charset=UTF-8"      pageEncoding="UTF-8"%>       <!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN"        "http://www.w3.org/TR/html4/loose.dtd">                    <meta http-equiv="Content-Type" content="text/html; charset=UTF-8">      <title> My first JSP </title>      <link href="http://localhost:8080/geoserver/openlayers/theme/default/style.css"type="text/css"/>

                But when I try to put this in tomcat nothing happens and I don't see a map. Im gonna do something with tomcat to run project which uses the Geoserver? PS im know my english is terrible.


Use a WMS raster in ArcGIS 10 and combine it with vector layer

I have connected to a WMS server which provides raster images, more specifically, LIDAR surface elevation maps. I also have a vector layer which contains the outline of buildings. I would like to measure the height of each building above the ground, and add that as an attribute to the vector layer. I have read around but can't find any information on how to work with the WMS layer beyond adding it to the map, and that you can't edit WMS layer. I have two questions:   1. Is it possible to download part/whole WMS layer so it can be edited?   2. Is it possible to extract the numerical value of the height from the raster? I'm using ArcGIS 10.1 Desktop Licence Type Advanced.