title: Are polygons created in Arcpy built in any order?

text: My script uses a search cursor to evaluate polyline coordinates and build polygons based on those coordinates. The script works well, but I am curious as to whether anyone knows the order (if there is any) that polygons are built. I know that the output polygon feature class reads the rows (and their coordinates) in order, but I do not know if they _build_ them in order. I could see the polygons being built according to which order they were read, or based on size, etc. Any input would be appreciated. Thanks in advance.


title: Is there any tool to reproject x,y coordinates, apart from Proj4?

text: I have a JEE project (Java, JBoss) and I'm looking for a library that allows x,y coordinates reprojection. I know about Proj4. Do you have any other option? Is it possible to find on the net a reprojection algorithm that I could use? The reprojection I'm interest in is the French NTF Lambert II Etendu to RGF Lambert 93.


title: Downloading multiple files from TIGER FTP site

text: ![TIGER FTP Site](http://i.stack.imgur.com/zpIOq.png) Hello everyone! I've been downloading TIGER data per county from the FTP site one at a time, with some states having hundreds of files. I was curious if there was a faster way to download maybe all the files for each county for each state at once? Thanks for your time and help! Attached is a picture of what i'm dealing with, would like to eventually download all the files for each state at once versus one at a time.


title: Joining features based on location and attribute at the same time

text: I have a task of assigning data from database that I represent as point layer, and containing information about characteristics of soils through the world, and also soil type. I also have polygon layers for each lattitude, that contain area of each type of soil (same ones as point layer, therefore can be ussed as join feature). Inside each lattitude layer, each soil type is only one, multipart polygon. I need to assign values from point layer in a way that the soil type always match, but also that it uses point closest geograpicly. If I use just join by attributes, it doesnt take into account spatial location, and vice versa, and I really have trouble figuring how I could approach that. Each layaer has around 40 soil types, and so doing it manually would take ages. I can also add that majority of points are way off the polygons that I intend them to join to, but as long as the soil types match and points is closest, its what I need. I have basicly no knowledge about python scripting. I tried making model to handle that, but its also beyond me. Thanks in advance


title: TileMill Exports Failing

text: I'm having trouble exporting my map in TileMill (v0.10.1). I'm using OSM Bright with a PostGIS database connection to import OpenStreetMap data. When exporting an SVG or PDF through the user interface in TileMill, the export fails more often than not with a SIGABRT error. The console shows the following error:               [tilemill] [tilemill-ui] Rendering file     [tilemill] [tilemill-ui] libc++abi.dylib: terminating     [tilemill] Export process failed with signal 'SIGABRT'      I then tried exporting through the command line by running `./index.js export OSMBright ~/Documents/MapBox/export/cityhall2.svg --format=svg --width=2000 --height=2300 --bbox="-79.8786965, 43.2604005, -79.8702965, 43.2516005"`, with a nearly identical error. It usually exports a blank file, although once in a while it does export nearly what I want (with zoom levels and symbology slightly off, but still usable). Has anyone had this problem before? Thanks for the help!


title: Ecoregions shapefiles: subsetting, extracting and summarizing data from raster

text: I am attempting to subset and project an ecoregions shapefile from : EPA ECOREGIONS my data set is found here: SHAPEFILES and RASTER When I run this code, two major problems occur: 1) I cannot get the subset of ecoregions to plot on top of my raster file and 2) I cannot extract the values (summary information) from the subset of ecoregions from the raster. Can anyone help with this?               require (raster) raster<-raster('FG_RCP26_2050.tif')      plot (raster)      ![Raster with Data](http://i.stack.imgur.com/5UQMj.jpg) ecoregions.all<\- readOGR("NA_CEC_Eco_Level2.shp",layer="NA_CEC_Eco_Level2") head (ecoregions.all) plot (ecoregions.all) ![ecoregions](http://i.stack.imgur.com/DthVm.jpg)                   #this is the subset of ecoregions that I need        NA_L2<-c('5.1','5.2','5.3','8.1','8.3','8.4','8.5')                     #subset the ecoregion.all to just those in the northeast      #this must be where the problem is      ne.ecoregions = coregions.all[match(toupper(NA_L2),toupper(ecoregions.all$NA_L2CODE)),]                    #plot ne.ecoregions on top of raster and extract summary statistics for each ecoregion      #summary statistics =distribution of values within each ecoregion     plot (raster) plot (ne.ecoregions, add=T)          ecoregions.raster<-mask (raster, ne.ecoregions)      summary (ecoregions.raster)


title: JSTS : Union on polygons with holes

text: ![JSTS union](http://i.stack.imgur.com/CcUPL.png) Hi all, I'm trying to union two polygons (see above) with JSTS. My first polygon contains an hole, and my second polygon fill partially this hole. The expected result is a unique polygon with two hole (or one hole with two points at the same location? I've tried with PostGIS and the result is two holes). The actual result is what you can see on C: my first polygon has no hole anymore. Is it a defect of the JSTS tool? Do you have an idea on how can I fix it? Or is it just some configuration to add somewhere? I'm using the JSTS union function like that: firstGeometry.union(secondGeometry); Please ask if you need more information. Thank you for time reading me! Sébastien EDIT: Thanks for your answers! @WhiteboxDev : I don't think I have topological errors. I'm easily able to reproduce this case by drawing polygons on OpenLayers, and my polygons seems to be correct. Please note that I'm using JSTS and not JTS. I've tried with the buffer(0,0) function, but I have unfortunately the same result. Here is an example of a union that fails: POLYGON((335179.0532258 328356.69669355,335178.1782258 328437.63419355,335285.3657258 328437.19669355,335284.9282258 328355.82169355,335179.0532258 328356.69669355),(335213.6157258 328373.32169355,335198.3032258 328400.88419355,335225.4282258 328426.69669355,335267.4282258 328415.75919355,335262.1782258 328371.13419355,335213.6157258 328373.32169355)) POLYGON((335213.6157258 328373.32169355,335225.4282258 328426.69669355,335267.4282258 328415.75919355,335213.6157258 328373.32169355)) And the result of a JSTS union() between these two geometries or a buffer(0,0) on a new geometry collection containing these geometries is the following polygon : POLYGON((335179.0532258 328356.69669355,335178.1782258 328437.63419355,335285.3657258 328437.19669355,335284.9282258 328355.82169355,335179.0532258 328356.69669355)) EDIT (2): It happens in each case where I have   * 1 to many hole(s)   * 1 to many polygon(s) within this/these hole(s), creating "sub-holes"   * Sub-holes are sharing a vertex Doesn't seems to happen in other cases. Geometries are topologically correct. EDIT (3): I've found a defect about that on GitHub: https://github.com/bjornharrtell/jsts/issues/135 I'm still looking for a way to fix it!


title: Editor and projections ArcGIS

text: When opening an edit session Arc is telling me that some of my shape files' coordinate systems do not match the data frame. Here is the data frame coordinate system: > Stereo_70 > > Projection: Double_Stereographic > > False_Easting: 500000.000000 > > False_Northing: 500000.000000 > > Central_Meridian: 25.000000 > > Scale_Factor: 0.999750 > > Latitude_Of_Origin: 46.000000 > > Linear Unit: Meter > > GCS_Dealul_Piscului_1970 > > Datum: D_Dealul_Piscului_1970 And here is the coordinate system that all the shape files in question have: > Projected Coordinate System: Stereo_70 > > Projection: Double_Stereographic > > False_Easting: 500000.00000000 > > False_Northing: 500000.00000000 > > Longitude_Of_Origin: 25.00000000 > > Scale_Factor: 0.99975000 > > Latitude_Of_Origin: 46.00000000 > > Linear Unit: Meter > > Geographic Coordinate System: GCS_Dealul_Piscului_1970 > > Datum: D_Dealul_Piscului_1970 > > Prime Meridian: Greenwich > > Angular Unit: Degree Now is it just me or is anyone else having difficulty seeing the difference? Doe this mean I'm going to have to reproject the files in question FROM their original projection, TO their original projection? I'm not even sure that's possible but will try now. And short of that I guess I just wipe the projections and define them again.


title: ArcGIS Server 9.3 - Javascript API - Selecting Feature Problem

text: I have ArcGIS Server 9.3, and I'm trying to use the JS API. I successfully displayed the map, and I can locate and zoom to a feature based on ParcelID value. The problem is I can't highlight this parcel. I mean I tried to configure the symbol, but its not being displayed. Here's the code:               <script>              dojo.require("esri.map");         dojo.require("esri.layers.FeatureLayer");         var map,featureLayer,query,sfs;            function init() {         map = new esri.Map("mapDiv", {            sliderOrientation: "horizontal"         });                    var imageParameters = new esri.layers.ImageParameters();         imageParameters.format = "PNG24";  //set the image type to PNG24, note default is PNG8.              //Takes a URL to a non cached map service.         var dynamicMapServiceLayer = new esri.layers.ArcGISDynamicMapServiceLayer("http://*******/*****/rest/services/****/MapServer", {           "opacity":1,            "imageParameters":imageParameters         });         dynamicMapServiceLayer.setVisibleLayers([4,10,11,16]);         map.addLayer(dynamicMapServiceLayer);                   //Locate Parcel--------         sfs = new esri.symbol.SimpleFillSymbol(esri.symbol.SimpleFillSymbol.STYLE_SOLID,new esri.symbol.SimpleLineSymbol(esri.symbol.SimpleLineSymbol.STYLE_DASHDOT,new dojo.Color([255,0,0]), 2),new dojo.Color([255,255,0,0.25]));              featureLayer  = new esri.layers.FeatureLayer("http://*****/*****/rest/services/*****/MapServer/4", {                 outFields: ["*"],                 mode: esri.layers.FeatureLayer.MODE_SELECTION             });         featureLayer.setSelectionSymbol(sfs);         map.addLayers([featureLayer]);              query = new esri.tasks.Query();                 query.where = "ParcelID_NS = '101010101'";              var deferred = featureLayer.selectFeatures(query, esri.layers.FeatureLayer.SELECTION_NEW, function (selection) {             var stateExtent = selection[0].geometry.getExtent().expand(3.0);             map.setExtent(stateExtent);         });         //----------------------           }                 dojo.ready(init);          </script>      Any help? Thanks.


title: Why does the attribute Table not maintain input data?

text: for a longer time I have a problem with the attribute table. I create the polygons on the map, I fill up the table with attributes. Then, if I click "save", all the data from the table are being deleted and in all fields of the table a word "error" does appear. After the restart of the program, only that columns containing the number values are being kept. That ones which contain the words are deleted and the word "null" does appear in all the fields. Is somebody able to help me?


title: Which SRS is computational more effective to store continental scale data? (given some conditions)

text: For a Precision Agriculture application I need to make calculations in UTM projection. The data will be stored in Postgis and is made of soil samples results (point data) and field boundaries of many farms. But the farms are scattered through many utm zones (ex. a farm is in Z23 and the other is in Z22, so I can't store the data in UTM, it needs to be stored in a SRS without zones. Since the calculations will be made in UTM the data will always need to be converted back and forth to the storage SRS. So which of those inputs is cheaper in terms of computations of SRS conversion when the output SRS is UTM Wgs84:   -Geometry vs. Geography Postgis types   -Geometry projected (ex World Mercator) vs. Unprojected (latlon)  I think that it involves the internal Proj4 math but I couldn't understand it myself...


title: Is 45005000 a recognisable location?

text: Hello and sorry if this is not a relevant question to this SE. Please if it is not, mitigate it to where you see fit. I have recently purchased an item online and I have been provided by the seller with a tracking number applicable to this site: http://track-chinapost.com/startairmail.php . So when I check where my item is, I get Location: 45005000 Any ideas what this could mean?


title: Geometry issues when creating table from existing spatial table

text: In short, I want to split an existing spatial table into two parts. The split itself is straight forward. I start with this (fully functioning) table:               CREATE TABLE table     (       my_identifier smallint,       geom geometry,       CONSTRAINT enforce_dims_geom CHECK (st_ndims(geom) = 2),       CONSTRAINT enforce_srid_geom CHECK (st_srid(geom) = 4326)     )     WITH (       OIDS=TRUE     );      and create the new tables by:               CREATE TABLE public.table_part1 AS     SELECT * FROM public.table WHERE idmy_identifier = 1;          CREATE TABLE public.table_part2 AS     SELECT * FROM public.table WHERE my_identifier = 2;      However, the problem is that although the geom column gets transferred to my new tables, the geometry (?) itself doesn't get transferred correctly:               CREATE TABLE table_part1     (       my_identifier smallint,       geom geometry     )     WITH (       OIDS=FALSE     );      I tried running `SELECT Probe_Geometry_Columns();` but it didn't help. Do I need to insert data into the geometry_columns table? If so, exactly what?


title: Delete polygon underlap?

text: I am editing a layer in QGIS with severe polygon overlaps. Is there a way to delete the underlapping portions (i.e. the invisible parts) using QGIS? If I import the layer into GRASS using the GRASS-plugin, I can get a very useful layer which contains the overlaps (polygons_2), but if I use the difference geoprocessing tool to slice the overlaps from the original file, it naturally cuts away the overlap as well as the underlap, leaving holes instead of overlaps. I have also tried the v.clean functions in grass, but the overlaps are really too big to use snapping without distorting the data.


title: How to transform Invalid Date format in Data Interoperability

text: My workbench is running into some date data it is rejecting. Here is the error I'm receiving... The invalid value '20100419000000' was supplied to the field 'install_date_orig' I believe the feature is rejected because of the '000000' for the time portion of the date. Any suggestions of which transformers may be useful to get around this issue. Here's an list of of what I've tried so far.   * DateFormatter with this format, FME Date (%Y%m%d)   * SubStringExtractor to split the field into the year, month and day parts, StringConcatenator to recombine those parts, the DateFormatter as per above, and finally AttributeRemover to exclude the bad field.   * TestFilter and Tester to filter out the bad records.  I've noticed that when it reaches the bad record it fails before it reaches the first transformer. Is there anyway to catch a record with a bad field value before it reaches any of the transformers? I appreciate any suggestions.


Is the "sdetable" command available with ArcGIS Workgroup, or only the Enterprise edition? We have ArcGIS Server Workgroup 9.3.1 and want to update some of the feature class tables with data from another system using SQL UPDATE commands. From taking a look at the documentation, we must first create a multi-versioned view, and then run the SQL UPDATE commands against this view. The problem is that we do not have the "sdetable" command required to create the multi- versioned view. We think that these command line tools may be available only in the Enterprise edition. Is there another way that we can perform basic data updates with SQL commands in the Workgroup Edition? I understand that the correct way to do this is not by editing the default version directly. EDIT:   Our data is versioned for different uses, including some field applications. We were planning on creating a "sql-updates" version just for making updates to attributes in a feature class. If I use SQL UPDATE commands to update this version, my concern is that the changes will not get replicated to the other versions, or the data could become corrupt. This seems like it would be bypassing delta table entries. From what I have read, we should use "sdetable" to create a multiversioned view, and then use stored procedures create_version, set_current_version, edit_version (1 - start session), run my updates, edit_version (2 - close session), and finally delete_version (after reconciled). I see these stored procedures in the database, but can't do the first step of creating the multiversioned view. Maybe it would be possible to perform these steps on the "sql-updates" version without creating the view. **EDIT 2:**   Since the original question has been answered, I opened a new question:   How can we use SQL to update feature class attributes in a versioned geodatabase in ArcGIS Workgroup


title: Displaying road outlines but not segments i.e. no lines at intersections?

text: ![enter image description here](http://i.stack.imgur.com/uuE0F.png) I have the roads that I want to just display at outlines but not segments (I.e not disconnected from each road). I am trying dissolve but its not working. Attached is a link to a package file of the data. Please help! https://www.dropbox.com/s/55hqr9471v9qxwf/Road%20Test%20package.mpk


title: QGIS composer - How can I eliminate striping across transparent raster layers?

text: I have downloaded QGIS dev version 1.7.0-114 and I have come across some weird behaviour in the composer. I create a map with a hillshade and a transparent DEM ontop of it, pretty standard stuff. When I export this map to an image format I get two strips across my output, this happens for a number of formats (png, jpg, bmp, tiff) and the stripes seem to be in the same place everytime. I have attached an example below (sorry for the size of the attachment). ![enter image description here](http://i.stack.imgur.com/RKFhO.jpg) The stripes appear in the same position on the page every time and this happens with different raster datasets, so I have ruled out the data. Also depending on the page size you may get 1 stripe (A5 as shown above), 2 stripes (A4) or no strips (custom dimensions h100 x w100). Is this something that can be overcome via some settings which I may have overlooked? Or might this be something funny in the software?


title: how to create Multiple layers?

text: I'd like to create multiple layers for asset management using openlayers, could anybody help?


title: Can ArcGIS layer be shown on Google Earth?

text: If we have a table in SQL Server which was used as layer in ArcGIS. Can the same thing be done in Google Earth? I mean does Google Earth support adding custom layers? If so how?


title: Open Layers 2.13.1 Image Layer Zoom Animation Issue

text: I updated my web app to use 2.13.1 from 2.12 and see that Image layers now animate zoom level changes? I don't believe it previously did this and investigating the code only briefly, I'm not sure how this even happens. I also do not see any obvious way to specify that zoom level changes should animate or not with Image Layers. The reason I care is that I'm using a Google base layer with an Image layer on top, and zoom changes animate for the Image Layer but not really the same way with the Google layer. The result is that the Image overlay does not appear to be properly attached to the Google layer through the zoom animation and it looks wrong. I experimented with the "animationEnabled" option on the Google layer but it doesn't improve the situation I'm not sure that I've explained this very well so I created a quick screen cap move of it in action. Unfortunately screen recorded slowed my computer down so the movie isn't as smooth as what I see without screen recording but it still shows the issue. http://youtu.be/i53UoxOmVQc I simply want all layers to have synchronized zoom animations. Any ideas how to fix?


title: Geoserver WPS - Tiled Barnes Surface with geospatial index

text: I'm using Barnes Surface to compute a set of points from my PostgreSQL db. The problem is that my whole table is processed for each tile. I added a geospatial index on my points, how can I tell geoserver to compute only the points present in the bounds of the tile?


title: PGRouting Layer plugin error in QGIS 2.2

text: I am trying to do some basic computations using the PGRouting plugin for QGIS, but I am getting this error all time:               An error has occured while executing Python code:          Traceback (most recent call last):       File "C:/Users/Gonzalo/.qgis2/python/plugins\pgRoutingLayer\pgRoutingLayer.py", line 407, in run         if db and db.con:     UnboundLocalError: local variable 'db' referenced before assignment          Python version:     2.7.5 (default, May 15 2013, 22:44:16) [MSC v.1500 64 bit (AMD64)]               QGIS version:     2.2.0-Valmiera Valmiera, c3a2817          Python path: ['C:/Users/Gonzalo/.qgis2/python/plugins\\processing', 'C:/PROGRA~1/QGISVA~1/apps/qgis/./python', u'C:/Users/Gonzalo/.qgis2/python', u'C:/Users/Gonzalo/.qgis2/python/plugins', 'C:/PROGRA~1/QGISVA~1/apps/qgis/./python/plugins', 'C:\\PROGRA~1\\QGISVA~1\\bin\\python27.zip', 'C:\\PROGRA~1\\QGISVA~1\\apps\\Python27\\DLLs', 'C:\\PROGRA~1\\QGISVA~1\\apps\\Python27\\lib', 'C:\\PROGRA~1\\QGISVA~1\\apps\\Python27\\lib\\plat-win', 'C:\\PROGRA~1\\QGISVA~1\\apps\\Python27\\lib\\lib-tk', 'C:\\PROGRA~1\\QGISVA~1\\bin', 'C:\\PROGRA~1\\QGISVA~1\\apps\\Python27', 'C:\\PROGRA~1\\QGISVA~1\\apps\\Python27\\lib\\site-packages', 'C:\\PROGRA~1\\QGISVA~1\\apps\\Python27\\lib\\site-packages\\PIL', 'C:\\PROGRA~1\\QGISVA~1\\apps\\Python27\\lib\\site-packages\\win32', 'C:\\PROGRA~1\\QGISVA~1\\apps\\Python27\\lib\\site-packages\\win32\\lib', 'C:\\PROGRA~1\\QGISVA~1\\apps\\Python27\\lib\\site-packages\\Pythonwin', 'C:\\PROGRA~1\\QGISVA~1\\apps\\Python27\\lib\\site-packages\\Shapely-1.2.18-py2.7-win-amd64.egg', 'C:\\PROGRA~1\\QGISVA~1\\apps\\Python27\\lib\\site-packages\\six-1.3.0-py2.7.egg', 'C:\\PROGRA~1\\QGISVA~1\\apps\\Python27\\lib\\site-packages\\wx-2.8-msw-unicode', 'C:\\PROGRA~1\\QGISVA~1\\apps\\Python27\\lib\\site-packages\\xlrd-0.9.2-py2.7.egg', 'C:\\PROGRA~1\\QGISVA~1\\apps\\Python27\\lib\\site-packages\\xlwt-0.7.5-py2.7.egg', '.', 'C:\\PROGRA~1\\QGISVA~1\\apps\\qgis\\python\\plugins\\fTools\\tools']      I made a local installation of postgres 9.3.4, later I added postgis using the stack builder, and finally I followed the instructions posted here: PGRouting 2.0 for windows. I also installed python-psycopg2 for my version of python. The SQL query runs well in pgAdminIII, but I can not make that QGIS plugin work, so any help would be really appreciated.


title: Shoehorning hierarchial data into MS SQL ArcGIS SDE

text: So, I've been tasked with creating a set of tables to store hierarchical data within an SDE. This is in reference to utility data, i.e, a set of valves, control panels, actuators, etc (each potentially having their own attributes) belong to a pump, the pump belongs to a process (also having its own attributes), which belongs to a pumpstation (which can have its own attributes) etc etc. The only way I can think of to do this is to basically have a table for each asset (a table each for valves, control panels, actuators, pumps, process, pumpstation) and a relationship class for each one to link to the table higher up in the hierarchy for it, all the way to the top, and then a final relationship class to link the pump station table to an already existing pump station point feature class. I've already googled some to find examples of using hierarchical data in ArcGIS and i'm not having much luck. Anyone have any suggestions if there's a better way of doing this?


title: small difference in projection between tilemill and postgis

text: I would like to create a world map with interaction on continents created directly from a postgis database. For the moment, I have a POSTgis database with the data. I created a background with TileMill. I get data from Postgis in SVG and put them on the map using Raphael. All seem well except that I have a small difference between both map :s I use default projection of TileMill (and I think we can't change it :)). And here is the POSTGIS command I'm executing to get the SVG:               SELECT st_asSVG(            st_scale(                st_translate(                    St_Intersection(                        ST_Transform(ST_SetSRID(ST_GeomFromText('POLYGON((-151.875 80.6470,175.7812 80.6470,175.7812 -63.2336,-151.875 -63.2336,-151.875 80.6470))'), 4326), 3857),                    (SELECT st_union(geom) FROM geo.continentboundaries))                , 20037508.72, -20037508.72)            , 0.0000199043047125,0.0000199043047125)         , 0, 0)      Here is the result: http://i.stack.imgur.com/7lS0F.png Do you know where I have a mistake? Edit: some precisions: I'm using the same data for Tilemill and the svg. I added in the tilemill project the example of tilemill there is no problem of projection. In the layer in tilemill I'm using 900913. st_srid give me 3857 in Postgis which seems ok? I'm a beginner in SIG


title: Exploding overlapping to new non-overlapping polygons?

text: Given multiple polygons that overlap in multiple ways, I would like to export from these features all polygons that don't overlap with others, iteratively. The product would be a number of features with no overlap that, when summed together, make up the original. The products could then be used as input to Zonal Statistics, and this would be much faster than iterating Zonal Statistics over each polygon. I have been trying to code this in ArcPy without success. Does code to do this already exist?


title: What is geom_way column in osm2po generated data?

text: Could anybody tell me the meaning of geom_way column generated by osm2po?


title: Looking for a functionality in QGIS to find all neigbours in a polygon layer

text: I am looking for a functionality in QGIS to find all neigbours for each polygon in a polygon layer. The result should be a table with 2 ID-Columns: one for the polygon and the second for the neighbours. I tried "spatial query" in QGIS with the layer as source and its copy as reference, but to no avail. In Mapbasic of Mapinfo it would be a simple expression:               select t1.id1, t2.id2 from t1, t2 where t1.obj intersects t2.obj       (where t1 and t2 are two copies of the polygon layer). Thanks for any help!


title: Is it possible to automatically set CRS of Rasters from PRJ file?

text: I am trying to learn my colleagues to display measured data (ASCII) in maps relatively easily. I can prepare project, color and symbology style for the point layer and background map. I also found good way to get background map (PNG + world file) for offline use via MOBAC and OSM. MOBAC uses EPSG:3395 WGS 84 / World Mercator and I have the appropriate PRJ file. However, I found out, that QGIS does not use the PRJ for rasters and load them always using default CRS (in my case EPSG:4326). Is it possible to make QGIS load the CRS settings from the PRJ (or other file) without need to set the CRS manually for every raster? thanks


title: How to produce "SpatialPixelsDataFrame" from stacked RasterLayer?

text: How to produce "SpatialPixelsDataFrame" from staccked RasterLayer ? I have a stacked rasterlayer and I want to convert them to SpatialPixelsDataFrame. RasterLayer information is               > RAD.all[[1]]     class       : RasterLayer      dimensions  : 180, 240, 43200  (nrow, ncol, ncell)     resolution  : 0.01, 0.01  (x, y)     extent      : -123.305, -120.905, 36.705, 38.505  (xmin, xmax, ymin, ymax)     coord. ref. : +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0      data source : C:\Users\Michael\Documents\1998bil\1998ASC5min_ppt_spas1214_0001_19980202_0810_UTC.asc.bil      names       : X1998ASC5min_ppt_spas1214_0001_19980202_0810_UTC.asc


title: Near analysis by matching attribute in two different layers

text: I need to perform a near analysis to assign cases in one feature class to their nearest match on an attribute in a different feature class (ex: Tutors and Tutees need to be paired by same language spoken). I have not been able to find a way to do this. The closest I have reached to simplifying this process is to perform a split by attribute on the values for the matching variable. So instead of having all tutors in one feature class, I now have n language number of feature classes containing tutors and n language number of feature classes containing tutees. I then plan on performing a near analysis between each feature class. However this is extremely tedious and may result in 100+ near analyses needing to be performed because there are other criteria I will need to match independently on in a similar way. Is there a better approach? Much thanks in advance. Edit: The near must be performed between two different layers to preserve the different roles of individuals being matched. I am new to learning modelbuilder so any ideas on conditionals, or iterators that could determine the "nearest" match on a shared attribute for points housed into two different layers would be much appreciated!


title: Mapping elements in a warehouse

text: I am new to GIS. For a small project I want to map elements in a warehouse and want to find the location of the elements later (e.g. nails are located at x,y). Now I am wondering what kind of position I should choose. Should I choose an absolute location coordinate system like `1,3`, or should I use real geocoordinates like `-7.023, 32.2423` for my use case? Where are the benefits/disadvantages with each approach?


title: Onclick event with polygons

text: When I click on any polygon in the openlayer I need to show the respective information in the popup. Right now I am using the below code.               eventListeners: {         getfeatureinfo: function(event) {             window.open("http://mydomain.com/polygoninfo.php?polygonid=" + event.feature.attributes.id, "window_name","location=1,status=1,scrollbars=1,resizable=no,width=650,height=650");         }     }      The issues with the above code are:   1. The popup is opening on complete layer, meaning that the popup is opening only on polygons if I click some other parts of the layer where the polygon is not present.   2. I need to capture the unique ID of the polygon when the polygon is clicked. I am using sld for labels, etc. Is there a way to create a uniqueId of polygon inside sld like:                       <sld:Label>                <ogc:PropertyName>CAGE_NO</ogc:PropertyName>              </sld:Label>      Sorry if any thing in my previous question is unclear or confusing. So, please let me elaborate on the issue as below:   1. As per my understanding, the actual requirement/task is to provide a layer which represents various states in the India map. Once the user clicks on any of the polygons/states, then the request should be redirected to the respective state's home page.   2. My concern is when a user clicks on any of the polygons/states, how to capture or send the respective state's ID/name so that I can redirect to the respective state's home page based on ID/name.


title: Is it possible to create map packages for ArcMap version 10 using 10.1?

text: Is it possible to create a map package in an older version? Currently I am using 10.1 and the client uses 10.0. Can you create a map package for different versions? I see you can create layer packages using Python for different versions but you cannot specify the version when creating a map package. As I don't have 10.0 I can't test whatever I create these. When you unpack the file using ArcCatalog does it create it v.10. and 10.1? Can you create the original file in 10.0?


title: Draw lines between two points in QGIS

text: I have an excel (Fig) which has both X & Y in projected coordinates (WGS84/UTM). Is it possible to draw lines between the X & Y in QGIS ? I've tried to create points based on **X,From_y** and **X,To_y**. I've created two points by just loading them using add delimiter layer (pasted To_Y below From_Y and loaded it). Then tried with Points 2 paths. But it didn't work out. Anyone know how to fix this ? All i need is traverse line between the points ![enter image description here](http://i.stack.imgur.com/hWmhd.png) Edit: This is what i am looking for. I've coordinates(X1, a1 and X1, a2) for n number of points which has to be connected by a vertical line ![enter image description here](http://i.stack.imgur.com/o7gRz.png)


title: Parse Field to Select Records

text: I wish to select records based on the NRN attribute in my dataset : `"NRN" IN ('1','2','23','31')` Some records have multiple values in the NRN field separated by the underscore character.   For example : `20_21, 14_23` etc How do I select any records that contain the values that I want ignoring the records that have part of the string but not the whole number? Example : `select 14_23 but not 20_21`.


title: How to select two points on a map using QGIS?

text: I need to select two points on a map in QGIS. I selected the element selection tool in menu view and I want to select two points. How can I accomplish this? Are there other tools or methods to select points in QGIS?


title: How to create one shapefile from multiple rasters?

text: Ok, I have about ten maps that connect along their edges to each other linking rivers and properties. When I link them all in Inkscape I end up with a large file that won't save (and a smaller one with 5 of the properties saves but is too big and slow in QGIS). So I am thinking that I will load them, about three linked together at a time, into QGIS and create my shapefiles. Is it possible in QGIS to then merge/link the shapefiles along their edges so that the rivers and properties line up and create one large shapefile? additionally, is it possible to then take this large shapefile and georeference it to a base map or can that only be done with rasters? thanks


title: On the usage of Proj4Js->tmerc.js

text: I am using "proj4js & open layers" for my mapping. I got projection & datum (which is based on "Transverse Mercator") of a local plane and this plane is part of a grid which returns coordinates of local points in the form of northings and eastings. Now, as we know that Open Layers accepts lattitude and longitude hence I am doing transformation of northings,eastings of a point to lat,long using transverse mercator (tmerc.js) script of proj4js...So my question is will this be enough for the right transformation or there is any thought or suggestion you would like to give regarding the transformation script I am using ? Thanks


title: Merge nearby items into a single representation

text: I'm investigating OpenLayers and I'm trying to see if it's possible to filter down the number of displayed feature items. Is it possible to display many nearby items as a single point item on a map when zoomed out. So instead of hundreds of nearly overlapping points I can filter them down a few dozen consolidated points? Or is this something that has to be in the Spatial server?


title: How do I set custom spatial reference/coordinate sys to my shapefile

text: I have a python code which converts JSON text file to a shapefile. I have to set coordinate system to this shapefile created. I do not want to access the .prj file from my computer as this code will be embedded in a web toolbox.


title: Google Maps - Use fitBounds on KML polygon fusion table layer

text: I have seen example here how to use map.fitBounds on fusiontablelayer here and here , but I think it works with markers table only. Its impossible to parse the lat and long when the vertex coordinates are in the geometry column as in a kml polygon layer. Is there a solution to use fitBounds or other function to zoom to the extent of a kml polygon fusion table layer?


title: ArcMap: Calculate Quarter using Date field

text: I am using ArcMap 10.0 and making a map where I symbolize on annual quarter (i.e. Q1 2013, Q2 2013, and so on by quarter through 2015). If a date field contains this: 01/15/2013, and I wanted to create a new field called QUARTER, that would contain the quarter in which the date falls, how could I do that? Right now I'm just definition querying a range of dates, but surely there is a way to just calculate Quarter, right? Thanks in advance!


title: How do I use st_dwithin with meters?

text: I am trying to follow PostGIS ST_Buffer Radius Help to use a 50 meter distance but postgis is returning an error. I am trying the following SQL,               SELECT id FROM table      WHERE ST_DWithin(ST_GeogFromText('SRID=4326;POINT(-3.165356 55.926665)'),            geography(geometry), 50);      The hints indicated at ST_GeogFromText are,   * function st_geogfromtext(unknown) does not exist at character 92   * No function matches the given name and argument types. You might need to add explicit type casts. Any help/alternatives would be appreciated.


title: Natural neighborhoods terminology

text: Sometimes the hardest part of an analysis is knowing what something is called. What are some R packages, but more importantly, terminology I should be searching for to define contiguous city neighborhoods based on home price changes - basically a topographic map treating these price changes as elevations (and ignoring how real estate people have defined neighborhoods). I'm not sure if this is called "natural neighborhood delination" or "spatially restricted clustering" or just "grouping analysis".


title: How to change data type in Excel file before creating feature class?

text: I am trying to import data from Excel spreadsheets where each point coordinate has a 'datetime' timestamp (e.g. "2012-06-01 14:45:32") associated to it. When I check the table properties ArcGIS 10 always recognizes the timestamp field as 'date' (i.e. "2012-06-01") and drops the time-of-the-day info. Is there a way to import the timestamp field as a full datetime type? And if yes, how do I go about it?


title: Oceans Baseline Data

text: I am working on a volunteer project and we are in desperate need of a baseline layer. I have been unable to find it anywhere (there are plenty of shoreline and coastline layers out there) does anyone have a source for this data? Thanks, Dave   **Added:**   The baseline from which the territorial sea is measured is mainly the low- water line along the coast as marked on large-scale charts officially recognized by the coastal state.


title: VB.Net Field Population Methods

text: Disclaimer: I am basically brand new to VB.NET and ArcObjects. I am computing some statistics as a component of a larger arcobject script. I write these components to an ITable. This ITable is added to the map display as a standalone table. I find that I am duplicating a fair bit of code to assign a value to each field in the table. This is not something I would do in Python and wonder if a better methods applies in the VB.Net world. VB.NET CODE:                  'Populate the stats table with the iteration.         Dim row As IRow = tTable.CreateRow         Dim fields As IFields = row.Fields              Dim fThres As Double = fields.FindField("ThreshDist")         row.Value(fThres) = lStats(0)         row.Store()              Dim fCountR As Double = fields.FindField("CountR")         row.Value(fCountR) = lStats(1)         row.Store()      where lStats is a list of doubles, and tTable is an ITable object. In python I would simply use a dictionary to map the index (as key) to the fields name, as value. A simple for loop would take care of population using an insert cursor. Does a similar paradigm exist using VB.Net? Am I missing something obvious in the documentation?


title: pgRouting on RHEL6: issues with driving distance library

text: I installed pgrouting on RHEL6 via yum and the ELGIS 6 repo (the only place where I could find pgrouting for RHEL6). The problem is that the `librouting_dd.so` library (for driving distance) is missing. In other words, I have `librouting.so`, `librouting_tsp.so`, but _not_ `librouting_dd.so`. _Does anyone know why this is or how to fix it?_ I've included the yum output below. I also tried compiling pgrouting from source (with the `-DWITH_DD=ON` flag so librouting_dd.so would be built) but this resulted in a waterfall of dependency issues (e.g., it needs CGAL-devel, which in turn needs things like lapack-devel and blas-devel which aren't installing because of other dependency issues, etc., etc.). Any advice would be greatly appreciated.               # yum install pgrouting     Loaded plugins: rhnplugin     This system is receiving updates from RHN Classic or RHN Satellite.     elgis                                                                                | 3.0 kB     00:00          elgis-plus                                                                                | 1.3 kB     00:00          elgis-source                                                                                | 3.0 kB     00:00          Setting up Install Process     Resolving Dependencies     --> Running transaction check     ---> Package pgrouting.x86_64 0:1.05-1_0.el6 will be installed     --> Finished Dependency Resolution          Dependencies Resolved          ============================================================================================      Package         Arch      Version         Repository                      Size     ============================================================================================     Installing:      pgrouting      x86_64     1.05-1_0.el6    elgis                            66 k          Transaction Summary     ============================================================================================     Install       1 Package(s)          Total download size: 66 k     Installed size: 226 k     Is this ok [y/N]: y     Downloading Packages:     pgrouting-1.05-1_0.el6.x86_64.rpm                                                                                |  66 kB     00:00          Running rpm_check_debug     Running Transaction Test     Transaction Test Succeeded     Running Transaction       Installing : pgrouting-1.05-1_0.el6.x86_64                                                                                1/1        Verifying  : pgrouting-1.05-1_0.el6.x86_64                                                                                1/1           Installed:       pgrouting.x86_64 0:1.05-1_0.el6                                                                                Complete!     [root@vm1 yum.repos.d]# rpm -sql pgrouting     normal        /usr/lib64/pgsql/librouting.so     normal        /usr/lib64/pgsql/librouting_tsp.so     normal        /usr/share/doc/pgrouting-1.05     normal        /usr/share/doc/pgrouting-1.05/README.routing     normal        /usr/share/postlbs     normal        /usr/share/postlbs/matching.sql     normal        /usr/share/postlbs/routing_core.sql     normal        /usr/share/postlbs/routing_core_wrappers.sql     normal        /usr/share/postlbs/routing_topology.sql     normal        /usr/share/postlbs/routing_tsp.sql     normal        /usr/share/postlbs/routing_tsp_wrappers.sql


title: Point to point distance over specified distance

text: I need to calculate the distance between neighbouring points which conform to a staggered line. I have two point datasets, one being the baseline and another the monitor. The points typically reside ontop of each other within a few metres. I need highlight the ones that fall outside a specified range. The specification is that one point must have a neighbouring point no further than +/- 9m away from it within 1km sections. The average cannot be over 12m as well. So I need to calculate the distance to its closest neighbour and output this into tabular form thus highlighting where the second feature is out of 'spec' with an appropriate colour. It is likely that only a few will be out over many 1000's. Any help appreciated.


title: Qgis - need to draw exact lines (mtrs) with fixed angle

text: I use GIS 1.8. I have the need to be able to draw exact lines and shapes on a bing map. (This is to project a yet to be build greenhous project on a map to show what i would look like). I have found the tool of Olivier Dalang (improved polygon capture) to be very usefull, since this tool enables me to draw lines while specifying the exact length and angle of the lines, and by doing so i build all walls of the greenhouse on the map). Here is my problem. The plugin enables me to enter the length in 'units', and drawing a line of 1000 units makes me come up with a line of about 622 meters. I cannot find the exact way to change or find out how to draw in meters. Any help highly appreciated! (If there are other tools for drawing 90 degrees angled lines to a map I welcome these too!) Rene van der Velde


title: QGis: Merge intersecting polygons into one which are part of the same feature

text: I'm not a GIS professional so I hope I use the correct wording, at least I hope you understand my question. Let's say there are 5 polygons in the feature A. In the shape file there are more features. In feature A 2 polygons from the 5 existing are intersecting. I want to merge the 2 so overlapping lines will be deleted and the result is one polygon. A possible solution would work like this   * I select the feature which 'activate' all 5 polygons    * I switch to edit mode    * I select the points of the 2 polygons    * I use function MERGE-SELECTED-POLYGONS So far I found only functions to merge objects (features?), merge shape files, and so on. Is there a function MERGE-SELECTED-POLYGONS?


title: ArcGis Trace Tool for QGis?

text: is there any tool in quantum gis that traces a line like the trace tool (Editor Toolbar) in arcgis? I want to draw a line and trace this line back so that i get the beginning and the end of the line at the same position. Any suggestions? Thanks! Kind regards Benjamin


title: Line segment connectivity analysis in FME

text: I have a set of line segments, all of them snapped at the endpoints: ![Groups of connected lines](http://i.stack.imgur.com/9tTz3.png) As seen in the above image, the line segments appear in groups (4 in this example). **Question:**   I would like to determine which line segments are connected to each other and assign an ID (in this example the values 1, 2, 3 or 4) to each line segment describing which group it belongs to. How can I do this in FME?


title: How do I prevent QGIS 1.8.0 changing my dbf table values when joining tables?

text: Im a newbie on GIS as well as QGIS. I have a table of data saved from Excel as dbf format. After dragging and dropping the file into QGIS and then opening the attribute table, all values less than 1 in my original table have been changed to 0. How can I prevent this?


title: violates check constraint

text: I try to add a constraint based on geometry type to a "street" table that prevents anything except Linestring geometry or Null.               alter table streets     add constraint streets_geom_linestring_chk     check(st_geometrytype(geom) = 'ST_LINESTRING'::text OR geom IS NULL);      When I try to store geometries in the table:               update streets set geom=st_geomfromtext('LINESTRING(0 8,12 8)') where id=1;      I got an error, that the new row for relation violates check constraint.


title: Can I change the symbology of selected features only?

text: Is there a way in ArcGIS (9.2, ArcMap) to change the symbology of selected features, over multiple field values? My example is, I've run an attribute query on the 'landscape categories' field of my shapefile, to find every feature in a category which includes the letters 'forest'. This includes categories which are essentially the same or similar enough that I want to display them as being the same, such as 'forestry', 'areas of forest', 'forest', 'forested hills'.... etc etc, there are loads. I can select them all using this attribute query quite easily, but I then want to make them all green for example. I thought I'd be able to do this in my Properties > Symbology > Categories > Unique Values tab, by just right- clicking and selecting 'Properties for Selected Symbols', but apparently selecting features through an attribute query doesn't count as selecting symbols in this tab, nothing is selected? Does anyone know an easy (code-free) way of doing this? I know I could export my selection as a new layer, but I want to keep all my info as one shapefile and I'll need to do this for quite a few categories. I'm also not keen on the idea of adding a field in the attribute table for 'colour' etc, I'm trying to keep the original shapefile as clean as possible. Any help much appreciated! Thanks.


title: map can only have one feature, after user click

text: I'm using openLayers 2.13.1. I have a map wih a vector layer and Editing tools. I want the map to have only one feature even if the user keeps adding. So I want **the new feature to automatically erase the previous one (if a previous exists)**. That way there is only one feature in the map, the new one that the user adds (using Editing tools) I tried some things , but I end up erasing the feature, as soon it is added. How can I do this? Thanks in advance Here is my code               myLayer.events.register('featureadded', this, feature_added);          map.addControl(new OpenLayers.Control.EditingToolbar(myLayer));          function feature_added(feature){          //erase previous features here (?)     pins.events.unregister('featureadded', this, feature_added);              }


title: How can I create random concave polygons for test purposes?

text: In the context of generating random test data using QGIS (in a shapefile): I can generate random, continuous (adjacent) polygons by first creating _n_ random points using ftools and then generating a Voronoi diagram from those. Using this technique I get _n_ convex polygons within some extent. Additionally, I know that all the polygons are contiguous and it is possible to build an adjacency matrix. Is it possible to generate random convex **and** concave polygons that are (1) bounded by some defined extent and (2) continuous? My initial thought is to generate _n_ x 3 polygons and then begin randomly merging them (using a Python script) until only _n_ remain. This _should_ generate at least a few convex polygons. Does QGIS (or GRASS) offer a means to do this? How about an external package? Edited for @whuber's comment:   1. It is possible to trace a path from any polygon to any other polygon without exiting a polygon in the set, e.g. the contiguous U.S. Adding the Hawaii multipolygon would fail the contiguity test because one has to exit California (for example), enter unbounded planar space, and then re-enter the polygon set.    2. The probability distribution of the random points is not important, in that, the density of polygons generated from the randomly distributed points does not matter. The goal is simply to vary the number of neighbors each polygon has - the distribution of the neighbor cardinality is not important, but it should not be equal across all polygons as with a regular lattice.   3. This is a test of regionalization algorithms using randomly generated data These tests have used lattices and Voronoi diagrams. I want to explore a particular algorithm on test data using convex and concave polygons. The generation of convex polygons that meet my constraints is handled using the technique I describe above.  Thanks!


title: How to rename attributes/fields in a shapefile attribute table?

text: I would like to rename some fields in my attribute table: ![enter image description here](http://i.stack.imgur.com/bXfkM.png) Like rename `XRF_N3`in something else. It seems quite a easy issue but I don't find any information about how can I do that... I've QGIS 1.8 in Mac Mountain Lion.


title: Batch processing Topojson

text: Can i convert more than one shape file .shp to topojson. Can this be performed using batch processing.?? I can convert shape file to topojson one by one but i have about 6000 files and i don't want to merge them i want them to be separately converted to topojson. I can also batch process ogr2ogr but don't know how to perform that on topojson. Is there any way i could perform this task using a folder of is there any alternate solution to this problem. Thanks Regards Saad


title: ArcGIS for Developers Free Development and Testing Account suspended/expired?

text: I signed up for an ArcGIS for Developers **Free Development and Testing Account** and when I signed into it today, it said it was suspended and or expired? I thought the account was free for life? Do I need to purchase a subscription like you need to do with the 30 day free trial account? I have a story map that I need access to. Is there any way to get my story map transferred to another account?


title: Connect point to closest feature in ArcGIS 10.2

text: I have some points and some features. These points fall outside the features. I want to connect each point with a line to the closest feature. How can this be done using ArcGIS 10.2? I have checked around toolbox but didn't manage to find the right tool for it. Many thanks in advance!DM


title: applyEdits method and File on arcgis feature service

text: I would like to know if it's possible to upload a file through the applyEdits method on a feature service? I have an Oracle database with a parameter "DOCUMENT" of type BLOB. I would like to upload the file in this parameter. Is this possible? IF so, how would I do it? Thanks for any help and sorry for my English, I'm French.


title: How to use “outputConfig” to configure LayerPanel from GXP

text: Can outputConfig be used for functions too - not just configuration? I have this code in geoexplorer.js (using GeoExplorer 2.4.3) which doesnt work:               var GeoExplorer = Ext.extend(gxp.Viewer, {     ...     config.viewerTools = [     ...     }, {         leaf: true,          text: gxp.plugins.Legend.prototype.tooltip,         outputConfig: {            filter: function(rec) {                alert("filter function");                if (record.get("title").indexOf("Europa -hvit bakgrunn") != -1) {                     alert("filter in outputConfig");                     return false;                 } else if ( record.get("title").indexOf("Europa") != -1) {                     return false;                 }                 return false;            }         },                         checked: true,          iconCls: "gxp-icon-legend",         ptype: "gxp_legend"                   }      Is it because filter() is a function and not a configuration parameter? I need to add this filter because some of the layers dont have getLegendGraphics and from the WMS spec you cannot see if a layer has getLegendGraphics or not.


title: QGIS 1.7.4 crashes when using interpolation tool

text: I am currently using QGIS with Lidar data to produce hillshade images. I have split the lidar data in 20 files (31kb). Some of the files are working some are not. I have checked the files xyz are correct and look normal. I can add a delimited vector layer, and everything is correct and save as shapefile. Then I load data for interpolation tool as it is building in may get to 7-10% then crashes. Why is this happening? Also I'm trying to add hillshade relief plugin; I can't find it and can't load load it. There is no link to this plugin at pyqgis.org.


title: How make .qgs files open in QGIS automatically on a Mac?

text: When saving a project my copy of Qgis used to assign file extension .qgis automatically. Now it assigns .qgs automatically, which means that I can't reopen it. All my old projects have also been given new extensions .qgs. Does anyone know what's wrong and how can I fix it? Thanks in advance


title: Are there WMS servers providing maps of UTM zones, actually in the zone projection?

text: There are lots of EPSG:4326 (lat lon) services out there but if one has to use an `UTM projected map` as base layer to display utm scenes, are there any WMS endpoints `publicly` available for that?


I am trying to find out if a certain function is available to link an excel sheet with a kmz. What i need to happen is have a Placemark that pulls information off an excel sheet row and display's the data in the description part of the Placemark. I know that there is a macro in excel to export KMZ file but I want the KMZ itself to perform the task of collecting the information and displaying it. I would guess that it would be an extension of the “ad link command,” but I can’t seem to get my head around it. This would be cut down on data duplication if there was a way of doing it. If you know of a way to do this please comment. Thanks Shane


title: Progress bar percentage is increasing and decreasing, never finishing

text: I am trying to clip a raster by shapefile perimeter using the Clip tool. The percentage has been jumping between 0% and 100% for over an hour, back and forth, climbing to 100% then dropping down to 35% or so all over again. Does anyone know how to fix this? I'm desperately trying to cut down this raster file so I can continue doing my work, since it has to re-draw the entire raster every time I add a point to a point shapefile. Arg! I'm using ArcGIS10 with SP4 installed. This occurred in ArcMap while running the Clip tool under Data Management Tools > Raster > Raster Processing > Clip. Thanks.


title: How can I get more verbose messages from a GP service?

text: When I run a tool I've published as a gp service, the only messages I get are these:               esriJobMessageTypeInformative: Submitted.      esriJobMessageTypeInformative: Executing...      esriJobMessageTypeInformative: Succeeded.       Is there some way I can get the same messages I see when I run the tool locally? I don't see anything obvious in the documentation. It is a custom tool written in arcobjects (implementing IGPFunction), so if there's some special way to call IGpMessages I'd consider that.


title: Create Vector Tiles for Polymaps

text: Ok, likely my bad for not making it to WhereCamp5280 and asking the FortiusOne geeks directly, but what's the process for creating _GeoJSON_ 'tiles' for use in Polymaps? The Population example over at polymaps.org notes that the data for the demo runs on Google AppEngine, and alludes to "rolling your own" vector tiles, but I can't seem to find any more information... Anyone got some insight they can share? Thanks!


title: How do I perform a proximity search with Postgis?

text: OK, so I downloaded the Geonames database (cities1000) and wrote a small Ruby program to import them into my table (`geo_cities`). I then added geography column called `geog`. So after all of that, I converted all of the lat/lon numbers into the geog column using:               update public.geo_cities set geog = st_GeogFromText('SRID=4326;POINT(' || longitude || ' ' || latitude || ')');      Things are looking good. Now what I want to do is find all cities within 100 miles of Prague. So I can get Prague like:               select * from geo_cities where asciiname = 'Prague' and countrycode = 'CZ';      I'm still learning gis and Postgres so could someone help me with the simple query? Thanks!


title: which map projection should I use to calculate distance in Poland

text: I'm using/learning geodjango + postgis, and I'm trying to calculate distance between 2 cities in Poland. First using WGS84/4326 I define 2 points, and then after transforming this points into spherical mercator 900913 i get distance which is about 40km too long. Which map projection should I use then? Where can I find proper srid's for different countries? Arek


title: arcpy layer issues

text: This should be a really simple process but I am having issues. I want to search a workspace for all files with a certain string and then add them to the current workspace with a selected symbology. Issues 1\. Error [[: Object: CreateObject cannot open map document]] I have to set the map doc to be the current workspace 2\. When run in the python window the symbology is not updated. 3\. It should only update the map frame/toc once everything is added Thanks,               import os, arcpy, arcpy.mapping     from arcpy import env          GDB = arcpy.GetParameterAsText(0)     Text = arcpy.GetParameterAsText(1)     symbologyLayer = arcpy.GetParameterAsText(2)     #OutputGDB = arcpy.GetParameterAsText(3)          arcpy.env.workspace=GDB     mxd=arcpy.mapping.MapDocument("CURRENT")     list=arcpy.ListFeatureClasses('*'+Text+'*')          for layer in list:         dataFrame = arcpy.mapping.ListDataFrames(mxd, "Layers")[0]         addLayer = arcpy.mapping.Layer(layer)         arcpy.mapping.AddLayer(dataFrame, addLayer, "BOTTOM")         arcpy.ApplySymbologyFromLayer_management (addLayer, symbologyLayer)          # Clean up the MapDocument object by deleting it     #del mxd          arcpy.RefreshActiveView()     arcpy.RefreshTOC()      ERROR DETAILS Messages > Executing: ProcessDatasets22 > P:\2012\Job_044_DM_Radio_Propogation\Working\FinalPropogation\cb_BEL07\working\bel07.gdb > _dow P:\2012\Job_044_DM_Radio_Propogation\Working\signal_strength_poly.lyr > Start Time: Mon Mar 19 14:11:27 2012 Running script ProcessDatasets22... : > Object: CreateObject cannot open map document Failed to execute > (ProcessDatasets22). Failed at Mon Mar 19 14:11:27 2012 (Elapsed Time: 0.00 > seconds)


title: How to query objects within 1000m epsg4326?

text: I have a postgis table having a column "geom" having a type geometry(MultiPolygon) and the projection epsg:4326 (degree decimal), i want to query all objects within 1000m and i use this query: SELECT ST_AsGeoJSON(geom) FROM table WHERE ST_DWithin(geom, 'POINT(long lat)', 1000.0); but i always ended up querying all objects.


title: What spatially-aware GUI administration softwares exist for PostgreSQL/PostGIS (free or paid)?

text: Most of my experience with PostgreSQL has been PGAdminIII and psql, but I'm interested to know what others use. Is there PostgreSQL equivalent to SQL Server Management Studio (SSMS) and its ability to display, in a limited way, results containing geometry. Or do most people just run a GIS client on the side for this purpose, as this question seems to indicate: GUI tools for PostgreSQL/PostGIS?? I looked down through the PostgreSQL software catalog, but I didn't find anything (other than broken links).


title: SQL expression to python expression

text: I have a SQL expression that I use in Make Feature Layer in Model Builder ("HOSPITAL" = '% HOSPITALNAME%'). The expression makes a selection in a OD Cost Matrix that a patient is patient only at one hospital and not all hospitals as an OD Cost Matrix would normally do. HOSPITAL is a column in orgins in add location and HOSPITALNAME is a column in the DESTINATION to add location. What I wonder is how this expression will be in a python script? I've tried something like this in python in the expression field in arcpy.MakeFeatureLayer "'HOSPITAL' = '" + HOSPITAL NAME + "'". But it does not works. I don't get any error message but it only runs OD Cost Matrix normally.


title: Interactive computation of spatial data on a web-map

text: All the GIS area is very new to me so I would like to ask about a guidance of my next steps. * * * ## This is the web-map I would like to make (if it is possible): Lets consider the following map: ![EGNOS Computation points](http://i.stack.imgur.com/7EBM0.png). Each cross represents spatial data in form of point. A number between 1-100 is associated to each point. Any other given point on the map is situated inside an imaginary square defined by the 4 closest spatial data points. For that given point, I would like to display a value witch is computed by linear interpolation of the 4 values of the square corners. So I would like to show for example: _lat: 43.9 deg; long: 18.1 deg, approx. value: 35.5%_ ![map with feature](http://i.stack.imgur.com/eIIqx.png) Is it possible to create such map? At which layer I should define the interpolation algorithm?


title: Uncaught SyntaxError: Unexpected token : when using ajax

text: I try to call json data from geoserver, but when i try to run it in chrome, it return this following error "Uncaught SyntaxError: Unexpected token :" This is my code,                   function initialize(){var urljson='http://localhost:1234/geoserver/Ven/ows?service=WFS&version=1.0.0&request=GetFeature&typeName=Ven:alamatpuskesmas&maxFeatures=50&outputFormat=json&format_options=callback:getJson';                       $.ajax({                         url :urljson,                         dataType: 'jsonp',                         jsonpCallback: 'getJson',                         success: handleJson                     });                  }             function handleJson(data)             {                 console.log(data);             }      i am using geoserver 2.4.2, whats wrong with my code?, when i try to paste the url in chrome/firefox, it return json. .


title: INSPIRE literature for the common man

text: As a GIS Coordinator I am constantly amazed at how difficult it is to comprehend the literature produced by the European Commission regarding the Inspire directive. If this is hard work for people who work full time in the profession how can anyone expect anyone working on the peripheries of GIS to understand it´s demands? With this in mind, does anyone know of any literature describing INSPIRE for the common man? I saw an "INSPIRE for dummies" pdf once which had great reviews but it was unfortunately in french. If anyone has something similar from a good source it would be very helpful. I am of the opinion that INSPIRE doesn´t have to be difficult. It has just become too complicated through bureaucracy.


title: How to convert EPSG:900913 or EPSG:4326 to EPSG:27700 using Proj4js (or to set displayProjection to EPSG:27700)

text: I'm using openlayers to create a vector feature (e.g. polygon, point) and i would like the geometry values of these vectors to be stored as EPSG:27700 coordinates (so when i send them off in WKT format they will be using EPSG:27700). I assume that this can be done by changing the displayProjection to EPSG:27700, but Proj4js doesn't have this projection, so how can i add it. If this is not the solution, how can i convert EPSG:900913 or EPSG:4326 (lat/long) to EPSG:27700? Update: I've tried using the EPSG:27700 projection (by adding to definition to Proj4js) but it's just reading as EPSG:900913 (unless they're extremely similar). Am i not converting properly? Code is here: http://pastebin.com/rV8NFgVh


title: ArcGIS: In the Legend, show legend items that are only displayed on the map

text: Using ArcGIS, is there a way for the legend, to only show the subitems that are displayed on the map?


title: Getting a "square" of lat/lon for database query

text: So I want to get the square formed by a single point of lat/lons, my plan was to take that point's latitude, divide by 360, and multiply by the circumference of earth. Giving us the distance from the equator. Now I will take 1/4 the circumference of earth, and subtract by the distance from the equator, giving us the height of triangle T. I can the find the base because I know the initial height, truncated height, and the original base. Now using the base and height I can use trig to find the minimum and maximum longitude, and the latitude is easy because lets are evenly spaced. Does this sound right? EDIT: Let me simplify this question. How can I find the distance between two points of longitude, at a single point of latitude?


title: Distance between two latitudes

text: I need to know the distance in feet between the coordinates N34-00.000 and N34-00.001 (Degrees, minutes and thousanths of minutes) I computed approx 6 feet but am not confidant in my math skills.


title: distance difference in Google Earth with exported GPX or KML compared to source shapefile

text: I've converted a shapefile (points) to KML with QGIS. Points in the shapefile are within 100m distance from each other. When I measure it in Google Earth, the distance changes to 80m. The same happens with the created GPX file. Why is this? My original shapefile was created in WGS84/UTM Zone 20S because I needed to work with the metric system. I'm working with data belonging to an area no larger than 6 square kilometers. As to avoid the problem I am having now, I reprojected it to WGS84 (EPSG 4326) by "saving as" and changing its SRC before I then created the KML and GPX. Both KML and GPX show the same surface, the difference comes along with the shapefile. Sample coordinates are:                   id  lat long     POINT(-61.82795607536099425 -37.83722374457502724)       POINT(-61.82861595101069696 -37.83772772712714527)       POINT(-61.82724009148368083 -37.83767676392982793)       POINT(-61.8278967620074269 -37.83817748597036257)


title: Allowing user to download shape files via geoserver

text: We provide layers via **shape files** and/or **PostGIS PostGreSQL Database** in `geoserver-2.x.x` with **GCS Geometry**. `geoserver-2.x.x` itself uses this shape files to render maps. The shape files are located on the server machine where `geoserver-2.x.x` is installed. I was just thinking to allow user to download specific shape file(s) of specific layer(s). I want user to download the shape file via the website developed using `Open Layers 3` and `Geo Ext` java script libraries with `Spring 4`,`Jsp`, `Hibernate`. As per this link WFS Output for shape file **Case 1:** Layer which is from PostGIS **csae 2:** Layer from shape file (.shp) (located somewhere in D:\shape\Layer.shp on server)   1. In both above cases will same code work for allowing user to download shape files?   2. If the files may too large or the data table used in `PostGIS` is too wide, still it will allow user to download the shape file?    3. For case 2 will it create a shape file by itself?   4. Because my layers are mixing of above two cases. Would I have same logic?    5. I want to understand how it works. Because I want to give bulk download also for group layers.   6. Is this possible? I also want to give the same feature in one of my project in world wind. Any help is invited.


title: get latitude & longitude as per address given for leaflet

text: In my project I am using leaflet to display points. I am getting address from my database & want to show point for that particular address. I dident found anything which will give me latitude & longitude as per address provided. Can any one help me on this.


title: Grouping of layer with GeoServer and OpenLayers

text: I need help to display layers group/non group 1\. with checkbox to select layer or group 2\. icon/symbol for each layer. I saw example codes from OpenLayers where it's adding the layers statically. I would like to populate it dynamically so in future if a layer is added/remove it will be displayed without changing any code. Sample layout of what I want to achieve is posted below, please note that checkboxes are displayed as dots. Groping of Layer   * Layer    *     * School      *     * Road      *       * Metal Road        * Unmetal Road      * Parecls      * Water Feature      * Cities      * Disulict      * India


title: Moran's I z-score in ArcPy

text: I have a ArcPy script to calculate Moran's I several times (50 to be precise). However, I am not interested in the Moran's I value, what I need is the z-score returned by the tool. Even though the Geoprocessing results window dispalays the z-score, I am not able to get the score through the script. I have been through the documentation on ArcGIS resource center and the example given by them is:               arcpy.SpatialAutocorrelation_stats("olsResults.shp", "Residual","NO_REPORT", "GET_SPATIAL_WEIGHTS_FROM_FILE","EUCLIDEAN DISTANCE", "NONE", "#","euclidean6Neighs.swm")      The result returned by this script is the Moran's I index value. Is there any way I can get the z-score besides looking at the Geoprocessing results window or the Reports?


title: Updating Attributes using Pyshp

text: Is it possible to update an attribute field in a shapefile using pyshp ? I went through their documentation and all they show is how to create new attribute rather than editing or updating an already existing one. If updating is not possible, is there any other recommended Python library for doing the same ?


title: How to do advanced Python labeling in Arcmap?

text: Does anyone have a good reference for advanced/complicated labeling expressions in Arcmap? I am attempting to write a Python label expression that would involve multiple label definitions. Here is what I am trying to do: (Using ArcMap 10.1) I have a point shapefile representing Monitoring Wells. For each Well I have taken concentration readings for various chemicals. What I need to do is label all of these concentrations for each well and color code them based on whether or not they exceed a certain limit. This is what I have:                    def FindLabel ([cis12DCE],[cis12DCE_Ex] ):             if int([cis12DCE_Ex]) == 2:              return "<clr red = '255'>"+"cis-1,2-Dichloroethylene  "+[cis12DCE]+"</clr>"            else:                 return "cis-1,2-Dichloroethylene  "+[cis12DCE]      My question is how do I repeat this so I get a label that looks somthing like this: ![enter image description here](http://i.stack.imgur.com/QmF2O.jpg)


title: How do I exclude raster layers from legend? (Silverlight, ArcGIS Server)

text: I'm trying to adjust my Print Widget so that raster layers do not show up in the legend. My thinking was that I could test for a layer type and not add the item to the legend if it is a raster. Something like:               Foreach (layer in List<Layers>)         If layer is a raster:             do nothing;         Else:              Add it to the legend;      I know there's a LayerInfo thing out there but I'm not really sure how to use it. Incidentally, I'm new to C#/Silverlight/ArcGIS for Server so this simple task is a bit challenging for me. I did notice there was a GetType() method(?) that I might be able to use but I'm not sure if how if that is indeed the answer. The foreach legend item code already exists and such -- I'm looking to wrap some conditions around the legendItems.Add bit. I found the following resources, but not sure how to use them if they are helpful... http://help.arcgis.com/en/webapi/sil...o_members.html \--> the "Type" property, perhaps? Alternatively, if I'm going about this in a difficult/wrong way, feel free to point that out! I'm a new intern with no prior C#/ArcGIS Server/Silverlight experience so it's all a bit overwhelming. Thanks in advance! Thread on arcgis forums: http://forums.arcgis.com/threads/83788-How-do-I- retrieve-layer-type additional thoughts: \-- I'll somehow have to create a list of all the layers, too >_<.


title: Is there a symbols level option/solution for SLD?

text: Is it possible to have an sld blend the symbology in the lines together just as QGIS and ArcGIS allow using symbol levels? In QGIS I can use symbol levels on the line layer and have this appearance: ![enter image description here](http://i.stack.imgur.com/T8XGB.png) currently in sld it has this appearance: ![enter image description here](http://i.stack.imgur.com/aHC1e.png) I was hoping to have the appearance of the individual lines blending together as one. The intent is to apply the sld to a geoserver layer. Currently I am only using rules like the one listed below.                   <Rule>           <MinScaleDenominator>2500</MinScaleDenominator>           <MaxScaleDenominator>12500</MaxScaleDenominator>           <LineSymbolizer>             <Stroke>               <CssParameter name="stroke">#ffe14c</CssParameter>               <CssParameter name="stroke-width">3</CssParameter>               <CssParameter name="stroke-linejoin">round</CssParameter>               <CssParameter name="stroke-linecap">round</CssParameter>             </Stroke>           </LineSymbolizer>         </Rule>


title: Using GeoXml3 to parse KML Markers loaded in Fusion Table Layer

text: I'm trying to apply the same code as this exmaple to zoom the google map to a KML Markers loaded into Fusion Table Layer. The example code works for KML Polygon loaded into Fusion Table Layer but not for markers. Does anyone know how to apply that to KML markers? kml example:               <Point><coordinates>-59.028645,50.772916,0.0</coordinates></Point>      script:               query = "SELECT * FROM  1eiOWq49sR8DXdymYALNdItE0uf-f4mNRbyJOVOU";     query.send(zoom2query);          function zoom2query(query) {        // zoom and center map on query results        //set the query using the parameter        var queryText = encodeURIComponent(query);        var query = new google.visualization.Query('http://www.google.com/fusiontables/gvizdata?tq='  + queryText);                   //set the callback function         query.send(zoomTo);                }            function zoomTo(response) {         if (!response) {           alert('no response');           return;         }         if (response.isError()) {           alert('Error in query: ' + response.getMessage() + ' ' ++ response.getDetailedMessage());           return;         }                 FTresponse = response;           //for more information on the response object, see the documentation           //http://code.google.com/apis/visualization/documentation/reference.html#QueryResponse           numRows = response.getDataTable().getNumberOfRows();           numCols = response.getDataTable().getNumberOfColumns();           var geoXml = new geoXML3.parser();           var bounds = new google.maps.LatLngBounds();           for (var i = 0; i < numCols; i++){             if (FTresponse.getDataTable().getColumnLabel(i) == 'geometry') {              var ColIndex = i;             }           }           if (!ColIndex){             alert('Geometry column "geometry" not found.')           }           for (var i = 0; i < numRows; i++){             var bounds = new google.maps.LatLngBounds();             for(i = 0; i < numRows; i++) {                 var kml = FTresponse.getDataTable().getValue(i,ColIndex);                 geoXml.parseKmlString("<Placemark>"+kml+"</Placemark>");                 bounds.union(geoXml.docs[i].bounds);             }             // zoom to the bounds             map.fitBounds(bounds);           }


title: How to delete a point with known coordinates in python?

text: I was wondering if there is a piece of code that does the opposite that setGeometry does and removes a point if there is one at x and y coordinates?


title: How to transpose reclassified field values from one polygon shapefile to identical polygons of another shapefile?

text: I used vector output from modeling software as clip features for ecoregion input. I reclassified and dissolved the clipped layers into 3 subunits, and now would like to transpose those simplified values to the source ecoregion layer in areas of intersection, such that the same 3 field values from the clipped layers are mapped onto the source ecoregion layer. I need to do this for about 200 pairs of shapefiles, and it seems like there must be some at least partially automated way to approach this. ![Top Layer displayed 50% transparent by the field value I need to transpose to bottom layer](http://i.stack.imgur.com/QnUkt.jpg) In the image below, I've symbolized the 3 field values semi-transparently over the source ecoregion layer. You can see that a single value polygon in the top layer overlaps multiple polygons in the source layer, but the outer boundaries of the individual polygons of the source layer exactly match the simplified polygons of the top layer. ![Field values from top layer transposed to bottom layer](http://i.stack.imgur.com/G8sFk.jpg) In this image, I've reclassified the source ecoregion layer to match the field values of the top layer, such that the field values from any polygons from the top layer that intersect with polygons in the source layer are copied to the attribute table of the source layer. I've been doing this via manual selection and the attribute table field calculator, but this will take forever for 200 sets of vector data. Anyone know of either a canned tool in arcGIS that could help with this, or python coding approach that could help speed this up? Thanks in advance for your help! **EDIT in hopes of clarifying:** For polygons in the reclassified Layer A and Layer B (ecoregions) that intersect, I'd like to populate the subunit field in Layer B with the same value as the subunit field in layer A. This is the attribute table of the reclassified Layer A: ![Reclassified attribute table](http://i.stack.imgur.com/mIwEa.jpg) This is the attribute table of the source ecoregion Layer B after I've manually populated the subunit fields based on intersection with Layer A: ![enter image description here](http://i.stack.imgur.com/Mnrro.jpg)


title: Extracting NetCDF raster values to time enabled shapefiles - what's wrong with my python script?

text: I'm trying to make a python script tool that extracts the values from a NetCDF raster to each data point from a time enabled shapefile at the correct time. My NetCDF raster consists of monthly sea surface temperature (SST) averages and I also have a plastic debris point shapefile that has been time enabled with time step interval of one month. I have checked with Time Slider that both the raster and the shapefile values change through time correctly. Can someone please help me figure out why I keep getting an error with this script:               mxd = arcpy.mapping.MapDocument('SST.mxd')     df = arcpy.mapping.ListDataFrames(mxd)[0]     rasterLayerName = "SST"     rasterLayer = arcpy.mapping.ListLayers(mxd, rasterLayerName, df)[0]     plasticLayerName = "Plastic"     plasticLayer = arcpy.mapping.ListLayers(mxd, plasticLayerName, df)[0]     df.time.currentTime = df.time.startTime     while df.time.currentTime <= df.time.endTime:        uniqueName = "Plastic_" + df.time.currentTime.strftime("%Y_%m")        outPointFeatures = "U:/gis/" + uniqueName        arcpy.sa.ExtractValuesToPoints(plasticLayer, rasterLayer,              outPointFeatures, "NONE", "VALUE_ONLY")     df.time.currentTime += df.time.timeStepInterval      Here is the error:               Runtime error       Traceback (most recent call last):        File "<string>", line 6, in <module>        File "c:\program files (x86)\arcgis\desktop10.1\arcpy\arcpy\sa\Functions.py", line 1311, in ExtractValuesToPoints          add_attributes)        File "c:\program files (x86)\arcgis\desktop10.1\arcpy\arcpy\sa\Utils.py", line 47, in swapper          result = wrapper(*args, **kwargs)        File "c:\program files (x86)\arcgis\desktop10.1\arcpy\arcpy\sa\Functions.py", line 1304, in wrapper          add_attributes)        File "c:\program files (x86)\arcgis\desktop10.1\arcpy\arcpy\geoprocessing\_base.py", line 498, in <lambda>          return lambda *args: val(*gp_fixargs(args, True))      ExecuteError: ERROR 001143: Background server threw an exception.      I'm using ArcGIS Desktop 10.1 and Python 2.7 on Windows 7.


title: Setting units in python using ArcGIS 10

text: I want to calculate the area of each polygon in a shapefile. How do I control what the units of measurement are?               def getArea(shapefile):         """Print area of polygons"""              Rows = ARCPY.SearchCursor( shapefile )               for row in Rows:                  tmp_area = row.area # this assumes a field name 'area' exists             print str(tmp_area)       I looked to see if there is any setting in the `arcpy.env.` module but could not find one. If I do set the units in my python script, does that mean that all aspects of the shapefile remain unchanged? For example, if the shapefile had units of feet, in my script I change the units to KM, will the shapefile still have units of feet afterwards?


title: Copy features between file geodatabase and ArcSde by using ArcObjects

text: I want to transfer all contents of a table from file geodatabase (*.gdb) to a table in ArcSde. In fact the data is meant to be replaced, so I need also to clear the table in ArcSde. Both tables have same structure. The target table is used in a map service. I started with ConvertFeatureClass and the effects were pretty good - it was working, but each time it was creating a new table. Every try to append the data to already existing table ended with errors. Here's my code - http://pastebin.com/RWhMTe4p Since this was not doing what I wanted, I tried something different. This time, I used IGeoDBDataTransfer interface. My efforts ended up with an error being thrown on the final Transfer method. That's my code - http://pastebin.com/MG3KUHdS The COMException is: "Error -2147216061: Exception from HRESULT: 0x80041543 Can someone please explain how to do this task?


title: How to draw some ellipses that can cover some categorized points automatically in QGIS?

text: Instead of using Rectangles Ovals Digitizing plugin to draw some ellipses manually, is there another way to draw them automatically that can cover some categorized points in QGIS?


title: import tile into PostGIS through java dll

text: i have used this command to import tiles from directory to postgis data base but it show error in command line please help me if any one know ?. thank advance to share knowledge C:>java -jar "C:\Program Files (x86)\GeoServer 2.4-RC2\webapps\geoserver\WEB- INF\lib\gt-imagemosaic-jdbc-10-RC2.jar" import -config "C:\Program Files (x86)\GeoServer 2. 4-RC2\data_dir\coverages\startimage\osm.postgis.xml" -spatialTNPrefix tileosm -tileTNPrefix tileosm -dir "C:\Program Files (x86)\GeoServer 2.4-RC2\data_dir\coverages\star timage\tiles" -ext png when i run the above command show this error java.lang.ClassNotFoundException: org.postgresql.Driver at java.net.URLClassLoader$1.run(Unknown Source) at java.net.URLClassLoader$1.run(Unknown Source) at java.security.AccessController.doPrivileged(Native Method) at java.net.URLClassLoader.findClass(Unknown Source) at java.lang.ClassLoader.loadClass(Unknown Source) at sun.misc.Launcher$AppClassLoader.loadClass(Unknown Source) at java.lang.ClassLoader.loadClass(Unknown Source) at java.lang.Class.forName0(Native Method) at java.lang.Class.forName(Unknown Source) at org.geotools.gce.imagemosaic.jdbc.Import.start(Import.java:366) at org.geotools.gce.imagemosaic.jdbc.Toolbox.main(Toolbox.java:46)


title: Add GPS tracks directly to a PostGIS table

text: Disclaimer: I am a total GIS newb and I have a use case I have not seen elsewhere. Background: I drive country roads looking for old building, cemeteries, school houses, etc and I record the locations and document them for a website. I want to use QGIS to keep track of the places in their categories (school, church, cemetery, etc) and also track the roads I have driven as a way to keep duplication of work to a minimum. I have successfully (I think) setup PostGIS, imported all my items (places + tracks from GPS) from a KML and setup a bluetooth GPS with a small netbook. I would like to use the built in GPS track maker to record where I have driven and append these lines to a PostGIS table which contains a list of other lines. I have seen examples of exporting the data to a GPX or SHP file but according to this post ( How to add features to a PostGIS layer using Quantum GIS? ) I should be able to just the line to PostGIS. According to the documentation I cannot add a GPS Track directly to the PostGIS DB, I must create a new layer and record it there then export / import the data to get it into the DB. I have had no luck with this up to this point. I get errors that there are missing columns, once those are added I get other errors about constraints. How do I accomplish this? I feel like I'm very close.


title: QGIS 1.8 map composition problem

text: I have recently installed QGIS 1.8 when I make a map with the adequate layer I switch to the map composition menu and try to insert a map in the white page. I draw the rectangle but nothing appears although I tried the 3 different option I also ticked on the show on map for all the layers (show in overview) can you tell me if this is an installation problem (I have windows XP professional and downloaded the french version) or a trivial problem I am not being able to resolve?


title: Access PostGIS Raster layer through mapserver

text: I am trying to use a cgi script calling mapserver to display a raster layer that is stored in a PostGIS database. However, only a blank white image is returned and I don't know where the problem is. The request seems to go through, but here is an error in the postgres log that I don't understand and the result is a white image. My questions are:   1. What does the error from the postgres server log mean (see below)?   2. How can I make sure the white image is not due to a misspecification of the map extent?   3. Is there something obvious that I am doing wrong? * * * _Here is some more detailed info for this problem._ The mapserver request is: `http://someurl.com/mapserv.cgi?map=../../usr/maps/test.map&mode=map&layers=kahoo&` The error from the postgres server log:               ERROR:  column "o_column" does not exist at character 39     STATEMENT:  select o_table_name, overview_factor, o_column, o_table_schema from raster_overviews where r_table_schema = 'public' and r_table_name = 'kahoo' and r_column = 'rast'      From the mapserver log:               .632473 msDrawMap(): rendering using outputformat named png (AGG/PNG).     .632541 msDrawMap(): WMS/WFS set-up and query, 0.000s     .633883 msDrawRasterLayerLow(kahoo): entering.     .859561 msDrawMap(): Layer 0 (kahoo), 0.227s     .859588 msDrawMap(): Drawing Label Cache, 0.000s     .859594 msDrawMap() total time: 0.227s     .862701 msSaveImage(stdout) total time: 0.003s     .862751 mapserv request processing time (loadmap not incl.): 0.231s     .862759 msFreeMap(): freeing map at 0x266bc60.      The mapfile used is:               MAP      CONFIG "MS_ERRORFILE" "../../tmp/ms_error.txt"      DEBUG 5           NAME "test"      SIZE 218 152      EXTENT 739335 2268075 7761335 2284075      UNITS METERS      STATUS ON      PROJECTION       "init=epsg:26704"      END           WEB       METADATA         "wms_title"           "WMS Demo Server"         "wms_onlineresource"  "http://someurl.com/mapserv.cgi?map=../../usr/maps/test.map&"         "wms_enable_reqest"   "*"       END      END           LAYER         NAME kahoo         TYPE RASTER         STATUS ON         DATA "PG:host=*** port=*** dbname='***' user='***' password='***' schema='public' table='kahoo' mode='2'"          PROCESSING "NODATA=0"         PROJECTION            "init=epsg:26704"         END      END     END      Here are the versions I am using: mapserver-6.0.3, postgresql-9.1.6, postgis-2.0.1


title: How to lock the legend in print composer?

text: I locked in the print composition (Druckzusammenstellung) the legend (attachment1). Then I saved the print composition and opened it (attachment2). But the locked legend is not locked? The changes are gone: polygon b is still part of the legend. What have I to do to lock the legend? Thank's for any advice. phil attachment1: ![enter image description here](http://i.stack.imgur.com/TIBnD.png) attachment2: ![enter image description here](http://i.stack.imgur.com/x6mQ6.png)


title: Maps API that look "fancy"

text: We have been using Google Maps V3 during the development of a web service. We've been struggling with getting the maps integrate good with the design, since the maps look kind of grey. Now when i look at Apple Maps, they look so much nicer (in standard mode that is). The colors are brighter, and the maps doesn't look "grey and boring" like Google Maps do. My questions are, 1\. Is there any way to get Google Maps to look like Apple Maps do? 2\. If not, is there any (preferably free) maps we could use that looks something like Apple Maps?


title: OpenLayers layer text popup override code not work

text: I have some problem with the popup code in OpenLayers, using Text layer (yes i know there already GML layer out there with text format, but i need to modify existing code), the problem was popup not removed when i open new popup (by clicking another marker)               OpenLayers.Feature.prototype.popupClass = OpenLayers.Popup.FramedCloud;                  OpenLayers.Layer.Text.prototype.markerClick = function(evt) {                 var sameMarkerClicked = (this == this.layer.selectedFeature);                 if (this.popup) {                     this.popup.toggle();                  } else {                     this.layer.map.addPopup(this.createPopup(true));                 }                 this.layer.selectedFeature = (!sameMarkerClicked) ? this :     null;                 OpenLayers.Event.stop(evt);                     }       code by Christopher Schmidt : http://www.mail- archive.com/users@openlayers.org/msg03869.html


title: How to get the First and Last Point coordinates of a Polyline with Field Calculator in ArcGIS10

text: I'm trying to get the First Point coordinates (x,y) and Final Point coordinates (x2,y2) from a polyline in ArcGIS 10 with the Field calculator. I've added 4 (double) fields like StartX,StartY,EndX and EndY,I tried with:               "!Shape!.positionAlongLine(0.0,True).firstPoint.X"     "!Shape!.positionAlongLine(0.0,True).firstPoint.Y"     "!Shape!.positionAlongLine(1.0,True).firstPoint.X"     "!Shape!.positionAlongLine(1.0,True).firstPoint.Y"      But no luck. **The error: "The row contains a bad value".** What am I doing wrong? Can anybody help me? \--- E D I T ---- Thanks to @radouxju for your answer because it was what I was looking for!!, **but NOW I don't know What's happening with my ArcGIS10.** Let's see: I could compute all these statements (4), to every field:   * "!Shape!.lastPoint.X"   * "!Shape!.lastPoint.Y"... But, after that I checked for errors and I found this! That there are times that the ArcGIS field calculator doesn't compute well. In the example, the yellow one, it's 16 meter long, but its coordinates are the same, but the the next polyline is well computed... What's happening? ![Error with ArcGIS10](http://i.stack.imgur.com/ri9G1.png) Thanks!


title: How to select features based upon multiple values in graphical modeler?

text: I'm trying to create a basic model using QGIS, the QGIS graphical modeler and the sextante toolbar. In part of the model, I need to select features based upon multiple attribute values (for example - A = 'value' OR A = 'another value'). The select by attribute algorithm provided by QGIS only allows you to enter a single value at a time - is there a simple tool/method to enter an expression to select multiple values from a dataset, and which can be incorporated in the graphical modeler? Thanks in advance.


title: Open QGIS project on different computers

text: I have multiple QGIS projects that need to go to different computers. I put all the files on a USB and loaded it under "Documents" on the other computers like it is saved on my own computer. However, when opening the projects on another computer the "Handle bad layers" window pops up. There are hundreds of files for each project. I am wondering if there is a way to get around this without having to click on each file individually and finding the new location? ![enter image description here](http://i.stack.imgur.com/MXT9A.png)


title: Disaggregating multiring polygon for editing, then recombining

text: **BUMP?** We have many rows of multi-ring polygons, displayed on a Flex front end. We need to be able to select a polygon (a row), split the rings for editing, edit them, then combine into one polygon again. I thought we could mulipart to single part, but that requires a whole FC and copies it to a new FC; whereas we want to be able to split, edit then recombine. I have built a spatial toolset class in Python and am now building the tools for this. But, I must admit, I am struggling to conceive how to do it; it being Friday doesn't help too ;). I also have another approach, which would be to select the id of the ring associated with the polygon; i.e., the polygon will have 1..n rings. I don't know the Flex API well, but I am sure you must be able to select the ring id of the polygon and removeRing()? Any ideas people? TIA. **EDIT** Anyone got any good resources online for Flex development?


title: Why is the 'straight line' path across continent so curved?

text: This is the result of mapping the straight line path from a point in US to Poland using Distance Measurement Tool. Also, planes from Asia to US would travel almost over North Pole. !['straight line' distance from Alberta to Poland](http://i.stack.imgur.com/v7RdT.png) Why is the path so curved? I agree that this is a flat representation of a sphere, so I do expect some arc, but I don't think earth has this much curvature. What am I missing here?


title: Is it possible to use GRASS GIS in Python stand alone scripts?

text: I am trying to find a least cost path script for Python. I came across r.drain from GRASS GIS. Is it possible to use GRASS GIS in a Python stand alone script? If so where can I download the bindings?


title: Can I alter font size using label expression?

text: I have two label classes displaying dynamically over my parcel polygons. The owner info is the major label with larger font size and bellow it I am displaying some property type info in smaller font. The problem is that the 2 label classes just will not display the way I want them to; Major label at the top, minor label just bellow it. I've used maplex and the regular labeling engine and played with the options but it just does not do what I want and the labels jump all over the parcels. The only way that works is to export the labels into a feature linked annotation size then adjust the label positions manually. This is OK but i need to change the scale often on this map. **Can I write the label expression so that I can have 2 different fonts and sizes within one label?** (Something similar to using text tags with text graphics)