ArcGIS Heads up digitizing - Rooms and doors

I am digitizing the floorplans of various schools in our county. This is for the Fire Department so they want to be able to see where doors are within the building. The doors do not necessarily have to be indicated; a gap in the wall is sufficient. It would be much more efficient for me to use the 90 degree tool and knock out these classrooms in full. Is there some sort of eraser feature that might split a line segment in two wherever it intersects? That way I could basically create a rectangle for the room and erase where the door might be to leave a gap. Or are there any other tools that could make this process easier? I've got a few dozen schools to digitize, and by far the most time consuming process is creating small segments around the doors just to leave that gap. Cheers!


OpenLayers crashes stock browser on Android 4.0.4 devices

* * * **Update:** The issue appears to be with OpenLayers Release 2.13 dev. I mistakenly mentioned 2.12 as being problematic but have retested with no problems. Also, no problems with 2.13rc3 ## **Update 2:** This issue appears to exist again in 2.13.1 We're developing a mapping tool using Openlayers ~~2.1.2 (issue also confirmed on~~ 2.1.3 ~~).~~ I've tested on 4 devices running Android 4.0.4. The map panning/zooming works for a few seconds before the browser completely locks up - no longer able to pan/zoom or navigate to another page. The browser then either crashes or I need to force stop in application settings. Devices running 4.1.1 have no problems. This is also an issue with OpenLayers own mobile example page: http://openlayers.org/dev/examples/mobile.html Again, after a few seconds, the browser freezes. I see no reports of this online. The OpenLayers documentation says it is supported on Android 2.1+ Has anyone else run into this issue?


What is difference between simultaneous users and simultaneous connections in ArcGIS for Server 10.1?

What is the difference between simultaneous users and simultaneous connections. We're looking at Workgroup standard vs Enterprise Standard. Here's what we want to do: Use Arc Gis Server to serve data to app. 30 vehicles using Ipads via Sierra Wireless 4g modems. This data would be maps for Right of Way, outage locations etc. We would also like to be able for users to send data to other users. Thanks, Russell


Is it possible to use UTF-8 characters for shapefile attributes?

I heard that according to a specification shapefile attributes stored in dbf- files cannot have UTF-8 as encoding. And then I found out that an additional file with the ending CPG can be used to define the code page to identify the character set used in the dbf file. For me this sounds contradictory. So my question is now - is it possible to use UTF-8 characters for shapefile attributes? and - additional question, since encoding and character sets are always an unpopular topic, do You know a page where the topic is described in detail?


Best way to get the overlapping segment of two lines

Could anyone help me figure out the best way to get the overlapping line segment of 2 line sdo_geometry objects? my query looks like this:               select SDO_GEOM.SDO_INTERSECTION(SHAPE, :SHAPE, 0.5) as geom     from MyTable      where MDSYS.SDO_OVERLAPS(SHAPE, :SHAPE) = 'TRUE'      the geometry objects in the table are all lines as is eth parameter being passed in, but the geometry being returned for all results is a point, not a line.


travel survey multi stop routing problem

I have a travel diary survey dataset that I am attempting to transform into more useful information. Currently, the data comes in the format of:               Survey# ,  PlaceVisited#,  TravelMode, TripPurpose, LocationType          1,      1,       NA,        NA,         home          1,      2,       auto,  doctor,     office          1,      3,       auto,  shop,       grocery_store          1,      4,       auto,  NA,     home                      where each place the survey respondent visits is a separate record (row) in the survey table, but this is really one multipurpose trip. What I would like to do with the data is separate each actual trip into a route with various points and purpose information included as (i think) a line using shortest path on a road network (I don't have the route actually taken). Currently, if a respondent leaves their house, and travels to several different stores by car, each store is entered as a separate record, as is their home starting location. What I would like to have is one "trip" for each time the respondent leaves their house and comes back, regardless of how many places they visit while preserving the information about each place they visit. I am wondering what the best way to tackle this problem would be. I have so far struck out finding prepackaged tools to perform such a task in QGIS, GRASS, ArcGIS, Maptitude, and PostGIS. Should I approach this as a scripting problem using SQL, scripting a shortest path tool in a GIS, or is there another more efficient way to do it. Many thanks in advance!


OSM XML: Should I generate contact:phone or phone?

I am generating OSM XML data for all POIs of Wikivoyage, the collaborative travel guide. When I want to express a phone number, should I generate `contact:phone` or `phone`? Both mean the same thing. According to a 2012 note on this page, the only difference is that `contact:phone` is newer but `phone` is more popular. Which one should I use, that so my data is readable by most 2013 apps?


how can the ogr2ogr help to merge a shp file?

I have over 30 shp files which is named "street.shp" into the sub-folder. There is over 30 sub-folders. In the sub-folder, there is another shp files which may be named a building.shp...etc. I only want to merge the street.shp files. The structure of all street.shp is exactly same! I tried this script -               for /R %f in (*.shp) do ogr2ogr -update -append merge.shp %f -f “esri shapefile” -nln merge       How can the script find the shp file named "street.shp" only? Thanks a lot.


Trigger a custom ArcToolbox tool outside Model Builder

I have created a model in a custom Arc Toolbox with Model Builder in ArcMap 10.1. Is it possible to trigger the model outside Model Builder and ArcMap? For example, use a Python script to trigger the model to operate?


ArcToolbox model run but warning output in Python

I have created an ArcToolbox and a Model in the toolbox. It uses a lot of Spatial Analyst tools and deals with raster images. When I imported the toolbox in a Python script and ran the model from Python in Windows command prompt, I got a few warning messages as below but the model got run successfully.               <string>:1: SyntaxWarning: import * only allowed at module level      It was shown several times. I am unable to tell where this error comes from. Any idea? Cheers


How to identify valleys on DEM in ArcGIS 10.1

First, What is the best and easiest way to identify valleys in a DEM raster? Second, how do you create a least-cost pathway for each valley in a region based on low angle slope and proximity to rivers? Any help would be massively appreciated!!!


SAGA raster calculator, ifelse command syntax problem?

I am trying to get the saga raster calculator in QGIS to spit out all values of a raster that are inferior to 10, and make all other values zero. I think the syntax is correct but I keep getting the original raster as an output which is definitely not what I want. Here is the code I am using:                   ifelse(lt(a,10),a,0)      I'm thinking there must be a problem with the way I am referencing the raster? I read in the SAGA help pdf that "if you have four grid map layers in your ‘Input’ list, they will be referred to by a, b, c, and d." So I'm not sure what I'm doing wrong.


Coordinate conversion (EPSG:3857 to EPSG:4326) using OpenGIS/JTS too slow

I am using the Java based libraries of OpenGIS and JTS to convert coordinates from EPSG:3857 to EPSG:4326 to proceed to compute the distance using the Haversine formula. The distance features in a constraint on whether to zoom in further on a map or not. We are using a tangible user interface and depend on the reactivity of the system. Unfortunately, the conversion is quite slow which induces even more delay than loading the map from a server. Would there be a less costly way to proceed with the conversion using local resources?                       CoordinateReferenceSystem targetCRS = CRS.decode("EPSG:4326");             CoordinateReferenceSystem sourceCRS = CRS.decode("EPSG:3857");             MathTransform transform = CRS.findMathTransform(sourceCRS, targetCRS, false);             GeometryFactory geometryFactory = new GeometryFactory(new PrecisionModel(), 3857);             Point minPoint = geometryFactory.createPoint(new Coordinate(xMin, yMin));             Point minTargetPoint = (Point) JTS.transform(minPoint, transform);             Point maxPoint = geometryFactory.createPoint(new Coordinate(xMax, yMax));             PointmaxTargetPoint = (Point) JTS.transform(maxPoint, transform);


Problem with precision when convert from WKT to shapefile

I have a problem with converting from WKT to shapefile in pygdal. In wkt I have fields with precision (for example real 7.2). When I create a shapefile I have field numeric with no precision (numeric 7.0). I can't find problem. Do you have any suggestions?


How to standardise colour-value graduations across multiple maps

In QGIS I have a vector dataset (imported csv vector + csvt) spanning multiple years, and a shape file shape which is essentially a grid. By joining the csv file with the shape file I can apply a colour to each grid square and hence a gradient scheme for the whole grid. I have created a series of maps - one map per year, each with its own colour gradient. However the colour gradients are not comparable between maps because the range of values is different among years. How can I make the gradient scales standardised across all maps? I've tried splitting the csv file into years, and with each year individually, it does not seem possible to create a 'custom' range of values for the gradient, say 0-100, even tho the values span say 20-70. I've also tried working with one csv file containing all years, thinking that I could somehow 'trick' QGIS into applying a gradient that spans the range of the whole dataset, while at the same time only map individual years. Categories doesn't seem to be an option as if I have 200 grid squares I get 200 categories each with own unique value. Any help would be much appreciated. It may be that there is a database-based solution which I am open to. I am also open to GRASS solutions, but I am less familiar with GRASS as this stage. Many thanks! Paul


Are there GIS tools for simulating landscape fragmentation?

I am currently looking for any GIS packages (or add-ons) that allow the simulation of landscape fragmentation on raster data. Maybe also other software than GIS tools? Does anyone of you have had experience with that issue already and can suggest a good tool for that task? Info on tutorials etc. are also welcome! ThanX! Jens


Map Refresh in Infrastructure Map Server 2015

I have a MapGuide application build in MapGuide 2011 where I refreshed my map using the following script               function MapRefresh() {            debugger;            ExecuteMapAction(20);             parent.parent.mapFrame.Refresh();              }          Now I am upgrading my application to Autodesk Infrastructure Map Server 2015 where the above script does not work, Please guide me as to how I can refresh my main map from a panel?


How to get elevation profile from .tif dsm on the web map?

I have my dsm .tiff file, pixels are floating point 32bit values giving the terrain elevation in metres. I need to create a web map where user can draw a path and get elevation profile for that path. Raster has attribute 'GRAY_INDEX' where values are stored. What are the options to achieve this and is it even possible? I thought about Geoserver `GetFeatureInfo`, but I'm not sure how that would work as I can get values for only one point at a time. I thought when user draws a path, I make a number of `GetFeatureInfo` requests exactly as points in that path, but I assume that would lead to slow performance. Is there a better option?


Custom style a particular road using CartoCSS/TileMill/OSM

For a particular city, I want to make the label for its main street labeled as prominently as Motorways. i.e. although this road would normally only be labeled from about zoom level 15 and up, I would like it to be visibly labeled from zoom level 12. Is there a syntax I can use in TileMill's CartoCSS to allow me to custom adjust label properties for this particular road? Ideally it could be based on road name (e.g. "Mary Street") but pretty much any solution could work such as using its ID number, because we only have to do it once.


GeoServer + ImagePyramid: Request is empty where data exists

I have set up an _image pyramid_ based on 5 levels of tiling. The data is served using the _ImagePyramid plugin_ for _GeoServer_. When I do a _layer preview_ the data is visualized just like I expect them to be. But when I zoom in a couple of levels, I get a _white_ tile in return. The _GeoServer log_ tells me the following: > 2012-05-14 08:07:07,495 DEBUG [gce.imagemosaic] - Reading mosaic from > file:/e:/GeoServer/data_dir/data/pyramid/2/2.shp   >  2012-05-14 08:07:07,495 DEBUG [gce.imagemosaic] - Highest res > 0.6400000000000005 0.6400000000000005   >  2012-05-14 08:07:07,495 DEBUG [gce.imagemosaic] - The computed > CropBoundingBox ReferencedEnvelope[534687.5 : 535156.25, 6225625.0 : > 6226093.75] Doesn't intersect the coverage BoundingBox > > ReferencedEnvelope[-1.0 : 0.0, -1.0 : 0.0] resulting in an empty request   >  2012-05-14 08:07:07,495 DEBUG [gce.imagemosaic] - RequestedBBox empty or > null   >  2012-05-14 08:07:07,495 DEBUG [gce.imagemosaic] - Request is empty: > RasterLayerRequest description:   >      RequestedBBox=ReferencedEnvelope[534687.5 : 535156.25, 6225625.0 : > 6226093.75]   >      RequestedRasterArea=java.awt.Rectangle[x=0,y=0,width=768,height=768]   >      RequestedGridToWorld=PARAM_MT["Affine",   >          PARAMETER["num_row", 3],   >          PARAMETER["num_col", 3],   >          PARAMETER["elt_0_0", 0.6103515625],   >          PARAMETER["elt_0_2", 534687.8051757812],   >          PARAMETER["elt_1_1", -0.6103515625],   >          PARAMETER["elt_1_2", 6226093.444824219]]   >      ReadType=JAI_IMAGEREAD   >  2012-05-14 08:07:07,495 DEBUG [gce.imagemosaic] - The response is empty. > ==> returning a null GridCoverage   > If I zoom in a couple of levels further, I get data again. The extent of the data is [530000, 6220000, 540000, 6230000] and the projection EPSG:25832, which the GeoServer admin interface agree on. **Then what could be the reason for GeoServer to look inside the BBOX [-1, 0, -1, 0], when my data is inside the requested map?**


osm2pgsql missing coordinates

I am quite new to GIS and OSM and am trying to get my head around exporting custom OSM region and uploading them into PostgreSQL. I have started with a very small map area to export with only two shapes - an island and a meadow - (http://www.openstreetmap.org/#map=19/43.86436/15.34266) which export the following .osm data:               <?xml version="1.0" encoding="UTF-8"?>     <osm version="0.6" generator="CGImap 0.3.3 (23630 thorn-02.openstreetmap.org)" copyright="OpenStreetMap and contributors" attribution="http://www.openstreetmap.org/copyright" license="http://opendatacommons.org/licenses/odbl/1-0/">         <bounds minlat="43.8635800" minlon="15.3411900" maxlat="43.8650000" maxlon="15.3445500"/>         <node id="63910084" visible="true" version="3" changeset="17709161" timestamp="2013-09-06T20:07:53Z" user="idelac" uid="183676" lat="43.8644519" lon="15.3432012"/>         <node id="63930603" visible="true" version="3" changeset="17709161" timestamp="2013-09-06T20:07:54Z" user="idelac" uid="183676" lat="43.8642369" lon="15.3431567"/>         <node id="63968158" visible="true" version="3" changeset="17709161" timestamp="2013-09-06T20:07:57Z" user="idelac" uid="183676" lat="43.8645365" lon="15.3424617"/>         <node id="1122490202" visible="true" version="2" changeset="17709161" timestamp="2013-09-06T20:07:44Z" user="idelac" uid="183676" lat="43.8640377" lon="15.3426171"/>         <node id="63879158" visible="true" version="3" changeset="17709161" timestamp="2013-09-06T20:07:52Z" user="idelac" uid="183676" lat="43.8639894" lon="15.3430113"/>         <node id="2448417703" visible="true" version="1" changeset="17709161" timestamp="2013-09-06T20:06:52Z" user="idelac" uid="183676" lat="43.8641440" lon="15.3427266"/>         <node id="2448417704" visible="true" version="1" changeset="17709161" timestamp="2013-09-06T20:06:52Z" user="idelac" uid="183676" lat="43.8641542" lon="15.3426165"/>         <node id="2448417706" visible="true" version="1" changeset="17709161" timestamp="2013-09-06T20:06:52Z" user="idelac" uid="183676" lat="43.8642112" lon="15.3428451"/>         <node id="2448417707" visible="true" version="1" changeset="17709161" timestamp="2013-09-06T20:06:52Z" user="idelac" uid="183676" lat="43.8642458" lon="15.3424979"/>         <node id="2448417708" visible="true" version="1" changeset="17709161" timestamp="2013-09-06T20:06:52Z" user="idelac" uid="183676" lat="43.8643394" lon="15.3428931"/>         <node id="2448417709" visible="true" version="1" changeset="17709161" timestamp="2013-09-06T20:06:52Z" user="idelac" uid="183676" lat="43.8644045" lon="15.3425008"/>         <node id="2448417711" visible="true" version="1" changeset="17709161" timestamp="2013-09-06T20:06:52Z" user="idelac" uid="183676" lat="43.8645022" lon="15.3428874"/>         <node id="2448417712" visible="true" version="1" changeset="17709161" timestamp="2013-09-06T20:06:52Z" user="idelac" uid="183676" lat="43.8645144" lon="15.3425403"/>         <node id="2448417713" visible="true" version="1" changeset="17709161" timestamp="2013-09-06T20:06:52Z" user="idelac" uid="183676" lat="43.8645937" lon="15.3427745"/>         <node id="2448417714" visible="true" version="1" changeset="17709161" timestamp="2013-09-06T20:06:52Z" user="idelac" uid="183676" lat="43.8645978" lon="15.3426673"/>         <node id="1122490153" visible="true" version="2" changeset="17709161" timestamp="2013-09-06T20:07:43Z" user="idelac" uid="183676" lat="43.8645802" lon="15.3431461"/>         <node id="1122490171" visible="true" version="2" changeset="17709161" timestamp="2013-09-06T20:07:43Z" user="idelac" uid="183676" lat="43.8646997" lon="15.3426923"/>         <node id="1122490322" visible="true" version="2" changeset="17709161" timestamp="2013-09-06T20:07:44Z" user="idelac" uid="183676" lat="43.8647101" lon="15.3429354"/>         <node id="1122490460" visible="true" version="2" changeset="17709161" timestamp="2013-09-06T20:07:45Z" user="idelac" uid="183676" lat="43.8641347" lon="15.3424630"/>         <node id="1122490524" visible="true" version="2" changeset="17709161" timestamp="2013-09-06T20:07:45Z" user="idelac" uid="183676" lat="43.8642836" lon="15.3424149"/>         <way id="8837893" visible="true" version="3" changeset="7098380" timestamp="2011-01-26T21:37:50Z" user="goranT" uid="217969">             <nd ref="63879158"/>             <nd ref="63930603"/>             <nd ref="63910084"/>             <nd ref="1122490153"/>             <nd ref="1122490322"/>             <nd ref="1122490171"/>             <nd ref="63968158"/>             <nd ref="1122490524"/>             <nd ref="1122490460"/>             <nd ref="1122490202"/>             <nd ref="63879158"/>             <tag k="natural" v="coastline"/>             <tag k="place" v="island"/>             <tag k="source" v="PGS"/>         </way>         <way id="236881851" visible="true" version="1" changeset="17709161" timestamp="2013-09-06T20:07:26Z" user="idelac" uid="183676">             <nd ref="2448417709"/>             <nd ref="2448417707"/>             <nd ref="2448417704"/>             <nd ref="2448417703"/>             <nd ref="2448417706"/>             <nd ref="2448417708"/>             <nd ref="2448417711"/>             <nd ref="2448417713"/>             <nd ref="2448417714"/>             <nd ref="2448417712"/>             <nd ref="2448417709"/>             <tag k="landuse" v="meadow"/>         </way>     </osm>      Running `osm2pgsql -S default.style -d gis -U postgres island.osm` uploads the data into my gis database with one row in the planet_osm_line table (for the island) and one row in the planet_osm_polygon table (for the meadow). Now, when issuing sql query `SELECT ST_AsText(way) FROM planet_osm_line`, pqsql returns               "LINESTRING(1707992.39 5444456.51,1707997.34 5444489.7,1707915.02 5444502.76,1707932.32 5444425.75)"      Here I don't understand two things- **Questions**   1. Why does the Linestring only contain four points when the original coastline of the island has 10 distinct nodes? Similarly when issuing the sql query on the planet_osm_polygon table to see the polygon for the meadow shape, the polygon has fewer coordinates than its original shape described by the nodes (7 instead of 11).   2. _(answered just below in update) How are the original lat,lng coordinates projected? I know osm2pgsql defaults to mercator projection. However, I don't understand how it computed the values in the Linestring above. Is the width of a tile (256px?) and some zoom (z=?) taken into account in the projection? Would LatLng[0,0] be projected as a Point[0,0]?_ **Update** The answer to my second question is that osm2pgsql projects from latitude and longitude coordinates to mercator in meters. LatLng[0,0] will indeed be projected as Point[0,0] where the point coordinates are the distance in meters from the origin on a mercator plane. Here is the projection function in case someone stumbles upon this question and looks for a solution               var degrees2meters = function(lon,lat) {             var x = lon * 20037508.34 / 180;             var y = Math.log(Math.tan((90 + lat) * Math.PI / 360)) / (Math.PI / 180);             y = y * 20037508.34 / 180;             return [x, y]     }      I still don't know why some of the nodes are missing in the resulting shapes in the database. Has anyone experienced this?


Viewing map tiles generated from Maperitive in leaflet or openlayers

for few days i have been working on developing my own osm tile server. I tried it on linux as well as on windows, no complete success till now. I Have to make tile server on both windows and linux. for windows i used Maperitive and succeeded in generating tiles and then i used this guys code               http://gis.stackexchange.com/questions/84697/viewing-map-tiles-from-maperitive-in-leaflet-or-openlayers      now the prob is tiles are not showing up Do i need to instal some Apache module for it ? My tiles are kept here               C:\Profiles\Maperitive\Tiles      Where as i am opening html file from server               http://localhost:8080/mwork/maperitive_test1/index_leaf.html      which is kept in here               C:\Profiles\ApachePHP\apache\www\mwork\maperitive_test1


OGR: How to best use GetFeature on a selection?

I think the following is very confusing               sel_str = "ND_14 = '108'"     l.SetAttributeFilter(sel_str)     l.GetFeatureCount()     >>> 1       In a for-loop this gives me               for f in l:        g = f.GetGeometryRef()        print g.GetEnvelope()     >>> (8.39746, 8.95788, 46.52757, 46.99331)      Whereas GetFeature() gives me               >>> feat0 = l.GetFeature(0)     >>> feat1 = l.GetFeature(1)      >>> g0 = feat0.GetGeometryRef()     >>> g1 = feat1.GetGeometryRef()     >>> g0.GetEnvelope()     (8.21206, 10.65409, 56.55045, 57.75181)     >>> g1.GetEnvelope()     (10.60014, 10.6172, 57.48933, 57.4979)      Apparently `GetFeature()` is based on the FID of the whole set, not the selection. But using a for-loop on a selection of length one is pretty ridiculous. How do I get the FID of my selected feature then? Or what are other ways of using `GetFeature()` with your current selection?


data range changes when warping

When the source band having pixel values in the range of 52 and 207 is warped using the following command:               gdal_warp -t_srs wgs84 -r near in_file out_file      The resulted band has the pixel values in the range of 0 and 207. How can I make the data range unchanged during warp?


QGIS 2.4: Basemap.at is displayed in poor quality

I have a huge problem displaying the austrian Basemap.at WMTS Service in QGIS 2.4. When I add the service, qgis logs multiple errors on failed tile requests. When I zoom into the map, it is displayed in a very poor quality (see screenshot).   I really don't know how to fix this problem - can anyone tell me how to increase the quality of the tiles? WMTS URL: http://maps.wien.gv.at/basemap/1.0.0/WMTSCapabilities.xml ![Basemap.at Scrennshot: Poor Quality](http://i.stack.imgur.com/gDW9M.png) UPDATE 1:   I tried setting the resampling methods in QGIS to Cubic/Average, as you guys advised me to and compared the result with the same zoomlevel in ArcGIS. In this screenshot you can see the result: At the same scale level QGIS just skips an entire tile level, whereas arcgis loads this level correctly. - Any other ideas? ![enter image description here](http://i.stack.imgur.com/YJJMs.jpg)


Does GeoServer's WorldImage plugin use UTM or WGS84 (EPSG:4326)?

I have some PNG images that I want to put on an OpenLayers map using GeoServer as a WMS, and it seems as though the easiest way to do that would be to create an accompanying `.pnw` file (ESRI world file) for each `.png` and load them using GeoServer's WorldImage plugin. I also have the longitude, latitude coordinates of nine points on each rectangular image: {top, middle, bottom} x {left, middle, right}, and I can use five of them to convert pixel x,y into longitude, latitude x',y': ![x' = \\left\(\\frac{mr_{lng} - ml_{lng}}{width}\\right\) x + \\left\(\\frac{tm_{lng} - bm_{lng}}{height}\\right\) y + \\left\(mm_{lng} + \\frac{mr_{lng} - ml_{lng}}{2} + \\frac{tm_{lng} - bm_{lng}}{2}\\right\)](https://chart.googleapis.com/chart?cht=tx&chl=x%27%20%3D%20%5Cleft\(%5Cfrac%7Bmr_%7Blng%7D%20-%20ml_%7Blng%7D%7D%7Bwidth%7D%5Cright\)%20x%20%2B%20%5Cleft\(%5Cfrac%7Btm_%7Blng%7D%20-%20bm_%7Blng%7D%7D%7Bheight%7D%5Cright\)%20y%20%2B%20%5Cleft\(mm_%7Blng%7D%20%2B%20%5Cfrac%7Bmr_%7Blng%7D%20-%20ml_%7Blng%7D%7D%7B2%7D%20%2B%20%5Cfrac%7Btm_%7Blng%7D%20-%20bm_%7Blng%7D%7D%7B2%7D%5Cright\)) ![y' = \\left\(\\frac{mr_{lat} - ml_{lat}}{width}\\right\) x + \\left\(\\frac{tm_{lat} - bm_{lat}}{height}\\right\) y + \\left\(mm_{lat} + \\frac{mr_{lat} - ml_{lat}}{2} + \\frac{tm_{lat} - bm_{lat}}{2}\\right\)](https://chart.googleapis.com/chart?cht=tx&chl=y%27%20%3D%20%5Cleft\(%5Cfrac%7Bmr_%7Blat%7D%20-%20ml_%7Blat%7D%7D%7Bwidth%7D%5Cright\)%20x%20%2B%20%5Cleft\(%5Cfrac%7Btm_%7Blat%7D%20-%20bm_%7Blat%7D%7D%7Bheight%7D%5Cright\)%20y%20%2B%20%5Cleft\(mm_%7Blat%7D%20%2B%20%5Cfrac%7Bmr_%7Blat%7D%20-%20ml_%7Blat%7D%7D%7B2%7D%20%2B%20%5Cfrac%7Btm_%7Blat%7D%20-%20bm_%7Blat%7D%7D%7B2%7D%5Cright\)) where (mm_lng, mm_lat) is the position of the middle-middle point, (mr_lng, mr_lat) is middle-right, (tm_lng, tm_lat) is top-middle, etc. for the middles of the four sides. This is a linear approximation to the full coordinate transformation, minimizing error at the center of the image. The World Image format describes an affine transformation in the same way: the six parenthesized terms above can be interpreted as the six lines of the `.pnw` file (in A, B, C, D, E, F order) if the map units are longitude, latitude coordinates (that is, WGS84 a.k.a. EPSG:4326). The World file Wikipedia page says "World files describing a map on the Universal Transverse Mercator coordinate system (UTM) use these conventions..." and goes on to describe the UTM application, though it also says, "World files do not specify a coordinate system; this information is generally stored somewhere else." If the result of the `.pnw` conversion is interpreted as UTM, then I would need to do additional work to transform from WGS84 into UTM, and GeoServer would also need more information to specify which of the 60 UTM zones that the `.pnw` converts into. I don't see how to provide that information to GeoServer. Am I right in assuming that GeoServer interprets the `.pnw`-converted coordinates as WGS84?


How to render a static map from CartoDB?

CartoDB is working great for dynamic maps, but I have a use case to render a totally static map, just an image online. Google offers its static map api. Afaik there is no similar API for CartoDB. But does anyone have a work around that would let me display a single, non-zooming non-panning image? It'd be ideal if I could just put it in an html img tag, but something with javascript could potentially work. It'd be ideal if it's lightweight, as many of these users have very slow connections.


CartoDB data not displaying on Google basemap

I'm new to CartoDB but hoping to use it to display various types of real estate data on my website. That said I've been playing around with some of the sample code (Google Maps V3 integration example - https://github.com/CartoDB/cartodb.js/blob/develop/examples/gmaps.html) and cannot get this data to display within my website (see test page - https://metropolisiq.com/MapTest). I've tried the same code within Dreamweaver and it works fine. Not sure what the issue is, but perhaps it relates to the differences between Dreamweaver and my website:   1. My site is encrypted (HTTPS)   2. My site was built using the Joomla CMS. Any thoughts on how to solve are greatly appreciated.


ArcGISDynamicMapServiceLayer is not refreshing properly

I have a DynamicMapServiceLayer that I am refreshing every 30 seconds. I want the data in the map to refresh every 30 seconds too but it doesn't appear to be working properly. The map refreshes but then the service goes away and I just see a blank map. The service will stay off for about 20 seconds then come back on again. My code is here too...Is there a way to keep the service on all the time and still refresh it or make it so it doesn't go away? Thanks in advance!               var map;         require(["esri/map", "dojo/domReady!", "esri/layers/ArcGISDynamicMapServiceLayer"], function (Map, domReady, ArcGISDynamicMapServiceLayer) {             map = new Map("mapDiv", {                 center: [-56.049, 38.485],                 zoom: 2,                 basemap: "streets"             });             var liveLayer = new ArcGISDynamicMapServiceLayer("http://server/MapServer",              {                  refreshInterval: .5,                 disableClientCaching: true             });             console.log(liveLayer);                  map.addLayers([liveLayer]);         });


Selecting features through a many-to-one relationship

I want to select multiple counties based on the occurrence of a species present within them. I have a County polygon layer with CountyName as an attribute, a flat table called tblSpCnty with the attributes SpeciesName and CountyName (multiple CountyNames for each SpeciesName, representing all the counties in which each species occurs), and a flat table called tblSpeciesList with the list of species (SpeciesName, one record for each species). What I'd like to do is select a single species from tblSpeciesList, resulting in the selection of all the county polygons in which the species occurs.   I used to do this by linking tblSpeciesList through tblSpCnty to the County layer in ArcView, but I haven't been able to replicate the functionality in ArcMap. Using relates doesn't seem to achieve the same result. I'm not tied to ESRI, in fact if there's an Open Source solution all the better.


What are lazy geometries

What are "Lazy Geometries"? How are they different from other types of geometries? I'm an experienced Python/Django developer but am very new to the world of GIS, and understand that this is a newbie question. I've been unable to find the answer on google (if anyone can direct me towards a good tutorial site, I'd be grateful). I encountered the term when reading up on GeoDjango (and it's use of the GEOS library). Thanks in advance.


How can I speed up my arcpy script?

I have written an application that does a lot of geoprocessing using arcpy. Once started, the application itself runs at an acceptable speed, but it is very slow to boot. It can take 30 seconds or so to start up the GUI, so I'd like to find ways to make it start faster. As far as I can tell, the main source of delay is when I set the dataframe as follows:               self.df = arcpy.mapping.ListDataFrames(self.mxd, "")[0]      Does anyone have any tips / tricks to make either the whole application or the ListDataFrames run quicker when the application is started? I use ArcGIS 10.1, python 2.7 and tkinter for the GUI. Feel free to ask if I've forgotten to mention anything. EDIT: As requested, here is the (very!) stripped down version of the start my code:               import arcpy     import sys     from Tkinter import *     import ttk     import tkMessageBox     import tkFileDialog     import webbrowser     import time     import csv     import os     import textwrap          arcpy.env.workspace = r".\..\Append\temp.gdb"          Title = "Automated Map Generator"          class Application(Frame):         def __init__(self, master=None):             Frame.__init__(self, master, background = "white")             self.grid()                  self.createWidgets(master)              def createWidgets(self, master):                  # setup map document             self.mxd = arcpy.mapping.MapDocument(r".\autolim_mapping.mxd")             arcpy.env.overwriteOutput = True             self.df = arcpy.mapping.ListDataFrames(self.mxd, "")[0]          #[Lots of missing code here]          root = Tk()     root.tkraise()     app = Application(master=root)     app.master.title(Title)     app.mainloop()      How do I know what it's the ListDataFrames that takes so long? I put a whole heap of print statements through the startup code and noted how long they took to come out... I bet there are better ways to do this, but I'm pretty new to python. EDIT 2: @dassouki - thanks! Here are the output times from after 'import time':               Start Class Application (Frame):  1.53200006485     Start def Create Widgets:         1.59400010109     Start  mapping.MapDocument:       1.6099998951     Start env.overwrite:              1.64100003242     Start mapping.ListDataFrames:     1.64100003242     Finish mapping.ListDataFrames:    12.2660000324


interpolate point values to line

This might be a stupid question and I apologize if so. I have count data for different trails. They are points and I would like to make an average of those counts for each trail so that each trail has a count value. I don't know if I am being clear... Later on I would like to have a probability of encoutering the different species that were counted along those trails depending on the average counts. I am fairly new to GIS and I have no idea which tool I could use to do that... I was thinking of interpolation but I want to model the trails only, not the whole layer. I hope I have been clear enough and you will not mind my ignorance, Thank you,


Using Value from Iterate Feature Selection to Calculate Field in ArcMap 10.1 ModelBuilder

I have two polygon feature classes, an alignment layer and a parcels layer. I am trying to use ModelBuilder to create a tool that will do the follow:   1. iterate through groups in a single field of a polygon alignment layer (ie select and iterate through certain attributes in a field with same names such as "tunnel", "portal", etc)   2. Clip a parcel layer by each group    3. Populate a new field in the parcel layer called "EsmtType" with the name of each group (group is called "Value" in model)   4. Merge the clipped groups into one feature class I have gotten steps 1 and 2 to function but I am stuck on step 3. I am not sure how I need to configure the "Calculate Field" tool in order to recognize the "Value" variant of the "Iterate Feature Selection" tool. I tried using the advice from this post Using a value from "iterate feature selection' to calculate a field but I couldn't figure out how to implement it for my situation. Here is the model I have thus far: http://imgur.com/ji2nAJJ I would also like advice on how to implement step 4, I am not sure if it would require integrating a model within a model. Thank you so much!!! -kkaszas


Part of text missing in legend and text boxes

I have had issues lately when exporting to pdf from ArcMap 10.2. Parts of the text sometimes fall off and goes missing. Mostly it's been some lines in legends, but now also text boxes. The text is also missing in legends displayed in Arcmap, but not always. ![enter image description here](http://i.stack.imgur.com/QkEUo.png) This is not how the text in the textbox supposed to look, obviously. ![enter image description here](http://i.stack.imgur.com/Vl9p7.png) This screenshot shows a legend, but some of the descriptions are missing. In Arcmap the legend sometimes looks like this: ![enter image description here](http://i.stack.imgur.com/MDoJD.png) When it's supposed to look like this: ![enter image description here](http://i.stack.imgur.com/axc9X.png) It happens that this gets solved when exporting a second (or third) time. Which is good, of course, but does not help me in trying to figure out why it's happening in the first place. Has anyone else encountered this problem? Is there a solution? This is really starting to get on my nerves, and not even creating a new legend helps every time (which is very tideous in the first place).


What if format for Getcapabilities document

I want to add capabilities URL for cascading a layer from remote WMS in Geoserver. What will be the exact format for Getcapabilities URL ?


Pick object by properties (and centre on it)

I have a Leaflet map with a Cloudmade base layer and a single layer of polygons. Each of these polygons has a reference number attached as a property. I would like to select a reference number from a list (or whatever) and then pick out the object by its reference number and zoom to the extents of that object. How can I pick an object by its properties? (Presumably once I have done that zooming to extents shouldn't be too hard)


reproject wms in geoserver

> **Possible Duplicate:**   >  Re-projection in Geoserver I've got a WMS service that publishes in WGS84 (EPSG:4326). I have several AutoDesk users that need to use the imagery in NAD83 (EPSG:2227). I have geoserver 2.1 setup to cascade the wms and reproject it. Problem is that as soon as I turn on SRS handling to reproject to declared SRS the reprojected wms fails. What am I doing wrong? Is Geoserver capable of doing the reprojection.


Attribute Table Decimal Places

As we all know, the attribute table field properties has number format to display number of decimal places to up to maximum 15. If we have many fields/columns in attribute table, how can we change all fields/columns decimal places roundoff/display at a time ( I mean in single operation. Not changing each and every field properties). Ex:               ColA:12.987654 ! ColB:34.20987766  ! ColC:89.654309877      In one operation my output should be (Rounding/display to 2 digits) For Rounding:               ColA:12.99 ! ColB:34.201  ! ColC:89.66      For Display:               ColA:12.98 ! ColB:34.20  ! ColC:89.65      is there any method/tool to do this in Arcgis?


Fix topology error tool- can't automatically fix multiple errors?

I am a ArcGIS.10. user. I'm want to solve multiple errors using the fix topology error. So I select all my features I click right but I don't have the correcting options (Snap, trim, extent), I only have Delete. Does anyone know why? Thank you very much,


Can someone help me to convert this Python to VBscript for labellling

I have two Python script snippets that need to convert to VBscript. Can anybody help me because I am really not familiar and know Python. the first script is               def FindLabel([BA2]):         name = str([BA2])         if "-" in name:             first, second = name.split("-",1)             for letter in first:                 if not letter.isalnum():                     first = first.strip(letter)             return first         else:         return name      the second script is:               def FindLabel([BA2]):         name = str([BA2])         if "-" in name:             first, second = name.split("-",1)             for letter in second:                 if not letter.isalnum():                     second = second.strip(letter)             return second            else:             return name


Problem with snap layer when in name is space

I'm trying to make a option, where user can choose layer to snap. All is okey when layers names are without space, but if in name is space (like 'something layer name') snapping does not work. How can I solve it? thank for all posts.


Pause ArcPy script until feature class lock goes away

My script needs to copy data into a feature class to an SDE geodatabase that might have a lock on it from another application. I shouldn't break the lock (that app has priority over a data update), but I also don't want the script to quit just because it tried at the wrong moment. How can I set it to keep trying until there is no lock? I'm looking for something like a try/except, but more like try/keep-trying-every-15-seconds. I considered making a separate function and adding a pause:               try:         funcCopyData()     except:         pause 10 seconds         funcCopyData()      but wouldn't this still be limit me to only two attempts?


QGIS : launch Atlas composer generation from command line

I am using QGIS 2.2 and have configured a composer with atlas generation. I would like to launch this atlas generation every day. The idea is to write a batch and launch it with a cron. Is it possible to launch atlas generation from command line ? Thanks


Export Stretch Function values - Arcpy

I'm looking for a way to export Stretch Function values from images that’s stored in a Mosaic dataset. Not values from the whole set, but specific stretch functions per image. The set is stored in either a filegeodatabase or SDE. The values I’m looking for is eg. Band[] Min, Max, Mean values and Gamma values.... EditRasterFunction_management() allows me to insert, update or delete which I’ve used to initially create the Stretch Function. But after some manual editing for the raster images I need a way to export the new values. I can’t manage to get access to this information using cursors since it’s a function? Any ideas how to maintain this function values? Preferably using Arcpy/python. Thanks!


arcgis silverlight - how to delete/cut a part of a multipart polygon

Say I cloned a single-part lake polygon that had a hole for an island. Then I used the Editor Reshape command to cut a piece off. That incidentally caused a second part to be created where the island was, as seen here. ![enter image description here](http://i.stack.imgur.com/v1J0v.png) Now I want to get rid of that bogus island. I've tried drawing a circle around it with the Cut command, tried the Reshape command again, no joy. Help a dummy out? (yeah I know my links are to ArcGIS Desktop editing tools, the better descriptions of the apparently more or less equivalent Silverlight API Editor commands) ;-)


Missing tiles when django application is hosted on an apache server

I'm building a GIS application using Django as my web framework and OpenLayers to implement a slippy map interface. I was initially using the django runserver for development and everything was working fine. However, I recently ported my app on an Apache server and I started seeing some missing tiles (see here): It's not always the same tile that is missing and usually there are two or three missing tiles. Sometimes, after refreshing the page several times the tiles appear properly. There are no missing tiles when I zoom in or out. Note that the page keeps on loading for an indefinite amount of time when I experience this issue. I'm using a TMS layer to display the map using data from the GSHHS database.The map is being rendered via a mapnik XML file which uses the GSHHS shapefile. Unfortunately, at this point I'm unable to paste any code over here since I'm unable to login to SE from work. I'll try to post some code later. I have no idea as to what might be causing this problem. There is nothing helpful in the Apache and Django log files. Any help would be greatly appreciated. Thanks in advance.


Problem with MIF/MID allignment

After creating a rather large map in MapInfo using US Census Tiger data and exporting it to the MIF/MID format, it appears the files are no longer properly aligned. When I try to open the MIF/MID in some other software it gives me an error. "Got Pline in the middle of reading other Pline" I assume this means that the MIF formatting is incorrect. I have recreated the MIF/MID files and the problem continues. Maybe a problem in the US census data? I am quite experienced with doing this process and I have never had this problem before. Help!


Why Add Join does not work in ArcMap and Arcpy?

I have problem with joining feature class with table from the same location (GDB) using Python script. When I use join Join and Relates -> Join... in ArcMap TOC it does not work because I do not get any joined records despite the fact there are equal values in fields used as a keys . The same problem occurs using arcpy. AddJoin_management(). When I use arcpy.JoinField_management () or the same tool from Toolbox in ArcMap I get the results. What is the reason that JoinField works and it gives me a good results, but in the case of AddJoin I do not get any joins? The field type in the Feature Class is the same as in the table to join. Any suggestions?


Find tunnel 'center line'?

I have some map files consisting of 'polylines' (each line is just a list of vertices) representing tunnels, and I want to try and find the tunnel 'center line' (shown, roughly, in red below). ![alt text](http://i.stack.imgur.com/jZdYb.png) I've had some success in the past using Delaunay triangulation but I'd like to avoid that method as it does not (in general) allow for easy/frequent modification of my map data. Any ideas on how I might be able to do this? Thanks. P.S. I'm working in fairly raw C++, and was referred here from Stack Overflow.


What are the most tangible benefits to getting your GISP or ESRI Certs?

> **Possible Duplicate:**   >  GISP - Value of being certified? Looks like I'm only a few presentations or other "contributions to the Profession" points away from being able to apply for my GISP. Has anybody found any actual benefits to getting the certification? Seems like it might be somewhat arbitrary. Of course it lets people know you have a high level understanding of GIS but beyond that are there any benefits? Salary bumps? Other? **ADDED** Additionally, what, if any are the benefits to getting the new ESRI Certifications? ESRI is obviously the most well known GIS software, but there are others... so if one is looking to move forward in their GIS career, what is the best way to do it as far as certs/continuing education go?


Point in polygon algorithm for multiple polygons

I've a Google map with bunch of polygons on it. Here is a problem I'm interested in: Given a lat,lng point, what's the best way to determine all the polygons that this point lies in? The obvious way is to run a "point in polygon" algorithm iteratively for each polygon, but I was wondering if there an efficient algorithm to answer such queries especially if you have thousands of polygons.


Is there a way to add turn restrictions in A* and Dijkstra?

We are currently working with pgRouting and found that there's no way of implementing turn restrictions (e.g. banned right or left turns). While it is possible to assign "to_cost" and "rule" in Shooting * algorithm... I couldn't find a way of implementing these restrictions using A star and Dijktra algorithms. Is there a way of implementing special turn restrictions in A star and Dijkstra methods ?


Command r.sunmask in GRASS

I have hypsometric model of relief, and I need to visualize shadows. I was suggested to use GRASS command r.sunmask, so I tried to set time values and run command, but it gives me an error               "ERROR: Neither sun position nor east/north, date/time/timezone definition are complete"      So then I tried to put values in Position tab, which I founded by SOLPOS Calculator (http://www.nrel.gov/midc/solpos/solpos.html)               r.sunmask elevation=model@PLANINA output=model_sunmask year=2011 month=1 day=15 hour=9 timezone=1 east=45.909361 north=15.965538      But it didn't work out, too, it says "Calculating shadows from DEM..." but nothing's going on for hours of waiting. So I ask if that command works at all in GRASS, or I'm doing something wrong? You can suggest some other application for same thing, except ArcGIS applications. Thank you.


Save spatial reference to application settings?

Having created a tool as a macro for several geoprocessing tasks, I'm saving user input to application settings so the user doesn't have to enter all required data on each use. I'm stuck on saving the spatial reference. ESRI code shows using the FactoryCode and an enumeration to create either a projected or a geometry spatial reference. What I'm trying to figure out is how to get the enumeration, and how to identify whether the factorycode/enumeration points to a projected or geometry spatial reference. I will greatly appreciate any answers/directions provided.


How to set the value of a field in all selected rows of a table using ArcGIS 10

I have an attribute table in ArcGIS (it's large - around 300,000 elements). I've selected some rows of the table using 'Select by attributes' and I want to set the value of a field for all of those rows. How should I go about doing this? I'm hoping there will be a nice, easy, GUI- based way - as I'm teaching someone with very little GIS experience.


What is future of Silverlight with ArcGIS Platform?

At the September 2011 Build Conference, Microsoft said a lot about Windows 8 and the Future of JavaScript over Silverlight. No Plugins wanted!! JavaScript with html5 ftw! What is the future of Silverlight with the ArcGIS Platform (ArcGIS for Server, ArcGIS Online, etc) in 2014 and beyond?


cad and geoserver

Hi Guys i have geoserver and openalyers running quite nicely and now need to overlay a cad drawing. I have looked at the geographic location part of cad and set the long / lat to the appropriate values. I saved the drawing and then converted to a shapefile and then imported in to postgis. When i add to geoserver the nativebounds and lon/lat generated in geoserver are no where near where i would expect them to be. Can anyone point out what i am doing wrong. When i set the long/lat in my drawing where does this point reference. is it the center, top and left, bottom and left of the drawing? Also how can i find out what projection the drawing is in. It doesnt seem that you can set this in cad. Ive asked my cad guys and no one seems to know. Thanks in advance


Is there any feature limit for Spatial data migration to Postgre SQL9.2...?

I am trying to migrate spatial data to Postgre SQL9.2 using 'Postgis Shapefile & DBF Loader 2.0'....Each time it returns a "Runtime error'. Primary data extension - .shp ; Data type - Multipolygon ; Feature Count - 1116930 . Is this a problem due to the lage amount of data involved....!!! Thanks in advance..


In the definition of **surface plots** in my book, (An introduction to data mining, chapter 3) it is said that: > Surface plots are often used to describe mathematical functions or physical > surfaces that vary in a relatively smooth manner. I don't understand what exactly is meant by the above sentence. What special property do mathematical functions and physical surfaces have that is appropriate for the case of surface plot?


maxScale or minScale dont working in vector (from wfs) layer

var options =          {             projection: "EPSG:900913",             displayProjection: "EPSG:4326",             units: 'm'         };         map = new OpenLayers.Map('map', options);                    ilceSinir = new OpenLayers.Layer.Vector("İlçe Sınırı",         {             minScale: 150000,             strategies: [new OpenLayers.Strategy.Fixed()],              projection: new OpenLayers.Projection("EPSG:4326"),             protocol: new OpenLayers.Protocol.WFS(             {                 version: "1.1.0",                 url: "http://localhost:8080/geoserver/wfs",                 featurePrefix: 'topp',                 featureType: "g_ilce_siniri",                  featureNS: "http://www.openplans.org/topp",                 geometryName: "bounds"              })         });      scale isn't working. layer is visible everyzoom.


Heat Map trouble: Is this a data formating problem?

I'm playing with the heat map plugin with QGIS. I've run through QGIS tutorial: http://docs.qgis.org/html/en/docs/user_manual/plugins/plugins_heatmap.html,   using the airports data set as explained and I reproduced tutorial results. I have a 600,000+ data set (five years worth) of police calls from Philadelphia, which I've reduces to 185 points (a day) for testing. The Meta data for both sets is below although I've used both the dbf file and saved shp file as input. For larger Radius settings (above about 0.1) I get a rainbow bulls Eye centered on but many time lager than my data set range. For smaller Radius settings I get a uniform single color covering the data set range and transparent holes opening over (but not centered on) each data point. I get just what I should with the Alaska data set. General:               Storage type of this layer: ESRI Shapefile     Source for this layer: C:\Users\Gallery of Tools\Documents\GIS DataBase\Alaska\shapefiles\airports.shp     Geometry type of the features in this layer: Point     The number of features in this layer: 76     Editing capabilities of this layer: Add Features, Delete Features, Change Attribute Values, Add Attributes, Delete Attributes, Create Spatial Index, Fast Access to Features at ID, Change Geometries     Extents:     In layer spatial reference system units : xMin,yMin -4480198.52,1433525.80 : xMax,yMax 4615124.98,6502586.83     Layer Spatial Reference System:     +proj=lcc +lat_1=40.96666666666667 +lat_2=39.93333333333333 +lat_0=39.33333333333334 +lon_0=-77.75 +x_0=600000 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=us-ft +no_defs               General:     Storage type of this layer: ESRI Shapefile     Source for this layer: C:\Users\Gallery of Tools\Recovered from the Dell600GB\Image Catalog\EngravingsNMaps\Philadelphia\Police Incident\PoliceIncJune302012.dbf     Geometry type of the features in this layer: Point     The number of features in this layer: 185     Editing capabilities of this layer: Add Features, Delete Features, Change Attribute Values, Add Attributes, Delete Attributes, Create Spatial Index, Fast Access to Features at ID, Change Geometries     Extents:     In layer spatial reference system units : xMin,yMin -75.2515,39.9047 : xMax,yMax -74.9594,40.1226     Layer Spatial Reference System:     +proj=lcc +lat_1=40.96666666666667 +lat_2=39.93333333333333 +lat_0=39.33333333333334 +lon_0=-77.75 +x_0=600000 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=us-ft +no_defs      P.S. And another thing. When I change the Radius setting, the cell size changes automatically. Is this supposed to happen?


How to intersect 2 polygon layers with a given tolerance for features in different layers only?

I can intersect 2 layers with polygons using ArcMap. However I can specify only overall tolerance no matter what layer a feature is in (even same layer). So I either end up with lots of slivers if I don't tolerate snapping, or I end up collapsing small features from the first layer that I'd want to keep. The problem is that I would like for features let's say less than 10 m apart in the first layer to stay as they are if they happen to be contained within some polygon of the second layer. However I'd like for edges of polygons of the 2nd layer to snap to those in the first layer if they are within 10 m. I'm open to a solution using any tool be it ArcMap or some python library/script. Here is the rough idea of what I want. ![enter image description here](http://i.stack.imgur.com/0XHhb.png)


Tools to calculate routing distance

I'm working to a web application, and I need to calculate the routing distance (not the distance "as the crow flies") between a starting address and several final addresses. I don't need to render the routes on the map. The application will have a lot of queries so I can't use Google APIs and probably the best solution is to perform the task not with an external service but on a local server. I think the only way is to use Open Street Maps (or not?). Which tool to perform the calculation do you recommend? Thank you in advance! * _*_ edit:   * At the moment geocoding (700 conversions daily) is performed by GoogleMaps and Nokia's HERE APIs. Do you recommend me some other tools also for this?   * The number of daily point to point routing calculations are approximately 100.000


ArcToolbox 9.3 converted to ArcToolbox 10 and now won't work in ArcMap 9.3

I added a custom toolbox created for 9.3.1 to a ArcGIS Desktop 10 ArcMap. Now it seemed to have disappeared. I can't add it to ArcMap 9.3.1. Is there a way to convert it back if that is even the issue? Or is there something else going on? Thanks.


grass gis- ERROR: nwalk (7000001) > maxw (7000000)!

i want to make r.sim.water (hydrological modelling) in Grass Gis, and this ERROR: nwalk (7000001) > maxw (7000000)! always appears, any ideas how to solve it? number of cells of dem is below 7000000.


Geoserver administration - defining different privileges

I need to create separate admin accounts for Geoserver. I'd like to have one superadmin account that would be able to read/edit/delete any style, layergroup, workspace, store or layer - that is what the default admin is already able to. I also need to have several accounts that would be able to create workspaces, layers and stores, would not be able to edit styles and would be able to edit and delete only the workspaces, layers and stores they created (e.g. if superadmin creates a workspace, he is the only one to delete it). Is this possible to achieve with Geoserver?


Convert Shape Polyline-M to Polyline

I have some pretty old software that I am forced to use (no option). It's old enough that it doesn't recognize shape files of type Polyline-M, it only recognizes shape files of Polyline. I have a shape file that is completely Polyline-M that I need to convert to plain old Polyline. I'm struggling to find a way to do this in ArcMap 10. Does anyone know how to do this? A helpful link, or even just a point in the right direction to where I could read how to do it would be really handy. I'd rather not have to create a tool to do it, surely something like this exists out there already! :) Thank you


Getting what seems to be incorrect data from GeoNames

I have a hiking website. I am trying to query for parks big and small around the world where people could potentially go walking or hiking. I was trying to use the GeoNames data and was using query arguments from this table: http://www.geonames.org/export/codes.html Here is an example query I made trying to get parks:               http://api.geonames.org/search?maxRows=10&username=demo&q=PRK      and it gives me some things which are I have no idea what. Here is also another thing where I try to query for arctic spaces:               http://api.geonames.org/search?maxRows=10&username=demo&q=LAND       And that gave places in the Carribean :) So I must be doing something totally wrong. My goal was to look up all parks in the world, and also cities (with state/country info) so that I can make that connection on my site for my users. Is my goal possible to achieve? What was I doing wrong in my queries? Thank you!


Is there a tool for automated quality assurance for vector data captured on raster data?

I wonder whether is any tool(or script) in ArcGIS or QGIS which allows to make an automated quality assurance for vector data captured on raster data? vector data on binary raster: ![enter image description here](http://i.stack.imgur.com/ewmav.jpg)


Import KML attributes into the correct schema

I have a kml file containing over 4000 polygons with attributes which must be retained (in 3 different fields). After using the standard `KML to Layer` tool, all the attributes are now in one `PopupInfo` field, along with the html code. I have the `Data Interoperability` extension, so after looking at the questions here and here, I came across a post on the ArcGIS forums: > The Data Interoperability extension will read KML (including folders and > subfolders) and use any schema object defining the attributes. Thus what it > writes out will have the correct field definitions. The same goes for > writing KML with a schema object, given geodatabase inputs. What I understand from that is that using the `Quick Import` tool will read the KML and export it to a gdb based on a schema which I defined. This seems to be exactly what I want, except I cannot find the option to specify the schema. I exposed the `kml_name` attribute as mentioned in the other question, but the field remained empty. I tried downloading the KML2shapefile plugin with MapWindow, but it is only for 32 bit and I am running 64 bit. I am unable to download the FME trial at work.


Background (sky/terrain) of a viewshed

Is it possible to determine whether a point feature of a particular height would be backed by sky or by terrain? I've done the viewsheed analysis of a wind turbine. Now I would like to make a raster map representing the background (sky/terrain) from each cell of the viewshed.


Nodata issues in map algebra

A dummy question. I need to change all values in my rasters into ones. I've been able to do this fine in the Rater Calculator just by multiplying the layer by zero and adding one (my_layer*0+1), but when I'm doing the same in arcpy, for some reason it works for a part of my raster but for some it changes all the nodata cells to ones as well. The nodata values in different rasters are different because my rasters have been generated from different source data. Does it not recognize some nodata values? If so, this is odd as it displays them correctly when I just import them on the mxd. Any ideas on why this works sometimes, but sometimes not? Here's the exact script:               bin_test = arcpy.Raster("F:/datasets/map1.img")*0+1


can't generate/display signed float values using python/gdal

I'm trying to create a raster with signed float values as the pixel values and display the results in MapWindow on windows. What is the proper gdal band type to use? I've tried both **gdal.GDT_Byte** and **gdal.GDT_Float32** and using both methods I can't seem to get the pixels displayed. If I change the pixel value to the abs() value they will display.               dataset = driver.Create(filename, cols, rows, number_of_bands, gdal.GDT_Byte)      **gdal.GDT_CFloat32** allows for negative numbers but doesn't properly display the floating point value past the point.


Integrate QGIS with ASP.NET MVC3

I'm completely new to GIS and i'm trying to work my way into this. My goal is to integrate a simple map created using QGIS with my asp.net web mvc3 application. The main objectives would be to show the map on the application as well as performing some operations with the different points and lines (points of interest and roads respectively). Can anyone point me in the right direction on how to achieve this?


Does the esri javascript api Identify Task Return a FeatureSet?

I am trying to get the results of an identify task using ArcGIS JS API 2.7 into a featureSet so I can utilize this sample http://arcscripts.esri.com/details.asp?dbid=16528 that takes the results of a datagrid store and exports them to an excel file. I have it working where one feature is being exported, I just need to have all selected features export. Does an Identify Task produce a FeatureSet? Thanks!                function executeIdentifyTask(geom) {       //clear the graphics layer        map.graphics.clear();            identifyParams.geometry = geom;       identifyParams.mapExtent = map.extent;       identifyTask.execute(identifyParams,function(response){           var polygonSymbol = new esri.symbol.SimpleFillSymbol(esri.symbol.SimpleFillSymbol.STYLE_SOLID, new esri.symbol.SimpleLineSymbol(esri.symbol.SimpleLineSymbol.STYLE_DOT, new dojo.Color([151, 249,0,.80]), 3), new dojo.Color([151, 249, 0, 0.45]));          var markerSymbol = new esri.symbol.SimpleMarkerSymbol(esri.symbol.SimpleMarkerSymbol.STYLE_CIRCLE, 20, new esri.symbol.SimpleLineSymbol(esri.symbol.SimpleLineSymbol.STYLE_SOLID, new dojo.Color([0, 0, 0]), 1), new dojo.Color([25,50,225,0.3]));               var controlItems = [];          var surveyItems = [];               dojo.forEach(response,function(result){                 var feature = result.feature;                 if (result.layerName =="surveys"){         showSurveysNameGrid();             feature.setSymbol(polygonSymbol);             searchType="selSurvey";         surveyItems.push(attributes);                   var attributes = feature.attributes;         var csvTextSurv = "DOCUMENT_N;TOWNSHIP_RANGE;SECTION\n";                     csvTextSurv += attributes["DOCUMENT_N"] + ";" +                                attributes["TOWNSHIP_RANGE"] + ";" +                                attributes["SECTION"] + "\n";                  document.getElementById("csvSurv").value = csvTextSurv;              map.graphics.add(feature);        }else{          showPointNameGrid();         feature.setSymbol(markerSymbol);         var attributes = feature.attributes;         controlItems.push(feature.attributes);          searchType="selControl";          var csvTextMons = "POINT_NAME;SECTION\n";                csvTextMons += attributes["POINT_NAME"] + ";" +                                attributes["TOWNSHIP_RANGE"] + ";" +                                attributes["SECTION"] + "\n";        }               document.getElementById("csvMons").value = csvTextMons;         //add selected feature to graphics layer                   map.graphics.add(feature);      });               if(surveyItems.length >0){          showSurveysNameGrid();          var surveysStore = new dojo.data.ItemFileReadStore({data:   {identifier:'DOCUMENT_N',items:surveyItems}});             var grid = dijit.byId('grid5');             grid.setStore(surveysStore);          }      if(controlItems.length >0){         showPointNameGrid();         var controlStore = new dojo.data.ItemFileReadStore({data:{identifier:'POINT_NAME',items:controlItems}});         var grid = dijit.byId('grid4');             grid.setStore(controlStore);                }                         });      }


Bearing/Course/Heading in WGS84 ellipsoid model

There is a common formula for computing the initial bearing for a course from `lat1, lon1` to `lat2, lon2`. (And yes, I'm aware that the course will change when traveling).               atan2( sin(lon2-lon1) * cos(lat2), cos(lat1)*sin(lat2)-sin(lat1)*cos(lat2)*cos(lon2-lon1))      The formulary where I took this from mostly uses a spherical model for simplicity. My questions about this formula are:   1. What is the appropriate formula for an ellipsoid model such as WGS84?   2. If it is the same formula, does this hold generally? I.e. is it the same formula for any ellipsoid model because we deformed the coordinate system, not the coordinates? Somewhat like Vincentry's formula, which uses a reduced latitude `u = atan((1 - f) * tan(lat))`. Also note that I'm looking for a global solution, including for the poles.


Confusing offsets in shapefile written by QGIS/OGR

I think I'm stuck here. A shapefile has been updated in QGIS. Records have been deleted and those records have been marked with an asterix in the first byte of the record in the dbf file. I open the shapefile in a C program to convert it to postGIS. The deleted records are skipped. Things happen to the next record after a deleted record, not after the first deleted record though. The shx file says that the offset for the next record after the deleted record is 30948 and the accumulated offset in the shp file says 28740. Neither is right, or rather when you start reading in the shp file at the two offsets you land up in the middle of a point stream and not at the start of a record header. The offset 28740 should be right because that is what you get when you add up the offset of the deleted record and the contentlength of the deleted record. The conclusion is that the shape is corrupt. But what is driving me nuts is that you can load the shape file into QGIS and it is shown correctly. What have I missed.How can QGIS handle this.


Google Maps API: how to get the street name in multiple languages

I'm trying to get the street name in multiple languages from the Google Maps API V3. For example in Belgium there are multilingual street names: ![enter image description here](http://i.stack.imgur.com/HU0JX.png) If i set the `language=FR` i will get "Rue Brogniez" and `language=NL` results "Brogniezstraat". How can i get these names together in one `geocoder.geocode` call?


How do I obtain layout coordinates of features using arcpy?

I just read the topic about Convert point XY to page units XY using arcpy? but I didn't find a solution to my problem. I have a point layer file, I know coordinates in map view, I'd like to know coordinates in layout view. I suppose the main difference with the case of a dataframe is that the point is not a graphic element itself in layout view, so the referencing of the point is more tricky. Is there anyone who solved this (or similar) issue? ps. the dataframe where I'm working is rotated, the solution proposed in the Convert point XY to page units XY using arcpy? post doesn't work for rotated dataframe. Does anybody knows what is the base point for rotation of dataframes in ArcGIS?


Oracle spatial. Cannot load data

Have just installed QGIS 1.9 with Oracle spatial abilities. We maintain geodata in an Oracle 9.2. spatial database. I get connected to the database but I can not load the data. The table shows up with no geometry and contains "nothing". Is there anyone else who has tried this ? Any suggestions ? Svein Inge Gausemel Lillehammer


c# How can I get the spatial reference of a ArcGIS layer?

I m using the following c# code to add my layer dynamically on a arcgis silverlight application:               ArcGISDynamicMapServiceLayer mylayer = new ArcGISDynamicMapServiceLayer()     {        ID = layerName,        Url = mapServiceUrl,      };      mylayer.Initialize();      mylayer.Initialized += mylayer _Initialized;      mylayer.InitializationFailed += mylayer _Failed;      myMap.Layers.Add(mylayer );      But when I try to get the spatial reference of the layer using:               SpatialReference sp = mylayer.SpatialReference;     int testww = sp.WKID;      I get an error "Object reference not set to an instance of an object." because mylayer.SpatialReference is Null. Where m I wrong? Thanks in advance!


Display GeoJSON layer on top a SVG image instead of a map.

I'm very new to GIS, so hope this is the correct place to ask. Not sure if this will even be possible. **This is what I want to achieve:** In short, I want to display a set of polygons, lines etc. in a browser, and then display an overlay which marks specific polygons and lines etc. **This is the way I think it should be approached:**   * Take a SVG which represents a set of lines and polygons.    * Display this SVG using something like leaflet in a browser, but without displaying a world map which seems to be the basis of leaflet, open layers etc.   * Overlay a GeoJSON layer on top of the SVG.   * If possible, I want the coordinates in the GeoJSON file to be relative to a reference point in the SVG. Is this possible at all? Also, is the last point possible, that is: specifying what the coordinates in the GeoJSON file should be relative to, or are they always interpreted as GIS coordinates? Thanks for any inputs, Jaco


Configure QGIS to show only the toolbars wanted.

I friend of mine who are using QGIS in their organization for non GIS- experienced people asked me a question the other day that I didn't knew the answer to. He would like to configure QGIS to only show the tools/toolbar that those non-experienced people need to see. A common complaint is that QGIS looks too so difficault with all toolbars. QGIS is installed to this users on a server and he would like to control this so that people don't need to change the toolbars on every computer that are using QGIS. Has anyone tried to solve this kind of problem (or that understand what the hell I'm talking about ;))


How to import a multiband image into IDL?

I want to run mathematical calculation on monthly NDVI images from april to october. I stacked monthly NDVI images through april to october into an image (with 7 band) using ENVI. Now I wanna import this image into IDL. How can I do this? What kind of function or procedure should be written? Please help me?


Which tool is recommended for a PostGIS design?

I need to create a design for a spatial database. I've found this list http://wiki.postgresql.org/wiki/GUI_Database_Design_Tools, but it seems to me that this software is orientated towards a regular (not spatial) DB design.


Creation raster image from interpolated image

I have interpolated temperature data that I have. This image's max and min values are 11.9 and -1.7, respectively. And then I exported this image to raster using 'export to raster'. But the max and min values of the output raster image are different from that of the interpolated temperature image. I'm wondering why values of both images (temperature and raster) are not same. They need to be same. Is there any way to make the values of the raster image be the same as the interpolated image?


Changing Buffer table geometry dynamically in PostGIS using plain SQL

I have two tables in PostGIS, a point and a point_buffer table. The point table has a buffer_distance field, with a default value say 200. Now I want to change the buffer table geometry whenever I change the buffer_distance value in my point table. I can do this for a single row in my point_buffer table using the following:               UPDATE point_buffer     SET the_geom = (SELECT ST_Buffer(the_geom,500) FROM point WHERE gid = 1)     FROM point     WHERE point.gid = point_buffer.gid      But whenever I am trying to change the whole of point_buffer table (droping WHERE clause in my sub-query) I am getting error message: > 'ERROR: more than one row returned by a subquery used as an expression'. My question is, can I change the entire point_buffer table in one go? I know one option is using a for lopp, with the upper bound being the count value of my point table and incrementing the point.gid value . But I want to do this in **plain SQL**.


Process to load elevation data into Oracle 11G

We have the National Elevation data set from USGS and considering loading it into Oracle (11G) database.   Can someone help with the process/steps?


How to store a password during one QGIS work session?

I'd like to get into qgis-plugin development, but am rather new to both python and QGIS development. I intend to write a plugin which connects data from a mysql database to shapefiles. As the data is related to persons, I need a kind of login before accessing the data. To avoid the user from typing his password every time he invokes the script, I wonder if it is possible to retain values between calls of the script. In the python cookbook, I found global project settings - but that is not quite what I'm looking for as the values should be forgotten on closing QGIS. I hope I made my intentions clear in my clumpy English and somebody could give me a hint with which way to proceed.


Updating and exporting XY positions in UTM and WGS84 to spreadsheet using ArcGIS ModelBuilder?

On a daily basis I'm editing the geographic positions of various points of a shapefile (corner_XY.shp) by dragging them on my map. When this editing task is done I'm updating the points coordinates in their attribute table column by column by right clicking on the column and calculate geometry. I have 4 columns: Latitude - Longitude - X_UTM - Y_UTM Latitude/longitude in WGS84 (format= DD MM SS,SSS) and X_UTM - Y_UTM in WGS84 UTM Zone 15 N. I'm then exporting the attribute table with the updated fields (Latitude - Longitude - X_UTM - Y_UTM) to an excel spreadsheet to share it with my team. Exemple of output:               X_UTM        Y_UTM      Latitude            Longitude     596697.989  2276072.78  20° 33' 58.529" N   91° 9' 16.333" W     599675.251  2123364.47  20° 32' 28.838" N   91° 12' 18.624" W     575317.464  2263171.52  20° 23' 51.409" N   91° 6' 35.489" W     ....        .....       ......              ......      My project is a projected one in WGS84 UTM Zone 15 N and my shapefile containing the points that I'm moving on a daily basis (corner_XY.shp) is in WGS84 not projected. I'm using ArcGis 10.1. GOAL: I'm trying to replicate those tasks using the model builder in order to automatically update and export the new position of my points in WGS84 and UTM to an excel spreadsheet.


How to simplify vector lines into an overlaid grid in QGIS?

I'm a relatively new user with GIS software (primarily QGIS) and I'm sure this is probably quite simple but I've been trying to find an answer for ages! I have transportation network data for an area (as vector lines) and I'm wanting to simplify these to then attempt to find a spatial correlation with another set of points that I'm interested in. I was told the best way to do this was to overlay a grid over the vector lines, so that each grid cell could be classified as to whether or not it contains one of the lines or not.


Shapefile to Rectangle Overlay/Intersect Test?

I have a shapefile and a repository of images that have corner Lat/Long points. I want to search that images who overlaps/intersects to shapefile. Here is the scenario that I am doing, creating a polygon from each image corner lat/long points and placing it in temporary memory. Now how can I check whether this created polygon intersects/overlaps the shapefile? Any solution recommended for this issue?


Digitizing a city map for thesis project; Transportaion

So, basically I need a digitized map of a city that can be used for analysis. My thesis will be focusing on land-use destination, so I'll be needing the fully digitized map that has all kind of land-use categories plus the road networks, precisely the measurement of accessibility. The biggest hurdle that I've came across is, as I belong to the third world country i.e. Pakistan, the nonavailability of data for GIS. Besides, I am not an expert in the GIS. However, I am major in urban planning. Things that I would like to ask are:   * How can I extract image for the study area?   * How can I digitize the map in the fastest way?


Why do ArcGIS Area Calculations not match Official Statistics?

I am trying to calculate area of polygons using a shapefile from the world bank where countries are defined. When however I am trying to calcuate the area, the numbers that emerge are different from the official statistics on area. Morever, whenever I use a different projection (e.g. equal area) and recalculate the area, I always get the same results. Can someone tell me where the mistake could come from? Thanks!


Precision for an FGDB raster dataset?

I would like to know the precision for an ArcGIS raster stored within a file geodatabase (FGDBR format); I can't seem to find this information anywhere in their documentation (apologies if I have missed the obvious...). By precision, I mean the number of decimal places that are used to store the value for each grid-cell (not the coordinate values for spatial referencing). I know this is 6 for the GRID format, but can't find the equivalent info for the FGDBR format. The reason for asking is that I get different results for exactly the same processing when I use FGDBR vs. file stored GRID raster datasets. I am interested in the cumulative sum over the entire raster grid; there is a difference of approx. 0.5% between the two methods (i.e. a result of 5301.316 when the analysis is carried out in a GDB, compared with 5323.825 in a file storage system). EDIT: In addition, students have reported getting vastly different answers when using 64 vs 32 bit machines; I don't have the values for these though. I will add when I can repeat the results myself.


Print Google and WMS layers with MapFishPrint

My application print maps using mapfish print with geoserver. But an problem, when print layer google and wms the two images are not properly overlaying, in other words my layer wms of Brazil not not overlaying territory of Brazil in Google Maps Street. app/pdf/print.pdf?spec= {"title": "Printing Demo","layout":"A4 portrait","srs":"EPSG:900913","units":"m","dpi":254, "pages": [ {                   "bbox": [-6227217.67059, -3780794.2499401,  -691301.11508219,  586512.2832727],         "scale": 28000000,         "rotation": 0,         "comment":"",         "mapTitle":"",          "geodetic": true     }], "layers":[      {          "baseURL": "http://maps.google.com/maps/api/staticmap",          "type" : "Google",         "maxExtent": [-20037508.3392,-20037508.3392,20037508.3392,20037508.3392],         "resolutions": [156543.03390625, 78271.516953125, 39135.7584765625, 19567.87923828125, 9783.939619140625, 4891.9698095703125, 2445.9849047851562, 1222.9924523925781, 611.4962261962891, 305.74811309814453, 152.87405654907226, 76.43702827453613, 38.218514137268066, 19.109257068634033, 9.554628534317017, 4.777314267158508, 2.388657133579254, 1.194328566789627, 0.5971642833948135, 0.29858214169740677],          "extension" : "png",          "format" : "png32",          "maptype": "street",          "sensor": "false",          "tileSize": [512,512],              },     {         "baseURL":"http://mygeoserver/geoserver/datageo/wms",         "opacity":0.5,         "singleTile":false,         "type":"WMS",         "layers":["BrasilUF"],         "format":"image/png",         "styles":[],         "overview":false     }         ]}


Arcmap [ArcView License only] 10 labels at different scales

Is it possible to have different labels at different scales (without duplicating layers)? For e.g. I have loads of point shapefiles and would like to label them with ID, Name and postcode. My idea was to give more information as we zoom in - meaning 1) at full extent - no labels 2) at 1: 1,000,000 - ID only 3) at 1: 500,000 - ID and Name 4) at 1: 100,000 - ID, Name and Postcode Any ideas or thoughts?


How can I get a single part dissolve using ST_Union or ST_Aggr_Union?

I am trying to get dissolve functionality using `ST_Union` using this query :               SELECT c.fid, ST_Union(c.boundaryshape) FROM c GROUP BY c.fid;       Than I found instead of using `ST_Union` I should use `ST_Aggr_Union` :               SELECT c.fid, ST_Aggr_Union(c.boundaryshape) FROM c GROUP BY c.fid;      I am getting results like dissolve but it looks like this is multi-part. Can somebody let me know how can I get the single-part dissolve?


ArcGIS Server Manager (10.2.2) - permissions applied to folder not "sticking"

I was able to do this before updating to 10.2.2, because I have several other folders with specific permissions applied via local role store. User store is Active Directory. Data store is flat files via UNC to shared drive. Each time I edit the permissions of a newly-created folder, I choose private, un-check allow access to all users, & add the role to allowed roles. Saving says that, "permissions applied successfully," but the lock icon remains open, & checking the permissions shows that the folder is public. I have tried deleting & re-creating the folder in the server manager.


Understanding URL that IdentityManager.findServerInfo() method expects?

Does anyone have experience working with the IdentityManager object in the JS API 3.9? I'm trying to understand what URL the findServerInfo() method expects. I've tried the base REST url, and individual service URLs, but so far no luck.


OpenLayers code for not showing intended area from Google Streets

Would anyone be able to help me with the following code? I have downloaded googleclient.htm world map and want to change it to display Australia only. So I have changed the max.Extent:new OpenLayers Bounds(-20037508, -20037508, 20037508, 20037508) from displaying the whole world to the coordinates I received when I zoomed in to Australia (12249491, 15302081, 17337140, 14910723). For some reason, I'm not getting Australia, and I don't understand what I'm doing wrong.                                  <title>Webmaps</title>

This is my map:


Is there a Targer-like tool for QGis?

I'm trying to make subsurface drillhole and I heard that Target (an AutoCad tool) can help me with this task but the issue is that I work with QGis.


QGis & GeoServer WPS error

When I try to use QGis WPS with GeoServer, I get the following error in QGis :               java.lang.IllegalArgumentException: argument type mismatch      Here are the GeoServer logs :               2013-09-19 15:18:50,934 WARN [geoserver.ows] - Could not get a ServiceInfo for service wps thus could not check if the service is enabled         2013-09-19 15:18:52,809 INFO [geoserver.wps] -          Request: getCapabilities         service = wps         baseUrl = http://localhost:8080/geoserver/     2013-09-19 15:18:58,121 WARN [geoserver.ows] - Could not get a ServiceInfo for service wps thus could not check if the service is enabled     2013-09-19 15:18:58,121 INFO [geoserver.wps] -      Request: describeProcess         service = wps         version = 1.0.0         baseUrl = http://localhost:8080/geoserver/         identifier[0]:             value = gs:HelloWPS     2013-09-19 15:19:00,824 WARN [geoserver.ows] - Could not get a ServiceInfo for service wps thus could not check if the service is enabled     2013-09-19 15:19:00,918 ERROR [geoserver.ows] -      java.lang.IllegalArgumentException: argument type mismatch      I do have a wps.xml file in geoserver2.4-RC2\data_dir :               <wps>       <id>WPSInfoImpl-15378ab6:140ed6a544e:-8000</id>       <enabled>true</enabled>       <versions>         <org.geotools.util.Version>           <version>1.0.0</version>         </org.geotools.util.Version>       </versions>       <citeCompliant>false</citeCompliant>       <schemaBaseURL>http://schemas.opengis.net</schemaBaseURL>       <verbose>false</verbose>       <connectionTimeout>30.0</connectionTimeout>       <resourceExpirationTimeout>300</resourceExpirationTimeout>       <maxSynchronousProcesses>2</maxSynchronousProcesses>       <maxAsynchronousProcesses>2</maxAsynchronousProcesses>       <processGroups>         <processGroup>           <factoryClass>org.geoserver.wps.DeprecatedProcessFactory</factoryClass>           <enabled>true</enabled>         </processGroup>         <processGroup>           <factoryClass>org.geoserver.wps.jts.SpringBeanProcessFactory</factoryClass>           <enabled>true</enabled>         </processGroup>         <processGroup>           <factoryClass>org.geotools.process.geometry.GeometryProcessFactory</factoryClass>           <enabled>true</enabled>         </processGroup>         <processGroup>           <factoryClass>org.geotools.process.raster.RasterProcessFactory</factoryClass>           <enabled>true</enabled>         </processGroup>         <processGroup>           <factoryClass>org.geotools.process.vector.VectorProcessFactory</factoryClass>           <enabled>true</enabled>         </processGroup>       </processGroups>     </wps>      Anyone knows where the problem comes from ? Thanks * * * Running on Win XP


How can I import data from Excel to SAGA?

I have just started using SAGA and I want to plot some geographic coordinates (latitudes and longitudes) and I want to import them from Excel. Please tell me what to do from Excel to when the points are plotted. Thank you.


How to Extract Data by Attributes (Clip/Zip/Ship) using Attribute Where Clause - Esri Server Tools

I want to use the out-of-the-box Esri geoprocessing task to extract data via ArcGIS Server (Server Tools --> Extract Data). Its pretty straight-forward, but it doesn't seem to support use of a "where clause" as an input parameter for the service. I could add this in using python, but would love if it just supported it without custom development. Has anybody seen a way to use ArcGIS Server to extract data using an attribute where clause WITHOUT custom development. Thanks ahead of time.


trsp route between two points ( not edge to edge )

I'm trying to solve common problem, how to generate route/navigate between two GPS coordinates and not between two edges using trsp algorithm. I did read many articles here but i'm not sure how should i create temporary edge on huge psql database where even now routing taking several seconds and insert/update can take forewer. Solution where i will navigate to nearest edge and then calculate distance between edge coordinates and destination coordinates is in my case not good b/c i'm trying to avoid one way streets and this solution is ignoring it. Another solution i have is to find edge where both start(x1,x2) and end(x2,y2) coordinates are nearest to the destination GPS coordinates ( which does not need to be nearest edge ) and navigate there but found this edge is also very slow.


Which bins should be used for demographic data for a choropleth map?

I'm currently working on a project which produces a choropleth map similar to a map like this: http://dl.dropboxusercontent.com/u/54512470/Telefondaten_VIZ/index.html That map has 2x5 classifications (5 for positive migration numbers, 5 for negative ones.) I have different underlying numbers than the one I linked. How do I decide the following things:   * Linear or Exponential increase of the break points (10, 20, 30, 40, 50 vs. 10, 20, 40, 80, 160)? Or even square-root-shaped?   * related: If linear(ish), when is it ok to make some bins larger than others?   * 'final break point', i.e. what would be the equivalent of -500 and 500 in the linked example Is there generic advice on these things which can be given without looking at the data? **Edit:** The reason I don't want to post the data here is that I'd like to use some kind of algorithm to automatically create bins and resulting choropleth maps. Unlike the linked project, I don't have the migration data for a fixed timespan, but I have data for 2007, 2008, and so on up to 2011. I'd like a radio button to switch between the years (can combine several years). I don't want the colors to be consitent between different different settings (i.e. + 5000 in year 2007 may have a different color than for 2010 + 2011), but to be informative for a fixed year. One approach I was thinking of is: Take the largest absolute migration number (e.g. 1200) and round down to the nearest 50, 100, 500, 1000 and so on. This is the lower border of the dark orange, let's call this threshold number T (i.e. 1000 in this example). The threshold for the dark blue is everything below -T (i.e. -1000). The binning points of the linked project are T/10, T/5 and T/2 (in the example 100, 200, 500). I'm unsure about these particular break points. Why not T/10, T/4, T/2 (100, 250, 500)?


How to convert 'RasterBrick' object to 'im' object in R?

During script execution I need to convert 'RasterBrick' to 'im' object. Package 'geostatsp' is supposed to do this job:               library(geostatsp)     myraster = raster(matrix(1:100,10,10), xmn=0,xmx=10,ymn=0,ymx=10)     class(myraster)     >[1] "RasterLayer"     >attr(,"package")     >[1] "raster"          myIm = as.im(myraster)     class(myIm)     >[1] "im"      But when I try it with real data I catch an error:               library(geostatsp)          b <- brick(some_file)     r <- raster(b)     class(r)     >[1] "RasterLayer"     >attr(,"package")     >[1] "raster"          bim <- as.im(r)     >Error in .local(x, ...) : 'x' has no values


Get current route segment from Route type using ArcGIS Runtime LocationListener?

I am trying to get the route segment on which I am currently navigating. i.e. I want to get the current route segment when `onLocationChanged` method of my `LocationListener` is called. A route is constituted from small segments. And we have the routing directions related to a segment but I do not know how to get the current segment in my `LocationListener`.


Cleaning up overlapping polylines

I am working on creating a project boundary for a specific area. To do this I merged several data sets including roads,rivers & railroads (polylines) and lakes, civil townships & congressional townships, (polygons converted to polylines) into a single polyline feature class. I then split the polylines at every intersect in the data. I then created a subset of the data using only the nearest features to the project area to create a boundary. My problem is my data has some issues with overlap. In some areas the civil townships are coincident with bodies of water or roads and the polyline for both types of feature are overlapping (or close to overlapping) in my subset of the data. I would like to clean up the data so there is only 1 line for each section of the boundary and duplicate boundary lines are deleted. after I get a clean polyline feature I was planning and using the feature to polygon tool to create the the final boundary. Is there a simple method of cleaning up this data? I already have used topology rules to clean up a few thing such as dangles, but I couldn't really find a rule that applied to my situation. any advice would be appreciated.


Using update cursor in a while loop to perform analysis on individual polygons gets stuck on first row

I am working through an open courseware project from PSU found at https://www.e-education.psu.edu/geog485/node/145 I am trying to use an update cursor in a while loop to perform an analysis on individual polygons within a multipart feature. However, the output is for the first row ad infinitum. Here is what i have so far:               import arcpy     arcpy.env.workspace = "c:/users/documents/arcgis/arcgis tutorials/psu tutorials/project3/policedata.gdb"     arcpy.env.overwriteOutput = True          # Declare variables for files and fields used in analysis          incidents = "c:/users/documents/arcgis/arcgis tutorials/psu tutorials/project3/policedata.gdb/GraffitiIncidents"     patrol_zones = "c:/users/documents/arcgis/arcgis tutorials/psu tutorials/project3/policedata.gdb/PatrolZones"     incidents_layer = "IncidentsLayer"     patrol_field = "NAME"          cursor = arcpy.UpdateCursor(patrol_zones)     zone = cursor.next()          # Loop through each patrol zone          while zone:         try:             arcpy.MakeFeatureLayer_management(incidents,"IncidentsLayer" )             zone_name = zone.getValue(patrol_field)             print zone_name                  sql = patrol_field + " = " + "'" + zone_name + "'"              # Make a feature layer for each zone             arcpy.MakeFeatureLayer_management(patrol_zones, "CurrentZone", sql)                   # Make selection             arcpy.SelectLayerByLocation_management(incidents_layer, "CONTAINED_BY", "CurrentZone")             num_incidents = arcpy.GetCount_management(incidents_layer)                       print "Graffitti incidents = " + str(num_incidents)     ##        arcpy.Delete_management("CurrentZone")     ##        zone         except:             print "Could not perform analysis."               del cursor      What in my code is causing the script to produce results only for the first row?


How to fix QGIS crashes when using GRASS toolbox?

I've installed QGIS 2.0.1 on Ubuntu 12.04. exactly as explained on qgis.org . Unfortunately QGIS crashes when trying to use any GRASS model (r.aspect as example): just opening QGIS, click on toolbox expand toolbox, click on GRASS function r.aspect (as example) and it hangs... no error message, no possibility to define data for input! It's probably a local problem but any help is appreciated. I've had QGIS 1.7.4 installed prior. Additionally I tried on a fresh ubuntu build (12.04) in a virtualbox. same results: when trying to us a function from Grass-toolbox it crashes QGIS :-(


How to repair a data source with joins to Excel?

I have a set of MXD files, that do not use relative paths, and related data that have been copied to a new computer. Normally it's not a big problem to repair the data sources in the project, but in this case some of the layers are joined to an Excel file. When I try to repair the data sources any layer that is joined to an Excel file looses it's symbolization. Has anyone run into this problem and figured out a fix?