GeoServer preventing overwriting html files

I have GeoServer installed on Windows XP as part of the OpenGeo bundle. It's working well and I have just started creating my own app using it. I have created an html file in \WebApps\apps and can view it via localhost:8080. My problem is that once I have viewed the file in a browser it becomes locked and I can't save over it. Has anyone else had this problem? Is there a fix?


Convert surveyor's calls in excel table to polygon shapefile

Do you know of a way to take surveyor's calls from an excel file (for example: N35-45-15W 60.0 ft, etc) and have them generate polygon shapefiles point by point in GIS? I know you can easily do this with point features that have coordinates, but I am trying to figure out how to build full polygons (without XY coordinate information, just surveyor's calls).


Working out parcel direction

I am looking at applying a direction indicator to land parcels and houses (in MasterMap). effectively, I am looking at applying a direction indicator to the land and the buildings associated to the land (NW-SE, N-S, E-W etc). Is there an easy way to do this? I guess you'd need to also find the direction of the road it was situated on for completeness?


Separate tilestache cache for each mapnik

I'm using tilestache with mapnik source and redis cache. I'm wondering if there is a possibility for defining seperate cache for each layer? I want to have unique "key prefix" for each layer, therefore redis administration would be simpler (for me). What is more I want to have all tilestache definitions in one config/cfg file.


Watershed Segmentation SAGA-GIS

I'm trying to segment a topographic map to vectorize it. I'm using the watershed algorithm of SAGA-GIS. It's working fine with polygons, but with roads or contour lines is not working properly. Here is an example about I refer... the countour lines have differents seed values...   * Anyone knows how could I improve it to obtain clear lines?   * Anyone knows how to remove the square appearance after vectorization? Thank you very much! ![Contour Lines after Watershed segmentation](http://i.stack.imgur.com/hoOBm.png)


How to get ESRI Metadata into GeoNetwork?

Has anyone had success in transferring metadata created in ArcGIS 10.2 into GeoNetwork? Currently I am creating metadata with ArcGIS using the Item Description menu item in ArcCatalog. I would like to store and publish that metadata in GeoNetwork but I have not been able to get this to work. I have tried to import that xml into GeoNetwork using a variety of formats and style sheets but I either get a validation error or incomplete information.


What causes an Out of Memory error in ArcServer?

I'm creating a web app in the ESRI JS API and have recently added a 2008 FIRM service to the app. I cached the service. When I would try and click and identify certain flood zones I get an error in the console that says `ERROR: Out of Memory` I thought maybe this was due to the cache since I believe its around 7gb. So I created a new service with the data without a cache. Same Error. Is there a good way to troubleshoot this? Is there a typical reason why this error occurs? The data is only for 1 county so its not an overly large dataset by any means. I'm using ArcServer 10.1 **SERVER SPECS** 4x 2.0GHz 4th Generation Intel Xeon E5-2650 CPUs w/ 1600MHz FSB, Hyper- Threading, and Virtualization Technology enabled 8GB DDR3 1600MHz quad-channel ECC RAM VMXNET3 Ethernet Adapter connected at 10Gbs Windows Server 2008R2 – Standard Edition Virtual Hardware version 9 (vmx-09) w/ the latest VMware tools running on ESXi 5.1 Enterprise Plus **UPDATE** The memory error only happens in IE. The app works fine in Chrome and Firefox.


Frequency analysis and overlay analysis

Am new to gis and I hope to find help here. I am working on a research that requires that I overlay conditions of drainage systems and quality of healthcare in a coastal state to get areas which have potential for increase in the incidence of malaria. Each of these two attributes of areas in the state have more than one value. For instance, under quality of health care, I have data on:presence of healthcare service, drug provision and attitude of healthcare personnel. Under conditions of drainage system, data on: presence of open drainage, effectiveness of waste management. The data is derived from answers from respondents in a household survey. The data is in numeric form Am using arcgis 9.2 To my question, can I use frequency analysis in gis to summarize the data and go on to overlay the two attributes


How to disassociate roughness from slope

How to measure topographic roughness irrespective of slope in ArcGIS? I have used many methods to generate roughness layers, including TRI, VRM (which is supposed to be less correlated with slope), slope aspect index, slope variance, slope standard deviation etc. So far, all of the output layers are extremely similar to each other and to my slope layer. They do not seem to be clearly differentiating between areas of steep and rugged terrain and those which are steep but smooth. How can I generate a roughness layer which is not so dependent on slope?


Is it possible to create a standalone model in ArcMap Modelbuilder?

Is it possible to create an algorithm in ArcMap Modelbuilder that could then become a fully standalone applet? I would like to be able to provide this processing algorithm even to users who do not own an ArcMap license.


Esri basemap render speed abnormally slow in Chrome

**_Summary of Issue:_** In my project, using the oceans basemap from Esri is very slow to render when zoomed in very close (for Chrome only). Topo is also slow at high zoom levels, but it is bearable in comparison to oceans. Using Chromes developer tools to look at the Network timing for the tiles reveals some are waiting 400+ ms, but this is also the case for topo so that's not too helpful. Profiling using the Timeline tool shows that painting takes considerably longer with the oceans basemap. You can try it yourself quickly in the sandbox by replacing "topo" with "oceans": http://developers.arcgis.com/javascript/sandbox/sandbox.html?sample=map_simple **_Question:_** Is there anything on my end that I can do to optimize this in Chrome?


Open Raster stored in ESRI V10.0 Geodatabase in QGIS V2.4

Is there a way to open a raster that is stored in ESRI V10.0 Geodatabase in QGIS V2.4. I can open the vector layer, but there does not seems to be a way to open a raster.


Set a split policy/rule on a field without domain in ArcMap

Can you define a split policy or split rule on a field without a domain? When a domain is involved, you can set the default value of the field and chose a rule for when a feature is split or merged. The setSplitRule exists, but the only two options are 'None' or 'Ratio'. 'None' sets the attributes to the copy of the original value. http://resources.arcgis.com/en/help/main/10.2/index.html#//018z0000004v000000 I am looking to set a field to the Default Value (null) instead of copying the value from the original feature.


How to query one shape file using another shape file?

How to query a shapefile using another shapefile in ArcGIS? For example: I had one shapefile with states, another shapefile with towns. Now I want to query towns using the states file like towns of particular states. Note: I can extract the data using _clip_ , _union_ or _intersect_ tools in ArcGIS. Here, I don't want to extract data, I need to view the data.


C++ parser for OGC Filter Encoding 1.0

I'm developing a C++ application and I have to parse a OGC Filter in order to implement querying functionality on a metadata repository. Does anyone know any C++ library to parse OGC Filters? Thanks.


Can I simplify polygons in QGIS by Size?

Just wondering if it's possible in QGIS to simplify polygons based on size? I've done it many times using the tolerance level to reduce the overall size however in this case each polygon cannot be bigger than a certain size? Or can I simplify only those polygons which have more than a certain amount of points (which would determine size). So the very small polygons wont be affected. I'm attempting to create a google map of every polygon so I can't exceed a certain size. Thanks


Using gdalwarp with multiple elevation models

I'm trying to use `gdal_warp` to slice and convert some SRTM data. The data I'm interested in are in S30E152.hgt and S30E153.hgt. What I'd like to do is:   1. Convert these DEM files from EPSG:4326 to EPSG:3857   2. Select a subsection of the data at coordinates 17008000,-3590694 17051000,-3503623   3. Write this subsection to a GeoTIFF file for further processing Firstly, I add the two DEM files to a VRT: `gdalbuildvrt combined.vrt S30E152.hgt S30E153.hgt` Then I check the VRT to ensure the height values are sane: `gdalinfo -mm combined.vrt [snip] Band 1 Block=128x128 Type=Int16, ColorInterp=Undefined Computed Min/Max=-14.000,1501.000 [snip]` Then I use `gdal_warp` to convert and slice the VRT: `gdalwarp -s_srs EPSG:4326 -t_srs EPSG:3857 -te 17008000 -3590694 17051000 -3503623 -r bilinear combined.vrt sliced.tiff` Then I check the GeoTIFF to ensure the height values are still sane: `gdalinfo -mm sliced.tiff [snip] Band 1 Block=449x9 Type=Int16, ColorInterp=Gray Computed Min/Max=0.000,0.000 [snip]` It appears that whenever I use `gdalwarp` with a VRT containing multiple DEMs, it always clobbers the height values. For comparison, I created a VRT with only a single DEM and tried the same operations (using a smaller subsection that fits within range). `gdalbuildvrt test.vrt S30E152.hgt gdalinfo -mm test.vrt (produces Min/Max=-14.000,1501.000) gdalwarp -s_srs EPSG:4326 -t_srs EPSG:3857 -te 16931600 -3388300 16976200 -3439400 -r bilinear test.vrt test.tiff gdalinfo -mm test.tiff [snip] Band 1 Block=446x9 Type=Int16, ColorInterp=Gray Computed Min/Max=53.000,1471.000 [snip]` This seems to work OK. Maybe the problem was with my coordinates. So I thought I'd warp the entire VRT, and then select the subsection: `gdalwarp -s_srs EPSG:4326 -t_srs EPSG:3857 combined.vrt converted.tiff gdalinfo -mm converted.tiff [snip] Band 1 Block=2328x1 Type=Int16, ColorInterp=Gray Computed Min/Max=-14.000,1501.000 [snip] gdaltranslate -projwin 17008000 -3503623 17051000 -3590694 converted.tiff sliced.tiff` Although I've double-checked the coordinates, I always get this warning: `Warning: Computed -srcwin 914 1337 449 910 falls partially outside raster extent. Going on however.` Now this gives me a file with reasonable height values: `gdalinfo -mm sliced.tiff [snip] Band 1 Block=449x9 Type=Int16, ColorInterp=Gray Computed Min/Max=0.000,369.000 [snip]` However when I check to see if the data is sane, I get an empty image: `gdal_translate -scale 0 369 0 255 -outsize 100 200 -of PNG sliced.tiff sliced.png` At this point, I'm totally lost. I think I have data, but am not sure.


"Definition Query on a versioned Geodatabase"

I would like to implement Specific views to the SDE users, based on Definition Queries of common attributes , like Date created > = October 1. But on first look, the definition queries are for "feature Display setting" of each featureclass - Our aim is to implement it for the Feature dataset in a particular version The below article did not help, as I could not locate , choose Production Cartography In the Categories list. http://resources.arcgis.com/en/help/main/10.1/index.html#//0103000000s5000000 What is the procedure to display features on a common attribute, across the feature dataset in a Version.


How to capture contour lines from OSM cycle map as vector?

In QGIS, there is a plugin to insert a OSM cycle map layer, that brings as an image from http. But I can't convert this OSM Plugin layer to a raster image. I can't even imagine how to get the contour lines that are showing by the cycle map as a vector. Is this possible? If is not possible through OSM, there is another open resource to get this?


ArcMap 9.3 Hangs on trying to save

I am having trouble with a machine running Arc GIS 9.3. When one tries to save an mxd, ArcMap stops responding and just hangs. I have even tried with an blank new mxd but the moment one clicks save it hangs without even opening the dialogue box that allows you to navigate to the location to save. Please help


Projection ETRS89 (Ellipsoid: GRS80) 6-Degree Zone 33 to UTM33N - WGS84

Hi I'm having trouble to reproject a shape file in Qgis. The projection is supposed to be "ETRS89 (Ellipsoid: GRS80) 6-Degree Zone 33" and I would like to project it to **UTM33N - WGS84 (EPSG 32633)**. When I set the Layer CRS to EPSG 25833 or EPSG 4258 and save it as UTM33N-WGS84 I don't get the result as intended... This is the file I would like to reproject: http://www.projektidee.org/images/extern/landschaftsschutz.zip It should fit to this area in the end: http://www.projektidee.org/images/extern/gebiet_UTM33.dbf.zip Maybe one of you could give me a hint what I'm doing wrong... Best regards desputin


How to transfer attributes to thousands of local points from ArcSDE feature class with millions of polygons?

I have a set of points and a huge feature class with topographic information (OS MasterMap for who knows it). I'm talking millions of polygons huge. I want to transfer attributes of the polygons to the points where they intersect. I've tried spatial join, but it ran out of memory. I've tried identity and it's giving me no clue at all about how long it might take, though I'm going to let it run over the weekend to see where it ends. I then thought of making a script to take one point at a time and use that to select the underlying polygons. Doing a select layer by location with the little dialog window in ArcMap works great, it takes less than a second. However doing just that from the python window takes many minutes. I've started writing a full ArcPy script to do a select by location one by one point and then update the point's attributes with the selected polygons' attributes. Updating one point took about 15 minutes. And I've got about 8000 points to do. Unfortunately my boss won't let me sit idle and wait for a script to finish if it takes 80 days to run :( Anyone has an idea about how this might be done more quickly, and why (why ESRI, why???) SelectByLocation works so much faster in ArcMap itself? edit: maybe should mention the points feature class is local and the polygons feature class is SDE ============================================================================= SOLUTION: Among many suggestions, the only feasible one for me, being unable to install anything or use multiple workstations at the same time, was Polygeo's solution, please see the answer below. The key thing to do is: set a processing extent to a very small area around each point (do this in a cursor), make a one-polygon copy of the large dataset using this processing extent and get the attributes from the small copy. I'm bound to try to do a query within the oracle database itself as well, but I'll need some assistance from a colleague with that. As well I'm trying to get PostGIS installed and use that. I'll report back on these as well as soon as I have results. Thanks for your help!


Count points inside polygons using PostGIS

I have 2 vector data in my Postgis DB: one is a point layer (punti) with about 200000 elements and the other is a polygon (province) with 4 elements. They cover the same area and I need to find the number of points for each polygon. This is the query:               province.gid, count(*) AS totale FROM punti, province WHERE      st_intersects(province.the_geom,punti.the_geom) GROUP BY province.gid;      This is the EXPLAIN ANALYZE > HashAggregate (cost=35.35..35.36 rows=1 width=4) (actual > time=23660.787..23660.789 rows=4 loops=1) -> Nested Loop (cost=0.00..35.34 > rows=1 width=4) (actual time=0.307..23532.479 rows=199977 loops=1) Join > Filter: _st_intersects(province.the_geom, punti.the_geom) -> Seq Scan on > province (cost=0.00..1.04 rows=4 width=36) (actual time=0.003..0.009 rows=4 > loops=1) -> Index Scan using punti_the_geom_idx on punti (cost=0.00..8.31 > rows=1 width=100) (actual time=0.038..76.254 rows=69541 loops=4) Index Cond: > (province.the_geom && punti.the_geom) Total runtime: 23660.850 ms What do you think about this total runtime? Thx!


Looking to combine temporal data from multiple shapefiles to produce a days per year file for each feature

I have many shapefiles in which there is a start date and an end-date per feature. I am hoping to produce a single shapefile that has each enclosed area with a "days per year" kind of calculation. Trying to visualize this... if feature a is [ __ ] and feature b is ( __ ) and they overlap... **[ feature a 10/1/2011-10/3/2011 _( overlap ]** feature b 10/2/2011-10/8/2011 )_ into **[ feature a = 3 days ]** [ feature c= 8 days ] _[ feature b = 6 days ]_ So, overlaps calculate properly and become their own feature. Not sure if this is making sense? Any help would be great! I've tried using unions, looked at rasters, checked out the toolboxes... Not sure what route to take. My final map I'm hoping to produce should be a map that shows the duration differences. Thanks!


Extract Data widget flex viewer 2.5

I am working on the data extract widget for the flex 2.5 viewer. It has installed properly and works fine but the output is not what I expected. The tool will not divide or split polygons. For example when I draw the AOI the return will omit any polygons that have been split or divided. I was expecting to get all data in the AOI. Is this the way the tool works? Or is there something I can do to fix this? My .xml file for the widget is below: -  --> - http://s2-cor- pangaea:8399/arcgis/rest/services/demo1/GPServer/Extract%20Data%20Task draw


ArcGIS 10 Raster Mosaicing versus Lizardtech Express Server

Trying to determine if I still need to use Lizardtech Express Server to mosaic rasters together or just use the Mosaic Raster Dataset function of ArcGIS 10? I know for earlier releases of Express Server and ArcGIS you could find info on tying the 2 together but I can't seem to find anything about it now. I've learned that when the silence is deafening on a topic it generally means that "they" want it to just die and go away. Help me StackExchange GIS, you're my only hope.


Best way to build my own simple web-based map server

I've got a requirement for a project where I need some relatively simple mapping functionality on a web page. I'm familiar with the Google Maps APIs but unfortunately this is a project that will be hosted on a closed network with no internet access. Is there a relatively straight-forward way to host my own map server with the following functionality (preferably on Windows but I would consider Linux):   1. I only need to cover a map area no bigger than the size of Texas   2. I don't need streets, terrain, or any other complex layers. At most I would needs outlines of the counties in the region with a few text markers for large cities.   3. I need to be able to display custom data driven markers on the map when it is rendered (I don't care if that is on the server side or if I need to do it in javascript on the client)   4. 2 or 3 levels of zooming would be great, but no necessary I've tried SUAS server v3 with XAMPP and ran into numerous issues, including low-level PHP errors.


Finding a route with desired cost not minimum cost

The goal of routing is usually finding the route with minimum cost. That cost can be distance, time, calories or some other value but still our desire is to find the route that costs least. But what if I want to find a route with desired cost? Not the fastest route but the route that will take around an hour (let's assume I have an hour to walk from point A to B and I want to see as much on my way there as I can). Or that I want a route that will be around 2 kilometres long. Using a standard routing algorithm and rejecting the proposed routes until one with desired cost is found will probably lead to something like "Head straight from A to B and then circle around B a few times". Is there any previous work on this? Does anyone have any ideas?


How to smooth the polygons in a contour map

Here is a contour map for which all the polygons of levels are available. Let ask **how to smooth the polygons keeping all vertices preserved in their exact locations?** Indeed the contour is made on top of a grid data, you may suggest then to smooth the grid data and hence the resulting contour will be smoother. Note that this is not working as my desire since the smoothing function such as Gaussian filter will remove small packs of data and will change the range of the third variable e.g., the height which are not allowed in my application. Actually I'm looking for a piece of code (preferably in **Python** ) which can do the smoothing of 2D polygons (any type: convex, concave, self-intersecting etc) reasonably painless (forget pages of codes) and accurate. FYI, there is a function in ArcGIS that does this perfectly, but using third party commercial applications is not my choice for this question. ![enter image description here](http://i.stack.imgur.com/nAeyI.png) **Updates**   **1)** Scipy.interpolate: ![enter image description here](http://i.stack.imgur.com/ZLMhT.jpg) As you see the resulting splines (red) are not satisfactory! **2)**   Here is the result using the code given in here. It is not working well!   ![enter image description here](http://i.stack.imgur.com/JL9aP.png) **3)**   To me the best solution should be something like the following figure in which a square is being smoothed gradually by changing only one value. I hope for similar concept for smoothing any form of polygons. ![enter image description here](http://i.stack.imgur.com/xHyxR.png) Satisfying the condition that spline passes the points:   ![enter image description here](http://i.stack.imgur.com/5cmVu.png) **4)**   Here is my implementation of "whuber's idea" line by line in **Python** on his data. There are possibly some bugs since the results are not good. ![enter image description here](http://i.stack.imgur.com/b0E2b.png)   K=2 is a disaster and so for k>=4. **5)**   I removed one point in the problematic location and the resulting spline is now identical to whuber's. But it is still a question that why the method does not work for all cases?   ![enter image description here](http://i.stack.imgur.com/3isPI.png) **6)**   A good smoothing for whuber's data can be as follows (drawn by vector graphics software) in which an extra point has been smoothly added (compare with update 4):   ![enter image description here](http://i.stack.imgur.com/4cVxi.png) **7)**   See the result from Python version of whuber's code for some iconic shapes:   ![enter image description here](http://i.stack.imgur.com/803BS.png)   **Note** that the method seems doesn't work for polylines. For the corner polyline (contour) green is what I want but got red. This needs to be addressed since contour maps are always polylines although closed polylines can be treated as polygons as in my examples. Also not that the problem emerged in update 4 has not been addressed yet. **8) [my last]**   Here is the final solution (not perfect!):   ![enter image description here](http://i.stack.imgur.com/YO8TF.png)   Remember that you will have to do something about the area pointed by stars.There is perhaps a bug in my code or the proposed method needs further development to consider all situations and to provide desired outputs. Thanks and Good luck!


Looking for market share of GIS desktop tools

What are the relative market shares of the various desktop GISes? I'm interested in this within the open-source and non-opensource realm. One article that I read, is that the market share of ESRI is estimated to be approximately 30% of the global GIS market in 2009 (Daratech, 2011 - dead link unfortunately) but that they have a much higher perceived 'mind-share' (Batty). I started thinking about this as I wondered which was the most popular open- source GIS program, and assumed it was QGIS but I'm finding it hard to locate any numbers. To clarify (in response to a comment), I am interested in GIS Desktop tools that allow quantitative analysis (which would exclude Google Earth).


OpenStreetMap or Google Map Maker: where should I donate and volunteer my geographical knowledge?

I definitely know a lot about the places I visit and want to share it with everyone. I want to know if donating and volunteering at OpenStreetMap—which is where I _have_ been helping—is a waste of time: it is possible that it is all for naught to help on such a relatively low-trafficked site. On the other hand, I am concerned about the closed, commercial nature of Google Maps, even though the traffic is high and my knowledge would likely get more reuse. Where should I help with collaborative mapping?


NDVI from QGIS is giving me NaN values

NDVI,TTVI maps are nan values. when I do the NDVI calculation using the raster calculator it looks fine. what's going on?


qgis-web-client stuck on Loading Map

I'm trying to publish a map created in qgis desktop via qgis-web-client. The left hand bar loads with the layers I would expect visible but I never actually see any mapping data just a grey screen with "Loading Map....". If I make the call to the qgis server manually via a web-browser e.g.               http://adddress/qgis/qgis_mapserv.fcgi?map=/var/www/qgis-web-client/projects/test/filename.qgs&&_dc=1394809120649&SERVICE=WMS&VERSION=1.3&REQUEST=GetCapabilities      I get what looks like valid XML returned. In the chrome dev/javascript console I'm seeing error: ` Uncaught TypeError: Cannot read property 'abstract' of undefined ` From webgisinit.js line 1680 ` var thisAbstract = wmsLoader.projectSettings.service.abstract; ` Anyone seen this before or got any clues how to fix?


Working with GDB in QGIS 2.0 - plain language instructions please!

I'm a QGIS user running v2.0 (Dufour) on Mac OSx10.8.5. I have some data in .gdb format and need to access this and work with it in QGIS. I'm a relatively savvy computer user, but not at all a programmer or developer. I've been reading the forums and have seen some steps outlined for GDAL installs etc., and compiling etc., but am a bit daunted and don't want to mess something up on my system. Can some patient soul please spell out for me the steps I need to take to get .gdb files functioning in QGIS? Alternately I could try to use GDB Flee or something like it, but not sure where to start. Thank you so much. I'm learning so much from reading these forums but still have a long ways to go!


One-To-Many Table Relate Viewed Through ArcGIS Server

I want to create a Table Relate with a one-to-many relationship between a layer and a non-spatial table based on a common field. My question is, if I publish this data to ArcGIS Server, will an identify operation on that layer return _all_ of the related values in the non-spatial table? I want to try this, but I don't have sample data that meets these requirements (if anyone knows where I can find some, it would be very helpful). Here is an example from ESRI of an Identify operation in ArcMap on two tables with a one-to-many relationship. I need to know if the results in ArcGIS Server would be consistant. I am using ArcGIS 10 and SQL Server 2008 R2. Thanks, Seth


Monitoring programmes

I am not sure if it is a right place to ask this one but let me try. Monitoring programmes are known for various targets and for various spatial levels. They may include regions, countries or even continents. It is widely accepted that monitoring systems are of increasing importance as economic activities and their impacts on nature become more and more global. I am looking for a various monitoring programmes. The best will be that related to forestry (but in fact i am seeking for anyone) like this one I am interested how do you evaluate the selected level I and level II indicators used in the ICP forest ecosystem monitoring programme? Do you consider them exhaustive to get a comprehensive view of the state of the forest ecosystems under the influence of air pollution? Also if you have some examples in mind it will be great to share it!


Find AB_CSRS.gsb for NAD_1983_TO_WGS_1984_8 Transformation

I am trying to change my WGS84 data to NAD83 for Alberta data. My research tells me that I need to use the NAD_1983_TO_WGS_1984_8 Datum Transformation which is in the AB_CSRS.gsb. I can not find this gsb file for download on the internet. I did manage to find and obtain the NTv2_0.GSB but that seems to only be for NAD27 To NAD83. Does anyone know where I can get the AB_CSRS.GSB file? Thanks


Split Lines at Intersection of Other Lines

I have two line layers (roads and contours in this example) and I would like to split the roads layer every time it crosses/intersects a feature on the contour layer. I can't seem to find a tool that does this, as the intersect tool is expecting a polygon for at least one of the layers. It also appears that GRASS' v.clean only works on a single layer. It seems like a simple task, but I'm at a loss!


How to make a different font and color of text in labels?

Is it possible to make a bold or italic, and change the color of individual words Label? Example: the first line is bold in red text, font 12, second row 10 green italic font ![enter image description here](http://i.stack.imgur.com/0y3Gr.png)


Misalignment when generating centroids with Qgis

I was wondering if anyone could shed a light on why this behavior is happening in QGIS ? ![alt text](http://i.stack.imgur.com/mziOb.jpg) I'm pretty confident that the issue here is not due to "centroids within polygon" **Edit** Added a detailed image ![alt text](http://i.stack.imgur.com/JndXW.jpg)


Issues using Arc 10.0 extension recompiled for 10.1

I have a written an ArcGIS extension in C# .Net for version 10.0. I have attempted to recompile this extension for 10.1 and that works. But running the extension in 10.1 is causing ArcGIS to crash. While debugging the code I'm getting an error saying "ExecutionEngineException was unhandled". This is not always happening on the same part of the code. One time its a cursor, another time its setting symbology, etc. its different every time. I have read through the documentation here on how to upgrade and none of these solutions are helping. Has anyone ever experienced this before? I have no idea what else to try.


Network Analyst hierarchy routing issues. (ArcMap)

My friend and I have created a walking network. We created a network without a hierarchy evaluation and the routing worked fine, but when we added a hierarchy to the equation many routes won't solve. We created a hierarchy class and attributed each of our walk types with a hierarchy of 1, 2, or 3. We thought that the cause of the routes not working was when a route would jump from a hierarchy of 3 to 1, or 1 to 3, skipping 2. But routes do work in some of those situations. Does anyone know why it breaks, or if there is a way to fix it? (We know that all of our paths are snapped together)


Query to select MIN and MAX values by groups

In ArcMap 10 i'm working in a Personal Geodatabase where the points (see screenshot below) represent the number of each house in the street. Each point has an ObjectID. This ObjectID is an unique value for each "link" of the street to which all the points are related.![enter image description here](http://i.stack.imgur.com/swdTz.jpg) I have a total of 35,000 points (house numbers), which are all related to 14,000 different links. What I want is to select both the MAX and MIN value in each side of the street. Each side of the street also has a unique ID, which appears in my attribute table as CVEFT. On every link, the CVEFT (for both sides) is different. This is what my Table looks like: ![enter image description here](http://i.stack.imgur.com/pBK0H.jpg) As you can see in the screenshot above, there are several points (NUMEXT) related to the same link (OBJECTID). You can also see that all the points related to a single link may have 2 different street-side codes (CVEFT). I want to build a SQL query that allows me to select the MIN and MAX values for every Link, for each side of the street. Does anyone have any idea of how to do this? Any ideas would be appreciated! Thanks!


Longitudinal and Latitudinal coordinates of interstate exits

I need to find the coordinates of all the interstate exits for some research. I'm not looking to code up a program that will do it for me, but I was hoping that someone had already done it and would be willing to share that information with me. Any help would be much appreciated.


Should Python Add-In installed in map which is packaged be part of *.mpk?

Using ArcGIS for Desktop 10.1 SP1, I just created a map package (.mpk) from a map (.mxd) which had a Python Add-In (.esriaddin) installed. When I unpacked the .mpk the .esriaddin was neither installed nor present in the folder that was created in C:\Users\Graeme\Documents\ArcGIS\Packages. Intuitively, I would expect it to be at least present in the package, so am wondering if it that capability has been omitted by accident or design?


R, GeoJSON and Leaflet

I recently learned about leafletjs.com from an R-Bloggers.com post. One such tutorial that I would like to achieve is to create interactive choropleth maps with leaflet. I have been using the rjson package for R to create the data.js file to be read by leaflet. Although I have had success with using the provided shape file as a readable JSON file in leaflet, I am unable to repeat the process when trying to merge additional properties from the data frame ("data.csv") to the JSON file. In this case, I have done rGIS to attach data on the number of cans in each school listed in the data frame. What I would please like to achieve is to create a choropleth map in leaflet that displayes high school district (as identified by the NAME variable), and the sum of "cans". The issue, I believe, is that writeOGR exports the information as points, rather than polygon?               {         "type": "Feature",         "properties": {             "name": "Alabama",             "density": 94.65         },         "geometry": ...         ...     }              ###load R scripts from dropbox         dropbox.eval <- function(x, noeval=F) {         require(RCurl)         intext <- getURL(paste0("https://dl.dropboxusercontent.com/",x), ssl.verifypeer = FALSE)         intext <- gsub("\r","", intext)         if (!noeval) eval(parse(text = intext), envir= .GlobalEnv)         return(intext)         }              ##pull scripts from dropbox          dropbox.eval("s/wgb3vtd9qfc9br9/pkg.load.r")             dropbox.eval("s/tf4ni48hf6oh2ou/dropbox.r")              ##load packages         pkg.load(c(ggplot2,plyr,gdata,sp,maptools,rgdal,reshape2,rjson))              ###setup data frames         dl_from_dropbox("data.csv","dx3qrcexmi9kagx")         data<-read.csv(file='data.csv',header=TRUE)              ###prepare GIS shape and data for plotting         dropbox.eval("s/y2jsx3dditjucxu/dlshape.r")              temp <- tempfile()         dlshape(shploc="http://files.hawaii.gov/dbedt/op/gis/data/highdist_n83.shp.zip", temp)         shape<- readOGR(".","highdist_n83") #HDOE high school districts           shape@proj4string               shape2<- spTransform(shape, CRS("+proj=longlat +datum=NAD83"))                   data.2<-ddply(data, .(year, schoolcode, longitude, latitude,NAME,HDist,SDist), summarise,                 total = sum(total),                 cans= sum(cans))                   ###merging back shape properties and data frame         coordinates(data.2) <-~longitude + latitude         shape2@data$id <- rownames(shape2@data)              sh.df <- as.data.frame(shape2)              sh.fort <- fortify(shape2 , region = "id" )              sh.line<- join(sh.fort, sh.df , by = "id" )                        mapdf <- merge( sh.line , data.2 , by.x= "NAME", by.y="NAME" , all=TRUE)              mapdf <- mapdf[ order( mapdf$order ) , ]                   ###exporting merged data frame as JSON     mapdf.sp <- mapdf     coordinates(mapdf.sp) <- c("long", "lat")     writeOGR(mapdf.sp, "hssra.geojson","mapdf", driver = "GeoJSON")      However, it appears that my features are repeating itself constantly. How can I aggregate the features information so that it looks more like the following:                 var statesData = {"type":"FeatureCollection","features":[     {"type":"Feature","id":"01","properties":{"name":"Alabama","density":94.65},"geometry":    {"type":"Polygon","coordinates":[[[-87.359296,35.00118],[-85.606675,34.984749],  [-85.431413,34.124869],[-85.184951,32.859696],[-85.069935,32.580372],[-84.960397,32.421541],[-85.004212,32.322956],[-84.889196,32.262709],[-85.058981,32.13674],[-85.053504,32.01077],[-85.141136,31.840985],[-85.042551,31.539753],[-85.113751,31.27686],[-85.004212,31.003013],[-85.497137,30.997536],[-87.600282,30.997536],[-87.633143,30.86609],[-87.408589,30.674397],[-87.446927,30.510088],[-87.37025,30.427934],[-87.518128,30.280057],[-87.655051,30.247195],[-87.90699,30.411504],[-87.934375,30.657966],[-88.011052,30.685351],[-88.10416,30.499135],[-88.137022,30.318396],[-88.394438,30.367688],[-88.471115,31.895754],[-88.241084,33.796253],[-88.098683,34.891641],[-88.202745,34.995703],[-87.359296,35.00118]]]}},          {"type":"Feature","id":"02","properties":{"name":"Alaska","density":1.264},"geometry":   {"type":"MultiPolygon","coordinates":[[[[-131.602021,55.117982],[-131.569159,55.28229],[-131.355558,55.183705],[-131.38842,55.01392],[-131.645836,55.035827],[-131.602021,55.117982]]],[[[-131.832052,55.42469],[-131.645836,55.304197],[-131.749898,55.128935],[-131.832052,55.189182],[-131.832052,55.42469]]],[[[-132.976733,56.437924],[-132.735747,56.459832],[-132.631685,56.421493],[-132.664547,56.273616],[-132.878148,56.240754],[-133.069841,56.333862],[-132.976733,56.437924]]],[[[-133.595627,56.350293],


How to select multiple polygons or points from WFS or WMS and getFeatureInfo

Is there a way to select multiple polygons and/or points from a WFS or WMS service and then get the feature information for each selected item? There are numerous examples on the OpenLayers website but none of them show how to multi-select AND get the attributes for the selected items. For instance I can get the feature info for a single selected item by using the getFeatureInfo option for a WMS service but this does not allow me to select multiple items. I absolutely have to have the ability to select multiple items! Also, I can select multiple items by drawing a box, a freehand polygon, or multi-selecting individual features using a number of different methods but none of those methods have the ability to get the feature info. What good is being able to select multiple items but not be able to get the attributes for the features? Am I missing something really obvious? My goal is to be able to select multiple features and then extract my unique identifier (point number or polygon number) from each of the features info in order to query data from a SQLServer database. I don't need to do anything fancy to the spatial data all I want to do is select multiple items and take the unique id from each of them so I can query my database for info that relates to the selected items. I am using OPENLAYERS to build my map interface and GEOSERVER to serve my WMS/WFS services. I can connect to my point and polygon shapefiles using either WMS or WFS but I can not find a way to select multiple items and then get the attribute info for the items. Is there a way in OpenLayers to select multiple items and then loop through the list of items? I have been stuck on this problem for four very long and frustrating days and if anyone has any ideas on how to do this I would be eternally grateful!! 10/11/2012 - HERE IS THE CODE I ENDED UP USING //THIS FUNCTION UPDATES THE TEXT BOX WITH THE LIST OF SELECTED ATTRIBUES WHEN THE USER SELECTS A NEW POLYGON OR MULTIPLE POLYGONS var update = function(feature) { names.push(feature.attributes['yourLayersAttributeNameHere']); showStatus(names.join(", ")); };               //THIS FUNCTION UPDATES THE TEXT BOX WITH THE LIST OF SELECTED ATTRIBUTES WHEN THE USER UNSELECTS A POLYGON     var clear = function(feature) {       name=feature.attributes['yourLayersAttributeNameHere'];          for(var y in names){            if (name === names[y]){              names.splice(y,1);             }         }         showStatus(names.join(", "));     };                         //INSTANTIATE A NEW CONTROL THAT ALLOWS THE USER TO SELECT FROM THE AVAILABLE WFS FEATURE      // var control = new OpenLayers.Control.SelectFeature( nameOfyourLayerHere, /{ clickout: true, toggle: true, multiple: true, hover: false, //if hover is set to true it takes away the selection functionality.   toggleKey: "ctrlKey", // ctrl key removes from selection /multipleKey: "shiftKey", // shift key adds to selection /onSelect: update, //ON SELECT CALL THE UPDATE FUNCTION onUnselect: clear, //ON UNSELECT CALL THE CLEAR FUNCTION // // } // );               //ADD THE CONTROL TO THE MAP CONTROLS COLLECTION      // map.addControl(control); //ACTIVATE THE NEW CONTROL //control.activate();


Select only features within the layout view

I want to create a tool layer in ArcGIS10 that will select only features within the current layout view so that I can add the attribute table to the layout and only show attributes for the selected features. Has anyone got any ideas on how to do this? I want the selection to be done automatically each time the mxd is opened, and not have to do it manually.


What is the Difference between Geometric and Geographic columns?

I am a beginner in GIS and Learning PostGIS. I came across this two terms Geom and Geog, Could anyone explain in brief. What is the best to use with map data and Do all functions support both Geom and geog.


Qgis; How can I add only mean values instead of coun, sum and mean from QgsZonalStatistics (zonal statistic)?

I am using Qgis 2.2.0 Valmiers. I have ASCII grid, total 864 and I have a polygon layer as a zone. after zonal static batch run, attribute fields stop at 2046, fields number shown from "Table Manage" plug-in. Question 1. Is there any way I can create only mean column field instead of default three fields? Question2: Does QGis have a limit on number of column fields which I don't think? Here is a my code and it runs in a Python console. from qgis import core from qgis.utils import iface from qgis.analysis import QgsZonalStatistics import os,sys,time print "Process started: ", time.ctime() s_tm = time.time() f_list = [] LOCTN = "C:\1SFO\Tye\Data\2011\ASC_QGIS\" polygonLayer = QgsVectorLayer(u'C:/1SFO/SHP/WS_SFO_epsg4326_2011.shp','WS EPSG4326_2011','ogr') QgsMapLayerRegistry.instance().addMapLayer(polygonLayer) iter = polygonLayer.getFeatures() for files in os.listdir(LOCTN): extension = os.path.splitext(files)[1][1:] if extension == "asc": f_list.append(files) for y in range(len(f_list)):                 FLNAME = f_list[y]            INFILE = LOCTN+f_list[y]            RSTRNAME = FLNAME.replace(".asc","")[29:36]            rasterFile=LOCTN+ f_list[y]            zoneStat = QgsZonalStatistics (polygonLayer, rasterFile, RSTRNAME, 1)            zoneStat.calculateStatistics(None)      print "Process Ended: ", time.ctime() Thanks a lot


error using the QGIS "terrain profile" plug in

I have searched the home page for the "terrain profile" tool but have not found why I am getting error when trying to utilize it. Any ideas? I have a polyline and 4 dems ready to run the tool.![enter image description here](http://i.stack.imgur.com/u40u2.png) Would like to export the profile as an image and excel file.


I am using the WGS 84/UTM coordinate sysytem. When I go to > Project Properties / General / Canvas Units to Degrees, Minutes, Seconds my coordinates are wrong. I am using QGIS 2.0.1 ![enter image description here](http://i.stack.imgur.com/fQetf.png)


Feature class to Feature class not working

I am trying to use the Featureclass to Featureclass conversion tool to make a copy of an existing layer (in a geodatabase). My primary motive is to copy over only a few (2 out of 18) fields to the new feature class. I have put a simple script to do this, but this copies all the fields even tough I specify the field names that should be copied. Can someone please let me know what am i doing wrong. I want the output feature class (Test) to have two fields A and B.               # Import arcpy module     import arcpy               # Local variables:     Input = "C:\\Work\\Projects\\0a_IL_Naperville\\Segmenter_v8\\1 00001 IL Naperville 2013.gdb\\Streets"     Out_location = "C:\\Work\\Projects\\0a_IL_Naperville\\Segmenter_v8\\1 00001 IL Naperville 2013.gdb"          keep = ['ST_NAME','FUNC_CLASS','Shape_Length','OBJECTID','Shape','A']          arcpy.FeatureClassToFeatureClass_conversion(Input,Out_location, "Test",field_mapping=keep)


How to reclassify raster data by their sign?

I want to do the following: I have two data raster sets with different units and values. One results from measured metric data, the other from modelled grid values. What I want to compare is simply whether the data sets differ in where they show positive and negative values. So - all negative data on both sets should be "-1", all data between -0.001 and 0.001 should be "0", and all positive values simply "1". Do you have an idea how to do this best on ArcGIS 10? I suspect raster calculator...? Any hints are appreciated.


Is it possible to grow QGIS symbol size over time?

I've got a data set of people movements within my municipality. (between points of places) I've managed it with the plugin for QGIS 'time manager'. Over time people come and go, and based on a count of movements of a specific point at a specific time (time manager), i want to let my point symbol grow. Is this possible with QGIS rule based symbol?


Changes made using wfs-t and GeoExt don't get saved to PostGIS

I'm using GeoExt to display and edit a vector layer from GeoServer (running with Apache Tomcat) but if I modify the shape of an existing polygon (or create a new one but I'm trying not to leap ahead of myself) the change doesn't get written back to PostGIS through GeoServer. I'm using code obtained from a tutorial at http://workshops.opengeo.org/geoext/wfs/wfst.html. I see a few people have had problems with the same code but can't find any solutions that solve my problem. The code for the save button action:                   handler: function() {       app.featureGrid.store.protocol.commit(        vectorLayer.features, {         callback: function() {             var layers = app.mapPanel.map.layers;             for (var i=layers.length-1; i>=0; --i) {               layers[i].redraw(true);             }             app.featureGrid.store.reload();         }       });     }      My map layers get redrawn and the polygon changes persist but if I reload the page, the original shape is back. There is no error in Chrome developer tools but if i look in the GeoServer log file, I get:               2012-05-10 14:51:31,311 INFO [geoserver.wfs] -      Request: getServiceInfo     2012-05-10 14:51:31,321 ERROR [geoserver.ows] -      org.geoserver.wfs.WFSTransactionException: Feature type 'dev_apps' is not available:          at org.geoserver.wfs.Transaction.execute(Transaction.java:222)         at org.geoserver.wfs.Transaction.transaction(Transaction.java:108)         at org.geoserver.wfs.DefaultWebFeatureService.transaction        (DefaultWebFeatureService.java:162)      I'd really appreciate some help as I don't know where to look anymore!


Import and Plot MS Excel file in ArcGIS (v9.3) ModelBuilder

I am trying to create a tool using ArcGIS ModelBuilder (v9.3) that will allow me to import an Excel file, and output it as a shapefile in a coordinate system of my choice - or, if choosing the CRS is not an option, in WGS84. I have set the input dataset (as part of the model), but can't seem to find a tool to facilitate the import. I have an ArcView licence.


How to migrate SVG-based styles for QGIS to other computers?

I'm working on a project that has involved creating many recolourable SVG icons and bringing them into QGIS. These are then used in layer symbologies, with these subsequent layer styles exported as .qmls and .slds. When I open these .qmls on another computer the svg points are not attached - rather there is just a blank spot. I've tried to get just the svgs to load by exporting the style library from the computer that all of the svgs were created on. However, when I imported this into another computers style library everything was once again blank. Putting the folder with all of the svgs in the C:\Program Files (x86)\Quantum GIS Lisboa\apps\qgis\svg directory is not working either...they show up in the option to select svg icon, however when a layers qml file is loaded it is once again blank. I'm trying to get my head around this, so if there's an obvious answer please let me know. How can I move these svg/style to other computers so that when I load a file with a set qml layer the proper svgs load as well? Thanks anyone for your help.


Querying intersecting areas areas - Oracle Locator or Spatial?

My spatial database is Oracle. We are wondering if we can use basic Oracle Locator or need to get the more expensive, full, Oracle Spatial. **My application's requirements include having to be able to query "which geometries in the spatial db intersect this geometry X".** The whole locator / spatial issue is a bit confusing since it looks like Oracle initially installs all Oracle Spatial features even if you are just using basic Oracle db and do not have license for full Oracle Spatial. And you may end up using the Spatial features by mistake: http://www.spatialdbadvisor.com/oracle_spatial_tips_tricks/55/oracle-locator- vs-enterprise-licensing-the-sdogeom-package What I understand of Oracle documentation http://docs.oracle.com/cd/E11882_01/appdev.112/e11830/sdo_locator.htm#SPATL340 is that geometry function SDO_GEOM.SDO_INTERSECTION is **not** supported by Oracle Locator. So does this mean that we really need full Oracle Spatial license to fulfill the requirements? Or is there some other way then "SDO_GEOM.SDO_INTERSECTION" of finding about whether areas intersect? We are using GeoServer on top of Oracle (and accessing the geometries with WFS-T). Does this change anything, will the intersection be somehow calculated in GeoServer (doubt it?) or will the intesection-calculation have to be done in Oracle (probably?)?


identify the fraction of raster cells and listing all the cell values in every polygon

In Arcgis 10 I have a netcdf file that maps the fractions of cropland in the world from 1700 to 2007. I extracted the years that i needed as raster layer format (continous data from 0 to 1) and now, using a shapefile of the world's countries, I would like to identify every raster's cell inside every country (possibly in a table format where I find id of the cell, country name, fraction cultivated). Is it even possible to do it? Moreover I would like to have also the fraction of cell belonging to a country and the fraction belonging to another one (i.e. in the boundary). Please I need your help, I am new in Arcgis and Gis in general!


How to view values from stacked raster maps of a particular pixel at the same time (ArcGIS)

Re: ArcGIS v10.1 Hi. I have a stack of raster maps. I would like to click on a pixel, and view ALL of the values for that particular pixel at the same time. The only tool that I could find is the pixel inspector, which only displays the value for one map at a time. I want to compare all values from my layered maps for this particular location at the same time. Has anyone done this with ESRI software? Cheers


How to add graphs - pie charts, bar diagrams in geoexplorer?

Does any one know how to add interactive graphs in Geoexplorer as shown in this video ? Do we have to incorporate Google fusion tables ? Any kind of help will be appreciated. https://www.youtube.com/watch?v=d9ADUHGUXL0


How to make elements size fixed (Arcobjects)?

I'm using arcobjects .Net 10 to add graphic elements to my map. My problem is that when the map scale changes, the element size change too. How to fix element size at all scales?


Save GPX from drawn feature in OpenLayers

I am trying to modify the Vector Formats example to only display the GPX in the textarea rather than all the options in the original example but I can't get it working. My code below currently displays [object Object] is the text area. The next stage would be to save it as a GPX file but I can't even get this part to work. Help appreciated.

Draw route to update the output below.

<textarea id="output"></textarea>


How to set custom fullextent?

I have a WPF Map control whose extent I would like to restrict to a certain area (e.g. to the full extent of some layer). That is, a user should not be able to zoom out too far, or shift the map to an area outside that full extent. How can I restrict Map display to a certain extent? This is my code:               <esri:Map x:Name="Map" SizeChanged="Map_SizeChanged">         <esri:ArcGISTiledMapServiceLayer x:ID="LB" ... />         ...     </esri:Map>      And here is the code-behind (this code does not work):               void Map_SizeChanged(object sender, SizeChangedEventArgs e)     {         if (currentExtent > Map.Zoom(Map.Layers["LB"].FullExtent)) // check the current extent of the map. If currentExtent is more than FullExtent         {             // jump back to full extent             Map.Zoom(Map.Layers["LB"].FullExtent);         }     }      I'd like to get something like that:               private void Map_SizeChanged(object sender, SizeChangedEventArgs e)     {         … // get current extent of the map              … // check the current extent of the map.           // if currentExtent is more than FullExtent jump back to full extent     }      I'm using the ESRI API for WPF and C# in Visual Studio 2010.


openlayers cluster strategy - exhibiting custom icons for clusters of 1 feature

i have a page that displays a **vector layer** imported from a **kml file** (basically like the the sundials example). it works perfectly well, both using fixed strategy ONLY, as well as using both fixed and CLUSTER strategies simultaneously. i would like to create a **hybrid display** , whereas **"clusters of 1 feature"** would be displayed with the original **custom icon** contained in the KML file under `<Style><IconStyle><Icon><href>img/arrowRed.png`, for example. right now, **if i use cluster strategy, clusters of 1 (feature) are displayed using the default icon** (yellow disk). i would rather not use any plugins or libraries that aren't standard openlayers. **any suggestions?** below a part of the original javascript code using clusters (remove the cluster strategy declaration `new OpenLayers.Strategy.Cluster()` and custom icons are displayed just fine):               var urlKMLClient = 'KMLClientsAll.kml';      var layerKMLClient = new OpenLayers.Layer.Vector("Clients", {             strategies: [new OpenLayers.Strategy.Fixed(), new OpenLayers.Strategy.Cluster(), refresh],             protocol: new OpenLayers.Protocol.HTTP({                 url: urlKMLClient,                 format: new OpenLayers.Format.KML({                     extractStyles: true,                      extractAttributes: true,                     maxDepth: 2                 })             })         });       ![WITHOUT CLUSTER](http://i.stack.imgur.com/3muMR.jpg) see Client custom icon (green marker) gives way to cluster default icon (yellow disk) below ![WITH CLUSTER](http://i.stack.imgur.com/OTnYY.jpg)


Geoserver + MongoDB: connection OK but layers list empty

I setted up Geoserver (ver. 2.5.2) with the MongoDB plugin (ver. 11.2 - http://repo.opengeo.org/org/geotools/gt-mongodb/11.2/). The connection between Geoserver and MongoDB is OK and when I create a new datastore in geoserver->stores->Add new Store, it creates the database on mongo as expected. The problem is that, if I go in layers->Add a new resource and select the correct resource, the layer's list remain empty. On the top i see only "You can create a new feature type by manually configuring the attribute names and types. Create new feature type..." but i cannot find documentation on how to use that. Any help? Thanks in advance **EDIT:** in the goserver's log i've just found this error:               ERROR [data.mongodb] - getLayers error; command failed [mapreduce]:          { "serverUsed" : "[ip:port]" ,          "errmsg" : "exception: ReferenceError: mapfields_recursive is not defined near 'ction()         { mapfields_recursive (\"\", this);' " , "code" : 16722 , "ok" : 0.0}       but I can't still figure out what is the problem and how to solve it. Any ideas?


How can I extract variable information from an OpenLayers 3 object?

I have tried to make a simple map with OL 3's Geolocation module. I succeeded with a resource wasting method. There must be a simpler way to extract variable information from an OpenLayers 3 object. My code is the following:               var view = new ol.View({         center: [0, 0],         zoom: 16         })          var map = new ol.Map({         target: 'map',         layers: [             new ol.layer.Tile({                 source: new ol.source.OSM()             })         ],         view: view         });          var geolocation = new ol.Geolocation({         projection: view.getProjection()     });          geolocation.setTracking(true);          var geolocation2 = new ol.Geolocation({         projection: view.getProjection()     });          geolocation2.setTracking(true);          var featureSource = new ol.source.Vector({     });          var feature = new ol.Feature();          feature.bindTo('geometry', geolocation, 'position')         .transform(function() {}, function(coordinates) {           return coordinates ? new ol.geom.Point(coordinates) : null;         });               featureSource.addFeature(feature);          var featureLayer = new ol.layer.Vector({         source: featureSource,         name: 'Your position'     });          view.bindTo('center', geolocation2, 'position');          map.addLayer(featureLayer);      The code works, but not the way it should operate. It uses two `Geolocation` instances and because the `bindTo` method the view periodically resets. I couldn't manage to extract the `position` variable of the `Geolocation` object. I've tried `Geolocation.getPosition()`, `Geolocation.get('position')` and `Geolocation.p.position`. Every one of them acted like a void method and returned an `undefined` value. I use Google Chrome version 35.0.1916.153 on a Debian Wheezy. To complete the story, if I use only one `Geolocation` instance and bind it to the `feature` and the `view` objects, the point follows the center of the map with every change in the view.


Non-spatial Data Driven Pages

Is it possible to utilise the Data Driven Pages (DDP) tool in ArcMap or ArcPy to create a set of maps with the same extent but showing different layers? To clarify I regularly make large sets of maps for projects all in the same location often with the same baselayers but modified to focus on different aspects for the relevant section of a report such as monitoring locations, land-use, key features etc. which all require separate .mxds to be created and maintained. A DDP type system would allow me to utilise the dynamic text and ensure a standard template / output throughout, it would also make life a lot easier when updating a general figure number or comment as a report is written in one .mxd rather than multiple for a client. I'm using ArcMap 10.1, but having read through the ArcGIS DDP online help I assume this would have to be done with ArcPy as I've been unable to figure out how to do it within ArcMap. I'm aware that QGIS has the Atlas Plugin, although I'm not as familiar with its capabilities but would happily use it as an alternative.


How can you create a depth raster from a DTM and a polygon in QGIS?

I have a DTM raster and a flood outline in polygon format - how do I create a depth raster for the flood outline (i.e. flood level minus DTM level)?


Looking for system specifications for a modern GIS workstation

My current employer has given me the opportunity to design my own GIS workstation. The problem is I have no idea what is important when it comes to hardware! I have been given a datasheet from Fujitsu which includes the latest CELSIUS M720 series workstations but the configuration includes about 30 different processors, graphics cards, and hard drives. I will be working a lot with the ESRI desktop suite (spatial analyst and 3D analyst) and doing a fair amount of raster processing (DEMs, visibility analysis, etc.), but I also use the GDAL library, SAGA-GIS, Quantum GIS and others. Does ArcGIS 10 actually use multiple cores? I have the choice from 4 cores with 8 threads up to 8 cores with 8 threads - 16 threads? Is SSD SATA III better than HDD SATA III? The Celsius R920 supports up to 512 GB RAM. But what do I need? It will probably be the last time I get a new workstation for the next 5 years, so I am curious what others have ordered or feel will be necessary for the coming years. Obviously the cloud isn't taking over as quick as they say!


System Requirements for Quantum GIS

> **Possible Duplicate:**   >  Looking for system specifications for a modern GIS workstation I have been using the program to make maps to include in my reports over the past year or so. It is excellent. However, my computer seems to not be able to handle all the data sometimes, and the program crashes. It happens when I do too much editing to data points that I have manipulated. My question is what system requirements are recommended to run this program with a lot of data points? Or does it simply depend on what I want to do with the program? I am probably going to get a faster computer primarily based on wanting to use this program. Thanks Marcus


Which processor for GIS Desktop?

> **Possible Duplicate:**   >  Looking for system specifications for a modern GIS workstation I am picking a new pc configuration for my GIS workstation. I use ESRI ArcMap and ArcCat geoprocessing tools, together with PG Admin III, SAGA-GIS. Quantum GIS, Adobe Photoshop/Illustrator as well as command line tools like FW-Tools, Putty and text editors....Not necessarily all at the same time, but more often than not there are between 5 and 10 programs running at the same time. I am configuring a new workstation and have a list of processors which are recommended by Fujitzu...namely the Intel Xeon E5 range. What I am not sure about is how many cores/threads I need. Does each application use a thread?...or even a core? Should I get as many cores as I can afford? any advice most welcome, Yours,


Covering a Polygon with Minimum number of circles with fixed locations

In QGIS, I have a vector polygon of a fixed size and shape. I also have a set of circles of a fixed radius and fixed location, which cover part of the polygon. What I would like to do is remove as many of the circles as possible while leaving as few gaps as possible, a gap being an area not covered by the circles. No circle should have overlap with another circle of more than 90% (or 80% or 70%- I am trying to run a few scenarios) of its area. I would also like to treat circles which go off the edge of the polygon in the same way as overlap, i.e. a circle with more than 90% of its area not covering the polygon would be eliminated. Inevitably there will be times when two or more circles could be chosen, I would like to prioritize them based on specific attributes. Here is a picture of the complete set of circles(pink) and the polygon I need them to cover (blue). ![enter image description here](http://i.stack.imgur.com/dkDlj.jpg) I am new to GIS and any extra elaboration on how to implement is greatly appreciated.


cast from IFeatureClass to ISchemaLock

I'm programming with the ARCGIS Java SDK, in Eclipse. I need to cast an **IFeatureClass** to an **ISchemalock**. This cast is supposed to be allowed. There is an example of this cast in: http://resources.esri.com/help/9.3/arcgisengine/java/doc/b22267cb-642c-11dc-9ca3-0b35f906bb2e.htm However, when I try to perform this cast in my code, an exception appears: _Exception in thread "main" java.lang.ClassCastException: com.esri.arcgis.geodatabase.IFeatureClassProxy cannot be cast to com.esri.arcgis.geodatabase.ISchemaLock_ Are you able to make this cast on an IFeatureClass ? Here is the code to perform this:               //Attempt to acquire an exclusive schema lock on the feature class.     IFeatureClass fc = fw.openFeatureClass("featureclass name");     ISchemaLock schemaLock = (ISchemaLock)fc ;      I would appreciate if you could try this simple line of code in one of your programs. Keep me posted... **EDIT:** I wanted to see if an IFeatureClass was, at least, recognized as an instance of ISchemaLock. Believe it or not, it was NOT recognized as an ISchemaLock! The execution of the code below did not enter the **if** statement.               if (fc instanceof ISchemaLock)      {       System.out.println("all ok");       ISchemaLock schemaLock = (ISchemaLock)fc;     }      Maybe is there something not correctly installed with Eclipse ?


Styling a WFS layer, how to control title in legend?

I have a problem with styling a WFS in OpenLayers. I tried different ways (from the examples) to style the WFS, but I could only get the default style in orange. Then I found one way, which worked for me. The problem is that in the legend is written "Untitled 1" next to the now correct symbol as with the default style. How can I change this title or remove it? Has anybody have a solution on how to style the WFS with an external SLD? I also tried different examples found on the web but I couldn't get it to work. Here's the code I am using to style the WFS but that displays the wrong title:               var StyleObjekt = new OpenLayers.StyleMap({                                 "default": new OpenLayers.Style({                                                               fillColor: "#fff000",                                                               strokeColor: "#666666",                                 }),                                 "select": new OpenLayers.Style({                                                                 fillColor: "#ffffff",                                                                 strokeColor: "#666666",                                                                 }),                                 "temporary": new OpenLayers.Style({                                                                 fillColor: "#ffffff",                                                                 strokeColor: "#666666",                                                                 })                 });            var gwg_f_wfs = new OpenLayers.Layer.Vector("Flurstuecke_WFS", {                 styleMap: StyleObjekt,                 strategies: [new OpenLayers.Strategy.BBOX()],                 protocol: new OpenLayers.Protocol.WFS({                     version: "1.1.0",                     url: "http://localhost:8080/geoserver/wfs",                     featurePrefix: "NewGeo",                     featureType: "Flurstuecke_GWG",                     featureNS: "http://www.newgeo.de/GeoServer",                     srsName: "EPSG:2398"                 })                           });


How to install PostGIS correctly with OpenGeo Suite 4.0 community edition on Windows 7?

I have installed the OpenGeo Suite community edition, version 3, onto Windows 7 multiple times in the past, with virtually no problems. With the latest version of OpenGeo Suite (4.0.1) (which includes PostGIS 2.1), I am able to install almost all of the components successfully, except PostGIS. With earlier editions, PostGIS was installed as part of the full Opengeo installation process, but now, although PostGreSQL 9.3 is successfully installed, there are no PostGIS or opengeo database templates installed, nor are the spatial functions installed. Nevertheless PostGIS 2.1 is installed. It seems that OpenGeo Suite 4 intalls all the component but no longer creates a default database or user account in Postgresql/PostGIS. What steps are necessary to configure OpenGeosuite (windows installation) in order to be able to fix this issue? thanks


Extract points from regular data.frame with coordinates based on SpatialPolygonsDataFrame in R

I have a regular data.frame called datos.merged > class(datosmerged)               [1] "data.frame"      The last two columns are the corresponding geographical coordinates > datosmerged[1,c(23,24)]                  XCOORD   YCOORD     1 2473520   668648.7      I have a file called shape > class(shape)               [1] "SpatialPolygonsDataFrame"     attr(,"package")     [1] "sp"      I would like to use this very simple shape to subset the previous data.frame. How can this be done?


How to create links between points and nearest vertex of polyline feature class?

I have a point feature class and for each point I want to find the nearest vertex on a polyline feature class. I then want to create a line between each point and it's nearest vertex on the other layer. How is this done? If it can't be done with a polyline layer, I can also use a point layer in lieu.


map() Function and IE8

I am using Javisantana example of Autocomplete an everything works inside my map, until I tested my map with IE8 and experienced this error:               Object doesn't support property or method 'map'       After some research it seems as if IE does not support that function some of the hacks/patches that I have tried fix that error but seem to cause errors within cartodb js. Has anyone ran into this problem? If so any tips or advice on a work around?


does cartodb works offline

i had installed cartodb on a local machine. It works fine. But when i tried to go off-line, and then i was able to login, nothing can be done. It seems page is looking for some on-line data's. Is it possible to completely work **cartodb in an off-line mode**


How to make a freehand polyline option from drawing toolbar enabled throughout in arcmap

Following code selects freehand polyline as default draw tool on first go: final ICommandBars commandBars = app.getDocument().getCommandBars(); final UID windowId = new UID(); windowId.setValue("{57E77868-D015-11D2-9F31-00C04F6BC8DD}"); final ICommandItem commandItem = commandBars.find(windowId, false, false); app.setCurrentToolByRef(commandItem); But once we finish with drawing first element, it goes back to its original settings i.e. it starts showing 'Polygon' as default draw option. Is there a way to set drawing tool to freehand polyline option throughout application execution in Arc map Regards, Viji


pg_escape_string() equivalent?

A client is hosting one of my maps on a server that does NOT have a postgres module installed, and will not install it (nor anything else... ) . We have been told also by the site owner that we can't move it due to wanting to retain the IP etc for seo (etc etc). so obviously my first call to escape a string for CartoDB results in "Call to undefined function pg_escape_string() in /usr/ etc etc... Is there a preferred method for circumventing this? it looks like there may be some remedy in PEAR , but it would be very preferable to not have to depend on being able to do /anything/ on this particular server.


NDVI time series in R -- how do I create a time series object from a data frame?

I am doing a time series analysis of NDVI using the bfast package in R. Bfast requires a time series (TS) object. I have successfully created a data frame with the NDVI at various point locations for tiffs in a given directory (code for this is at the bottom of the post). Problem 1: the code does not extract the date of the images. Landsat file names give year and Julian day starting at the 10th character (YYYYDDD) so this should be possible. But how? And will the Julian day be recognized as a date? Problem 2: once the date is extracted, how do I create TS object from the data frame? Thanks for any tips! I'm new at R, so be gentle. Code for extracting values at point locations:               require(raster)     require(rdgal)     require(sp)     setwd("C:\\Users\\Work\\LSTest")     pts <- readOGR(getwd(), "MyPoints")     LS.files <- list.files(getwd(),pattern="\\.tif$", full.names=FALSE)     India.LS <- stack(LS.files)     pts@data <- data.frame(pts@data, extract(India.LS, pts))


Convert an .xlsx datasheet with point feature data into a POINT feature class

I am trying to create a script tool to automate this process, eventually it will be used for batch processing but at the moment I can't even get test data to work. This is pretty basic at the moment, and I know it can be done through ArcGIS, but I would like to develop a python script for processing. My code is as below: It completes without error, but I am expecting the creation of a feature layer (even better would be a feature class with POINT geometry). I can't find or access the output of the SaveToLayerFile step. I apologise for the messiness of the code, I've been trying to beat this for 2 days now.. Any help would be appreciated.               #Import libraries     import arcpy     import os     import sys     arcpy.env.overwriteOutput = True          #Set workspace     home = arcpy.env.workspace     arcpy.env.workspace = arcpy.GetParameterAsText(0)     mxd = arcpy.mapping.MapDocument("Current")     df = arcpy.mapping.ListDataFrames(mxd, "*")[0]               #Create spatial reference object     prjFile = arcpy.GetParameterAsText(2)     spatialRef = arcpy.SpatialReference(prjFile)     arcpy.AddMessage('Setting Spatial Reference')          #Set Excel Datasheet for Migration     TEST = arcpy.GetParameterAsText(1)          try:     # Process: Table to Table Conversion with destination variable         FarmerTable = 'Farmers1'         arcpy.AddMessage('Migrating Data...')              arcpy.TableToTable_conversion (TEST, home, FarmerTable)     except:         arcpy.GetMessages()               # Set local variables     in_layer = "FarmerLyr"     out_layer = "FarmerLyr.lyr"     out_layer0 = "FarmerLyr"               try:        #Make XY Events of each Farmer         arcpy.AddMessage('Creating Features...')              result = arcpy.MakeXYEventLayer_management(FarmerTable, "X", "Y", out_layer0, spatialRef, "")         result.getOutput(0) = out_layer0         arcpy.AddMessage(arcpy.GetCount_management(out_layer0))         arcpy.AddMessage('Successfully Created Features')              # Execute SaveToLayerFile         arcpy.SaveToLayerFile_management(out_layer0, out_layer, "ABSOLUTE")     except:         print arcpy.GetMessages()


Is there demand for a higher resolution alternative to GPS tracking in the GIS market?

If there was an alternative to GPS that provided higher resolution tracking, what applications would there be for it in GIS? Note that this means higher resolution, but not necessarily higher accuracy. For instance, the absolute position (lat/long) would not be more accurate than GPS, but the relative movement during tracking could be within a few inches. This means you could get a very detailed track recording of a device's movement, allowing you to see small changes in position or elevation. If you were to carry this device while hiking a trail, you could see elevation changes with 1-5 inch resolution. I imagine this might be useful for GIS, but I am not familiar with the market. Any insight would be appreciated.


Thermal Map Data

I am curious about making a thermal map to display where the hot spots and where the cool spots in a community are. Obviously the roads will be hot and the shaded areas with trees will be cool but this is more to brush up on my skills. My main question is how would I go about collecting this data for an area where none currently exists? It will be some form of raster data but I am lost from there on out. Does this relate more to remote sensing? If anyone could help me get the ball rolling I would greatly appreciate it.


Iterating using search cursor and layer selection performance?

I'm trying to use the results of a search cursor to iteratively make selections on a feature class and then create new feature classes (so you get a new feature class or shp for each unique value in a field... kind of an explode by attribute tool). _As an example - if you had a feature class containing 14763 records and one of the fields had 18 unique values in it, after running the script you would have 18 new feature classes in your gdb._ So far I have this:               def unique_values(table, field):     with arcpy.da.SearchCursor(table, [field]) as cursor:         dict = sorted({x[0] for x in cursor})         print type(dict)         for i in range(len(dict)):                  whereD = dict[i]             fname = whereD[0:5] + "_" + whereD[-4:]             arcpy.MakeFeatureLayer_management(table, fname)             query = "relatedBirds_Species = '" + whereD + "'"             arcpy.SelectLayerByAttribute_management(fname, "NEW_SELECTION", query)             arcpy.CopyFeatures_management(fname,fname+ "_c")      This works but it's incredibly slow when it runs. Does anyone know of a quicker way to get the same results? **Requested edits below** _The Data_ Geodatabase workspace The feature class I'm running this with has 18 unique values in the field in question. Total number of records in the original feature class is 14,673 _Timings_ As the code is above = 2' 35.04" With the {where clause} added to MakeFeatureLayer (w/out copy feature either)= 1' 27.75" (thanks Stephen) With Select_analysis = 2' 10.58" I think most of the time was being taken to actually draw these in Arc as I was running this directly from the python window so I ran it again from cmd with arc closed. Stephen wins at 57.93" (with copy feature put back in) running from cmd.


How to parse xsd:schema file using openlayers?

I got layer column names using this thread. But I am stuck parsing it. How to parse this xml using openlayers?               <xsd:schema elementFormDefault="qualified" targetNamespace="http://www.openplans.org/topp">         <xsd:import namespace="http://www.opengis.net/gml" schemaLocation="http://localhost:11032/geoserver/schemas/gml/3.1.1/base/gml.xsd"/>         <xsd:complexType name="statesType">             <xsd:complexContent>                 <xsd:extension base="gml:AbstractFeatureType">                     <xsd:sequence>                         <xsd:element maxOccurs="1" minOccurs="0" name="the_geom" nillable="true" type="gml:MultiSurfacePropertyType"/>                         <xsd:element maxOccurs="1" minOccurs="0" name="STATE_NAME" nillable="true" type="xsd:string"/>                         <xsd:element maxOccurs="1" minOccurs="0" name="STATE_FIPS" nillable="true" type="xsd:string"/>                         <xsd:element maxOccurs="1" minOccurs="0" name="SUB_REGION" nillable="true" type="xsd:string"/>                         <xsd:element maxOccurs="1" minOccurs="0" name="STATE_ABBR" nillable="true" type="xsd:string"/>                         <xsd:element maxOccurs="1" minOccurs="0" name="LAND_KM" nillable="true" type="xsd:double"/>                         <xsd:element maxOccurs="1" minOccurs="0" name="WATER_KM" nillable="true" type="xsd:double"/>                         <xsd:element maxOccurs="1" minOccurs="0" name="PERSONS" nillable="true" type="xsd:double"/>                         <xsd:element maxOccurs="1" minOccurs="0" name="FAMILIES" nillable="true" type="xsd:double"/>                         <xsd:element maxOccurs="1" minOccurs="0" name="HOUSHOLD" nillable="true" type="xsd:double"/>                         <xsd:element maxOccurs="1" minOccurs="0" name="MALE" nillable="true" type="xsd:double"/>                         <xsd:element maxOccurs="1" minOccurs="0" name="FEMALE" nillable="true" type="xsd:double"/>                         <xsd:element maxOccurs="1" minOccurs="0" name="WORKERS" nillable="true" type="xsd:double"/>                         <xsd:element maxOccurs="1" minOccurs="0" name="DRVALONE" nillable="true" type="xsd:double"/>                         <xsd:element maxOccurs="1" minOccurs="0" name="CARPOOL" nillable="true" type="xsd:double"/>                         <xsd:element maxOccurs="1" minOccurs="0" name="PUBTRANS" nillable="true" type="xsd:double"/>                         <xsd:element maxOccurs="1" minOccurs="0" name="EMPLOYED" nillable="true" type="xsd:double"/>                         <xsd:element maxOccurs="1" minOccurs="0" name="UNEMPLOY" nillable="true" type="xsd:double"/>                         <xsd:element maxOccurs="1" minOccurs="0" name="SERVICE" nillable="true" type="xsd:double"/>                         <xsd:element maxOccurs="1" minOccurs="0" name="MANUAL" nillable="true" type="xsd:double"/>                         <xsd:element maxOccurs="1" minOccurs="0" name="P_MALE" nillable="true" type="xsd:double"/>                         <xsd:element maxOccurs="1" minOccurs="0" name="P_FEMALE" nillable="true" type="xsd:double"/>                         <xsd:element maxOccurs="1" minOccurs="0" name="SAMP_POP" nillable="true" type="xsd:double"/>                     </xsd:sequence>                 </xsd:extension>             </xsd:complexContent>         </xsd:complexType>         <xsd:element name="states" substitutionGroup="gml:_Feature" type="topp:statesType"/>     </xsd:schema>


Problem in Spatialite database creation

When I try create a new DataBase, Spatialite plugin return the dialog error box:               The SQL query seems to be invalid. no suck table: views_goemetry_columns      When I click in OK button, a second dialog box (Transform DataBase into SpatiaLite) is open with the text :               This DataBase seems to be a valid SQLite DataBase but not a valid SpatiaLite One     Would you like QSpatiaLite to automatically transform it to SpatiaLite (information won't be lost) ?      When I click in the YES button, a new dialog error box open with the text:               Unable to convert Database to SpatiaLite      When I click in OK button, another dialog error box open with the text:               The SQL query seems to be invalid.      The DataBase was set up, but I cant import spatial layers. How resolve this problem?


How to convert a point grid to a raster?

If I load the original .csv file, I then have a map of the African continent represented by 14069 data points covering it at 40x40km resolution. Each data point shows solar radiation means across the year, (jan, feb, etc). What I want to be able to do is convert that vector into a raster using GRASS. What I have managed to do already, is import the .csv (tab delimited) file into the GRASS mapset using v.in.ascii. However, when I try use v.to.rast.constant to convert that vector into a raster, I get a solid square covering Africa, that cannot be edited. I.e. there is no data from the vector within the created raster. What am I doing wrong? This work is for my masters thesis and I'm constantly working on it. Thanks!


Repaint seems not working when called multiple times from the same method

I have the following scenario:   `   displaySelectedFeatures(IDs);   then some code... and then:   displaySelectedFeatures(features);    Display selected features looks like this:               public void displaySelectedFeatures(Map<Color, Set<FeatureId>> features) {              Style style;              style = createSelectedStyle(features);              Layer layer = mapFrame.getMapContent().layers().get(0);         ((FeatureLayer) layer).setStyle(style);              mapFrame.getMapPane().repaint();         }         What I need is the following. First IDs object should be drawn and then after some time e.g. 5 seconds, object - features should be repainted. What happens it that only the IDs is drawn while features(second repaint) is omitted. How to solve this?


How to install PostGIS 1.5 on Mac OS 10.5?

I am running PostgreSQL 9.0 on Mac OS 10.5. I've installed postgis- pg90-1.5.3-1-osx through "Stack Builder" but still don't see postgis in my list of databases when I open PgAdmin. When I open my Qgis 1.7, I can connect to my postgres database and display all tables with no spatial geometry. However, when I try to query a non-spatial table, I get the following message "Your database has no working PostGIS support". I appreciate your ideas to get postgis to integrate into postgresql. Thanks! * * * Thanks for your ideas! Unfortunately, the first link (http://www.kyngchaos.com/software/postgres) is limited to "...Lion and Snow Leopard, 64bit only." I'm using Leopard (Mac OS 10.5). The second link (http://trac.osgeo.org/postgis/wiki/UsersWikiPackages) provides two options. Both of these (http://yum.postgresql.org/repopackages.php) & (http://www.postgresonline.com/journal/archives/204-postgis15-install- yum.html) point to Linux installations. So, I'm still stuck. I've downloaded PostGIS 1.5 using Application Stack Builder 3.0.0 from PostgreSQL 9.0. It installs, but does not run. I appreciate any other suggestions.


How to convert WMS to SHP?

I'm looking for tool that will be convert WMS to vector layer. Are they ready tool?(for example in python)


Which parameter values should I use with QGIS Heatmap plugin?

I've read all the posts and docs about the heatmap plugin but can't seem to get it to work. The map is the size of the United States and has a little over 2000 points. I've tried small radi (1 meter) and large radi (5 million meters), but still get a red hot map. I'm using QGIS 1.8 and WGS84. The data was initially in a spreadsheet, but I saved it as a .shp vector layer and deleted the old csv layer. The goal is to convert that to contour lines to map density of these points like in this Mapbox blog post.


looking for an esri compatible javascript engine

I work on an esri heavy team, i'm looking for a mapping engine that will give me flexibility to build on top of and leverage only the pieces I need and or add more as the app grows. I've been using ESRI's API to great success for my indiviual code but leveraging their codebase into a reusable solution grows too large. We have arcgis server 10.1 and create featureserver and mapserver services. Along with those I will consume data thats created from various tabular locations and sent to my front-end via JSON/Ajax. I've started on leaflet, and have began testing on esri-leaflet to some success, dynamic layers are drawing slowly, but we are currently upgrading to 10.1 so that may wash out. Goal: find a lightweight javascript based mapping engine that will handle ArcGIS feature&mapserver services. Along with the ArcGIS services the engine needs to handle JSON data. Geocoding with the service of my choice would also be great(bing, google, esri, or internally created) Thx


Is there a way to reference a selection in model builder or python?

Is there a way to reference a selection in model builder or python so that I can edit my python script to only run if there is a selection. And even better if there is only one row selected? I made a model to make a selection with python script to ask for variables, make a selection based on those variables, and then zoom to the selection. When the tool does not make a selection for whatever reason it takes forever to finish and zooms to the entire extent of the map or crashes ARCmap entirely.


How to create a text box with inverted rounded corners

I just saw this post on the tex SE site and wondered how I would do the same in ArcGIS. I created a text box in a layout and drilled down in all of the dialog boxes but couldn't see anything similar. Sure, it's easy to create a rounded box, but how do you invert the corners? One solution I had was to manually draw a rectangle with circles on the corners. ![enter image description here](http://i.stack.imgur.com/td87c.png) Then, I used the _subtract_ tool on the Graphics Operations context menu for each of the corners. ![enter image description here](http://i.stack.imgur.com/eUREn.png) The result was satisfactory: ![enter image description here](http://i.stack.imgur.com/5icEm.png) I could easily put some text in there and it would look fine. However, if the amount of text changes and I need to stretch the box, the corners stretch as well: ![enter image description here](http://i.stack.imgur.com/BB8Ly.png) Does anyone know of a way to create such a frame around text in the layout in ArcGIS? I'd be happy to know if it's possible to do this using Python.


In ArcGis 10, spatial join of a hot street value to a fishnet cell is counting multiple values for same street segment

So, I am working on a fishnet, so that we can attribute multiple values to points within different cells (i.e., ranking light poles as optimum for traffic cameras based on the number of speeding tickets, accidents, etc) so that all poles in a high value cell would be ranked as prime candidates. To do this, we have assigned numerical counts of accidents to street segments (sections between intersections) and are attempting to join those values to the fishnet through a spatial join. The issue is, if a street segment runs through the cell, and has a value of 5, it is counted 15 times (or whatever for each cell) and the cell then has a value of 75 instead of 5). Any way to fix this and get one value from each unique street segment in the cell? (averaging as a run around doesn't work, as some cells have different street segments and if one street has a 0, we don't want that to discount the value of the other street that may have a higher value.) We've tried dissolving, interection, multi to single, and none of them seem to make each street segment unique or with only the requisite value display once. Thanks. Chris NYC as requested, here is a shot of the image. the issue is, and this example actually accidentally illustrated it clearer than I'd seen, is that for the horizontal street near the middle, with a "3" (there are only two dots because one has 2 values) looks correct, but if one tried to join that "3" value to a cell intersecting that street segment, it would likely join the value "45." This is demonstrated a street over where one street has a "3" a "0" and a "7." This seems to demonstrate the single street is actually 3 segments to arcgis. I can't figure out how to make, for example, that street equal "10" and then join it's value to a cell so the cell would only have "10" as the joined value, not "30" or "150" or however many times it joined. ![map example](http://i.stack.imgur.com/GOIX9.png)


Custom python scripts do not appear in SEXTANTE QGIS

Custom python scripts do not appear in Scripts menu in SEXTANTE plugin unless they are either empty or are as simple (and useless) as:               ##foo=group     ##input=raster      OS is openSUSE 12.3; QGIS 1.8 and Master both demonstrate the same behaviour. I was told that the issue may be caused by insufficient permission rights for the scripts folder. The path to the script's (for Master in this case) is that is defined by default: `/home/username/.qgis2/processing/scripts`. For both 'processing' and 'scripts' folders I have the following permissions: `drwxr- xr-x`. For the scripts in 'scripts' folder the permissions by default are `-rw-r--r--`, but even if I enable execution rights manually to be `-rwxr- xr-x` this doesn't affect their appearance in SEXTANTE scripts. So I suppose that permissions doesn't seem to cause trouble here... Please advise me how to overcome this issue.


ArcObjects SDK: Textured MultiPatch: problems with large texture images

I am using IGeneralMultiPatchCreator to create a textured MultiPatch. This works fine for small texture images, but fails for larger images (15 MB .png). In ArcScene the MultiPatch is displayed in grey, the same way it would be if the "Disable material textures" checkbox on the Rendering tab of the Layer Properties was checked. I have tried to increase the max texture memory for the layer:               I3DProperties2Ptr prop3D = layer;     long maxTextureMemory = 128*1024*1024;     prop3D->put_MaxTextureMemory(maxTextureMemory);     prop3D->Apply3DProperties(layer);      but this does not seem to have any effect :-(. I have also tried to split the MultiPatch into parts where each part have a smaller texture image, but this does not seem to work either.


Creating layer from in-memory shape data using ArcObjects?

In my current project, a feature layer is created based on shape data generated after some complicated surface analysis procedure. The sequence right now is like [Surface Analysis Algorithm] -> [array of xyz data] -> [shp file written to disk] -> [create feature class from shape] -> [create feature layer from feature class] I'm on the process of optimizing the project, so I wondered if there was any way to display directly the in memory shape data in a layer, bypassing the read/write to hard disk?


Is it possible to show Pie charts in QGIS Legends?

I would like to display some pie charts on my map. They show in the map frame perfectly but in legend on the print composer I am unable to show the legend for the pie charts. QGIS just shows the style of the point underneath the pie. Is it possible to get these shown? Using QGIS v1.8 from OSGEO4W on Windows XP


Geometric Network error in File Geodatabase<Network Dataset

I need to create a geometric network in my network dataset which resides within my file geodatabase. I have been receiving the following error: Cannot create geometric network. _The selected feature dataset does not contain any feature classes which can participate in a new geometric network._ This error should not come up because there are tons of functional feature classes in the network data set. Any workarounds or solutions?


How to extract vegetation polygons from Imagery using Python?

I have been trying to roughly polygonize vegetation areas from Imagery (RGB). I tried using Iso Cluster Unsupervised Classification (Spatial Analyst) but it is still too detailed for my purpose. I just need to roughly draw polygons wherever I have dense vegetation as can be seen here:![enter image description here](http://i.stack.imgur.com/6ac4u.jpg) Any idea/ thought to put me on the right track is much appreciated.


Plot the nodes of the shortest path found by dijkstra algorithm

I am using qgis and pgrouting plugin to find the shortest path between two nodes and display the route. Currently, the plot shows the shortest path. However, I also want to know which nodes the path passes through. So is it possible to actually get this information?


QGIS 2.4 plugin PyQT qWebview javascript interface problem

I have a QGIS plugin which contains a `qtwebview` and Javascript communication between the page and the plugin. It is working perfectly with QGIS 2.2 but not with 2.4 (both QGIS versions share the same plugin directory), The javascript side never gets called with `evaluateJavaScript()` and the python side never gets invoked from Javascript either. What could explain that it is working in 2.2 but not in 2.4?               def on_loadFinished(self):          ...               def init_Clicked(self,iface):         self.ui.webViewResult.page().mainFrame().addToJavaScriptWindowObject("QGISInterface", self.moTAZJsInterface)         self.ui.webViewResult.loadFinished.connect(self.on_loadFinished)         self.ui.webViewResult.load( QUrl(config["TAZ_URL_INTERFACE"]))         self.ui.webViewResult.show();


How to use pgsql2shp when PostGIS extensions installed in a separate schema as geometry

I have the PostGIS extension (postgis functions and spatial_ref_sys table) added to a specific schema called postgis. I'm wanting to run the pgsql2shp utility with a select query on geometries stored in other schemas that do not have postgis functions available within their respective schemas. Normally from psql I would first set the search path to contain both the postgis schema and the schema with the tables I'm querying on, but I'm having trouble figuring out how to do the same thing with the pgsql2shp utility. How would I go about setting a search path or having it look in multiple schemas so it picks up the postgis functions? Example Query (field_polygons is in different_schema):               pgsql2shp -f shapefile -h localhost -u user -P password gisdb "SET search_path TO postgis, different_schema; SELECT the_geom FROM field_polygons WHERE id = 1"      Results in error:               ERROR:  function postgis_version() does not exist     LINE 1: SELECT postgis_version()


How to query a PostGIS database through Geoserver and OpenLayers

I've built a web application using OpenLayers, which is connected to my PostGIS database through Geoserver. Now that I'm able to successfully display all of the layers, I'd like to add querying/interactivity capabilities to the map. For instance, I'd like to create a dropdown menu of checkboxes so the user could "Choose X who have done Y", another dropdown to "Choose X who haven't done Y", etc. I'm having trouble figuring out how exactly to add in this functionality. Does anybody know any good resources/examples for how I might go about doing this? Thanks for your help.


How to download hospitals and restaurants from OSM in shp format?

I need information on hospitals and restaurants in Guatemala. How can I download this information from OpenStreetMap? I want to use the data in Arcgis 10.1


Value of a Python programming course

I am considering upping my programming skills by taking a python geoprocessing course through ESRI. Before I part with the hefty registration fee, I would like to get your opinion on the value of 24 hrs of instructor led python programming. I am primarily interested in learning the following:   * For those of you who have taken a course such as this, did you walk away with applicable skills or more theoretical knowledge?   * For those of you who have recently been in the job market, did you find employers looked favorably on having a class such as this on your resume?   * For employers, does a specialty course like this really matter, or do most GIS analysts/specialist learn python programming on their own? Thanks for the advice.


Mapping zipcode to legislative district

Did anyone work with SunLight foundation data to map the zipcodes with legislative district data. I have to work with Illinois zipcodes only.


Do Different ArcGIS for Desktop License Levels Process Faster?

At my job, I have access to all three levels of ArcGIS licensing on my computer. My question is: Does the license you are using in ArcMap make a difference on the speed of ArcMap? I'm pretty sure that the drawing of the map isn't affected, based on tests I've done myself. But, it really does seem like some of the ArcToolbox tools do work faster in ArcInfo. I'm just trying to either validate my suspicious or prove I am going a bit crazy. Thanks for any replies!