title: Is it possible to setup ArcGIS Online edits to post to ArcGIS Server?

text: I've got ArcGIS Online setup and have all relevant layers on ArcGIS Server shared out as map services and added to ArcGIS Online. But what I'd like to find out is if they can then be edited in Online and those changes would reflect in Server.


title: Mississippi Township Range Section

text: Does anyone know where I can get a full complete T/R/S for Amita County in MS? Geocomm shp files looks incomplete. The counties in Mississippi are not helpful. I ended up finding the TRS from this website http://www.maris.state.ms.us/HTM/DownloadData/Statewide.html


title: find and replace, field calculator

text: Can someone show me a very simple find and replace in field calculator for just one or two characters in a field based hyperlink? Thanks Dustin


title: How can I use Oracle views with Geoserver and WFS-T?

text: What goes on when Geoserver handles updates to my table using WFS-T? I have a table               AREA (     ID NUMBER(10) NOT NULL PRIMARY KEY,     SOME_OTHER_ID NUMBER(10) NOT NULL FOREIGN KEY REFERENCES FOO(ID),     GEOMETRY MDSYS.SDO_GEOMETRY NOT NULL     );      and a sequence SEQ_AREA for the primary key values. When I publish this table as a Geoserver layer, it seems that it automatically detects that ID is the primary key and does not show it in the layer's properties. When a new row is inserted, Geoserver does not use the index - how does it determine the next value? Does it do SELECT MAX(id) FROM AREA or something like that? Could I make it use the sequence? And - if I create a view and a layer from that - the view does not have primary key as such - how can I make Geoserver understand it so that it is updatable? When I create a view such as:               CREATE VIEW TESTVIEW AS       SELECT a.ID, foo.somefield, a.some_other_id, a.GEOMETRY        FROM FOO foo, AREA a       WHERE foo.ID = a.SOME_OTHER_ID;            insert into USER_SDO_GEOM_METADATA (table_name, column_name, diminfo, srid)     values ('TESTVIEW ',     'GEOMETRY',     SDO_DIM_ARRAY     (       SDO_DIM_ELEMENT('X',50199, 761274,0.5),       SDO_DIM_ELEMENT('Y',6582464, 7799839,0.5)     ),     3047);       The layer ends up being viewable OK, but cannot be updated / inserted. When updates are attempted through wfs-t I get: 2013-02-26 19:00:33,636 ERROR [org.geoserver.ows] - org.geoserver.wfs.WFSTransactionException: {http://www.opengeospatial.net/cite}TESTVIEW is read-only


title: Map Coordinates to Neighborhoods from MySQL Table

text: I have a MySQL database with a set of coordinates and a shape file (Zillow Neighborhood Database for Chicago) containing the boundaries for each neighborhood which I have converted to a spatial database. How do I map each set of coordinates to the neighborhood that it is in? UbuntuGIS is installed.


title: Deployment options of arcgis desktop add-in

text: I made an addin for arcmap and want to deploy it using setup file. By double clicking the setup file anyone can install it in his machine and through installation it will check all dependencies and license issues. After the installation, the addin can be accessed from arcmap. I used visual studio setup file to do the first portion of the setup, but I couldn't access the addin from arcmap. Is there any other procedure or any special step to do this?


title: Percent cover based equation for horizontal ground visibility?

text: I'm trying to put together a layer representing horizontal visibility at ground level. I'm thinking of using the EVC (vegetation cover) layer from the landfire project (see http://www.landfire.gov/vegetation.php). However, I know that, for instance, 70% cover does not correspond directly to 30% horizontal visibility...I imagine there is a threshold at which the cover saturates any visibility, with visibility increasing in a non-linear way as percent cover decreases - probably calculated over some neighborhood (100 meters, maybe?). I'm thinking of cooking something up myself along these lines. However, I'm wondering if there hasn't been work done in this realm, or some theory to exploit that anyone may be aware of. A lit/google search hasn't turned up much of interest.


title: Instalación masiva de QGIS - Automated installation of QGIS

text: Solicito lo siguiente requiere la línea de comandos que permitan realizar la instalación automatizada o desatendida para poder enviarla desde SCCM para la Instalación de QGIS Google Translate - Spanish to English I need command line syntax that allows for automated/unattended installation of QGIS via System Center Configuration Manager SCCM.


title: Speeding Up ArcGIS, redrawing sooooo slow!

text: I just updated ArcGIS to 10.1 and it is so slow. When I load my template mxds they can take 30 plus seconds. And if I have to change the view it's another 30 seconds. Does anyone have any ideas or methods to speed up the loading and redrawing purpose? Thanks


title: QGis Grid statistics for polygons error - Module caused access violation!

text: I'm trying to calculate statistics for a raster layer using an overlaying polygon layer, I'm using SAGA "Grid Statistics for polygons" function, but it returns an error, sadly.. So in the sextante log it says :               Parameters          Grids: No objects     Files: "Y:/.../clipper_grassexport_0to11"     loading: Y:/.../clipper_grassexport_0to11     Driver: Y:/.../clipper_grassexport_0to11/     Cells: x 44069, y 42309     Bands: 1     Transformation x' = 2635889.298915 + x * 100.001118 + y * 0.000000     Transformation y' = 5421705.687421 + x * 0.000000 + y * -100.000407          Module caused access violation!          GDAL: Import Raster: Module caused access violation!     error: executing module [GDAL: Import Raster]          ***          go...     Load grid: C:\Users\...\AppData\Local\Temp\sextante\1368631882.798.sgrd...     failed     error: Grid file could not be opened.     error: input file [C:\Users\...\AppData\Local\Temp\sextante\1368631882.798.sgrd]     error: empty input list [GRIDS]     Usage: Grid Statistics for Polygons -GRIDS <str> -POLYGONS <str> [-RESULT <str>] [-COUNT] [-MIN] [-MAX] [-RANGE] [-SUM] [-MEAN] [-VAR] [-STDDEV] [-QUANTILE <num>]     -GRIDS:<str>    Grids     Grid list (input)     -POLYGONS:<str> Polygons     Shapes (input)     -RESULT:<str>   Statistics     Shapes (optional output)     -COUNT          Number of Cells     Boolean     -MIN            Minimum     Boolean     -MAX            Maximum     Boolean     -RANGE          Range     Boolean     -SUM            Sum     Boolean     -MEAN           Mean     Boolean     -VAR            Variance     Boolean     -STDDEV         Standard Deviation     Boolean     -QUANTILE:<num> Quantiles     Integer     Value Range: 0.000000 - 50.000000          error: executing module [Grid Statistics for Polygons]      I have tried to load my raster and vector layers into SAGA directly, but the raster layer wouldn't load, and same error appears - "Module caused access violation!" Checked if SAGA is properly installed, and is 2.0.8 version (QGis 1.8.0) Raster is an GEOTIFF exported from GRASS Tools in QGis in Byte format (Values 0-255), CRS's also match, I have run out of ideas, what could possibly be wrong? It would be great to have some help! Thanks! * * * Update - June 17 2013 So to sum it all up - with a large dataset it still gives an error, but if I want to divide the raster in say 4 smaller parts, raster Clipper Tool does it for ages and freezes with an approximate filesize of 100-140 MB, that means I have found no solution yet, that does the job efficiently. Suggestions?


title: adding kml to google maps

text: I am trying to overlay the KML file over google maps and want to know where I can host the kml file which I will then be able to add to the google maps. I tried using google drive and dropbox but they do not allow the file to be read from these sites. Can somebody tell me where I can host it? **UPDATE** According to this google-groups discussion, I cannot overlay a KML file from my local system.


title: Spatialite, to speed performance, when should I create Views, Indexes or additional spatial tables?

text: I’ve got a **large Spatialite database** with a table of lines (roads, rail and water). I created a spatial index of the geometry to speed things up. This works amazingly well when I’m zoomed in close. However, when I’m zoomed out in QGIS 1.8, the display is slower than I’d hope. And so I created a separate Table-View of just the major roads, intending to speed things up. That didn’t work—I suppose because the display is still plowing through all the roads just to find the major roads. Is it justified for me to create another spatial table of “just the major roads” to use when I’m zoomed out. Or, _is there some way to index the queried column for major roads, or some other trick, and thus create a very speedy display of the major roads without creating a brand new spatial table?_ I’m trying to keep my Spatialite database as small and practical as possible, without filling it with a lot of redundant tables of data. I’m finding that **Spatialite is extremely** quick, but perhaps I’m expecting too much performance, though I’d still like to get as much as possible.


title: How to query if a set of polygons intersect another set of polygons

text: Is there a way to leverage .STIntersect to check if a set of polygons intersect another set of polygons? I have a MasterTable of polygons that represent all waterbodies in a florida County. These polygons are used in an inspection process and I need to query the data to see how many inspections were completed in a given region of polygons. Ive been working on this all day and have got no where. HELP!


title: Network Analyst provides different values for Total length and Shape length

text: I'm using the Network Analyst extension in ArcMap 10.2 to calculate the shortest routes between ~50000 pairs of stops. As the impedance is set to "Length (Meters)" the resulting Routes Featureclass has a field for "Total_Length" and "Shape_Length". What is confusing me, is that the values in these two fields are not the same, they differ by +15m to -0.5m. At first I assumed that this had something to do with the distance of the stops to the network, but looking more closely at some of the routes, the distance to the network doesn't match the difference in length. I've also re-run my analysis with a 0m offset of the Stops from the network, but the length differences stayed the same. Has anybody else come across this before or does anybody have an idea of why this is happening?


title: building a water mask

text: I'm working with multi-temporal datasets from modis and had to do various vegetation indices. I ran into a problem now because there is a significant amount of water too, which skews the resulting map. How can I build a mask to just blend out the water bodies?


title: Mapserver file to localhost

text: I have installed OSGeo live 7 - Xubuntu I made a nice map using qgis and I have exported the project using RT mapserver exporter. So I created a *.map file in home/user/documents/ How can I load this *.map file to localhost?


title: Principal Component Analysis in QGIS

text: How do I go about to perform a band ratio and Principal Component analysis in QGIS? I need to derive principal components from ratios of landsat TM bands .


title: How to achieve point to point routing?

text: Currently, the pgrouting 1. 5 version has a shortest path function which actually traces a path from nearest vertex of source to nearest vertex of destination but. As, Dijkstra's algorithm works on the principle of vertex to vertex routing. But, if the location (point geometry) lies at somewhere in the mid of the road geometry and I want the path to be traced until that location exactly, is there a special function or way to do the same in pgrouting? In short, currently node to node routing happens. Whereas, I am looking for point to point routing. Please refer: http://www.mapgears.com/downloads/foss4g2013/foss4g2013-pgrouting-for- dummies.pdf As per this document as well as issues raised in the mailing list, this script helps in achieving that. But, I figured out that there are multiple issues in the routing_core_smart.sql file. Some tweaking needs to be done to make it work. And that is what is unknown.


title: How to find the projection of an image?

text: Instead of georeferencing the image and deforming it, I'm trying to find the projection used by the author of the image. ![enter image description here](http://i.stack.imgur.com/XKLDI.png) I've got tens of these distribution maps to digitize. The idea is to avoid image distorsions and to find only the parameters of pixel dimension and top- left corner coordinates to move maps in the correct position. To do that I need the original coordinate system. I used ArcGis 10.2.1 to compare and customize different coordinate systems, but no one matched and the variable parameters for every projection were few. I used Flex Projector, but it's not possible to import and overlap an image for a comparison, moreover the software exports many parameters I cannot use in a desktop GIS. Using EPSG codes in Qgis would take me months to find the proper one. Any ideas how to solve it? or which GCS could fit better?


title: Postgresql Connection issue "unexpected EOF on client connection"

text: I'm using Tomcat, Apache (for Load Balencing),OpenLayers - 2.11,GeoServer 2.13, ExtJS 3.4, PostgreSQL 8.4,PostGIS 1.5 and Java for developing GIS based web applications. I have 2 Server. Server1 contains   1. Apache(for Load balancing. This will divert all GIS links to ServerIP1:8081 and ServerIP2:8080)   2. Tomcat(ServerIP1:8080)   3. GIS Application   4. Geoserver Files   5. GIS Database(PostgreSQL 8.4,PostGIS 1.5) Server2 contains   1. Tomcat(ServerIP2:8081)   2. GIS Application   3. Geoserver Files All Map link will directly go to Apache from there after the load balancing process that will divert into Server1 or Server2. ## LoadBalancer code in Apache               <Proxy balancer://clusterABCD>         BalancerMember http://ServerIP1:8080         BalancerMember http://ServerIP2:8081              Order allow,deny         Allow from all     </Proxy>     ProxyPass / balancer://clusterABCD/      ServerIP1 and ServerIP2 contains same GIS applications and the same Geoserver files. Both Server accessing datas from ServerIP1:5433 Postgres Database Files. While running `select count(client_addr),client_addr from pg_stat_activity group by client_addr;` this query i'm getting from ServerIP1 Postgres counts are increasing continuously(count upto 50 or 100 or 200 like) and after some time intervals(not a fixed time) that will come down to 1 then 2 then starting to increase the count. If i stop ServerIP1:8080 tomcat then there is no problem. Then i have checked which query is taking this much of connections in PostgreSQL But, i'm getting results only **IDLE**. And getting only               2013-10-24 18:23:43 IST ERROR:  relation "public.gt_pk_metadata" does not exist at character 15          2013-10-25 06:15:31 IST LOG:  unexpected EOF on client connection      these errors in postgresql logs. Geoserver Store: ![enter image description here](http://i.stack.imgur.com/91rxj.gif) If anyone faced a problem like this please guide me to solve this problem.


title: Releasing GeoServer WFS-T Feature Lock

text: For some reason I'm having trouble finding any documentation about releasing feature locks in a GeoServer WFS (or any WFS for that matter). I can easily request a FeatureLock with a request like this `http://localhost:8080/geoserver/wfs?service=WFS&version=1.1.0&typeName=namespace:layer_name&maxFeatures=1&outputFormat=json&request=GetFeatureWithLock&cql_filter=Id%3D0` I have read that an update will release the feature lock once completed, but what if my user cancels? Should I just update the feature with its original data, or wait for the lock to time-out? Looking through the WFS GetCapabilities response I thought I could make a LockFeature request and pass a ReleaseAction parameter, but this didn't work out. The URL was `http://localhost:8080/geoserver/wfs?service=WFS&version=1.1.0&typeName=namespace:layer_name&maxFeatures=1&outputFormat=json&request=LockFeature&releaseAction=ALL&cql_filter=Id%3D0` but the result was 'A LockFeature request must contain at least one LOCK element'. The GetFeatureWithLock response didn't include any new lock identifier, and there are no cookies in the response headers, so I don't really understand how I'm supposed to release the lock I just created. Any thoughts on this would be very much appreciated!


title: can grass modules work separately

text: Grass is very useful tool for pre processing of imagery.i want to use i.atcorr module separately.will it be possible. i have compiled using .cpp files(AerosolConcentration.cpp, AerosolModel.cpp,Altitude.cpp,AtmosModel.cpp,etc) given but its asking grass functions to be declared.i dont want to use entire grass, only this i.atcorr is required. how can i build that single module?please suggest me thank you


title: Testing Tools For Simulating User Interaction - Especially Map Panning

text: Can anyone suggest any tools that can simulate users panning around in a map? I have used tools such as SoapUI to do load testing, but never anything to simulate users. In case it matters, the app I will be testing is based on ESRI Javascript API.


title: QGIS not reading colormap from hand-edited .qgs file

text: I have a set of rasters which store solar radiation incidence over the course of a year. In QGIS, I want to map the same extent, and with the same reference features, varying only the solar raster within the set. My idea was to create a master project file, symbolize and lay out the map in print composer using the first raster of the set, save the project, and then edit the .qgs file as XML. For each map in the series, I did a find and replace for all instances of the raster name, then saved as a new .qgs file. When I opened the projects again in QGIS, each raster showed in the proper file, but it had lost its symbology and reverted to greyscale. Each file contains this code:               <customColorRamp>       <colorRampType>INTERPOLATED</colorRampType>       <colorRampEntry red="255" blue="236" value="0.000000" green="247" label=""/>       <colorRampEntry red="254" blue="200" value="5000.000000" green="232" label=""/>       <colorRampEntry red="253" blue="158" value="10000.000000" green="212" label=""/>       <colorRampEntry red="253" blue="132" value="15000.000000" green="187" label=""/>       <colorRampEntry red="252" blue="89" value="20000.000000" green="141" label=""/>       <colorRampEntry red="239" blue="72" value="25000.000000" green="101" label=""/>       <colorRampEntry red="215" blue="31" value="30000.000000" green="48" label=""/>       <colorRampEntry red="153" blue="0" value="35000.000000" green="0" label=""/>     </customColorRamp>      Furthermore, a diff check on the files shows that resetting the symbology in QGIS and saving changes only the timestamp property within the file, nothing else. So if the colormap isn't being saved in the .qgs, where is it being saved?


title: How to create a worldfile by bounding box?

text: I need to create a worldfile from Bigmap Application. I've successfully to grab a pack of PNG tiles and the pearl app mosaicking all in one. The BigMap site shows the raster bounding box and util metadata: Map is 12x12 tiles (3072x3072 px) at zoom 18, aspect 1:1 -43.475647,-22.886297,-43.459167,-22.871114 (l,b,r,t) This is my local png map | http://goo.gl/oY8vA This is my pgw configuration: 3072----------X?   0   0   -3072---------Y?   -43.475647----left?   -22.871114----top?   How I create a PGW file in Quantum GIS or another app? Jorge Santos


title: Python Toolbox (.pyt) Tool Output Layer Name

text: How can I manipulate the TOC display name of my output parameters when using a Python toolbox (.pyt) Tool class. I tried adding a symbology layer, however this has no affect on the TOC name. Is it possible to reference the layer(s) added as output from the execute method of a tool. The expected results are added to map and TOC, but I can't figure out how to manipulate the TOC display name. I could loop through the layers in the TOC after exectution and check each layers datasource, but this seems problematic. I've searched the help and applied the extent of my Google Foo. In my tool class I define my parameters as follows:               def getParameterInfo(self):         extent = arcpy.Parameter(                 displayName='Analysis Extent', name='extent_features'                 , datatype='GPFeatureLayer', parameterType='Required'                 , direction='Input')         extent.filter.list = ['Polygon', ]              tilelayer = arcpy.Parameter(                 displayName='Analysis Tiles', name='tile_layer'                 , datatype='GPFeatureLayer', parameterType='Required'                 , direction='Output')         tilelayer.symbology = os.path.join(os.path.dirname(__file__), 'pyltk/analysis_tiles.lyr')      **Follow Up:** This normally would not be a problem as the output would reflect the name of the feature class. However, I've been evaluating the use of the new SQLite workspace for integration with external libraries, etc. ArcGIS prepends the internal SQLite database name 'main' to all my feature classes and tables. When ArcPy adds the results to the map it only adds the database name, and not the whole name. For instance if the output feature class is called "..\analysis_tiles" ArcGIS will refer to this as "..\main.analysis_tiles", but the resulting TOC name is just "main".


title: Unable to get processing extent to match desired raster output

text: I have a few thousand ortho mosaic tiles that I've resampled to a lower GSD and am trying to cut into larger tiles. I developed a script that has worked perfectly for thousands and thousands of tiles over the course of many projects. I am having a problem now where at the upper right corner of my current project area the tiles are not cutting into the proper image extent, but only for two of the five tiles I'm attempting to cut. I'm using ArcGIS 10.2 for the clipping process. The script utilizes arcpy.Clip_management. The images are TIFFs. I've played with dozens of environment settings, including processing extent, and I get the same thing every time. The input data has RGB values of 0,0,0 outside of the bad extent it is clipping to. The script itself clips to a shapefile defining the area of interest, but has always matched the input extent regardless of whether or not there was actually image there. The process clips to these bounds: ![Bad extent](http://i.stack.imgur.com/hUe0u.jpg) When I need it to clip to these bounds (and yes, I need/want all of those black pixels so image size will match all other tiles): ![Good extent](http://i.stack.imgur.com/Ul4vn.jpg)


title: NWS Grib Files scheduled download script

text: I am downloading grib2 files from the NWS NOMADS GFS service details here The files are provided on a 4x per day cycle and the date format is YYYYMMDDHH. I believe that the processing is actually done 4 times per day. I am a little confused about the naming convention and what time to actually go and get the files. An example url using the filter service is: http://nomads.ncep.noaa.gov/cgi- bin/filter_gfs.pl?file=gfs.t00z.pgrbf00.grib2&lev_surface=on&var_PRES=on&leftlon=0&rightlon=360&toplat=90&bottomlat=-90&dir=%%2Fgfs.YYYYMMDDHH My understanding is that he HH can be 00, 06, 12, 18.   * How does this file naming system work?   * Does any one have a script that downloads this data on a scheduled basis and figures out what the file names are?   * does this file have forecast weather available in advance?


title: How to get 25 points within 1 square kilometer in Arcmap

text: I tried to generate 25 sampling points within 1sqkm using distance point to point is 250m but i got only 15 points not 25 points as i expected because starting point was not at the corner. The tools I used was Creat fishnet.


title: openlayers Googlemap and custom wms

text: I have a custom wms server the transform for this wms is UTM and i cannot set googlemap layr for background for example my feature in wms is Iran country but googlemap background is Algeria this is my code :               var options = {                  //maxExtent: new OpenLayers.Bounds(35.7008418592904, 51.3777884858342, 35.7520557036613, 51.4272609453396),             //maxExtent:new OpenLayers.Bounds(506968.03125, 3936165,557431.4375, 3965082.5),             units: "meters",             format: "image/jpg",             projection: new OpenLayers.Projection("EPSG:2058"),             displayProjection: new OpenLayers.Projection("EPSG:900913"),             numZoomLevels: 50         };     function init() {         var capaGoogle = new OpenLayers.Layer.Google(                    "Google Satellite",                    { sphericalMercator: true,projection:"EPSG:900913",                    transparent: true                    }                );                   map = new OpenLayers.Map('MapDiv', options);         wms = new OpenLayers.Layer.WMS(                     'WMS Layer Title',                     '../wms.ashx?dis=t',                     { layers: '', version: '1.3.0', srs: "EPSG:2058", transparent: true },                     { singleTile: true, ratio: 1, displayOutsideMaxExtent: true }                 );          map.addLayers([wms, capaGoogle]);      var lon = 3943679.05500279;         var lat = 538414.111109505;              var zoom = 11.3;              center = new OpenLayers.LonLat(lat, lon);         alert(center);         map.setCenter(center, zoom);               }      i used transform like               center = new OpenLayers.LonLat(lat, lon).transform('EPSG:2058', 'EPSG:900913');      but that is not worked :(


title: It is possible to invoke PostGIS internal function within new PostgreSQL user-defined function in C?

text: I have fuction f() and function ST_Area() for example. Can I invoke ST_Area() in the body of my function f()? Function f() is written in C. The code might look like this:               #include "postgres.h"     #include "fmgr.h"          #ifdef PG_MODULE_MAGIC     PG_MODULE_MAGIC;     #endif          PG_FUNCTION_INFO_V1(f);          Datum f(PG_FUNCTION_ARGS)     {         PG_RETURN_FLOAT8(ST_Area(PG_GETARG_GEOMETRY(0)));     }      Thanks.


title: How to restore the table window in ArcMap 10.2 (stuck in maximized view)

text: I've somehow managed to maximize the table window in ArcMap 10.2 (after choosing > Open Attribute Table from the TOC). Whenever this window opens, it's maximized, and I can't see an option to restore it so it doesn't take up the full screen - the only option is to close it. (A possible complication is that I'm running this in Windows 7 via Bootcamp on a Mac).


title: Calculating watershed runoff with landuse and rainfall data?

text: What is the method for computing rasters of total rainfall of a given month, flow accumulation, and land use (runoff coefficient) or soil K factor? I am using ArcMap 10.


title: opengeo suite and pg_hba.conf

text: I need help with how to add access to my computer. When I edit the pg_hba.conf i still get the access denied error. I just want to give a specific ip access to all db and all users from that ip... no encryption. I am using opengeo suite 3.0.2 and pgsql 9.1.


title: Embedded spatial database for Java

text: Are there any embedded spatial databases for Java? I've come across H2 Spatial, but it was difficult to tell the status of that and if it can be used as an embedded spatial database. thanks, Jeff


title: insert maximal number of circles into polygons

text: > **Possible Duplicate:**   >  algorithm to place maximum number of points within constrained area at a > minimum spacing I'm trying to fit as many 30ft diameter circles into polygons of certain land cover classes without any overlap, neither with one another nor with adjacent classes. I'm trying to fit the maximum number possible. I'm having two problems, I can't figure out how to drop circles and I can't figure out how to maximize the number of them. Any thoughts?


title: Placing label designs in the legend - Is it possible in ArcGIS 9.3.1?

text: I'm working in ArcGIS 9.3.1 and apart from "regular" slayers with symbology, I have a few layers which I do not show their symbol on the map, but only their label (for example: cities layer in which I've made the point symbols transparent and their labels printed in blue). Question: Is there any way I can show the labels of the layers on the legend instead of their symbols (which is the default)??


title: ArcGIS: Scalable labels in x direction only

text: I'm trying to label a road layer that is named in 2 different ways:   * XYZ    * String between 5 to 20 characters Is there a way using ArcGIS 9.3 to produce a scalable ~~layer~~ **label** only in the x direction? Specifically I'm looking for a rectangle with semi circles on either end.


title: How is Mapquest Nominatim Service's "routeWidth" search param implemented?

text: I was looking at MapQuest's Nominatim service, and it seems like they have a non-standard search parameter (doesn't seem to exist in standard Nominatim's list of params) that is the routeWidth. Any guess how this is implemented in their service?


title: ArcObjects (vba) update selected feature/record

text: I would like to edit a selected feature/record via ArcMap vba, but I'm not having much luck with this code sample (see code below). I'm not getting an error, but the record selected is not updating.               Private Sub Update_Click()          Dim pMxDoc As IMxDocument          Dim pLayer As ILayer          Dim pInFtrLyr As IFeatureLayer          Dim pFtrSel As IFeatureSelection          Dim pFeatureClass As IFeatureClass          Dim pFields As IFields          'Dim ii As Integer          Dim pPoint As IPoint          'Dim pSegColl As ISegmentCollection          Dim pEnumLayer As IEnumLayer          Dim pMap As IMap          Set pMxDoc = ThisDocument          Set pMap = pMxDoc.FocusMap          Set pEnumLayer = pMap.Layers               Set pLayer = pEnumLayer.Next          Do Until pLayer Is Nothing          If pLayer.Name = "Copy of Wells_test" Then          Set pInFtrLyr = pLayer          End If          Set pLayer = pEnumLayer.Next          Loop               Set pFtrSel = pInFtrLyr               If pFtrSel.SelectionSet.Count = 0 Then          MsgBox "Please select one Well record", vbOKOnly          Exit Sub          ElseIf pFtrSel.SelectionSet.Count > 1 Then          MsgBox "Please select one Well record", vbOKOnly          Exit Sub          Else          If pFtrSel.SelectionSet.Count = 1 Then          Set pPoint = pInFtrLyr.FeatureClass.GetFeature(pFtrSel.Selectio nSet.IDs.Next).Shape          End If          End If               Set pFeatureClass = pInFtrLyr.FeatureClass               Dim pFSelection As IFeatureSelection          Set pFSelection = pInFtrLyr               Dim pFClass As IFeatureClass          Set pFClass = pInFtrLyr.FeatureClass               Dim pFCursor As IFeatureCursor          Set pFCursor = pFClass.Update(Nothing, True)               Dim pID As New esriSystem.UID          pID.value = "esriEditor.Editor"          Dim pEditor As IEditor          Set pEditor = Application.FindExtensionByCLSID(pID)          Dim pFeature As esriGeoDatabase.IFeature          Dim pEnumFeature As esriGeoDatabase.IEnumFeature          Set pEnumFeature = pEditor.EditSelection          Set pFeature = pEnumFeature.Next          Do Until pFeature Is Nothing          pFeature.value(pFeature.Fields.FindField("Type")) = "update"          Set pFeature = pEnumFeature.Next          Loop               End Sub      Thank you for your help.


title: Designing Better Maps

text: Using this chart, enhanced my presentations and presentation skills: ![Extreme Presentation Method - Ten Steps to Persuasive Presentation](http://i.stack.imgur.com/Zmp6L.jpg) found on Extreme Presentation Method. How would you modify this chart to be directly related to map design?


title: Style Manager in ArcMap

text: On my desktop PC, I am using Ministry of Transportation Ontario road shields to label major highways on my maps. I often need to work with a laptop when I'm visiting a client at their location, so I copied the appropriate style files to the same folder on my laptop. When I try to apply the roads shield markers on the laptop, it gives me weird symbols and text instead of the appropriate shield. I had no problem doing this on my desktop, so I must be missing something. On the left (in the image below), I have an example of the road shield marker style that displays on maps on my PC. On the right (column highlighted in yellow) are the marker symbol options I get in the style manager menu when I try to access these shields when using the laptop. I assume this is a simple fix and I'm overlooking something. Any assistance is appreciated. ![enter image description here](http://i.stack.imgur.com/9Cl0n.png)


title: Is there a standard projection for Sub-Saharan Africa?

text: With the help of d3.js, TopoJSON, and data from Glottolog I have made this map of Bantu languages that exist in Sub-Saharan Africa. I've used the standard Mercator projection so far. I know that there are always tradeoffs between different projection properties but I'd like to know if there is some standard projection that is recommended for Sub-Saharan Africa other than Mercator.


title: Connecting to SQL Server Geodatabase using ArcObject in C#

text: I am going to save a polygon that is created from points to remote geodatabase.I am new with arcobject and cant undrestand how it works..at this point I want to open a connection to database using arcobject and save geometry in it.I am using this code:                   public IWorkspace open_ArcSDE_Workspace(string server, string instance, string user,     string password, string database, string version)     {                        // Create the workspace factory.         Type factoryType = Type.GetTypeFromProgID("esriDataSourcesGDB.SqlWorkspaceFactory");         IWorkspaceFactory workspaceFactory = (IWorkspaceFactory)Activator.CreateInstance             (factoryType);              // Create the connection properties.         IPropertySet connectionProps = new PropertySetClass();         connectionProps.SetProperty("dbclient", server);         connectionProps.SetProperty("serverinstance", instance);         connectionProps.SetProperty("authentication_mode", "OSA");//The type of authentication to use. Valid values are DBMS and operating systems authentication (OSA). DBMS is the default mode and is not required.         connectionProps.SetProperty("user", user);         connectionProps.SetProperty("password", password);              // Open the workspace.         //IWorkspace workspace = workspaceFactory.Open(connectionProps, 0);         return workspaceFactory.Open(connectionProps, 0);               }      But it trows errors               enter code hereError HRESULT E_FAIL has been returned from a call to a COM component.      Can you help me and give me some examples that I can undrestand..The only source I have found for arcobjects is esri and I couldent find ant solutions for it. Thanks in advance


title: Issue using Matplotlib in QGIS

text: I am developing a QGIS-plugin and for graphing I want to use matplotlib, like many other QGIS -plugins. This worked fine under linux, but under windows with QGIS2.4 I am experiencing problems. I checked and matplotlib seams to be included in the standard windows-install of QGIS. But when people try to install the plugin, they get this report: **The plugin is broken: No Module named tkinter** Turns out even though I'm using pyqt, matplotlib requires tkinter and this seams be to be excluded from the QGIS installation. Is this a Bug? (should i report it?) Or is there a way to use matplotlib without loading tkinter? If I really need tkinter, how can I solve this? As tkinter is a binary python package, i can't include it in my plugin or is there an alternative (pyQWT is also binary, so also not a solution) I call matplotlib this way:               from matplotlib.backends.backend_qt4agg import FigureCanvasQTAgg as FigureCanvas     from matplotlib.backends.backend_qt4agg import NavigationToolbar2QTAgg as NavigationToolbar     import matplotlib.pyplot as plt      its the experimental version (1.2.0) of this plugin, if you want to try it: http://plugins.qgis.org/plugins/geopunt4Qgis/


title: How to copy attribute table values from one feature class to another

text: I am looking for a way to copy attribute values from one attribute table to another. I am mainly looking to do this using python, so would prefer a way to do this with a tool or a a stand-alone python script. A little background about my problem. I have a collection of polylines situated within a collection of polygons (a watercourse network). The crux of what I want to do is identify, spatially, which polygons hold which polylines and record this in the relevant field. The script below shows how far I have got. I have used python to select the relevant polygon and highlight the polylines that are situated inside it (selected features). However, I have so far been unable to copy these selected features into a particular field of another feature class. I have tried using both the Append (got 'warning' 000597) and Join tools but neither have worked.               import arcpy          reachesBoundaries = "C:\Users\maureen\Documents\ArcGIS\EDRN\Catchments\MainRiversCatchmentsPolygonMarch2013.shp"     edrnLink = "C:\Users\maureen\Documents\ArcGIS\EDRN\EDRN_LINK.shp"     reachField = "CATCHNAME"     reachName = "Aire - Lady Beck Catchment"          arcpy.MakeFeatureLayer_management(edrnLink, "edrnLayers9")     arcpy.MakeFeatureLayer_management(reachesBoundaries, "boundaryLayer8", '"' + str(reachField) + '" =' + "'" + str(reachName) + "'")          arcpy.SelectLayerByLocation_management("edrnLayer9", "WITHIN", "boundaryLayers8")      The script works up the point mentioned above, if anyone has any ideas on how to progress it, please let me know.


title: Convert from Julian Date to Calender Date in Field Calculator

text: Does anyone know of a way to convert from a Julian date to a regular calender date in either QGIS or ArcGIS field calculator? Thanks...


title: osm2pgsql import crashed!

text: I am currently importing a bz2 country file into PostgreSQL 9.1 database for Nominatim 2 using its `setup.php` that runs `osm2pgsql`. After processing the nodes. it seems to be stuck at a certain spot for hours without any changes and both CPU/memory usage for the server went down to idle levels. Any idea what happened? The full output can be found here. **Output Snippet**                                           addgeometrycolumn     --------------------------------------------------------------------------      public.import_polygon_error.prevgeometry SRID:4326 TYPE:GEOMETRY DIMS:2     (1 row)                                      addgeometrycolumn     -------------------------------------------------------------------------      public.import_polygon_error.newgeometry SRID:4326 TYPE:GEOMETRY DIMS:2     (1 row)          CREATE INDEX     CREATE TABLE     CREATE INDEX     CREATE SEQUENCE     CREATE TABLE     ALTER TABLE     CREATE INDEX     CREATE TABLE     ALTER TABLE     Partitions


title: Mapnik, Mapserver or Geoserver

text: I'm currently developping an ArcGIS based flex webviewer. ArcGIS is great for most of the things but there are some issues like performance and license costs. That's why I'm looking into possible alternatives. At the moment I found 3 possible alternatives, but I'm not sure what they are capable of so I'm wondering if some people on here have more experience with them. I have been reading the docs and will try to experiment with some of these but it will take a lot of time to find out all these requirements. So it would be of great help if you could point me in the direction of software I should try out. These are our requirements:   * fast   * easy scalability   * advanced labeling (different label colors based on attributes, labels placed in a rectangle)   * render cached data   * network tracing possibilities (for sewers, so trace upstream and downstream)   * printing (generating pdf's) nice to have:   * good flex or javascript viewer (I know this is different software, but I'm not sure if all 3 packages support those equally well)   * easy to start with


title: Comparison of the different open source GIS servers

text: I currently know the following open source GIS servers:   * Geoserver: java based, vector and bitmap support, also supports. Is now on par with mapserver (performance-wise).   * mapserver: c/c++, used to be the fastest server. Better for raster than vector?   * mapguide: do not know this   * mapnik: the new kid on the block? do not know much about it, but it looks appealing   * mapfish: as far as i know, only retrieves vector-data for display in openlayers. There does exist a rails implementation. If possible i would like some kind of comparison, which did you choose or prefer and why? I am looking at building a rails website, and need some kind of GIS server. I will need raster and vector data (clickable). Is this doable with one server? Backend database will be Postgis.


title: UMN MapServer alternatives

text: I am new to GIS application world. I would like to ask what are the other alternatives to **UMN MapServer** in the market preferably open source?


title: postgis distance calculation between points with lat, long coordinates

text: Team, I am new to postgis, but a postgres user for quite a while. I would like to compare the geospatial capabilities between postgresql and mongodb. I have the lat.lon, zip code of US locations; I ran across this blog entry at http://unserializableone.blogspot.com/2007/02/using-postgis-to-find-points- of.html Also, I see this the function ST_DWithin               SELECT * FROM geotable      WHERE ST_DWithin(geocolumn, 'POINT(1000 1000)', 100.0);      Please advise, I see that the lat/lon must be converted into points; but I am having issues with the best mechanism to perform indexed distance calculations. Update.. I have made progress and loaded a spatial aware table .. see below               create table public.disttest     (              zipcode varchar(5),         state varchar(2),         city varchar(55),         long numeric,         lat numeric          );          SELECT AddGeometryColumn( 'public', 'disttest', 'geom', 32661, 'POINT', 2 );          COPY public.disttest FROM 'c:\temp\zips.txt'          select * from public.disttest limit 10;               UPDATE public.disttest SET geom = ST_Transform(ST_PointFromText('POINT(' || long || ' ' || lat || ')',4269),32661) ;      Now the question is how to use the below type query for a 10 mile radiuus               SELECT * FROM public.disttest WHERE ST_DWithin(geom, geom, 16093);


title: Can I do Intersection queries in ECQL?

text: I am building a statistical analysis project using ASP.NET, Geoserver, and OpenLayers. I need to let the user draw a polygon on OpenLayers then find geometries(villages etc.) that intersects this polygon. Then get their attributes like temperature, population etc., make some mathematical calculations and present this information to user. My question is can I do the following operation: (SELECT population,temperature,humidity,height FROM villages WHERE geometry INTERSECTS userDrawnPolygon, SELECT data1,data2,data3 FROM tableX WHERE geometry INTERSECTS userDrawnPolygon .... goes on for multiple tables) with ECQL or CQL or some filter on Geoserver or do I need to develop this algorithm in my ASP.NET code - connecting to MSSQL? Thanks already.


title: Coordinate System Alignment

text: I have four shapefiles with different original projections that I am trying to bring into ArcMap 10.2. How do I set these four shapefiles to a single projection that will display the layers correctly? Any help would be appreciative.               Layer         Current Coordinate Ref Sys   Lines Up     Street        NAD83 UTM 16                 YES     Major Roads   NAD83 UTM 16                 YES     County        NAD83                        YES (probably)     Railroads     ??


title: FeatureClass to IFeatureLayer

text: I am trying to covert my feature class to feature layer so that i can pass "IQueryByLayer" over my feature class. Till now i have my feature class ---               FeatureClass className = new FeatureClass(cc);      want to do something like               IFeatureLayer newLayer = className.as---(FeatureLayer)      Can anybody point out a correct way to impliment it.


title: Layer is refreshed after an ajax request with openlayers

text: I can't understand why my vector is empy after this ajax request :               function loadKMLTMINT(terminalcode) {         MapPanel.el.mask();         Ext.Ajax.request({             url: 'php/GetKMLTM.php',             method: 'POST',             params: {                 terminalcode: terminalcode             },             success: function (r) {                  if (r.responseText.length > 0) {                     var data_tmp = JSON.parse(JSON.stringify(r.responseText));                     vector_term.parseFeatures(data_tmp, new ol.parser.GeoJSON(), ol.proj.get("EPSG:3857"));                     var first = true, precExtent = [];                     for (var vec in vector_term.featureCache_.idLookup_) {                         if (vector_term.featureCache_.idLookup_[vec].values_.Intermodal == 'Y' && vector_term.featureCache_.idLookup_[vec].values_.NAME.substring(0, 5) === terminalcode.substring(3, 8)){                             if (first) {                                 first = false;                                 precExtent = vector_term.featureCache_.idLookup_[vec].values_.geometry.rings_[0].bounds_;                             } else {                                 ol.extent.extend(precExtent, vector_term.featureCache_.idLookup_[vec].values_.geometry.rings_[0].bounds_);                             }                         }                     }                     olmap.getView().fitExtent(precExtent, olmap.getSize());                 }                  MapPanel.el.unmask();             }         });      }      I'm calling this function many time, in order to add some values to the layer vector_term. With the debugging I can see that my vector_term is full when I call the function, but when it reaches the `success: function(r)` it becomes empty ! Do you know if I'm doing something wrong? Should I use `addFeatures` instead of `parseFeatures`? Thank you EDIT : vector_term definition               var vector_term = new ol.layer.Vector({         source: new ol.source.Vector({             parser: new ol.parser.GeoJSON()         }),         style: new ol.style.Style({             rules: [                 new ol.style.Rule({                     filter : 'geometryType("polygon")',                     symbolizers: [                         new ol.style.Fill({                             color: ol.expr.parse('color'),                             opacity: 0.4                         }),                         new ol.style.Stroke({                             color: ol.expr.parse('color'),                             opacity: 0.8                         }),                         new ol.style.Text({                             color: '#0A3ABF', //'#165EF0',                             text: ol.expr.parse('NAME'),                             fontFamily: 'Arial',                             fontSize: 14,                             zIndex: 7                         })                     ]                 })             ]         })     });


title: Make Fields Required in the Relationship Inspector - Flex

text: I need to make some fields required in the attribute inspector of my relationship inspector. The code I use works with the editor attribute inspector but not with the relationship inspector attribute inspector. I get an error on this line: attributeInspector.fieldInspectors = [fieldInspector1, fieldInspector2] TypeError: Error #1009: Cannot access a property or method of a null object reference. When I use the code assist in the relationship I can see that there is a class called Attribute Inspector - Relationship Inspector (please see image), but I can't find this class to instantiate the attribute inspector in action script. Anybody know the right way to do this? Thanks               myEditor.featureLayers = [ myFeatureLayerstewardship];                     var attributeInspector: AttributeInspector = relationshipInspector.attributeInspector                     var fieldInspector1:FieldInspector = new FieldInspector();                     var fieldInspector2:FieldInspector = new FieldInspector();                          fieldInspector1.featureLayer = myFeatureLayercounty                     fieldInspector2.featureLayer = myFeatureLayercounty                     attributeInspector.fieldInspectors = [fieldInspector1, fieldInspector2]                     fieldInspector1.fieldName = "Acres"                     fieldInspector1.required = true      ![Image 1](http://i.stack.imgur.com/4nvEt.png)


title: Pinning ArcMap Table of Contents but not Attribute table window?

text: I have pinned by table of contents but I would like to view my attribute table as floating - i.e. I can make it much larger so I can view the data easily.


title: How to make a loop to MakeNetCDF Table View?

text: I was trying to write a code to extract grid point data from the NetCDF file. I am able to view and save the table data extracted as follow               #Import system modules     import arcpy          # Set local variables     inNetCDFFile = "F:/7071.nc"     variable = "prec"     outTableView = "F:/precipm"     rowDimension = "time"     dimensionValue = "lat 64;lon -146"     valueSelectionMethod = ""          # Execute MakeNetCDFTableView     arcpy.MakeNetCDFTableView_md(inNetCDFFile, variable, outTableView, rowDimension, dimensionValue,valueSelectionMethod)          outTable = arcpy.CopyRows_management(outTableView)     print      The problem is I have about 1000 points that I am going to extract the time series data from the netCDF file. I want a to write a script that loops through the csv grid point and save each table using the point ID. The grid point is given as Grid.csv               ID lat lon     1 10 20     2 30 40      etc


title: Is it possible to calculate statistics across multiple fields in arcpy?

text: The summary statistics command can be used in arcpy to calculate statistics - this tool runs against a single field, eg: > SUM—Adds the total value _for the specified field._ Is there an out-of-the-box method to calculate statistics _across multiple fields_? eg I'd like to calculate the MIN, MAX, MEAN and STDEV of Field1, Field2 and Field3, combined:               OID  Field1  Field2 Field3 Field4      1     23      24     25     26      2     12      13     14     15      3     17      18     19     20      One approach is to read through the table with a cursor, and calculate these values manually in a script. Is there a better/faster way? Edit: I need this to run on only the selected records, and these will change frequently.


title: How to zoom to point and place a pin?

text: Is there any means to combine the function of zoom to point and place a pin plugins in QGIS? Would it be possible to enter a xy co-ordinate and a pinpoint is placed at the point of the grid reference. My apologies if I have missed something obvious on other forums. GSherman maintains these both so I am hoping this query reaches him!? Thanks


title: Mosaic Dataset would not show raster unless you zoom in further to around 1:2,500 scale?

text: So when I make a copy the Mosaic Dataset File GeoDatabase to my external hard drive and then copy it to my computer. It would show the checkerboard at the full extent but when you zoom in it would show the imagery. I repair the Moasic Dataset to set it correct the path files. How do I make the imagery to show it at the full extent ? At the Full extent ---> ![enter image description here](http://i.stack.imgur.com/BODh6.jpg) At the 1:2,500 scale ---- > ![enter image description here](http://i.stack.imgur.com/IlJhs.jpg)


title: How to migrate Network Dataset in WGS84 to Web Mercator projection using ArcGIS 10.1?

text: regards I'm using Arcgis 10.1 I need to migrate a Network Dataset developed in WGS84 to Web Mercator projection, I could tell the steps to do this I want to automate this process in model builder, but I need to know the steps to make the change from WGS84 to Web Mercator projection could help me please


title: Problem with pgRouting and pgr_createTopology

text: I follow this tutorial: http://anitagraser.com/2013/07/06/pgrouting-2-0-for-windows-quick-guide/ and create a graph using OSM data (planet_osm_roads table). All goes with no error, but searching routes give me no results. My steps: 1) Create the graph (step3 of tutorial) 2) Check for a highway of my country:               select * from routes where name like '%BR-040%'          31575201    "BR-040/BR-050" 1728    37342     76240169    "BR-040"        41109   2174     31574903    "BR-040"        60030   35623     31574903    "BR-040"        60030   35623     31574903    "BR-040"        60030   35623     169184692   "BR-040/BR-050" 37342   41109     32829121    "BR-040"        28920   28921     169184696   "BR-040/BR-050" 37342   37342      Note is the same road, with a few segments. Lets try a small search: (I just look for a small source/target sequence in the result set)               SELECT seq, id1 AS node, id2 AS edge, cost, name       FROM pgr_dijkstra('SELECT id, source, target, st_length(the_geom) as cost FROM routes',  1728, 2174, false, false  ) as di        JOIN routes pt ON di.id2 = pt.id ;      Results:               0   1728    31575201    0.332977652802372   "BR-040/BR-050"     1   37342   169184692   0.348828028927443   "BR-040/BR-050"     2   41109   76240169    0.0616867326181856  "BR-040"     2   41109   76240169    0.0616867326181856  "BR-040"      Fine! The issue: Take a look in road segments result set: The record 31575201 have source = 1728 and target 37342 wich leads me to record 169184692 ( source 37342 and target 41109 ) in your turn, leads me to record 76240169 ( source 41109 and target 2174 ). At this point, the route is broken, because I can't find any record with source = 2174:               select * from routes where source = 2174      Results no records! What I'm doing wrong? Why I have a TARGET pointing to an inexistent SOURCE? I have many records pointing to itself, like the record number 169184696.               select count(*) from routes where source = target      Results 86.904 rows !! I've already play with the tolerance parameter in               select pgr_createTopology('routes', 0.0001, 'the_geom', 'id');      Changing from 0.01 to 0.00000001, give me some variations in self-pointing rows, but still having broken routes and can't go far away from 4, 6 steps in same road. Any real route searcing give me no results. Another test: Get the lower source of a road segment:               select  name, source, target from routes where name like '%BR-040%' order by source asc      Get the higher target of same road:               select  name, source, target from routes where name like '%BR-040%' order by target desc      By using this parameters in pgr_dijkstra function, I can see this road is broken! ( no results). Thanks!


title: toggling a list of paths in Leaflet

text: I drew a bus map of San Juan, Puerto Rico in leaflet.js. It's on dropbox. Using the layers control I let users toggle the bus routes they want to see. Unfortunately, the list of 45 or so bus routes is too long to display on the page. Route 21 is hidden behind the Leaflet link, and there are more routes beyond that. How can I let users see all their options in the layer control popup? ![enter image description here](http://i.stack.imgur.com/sg6GV.png)


title: problems using the Split Features tool

text: i'm using the split features tool to create two polygons from a bigger one. The problem is, when i save the split i have done, half of the polygon disappears. But sometimes, after doing some more polygons they appear again, and the last ones in with i was working disappear. Do you know what is happening and what is the solution? PS: This Shape was made by other people and i don't know how it was made, neither in which program has been created.


title: GRASS v.Overlay/QGIS Intersect Speed intersecting line and polygon

text: Im currently trying to convert my scripts from ArcGIS to QGIS but one aspect that is extremely slow (imo) is trying to intersect a polygon with a line feature. I feel as though there must be a better way to accomplish trying to clip the voronoi polylines with the original polygon itself i.e. : ![Voronoi_Lines](http://i.stack.imgur.com/lpt3Q.png) **by** ![Original Polygon](http://i.stack.imgur.com/yFmIx.png) These are some of the bench tests that i've done with these two layers. **ArcGIS 10.1 Intersect** 2min 24sec **Grass 7 v.Overlay** quit after ~1hour (~70% Selecting lines...) - Currently not supported in processing (Sextante) **Grass 7 v.split + v.Overlay** ~2min + 45min **QGIS 1.9** quit after ~3hours **ogr2ogr -clipsrc** quit after > 1hour **SAGA GIS Intersect Line-Polygon** ~15min -- acceptable but still 7x slower than ArcGIS and produces wrong results? The images below show SAGA GIS Intersect(with dangles removed) and same area with ArcGIS Intersect - Clearly SAGA did not intersect the line and polygon correctly as there should be no connection of the lines. These errors are scattered throughout the new output!!! In addition Line-Polygon does not retain the polygon intersect attributes to the line unlike the other methods mentioned above... ![SAGA GIS Intersect](http://i.stack.imgur.com/hbuEu.png) I have placed the **original dataset** in a GitHub repository here so hopefully others can make their own benchmarks and thereby provide suggestions? I notice that ArcGIS splits (cracking) the featureclass so I do not know if that is another possibility to improve intersection speed within one of the opensource alternatives mentioned above? ***Note that the reason I mention QGIS here is that I have made some scripts within the processing toolbox (formally known as Sextante) and i'd like to remain in that environment if at all possible. **System** - Intel Core i7-2760QM CPU @2.40GHz and 8GB ram Any help on increasing performance would be most appreciated! Cheers


title: Looking for USDA plant hardiness zone shape file

text: given that this is a USDA resource I really thought that it would be available online somewhere, but for the life of me I cannot find a shape file for the USDA plant hardiness zones anywhere. anyone know where there might be one? This link is the closest I found but it doesn't seem to offer any sort of shape file download: http://www.usna.usda.gov/Hardzone/ushzmap.html


title: Where are imported Mosaic Dataset images stored on ArcSDE server?

text: I am trying to learn more about Mosaic Datasets and specifically how they interact with ArcSDE (SQL Server 2008). I read in How raster data is stored in a geodatabase that "...In the case of the mosaic dataset, this table remains empty, since the raster data for a mosaic dataset is not stored in the DBMS; instead, the raster data is referenced from an image file..." Presumably this means that if I perform the Mosaic Dataset -> Add Rasters operation using files on the same machine as SDE then a file location is inserted into the murky depths of SDE. However, what if I connect to SDE remotely and Add Rasters from a different machine? I can only assume that SDE has a default directory to which the files are copied but I can't find anything about it, even in Mosaic datasets in a geodatabase in SQL Server. I did find a table named `<schema>.AMD_<mosaic name>_CAT` which includes a field `Uri` ("A serialized Uniform Resource Identifier (URI) that is used to re-create the raster dataset") but it means nothing to me. The test server in this case just has a single partition, but what if I wanted SDE to store mosaic dataset images on a separate drive for performance? How would I go about changing the default location, or even finding out what it is? Any input very much appreciated.


title: Google Maps legend ? Where do I find it for more information about them?

text: Can anyone point me where can I found information about Google Maps and their legend ? Since Google maps does not provide their legend, they display their colors like by showing the tan colors as schools, blue as lakes, green as could be trees, football fields and so on


title: How to work with GeoSMS?

text: Im looking for any libraries to work with GeoSMS standart on Android.   In standart says that geosms it just:               geo:lat,lon[,alt][;u=accuracy],       So maybe everyone work with it without any libs?


title: Quantum GIS Reproject a projected coordinate system with rotation

text: I have two shapefiles:- one in WGS84 Z30S and one in a custom, system based, projection. I need to convert the custom projection into WGS84 Z30S. I have a set of parameters including a shift in X and Y rotation, scale, and origins for rotation and scale. How do I format a custom coordinate system in QGIS to achieve this? Is there a better way of doing it? The syntax is foxing me at the moment.


title: Tiger Census data of 2013

text: I'm looking into geocoding U.S. locations without using Google or Bing. I found a service that geocodes according to the Tiger Census data of 2013. I have little knowledge about this Census and didn't get extremely far poking it around. When I use the api of the service that geolocates based on that census, I get several possibilities for the same full address in some cases. E.g. I got that one address in San Francisco is related to 3 different zip codes (that address, randomly picked, is "1354 Market St, San Francisco, CA") Is that at all possible in the U.S. zip code system? Do you think the Tiger Census data of 2013, is a very good source for geocoding? Thanks, Matan


title: Can ArcMap text formatting support mailto in text elements (ready for Export to PDF)?

text: I currently include an email address in a text element created using the Rectangle Text tool and its map gets exported to a PDF file. Is it possible to use ArcMap text formatting (or any other means) to assign a mailto: hyperlink to that email address so that it persists into the PDF file as something that can be clicked on to start an email?


title: Can any non-Terrago products produce geopdfs recognized by the Terrago toolbar in Acrobat?

text: Neither ArcMap nor GlobalMapper seem to be produce geopdfs that are recognized by the Terrago toolbar. Can anything besides Terrago products produce geopdfs that the Toolbar will accept?


title: How to get field names in pyqgis 2.0

text: I am trying to change an plugin from QGIS 1.8 to 2.0. In the 1.8 version I create a table of the field names like this:               layer = vlayer     provider = layer.dataProvider()     feat = QGSFeatrue()     allAttrs = provider.attributeIndexes()     provider.select(allAttrs)     fieldNames = []          while provider.nextFeature(feat)         fields = layer.pendingFields()         for (k, field) in fields.iteritems():              fieldNames.append(field.name())      I tried this for the 2.0 version with               layer = vlayer     fieldNames = []          for feature in layer.getFeatures():         fields = layer.pendingFields()         for (k, field) in fields.iteritems():              fieldNames.append(field.name())      But pendingFields() does not work for me. Is there now a new method to get the names of the fields?


title: How to get geolocation of a point in Autocad 2013?

text: I would like to get latitude and longitude of a point in Autocad 2013. Could you please help me with this? Thanks in advance. Regards, madhu


title: How can I import only a subset of the polygons in a Shapefile using QGIS?

text: I have a shapefile that is covering a huge area, with lots of polygons. The area I'm insteresting in is small and limited. Is there any way to only import (using "add vector layer") the polygons from the shapefile that covers my smaller area? How can I do this in QGIS? How do I specify the smaller area for my import?


title: arcpy fieldmappings: can't set output field name for fieldmap

text: I am testing out some code in the Python window of ArcGIS so that I can insert it into a standalone script later. I want to create a copy of a table and add a new field that is the result of concatenating two of the existing fields using a space as a delimiter. I am sort of following the suggestion below to use field mapping which is from this page http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//001700000065000000 because I ultimately want to use the JoinField command to join two tables based on two attributes: "To merge two or more fields in the join table before the join, first export the table or feature class using the Table To Table tool and merge using the tool's field map." My code does the concatenation or merging just fine, but I cannot get it to rename the output field. In my resulting table, the name of the new field is the same as the first field I add to the fieldmap (ie, "COC_NAME" in the code below.)               fieldMappings = arcpy.FieldMappings()     fieldmap_join = arcpy.FieldMap()     fieldmap_join.joinDelimiter = " "     fieldmap_join.mergeRule = "Join"     fieldmap_join.addInputField("COC_DATA", "COC_NAME")     fieldmap_join.addInputField("COC_DATA", "GRID_ID")     fieldjoin = fieldmap_join.outputField     fieldjoin.name = "JOINID"     fieldmap_join.outputField = fieldjoin     fieldMappings.addTable("COC_DATA")     fieldMappings.addFieldMap(fieldmap_join)     arcpy.TableToTable_conversion("COC_DATA", "C:\\HEA_GIS.mdb", "COC_MOD", "", fieldMappings)      Can someone point out what am I doing incorrectly or perhaps provide another solution?


title: LiDAR data density, returns per square meter

text: How can I compute a raster showing the number of return points per square meter using SAGA GIS or other open source software? My goal is to obtain a mask which would show me the places (cells) where there were no hits and would help me in interpreting LiDAR-derived data.


title: How to pass a variable (FilesToAdd) from VBA script (interface) to a Mapbasic program?

text: Option Explicit          Public FilesToAdd As Integer          Private Sub TEST_Click()          FilesToAdd = 1          Dim b As Variant     Dim Prog As String     Dim File As String          Prog = "C:\Program Files\Mapinfo\Professional\Mapinfow.exe"     File = "C:\Mapinfo\Combine\SelectCarrier.MBX"          b = Shell(Path + " " + File, vbNormalFocus)          End Sub


title: VBScript for Advanced Labeling in ArcGIS

text: I found a template I thought might be useful to create a table of sorts for a label. I modified it with my variables, but it still doesn't seem to be working. I'm not very familiar with VBScript, but I didn't really change anything from the original other than adding my variables. I don't see anything in the code that would be dependent on the number of variables, original had 4, so I don't think that's a problem. Any ideas how I can make this work. I want the result to look like the attached. And here is the link to the Template ![Table](http://i.stack.imgur.com/yCCaj.jpg)                   Function FindLabel([SWSID], [Stn_Name], [First_Date], [Last_Date], [Percentage], [First_Da_1], [Last_Dat_1], [Percenta_1], [First_Da_2], [Last_Dat_2], [Percenta_2], [First_Da_3], [Last_Dat_3], [Percenta_3])         iMxLblSz = 0          iMxValSz = 0          iSpace = 5          m = Array(Array(“SWSID”, [SWSID]), Array(“Station Name”, [Stn_Name]), Array(“First Date in DB”, [First_Date]), Array(“Last Date in DB”, [Last_Date]), Array(“Percentage Complete”, [Percentage])_      Array(“First Date in DB”, [First_Da_1]), Array(“Last Date in DB”, [Last_Dat_1]), Array(“Percentage Complete”, [Percenta_1])_      Array(“First Date in DB”, [First_Da_2]), Array(“Last Date in DB”, [Last_Dat_2]), Array(“Percentage Complete”, [Percenta_2])_      Array(“First Date in DB”, [First_Da_3]), Array(“Last Date in DB”, [Last_Dat_3]), Array(“Percentage Complete”, [Percenta_3]))              For i = 0 To UBound(m)         j = m(i)         If (Len(j(0)) > iMxLblSz) Then          iMxLblSz = Len(j(0))         End If         If (Len(FormatNumber(j(1), 0, 0, 0, -1)) > iMxValSz) Then          iMxValSz = Len(FormatNumber(j(1), 0, 0, 0, -1))         End If         Next              strGT10M2 = “</CLR>”              FindLabel = strGT10M1 + “<ACP><BOL><FNT name=’Arial’ size=’10’>” & _         [SWSID]+ "-" + [Stn_Name] & “</FNT></BOL>” & vbNewLine & _         “<FNT name=’Lucida Console’ size=’8’>” & vbNewLine               FindLabel = FindLabel & String(iMxLblSz + iMxValSz + iSpace, “_”) & _          vbNewLine               For i = 0 To UBound(m)          j = m(i)               k = (iMxLblSz - Len(j(0))) + iSpace + _         (iMxValSz - Len(FormatNumber(j(1), 0, 0, 0, -1)))               FindLabel = FindLabel + j(0) + String(k, “.”) + _          FormatNumber(j(1), 0, 0, 0, -1) + vbNewLine          Next               FindLabel = FindLabel + “</FNT></ACP>” + strGT10M2         End Function


title: Link to pdf file through mouse clicking

text: I need help in identifying the feature attributes displaying a pdf file when click in that feature. I see many examples but I don't have success. I try this code:               command = "firefox";     imagerelpath = "images_test/test_image.jpg";     layer = qgis.utils.iface.activeLayer();     import os.path;     layerpath = layer.source() if layer.providerType() == 'ogr'       else (qgis.core.QgsDataSourceURI(layer.source()).database()       if layer.providerType() == 'spatialite' else None);     path = os.path.dirname(str(layerpath));     image = os.path.join(path,imagerelpath);     import subprocess;     subprocess.Popen( [command, image ] );      but it gives me always an error:               WindowsError: [Error 2] The system cannot find the file specified      Basically, I have a shapefile and I want to click in a feature and open a pdf file. I want to develop that in a python plugin.


title: How to compare a set of images with a shapefile?

text: I have around 150 images (rasters), each one is just a map of certain "classes", e.g. pixel value 1=forest. I also have a shapefile that consists of roughly the same area as the 150 images constitute when put together. The shapefile also includes the clasee 1=forest. How could I go around to compare if the "classes" match, and how well do they. Any ideas? I'm using python and Qgis. Any idea would be appreciated. Greetings. Julián.


title: How to repair topology of a shapefile?

text: I performed an object-based image analysis using ENVI EX, and each of the classified features were exported as a shapefile. I transferred those shapefiles into ArcMap 10, and tried to intersect them with other shapefiles of the same area that were acquired from GeoBase. When I try to run the intersect tool, I get an error message stating that the topology of the shapefiles extracted from ENVI EX are invalid. I tried to turn on the topology tool to repair the topology of those shapefiles, but the buttons are greyed out, and I can't seem to click on them. The topology tool does not work for any of the shapefiles that I have. Any solutions? How I can I fix the topology and run the intersect?


title: GRASS GIS: How do I run a batch (bash?) script?

text: I am in a GRASS database (in command shell), where I want to run a loop function. For example, I want to do an operation on many raster files that are in the database. I know this must be done through a txt file. What i am confused about is where I should put this file, and how exactly I should call it. I am also a bit confused if I should be writing a batch or a bash file. GRASS is UNIX based which makes me think I should be looking to bash, but I am in windows command shell, where batch scripting applies. Thanks for your help!


title: Zooming into landsat-8 produces pixellated view

text: I am working on a project were I am using LANDSAT 8 images for remote sensing. The issue that I have is that when I zoom into them the photo is pixellated. I can not see a thing. When I zoom out and view it from a distance it works, but I need to be able to view it closely. What I need is to have a clear view at say a scale of 1:1,250. Is this possible? I have taken and looked at pan- sharpening the bands, is there anything else that I may be missing? The program that I am using is ArcMap.


title: Python add-in won't buffer from output layer

text: I have a simple add-in script for Arcmap 10.1 to create a toolbar, select a layer from a drop down list of the current layers in the project, buffer that layer four-miles, then use the resulting four-mile buffer output as an input for a six-mile outside buffer. I have a work flow where I repeatedly need to create a four-mile buffer with a six-mile outside buffer surrounding this. Everything functions up to and including the four-mile buffer layer creation, but I get no six-mile outside buffer being created. I have been working with ESRI support but having trouble resolving this. Any ideas?               import arcpy     import pythonaddins     import os          class ComboBoxClass1(object):         """Implementation for DataRequestBuffer_addin.combobox (ComboBox)"""         def __init__(self):              self.editable = True             self.enabled = True             self.dropdownWidth = 'WWWWWWWWWWWWWWWWWWWW'             self.width = 'WWWWWWWWWWWWWWWWWWWW'         def onSelChange(self, selection):             arcpy.env.overwriteOutput = True             # When a new layer is selected, create a 4 mile buffer around that layer then create a 6 mile buffer surrounding the outside of the 4 mile buffer.             layer = arcpy.mapping.ListLayers(self.mxd, selection)[0]             outPath = os.path.dirname(layer.dataSource)             print outPath                  # Buffer 4 miles             print 'Buffering %s with 4 mile buffer' % layer             arcpy.Buffer_analysis(layer, outPath + os.sep + str(layer) + "_4milebuff", "4 Miles", "FULL", "ROUND", "ALL", "")             buff6 = outPath + os.sep + str(layer) + "_4milebuff"             # Buffer 6 miles surrounding outside of 4 mile buffer             print 'Buffering %s with 6 mile outside buffer' % buff6             arcpy.Buffer_analysis(buff6, outPath + os.sep + str(layer) + "_6milebuffoutside", "6 Miles", "OUTSIDE_ONLY", "ROUND", "ALL", "")             arcpy.RefreshActiveView()                   def onEditChange(self, text):             pass         def onFocus(self, focused):           # When the combo box has focus, update the combo box with the list of layer names.             if focused:                     self.mxd = arcpy.mapping.MapDocument('current')                     layers = arcpy.mapping.ListLayers(self.mxd)                     self.items = []                     for layer in layers:                             self.items.append(layer.name)              def onEnter(self):             pass         def refresh(self):             pass


title: PostGIS: Calculate the length of just the portion of linestring that intersect with a polygon

text: I want to calculate the portion of the linestring that lies within the limits of a polygon. Better with an image ![Distance within province](http://i.stack.imgur.com/Phx8Y.png) What is the length of the portion of linestring within Province A? and within Province B? I'm looking for a query that returns something like:               province_name | length   |     --------------------------     Province A    | 30.5342  |     Province B    | 19.4321  |      My approach is a PL/pgSQL script that loops over all the points in the linestring and set a mark in the first point that is in province B (point N). Then, calculate the length from point 0 to point N (distance within province A) and the length from point N+1 to last point (distance within province B). But sounds rusty, and difficult to extend to an indefinite number of provinces. Is there any better approach? Suggestions welcome. I'm using PostGIS 2.0, by the way. **UPDATE** : @HeyOverThere suggests a query that almost works, with a little modification:               SELECT a.name, sum(ST_LENGTH(ST_Intersection(a.geom, b.geom))/1000)     FROM polygons a, lines b     WHERE ST_Intersects(a.geom, b.geom) group by a.name order by a.name      That query returns 2 columns, with the province name and a length. But the length is **exactly twice** the value it should be. Does it make sense? More info: the provinces are really multipolygons (but all share a the province name as common field), and the line has been built from gps tracks (points) and st_makeline. **UPDATE 2:** As @GetSpatial suggested, I've tried using subqueries to sum and group the results of the first query. Same result. The length of each intersection is correct. So, maybe the intersection is two times as long as I expect (?). Here, the 3 things I've tried. All of them get a distance 2x the expected one (link to real data at the end): First try, with one query:               SELECT a.province, sum(ST_LENGTH(ST_Intersection(a.geom, b.the_geom))/1000) as length_per_province     FROM provinces a, line_tracks b     WHERE ST_Intersects(a.geom, b.the_geom) group by a.province order by a.province;      Second try, with CTE:               with line_segments as (     SELECT a.province, ST_LENGTH(ST_Intersection(a.geom, b.the_geom))/1000 as length     FROM provinces a, line_tracks b     WHERE ST_Intersects(a.geom, b.the_geom) order by a.province     ) select province, sum(length) as length_per_province from line_segments group by province;      Third one, with subquery, as suggested:               select province, sum(length) as length_per_province from(     SELECT a.province, ST_LENGTH(ST_Intersection(a.geom, b.the_geom))/1000 as length     FROM provinces a, line_tracks b     WHERE ST_Intersects(a.geom, b.the_geom)) as foo group by province;      I get this, in all cases:               "province";"total_length"     "p1";62.8687677452213     "p2";40.254207910339      The total length of the line:               select sum(st_length(the_geom))/1000 as total_length from line_tracks      The result               "total_length"     51.5614878277802      And the 2 tables here (SET client_encoding = 'LATIN1'): https://dl.dropboxusercontent.com/u/6599273/gis_data/data.sql.tar.bz2 If anyone could try, that would be great. Anyway, many thanks.


title: In ArcGIS 9.3, how do I split a 3-band raster (.tiff) file and summarize the band values?

text: Using ArcGIS 9.3, I would like to make a new, single-band raster file that includes only one value: a summary of the three existing bands. The composite band is currently a composite .tiff file. Something like Band1 + band2 + band3 = newValue , for each pixel in the raster. The new raster would contain only one output band, with the values for that band = newValue.


title: What is the structure of Standard Interchange Format (SIF) used by spatial data providers?

text: What is the structure of Standard Interchange Format(SIF) used by navteq and other spatial data providers? Any documentation/sample files? Is there a SIF parser available in C/C++/Java/PHP ?


title: How to add GDAL to existing Python installation?

text: I would like to download and install GDAL for Python version 2.6 [MCS v.1500 32 bit(Intel)] on win32 that supports reading FileGDB's with Python. But I am not sure how to go about this. I found one exe at this link but it was not built with the FileGDB included. I also tried the OSGeo4w Installer but it currently installs its own 2.7 version of Python rather than allowing the GDAL to be installed with the existing version of Python. I must work with the 2.6 version installed on our systems. Can anyone be of assistance in locating a GDAL build that supports FileGDB's, can be installed with my existing Python version 2.6 and provide detailed instructions?


title: ArcGIS Server 10.1 and Postgres Geoprocessing Service issue

text: I'm hitting an error when publishing a model from ArcGIS Desktop 10.1 to ArcGIS Server 10.1. The model I want to publish uses a layer from a Postgres DB defined via a Query Layer in my MXD. The model is very simple for now, run the Select function on the Query Layer and return the result as a shapefile. Publishing fails with a _"Failed to create service definition"_ message and a _"ERROR 001270: Consolidating the data failed"_ message in the results window. This error links to http://resources.arcgis.com/en/help/main/10.1/index.html#//00vp0000002r001270 I have taken all of the steps mentioned therein - shortening service names, paths etc, and ensuring that the DB connection is registered in my ArcGIS Server instance. I can add the PG layer as a standalone map service which works fine. Thanks for any help or suggestions! EDIT: I should add that my Postgres instance (9.1.3) is not configured as an ESRI Enterprise Geodatabase, so that could be an issue? PostGIS 1.5.3 is installed on the DB, which is currently used for web mapping (MapServer) and hopefully data connections into ESRI GP services.


title: On the fly projection not working in QuantumGIS

text: I am trying to put a point layer of cities on top of a shapefile of counties in QuantumGIS. I have on the 'On the Fly' turned on and both files are in the same projection system. However, the point file is not showing up on the county file and both have different coordinate systems for some reason. How to I get the layers aligned?


title: Why is my .csvt not recognized?

text: I am trying to associate a .csvt file with a .csv file, so QGIS doesn't insist on calling all the fields text strings (thereby preventing the use of graduated symbols). I used the syntax described by other answers here and on other sites in the .csvt file, ie: "String","Integer","Integer" etc. I saved it in Text Editor on my Mac with the same name as the .csv file and in the same folder. For whatever reason, QGIS is not recognizing the .csvt file's information and continues to deem every field of the .csv file a text string. I do not know if this is the fault of the Text Editor not recognizing the file format, or of QGIS not recognizing the .csvt when it ought to. Anyone have troubleshooting ideas or work-arounds?


title: Polygon to Raster Tool Problem

text: I am trying to create a raster data set from county level polygon data for the US. I have successfully created one raster from a dataset, but when attempting another, it does not work. I am wondering if I am missing something. It's in ArcMap 10, I have tried moving the cell size up to 100 (it took the recommended cell size for the other data set) and have also changed any NULL values to zero. I get the "999999 : Error executing function". I have also tried the Feature to Raster tool. Same error message. Any ideas or troubleshooting tips would be great! Thanks.


title: Openlayers: How to highlight a point on mouse over a line

text: Any suggestions on how to highlight a point on mouse over a line, I've read points from a javascript array, lines are drawn by using them and I want to highlight some points in that array by moving the mouse pointer through the line and near to that point. However, i'm stuck on actually how to get started with this.   ![This image contains the lines:](http://i.stack.imgur.com/y7PEu.jpg) Any guidance and help is appreciated..!


title: geoxml3 parseKMLString function not working?

text: I am unable to use Parser.parseKMLString(). Its giving me a error saying **Parser.parseKMLString() is not a function**. What i am trying to do is, read the KML file on the server and return the entire file as a json string which is then passed to the parseKMLString function. But i am getting the above error.


title: ArcGis 10 How judge ArcMap dock window ArcCatalog is open?

text: 1. From ArcMap 10 open ArcCatalog window, how to judge ArcCatalog isOpen?(Using C#)


title: What would be 'normal' processing time for a large postgis query

text: I have two datasets which I want to spatially join (with an intersect). The first dataset consists of 9 million addresses. The second dataset consists of 12.000 neighbourhoods. I want the add the neighbourhood to the address. Some information about the system: Microsoft Windows Server 2003 (SP-2) Inter(R) Xeon(R) CPU X5650 @ 2.67 GHz 2.66 GHz, 4,00 GB of RAM I log in via Remote Desktop The processing time for this query is about 50 minutes. Quite long I would say. Can someone tell what would be the general processing time? Both geometry columns have an index. First I add a column               ALTER TABLE temp.test ADD COLUMN buurtnr INTEGER;       Then I want to fill the column by doing a spatial join.               UPDATE temp.adresses a SET buurtnr = n.buurtnr FROM temp.neighb n WHERE st_intersects (a.geometrie, n.the_geom)       I have to run the query monthly, together with similar query. The whole process now takes about a day. I want to make it faster.


title: How to generalize and vectorize a LandUseClassification raster file?

text: I want to vectorize that rasterfile with its 6 landuse classes. The classification is very noisy and have to be generalized before. The result should be a polygone vector file. I use SAGA, (GRASS) and QGIS. ![enter image description here](http://i.stack.imgur.com/igKAf.png)


title: Why is this .prj file causing the shapefile to render in meters instead degrees?

text: So I've got a shapefile that has polygons of countries, with coordinates in decimal degrees:                Country (String) = Algeria      POLYGON ((6.928073999999953 36.88360652,                6.949973999999969 36.883306519999991, ...      `gdalsrsinfo` on the .prj file outputs:               PROJ.4 : '+proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +a=6378137 +b=6378137 +units=m +no_defs '          PROJCS["Mercator",         GEOGCS["GCS_unnamed ellipse",             DATUM["unknown",                 SPHEROID["Unknown",6378137,0]],             PRIMEM["Greenwich",0],             UNIT["Degree",0.017453292519943295]],         PROJECTION["Mercator_1SP"],         PARAMETER["central_meridian",0],         PARAMETER["false_easting",0],         PARAMETER["false_northing",0],         UNIT["Meter",1],         PARAMETER["latitude_of_origin",0.0]]      To my untrained eye, the `GEOGCS[]` bit looks like it's defining the coordinate system in degrees, so this should work. However, if I include the .prj file, Google Maps Engine thinks the coordinates are in meters and renders everything as a tiny blob around (0,0). If I drop the .prj file, it defaults to degrees and returns a much more sensible rendering. What's wrong, and how do I fix it? (Note the the shapefile is from a third party and produced with an unknown toolset.)


title: Routes Consult with pgrouting / openlayers

text: I'm starting with pgrouting. After suffering a bit unable to put the data in osm postgis using osm2po with the excellent tutorial http://anitagraser.com/2011/12/15/an-osm2po-quickstart/ ran the sql generated by osm2po visualized data with http: / / localhost: 8888/Osm2poService routing is fantastic and has also insert the data in my town in postgis successfully. Then I started reading the manual and pgrouting in http://docs.pgrouting.org/2.0/en/doc/index.html me questions arose: 1 - The data provided by osm2po are already topologically arranged just so I do consultations? 2 - The fields of my table are:               id     osm_id     osm_name     osm_meta     osm_source_id     osm_target_id     clazz     flags     source     target     km     kmh     cost     reverse_cost     x1     y1     x2     y2     geom_way      Can I get the total distance in miles doing a "sum" of the KM field, is this correct? Where can I learn what each field, for example, do not know what the field clazz, cost, reverse_cost ....? 3 - The flow of my application the user is typing two addresses, I'll get the latitude and longitude of the point of origin and destination and select the launch in postgis, how do this if the source and target parameter settings are whole? 4- Does anyone have any examples of how to do this with openlayers? * * * Edited here I'm getting the distance of the route like this:               SELECT                    sum (cost) cost AS      FROM pgr_dijkstra (                      'SELECT id, source, target, the cost km FROM at_2po_4pgr'                       50347, 32988, false, false               );       Where 50347 and 32988 are my origin and destination. I do not understand what this format is, are coordinated? The coordinates of my geocode comes in decimal degrees are respectively 51.13, -29.77 and -51.21, -30.04 if I use them so do not run the sql. Could you explain me I can pass the coordinates for the sql?


title: Embedding a open source topographic Earth (maybe Google Earth) onto a website?

text: For a project, I want to create a website that includes something that resembles Google Earth. For this website I would like the starry night background with the Earth in the middle (like that on Google Earth). I do not need multiple layers just a simple topographic Earth and maybe a layer that includes the clouds. Is there anyway that I could embed the Google Earth map into the website and reference them for credit? If the above is not possible, is there anyway to get a 3D topographic Earth layer open source that I can embed into the website? It needs to also be interactive like the one on google earth (as in you can zoom in and out and spin the Earth to view different parts) Thank you


title: Print and Export Functionality on web GIS application built with OpenLayers+GeoExt+Extjs+Heron Configuration

text: I have used **OpenLayers+Extjs+GeoExt+Heron** configuration to create a thin mapping client. I have added print icon to the client using `heron.options.map.toolbar` (`type: printdialog, options: {printURL etc}`). I have also the mapfish printing e.g. `config.yaml` in the geoserver data directory and validation of this installation via `xxx/geoserver/pdf/info.json` shows the print options as expected. However back to the app, whenever I click on the print icon, I get an error saying:               Error getting print options from server: undefined.        Has anyone experienced this and have solved this problem. I would appreciate any advice on how to fix this error. Thank you in advance. Any help would be very much appreciated!! Barbara


title: Recovering shape layer which was accidentally overwritten?

text: I duplicated a layer set, as it was more convenient (same style and attributes). It was working fine and both the original and new layers were showing. Then I open it up another day and the original layer looks to have vanished, and only shows the new shape layer (polygons/lines). Upon further investigation the new layer in now referencing the original layer shape files, the new ones are nowhere to be found. The big problem is that it resulted in the deletion of all the original old lines/polygons. Is there a way to retrieve the original lost shape data and to make them reappear on the layer? E.g. restoring a file to an earlier date? I'm new at this and would appreciate a step-by-step answer, links can be confusing without directions.


title: Spatial join including attributes - kind of like right join

text: I've run in to a little problem. I have two sets of polygons with important attributes. I need to join over the attributes from the first set to the other, but I can't seem to figure out how. Theoretically it seems like an SQL right join, I think. Can you help? The only way I can get this to work, I think, is the somehow copy the blue polygon. The first problem is that not all of the red polygons are completely inside the blue, but I have selected by centroid (I need to do this). The second problem is that there could be 1, 2, 3 or more red polygons per blue. The example below would be Blue1+Red1 and Blue1+Red2 and so on. I have ArcEditor 10.1, QGis and FME. ![I need to join the attributes from the red polygons over to the blue polygon. ](http://i.stack.imgur.com/GJfTJ.jpg)


title: Python syntax error -can't figure out why?

text: This should be simple. I hope I don't get down votes! I just got the following code from an article about using Google Maps API and ArcGIS to calculate routes without network analysis and just can't get it to run. I am getting a syntax error on line 7. I have tried to type in other values before it and they are fine, also tried the to change the value "fromFile" to something else and it still gives the same error. Any ideas?               #Estimating O–D travel time matrix by Google Maps API: implementation, advantages, and implications     #http://dx.doi.org/10.1080/19475683.2011.625977     import arcpy     import urllib      import time      from xml.etree.ElementTree import XML, fromstring, tostring     fromFile = arcpy.GetParameterAsText(0)     toFile = arcpy.GetParameterAsText(1)     Resultfile = arcpy.GetParameterAsText(2)      Result = open(Resultfile,”w”)      Result.write(“FromFID,toFID,TravelTime”)      Result.write(“\n”)      fromCursor = arcpy.SearchCursor(fromFile)      toCursor = arcpy.SearchCursor(toFile)      fromRow = fromCursor.reset()      fromRow = fromCursor.next()      while (fromRow!=None):      fromX = fromRow.shape.centroid.X      fromY = fromRow.shape.centroid.Y      fromFID = fromRow.FID      arcpy.AddMessage(str(fromFID))      toCursor = arcpy.SearchCursor(toFile)      toRow = toCursor.reset()      toRow = toCursor.next()      while (toRow!=None):      toX = toRow.shape.centroid.X      toY = toRow.shape.centroid.Y      toFID = toRow.FID      arcpy.AddMessage(str(toFID))      googletext = “http://maps.googleapis.com/maps/api/directions/xml?origin=(“ + str(fromY) + “,” + str(fromX) + “) &destination=(“ + str(toY) + “,” + str(toX) + “) &sensor=false”      time.sleep(3)      xmlfile = urllib.urlopen(googletext)      xml = xmlfile.read()      value = “NA”      dom = fromstring(xml)      nodelist = dom.getchildren()      if (nodelist[0].text =  = “OK”):      arcpy.AddMessage(nodelist[0].text)      route=nodelist[1]      leg=route.getchildren()[1]      duration = leg.find(“duration”)      value = duration.getchildren()[0].text      else:      arcpy.AddError(nodelist[0].text)      Result.write(str(fromFID))      Result.write(“,”)      Result.write(str(toFID))      Result.write(“,”)      Result.write(value)      Result.write(“|n”)      toRow = toCursor.next()      fromRow = fromCursor.next()      Result.close()      del fromCursor      del toCursor


title: Adding QGIS Server to an existing Windows web server

text: What is the best way to add the QGIS Server capabilities to an existing web server running under Windows? Specifically, I have WAMPServer 2.2 installed and running on Windows Server 2012, and would now like to add the capability to serve WFS layers. On a different machine (running Windows 7) I have successfully set up QGIS Desktop and Server via the OSGeo4W installer, but now I want to achieve the same result without installing a second version of Apache. The tutorials I have found are either for linux/unix-type systems, are several years old and use pre-2.0 versions of QGIS, or use the OSGeo4W installer to create a new Apache installation. I haven't been able to find an individual installer for a Windows version of QGIS Server, and it doesn't appear to be part of the QGIS Standalone Installer.


title: MapInfo MI_STYLE Column Definition

text: I have been looking for some documentation or post which tells me what the column definition for a MapInfo MI_STYLE column should be in Oracle. My guess is that it's VARCHAR2(128) or something like that since you have to give Pen, Brush, Line and Symbol style information. I just don't want to be sitting around debugging problems created from an incorrect field definition so if anyone can give me some info or point me to the documentation that would be great.


title: json to shapefile in javascript

text: I know, there are some libraries to convert .shp to .json even with javascript, but haven't found a way to convert .json back into .shp that could be implemented in web enviroment. Do you know or have experience with something? Or what are my possibilities? Thanks


title: What are the most efficient Geographic Weighted Regression packages?

text: What are the most memory efficient open source packages for calculating a geographically weighted regression (GWR)? I am in a situation where I need to do a geographically weighted regression on a set of points where training data consists of about 40,000 observations and each observation has about 20,000 variables. I have implemented a version of GWR myself using a combination of Python Numpy/SciPy and PostGIS. I solve the regression using a matrix algebra approach, but this fails due to memory issues when I have dense, feature rich systems with many observations. One way to get around the memory issue is use an iterative approach for finding a line of best fit, such as an incremental gradient descent. I'm thinking it should work something like (http://www.eecs.wsu.edu/~cook/dm/lectures/l5/node14.html). Incremental Gradient Descent is described pretty well here in pages 4-7 (http://cs229.stanford.edu/notes/cs229-notes1.pdf). Obviously I could implement this myself, but I was hoping maybe someone else had already coded something similar.


title: Is it possible to update the arcgis bundle file dynamically

text: We are using Arcgis Server to generate the tiles with Compact format,then we depoly the tiles by WMTS or something else. Now we met some problem when we update the source data,then we will re- generate the bundle file which contain the update area. But we have to generate the whole bundle even we have only a single tile updated, this is a waste of time. Then I wonder if we can modify the bundle file dynamically by program to put the replace the single or multiple tiles? Example: Suppose we have a bundle file named `R21212C2112.bundle`,and it contains 2000+ tiles, and we suppose there is a tile named 'tile1'. Now our source data have some update, and all the change are inside the `tile1`, now we have to re-generate the whole bundle `R21212C2112.bundle`,we wonder if we can only re-generate `tile1` and put it to the original bundle?