title: Spatial Indexes on OS MasterMap Data

text: What are the optimal settings for setting up a SQL Server 2008 Spatial database with full UK OS MasterMap data in?   * Seperate tables for the different feature types, or a single table with everything in it?   * What spatial index settings have you found work best? Any "best practice" or experience gained would be appreciated. For reference, I'm using Cadcorp to display the data.


title: After two buffers operations (+/-), why is the result file bigger than the original?

text: First file ( .shp ) had 1000KB, after 2buffers ( buffers on polygon, 1st buffer +100, 2nd buffer -100 ( small things near polygon should be removed )) in quantum gis file grown to 1500KB, and after the same 2 buffers in postgis file has 20000KB. Anyone can tell me why?


title: How can I merge the attribute tables from feature classes into a single table, easily?

text: I have a large number (approx. 1,100) polygon feature classes. What I would like to do is to merge the attribute tables from each of these feature classes into a single **table** , in the simplest and fastest way possible. This will be the first step in a script or model that then uses that table with all polygon values from all layers to perform some tabular joins and selections. Each time I run this script, I will want to re-do the table merge, as the underlying polygon feature classes may change. If it's possible to write the merged table into memory, for speed, that would be best. I'm using ArcGIS 10.1. I think I can do this easily enough in VBA/VB.net, but it might be simplest to set it up as a ModelBuilder model, so that it's easy to share with colleagues.


title: Simplifying polygons made from binary rasters

text: I'm using MaxEnt to produce a species presence / absence raster layer. It works fine but I have to choose a cut off point for what constitutes presence or absence. As a result I end up with a lot of single cell polygons. Using the "simplify polygons" check box in the raster to polygon tool helps a little but I still have a lot of single cell polygons (that come as triangles left). I've tried using the exclude and dissolve tools but they seem to leave the outlines there. Really I want to change all polygons under a certain area from 0 (absence) to 1 (presence) and then dissolve those boundaries and finally smooth out the boundaries of the polygons I have left. Any help appreciated.


title: Pulling CSV file directly from website into QGIS?

text: I want to be able to automatically update my csv files in a QGIS project. My organization has an online database that can easily export csv files. These csv files are constantly updated with new information. I want to be able to link my QGIS project layers directly to the updated csv files so I do not have to constantly be uploading and rejoining layers to have current maps. Is there an easy way to link the csv file? I am fairly new to QGIS.


title: How to get a list of adjacent lots to a road in ArcMAP and PostGIS/QGIS?

text: In the following map, I have:   * A line / road layer   * A parcel / lot layer The Task:   * Find and count the list of lots adjacent to the road with their types (parcel_type is a field in the parcel layer)   * The data is both stored as shp files and as a postgis tables The Issues:   * The line layer and the lot layer do not intersect i.e. the lot layer has the right of way specified in it of the actual width of the road; whereas, the road layer is just a line layer.    * a parcel might be adjacent to a road on more than one location. think of it as a large lot, with a lot within it. or lots with double frontages. (there is an example of such a lot in the picture below   * This has to be automated on all roads of the study area (roughly 20,000 segments) **Image Sample** ![alt text](http://97.107.136.148/sample_map.jpg)


title: Is there a Raster2pgsql GUI?

text: Has anyone ever attempted to create a raster2pgsql gui (like the shp2pgsql gui). As I am sometimes lazy to use command line, this could be a very nice to have.


title: How can I get ArcGIS Server to give a GML GetFeatureInfo response?

text: Although I have read the following document (http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//005300000228000000) I am still totally unsure about how to enable a gml GetFeatureInfo response from ArcGIS Server 10.0 . Is this something on the client-side that needs to be implemented or something in my GetCapabilities Document?


title: Run ArcGIS 10.1 application on 10.2

text: I support a .NET application that has been compiled on a PC with ArcGIS 10.1 installed. The application comprises a user interface and an ArcGIS desktop extension, plus an installer assembly that is called by Windows Installer to register the desktop extension during installation. The problem I have is that the application won't run on a PC with ArcGIS 10.2 because the 10.1 ESRI assemblies are not found so it crashes. It seems that since v9.3 ESRI no longer supply publisher policy files any more to redirect the program at runtime to the latest versions. I have managed to redirect the program to load the 10.2 assemblies by adding runtime assembly binding into the application's .exe.config. However, this doesn't redirect the desktop extension and won't solve the problem with the installer looking for the 10.1 libraries when installing the application. ESRI's solution seems to be to re-compile the application every time a new version of ArcGIS is released. However, this means I would need to have a copy of every release of ArcGIS that the end users are using. And I'd have to uninstall and re-install different versions every time I need to change the application code in order to re-compile it. Does anyone have any suggestions? Can the libraries be loaded dynamically at runtime to avoid looking for a specific version? Is there any other solution?


title: Volume calculation analysis with Quick Terrain Modeler 8.0.2. (QTM)

text: I have been unsuccessful in finding a tutorial, guide or example of the 'Volume Calculation' analysis tool in Quick Terrain Modeler. http://appliedimagery.com/ I was wondering if anyone else has any success stories they would wish to share and/or how (in)/accurate they have found it to be when comparing (for example) a building while empty versus a building while occupied. How accurate were the measurement(s) between an empty structure and one that is filled, and how long does the calculation generally take (obviously it depends on the size though)? Unfortunately, I'm currently having issues with my LiDAR datasets, so I have yet to test it out myself, nor do I currently have a 'bare' dataset in which to reference to one that is occupied/filled. Is anyone aware of better (or free) software that also does volume calculation analysis? I'm using Quick Terrain Modeler v8.0.2. (x64) Build 80777.


title: Setting up an application in QT Creator that uses QGIS API

text: I am a 15 year veteran of Java making my first entry in to C++ applications. I have a specific need to develop a desktop application that produces graphics type output as well as displaying imbedded maps. It looks like Qt and QGIS API are exactly what I need to get this going. I have been playing around with Qt Creator and generated several of the tutorials without too much problems. Now I am trying my first Qt project that incorporates the QGIS api. I am really stuck. I can't figure out where the libraries are that contains the QGIS modules. I think I have an idea how to configure the Qt project if I could find QGIS. I found all the .h files but some of what I need just aren't there. There are some good examples of code that shows how to code with QGIS classes. But I can't find anything that explains how to set the project up. Could someone please help guide me through setting up a project in Qt Creator that will allow me to incorporate the QGIS modules? I have been searching everywhere. Google can be a strange beast. Sometimes you find what you need on the first search. Other times it just never happens. Thanks!


title: Natural Earth: How to get the position of each country?

text: I've just split the 1:10m Natural Earth countries into individual JSON files which I plan to render and partially stitch back together using RaphaelJS. The 1:10m map is very big (10mb or so) so loading it all in one go in the browser is not really feasible. But how do I get the position of each country? Opening the 1:10m map in ArcMap's Attribute Table shows no signs of any latitude, longitude, X or Y or whatever. Maybe this info needs to be extracted from an exported JSON of the full map using some JS virtual canvas / positioning system? I.e. for each of the individual countries that need to be stitched together, some server-side script would query the full JSON to get their paths? Thanks!


title: Is there a website that will let me see add markers and bounding boxes on top of a street map?

text: Is there a website out there that will let me markers (lat/lon) and bounding boxes (minimum / maximum lat/lon) and let me see it overlayed on top of Google Maps (or OSM, or Bing, etc)? I'm not really looking to create my own webpage for this... I'm hoping someone's implemented something like this for a general use case or a demo and it's online somewhere.


title: How to manipulate OpenStreetMap schema to be useful in ArcGIS applications?

text: I have a shapefile of OSM data. I plan to use the data for an ArcEngine application, and I was wondering if there are any existing solutions to getting the data's schema into your typical GIS format. What I'm talking about is the "tags" attribute. It contains a comma-separated list of categories and their values, whereas typically you would see these values as separate fields in a feature class. **Are there any tools available that separate the tag values out into separate fields?** It would probably be straightforward to create a Python script that separates these out, but I'd like to know if a solution already exists before I create one myself.


title: Improve performances on distance coverage routing - PostgreSQL

text: I have to create layers of land coverage from one point to anywhere around on a given street network for a given distance (ex: FROM house, DISTANCE 30000m) An important factor to understand; I need the result to be exactly superposed over my network, and the coverage to represent only the parts of the road (let say, the addresses) that can be reached from the network for the given distance. Nothing more. I need to be precise, because it will be used for emergency services and efficiency analysis over the covered network. And to add to the fun, the tool must not only give the right answer, but do it in a decent length of time, as it will ultimately be open to other uses on a public web interface (so no, I don't only need an application that will generate the right routing over a night of work!) The point is not to find the shortest way to get to a known point, so basic pgrouting tools can't get me a direct answer. I developped a working algorithm, and it works well on long routes, but as soon a I enter a "real" network (as for a city, or even a village), it crashes... well, it simply run for hours. Basically, my algorithm goes like this: First, we link our network table to our previous results to make sure we don't take a route that as been looked uppon, unless we have more "remaining distance" to cover at this point, which mean that our new route is better. > FOR a in EXECUTE 'SELECT aa.ogc_fid, aa.'||col_contr||' as cost, > aa.the_geom, aa.source, aa.target > >> >>        FROM '||table_name||' aa LEFT JOIN temp res ON res.ogc_fid = aa.ogc_fid AND res.dist > '||distance||' >>   >> >>> >>>         WHERE (aa.source = '||debut||' >>>                 OR aa.target = '||debut||') >>>                 AND (ARRAY[aa.ogc_fid] && ARRAY['||array_to_string(fromSegment,',')||']) = false >>>                 AND res.ogc_fid is null' >>>   >>> >>> LOOP After that (I'll short the code), IF(distance > a.cost) THEN We insert the informations in the result table, and launch back the function from this node, which developps a tree. ELSE (we don't have enough "remaining distance" to cover the line) We return a st_line_substring of the lenght we can cover, and the function stops So it's a reccursive function that tries to cover every possible routes from the startpoint to the given distance... Any idea about how to improve it's speed? I think the answer would be to modify the fist FOR statement... but after hours of brain boiling, I have to admit that I can't do it alone! Thanks **Some new details** ![Beautiful network](http://i.stack.imgur.com/PAMYb.jpg) Suppose this fictive (and horrible, sorry, not on my computer) network, with the blue point as a startpoint for a given distance. ![As wonderful, but with colors!](http://i.stack.imgur.com/ysIkV.jpg) Driving_distance would give a result such as this (red dots being end vertex for the distance, and green ones intermediate vertex on the path). The problem I face here (but still, not on my computer, can't try it tonight, maybe I'm wrong) is that I need a very precise and stuck-on-the-network polyline result. With these dots, I guess I could recreate the path, using some kind of point_to_line function, but it would not work for complexe polylines that are not straight from startpoint to endpoint (like, if the road makes an "S"...). Same problem with isocrone: ![In red dress?](http://i.stack.imgur.com/A6nUF.jpg) The coverage here is not what I need, simply because it covers spaces that cannot be reached on the real network (like the segment in the pink ellipse, which is pass the terminal red dot). (Sorry about the pictures, I'll probably give you some real screenshots of the problem tomorrow)


title: What is the origin of the Degrees, Decimal Minutes format?

text: In my experience, for most non-GIS folks, the most commonly used and understood coordinate system is good old Latitude and Longitude, displayed as degrees, minutes, and seconds (DDMMSS.ss). In GIS, it's of course much easier to use decimal degrees (DD.dd), and with the growth of everyday GIS apps like Google Maps, the non-GIS user is much more familiar with decimal degrees these days than they were 10 years ago. But every now and again, someone brings me a set of coordinates in degrees and decimal minutes. This seems like the worst of both worlds: it can't be used in GIS without a little work, and the format is not immediately obvious to most non-GIS people (many think that it's DDMMSS at first). So why does it exist when we already have decimal degrees, which is easily understood by computers, and degrees-minutes-seconds, which is easily understood by humans?


title: An unhandled exception of type System.AccessViolationException

text: I am trying to upgrade a desktop extension from ArcGIS 10.0 to ArcGIS 10.2. I am getting this error:               An unhandled exception of type 'System.AccessViolationException' occurred in GIS.Metadata.dll.       Not sure what may be causing this. My Environment: VS2013, ArcGIS 10.2, VB.NET, Windows 7 (64 bit). I have also added code where it’s failing. Appreciate any feedback.               pGxCatalog = New ESRI.ArcGIS.Catalog.GxCatalog     '// connect to root ...      pGxCatalog.ConnectFolder(IO.Path.GetPathRoot(pPath))          '// connect to folder via Location property     pGxCatalog.Location = IO.Path.GetDirectoryName(pPath) '       **Fails here**               pCatObject = pGxCatalog.GetObjectFromFullName(pPath, numFound)       IF I set the location manually (say C:) it errors on the this line.


title: Geographic constraints on tables in MapInfo

text: In MapInfo, how can I remove polygons from a table according to whether or not points from another table fall within their boundaries?


title: in seach of gps tracking (web) service to use with rails web app

text: hey out there GISers, i am looking for a gps tracking service i can hook gps devices into and request their infos via a ruby on rails webapp. does anybody know some? any help is appreciated! thanks, martin a site which does similar things is eg.: www.cabulous.com


title: How to find the country in which a point is located?

text: I have a PostGis database created with osm2pgsql. I am looking for a way to add the country name in which a point is located to the point row. It must work also for countries with enclaves, such as France (Llívia). In this case, osm2pgsql creates several polygons from the border relation: a polygon with the outer border and and the enclave, and for each of the both countries a polygon with the enclave only. A simple ST_Within or ST_Contains does not work in this case. I found I way to do the job:               update planet_osm_point p set country=            (SELECT distinct poly.name               from planet_osm_point point inner join planet_osm_polygon poly               on (st_within(point.way, poly.way)                  and (st_nrings(poly.way)=1 or not st_within(point.way, st_makepolygon(ST_InteriorRingN(poly.way,1))))                 and not exists (select * from planet_osm_polygon other                                    where other.admin_level='2'                                          and other.name=poly.name                                          and st_contains(st_makepolygon(st_exteriorring(other.way)), poly.way)                                         and other.way_area>poly.way_area))                 where  poly.boundary='administrative' and poly.admin_level='2' and point.osm_id=p.osm_id)            where p.name is not null;      This seems a bit complex to me and takes a lot of processing time. Is there any simpler and faster way to do this?


title: Incorrectly Projected PostGIS data in TileMill

text: When I bring `PostGIS` data in `EPSG:900913` into `TileMill`, my layers end up on other continents. This is my workflow: (1) Save `shapefiles` in `QGIS` to `EPSG:900913` (2) Import the `shapefiles` using `PostGIS Shapefile and DBF loader 2.0` in `pgadminIII` ... specify SRID as `900913` in the loader (3) Add `PostGIS` layer in `TileMill`... specify SRS as `900913`. I get this result: https://www.dropbox.com/s/ztw2c03o9r7gv30/tileMill_projection.png Some notes: I have no problems loading the `PostGIS` data into `QGIS` afterwards in the same projection (on-the-fly disabled). I have tried the same process using WGS84 `EPSG:4326` & WGS84 / Pseudo Mercator `EPSG:3857` with no luck. My guesses are that either (1) the `Extent` parameters when Adding PostGIS Layer in `TileMill` have something to do with this (2) TileMill is not reading my PostGIS spatial_ref_sys table correctly. Any thoughts and/or suggestions? I'm going to try inserting a new spatial_ref_sys table like here. Thanks again!


title: Nearest Neighbor Focal Statistics and Map Algebra

text: I would like to use a nearest-neighbor type analysis to iterate over a raster. The nearest neighbor analysis will give a cell value of 1 IF AND ONLY IF the cell being analyzed has a cell value larger than all adjacent cells, and 0 otherwise. Is this type of analysis possible with Arcobjects or with the raster calculator tool using map algebra?


title: Creating spatial tables with PostGIS

text: In the PostGIS documentation it says that there are two steps to creating a spatial table with SQL: >   1. Create a normal non-spatial table. >   2. Add a spatial column to the table using the OpenGIS "AddGeometryColumn" > function. > If I followed the examples, I would create a table called `terrain_points` like this:               CREATE TABLE terrain_points (        ogc_fid serial NOT NULL,        elevation double precision,     );          SELECT AddGeometryColumn('terrain_points', 'wkb_geometry', 3725, 'POINT', 3 );      Alternatively, if I look at existing tables in **pgAdmin III** , it seems like I could create the same table like this:               CREATE TABLE terrain_points     (       ogc_fid serial NOT NULL,       wkb_geometry geometry,       elevation double precision,       CONSTRAINT terrain_points_pk PRIMARY KEY (ogc_fid),       CONSTRAINT enforce_dims_wkb_geometry CHECK (st_ndims(wkb_geometry) = 3),       CONSTRAINT enforce_geotype_wkb_geometry CHECK (geometrytype(wkb_geometry) = 'POINT'::text OR wkb_geometry IS NULL),       CONSTRAINT enforce_srid_wkb_geometry CHECK (st_srid(wkb_geometry) = 3725)     )     WITH (       OIDS=FALSE     );     ALTER TABLE terrain_points OWNER TO postgres;          -- Index: terrain_points_geom_idx          -- DROP INDEX terrain_points_geom_idx;          CREATE INDEX terrain_points_geom_idx       ON terrain_points       USING gist       (wkb_geometry);      Do these two methods produce the same result? Is the version based on **pgAdmin III** simply more verbose, and doing things that `AddGeometryColumn` would do by default? Thanks.


title: Error adding buffer to polyline in ArcMap 10

text: As the next step to what I did in this previous question, I would like to add a buffer to a polyline, with 5 meters on either side. I tried adding buffer to two versions of the line, one in geographic WGS 1984, and the other in projected UTM 16N, but in either case I get the error "spatial reference does not match data frame". Does ArcGIS have specific requirements as to whether I need to have the project/data frame in geographical vs. projected coordinate system? And how do I alleviate this problem? Thanks. P.S. I am running ArcMap 10 in Windows.


title: Basic Field Calculator Python Question

text: I'm trying to used the Python syntax in ArcGIS 10's field calculator and am thoroughly confused. So far the most basic calculations are failing me, and I don't understand why. I'm trying to create a function to just set the values of all cells in a column to be a multiple of the argument that I pass in ( I wanted to try something simple to test it out, as anything more complicated I tried also failed). ![enter image description here](http://i.stack.imgur.com/cFiCc.png) However, the error messages are telling me that my syntax is incorrect: ![enter image description here](http://i.stack.imgur.com/HYYA3.png) but I can't understand what is incorrect. * * * As an aside, I'm also horrified at how hard it is to write readable python in the field calculator box as there is no syntax highlighting, or clear indentation. It would be quicker, and easier to write an external Python script to do this calculation, but I want to try to figure out how this new function worked within the program.


title: Google Earth portable solution and imagery data

text: I wish to clarify with users who are using Google Earth portable solution. As per portable solution link http://www.google.com/enterprise/earthmaps/portable.html, we can create part of globe and create offline map to use it. However FAQ section in : http://www.google.com/enterprise/earthmaps/earth_enterprise_faq.html mentioned as below: Does the Google Earth Portable solution allow you to use Google's imagery? No, the Google Earth Enterprise portable solution is intended to be used in situations where your employees need to access your internal imagery and mapping data. Google does not provide customers with Google licensed data for portable use. Google has many partners that provide custom Google Earth Enterprise portable solutions with high resolution imagery. Then i am not sure how to use portable solution in offline without Google's imagery data. Any experiences or details on this solution are welcome. Thanks


title: Is there Open Source equivalent to ArcGIS Integrate (for topology correction)?

text: I need to transform a set of polygons so that it forms an exact spatial partition/tesselation (no gap - no overlap). What are the existing open-source solutions (software or library) to perform this correction? Ideally, I am looking for open-source versions of the ArcGIS integrate. > Integrate is used to maintain the integrity of shared feature boundaries by > making features coincident if they fall within the specified x,y tolerance. > Features that fall within the specified x,y tolerance are considered > identical or coincident.


title: Looking for mobile web app map control that can handle layer redraw with good performance

text: I want to develop mobile web app to track GPS coords in real time (refresh every 10 seconds) from geoserver. Currently I have desktop site I wrote in ExtJS 4 and it works pretty well (Openlayers + GeoJson from WFS) - the layer redraws are not noticable and it gives a feeling that the points are moving (because there are no flickering when layer redraw). I've tried doing the same using sencha touch 2 (found leaflet and openlayers plugins on github) . here what I achieved :   1. Tried google map v3 api - the map renders fast but when drawing the geojson layer or kml layer , in order to update i need to delete the previous features and redraw which cause all the markers to flicker.   2. Openlayers - The map renders very slow when panning / zooming.   3. Leaflet - Also , the map renders very slow (and sometimes idle...) when trying to pan fast or zoom multiple times. so , now I'm stuck with nice UI but poor performance map control (or maybe I don't use the google api correct) This is how I implemented it using google api:               var gMap = Ext.getCmp('map').getMap();     georssLayer.setMap(null);     georssLayer = new google.maps.KmlLayer({        url: 'http://apps.utilis.co.il/geoserver2/wms/reflect?layers=ClearcutCB:DeviceInfo& format=rss&rnd=' + Math.random()});     georssLayer.setMap(gMap);      Will be more than happy having your advise on this one.


title: How to create a many to one spatial join?

text: I have a polygon feature class of territories and a point feature class of trees. The points have many attributes, one of them is the tree species name. I need to get the list of all the tree species present within each polygons. I am trying to perform a many to one spatial join (many points to one polygon), but I am unable to make it run successfully. Ideally, I would have as ouput the same polygon layer, with all the attributes of the trees (including names) attached to it. Any ideas ?


title: Clip Spatial object to bounding box in R

text: Given a Spatial object in R, how do I clip all of its elements to lie within a bounding box? There are two things I'd like to do (ideally I'd know how to do both, but either is an acceptable solution to my current problem--restricting a polygon shapefile to the continental U.S.).   1. **Drop each element not fully within the bounding box.** This seems like `bbox()<-` would be the logical way, but no such method exists.   2. **Do a true clip operation, such that non-infinitesimal elements (e.g. polygons, lines) are cut off at the boundary**. `sp::bbox` lacks an assignment method, so the only way I've come up with would be to use `over` or `gContains`/`gCrosses` in conjunction with a SpatialPolygons object containing a box with the new bounding box's coordinates. Then when clipping a polygon object, you'd have to figure out which are contained vs. cross, and alter the coordinates of those polygons so that they don't exceed the box. Or something like `gIntersection`. But surely there's a simpler way? While I know that there are many problems with bounding boxes, and that a spatial overlay to a polygon that defines the region of interest is typically preferable, in many situations, bounding boxes work fine and are simpler.


title: Tool to search for repeated codes in ArcGIS 10.1

text: I am using ArcGIS 10.1. I have a shapefile that has a field called "Code". I want to make a Modelbuilder model to check if there are "duplicate" codes. Is there a tool that could do this?


title: Does a road connect one grid cell to an adjacent grid cell?

text: So here is the problem I am working on. I have an OpenStreetMap roadmap which I overlaid with a 2-D cartesian grid. So for a given grid square, I want to find out if a road segment continues from the original square to an adjacent square--meaning a square that touches the original square. If I can find out if that road segment continues into an adjacent square, I want to find its direction as well. The goal is to understand whether it is easier for a person residing in one square to get to an adjacent square. So I am saying that if a road exists between squares, then it is easier for a person to move between squares. I am not sure if there is a way to do this in python, GRASS, PostGIS, or other GIS software. Has anyone tried to do this before?


title: Crossover Anaylsis

text: I am new to GIS. I have a text file with (lat,long) of flight lines. If I view the coordinates, they are points, but I am interested in drawing the path of those points, basically connecting them and determining the intersection (crossover) of flight lines (lat, lon). Is there a way to do this in QGIS? Thanks in advance


title: Tool for exporting las data

text: Is there any tool available for extracting the X, Y and Z channels from a .las file and exporting them to a .csv or similar file?


title: PostgreSQL/GIS table relations in QGIS 1.8.0

text: I am building a database in PostgreSQL 9.2 with the PostGIS 2.0 extension. The database consists of a schema called Water which contains water related geographic features (as example, boreholes and water filters as points with coordinates). Each feature type, eg. water filters or boreholes, is in a separate table. In other words, water filters are in a table named water_filters and boreholes are in a table named Boreholes and the tables are in a schema called Water. In the same database I have a schema called maintenance with non-geographic tables. The maintenance tables houses the maintenance records of the respective geographic features. The maintenance table for water_filters will house, for example, a unique number/name/identifier for a water filter and the latest service date. Thus far, I have selected for the geographic features (eg water filters) a unique number/name/identifier as the primary key. For the corresponding maintenance table the number/name/identifier of the geographic features was selected as a foreign key. **My questions** :   1. Is this the correct way to define relations between two tables in PostgrreSQL/GIS?   2. How do I view the relations in QGIS 1.8.0? Thanks in advance


title: Update Attribute on location

text: How do I add a location attribute to a point file, where the location is an ID from an underlying shape file, using ArcMap 10.1? * * * 2 layers, I have polygones with a location ID, 1,2,3,4,5 etc. I have point layer which I would like to add a new column to, and update that column with the ID from the polygon file


title: How to set proportional value in legend?

text: I have a GIS table with traffic flows.   The thickness is proportional to the traffic value and I'd like to have a little legend with 3 links proportional to setted value (eg: link1 proportional to 100 vehicles, link2 proportional to 500 vehicles and link3 proportional to 1000 vehicles).   Any suggestions?   I tried to use the "Legend" but it only represents a coloured link with no thickness.   I am using QGIS 1.8 Lisboa on Win7.   Thanks


title: geoserver configuration 'Could not initialize class org.geotools.image.ImageWorker'

text: I installed geoserver 2.2.4 on windows 7. I then tried to open geoserver demo service with quantum gis desktop 1.8.0. First I tried the wfs service with ghttp://localhos..:8080/geoserver/wfs?version=1.0.0 url and it worked perfectly. Second I tried the same service (tobb:states, tasmania:roads) with ghttp://localho..:8080/geoserver/wms and the service catalog loaded. When I choose one of them I recieve this error. I think I have to add some configuration for wms but I'm not sure.                    Could not initialize class org.geotools.image.ImageWorker           * * * Is it enough?               at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:68)         at org.springframework.security.web.authentication.www.BasicAuthenticationFilter.doFilter(BasicAuthenticationFilter.java:150)         at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:72)         at org.geoserver.security.filter.GeoServerCompositeFilter.doFilter(GeoServerCompositeFilter.java:91)         at org.geoserver.security.filter.GeoServerBasicAuthenticationFilter.doFilter(GeoServerBasicAuthenticationFilter.java:82)         at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323)         at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:68)         at org.springframework.security.web.context.SecurityContextPersistenceFilter.doFilter(SecurityContextPersistenceFilter.java:87)         at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:72)         at org.geoserver.security.filter.GeoServerCompositeFilter.doFilter(GeoServerCompositeFilter.java:91)         at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323)         at org.springframework.security.web.FilterChainProxy.doFilter(FilterChainProxy.java:173)         at org.geoserver.security.GeoServerSecurityFilterChainProxy.doFilter(GeoServerSecurityFilterChainProxy.java:97)         at org.springframework.web.filter.DelegatingFilterProxy.invokeDelegate(DelegatingFilterProxy.java:346)         at org.springframework.web.filter.DelegatingFilterProxy.doFilter(DelegatingFilterProxy.java:259)         at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084)         at org.geoserver.filters.LoggingFilter.doFilter(LoggingFilter.java:71)         at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084)         at org.geoserver.filters.GZIPFilter.doFilter(GZIPFilter.java:41)         at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084)         at org.geoserver.filters.SessionDebugFilter.doFilter(SessionDebugFilter.java:46)         at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084)         at org.mortbay.jetty.servlet.ServletHandler.handle(ServletHandler.java:360)         at org.mortbay.jetty.security.SecurityHandler.handle(SecurityHandler.java:216)         at org.mortbay.jetty.servlet.SessionHandler.handle(SessionHandler.java:181)         at org.mortbay.jetty.handler.ContextHandler.handle(ContextHandler.java:726)         at org.mortbay.jetty.webapp.WebAppContext.handle(WebAppContext.java:405)         at org.mortbay.jetty.handler.ContextHandlerCollection.handle(ContextHandlerCollection.java:206)         at org.mortbay.jetty.handler.HandlerCollection.handle(HandlerCollection.java:114)         at org.mortbay.jetty.handler.HandlerWrapper.handle(HandlerWrapper.java:152)         at org.mortbay.jetty.Server.handle(Server.java:324)         at org.mortbay.jetty.HttpConnection.handleRequest(HttpConnection.java:505)         at org.mortbay.jetty.HttpConnection$RequestHandler.headerComplete(HttpConnection.java:828)         at org.mortbay.jetty.HttpParser.parseNext(HttpParser.java:514)         at org.mortbay.jetty.HttpParser.parseAvailable(HttpParser.java:211)         at org.mortbay.jetty.HttpConnection.handle(HttpConnection.java:380)         at org.mortbay.io.nio.SelectChannelEndPoint.run(SelectChannelEndPoint.java:395)         at org.mortbay.thread.BoundedThreadPool$PoolThread.run(BoundedThreadPool.java:450)     2013-02-13 10:54:34,220 ERROR [geoserver.ows] -      org.geoserver.platform.ServiceException: No service: ( ows )         at org.geoserver.ows.Dispatcher.service(Dispatcher.java:518)         at org.geoserver.ows.Dispatcher.handleRequestInternal(Dispatcher.java:250)         at org.springframework.web.servlet.mvc.AbstractController.handleRequest(AbstractController.java:153)         at org.springframework.web.servlet.mvc.SimpleControllerHandlerAdapter.handle(SimpleControllerHandlerAdapter.java:48)         at org.springframework.web.servlet.DispatcherServlet.doDispatch(DispatcherServlet.java:923)         at org.springframework.web.servlet.DispatcherServlet.doService(DispatcherServlet.java:852)         at org.springframework.web.servlet.FrameworkServlet.processRequest(FrameworkServlet.java:882)         at org.springframework.web.servlet.FrameworkServlet.doGet(FrameworkServlet.java:778)         at javax.servlet.http.HttpServlet.service(HttpServlet.java:707)         at javax.servlet.http.HttpServlet.service(HttpServlet.java:820)         at org.mortbay.jetty.servlet.ServletHolder.handle(ServletHolder.java:487)         at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1093)         at org.geoserver.filters.ThreadLocalsCleanupFilter.doFilter(ThreadLocalsCleanupFilter.java:23)         at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084)         at org.geoserver.filters.SpringDelegatingFilter$Chain.doFilter(SpringDelegatingFilter.java:74)         at org.geoserver.filters.SpringDelegatingFilter.doFilter(SpringDelegatingFilter.java:45)         at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084)         at org.geoserver.platform.AdvancedDispatchFilter.doFilter(AdvancedDispatchFilter.java:49)         at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084)         at org.vfny.geoserver.filters.SetCharacterEncodingFilter.doFilter(SetCharacterEncodingFilter.java:109)         at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084)         at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:311)         at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:68)         at org.springframework.security.web.access.intercept.FilterSecurityInterceptor.invoke(FilterSecurityInterceptor.java:116)         at org.springframework.security.web.access.intercept.FilterSecurityInterceptor.doFilter(FilterSecurityInterceptor.java:83)         at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:72)         at org.geoserver.security.filter.GeoServerCompositeFilter.doFilter(GeoServerCompositeFilter.java:91)         at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323)         at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:68)         at org.springframework.security.web.access.ExceptionTranslationFilter.doFilter(ExceptionTranslationFilter.java:113)         at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:72)         at org.geoserver.security.filter.GeoServerCompositeFilter.doFilter(GeoServerCompositeFilter.java:91)         at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323)         at org.geoserver.security.filter.GeoServerAnonymousAuthenticationFilter.doFilter(GeoServerAnonymousAuthenticationFilter.java:53)         at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323)         at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:68)         at org.springframework.security.web.authentication.www.BasicAuthenticationFilter.doFilter(BasicAuthenticationFilter.java:150)         at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:72)         at org.geoserver.security.filter.GeoServerCompositeFilter.doFilter(GeoServerCompositeFilter.java:91)         at org.geoserver.security.filter.GeoServerBasicAuthenticationFilter.doFilter(GeoServerBasicAuthenticationFilter.java:82)         at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323)         at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:68)         at org.springframework.security.web.context.SecurityContextPersistenceFilter.doFilter(SecurityContextPersistenceFilter.java:87)         at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:72)         at org.geoserver.security.filter.GeoServerCompositeFilter.doFilter(GeoServerCompositeFilter.java:91)         at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323)         at org.springframework.security.web.FilterChainProxy.doFilter(FilterChainProxy.java:173)         at org.geoserver.security.GeoServerSecurityFilterChainProxy.doFilter(GeoServerSecurityFilterChainProxy.java:97)         at org.springframework.web.filter.DelegatingFilterProxy.invokeDelegate(DelegatingFilterProxy.java:346)         at org.springframework.web.filter.DelegatingFilterProxy.doFilter(DelegatingFilterProxy.java:259)         at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084)         at org.geoserver.filters.LoggingFilter.doFilter(LoggingFilter.java:71)         at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084)         at org.geoserver.filters.GZIPFilter.doFilter(GZIPFilter.java:41)         at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084)         at org.geoserver.filters.SessionDebugFilter.doFilter(SessionDebugFilter.java:46)         at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084)         at org.mortbay.jetty.servlet.ServletHandler.handle(ServletHandler.java:360)         at org.mortbay.jetty.security.SecurityHandler.handle(SecurityHandler.java:216)         at org.mortbay.jetty.servlet.SessionHandler.handle(SessionHandler.java:181)         at org.mortbay.jetty.handler.ContextHandler.handle(ContextHandler.java:726)         at org.mortbay.jetty.webapp.WebAppContext.handle(WebAppContext.java:405)         at org.mortbay.jetty.handler.ContextHandlerCollection.handle(ContextHandlerCollection.java:206)         at org.mortbay.jetty.handler.HandlerCollection.handle(HandlerCollection.java:114)         at org.mortbay.jetty.handler.HandlerWrapper.handle(HandlerWrapper.java:152)         at org.mortbay.jetty.Server.handle(Server.java:324)         at org.mortbay.jetty.HttpConnection.handleRequest(HttpConnection.java:505)         at org.mortbay.jetty.HttpConnection$RequestHandler.headerComplete(HttpConnection.java:828)         at org.mortbay.jetty.HttpParser.parseNext(HttpParser.java:514)         at org.mortbay.jetty.HttpParser.parseAvailable(HttpParser.java:211)         at org.mortbay.jetty.HttpConnection.handle(HttpConnection.java:380)         at org.mortbay.io.nio.SelectChannelEndPoint.run(SelectChannelEndPoint.java:395)         at org.mortbay.thread.BoundedThreadPool$PoolThread.run(BoundedThreadPool.java:450)     2013-02-13 11:27:19,877 WARN [config.CustomEditorConfigurer] - Passing PropertyEditor instances into CustomEditorConfigurer is deprecated: use PropertyEditorRegistrars or PropertyEditor class names instead. Offending key [org.geotools.util.Version; offending editor instance: org.geoserver.platform.util.VersionPropertyEditor@18eb4a8     2013-02-13 11:27:20,047 ERROR [geoserver.global] -      ----------------------------------     - GEOSERVER_DATA_DIR: C:\Program Files (x86)\GeoServer 2.2.4\data_dir     ----------------------------------     2013-02-13 11:27:20,113 INFO [org.geoserver] - Loaded style 'burg'     2013-02-13 11:27:20,114 INFO [org.geoserver] - Loaded style 'capitals'     2013-02-13 11:27:20,115 INFO [org.geoserver] - Loaded style 'cite_lakes'     2013-02-13 11:27:20,116 INFO [org.geoserver] - Loaded style 'dem'     2013-02-13 11:27:20,116 INFO [org.geoserver] - Loaded style 'giant_polygon'     2013-02-13 11:27:20,117 INFO [org.geoserver] - Loaded style 'grass'     2013-02-13 11:27:20,118 INFO [org.geoserver] - Loaded style 'green'     2013-02-13 11:27:20,118 INFO [org.geoserver] - Loaded style 'line'     2013-02-13 11:27:20,119 INFO [org.geoserver] - Loaded style 'poi'     2013-02-13 11:27:20,119 INFO [org.geoserver] - Loaded style 'point'     2013-02-13 11:27:20,120 INFO [org.geoserver] - Loaded style 'polygon'     2013-02-13 11:27:20,121 INFO [org.geoserver] - Loaded style 'poly_landmarks'     2013-02-13 11:27:20,121 INFO [org.geoserver] - Loaded style 'pophatch'     2013-02-13 11:27:20,122 INFO [org.geoserver] - Loaded style 'population'     2013-02-13 11:27:20,122 INFO [org.geoserver] - Loaded style 'rain'     2013-02-13 11:27:20,123 INFO [org.geoserver] - Loaded style 'raster'     2013-02-13 11:27:20,124 INFO [org.geoserver] - Loaded style 'restricted'     2013-02-13 11:27:20,124 INFO [org.geoserver] - Loaded style 'simple_roads'     2013-02-13 11:27:20,125 INFO [org.geoserver] - Loaded style 'simple_streams'     2013-02-13 11:27:20,126 INFO [org.geoserver] - Loaded style 'tiger_roads'     2013-02-13 11:27:20,127 INFO [org.geoserver] - Loaded default workspace cite     2013-02-13 11:27:20,130 INFO [org.geoserver] - Loaded workspace 'cite'     2013-02-13 11:27:20,132 INFO [org.geoserver] - Loaded workspace 'it.geosolutions'     2013-02-13 11:27:20,133 INFO [org.geoserver] - Loaded workspace 'nurc'     2013-02-13 11:27:20,134 INFO [org.geoserver] - Loaded workspace 'sde'     2013-02-13 11:27:20,136 INFO [org.geoserver] - Loaded workspace 'sf'     2013-02-13 11:27:20,137 INFO [org.geoserver] - Loaded workspace 'tiger'     2013-02-13 11:27:20,138 INFO [org.geoserver] - Loaded workspace 'topp'     2013-02-13 11:27:20,561 INFO [org.geoserver] - Loaded store 'ch11', enabled     2013-02-13 11:27:20,564 INFO [org.geoserver] - Loaded data store 'ch11'     2013-02-13 11:27:21,248 WARN [referencing.factory] - Axis elements found in a wkt definition, the force longitude first axis order hint might not be respected:     PROJCS["WGS84 / Simple Mercator", GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS_1984", 6378137.0, 298.257223563]], PRIMEM["Greenwich", 0.0], UNIT["degree", 0.017453292519943295]], PROJECTION["Mercator_1SP_Google"], PARAMETER["latitude_of_origin", 0.0], PARAMETER["central_meridian", 0.0], PARAMETER["scale_factor", 1.0], PARAMETER["false_easting", 0.0], PARAMETER["false_northing", 0.0], UNIT["m", 1.0], AXIS["x", EAST], AXIS["y", NORTH], AUTHORITY["EPSG","54004"]]     2013-02-13 11:27:21,302 INFO [org.geoserver] - Loaded feature type 'DNM_DURAK', enabled     2013-02-13 11:27:21,307 INFO [org.geoserver] - Loaded feature type 'ch11'     2013-02-13 11:27:21,317 INFO [org.geoserver] - Loaded layer 'DNM_DURAK'     2013-02-13 11:27:21,486 INFO [org.geoserver] - Loaded feature type 'zonsinir', enabled     2013-02-13 11:27:21,487 INFO [org.geoserver] - Loaded feature type 'ch11'     2013-02-13 11:27:21,488 INFO [org.geoserver] - Loaded layer 'zonsinir'     2013-02-13 11:27:21,495 INFO [org.geoserver] - Loaded store 'arcGridSample', enabled     2013-02-13 11:27:21,497 INFO [org.geoserver] - Loaded coverage store 'arcGridSample'     2013-02-13 11:27:21,517 INFO [org.geoserver] - Loaded coverage 'Arc_Sample', enabled     2013-02-13 11:27:21,523 INFO [org.geoserver] - Loaded coverage 'arcGridSample'     2013-02-13 11:27:21,885 INFO [org.geoserver] - Loaded layer 'Arc_Sample'     2013-02-13 11:27:21,895 INFO [org.geoserver] - Loaded store 'img_sample2', disabled     2013-02-13 11:27:21,895 INFO [org.geoserver] - Loaded coverage store 'img_sample2'     2013-02-13 11:27:21,909 INFO [org.geoserver] - Loaded coverage 'Pk50095', disabled     2013-02-13 11:27:21,910 INFO [org.geoserver] - Loaded coverage 'img_sample2'     2013-02-13 11:27:21,911 INFO [org.geoserver] - Loaded layer 'Pk50095'     2013-02-13 11:27:21,916 INFO [org.geoserver] - Loaded store 'mosaic', enabled     2013-02-13 11:27:21,916 INFO [org.geoserver] - Loaded coverage store 'mosaic'     2013-02-13 11:27:21,919 INFO [org.geoserver] - Loaded coverage 'mosaic', enabled     2013-02-13 11:27:21,920 INFO [org.geoserver] - Loaded coverage 'mosaic'     2013-02-13 11:27:21,949 WARN [imagemosaic.catalog] - BBOXFilterExtractor::extractBasicProperties(): passed typename is null, using: mosaic     2013-02-13 11:27:22,101 WARN [imagemosaic.Utils] - Unable to find sample image for path C:\Program Files (x86)\GeoServer 2.2.4\data_dir\coverages\mosaic_sample\sample_image     2013-02-13 11:27:22,102 WARN [gce.imagemosaic] - Unable to find sample image for path file:/C:/Program%20Files%20(x86)/GeoServer%202.2.4/data_dir/coverages/mosaic_sample/mosaic.shp     2013-02-13 11:27:22,103 INFO [org.geoserver] - Loaded layer 'mosaic'     2013-02-13 11:27:22,182 INFO [org.geoserver] - Loaded store 'worldImageSample', enabled     2013-02-13 11:27:22,182 INFO [org.geoserver] - Loaded coverage store 'worldImageSample'     2013-02-13 11:27:22,185 INFO [org.geoserver] - Loaded coverage 'Img_Sample', enabled     2013-02-13 11:27:22,186 INFO [org.geoserver] - Loaded coverage 'worldImageSample'     2013-02-13 11:27:22,196 INFO [org.geoserver] - Loaded layer 'Img_Sample'     2013-02-13 11:27:22,201 INFO [org.geoserver] - Loaded store 'sf', enabled     2013-02-13 11:27:22,202 INFO [org.geoserver] - Loaded data store 'sf'     2013-02-13 11:27:22,230 INFO [org.geoserver] - Loaded feature type 'archsites', enabled     2013-02-13 11:27:22,230 INFO [org.geoserver] - Loaded feature type 'sf'     2013-02-13 11:27:22,232 INFO [org.geoserver] - Loaded layer 'archsites'     2013-02-13 11:27:22,234 INFO [org.geoserver] - Loaded feature type 'bugsites', enabled     2013-02-13 11:27:22,235 INFO [org.geoserver] - Loaded feature type 'sf'     2013-02-13 11:27:22,236 INFO [org.geoserver] - Loaded layer 'bugsites'     2013-02-13 11:27:22,238 INFO [org.geoserver] - Loaded feature type 'restricted', enabled     2013-02-13 11:27:22,238 INFO [org.geoserver] - Loaded feature type 'sf'     2013-02-13 11:27:22,239 INFO [org.geoserver] - Loaded layer 'restricted'     2013-02-13 11:27:22,241 INFO [org.geoserver] - Loaded feature type 'roads', enabled     2013-02-13 11:27:22,242 INFO [org.geoserver] - Loaded feature type 'sf'     2013-02-13 11:27:22,242 INFO [org.geoserver] - Loaded layer 'roads'     2013-02-13 11:27:22,245 INFO [org.geoserver] - Loaded feature type 'streams', enabled     2013-02-13 11:27:22,245 INFO [org.geoserver] - Loaded feature type 'sf'     2013-02-13 11:27:22,246 INFO [org.geoserver] - Loaded layer 'streams'     2013-02-13 11:27:22,250 INFO [org.geoserver] - Loaded store 'sfdem', enabled     2013-02-13 11:27:22,250 INFO [org.geoserver] - Loaded coverage store 'sfdem'     2013-02-13 11:27:22,254 INFO [org.geoserver] - Loaded coverage 'sfdem', enabled     2013-02-13 11:27:22,254 INFO [org.geoserver] - Loaded coverage 'sfdem'     2013-02-13 11:27:22,358 INFO [org.geoserver] - Loaded layer 'sfdem'     2013-02-13 11:27:22,363 INFO [org.geoserver] - Loaded store 'nyc', enabled     2013-02-13 11:27:22,363 INFO [org.geoserver] - Loaded data store 'nyc'     2013-02-13 11:27:22,377 INFO [org.geoserver] - Loaded feature type 'giant_polygon', enabled     2013-02-13 11:27:22,378 INFO [org.geoserver] - Loaded feature type 'nyc'     2013-02-13 11:27:22,379 INFO [org.geoserver] - Loaded layer 'giant_polygon'     2013-02-13 11:27:22,380 INFO [org.geoserver] - Loaded feature type 'poi', enabled     2013-02-13 11:27:22,381 INFO [org.geoserver] - Loaded feature type 'nyc'     2013-02-13 11:27:22,381 INFO [org.geoserver] - Loaded layer 'poi'     2013-02-13 11:27:22,383 INFO [org.geoserver] - Loaded feature type 'poly_landmarks', enabled     2013-02-13 11:27:22,383 INFO [org.geoserver] - Loaded feature type 'nyc'     2013-02-13 11:27:22,384 INFO [org.geoserver] - Loaded layer 'poly_landmarks'     2013-02-13 11:27:22,386 INFO [org.geoserver] - Loaded feature type 'tiger_roads', enabled     2013-02-13 11:27:22,387 INFO [org.geoserver] - Loaded feature type 'nyc'     2013-02-13 11:27:22,389 INFO [org.geoserver] - Loaded layer 'tiger_roads'     2013-02-13 11:27:22,394 INFO [org.geoserver] - Loaded store 'states_shapefile', enabled     2013-02-13 11:27:22,395 INFO [org.geoserver] - Loaded data store 'states_shapefile'     2013-02-13 11:27:22,401 INFO [org.geoserver] - Loaded feature type 'states', enabled     2013-02-13 11:27:22,401 INFO [org.geoserver] - Loaded feature type 'states_shapefile'     2013-02-13 11:27:22,403 INFO [org.geoserver] - Loaded layer 'states'     2013-02-13 11:27:22,412 INFO [org.geoserver] - Loaded store 'taz_shapes', enabled     2013-02-13 11:27:22,412 INFO [org.geoserver] - Loaded data store 'taz_shapes'     2013-02-13 11:27:22,430 INFO [org.geoserver] - Loaded feature type 'tasmania_cities', enabled     2013-02-13 11:27:22,431 INFO [org.geoserver] - Loaded feature type 'taz_shapes'     2013-02-13 11:27:22,432 INFO [org.geoserver] - Loaded layer 'tasmania_cities'     2013-02-13 11:27:22,433 INFO [org.geoserver] - Loaded feature type 'tasmania_roads', enabled     2013-02-13 11:27:22,434 INFO [org.geoserver] - Loaded feature type 'taz_shapes'     2013-02-13 11:27:22,434 INFO [org.geoserver] - Loaded layer 'tasmania_roads'     2013-02-13 11:27:22,436 INFO [org.geoserver] - Loaded feature type 'tasmania_state_boundaries', enabled     2013-02-13 11:27:22,436 INFO [org.geoserver] - Loaded feature type 'taz_shapes'     2013-02-13 11:27:22,437 INFO [org.geoserver] - Loaded layer 'tasmania_state_boundaries'     2013-02-13 11:27:22,438 INFO [org.geoserver] - Loaded feature type 'tasmania_water_bodies', enabled     2013-02-13 11:27:22,439 INFO [org.geoserver] - Loaded feature type 'taz_shapes'     2013-02-13 11:27:22,440 INFO [org.geoserver] - Loaded layer 'tasmania_water_bodies'     2013-02-13 11:27:22,443 INFO [org.geoserver] - Loaded layer group 'spearfish'     2013-02-13 11:27:22,444 INFO [org.geoserver] - Loaded layer group 'tasmania'     2013-02-13 11:27:22,445 INFO [org.geoserver] - Loaded layer group 'tiger-ny'     2013-02-13 11:27:22,479 INFO [storage.DefaultStorageFinder] - **************************************************************************************************************************************************     2013-02-13 11:27:22,479 INFO [storage.DefaultStorageFinder] - *** Found Java environment variable GEOSERVER_DATA_DIR set to C:\Program Files (x86)\GeoServer 2.2.4\data_dir, using it as the default prefix. ***     2013-02-13 11:27:22,480 INFO [storage.DefaultStorageFinder] - **************************************************************************************************************************************************


title: Problem with deployment of GeoProcessors

text: I am developing an AddIn which requires lots of Geoprocessing but whenever I deploy any Geoprocessor, the ArcMap crashes and its always with the Geoprocessors only. All other commands are working fine But whenever I try to run Geoprocessing, ArcMap crashes; whether its Geoprocessing related to Densify, Interpolate Shape, Split etc. Please Help Thanks in Advance


title: What is the best ETL tool within ArcGIS standard to bring CSV to an existing schema with domains and other constraints?

text: We will be bringing some GPS related data from a business partner to our existing schema. We have provided them the sample data and data model from ArcGIS Diagrammer. There are bunch of data import options available in ArcMap and ArcCatalog. I am looking for an efficient approach to finish the ETL process with proper messaging about the exceptions of domains, datatype and mandatory/non-mandatory - so that any novice user of ArcMap can follow it. I don't want any third party tool or code development, but I can build a model, if there are 2-3 tools involved. I am new to data interoperability - has anyone used it in similar situation. -Thanks


title: Most optimal algorithm for contour correctness?

text: Am creating an app which correctly calculates contours from a DEM file or any height-map given to it. The correctness of the contour holding the most importance, since a person viewing the contours shouldn't end up taking wrong planning decisions and maybe losing their life by the route they choose after looking at the map. I have already tried using GDAL 1.7 source code, compiled it and added some features of my own to it for some other objective. Now I want to add-to/modify the source to be able to use the contour generation function of GDAL in C++. **But before proceeding, I wanted to know:**   **1.** If GDAL offers the best free contour generation algorithm available or is there any other library which I can link to my code and feed my DEM map to it and receive contour vector lines (or go thru research papers and implement the contour generation algorithm myself in C++)? If GDAL, then why is it the best?   **2.** A generated contour will basically be points connected with straight lines. Is it more correct for a contour to be depicted with connected straight lines or will it give a more correct result to use a spline to smoothen out the lines? I won't be interpolating any contours or the DEM because I've read on this website that interpolation leads to incorrect results.


title: Merging OSM files with OSMOSIS and import to Postgres with OSM2PGSQL

text: I have a problem merging and then importing two OSM files. I downloaded to distinct .osm-files from cloudmade. They are not sharing any geographical attribute (as the first is the country Monaco, the 2nd is the county Bremen in Germany). I merged the with OSMOSIS and it gives me an output file. Then I import it via OSM2PGSQL. Looking at the dataset using QGIS it either shows me the data of Bremen or Monaco (depending the order in the OSMOSIS merge statement). Does anyone has a solution for this problem? Thanks a lot Carsten


title: How to change Split Vector by Attribute to produce files with attribute values as file names?

text: Using: QGIS 2.0.1 I am trying to separate a large fauna point .shp file into many shapefiles based on each species' unique name attribute (i.e. by "SCI_NAME" field)??? I am currently using "Split vector layer by attribute" from the toolbox but am unable to name the resulting .shp files with the attribute name (It is currently naming them with say: "SCI_NAME_1", "SCI_NAME_2" ...) The correct naming issue is the main concern. Unless a different methodology is suggested. I have seen this post: > The SAGA algorithm has been removed and there is a native one implemented as > an example script in SEXTANTE. It is called "Split vector layer by > attribute". > > answered Jan 28 at 17:32 Victor Olaya > > The resulting files are named, appending a number. You you edit the script > and edit the following line: outputFile = output + "_" + str(len(writers)) + > ".shp" to suit your needs in case you want a differnt naming – Victor Olaya > Jan 28 at 17:35 ... but have been unable to work out how to "edit the script". (hehe: have not ventured into the world of scripting/programming {or whatever all that text means} yet) ... but will be able to work out copy/paste if you transcribe(?) how it should look. Cheers. **EDIT:** i have been informed how to change a line in the script to this: > outputFile = output + "_" + str(clazz) + ".shp" which get's the attribute value into the filename, however, it leaves the field/column title in the front of the filename which is not good for organising/searching my .shp files. when i simply remove this part (in bold): outputFile = **output + "_" +** str(clazz) + ".shp" the script no longer runs (or no longer produces output files). can anyone help me overcome this problem? (i.e. how to hide the field/column heading presence in the filename)


title: Attribute manipulation on shapefile

text: I have a shapefile which contains landuse types. The Landuse types are written in a field called "Type" and contain also the area in m² and a landuse code. Some polygons contain more than one landuse type separated with a tab An example or such a field would be (43001-1200) Brachland 1185 m² (43002-1300) Laub- und Nadelholz 14 m² What I would like to do is make a list of the landuse types that interest me, and then search the field for these types. If the word is contained within the field, that polygon should be selected. I then need to remove the code and the area leaving just the Landuse type. IE int he above example the resulting polygon would simply have Brachland, Laub- und Nadelholz in the field. I have ArcGIS (Editor License) and FME Desktop. Does anyone have an idea which approach would work best? Thanks for any replies, Rob


title: Apply styles to WMS layer

text: I am working with OpenLayers (and GeoServer) and I am showing a WMS layer with countries in a map. Each object (country) in the layer has an attribute called NAME (the name of the country). I created a filter like the following in order to see only Italy (it's an example)               myFilter = new OpenLayers.Filter.Comparison({         type: OpenLayers.Filter.Comparison.EQUAL_TO,         matchCase: false,          property: 'NAME',         value: 'Italy'     });      This work perfectly and if I apply it to the layer and call redraw() method I can see only one country. Now I would like to do the same thing with sytles, in order to highlight Italy but keep visible all the other countries in the map. I created a style like the one below, but now I don't know how I can apply this to the layer.               var style = new OpenLayers.Style();     var ruleLow = new OpenLayers.Rule({         filter: myFilter,         symbolizer: {fillColor: "green"}     });     style.addRules([ruleLow]);                Can I apply the style to a WMS layer? Or only to vector layer? Thank you


title: How to change value of an Line edit dynamically

text: I create a form in QT4 and compile with pyuic4. it create an form.py I want to modify the value of the line edit dynamiclly with code. I put this code in the form.py but never happen:               self.pointStartId.setText(QtGui.QApplication.translate("Dialog",     "example", None, QtGui.QApplication.UnicodeUTF8))      how can I do that?


title: In MapBasic, how do I display a dialog box as the result of a menu item being clicked on?

text: I have created a menu item in MapBasic. I need to call a procedure to open a dialog box as the result of one of the menu items being selected. I believe I am having trouble declaring the sub procedure that will be called. I have not included the sub procedure, because it is not right. Here is my code:               Declare Sub Main     Sub Main         Create Menu "Test" ID "12" As             "Preferences",             "Leads",             "LeadsRemain",             "Accts",             "AcctsRemain"         Alter Menu Bar Add "Test"     End Sub


title: Adding feature layer to ArcMap from url

text: Is it possible to somehow load feature layer from REST (like: ArcGIS Example Service) to ArcMap? No need to edit that layer, we just need to select data from it and then create a local FeatureClass.


title: Add layer to open mxd using ArcObjects from outside the map document crashes ArcMap

text: I have an Access database (MS Access 2010) that allows users to enter a variety of land treatment information including associated points, lines, or polygons which users can select and the database would then load as part of the project, to a master geodatabase, and allow the user to view the features in an ArcMap. In Arc 9.3 this was all accomplished with VBA and worked well. Because of the upgrade to Arc 10, VBA is no longer an option and I have been learning how to program in python and call the commands through Access using the Shell feature. I have nearly completed everything except for the ability to allow users to add their feature to the current ArcMap document. Here are the steps that should occur. 1) User clicks a button in the Access program which identifes the feature(s) of interest. 2) Access determines if the correct mxd is open and if not opens it (the document is usually open so this step is often uneccessary) 3) The Python script opens using the Python shell command, identifes the correct and open mxd, and adds feature to the open ArcMap document I have tried both arcpy and ArcObjects. I can run arcpy within the open map but opening a Python file via VBA in MS Access does not manipulate the open mxd. The ESRI forum suggested that I work with ArcObjects. I have found and used code to identify and manipulate layer already loaded in the ArcMap document of interest but am completely stuck on how to identify and add a new layer. Here is the code so far, I think it will work if I can correctly identify the layer file in the last section of code. At least it all runs without the pMap.AddLayer(FC) near the bottom of the code, but I am an ecologist and not a programmer by training so there may be something else I am missing. Any help would be appreciated. Update: I am now very close. I can identify the feature in the geodatabase but now ArcMap crashes when attempting to add the data to the map and I generate the error: COMError: (-2147467259, 'Unspecified error', (None, None, None, 0, None)). However, my output appear to identify the correct feature. I am unsure as to why ArcMap is crashing. Please help, I am so close.               import arcpy     from arcpy import env     import os     import sys     import comtypes.gen.esriArcMapUI as esriArcMapUI     import comtypes.gen.esriCarto as esriCarto     import comtypes.gen.esriGeoDatabase as esriGeoDatabase     import comtypes.gen.esriFramework as esriFramework     import comtypes.gen.esriSystem as esriSystem     import comtypes.gen.esriDataSourcesFile as DataSourcesFile          #Eventaully these will be the arguments that call the feature     #test = sys.argv[0]     #Name = sys.argv[1]     #Alias = sys.argv[2]          def NewObj(MyClass, MyInterface):         from comtypes.client import CreateObject         try:             ptr = CreateObject(MyClass, interface=MyInterface)             return ptr         except:             return None          def CType(obj, interface):         try:             newobj = obj.QueryInterface(interface)             return newobj         except:             return None          def GetApp():         """Get a hook into the current session of ArcMap"""          #from comtypes.gen import esriFramework         pAppROT = NewObj(esriFramework.AppROT, esriFramework.IAppROT)         if pAppROT is not None:             iCount = pAppROT.Count              if iCount == 0:                  print 'No ArcGIS application currently running.  Terminating ...'                  return None              for i in range(iCount):                  pApp = pAppROT.Item(i)  #returns IApplication on AppRef                 print pApp.Name                 if pApp.Name == 'ArcMap':                      print "ArcMap found"                     pDoc = pApp.Document                     print pDoc.Title                     if pDoc.Title == "CreateShapefile10x.mxd" or pDoc.Title == "CreateShapefile10x":                         return pApp          print 'No ArcMap session is running at this time.'      print "No AppROT found"     print "Failed"     return None          pApp = GetApp()          try:         pDoc = pApp.Document         pMxDoc = CType(pDoc, esriArcMapUI.IMxDocument)         pMap = pMxDoc.FocusMap         if pMap.Name == "LTDL Layers":                  ###New Code###                  sPath = r"Y:\LTDL_Data\LTDL_Files\Abate_Seeding_1965\Abate_Seeding_1965.gdb"             sPath1 = r"Abate_Seeding_1965"             pWSF = NewObj(esriDataSourcesGDB.FileGDBWorkspaceFactory, esriGeoDatabase.IWorkspaceFactory)             print pWSF             pWS = pWSF.OpenFromFile(sPath, 0)             pDS = CType(pWS, esriGeoDatabase.IDataset)                  print "Workspace name: " + pDS.BrowseName             print "Workspace category: " + pDS.Category                  FWSS = CType(pWS, esriGeoDatabase.IFeatureWorkspace)             FL = NewObj(esriCarto.FeatureLayer, esriCarto.IFeatureLayer)             FL.FeatureClass = FWSS.OpenFeatureClass(sPath1)             print FL             print FL.FeatureClass.AliasName             LAN = CType(FL, esriCarto.ILayer)             LAN.Name = FL.FeatureClass.AliasName + ": Project Boundary"             print LAN.Name                  ###New Code###                  pMap.AddLayer(FC)             pLayout = pMxDoc.PageLayout             pActiveView = CType(pLayout, esriCarto.IActiveView)             pActiveView.Refresh()             pMxDoc.UpdateContents         print "Done"     except:         print "Failed"


title: set anchor for a PictureMarkerSymbol in arcgis android sdk

text: I am trying to create a `PictureMarkerSymbol`in arcgis android, however I found that the marker use the center of the drawable as the anchor point, is it possible to set the anchor? * * * update: This is the drawable(with 40x56 in pixel): ![enter image description here](http://i.stack.imgur.com/4H0wU.png) This is my code for adding marker:               protected void onCreate(Bundle savedInstanceState) {         final MapView mapView = (MapView) findViewById(R.id.map_imp_esri);         mapView.addLayer(new ArcGISTiledMapServiceLayer("http://services.arcgisonline.com/ArcGIS/rest/services/World_Topo_Map/MapServer"));         Point p = new Point(-9664114.480484284,3962469.970217699); //Birmingham         mapView.setResolution(mapView.getMinResolution());              GraphicsLayer graphicsLayer = new GraphicsLayer(mapView.getSpatialReference(), new Envelope(-180, -90, 180, 90));         mapView.addLayer(graphicsLayer);              //add marker         Drawable d = getResources().getDrawable(R.drawable.marker_a_normal);         PictureMarkerSymbol sym = new PictureMarkerSymbol(d);         sym.setOffsetX(-convertToDp(d.getIntrinsicWidth()) / 2);         sym.setOffsetY(-convertToDp(d.getIntrinsicHeight()) / 2);         Graphic g = new Graphic(p, sym);         graphicsLayer.addGraphic(g);     }                    private int convertToDp(double input) {         final float scale = getResources().getDisplayMetrics().density;         return (int) (input * scale + 0.5f);     }      But this is what I got:![enter image description here](http://i.stack.imgur.com/4V9TS.png): The four screenshot are captured during zooming in the map. Note, at the image 1, the marker event does not at the united states. Do I miss anything?


Is there a way to build topology and work with topology relations in PyQGIS? I'm asking because I need to run some custom generalization algorithms on my data, and I need to preserve topology. Presumably this would be easier if there is already some functionality for handling topology in QGIS. The data is basically shapefiles with multipart-polygons.


title: Combining polygon geometries with python and QGIS segfaults

text: I am trying to combine polygons, using Qgis 1.8, and every time I run this script in the python console, I get a segmentation fault. Here is the script:               from PyQt4.QtCore import *          iface = qgis.utils.iface;     legend = iface.legendInterface()     poly_layer = iface.activeLayer()          MapLayer = QgsMapLayerRegistry.instance()          #######NEW POLY LAYER TO CREATE#########     poly_provider = poly_layer.dataProvider()          fields = { 0 : QgsField("STAT_NO", QVariant.String) }          writer = QgsVectorFileWriter("/localstore/data/shps/contrib_area.shp",     "CP1250", fields, poly_provider.geometryType(), poly_layer.crs(), "ESRI Shapefile")          if writer.hasError() != QgsVectorFileWriter.NoError:       print "Error when creating shapefile: ", writer.hasError()                    #######POINT PROCESSING########     #get the point layer that contains the points contained by the polygons     point_layer = None          for layer_id, layer in MapLayer.mapLayers().iteritems():       if(layer.name() == "slope_and_catchment"):         point_layer = MapLayer.mapLayer(layer_id)         break          point_provider = point_layer.dataProvider()     point_id_index = point_provider.fieldNameIndex("SITE_NO")     point_link_index = point_provider.fieldNameIndex("pw_comb")     point_provider.select([point_id_index,point_link_index])     point_feat = QgsFeature()          intersects = []          point_feat_dict = {}     while point_provider.nextFeature(point_feat):       point_attrmap = point_feat.attributeMap()       point = point_attrmap[point_link_index].toString()       intersects.append(point)       point_feat_dict[point] = point_attrmap[point_id_index].toString()          for point in intersects:       print point               ########LINE PROCESSING######     #get the line layer that contains the stream network     line_layer = None          for layer_id, layer in MapLayer.mapLayers().iteritems():       if(layer.name() == "projectednet"):         line_layer = MapLayer.mapLayer(layer_id)         break          line_provider = line_layer.dataProvider()     line_index = line_provider.fieldNameIndex("LINKNO")     line_US1_index = line_provider.fieldNameIndex("USLINKNO1")     line_US2_index = line_provider.fieldNameIndex("USLINKNO2")          line_provider.select([line_index,line_US1_index,line_US2_index])     line_feat = QgsFeature()               line_feat_dict = {}     #loop over all features and store relevant information     while line_provider.nextFeature(line_feat):       line_feat_id = line_feat.id()       line_attrmap = line_feat.attributeMap()       link = line_attrmap[line_index].toString()       us1 = line_attrmap[line_US1_index].toString()       us2 = line_attrmap[line_US2_index].toString()       line_feat_dict[link] = [line_feat_id, link, us1, us2]           #all_lines = line_feat_dict.keys()          #line_feat_dict is now a map from link names to feature ids, and the ustream links               ########POLYGON PROCESSING#######     #set up attribute information     poly_id_index = poly_provider.fieldNameIndex("DN")     #get all features allong with the attributes     poly_provider.select([poly_id_index])     #make a feature variable to store features in     poly_feat = QgsFeature()          #######NOW WE NEED TO GET ALL CATCHMENTS THAT SHARE LINKS OUTLINED IN STREAM NETWORK########     poly_feat_dict = {}     #get a list of all related polygons:     while poly_provider.nextFeature(poly_feat):       poly_id = poly_feat.id()       poly_attrmap = poly_feat.attributeMap()       poly = poly_attrmap[poly_id_index].toString()       poly_feat_dict[poly] = poly_id          #poly_feat_dict now maps the link number to the feature id          #start traversing the network and getting geometries     geoms = []          def traverse(link1, link2):       #work down one side of the "tree"       if link1 != "-1":         feat = QgsFeature()         poly_provider.featureAtId( poly_feat_dict[link1], feat )         geoms.append( feat.geometry() )         #traverse this links uplinks         traverse( line_feat_dict[link1][2], line_feat_dict[link1][3] )       else:         #have hit headwater         return       if link2 != "-1":         feat = QgsFeature()         poly_provider.featureAtId( poly_feat_dict[link2], feat )         geoms.append( feat.geometry() )         #traverse this links uplinks         traverse( line_feat_dict[link2][2], line_feat_dict[link2][3] )       else:         #have hit headwater         return               #process the catchment up to station       for point in intersects:       #make sure the geoms list is clear for each point       geoms = []       #point is the starting link number       #get this initial geometry from poly_feat_dict       feat = QgsFeature()       poly_provider.featureAtId( poly_feat_dict[point], feat )       geoms.append( feat.geometry() )            #get the corresponding upstream links       up1 = line_feat_dict[point][2]       up2 = line_feat_dict[point][3]       #traverse the stream network links till we can get all geometries       traverse(up2, up2)       #have all geometries, join them       union_geoms = geoms[0]       for geometry in geoms:         union_geoms = union_geoms.combine(geometry)            fet = QgsFeature()       fet.setGeometry(union_geoms)       fet.addAttribute(0, QVariant( point_feat_dict[point] ))       writer.addFeature(fet)               del writer          layer = QgsVectorLayer("/localstore/data/shps/contrib_area.shp", "contrib_area", "ogr")          MapLayer.addMapLayer(layer)          if iface.setActiveLayer(layer):       print "New Layer is Active"     else:       print "Couldn't set active layer"      The following is the GDB backtrace from the coredump:               Thread 1 (Thread 0x7ffff2f27820 (LWP 29411)):     #0  0x00007ffff3645cb0 in geos::geom::Geometry::getPrecisionModel() const () from /usr/lib64/libgeos-3.3.2.so     No symbol table info available.     #1  0x00007ffff369f876 in   geos::operation::GeometryGraphOperation::GeometryGraphOperation(geos::geom::Geometry const*, geos::geom::Geometry const*) () from /usr/lib64/libgeos-3.3.2.so     No symbol table info available.     #2  0x00007ffff36baea0 in geos::operation::overlay::OverlayOp::OverlayOp(geos::geom::Geometry const*, geos::geom::Geometry const*) () from /usr/lib64/libgeos-3.3.2.so     No symbol table info available.     #3  0x00007ffff36bb901 in geos::operation::overlay::OverlayOp::overlayOp(geos::geom::Geometry const*, geos::geom::Geometry const*, geos::operation::overlay::OverlayOp::OpCode) () from /usr/lib64/libgeos-3.3.2.so     No symbol table info available.     #4  0x00007ffff46165da in std::auto_ptr<geos::geom::Geometry> geos::geom::BinaryOp<geos::operation::overlay::overlayOp>(geos::geom::Geometry const*, geos::geom::Geometry const*, geos::operation::overlay::overlayOp) () from /usr/lib64/libgeos_c.so.1     No symbol table info available.     #5  0x00007ffff46109ae in GEOSUnion_r () from /usr/lib64/libgeos_c.so.1     No symbol table info available.     #6  0x00007ffff56b1294 in QgsGeometry::combine(QgsGeometry*) () from /opt/lib64/libqgis_core.so.1.8.0     No symbol table info available.     #7  0x00007fffe346185e in meth_QgsGeometry_combine () from /opt/share/qgis/python/qgis/core.so     No symbol table info available.     #8  0x00000034c42d59e4 in PyEval_EvalFrameEx () from /usr/lib64/libpython2.6.so.1.0     No symbol table info available.     #9  0x00000034c42d7657 in PyEval_EvalCodeEx () from /usr/lib64/libpython2.6.so.1.0     No symbol table info available.     #10 0x00000034c42d7732 in PyEval_EvalCode () from /usr/lib64/libpython2.6.so.1.0     No symbol table info available.     #11 0x00000034c42d5c92 in PyEval_EvalFrameEx () from /usr/lib64/libpython2.6.so.1.0     No symbol table info available.     #12 0x00000034c42d7657 in PyEval_EvalCodeEx () from /usr/lib64/libpython2.6.so.1.0     No symbol table info available.     #13 0x00000034c42d5aa4 in PyEval_EvalFrameEx () from /usr/lib64/libpython2.6.so.1.0     No symbol table info available.     #14 0x00000034c42d7657 in PyEval_EvalCodeEx () from /usr/lib64/libpython2.6.so.1.0     No symbol table info available.     #15 0x00000034c42d5aa4 in PyEval_EvalFrameEx () from /usr/lib64/libpython2.6.so.1.0     No symbol table info available.     #16 0x00000034c42d6b8f in PyEval_EvalFrameEx () from /usr/lib64/libpython2.6.so.1.0     No symbol table info available.     #17 0x00000034c42d7657 in PyEval_EvalCodeEx () from /usr/lib64/libpython2.6.so.1.0     No symbol table info available.     #18 0x00000034c426acb0 in ?? () from /usr/lib64/libpython2.6.so.1.0     No symbol table info available.     #19 0x00000034c4243c63 in PyObject_Call () from /usr/lib64/libpython2.6.so.1.0     No symbol table info available.     #20 0x00000034c42566af in ?? () from /usr/lib64/libpython2.6.so.1.0     No symbol table info available.     #21 0x00000034c4243c63 in PyObject_Call () from /usr/lib64/libpython2.6.so.1.0     No symbol table info available.     #22 0x00000034c42cfc93 in PyEval_CallObjectWithKeywords () from /usr/lib64/libpython2.6.so.1.0     No symbol table info available.     #23 0x00007fffe45bf3db in ?? () from /usr/lib64/python2.6/site-packages/sip.so     No symbol table info available.     #24 0x00007fffe39927c3 in ?? () from /usr/lib64/python2.6/site-packages/PyQt4/QtGui.so     No symbol table info available.      My question is am I doing something wrong, or is this a bug? Any thoughts would be appreciated, Thanks! EDIT: More stack trace is available if needed. I only included the last 25 steps.


title: Join Tables related by a composite primary key in qgis

text: I have two tables, one has the geometry and other has attributes. The table with geometry has unique number primary key, and 4 columns `admin1`,`admin2`,`admin3`,`featureid`. Table two is a csv file has the same columns mentioned above and their combination forms a composite primary key. How can i join the two tables?


title: Displaying extruded polygon outlines in ArcScene 10.2.1

text: I am having a similar problem to the question posted on this site: How to display Arcscene polygon extrusion edges? When I extrude polygons in Arcscene 10.2.1, the polygons rise to the specified height, but the outlines of each feature are not displayed (Mapping Problem- top image). In particular, I'd like to see the outlines for each country, regardless of their extrusion height. I've tried the suggestions listed in the aforementioned post, but do not get the desired effect. I even tried the user's solution to use the Feature to 3D by Attribute tool on the polygon layer. I then added that output layer to Arcscene without any extrusion, symbolized with no fill and a black outline, but again fail (Mapping Problem- bottom image). I'm very new to ArcScene (but well-versed in ArcMap), so any help would be appreciated.


title: Retrieval of only new values from the database

text: I am using Postgis as my backend for my project along with geoserver and openlayers. I have a table in my DB which gets filled up very frequently. Now when a user gets connected to a server , he can see the map with all the markers whose points are in the DB. Imagine the table keeps on filling up and there is no relation between user viewing the map and tables getting filled by values. Now suppose after some time when user clicks on certain button on the page, he has to see the map with markers of all the values present in the DB(along with new points inserted in the table since last time he clicked that button). This task can be made simple by making that button as submit button, but when the user clicks it, it reads the entire DB again(reloads the page) along with old points and projects on the map, instead of reading only the new points. This method is not at all good when I have 10000 points already in DB and only say 5 points getting added in the meantime. In this case it reads all the 10005 points instead of only 5 points. So can anybody give me a solution on how to retrieve only the new points, still projecting all the points present in DB. To my knowledge Ajax can be used to do the same(which doesnt need reloading of the page) but I want a clear explanation of how to carry this work using postgis and geoserver. Are there any special functions which can handle this?? If not, what is the better way to do this? Thanking you in advance.


title: Average frequency of zip code change for a given address over the course of year

text: This is a difficult question to answer with a high degree of accuracy, but could someone estimate to within an order of magnitude for me the chances of a zip code for a given address changing over the course of a year? For example, would you say the chances of this are: a) ~10% b) ~1% c) ~.1% d) ~.01% Just trying to get a rough sense of the constancy of these things for some code I'm writing. If the answer were ~10% I would proceed very differently than if it were ~.01%.


title: Making multiple legend boxes all same size using ArcGIS for Desktop?

text: How do I make all the rectangle or (legend box) all the same ? Look at the circles I marked in the pictures. It should be very simple thing to do. Anyone out there have a good workaround ? ![Example ](http://i.stack.imgur.com/St7lB.png)


title: Can you use GetLegendGraphic for a ArcGIS Server 9.3 WMS?

text: The following request should return a legend graphic for a particular layer: http://myserver/arcgis/services/default/MapServer/WMSServer?Service=WMS&Request=GetLegendGraphic&Version=1.1.1&Layer=3&FORMAT=image/png However instead I get the following error:               <?xml version="1.0" encoding="UTF-8" standalone="yes" ?>     <!DOCTYPE ServiceExceptionReport SYSTEM "http://schemas.opengis.net/wms/1.1.1/exception_1_1_1.dtd">     <ServiceExceptionReport version="1.1.1">       <ServiceException code="InvalidFormat">     Can't parse XML request.       </ServiceException>     </ServiceExceptionReport>      The GetCapabilities request for the service for the layer contains the following:               <Style>     <Name>default</Name>     <Title>Layer Name</Title>     <LegendURL width="100" height="62">     <Format>image/png</Format>     <OnlineResource xlink:href="http://server/arcgisoutput/wms/default1.png" xlink:type="simple"/>     </LegendURL>     </Style>      The OnlineResource href is a valid PNG file that is accessible externally. It seems others have had issues, and recommend using external capabilities files, but when the legend already has a valid image file is this necessary? Has anyone successfully configured a WMS to support the GetLegendGraphic request?


title: How to view and change the linear units of a feature class in ArcGIS 10.0

text: According to the documentation for Near Analysis the output units for the tool are in the linear units of the feature. When accessing the properties of the feature class I can see where angular units are defined (as a result of the coordinate system), however I do not see where I can view or change the linear units of a feature. The Data Source textbox from the Source tab for the Layer Properties of my point feature class does **not** contain a Linear Unit specification. How do I access the linear units of a feature class?


title: Locating and using the FMV toolbox

text: I've just installed version 1.2 of the Full Motion Video Add-In, along with the new geoprocessing toolbox. After running the installation for the tools, it creates a folder on the C drive but does not contain a toolbox. There is no link from the FMV Video Manager pane to the toolbox, nor does it appear in the system toolboxes. I tried to "Add Toolbox" and looked in all the normal places (Documents, Users, C drive etc) and came up empty. Where is the toolbox located? I would like to test out the new tools as mentioned in the blog post. I'm running ArcGIS 10.2.1.


title: analysis of srtm data for gaps

text: I need to analyse the SRTM DATA (1025 X 1025) of a particular terrain which has gaps or voids. Could anybody suggest a solution to find the percentage of gaps or voids in the SRTM data? I would like to know if there is a tool for such data analysis? It is also possible for me to share the data if needed. The file format is (.geotiff)


title: Finding X, Y of nodes in feature layer

text: I have a map document in ArcMap and I want to find X and Y of my nodes in a special feature layer. I look through my access database and I just find the Shape Index table which has MinGX, MinGY, MaxGX, MaxGY. I will be thankful if you help me how to find the X , Y of my nodes. Generally, from what source is ArcMap locating the nodes on the map?


title: How do I create a map that is not to scale in QGIS?

text: I want to make a map that isn't to scale and I'm so used to being pedantic about CRS and scale that I don't know where to start. The map that I need to make is a simple guide for people to get to a site that is out of town. I'm addicted to QGIS, of course.


title: problem with open popup onclick marker in openlayers?

text: I am very new to openlayers. I have done showing map and whenever we click on the map it is showing marker. But i need to open the popup when i click on the marker. The following code has some problem. When i click on marker popup is displaying but at the corner it is displaying. I want to display the popup on the marker it self


title: What is the best way to prepare digital orthophotos for the web?

text: My system: Linux Ubuntu 10.04 server (100gb: 4gb RAM: 32bit) /tomcat6/geoserver/geowebcache/FWTools My Objective: Create Aerial Photo coverage in ETRS89 N32 My Data: 1300 3 band Orthophotos in Tif Formal with a resolution of 12,5m @ 75MB Question: How should I do this and at the same time save space on the server? I am thinking of reducing 3 bands to 1 , reduce resolution by half then image- mosaic in geoserver and serve with geowebcache. Has anyone every done this? I would be pleased to hear from anyone who has done this successfully. yours, Rob


title: Is there an XML equivalent to GeoJSON?

text: In my web application, I am retrieving geospatial data from a server backend (Java) to the browser frontend (JavaScript). I am considering two options of how to serialize the data in this communication - `XML` and `JSON`. `XML` seems to be a faster option as it can be deserialized by browser itself (into a documentElement), while `JSON` must be deserialized by a JavaScript library. On the other hand, there is a standard for storing geospatial data in `JSON` (`GeoJSON`), so using `JSON` would make the application more standard- compliant. So my question is, is there an `XML` equivalent to `GeoJSON`? I found out that the term `GeoXML` has already been coined, but it is the name of a library, not a format standard.


title: hyperlink in a package

text: I want to give a map that contains hyperlinks as a package of adobe or an exe file to someone who does not have access to arcmap. this package should have just some points and polygons that when he click on each of them, a file which is linked to that point or polygon opens. how can i do that? is there any tutorial file about how making this simple package? could you help me please? yours faithfully hamid azari


title: ArcGIS 10: How to clip point buffer to features in another layer?

text: I'm on a steep learnig curve here, so aplogies if my question is a one-banana problem. In ArcGIS 10 I have a contour layer that was generated from an ESR ASCII raster file (batyhmetry of the South Pacific). I also have a point layer (indicating certain seabird breeding locations). For each point I produced 20km buffer zones. So far so good. What I need to do now, is clip the circular buffer zones so that only the buffer regions between two contour lines remain (in my the case 0m and 150m isolines). Is there a way in ArcGIS to clip/remove the buffer zones that are outside the 0-150m contours?


title: Reproject script in Python with GDAL

text: I'm having a lot of trouble with GDAL. Besides just the at times lacking documentation, there seems to be little support in Python. Anyways this is just a batch reproject script that goes thusly:               source_file = gdal.Open(filepath.encode("ascii"))     source_wkt = source_file.GetProjectionRef()     source_srs = osr.SpatialReference()     source_srs.ImportFromWkt(source_wkt)                                      reproj_file = gdal.AutoCreateWarpedVRT(source_file, source_wkt, dest_wkt)     gdal.ReprojectImage(source_file, reproj_file, source_wkt, dest_wkt)     reproj_attributes = reproj_file.GetGeoTransform()          driver = gdal.GetDriverByName("GTiff")     dest_file = driver.CreateCopy(outputpath.encode("ascii"), reproj_file, 0)      The dest_srs and des_wkt part are not defined in this chunk of code but it's somewhere outside the loop (since it's only needs to be defined once). It seems to work once, I can get one nice looking tif out of it, then give me an 'ERROR 6 WriteBlock() not supported' and python crashes. They are all GeoTIFFs created in the same way with the same basic data (just different times). Also due to the nature of the reprojection (from GCS to PCS), the AutoCreateWarpedVRT tends to create a lot of blank space, but gives it a value of 0, which is an issue since that could be a real data value. Is there someway to set the nodata value to -99 instead?


title: Python SQL connection causes crash in ArcMap

text: I have the following code which I have striped down to just included the part that is causing a crash in ArcMap.               import pymssql          # Connect to the Research db using Windows Integrated Authentication.     conn = pymssql.connect(host='000.000.0.0', database='Research', trusted=True)          conn.close()     del conn      The first time this is run as a script from the toolbar works just fine. However if I run the script again it crashes ArcMap to crash. "ESRI has encounter a serious error and will now exit." Then asks if I want to submit error report... Thoughts on what's happening? I thought closing the connection would be enough. Do I need a Try/Exception statework work here. To handle if the connection is already made then don't override it? Also, If someone knows where to find the ESRI report log on the computer I will upload mine.


title: One to Many Relationship in QGIS with PostGIS

text: Does anyone know of any tool that can do a One-to-Many relationship in QGIS. I have my data in a PostGIS database. I would like to be able to query spatial layers that are related to spatial or non spatial tables and vice versa. I believe there was a way in ArcGIS 9.x that allowed to do something similar. If there isn't an existing tool, would anyone be willing to build such a tool? Thanks,


title: Digitize by User Input by Distance and Length in QGIS

text: I was wondering if anyone knows a way to digitize a line or polygon by user entered distance and length. This is an option in ArcMap, but I would like to be able to do this in QGIS. Thanks, Ryan


title: OpenLayers is not defined while accessing over the html and js

text:
![enter image description here](http://i.stack.imgur.com/IHSAA.png) The above is my working directory, all html files are under WebContent folder and the problem is that when all the code is written in a single html file working fine. If I split the OpenLayers code into 'js' files under 'js' folder. The files within 'js' folder couldn't refer OpenLayers library(OpenLayers.js). **ReferenceError: OpenLayers is not defined** So that var's not getting initialized.               {                                <meta name="viewport" content="width=device-width, initial-scale=1.0, maximum-scale=1.0, user-scalable=0">             <link rel="stylesheet" type="text/css" href="ext-321/resources/css/ext-all.css" />                  <link rel="stylesheet" type="text/css" href="css/style.css" />                                                                                }      working with localhost. and calling init() method explicitly from Ext.onLoad() while I am calling with in the same page it works fine and there is no problem with map it is displaying               Ext.onReady(function(){             Ext.QuickTips.init();             Ext.BLANK_IMAGE_URL = "ext-321/resources/images/default/s.gif";             // turn on validation errors beside the field globally             Ext.form.Field.prototype.msgTarget = 'side';             new Ext.Viewport({             layout      :   'border',             renderTo    :   Ext.getBody(),             region  :   "center",             //margins   :   '0 5 5 0',             defaults:   {border:false},             items       :   [{                 region  :   "center",                 layout  :   "fit",                 id      :   "center_region",                 /*//tbar    :   new Ext.Toolbar({//551                                 id: 'map-toolbar'                             }),                 html    :   "

"                 *///555//                 items   :[{                     xtype   :"tabpanel",                     id      :"tab_map_status",                     activeTab: 0,                     defaults    :   {border:false},                     items:[{                                 title   :   'Map',                                 id      :   'tab-map',                                 tbar    :   new Ext.Toolbar({                                                 id: 'map-toolbar',                                                 items: []                                             }),                                 html    :   "

"                              }]}]});init();});'      Please help me.


title: How to display a image from Postgres database

text: I have Images in bytea form in PostgreSQL/PostGIS database, i want to display them in a webpage using PHP - PostgreSQL. I got lot of examples and I worked out but not getting output from them... Few codes I have used.               $conn = pg_connect("user = postgres password = nkr@123 host = localhost dbname = ARKREJ");              $res = pg_query($conn, "SELECT photo FROM public.zdata_4326, public.zdata   WHERE zdata_4326.gid = zdata.gid AND zdata.gid = 1;");                  if (!$res) {                 echo "An error occured at result.\n";                 exit;             }                  if($res) {                 $Row = pg_fetch_row($res, '0');                      $Object = pg_locreate($conn);                 header("Content-Type: image/jpg");                 $File = pg_loopen($conn, $Object, "r");                 pg_loreadall($File);                 }             else {echo "The query failed!?! There is no row with id=1!";}      Another one display.php               $conn = pg_connect("user = postgres password = nkr@123 host = localhost dbname = ARKREJ");              if (!$conn) {             echo "An error occured.\n";             exit;         }              $temp = "C:/Program Files/OSGeo/MapGuide/Web/Apache2/htdocs/tmp.jpg";                   if (!$temp) {             echo "An error occured at temp.\n";             exit;         }                   $result = pg_query($conn, "SELECT (photo, '$temp') FROM public.zdata_4326, public.zdata WHERE zdata_4326.gid = zdata.gid AND zdata.gid = 1;");              if (!$result) {             echo "An error occured at result.\n";             exit;         }                  if ($result)                 {                     While ($line = pg_fetch_array($result))                     {                         $ctobj = $line["photo"];                         echo "<IMG SRC = show.php> </br>";                     }                 }             else {echo "File does not exists"; }      show.php               <?     header("Content-type: image/jpeg");             $jpeg = fopen("C:/Program Files/OSGeo/MapGuide/Web/Apache2/htdocs/tmp.jpg", "r");                  $photo = fread($jpeg, filesize("C:/Program Files/OSGeo/MapGuide/Web/Apache2/htdocs/tmp.jpg"));                  echo pg_unescape_bytea($photo);         ?>      Can anyone please help me from this... Thanks in advance..


title: How to get x and y coordinate values of an annotational featureclass

text: I am able to get the x and y coordinate values from point, polyline, polygon and multipoint type featureclasses. However, I am unable to export x and y coordinate values of an annotational featureclass programmatically. Please give me any code example. Thanks.