title: Visual Studio 2010 templates for ArcGIS missing

text: It seems that I do not have ArcGIS templates in VS 2010 templates folder - > \Program Files (x86)\Microsoft Visual Studio > 10.0\Common7\IDE\ProjectTemplates\CSharp\ and as a result ArcGIS templates do not show up as installed templates for New Project. To get these templates will I need to reinstall VS or can I get the zipped ArcGIS templates from somewhere? Thank you, -Dee


title: Color matching between Google Earth and Google Maps (same KML)

text: I'm working on a custom Google Map using API v3. I'm creating my KML files in Google Earth, uploading to Maps and then accessing the KML from there. I'm having problems matching polygon colours to the defaults in Maps - everything looks fine in Earth, but when viewing it on google maps, the colours have shifted. **How the KML looks on Google Earth:** ![enter image description here](http://i.stack.imgur.com/ewj0g.jpg) **How the KML looks on Google Maps:** ![enter image description here](http://i.stack.imgur.com/NX3j3.jpg) You can download the KML here.   KML file. Is anyone can help me with this issue?


title: Fusion Tables API advanced visualizations

text: From the google example; i want to replace the column bar chart with a dynamic pie chart. How do i do that. What codes do I need to change?


title: ArcGIS 10: How to programatically project coordinates?

text: Let's say I have coordinates defined for SRID 4686, is there a way to programatically transform them to SRID 4682?


title: Is it possible to simplify orthogonal polygons with OpenCarto Java Library?

text: I want to know if simplification of orthogonal polygons can be done with OpenCarto Java Library. If yes, does it have any API documentation? Cheers, Sanjeewa


title: Why do FusionTable layers sometimes not work / slow down?

text: Fusion tables layers are great. I have uploaded some data and polygons of areas in Europe and most of the time they respond very fast. First thing in the morning (after prolonged inactivity I assume), the layer will sometimes take some time (about 5 minutes) before it appears. Occasionally (about once a day), even whilst I am looking at and using the map layer, it drastically slows down, or stops working all together. This lasts for roughly half an hour. Are these common and unavoidable problem of using a service I neither control nor paid for, or is there something I might have done wrong? Is anyone else experiencing issues? Thanks


title: How to draw Sprite on Graphics layer?

text: I want to draw sprite on Graphics layer like this:                           var sp:SpriteVisualElement = new SpriteVisualElement();                 sp.graphics.beginFill(0x00ff00,1);                 sp.graphics.drawRoundRect(10,10,100,100,150,150);                 sp.graphics.endFill();                 var myGraphicMarker:Graphic = new Graphic();                 myGraphicMarker.symbol.draw(sp,geom,null,map);                 graphicsLayer.add(myGraphicMarker);      It complaints at line myGraphicMarker.symbol.draw(sp,geom,null,map); about that TypeError: Error #1009: Cannot access a property or method of a null object reference. Can anyone please guide me what to pass instead of null to make this sprite drawn. Thanks


title: Why do we need Mapserver/Geoserver to present data from spatial database to web browser?

text: I was trying "How I set up GeoServer on an old machine" and By @iant I was setting up Tomcat and all. I suddenly had a question why we need Mapsever and Geoserver? Can these server pull data from database forming layer as user asks and only presenting what user has asked for from database?


title: How to identify Polygon within another Polygon layer

text: I have two polygons files. The first one has several smaller polygons within it eg.a shapefile with all US states, the other polygon file has a single polygon eg. a river basin. I am trying to get all the polygon (states) within the second polygon(river basin). Intersection tool in Arcgis10 clips the first polygon whereas I want all the intersecting polygon to maintain its original boundaries.In this way I will have the exact state boundaries. The basic idea here is to identify all the polygons intersecting another polygon layer. thanks in advance.


title: Where can I find a high quality real globe of a decent size without axis?

text: I'm trying to find a physical globe of a high map quality that _doesn't have any axis or bearings_ as the axis mounted on standard desk globes completely obscures the arctic regions. High-quality map globes seem to be available but repeated searches failed me to find a globe model of a decent size (say 12'' and up) with an Atlas quality map and _without_ axis, so that one can view the earth from any angle and the arctic regions are fully mapped and easily visible. Do you know whether/where such a product exists? * * * Edit: Magnetic globes won't do, as they are too small and the one example from an answer also obscures the polar regions through its magnets. What I'd really like is a complete sphere without any knobs or visible seams (I don't mind seams if they are well covered) that can be placed in a holder. Inflatable globes seem quite OK, but their map-quality is very low as far as I could find. Hmmm ... it seems these would go in the right direction (as far as I can trust the images): http://www.worldglobeu.com/movaglobes \-- but they only go up to 8,5'' the maps look so-so and they are insanely expensive for the size and map-quality they offer. This one seems to get the axis/seam thing halfway right, but the mapwork is more artistic that accurate: ![Blue Green Earthsphere](http://i.stack.imgur.com/n3BWM.jpg). (Note that if that globe offered an Atlas quality map it's price tag _might_ even be acceptable at 30''.)


title: 3D roof polygon creation from polylines in ARCGIS 10.1

text: I am Working in arcgis 10.1,I would like to create roof polygons from lines, the lines are captured with photogrammetry technique. As per my requirement I have found tool "3D Building Generator Tool", but tool is based on VBA language. Please can anyone guide me in arcgis 10.1 to create Buildings roof polygon using with Arctoolbox or is there any alternative process. Thanks Venu


title: Missing lxml in QGIS 2

text: I made a QGIS plugin in QGIS 1.8 doing a lot of XSLT transformation and other XML handling with the library lxml. After upgrading to QGIS 2.0 this plugin does not work anymore:               "ImportError: No module named lxml.etree"      What libraries do one use to handle XML and XSLT in Python in QGIS 2.0? Would it be a mess to try and get lxml working, in regard to distributing the plugin to other computers?


title: How to acquire schema lock is one already exists?

text: I'm trying to delete a feature class from a feature dataset or to add a field to it, but I keep getting an error message (pics attached below). I have imported a geometric network I've created from one PGDB to another. However couldn't add any fields. I'm receiving the following error: ![enter image description here](http://i.stack.imgur.com/04ewc.png) Than I planned to delete the feature classes that were imported - I've received the following error: ![enter image description here](http://i.stack.imgur.com/ujWQD.png) I have the feeling that both errors are related but couldn't figure out the problem from ESRI help or from other posts in gis.stackexchange.


title: How can I improve the performance and start-up time of arcgis 10.x?

text: ArcGIS 10.0 take about 3 minutes to start. ESRI recommends doing these steps (http://support.esri.com/cn/knowledgebase/techarticles/detail/31672) but does not explain the consequences or adverse side-affects with respect to other installed applications. I have also heard from a recent ESRI Seminar that decreasing the number of installed printers on the System will also improve start-up performance. Can anyone comment. My system is Win7 64bit 4gb Ram. No Licence server.


title: Postgresql Trigger to insert records in my table from another one

text: I spent few days trying to resolve this question but it was impossible. I have two tables: Table1.               CREATE TABLE point     (       id_point integer NOT NULL DEFAULT,       layer text,       elevation numeric,       geom geometry(POINT,25830)     )      Table2.               CREATE TABLE point2     (     id_point2 integer NOT NULL DEFAULT,     layer text,     elevation numeric,     geom geometry (POINTZ, 25830)     )      In this situation I want get the values from **point.elevation,point.geom** and insert to **point2.elevation, point2.geom** each time I insert a new row in point table. Also to say I need transform the point.geom value before to insert into point2.geom like this:               geom=st_setsrid((st_makepoint(st_x(point.geom),st_y(point.geom),point.elevation)), 25830)      I tried to create this trigger but I did not get the expected result, my trigger inserted all my records from point to point2, but I need do it for each new record from point table: The function:               CREATE OR REPLACE FUNCTION function_prueba() RETURNS TRIGGER AS $$     DECLARE     BEGIN          INSERT INTO              point2(elevation,geom)             VALUES                 NEW.point.elevation,                  st_setsrid                     (                         (                                st_makepoint(st_x(NEW.point.geom),st_y(NEW.point.geom),NEW.point.elevation)                         ), 25830                     )                     FROM                              point                     ;          RETURN new;     END;     $$ language plpgsql;      The trigger:               CREATE TRIGGER trig_geotrans          AFTER INSERT ON point          FOR EACH ROW          EXECUTE PROCEDURE function_prueba();      How can I resolve this question? Thanks.


title: how create a dyamic layer in openlayer map

text: I am using geoserver, Postgis(database) and OpenLayers for map. Now I want to create a layer (WFS, WMS) dynamically. The goal is to allow the user to create any number of layers from front end without hardcoding in (map.js) file. Is it possible?


title: How to write script for Line of sight and Viewshed analysis in Global Mapper

text: I wanted to write script in Global Mapper for Line of Sight and Viewshed analysis. I tried to calculate using Viewshed and Global Mapper Tools but the issue I am facing is that analysis is user specific, meaning we have to give observer point manually. Every time it is very difficult to give observer point. I wanted to know whether there are any other methods to give multiple points or point to a .shp file as input point and DEM as input raster. I am very new to Global Mapper. I am reading about scripting in Global Mapper but I didn't understand how to use it. My thinking is to write script which will give multiple points as input (similar to batch processing) which will give me output of the viewshed of each points. Is it possible? Does anyone know about it? Thanks in Advance


title: Geoserver Caching does not save tiles

text: I want to cache a layer from Geoserver. I use Geoserver 2.5.1 on Ubuntu 13.10.   * In the `web.xml` file under `tomcat7/webapps/geoserver/WEB-INF` I added the following part: ` <context-param> <param-name>GEOWEBCACHE_CACHE_DIR</param-name> <param-value>/opt/gwc</param-value> </context-param> `   * I enabled `Enable direct integration with GeoServer WMS` under Caching defaults.    * I use the sample data `tiger:giant_polygon` where `Create a cached layer for this layer` is enabled After opening the tile preview with `EPSG:4326 /jpeg` for the layer, the folder where my tiles are supposed to be saved is still empty. Also the `DISK USED` colum shows `N/A`. ![enter image description here](http://i.stack.imgur.com/BDZbb.png) What is missing?


title: Keeping up to date with OSM changes in pgRouting

text: I'm working on a routing application using OSM and pgRouting. I'm only interested in a subset of data, mostly North America (performance was difficult for large amounts of data, so don't want to use the whole world). I want to be able to keep up with updates to OSM data, but I haven't seen any place that provides a smaller changeset then for the whole world. Is reimporting the data every time I want to update the only way?


title: Why does QGIS 2.4 crash at every shutdown?

text: I am using QGIS on two Win7, 64 bit machines for years. After upgrading to QGIS 2.4 it crashes on both machines at every shutdown with the message "Minidump...". The crashes occur independant of what I have done before or whatever project I am working on. Until now there seems to be no loss of data, but it is annoying. (I do not have installed the LecoS-plugin on either machine.) Any ideas, solutions or similar problems?


title: How to add a new permanent column that contains result from ST_AREA select queries

text: I am a total newbie here, I have a bunch of polygons in my cartodb map, and I want to create a new column called "AREA" in the table MY_TABLE that contains values from my area calculation query: select ST_Area(the_geom::geography) as area_sqm FROM My_TABLE this statement works, but the column area_sqm that contains the area vlue clears every time I clear the query, how can I have this permanently as part of my table so I can display the area value in the info table of my map?


title: Go to XY tool in an ArcGIS online web map and/or application

text: Has anybody created any language to attach to an ArcGIS online web map and/or application for a "Go to XY" tool?


title: postgresql export database with command or pgAdmin III

text: I want to export my postgresql db with pgadmin.   I search all of its GUI, I also searched the web and I don't find some sueful information.   I know how to import the db into postgresql: create a db with `template0' and then:                c:\program files\postgresql\9.3\bin> psql -h localhost -p 5432                                      -U postgres newDBname < D:\backup.sql      But I don't know how to export it. please help me to export my db.


title: Bootstrap and Leaflet - Grid Issue

text:
Trying to use integrate a simple leaflet map within Bootstrap grid / column. While I am trying to follow tutorials and templates, for some reason the `<div> id="map"</div>`does not seem to be nesting within the column. The column only seems to be 1px in height. In fact, when looking at Firefox Inspector, the full extent of the map window doesn't even appear to be registering as part of the body (apologies for any dubious terminology - very much starting out and probably running before I can walk). My simple code is as follows:

Simple Mapping Page

Some witty lead in text...

Some test text

      I also apply the following CSS:               #map {     width: 100%;     position: absolute;     height: 600px;     margin:0;     padding:0;      border: 1px solid black;     border-radius: 8px;     }      While the second column appears to the right of the map, there is no padding etc.


title: German coordinate system from Lat/Long

text: I have WGS latitude and longitude for data points near Munich and I'd like to use a grid-type coordinate system, like the Ordnance Survey grids in Ireland or the UK. What coordinate system do you recommend? Does Germany have a similar national grid? Then a bigger question, for any location how does one find out what the recommended local grid is? apart from asking here of course!


title: WFS Layer - capabilities document contained no layers

text: I use Quantum GIS 1.8.0 on Windows 7 Professional, 64 bit. I have problem with WFS Layer. I tried tree WFS: 1) Download services - Parcels (Czech Republic) http://services.cuzk.cz/wfs/inspire-cp- wfs.asp?service=WFS&request=GetCapabilities 2) Geonames (Czech Republic) http://geoportal.cuzk.cz/WFS_GN/WFService.aspx?SID=27082012_03561414729376000_PM_021892838861&service=WFS&request=GetCapabilities 3) Example geoserv.weichand.de:8080/geoserver/wfs?service=WFS&request=GetCapabilities Example #3 is OK. Example #2 shows Error "The layer geoportal.cuzk.cz/WFS_GN/WFService.aspx?SID=27082012_03561414729376000_PM_021892838861&SERVICE=WFS&VERSION=1.0.0&REQUEST=GetFeature&TYPENAME=gmgml:GEONAMES&SRSNAME=EPSG:102067 is not a valid layer and can not be added to the map" Example #1 shows Error "capabilities document contained no layers." I don't know, where are the differences between services. Could you help me? Thanks Petr


title: Quantum GIS Standalone Installer Windows folder

text: I have a very simple question for you: what is the "Quantum GIS Standalone Installer Windows" folder? Is it, as usual, something like "c:\program files\quantum gis"? I cannot make a test (at this moment) in a Windows machine. Thank you very much.


title: Retrieving XY from point graphics set in WPF network (arcgis 10.0)

text: I have a WPF application using ArcGIS 10. In the application I can set two point by clicking on map and calculate route between them. every time I click on the map, The application creates a GraphicsLayer and a Graphic element using this code:               GraphicsLayer stopsGraphicsLayer = myMap.Layers["myStops"] as GraphicsLayer;     Graphic stop = new Graphic() [Geometry = e.MapPoint };      My question: I need to get the XY Coordinates of the points clicked on the maps. Does anybody know what elements\attributes\methods to get this information? Thanks.


title: Does QGIS support MT83F

text: Does QGIS support NAD83 Montana State Plane Zone, US Foot? I have vector data in a CSV file that I would like to display, but the coordinate system is MT83F. I don't see a matching CRS. Is there a user defined CRS that is appropriate?


title: How to round a column in attribute table (ArcMap)

text: I want to use the Field Calculator in ArcMap to round an existing column to two decimals. Currently I have a column that is 6 decimal places long and would like to simply round it down to 2 decimals. I had planned on using the Field Calculator (possibly using Python) to do this, but maybe there is an easier way? **Edit** The accepted answer is probably the easiest way to change a single field, however, here is how to do it with the field calculator for both python and VB. VB:               round([column], 2)      Python:               round(!column!, 2)


title: Best Resampling Method for DEM reprojection in QGis

text: I need to reproject a DEM (aster file) in order to have the same CRS as all my other project data - I have different SHP for the same DEM's location but they have a different CRS. Which resampling method would you recommend? I am using QGis.


title: Where is the Road Graph panel?

text: I am trying to use the Road Graph plugin within QGIS 1.8 (running under Windows 7 64-bit). When I click Vector -> Road Graph -> Settings I get the settings-panel: ![enter image description here](http://i.stack.imgur.com/oyeBK.png) But, I don't know how to enable the "working-panel", which should appear on the left side of the main window, according to the producers of the plugin (screenshot of pic #2 is taken from their site): ![enter image description here](http://i.stack.imgur.com/PqLW2.png) Any hints? What am I doing wrong / what's missing?


title: looping the mask in python for arcgis 10

text: I am working with raster data on a worldwide level with a resolution of 30 arc seconds. I need to export the data into a table format but given the high resolution you can not run the "Sample" (Spatial Analyst) command for the entire dataset at once. Therefore, I have created a mask for each country(with The Set Null command). I was trying to write a code in Python, using a loop where each time the mask will become a different zone, i.e. env.mask = zone1, then zone2 and so on. However, I suppose that within a python session one cannot keep on changing the mask using a loop. Could you kindly suggest a way loop the mask (so that I can just have one code for all the countries but my program still runs)?


title: QGIS Attribute Table Not Displaying Data

text: Using QGIS and one minute I could see the attribute table, the next minute I couldn't. Using QGIS on Linux Mint 13. I tried purge/autoclean/reinstall to get a completely new download of all dependencies and no saved config files, but the problem persists. I can't see the attribute tables for either shapefiles or SpatiaLite data (haven't tested PostGIS yet). Note, the underlying data _is_ there. I can open the DBF in LibreOffice and see the expected data, and if I click on a feature with the Identify Features tool, I do see the attributes associated with that feature. Symbolization of the data (e.g. a simple choropleth map) also works based on the underlying data. The attribute table displays all expected columns, but no rows, as if it was a newly created completely empty layer. This was a recent install on a new computer, but for the last week I have been able to add data and view attribute tables as expected. Today I set the configuration option to "Open attribute table in a dock window (QGIS restart required)" and then when I tried look at the attribute table, no data was displayed. I restarted QGIS, still the same. I then turned off the docking, no effect, restarted the computer, no effect. Completely reinstalled QGIS as indicated above, no effect. One obvious gotcha is that "Show selected only" could be checked, but that is not the case here. Has anyone run into this or know how to fix it?


title: Distributed versioning of geospatial data

text: I'm getting started on the design of a new system that would reconcile, merge, and synchronize changes made across numerous geospatial databases of common ancestry, on an ongoing basis. These databases are distributed across numerous different organizations that have an interest in eliminating duplication of effort. I read that OpenGeo is working on a project called GeoGit, which is interesting but seems very much a work in progress. Are there any other tools, projects, or research that would be worth looking into? I don't want to reinvent the wheel, but can't seem to find much out there on this topic.


title: Where can I get road data that includes max speed?

text: Does anyone know how to detect the road speed I. E. If one of our vehicles was driving on a certain road above the speed limit, we would want to know But no idea where to get this data from Any help appreciated Thanks


title: Any open source alternative for ESRI maps for office

text: I am looking for an alternative for ESRI maps for office. I know it is very good product for analyzing data within MS office environment. My requirement: -I have few stores data contains location details. -Trying to analyzing the trend on the map itself in the excel document +Statistical analysis  -Is there any tool/app in the Open office(or MS Office & going through EuroOffice tool but it looks like its paid extension)? -Researching on Google. Any help regarding this will be great thanks in advance !!!


title: How to show number of patches in a neighbourhood?

text: Given one raster file containing only one zone that indicates patches of vegetation in a city and a shape file containing neighbourhoods for the city, explain how you could create a map that shows the number of patches in each neighbourhood. Do you use topology?


title: How does manual Label Placement in QGIS 1.9 work?

text: There are a bunch of questions here related to labeling in QGIS, but it's not clear which version they are referring to, and I know this is something that has changed recently. Hopefully it makes sense to start a new question due to the new labeling functionality. I'm trying to label a road network so that every street is labelled. I have the 'auto' labeling doing a decent job, but it is not placing a label on every street. How can I (with the minimum of effort) _keep the existing auto-placed labels_ (possibly with a few tweaks/manual moves,) but then also go in and manual 'force' missing labels to appear, and position them myself and/or change font size to make them fit? The new label tools in the toolbar are not very easy to figure out, and the 'pinning' concept is lost on me. Can anyone explain the steps?


title: Quick Way to Finish Area of Interest GIS Data Export Requests

text: I would say that after custom map creation, data requests are the most frequent GIS requests I have on a daily basis. Now they don't take long, but they add up and I have much more important work to do sometimes. So I am wondering if there is a quick and easy way finish a request like this. Usually these requests involve a picture like this in my inbox: ![Area of Interest Request](http://i.stack.imgur.com/6jPLm.png) After I figure out where the heck the person is talking about I usually need to clip about 1 to 6 layers in the area of interest for them. Then I need to zip it up and send it back. Like I said, it's not hard whatsoever, but it is time consuming. I've tried using batch Clip, but it seems to take forever and takes a while to set up. I'm open to any suggestions: programming, ArcToolbox, Open-Source, ect.


title: Postgres SSL SYSCALL error: EOF detected with python and psycopg

text: Using psycopg2 package with python 2.7 I keep getting the titled error: `psycopg2.DatabaseError: SSL SYSCALL error: EOF detected` It only occurs when I add a `WHERE column LIKE ''%X%''` clause to my pgrouting query. An example:               SELECT id1 as node, cost FROM PGR_Driving_Distance(       'SELECT id, source, target, cost           FROM edge_table          WHERE cost IS NOT NULL and column LIKE ''%x%'' ',       1, 10, false, false)      Threads on the internet suggest it is an issue with SSL intuitively, but whenever I comment out the pattern matching side of things the query and connect works fine. I understand this may not be the best 'Stack' for this query given that it is not really GIS orientated so just let me know if I should migrate it to where it best fits.


title: VB ArcGIS 9.3 Field Calculator

text: I'm trying to join one column into another. one is "acct id" and the other is "taxid number". Acct is from one county whereas the tax is from another, so there are blanks in the number of parcels/rows. I did it before last week and I had it written down how to do this. What I need to do is move 'acct id' into the 'taxid' column. I remember having " " as part of the vb script. What are the steps involved as I keep getting a blank field when I attempt it. Thank You, Z


title: Querying data linked to a postGIS database

text: My project is coming along nicely! A brief background on my project. I have a world map that is colored to show the country of origin of the students at a college. Each color represents the number students (a range). E.g. any country colored red has 1 to 20 students at the college. When you click on a country, I was successful in displaying data specific to that country in a pop-up. Hopefully I can describe what I want to do next.... I want to add the ability to query the map such that the countries that satisfy the query will be highlighted on the map.   So for example, in my html form I enter > 300... click on submit.... England, India and USA are highlighted. I've figured out that I need to use WFS getFeature but as you can see it's the how that I need help with... Can anyone direct me to an example of the above or give me some suggestions, I would be most appreciative. I am using MapSever, PostGIS and Openlayers Chris


title: Can GeoNetwork store metadata for a PostGIS table?

text: Can GeoNetwork be used to store Metadata for a PostGIS table? Or is it only for WMS/WFS services? If it can be used to create a store metadata for a PostGIS table, how is it done? I can't seem to find this in the GeoNetwork documentation. Thanks. DBlack


title: Write to FGDB in C++ code?

text: Anyone know the correct syntax for this? I'm a C++ beginner but experienced in GIS.. I'm guessing ofstream won't do it and am not sure what else to try. Never written to any sort of database in C++.


title: Converting/clipping KML to shapefile in QGIS without truncating field information (data)?

text: While I am clipping KML file by using polygon shapefile in QGIS, filed information (data) are truncating .Anybody can suggest a solution for this one.The same data losing issue occurring, at the time of KML to shape file conversion also


title: QGIS, CRS's and metric measurements of distance matrix and nearest neighbour

text: I have thousands of points on a multiple layers and I need to get the distances between each point and every other point within a layer (so i can get min, max, mean and median distances) and the nearest neighbour values. I have used the distance matrix to achieve distances between each point, and the nearest neighbour tool to get the nearest neighbour value. It all seemed to be working fine until I realised that the distance values are wrong. I assumed they were in km as the measuring tool displayed the distances of known distances correctly. However, when running it through the distance matrix and nearest neighbour, the values are massively different. In a similar vein to this question: Distance Matrix tool in QGIS. How to get result in metric?   * **Example:**   * **Point 1 to Point 2** is 1.12km in reality   * **Manual measuring tool on QGIS** says 1.12km   * **Distance matrix** says 0.010682   * **Nearest neighbour mean observed distance** says 0.010682 I have tried changing CRS and reprojecting and saving as a new shapefile, but it still doesn't work. I can get the measuring tool to display the correct measurements but as soon as i run it through the distance matrix or the nearest neighbour i get the same wrong values no matter what the CRS is. I thought it might be in degrees, but how do I then convert it to km's so when I am writing up my results I can say that 'the mean value was Xkm' rather than some obscure number which doesn't readily mean anything to anyone? Can anyone offer any advice? Is there a projection I should be using (my points are plotted from long/lat)? Is there a setting I need to change for the plugins? Is there another tool which will do this for larger numbers of points? Is there a formula I can use to change all the 'wrong' values (possibly degrees) into m/km's? I don't know python, but is there a clever bit of code I can use to do this? I'm lost and this is driving me insane! Help gratefully appreciated! **UPDATE** This worked for a while, but now it's reverted again. I managed to batch execute the reprojection, and then batch distance matrix and nearest neighbour and it seemed to work. For all except one layer. I've tried reprojecting just one file using save as manually, but it won't work. I've installed it on a second computer and tried it won't work. I went back to the original data, imported it into QGIS from csv, straight into ESPG:27700 with the project set to ESPG:27700 and it's still measuring in degrees. It's not the coord's, i've checked and know what the measurements should be. I then ran it with the second year of data ready to disregard one bad file, and it is refusing to work. I really don't understand the problem. It says it's reprojected the file, but it isn't, or it is but something is making QGIS continue to measure in degrees. I didn't changed any settings from the set which worked, to the set which didn't. Is this a glitch in the program? Is there another free GIS program that can do distance matrix and nearest neighbour which might save me this continuous circling? **UPDATE2** I have just taken an original unsplit file and reprojected the entire year as one file. I saved and then split it by date, before using distance matrix. It is still running and outputting files and this is likely to take some time but the first couple of outputs so far appear to be in metres. I believe the problems I am experiencing are related to a bug. I will try to report it as such, but for anyone experiencing similar problems all I can suggest is that you disregard any files which aren't working and start right from the very beginning with the data, reprojecting before any other analysis as it seems something lodges itself in the file and stops it working properly when you try to run distance related analysis at some point.


title: Reduce shapefile to outline

text: I have a shapefile consisting of 10 contiguous polygons. I want to reduce this to a single polygon encompassing the entire area of the existing 10. I don't need to retain data. I've googled but without the proper terminology I get pointed to Merge and Clip functions which I don't think are what I am looking for.


title: Pull Data from Specific Excel Cell into ArcPy

text: I am new to Python and ArcPy, and wish to create a short script that pulls in a value from one specific Excel worksheet cell (namely, A2) and then performs an 'if' statement, whereby if the value of A2 is greater than 1, then a certain text string is returned, but if it is less than or equal to 1, then a different text string is returned. So far, my code looks like this:               import arcpy     from arcpy import env     arcpy.env.workspace = r"C:\MMO_Model_Test\testdata.xlsx"     inputData = "data"      So I have defined my workspace, and set the relevant worksheet as a variable. I now need to reference the cell 'A2' as part of an 'if' statement. This is probably simple for someone who knows how, but I haven't managed to find a relevant example online. Any help is hugely appreciated.


title: Delete files that were used in python code

text: I have a code that creates some tif images in .qgis2 folder, but I want, in the end, delete these files. How do I do that? I try with os.remove but it gives me the error32, that said that the process can´t access to file because this is used by another process. Hope you help me.


title: selection by attributes with python

text: I have a problem concerning the selection by attribute with python on arcmap. I use this code:               import arcpy     from arcpy import env     arcpy.env.workspace= "E:\cities"     Mi=(r"E:\cities\rivers.shp")     rows = arcpy.SearchCursor(Mi,'"NAME" = \'Mississippi\'')     row = rows.next()     largeur = row.getValue("MILES")     arcpy.selectLayerByAttribute_management("rivers","New_SELECTION",'"MILES" >' + str(largeur))      And i got this error: > Traceback (most recent call last): File "E:\cities\scripts\Selection.py", > line 8, in > arcpy.selectLayerByAttribute_management("rivers","New_SELECTION",'"MILES" > >> ' + str(largeur)) AttributeError: 'module' object has no attribute 'selectLayerByAttribute_management' > > Échec de l’exécution de (SelectionByAttributes) Can u help me please??


title: Failure to Access DBMS Server Error - ArcGIS/Postgres

text: If connected to an SDE enabled Postgres database via ArcGIS (10.1) for a long period of time, maybe an hour +, I will often receive an error from ArcGIS stating "Failure to access the DBMS Server[my schema name.sde.GDB_items]. I can't pin down any specific event causing the error, and my fix is to simply restart ArcGIS. Disconnecting and re-connecting the database connection without an ArcGIS restart results in the same error message. Its a minor inconvenience, but one that would save myself and others daily grief if corrected. Any ideas?


title: How to get the width and height of a rectangular polygon in PostGIS

text: I have a table with point features saved in the my PostgreSQL database, and I can run the following query to get the extent of the point features and I turned the extent back into a geometry:               SELECT ST_AsText(ST_SetSRID(ST_Extent(pts.geom), 4026)) AS Extent_Geom           FROM Table_w_Point_Features pts      Here is how the output geometry of a rectangle polygon that represents the extent of the point features looks like:               POLYGON((-117.536012078 33.7130048820001,              -117.536012078 34.3860049530001,              -116.825011941 34.3860049530001,              -116.825011941 33.7130048820001,              -117.536012078 33.7130048820001))      And here is my question: Is there any way in PostGIS that I can get the Width and Height of the extent that I created above? Thanks a lot!


title: Garmin conversion track

text: I have problem when import garmin GPS tracks with software for Android. I found this code:               public static int importTracks(String path, Context context) throws ClassNotFoundException, IOException     {         TracksDB tracksDB = new TracksDB(context);         TrackPointsDB trackPointsDB = new TrackPointsDB(context);         tracksDB.open();         trackPointsDB.open();              File file = new File(path);         ObjectInputStream in = new ObjectInputStream(new FileInputStream(file));              int count = 0;              try         {             Track track;             while ((track = (Track) in.readObject()) != null)             {                 long id = tracksDB.insert(track.getName());                      ArrayList<TrackPoint> points = track.getTrackPoints();                 for (TrackPoint point : points)                 {                     trackPointsDB.insert(id, point.getLongitude(), point.getLatitude(), point.getAltitude(), point.getAccuracy(),                             Long.toString(point.getTime()));                 }                      count++;             }         }         catch (EOFException eof)         {              }              // !// finally?         in.close();         tracksDB.close();         trackPointsDB.close();              return count;     }      Can someone help me?


title: Problems converting csv to shp in QGIS 2.01 on a Mac

text: When converting a csv to shp, I have been using the 'Add delimited text layer', and get one of two errors. With 'First record has field names' clicked, I am told that there is 'no data found in file'.   With 'First record has field names' unclicked, I am able to save the file, but get errors that there are 'invalid X or Y fields at line 1.' If I drag in the csv, and do 'Layer' 'Save as''shapefile', it will create a .dbf, .prj, and .qpj, but no .shp. Any help?


title: Cannot Save Edits to File Geodatabase Replica in QGIS 2.2

text: When trying to edit a File Geodatabase Replica in QGIS 2.2, everything works as expected until the user saves their edits. When pressing save, all edits revert back to their original state with no errors. If I copy the Feature Class into a new FGDB (thus breaking it out of the replica), I can edit it without any problems. Anyone have any thoughts on a solution for this? FGDB Replica created in ArcGIS Desktop 10.2 from a SQL SDE Geodatabase.


title: Update multiple features attributes in loop - JavaScript

text: I am trying to update multiple features in a feature layer through a loop. I have a feature layer that I would like to select and then assign someone to the features (through an attribute update. I can select the features and see all of the attributes but I can't figure out how to update the the features. See below for a sample of what I am trying to do. Right now it is a bit messy where I am trying to set the visibility of a form and button on after the selection so the user can choose an assignment. Has anyone updated multiple features without using the attribute inspector or editor widget which I have only seen examples for editing a feature at a time.                var map, geocoder, currentBaseMap, assignWindow, featureLayer;      var locatorUrl ="http://gisaprd/ArcGIS/rest/services/SDE.AddressLocator_Comp/GeocodeServer";      var selectedPoints = new Array();          require([     "esri/InfoTemplate",     "esri/map",     "esri/layers/FeatureLayer",     "esri/layers/ArcGISTiledMapServiceLayer",     "esri/layers/ArcGISDynamicMapServiceLayer",     "esri/dijit/Geocoder",     "esri/symbols/SimpleFillSymbol",     "esri/symbols/SimpleLineSymbol",     "esri/tasks/query",     "esri/toolbars/draw",     "dojo/dom",     "dojo/on",     "dojo/parser",     "dojo/_base/array",     "dojo/_base/Color",     "dijit/form/Button",     "dojo/domReady!"      ],      function (      InfoTemplate, Map, FeatureLayer, ArcGISTiledMapServiceLayer, ArcGISDynamicMapServiceLayer,         Geocoder, SimpleFillSymbol, SimpleLineSymbol,     Query, Draw, dom, on, parser, arrayUtil, Color      ) {            parser.parse();            var selectionToolbar, featureLayer;            // Create the map       map = new Map("map")            aerialsLayer = new esri.layers.ArcGISTiledMapServiceLayer("http://gisaprd/ArcGIS/rest/services/Aerials_4inch_12inch_Cached/MapServer", { "id": "Aerials", "visible": false });       map.addLayer(aerialsLayer);            streetsLayer = new esri.layers.ArcGISDynamicMapServiceLayer("http://gisaprd/ArcGIS/rest/services/SRSLandbase/MapServer", { "id": "Streets" });       map.addLayer(streetsLayer);       currentBaseMap = "Streets";            map.on("load", initSelectToolbar);            document.getElementById("AssignFrm").style.visibility = "hidden";       document.getElementById("AssignBTN").style.visibility = "hidden";            highlightSymbol = new esri.symbol.SimpleMarkerSymbol().setColor(new dojo.Color([255, 255, 0]));            var content = "Address: ${SDE.GasServiceLocation.LOCATIONDESCRIPTION}";       var infoTemplate = new InfoTemplate("${SDE.GasServiceLocation.CU_ID}", content);            featureLayer = new FeatureLayer("http://gisadev:6080/arcgis/rest/services/ISS_Inspections/FeatureServer/0",         {             mode: FeatureLayer.MODE_ONDEMAND,             infoTemplate: infoTemplate,             outFields: ["*"]         });            featureLayer.setSelectionSymbol(highlightSymbol);       featureLayer.on("selection-complete", createListForAssignment);       //featureLayer.on("selection-clear", function () {       //    dom.byId('messages').innerHTML = "No Selected Fields";       //});       map.addLayer(featureLayer);            var myGeocoders = [{           url: locatorUrl,           name: "SDE.AddressLocator_Comp"       }];       geocoder = new Geocoder({           map: map,           autoComplete: true,           arcgisGeocoder: false,           geocoders: myGeocoders,           value: "",           zoomScale: 1200            }, "search");       geocoder.startup();            on(dom.byId("selectFieldsButton"), "click", function () {           selectionToolbar.activate(Draw.EXTENT);       });       on(dom.byId("clearSelectionButton"), "click", function () {           featureLayer.clearSelection();       });       on(dom.byId("BaseMapToggle"), "click", ChangeBaseMap);       //on(dom.byId("assignCrew"), "click", AssignInspectionCrew);            function initSelectToolbar(event) {           selectionToolbar = new Draw(event.map);           var selectQuery = new Query();                on(selectionToolbar, "DrawEnd", function (geometry) {               selectionToolbar.deactivate();               selectQuery.geometry = geometry;               featureLayer.selectFeatures(selectQuery,                 FeatureLayer.SELECTION_NEW);                    });       }       function ChangeBaseMap() {           if (currentBaseMap == "Streets") {               currentBaseMap = "Aerials";                    aerialsLayer.setVisibility(true);           }           else if (currentBaseMap == "Aerials") {               currentBaseMap = "Streets";               aerialsLayer.setVisibility(false);           }       }       function createListForAssignment(event) {           document.getElementById("AssignFrm").style.visibility = "visible";           document.getElementById("AssignBTN").style.visibility = "visible";                arrayUtil.forEach(event.features, function (feature, index) {               selectedPoints.push(feature);           });        }       });       function assign() {       if (document.getElementById('Jordan').checked) {         x = document.getElementById('Jordan').value;         for (var i = 0; i < selectedPoints.length; i++) {             selectedPoints[i].attributes.ASSIGNED_TO = x;         }       }       else {         y = document.getElementById('Trenton').value;         alert(y);       }       }


title: Extracting z values from shapefile using OGR Api

text: Is there an example of how we can extract the z value from a shp file given the x,y coordinate value (namely utm) using gdal/ogr c++ lib? alternatively, would it be easier if we convert the shp file into raster layer format and read off value from there? Woukd it be easier? Is there some examples of how this could be done? Thanks


title: Map accuracy assessment by moving window in R

text: I would like to make an **accuracy assessment** using a confusion matrix between classified Landsat image and reference dataset. Maps have the same resolution and extent. I would like to evaluate the agreement by **pixel-by- pixel**. In many studies I have found that they used **moving-kernel window** (ideal 3x3) to "deal" with the Landsat pixels misregistration. However I couldnt find any approach to do use this moving window for confusion matrix evaluation in R, normally it is used to values interpolation. Do you have any ideas how to **implement moving window into classification accuracy assessment**? Or am I misunderstanding this approach? Thanks a lot, Example:               library("raster")     clas <- raster(ncols=5, nrows=5)     values(clas)<-c(2,1,0,0,1,              2,1,0,0,0,              1,0,0,0,0,              0,0,0,0,0,              2,1,0,0,2)     reference <- raster(ncols=5, nrows=5)     values(reference)<-c(2,1,0,0,0,              2,1,0,0,1,              0,0,0,0,0,              1,2,0,0,2,              0,0,0,0,0)      EDIT: After utilisation of movingFun (raster) I think this is mostly to get information in one cell taking into consideration the cells around. so if I have a               ref<-c(0,0,0,0,1,0,0,0,0)     class<-c(0,0,0,0,1,1,1,0,1)      and I apply movingFun               mov_ref<-movingFun(ref, n=3)     mov_clas<-movingFun(class, n=3)      I will obtain               mov_ref     NA, 0.0000000 0.0000000 0.3333333 0.3333333 0.3333333 0.0000000     0.0000000        NA          mov_clas     NA 0.0000000 0.0000000 0.3333333 0.6666667 1.0000000 0.6666667     0.6666667        NA      For **accuracy assessment needs** could I directly compare mov_ref and mov_class? Or do I have maybe to reclass it into 0,1 (by threshold 0.5) to compare them? If I leave the numbers as 0.3333 I dont reach a good result in classification accuracy. Thanks for your help.


title: OpenLayers vector feature not being highlighted

text: I'm trying to highlight a clicked area of a layer using a vector:               var highlightLayer = new OpenLayers.Layer.Vector("Highlighted Features", {         displayInLayerSwitcher: false,         styleMap: mystyle     });     // the layer is added to the map, GetFeatureInfo controls setup and added etc.      ...and this is the callback for the GetFeatureInfo call:               function showInfo(event) {         console.log( event );         highlightLayer.addFeatures( event.features[0] );         highlightLayer.redraw(); // i've tried without this too         console.log( highlightLayer );     }      No highlighting happens though. The object (event) passed to the callback looks alright, it has a features property containing an array with one vector feature element. But while inspecting the highlightLayer vector object printed to the console in the callback, I noticed that the returned vector feature object has been added to the `unrenderedfeatures` property. OpenLayers obviously could not render the feature, but I don't understand the reason. There doesn't seem to be a callback for detecting problems with rendering. How could I debug this further? I should also mention that the `bounds` property of the geometry object is a valid `OpenLayers.Bounds` object with values. **Edit** :   I dug into the OpenLayers source and added some debug logging. It turns out that the redrawFeature function belonging to the renderer is returning false, resulting in the feature to be added to the unrenderedfeatures property of the vector layer. The reason it's returning false is because of this:               var bounds = feature.geometry.getBounds();     var worldBounds;     if (this.map.baseLayer && this.map.baseLayer.wrapDateLine) {         worldBounds = this.map.getMaxExtent();     }     var intersects = bounds && bounds.intersectsBounds(this.extent, {worldBounds: worldBounds});     rendered = (style.display !== "none") && !!bounds && intersects;      The `intersects` variable is **false**. Why might that be? Could it be that some sort of conversion is needed on the bounds values returned from the GetFeatureInfo call? **Edit2:**   Here's some logging of the points that are compared in the `intersectsBounds` function: Feature / Renderer-Canvas-Extent   Bottom: 719322.39 / 7022915.0167624   Left: 689914.61 / -740991.87193189   Right: 689947.66 / -740514.14050517   Top: 719349.21 / 7023273.3153324


title: How to calculate Median of a table field in the Dissolve Tool of ArcGIS?

text: The Dissolve Tool of ArcGIS 10.0 has some statistical parameters as a type of calculation new values of each field in a process of dissolving of identical parameters . It has SUM, MAX, MIN, RANGE, STD and a bit more Statistical types , but it doesn’t have MEDIAN. It’s very strange, why it doesn’t have Median (and some others important statistical parameters). But maybe it’s possible to “turn on” calculation “Median” in the Dissolve Tool, or solve this task in some other way?


title: Is it possible to use two layer.bindPopup functions in Leaflet

text:
I have two different layers in a leaflet map. I have one layer in which I would like the popup to return some text from the layer's attributes. They other layer I would like the popup to return an image which I have embedded the link to in the layers attributes. I have successfully implemented both of the layer.bindpopup functions successfully individually but I'm having no luck implementing them together.               function onPlot(feature, layer) {     if (feature.properties) {         layer.bindPopup('');     }     };     function onEachFeature(feature, layer) {     if (feature.properties) {         layer.bindPopup( '' +feature.properties.Anum +'
 Hunting District: ' +feature.properties.DISTRICT+ "
 Elk Hunter Success " +feature.properties.Success +" %");     }     };     function style (feature,layer){     if (feature.properties) {      return  {color: feature.properties.color,fillColor: feature.properties.fillColor,fillOpacity: feature.properties.fillOpacity};     }     };          var Plot= new L.geoJson(Plots, {         onPlot:onPlot         style:style     }).addTo(map);          var S_2006= new L.geoJson(Success_2006, {         onEachFeature: onEachFeature,         style:style     }).addTo(map);


title: qgis pdf doesn't show google street

text: Apologies I've created a map? using location data imported from a csv file and labelled nicely and I've added (a very inaccurate) Google Street map. When I use Print Composer to create an annotated map with title, legend etc. and then export to pdf - everything is exported except the Google Street background map. I must be doing something wrong, can some kind person lead me on the path to righteousness please?


title: jQuery element value to js object ID

text: I am trying to set the visibility of a layer in Openlayers map according to a checkbox status with a class ".layercheck" and value is the name of the Openlayers layer. So here is my HTML:               <input type="checkbox" class="layercheck" value="layer134431245">text      And here is JS:                  var layer134431245 = new OpenLayers.Layer.WMS(...);         $(".layercheck").on( "click", function() {         var layer = $(this).val();         if ($(this).is(":checked")) {             $("#"+layer).setVisibility(true);         } else {             $("#"+layer).setVisibility(false);              }        });      This function and HTML does completely nothing. It also does not give an error in the js console. So how do I make this work, and switch the layer visibility with any of the checkboxes which have the class "layercheck" and value is the layer ID? I also understand that the $("#"+layer) is not the layer, but a jQuery object which maybe does not have the properties setVisibility. So then how do I convert a value (string) into a Openlayers layer name?


title: transform bound for google map

text: i have wms layer with coming geoserver and it's bounds is                bounds = new OpenLayers.Bounds(417138.6205791628, 4243634.852681775, 486211.1950145211, 4331542.555877601).transform(WGS84, WGS84_google_mercator);      but when i use google baselayer, in the opening map it shows me full extends how can i transform my bound with harmonic google. ?


title: Error when using spatial query in arcgis silverlight api

text: I'm developing a webgis using ArcGIS Silverlight API 3.1 . I used esri interactive SDK sample for spatial query in the webgis .The problem is when i draw a freehand polygon and the polygon intersect itself (fig1), the error "Query failed: Error code '400':unable to complete operation." raised (fig2) . How can i solve the problem ? ![fig1](http://i.stack.imgur.com/toGaa.png) ![fig2](http://i.stack.imgur.com/ztqJV.png)


title: SST process using just GDAL and python

text: I am trying to process AVHRR 4km SST using a module that uses only GDAL and Python. OSGEO4W, GDAL, Python 3.4.1. I borrowed the download code from http://grasswiki.osgeo.org/wiki/AVHRR. The rest I worked on. The code is highly inefficient but it works. My concerns are: 1) Being able to adjust quality flags. I have it set for >3 but it is easy to change. 2) Filling small Nodata holes. 3) Keeping the integrity of the land mask. **I need help automating some of this such as reading the scale automatically from the .txt. Running the multiply of the SST and the NoData,1 Binary quality data in an automated manner. This code works but is piecemeal and inefficient. I would like to make it efficient** I would like to streamline and improve the process but we must stay in GDAL python. V5.0 Eight day. Make our own month with high quality flags of 4 or grater only and fill small hole via a 5 by 5 kernel. http://www.nodc.noaa.gov/SatelliteData/pathfinder4km/userguide.html 1) Get The HDFs Via Filezilla or a web browser are two options. FTP http://data.nodc.noaa.gov/pathfinder/Version5.0/8day/2000/ 2) Downlaods and create folders, be sure to note 32 bit or 64 bit. 2a) CYgwin with Wget (https://www.cygwin.com/) 2b) OTSGEO4w From VTECH or similar (http://trac.osgeo.org/osgeo4w/) 2c) Install Python (https://www.python.org/downloads/release/python-341/) 2d) Bind it to GDAL http://www.lfd.uci.edu/~gohlke/pythonlibs/%20(GDAL%E2%80%911.11.0.win%E2%80%91amd64%E2%80%91py3.4.exe)%20or%20(GDAL%E2%80%911.11.0.win32%E2%80%91py3.4.exe) 2e) Create 4 Folder in C:\temp QUAL TXT ARCHIVE LANDMASK and download the landmask http://www.nodc.noaa.gov/SatelliteData/pathfinder4km/userguide.html 3) Get Data using Cygwin 8-day for January 2000 ---change drive to C:\temp Point out how to change night or day . Download just night (clouds) for TIME in 1 = night for TIME in 3 = day for TIME in 1 3 = night and day               DAYS="\     2000001-2000008     2000009-2000016     2000017-2000024     2000025-2000032"          BASE_URL="http://data.nodc.noaa.gov/pathfinder/Version5.0/8day/2000/"     for TIME in 1 ; do       for JULIAN in $DAYS ; do         wget "${BASE_URL}${JULIAN}.s048${TIME}pfv50-sst-16b.hdf"         wget "${BASE_URL}${JULIAN}.m048${TIME}pfv50-qual.hdf"       done     done      4) Convert to TIFF - HDF to Tiff, Set the projection, set the extent, compress the data, convert to GTiff. then move old HDF manually to Archive and move the quality files to QUAL. -ot Float32 to use decimals but we are using Ints               for %A in ("C:\temp\*.hdf") do gdal_translate -of GTiff -a_srs "+init=epsg:4326" -a_ullr -180 90 180 -90 -co "COMPRESS=PACKBITS" -sds "%A" "%A.tiff"      5) Set NoData, a_nodata 0 will convert the 0 to NoData if required (not required for this method but usually is) 6) Sacle Factor - Check for Metadata Scale factor (hopefully 0.07500000298) and Move Quality ito own folder (QUAL) and do same for Text (TXT)               for %A in ("C:\temp\*.tiff") do gdalinfo "%A" >>"%A.txt"      7) Quality 1 - In quality field- Convert Quality <3 to 0 and Quality >3 to 1. Manually move old quality files to Archive.               for %A in ("C:\temp\QUAL\*.tiff") do gdal_calc.py -A "%A" --outfile="%ABINARY.tiff" --calc="0*(A<3)" --calc="1*(A>3)"      8) Scale data from 6, be sure to first to cd back to temp cd .. Move old files without_C to Archive NOTE THis is for whole degrees if decimals are required set output raster to float. in step 4.               for %A in ("C:\temp\*.tiff") do gdal_calc.py -A "%A" --outfile="%A_C.tiff" --calc="A*(0.07500000298)"      8) Transfer Quality to SST files, so if quality is 4,5,6,7 then times SST pixels by 1 (QUAL) so they stay the same. If Quality if 0,1,2,3 times by zero (QUAL) so they become zero.               gdal_calc.py -A 2000001-2000008.s0481pfv50-sst-16b.hdf.tiff_C.tiff -B C:\temp\QUAL\2000001-2000008.m0481pfv50-qual.hdf.tiffBINARY.tiff  --outfile=A1_8_SSTQUAL4.tiff --calc="A*B"     gdal_calc.py -A 2000009-2000016.s0481pfv50-sst-16b.hdf.tiff_C.tiff -B C:\temp\QUAL\2000009-2000016.m0481pfv50-qual.hdf.tiffBINARY.tiff  --outfile=A9_16_SSTQUAL4.tiff --calc="A*B"     gdal_calc.py -A 2000017-2000024.s0481pfv50-sst-16b.hdf.tiff_C.tiff -B C:\temp\QUAL\2000017-2000024.m0481pfv50-qual.hdf.tiffBINARY.tiff  --outfile=A17_24_SSTQUAL4.tiff --calc="A*B"     gdal_calc.py -A 2000025-2000032.s0481pfv50-sst-16b.hdf.tiff_C.tiff -B C:\temp\QUAL\2000025-2000032.m0481pfv50-qual.hdf.tiffBINARY.tiff  --outfile=A25_32_SSTQUAL4.tiff --calc="A*B"      9) Combine SST - Combine all the SST QUALT 4 rasters together               gdal_calc.py -A A1_8_SSTQUAL4.tiff -B A9_16_SSTQUAL4.tiff -C A17_24_SSTQUAL4.tiff -D A25_32_SSTQUAL4.tiff --outfile=SST_SUM_C.tiff --calc="A+B+C+D"      10) Combine Quality - Combine all Quality Binaries together (0-4) 0 = no good pixels and 4 = Four Good Pixels. Be sure to cd into QUAL Folder               gdal_calc.py -A 2000001-2000008.m0481pfv50-qual.hdf.tiffBINARY.tiff -B 2000009-2000016.m0481pfv50-qual.hdf.tiffBINARY.tiff -C 2000017-2000024.m0481pfv50-qual.hdf.tiffBINARY.tiff -D 2000025-2000032.m0481pfv50-qual.hdf.tiffBINARY.tiff  --outfile=QUAL_final.tiff --calc="A+B+C+D"      11) Get to Monthly Mean - Divide SST by Pixel Count (Quality greate than four 0,1,2,3,4) Be sure to cd into temp Folder               gdal_calc.py -A SST_SUM_C.tiff -B C:\temp\QUAL\QUAL_final.tiff  --outfile=SST_Final_C.tiff --calc="A/B"      11) Convert 0 to NoData, calculate stats, Float . Move final to archive               gdal_translate SST_Final_C.tiff -of GTiff -ot Float32 -a_nodata 0 SST_Final_C_NoData.tiff      12) ALthough this has Stats for QGIS ir GRass ArcGIS requires you to build Sstats in ArcCatalog (Clear first then make new in arcatalog). 13a) In ArcMap do this to fill small holes. Watch out this will shrink the land. focal_sst_C (do it to a geodatabase) NOT USING               Con(IsNull("SST_Final_C_NoData.tiff"), FocalStatistics("SST_Final_C_NoData.tiff", NbrRectangle(5,5, "CELL"), "MEAN"), "SST_Final_C_NoData.tiff")      13b) Use NoData filler in QGIS 13c) gdal_fillnodata SST_Final_C_NoData.tiff -md 5 -si 2 SST_Final_C_NoData_Foc.tiff 14) CHange to Mask to tiff, as step 4 but this time add change NoDta to 2 (delete old HDF file) cd to landmask               for %A in ("C:\temp\landmask\*.hdf") do gdal_translate -of GTiff -a_srs "+init=epsg:4326" -a_nodata 2 -a_ullr -180 90 180 -90 -co "COMPRESS=PACKBITS" -sds "%A" "%A.tiff"      15) Change 2 to 1 in landmask `for %A in ("C:\temp\landmask\*.tiff") do gdal_calc.py -A "%A" --outfile="%A_FINAL.tiff" --calc="1*(A<`5)" NOW MULTIPLY BY FINAL 16) Burn landmask back into final SST, cd ..               gdal_calc.py -A SST_Final_C_NoData_Foc.tiff -B C:\temp\landmask\pfv50_land.m04.hdf.tiff_FINAL.tiff  --outfile=SST.tiff --calc="A*B"      17) General Clean up Build sats, compress, 0 to No Data,               gdal_translate SST.tiff -of GTiff -ot UInt16 -a_nodata 0 -stats -co "COMPRESS=LZW" SST_ALLDONE.tiff


title: Can we manually select which satellites to use in solution?

text: For example, I have used below satellites ID in current solution:               ID    Signal Strength     ______________________      1         21     22         23     19         12     36         26     20         18      3         29     32         14     13         15      But I want to filter and use only those satellites with **higher signal strength** in my solution, let say `satellite ID: 1, 22, 36, 3`. Is it possible?


title: OpenLayers with MVC and SQL Server

text: I have imported some shape files into SQL Server 2008. Now i want to use OpenLayers with my MVC 4 application to render those shape files from SQL Server Database. Any help in this regard?


I have successfully use tilestache-seed.py to get tiles in different zoom levels, and I use processing to render the geojson tiles. Now I want to use tilestache-server.py to get local tiles in a way like local/roads/x/y.geojson, maybe something like this. But I don't know how to write the cfg files.


title: Draw point on google map

text: i try to add google map on my own wms and wfs layer. everything works but when i create point, i intersect point geometry and my geojson data. So without googlemap it was run. But when add google map, points cordinates will be different and intersect result is wrong.you can see picture what i said. when i click on africa, intersect results is correct but i must click on my overlays layers. http://i.hizliresim.com/xVyDJj.png WFS_LAYER                wfs_layer = new OpenLayers.Layer.Vector("Point Layer",             {                      styleMap: new OpenLayers.StyleMap({                     "default": new OpenLayers.Style({                         pointRadius: 10, // sized according to type attribute                         fillColor: "#00FFFF",                         strokeColor: "#33CCCC",                         strokeWidth: 3,                         graphicZIndex: 1                     }),                 }),                      strategies: [new OpenLayers.Strategy.BBOX(), saveStrategy],                      protocol: new OpenLayers.Protocol.WFS({                     version: version,//1.1.0                     // loading data through localhost url path                     url: wfsUrl,                     featureNS: featureNS,                     visibility: true,                     // layer name                     featureType: featureType,                     // geometry column name                     srsName: "EPSG:900913",                     geometryName: geometryName,                     schema: schema,                      })             });      point intersect :                function mahalleBul(request) {                      var geojson_format = new OpenLayers.Format.GeoJSON({                     internalProjection: new OpenLayers.Projection("EPSG:900913"),                     externalProjection: new OpenLayers.Projection("EPSG:4326")                 });                 var getMahalle = geojson_format.read(request.responseText);                 if (getMahalle.length > 0) {                     mahalleadi = getMahalle[0].attributes.ADI;                     mahallegisID = getMahalle[0].attributes.GIS_ID;                 }     }      var request = OpenLayers.Request.GET({                 url: rootURL + "ows?service=WFS&version=1.0.0&request=GetFeature&typeName=Kaski:ADRES_MAHALLE&outputFormat=json&cql_filter=INTERSECTS(the_geom," + feature.geometry + ")",                           callback: mahalleBul             });


title: How to convert WGS84 to UTM when the study area covers two zones: UTM 51N and 52N

text: I want to use ArcMap to convert the WGS84 to UTM. But the study area covers two zones: UTM 51N and 52N. Which zone should I choose to output the result?


title: Is there maximum number of .DBFs that can be joined using ModelBuilder?

text: Is there a maximum number of .DBFs that can be joined using the Join Tool in ModelBuilder? I tried joining 480 .DBFs but results generated only 250 fields in the input table. Is this a glitch? I am using ArcGIS for Desktop 10.2.


title: OGC Compliance Test for GeoServer WMS 1.1

text: I have installed OpenGeo Suite and using GeoServer installed under this suite. Now i want to get OGC Compliance Certificate for WMS.   For the same.. I have also import the cite data in GeoServer. Now when i am running the WMS1.1.1 test from OGC portal Team Engine page it is failed for following test results:   1. Test wms:wmsops-getmap-params-bgcolor-1   2. Test wms:wmsops-getmap-params-bgcolor-2   3. Test wms:wmsops-getmap-params-transparent-1   4. Test wms:wmsops-getmap-params-transparent-2   5. Test wms:wmsops-getmap-each-style   6. Test wms:wmsops-getmap-each-style The log of test 1 is:               Test wms:wmsops-getmap-params-bgcolor-1 type Mandatory default result Passed (s0028/d1e1515_1)          Assertion: When a GetMap request is made with no BGCOLOR parameter, then areas which are not covered contain white (0xFFFFFF) pixels.          Request d1e7670_1:        Method: get        URL: `http://182.71.192.178:8090/geoserver/wms?SERVICE=WMS&&HeIgHt=100&SrS=EPSG:4326&BbOx=0,-0.0020,0.0040,0&VeRsIoN=1.1.1&ReQuEsT=GetMap&LaYeRs=cite:Lakes&WiDtH=200&FoRmAt=image/png&StYlEs=`        Response from parser example:ImageStats:           <image_stats>                                                             <model value="RGB">                                                                     <subimage height="100" width="25" x="0" y="0">                                                                             <count bands="RGB" sample="0xffffff">25</count>                                                                     </subimage>                                                                     <subimage height="50" width="50" x="150" y="0">                                                                             <count bands="RGB" sample="0xffffff">50</count>                                                                     </subimage>                                                             </model>                                                     <!--Response received in [477] milliseconds--></image_stats>          Result: Failed      But when i am browsing the same url:               http://182.71.192.178:8090/geoserver/wms?SERVICE=WMS&&HeIgHt=100&SrS=EPSG:4326&BbOx=0,-0.0020,0.0040,0&VeRsIoN=1.1.1&ReQuEsT=GetMap&LaYeRs=cite:Lakes&WiDtH=200&FoRmAt=image/png&StYlEs=      it is giving the showing the same image according to test condition. Same is for tests 2, 3 and 4. Please help me to resolve the problems. Thanks & Regards, Manish Gupta


title: What is the distinction between Difference and Clip in QGIS?

text: What is the primary distinction between the **Difference** and **Clip** geoprocessing tools in QGIS?


title: Make Google Analytics city data showed on the Google map

text: I use the Google fusion table to table to mark city on google map, but there are about 20~25% which quality is ambiguous. How to make the ambiguous data precise?   Here is my finished map (scroll down to see map)


title: How can I use data from PDF map?

text: I have a map over Sweden with plenty of data represented as circles of various sizes. The map is saved as PDF from Illustrator, fully editable. I wonder if there is any way I could get these circles into Tilemill so that I can make a searchable georeferenced map? I don't know the projection but I would guess WGS84 or SWEREF99. One way that comes to mind is to extrakt the geodata from a geotiff of Sweden (if I find one with the same projection), place PDF-map on top of the geotiff in Illustrator/Photoshop and then use gdal to insert the geodata again. Any ideas if this is a good way, or if there are better methods?![enter image description here](http://i.stack.imgur.com/sL2Sf.png)


title: looking for a GPS device which can transmit the location directly

text: I am looking for a `GPS` device or receiver that can transmit its own location to a specified configured server directly. there will be no medium like mobile phone or computer in between GPS and Server which receives the data. a dongle sort of thing which has a sim card in it and with the use of that it can transfer data directly to server


title: Capturing/overriding mousewheel zoom in openlayers

text: So I've created a map in openlayers and I am setting it up so that there are only certain zoom levels a user can see. I have the code created to handle forcing a user to the next acceptable zoom level when they attempt to zoom in/out. I have been able to create listener for handling double-clicks (see code below) but for the life of me I cannot figure out how to capture mousewheel events. The documentation makes it seem like it should be straightforward but I've been banging my head against a wall for more than a day trying to figure this out. All my googling has turned up plenty of results for disabling the mousewheel but nothing for capturing wheel events.               function init() {         map = new OpenLayers.Map('map_canvas',mapSettings);         ....         ....         );              var click = new OpenLayers.Control.Click();         map.addControl(click);         click.activate();     }          OpenLayers.Control.Click = OpenLayers.Class(OpenLayers.Control, {         defaultHandlerOptions: {         'single': false,         'double': true,         'pixelTolerance': 0,         'stopSingle': false,         'stopDouble': false         },              initialize: function(options) {             this.handlerOptions = OpenLayers.Util.extend(             {}, this.defaultHandlerOptions             );             OpenLayers.Control.prototype.initialize.apply(             this, arguments             );             this.handler = new OpenLayers.Handler.Click(             this, {             'dblclick': this.trigger             }, this.handlerOptions             );         },              trigger: function(e) {             var lonlat = map.getLonLatFromViewPortPx(e.xy);             curZoom = map.getZoom();             var newZoom = stepZoomIn(map.getZoom());             updateMapZoom(newZoom);             map.setCenter(event.lonlat);         }          });


title: How to filter out multipart polygons and then buffer?

text: I've written a small script trying to accomplish the task described in the question title. It basically looks like:               in_fc = '...' # Some shapefile goes here          # Get the fields (column names) of the attribute table     desc = arcpy.Describe(in_fc)     shape_field = desc.ShapeFieldName          # Initialize an empty polygon list     singlepolys = []          # Open the attribute table of the shapefile     rows = arcpy.SearchCursor(in_fc)     for row in rows:         poly = row.getValue(shape_field)              # Check if single part (not multipart)         if not poly.isMultipart:             singlepolys.append(poly)          # Print areas of single polygons - just for debugging     for p in singlepolys:         print p.area          # Buffer code would go here...      The final printing of the areas shows that only multiples of the last polygon have been copied to the list singlepolys. I really don't understand why. Any help appreciated! **Edit** A clumsy work-around is possible (I've found a faster one but this is closer to the original thought). Here goes:               # ... Skipping initialization          # Open the attribute table of the shapefile     rows = arcpy.SearchCursor(in_fc)     for row in rows:         poly = row.getValue(shape_field)              # Check if single part (not multipart)         if not poly.isMultipart:             # Buffer this polygon only and add to list             geom = arcpy.Geometry()             result = arcpy.Buffer_analysis(poly, geom, '1000 Meters')             singlepoly_bufs.append(result[0])          # Now union the buffered features     out_geom = arcpy.Geometry()     union = arcpy.Union_analysis(singlepoly_bufs, out_geom)          # ... Skipping the end


title: What is the best way to change a data type in the ArcGIS attribute table?

text: What is the easiest way to change the data type in an attribute table in ArcGIS for multiple fields? Say if the fields have been exported as a string from Excel and I would like to use them as numerical values?


title: ArcGIS JS API - Can't draw results of Query Task

text: I have an ArcGIS JavaScript map that I'm trying to use with a Query Task. The query runs fine on the server side and the proper geometry is returned. Then, I draw the feature as follows:               var features = featureSet.features;          for (var i = 0; i < features.length; i++) {         var feature = features[i];         feature.setSymbol(lbSymbol);         feature.setInfoTemplate(lbInfoTemplate);         map.graphics.add(feature);     }      This was working for a while but it randomly broke and I can't figure out why. No errors -- the map just doesn't draw anything. If I listen for the `graphics-add` event on the graphics layer it fires and all the event data looks fine (including spatial reference, geometry, etc.) Can anyone think of why this wouldn't work or how else to troubleshoot? For what it's worth - this code to zoom to the extent of the feature is also broken:               var extent = esri.graphicsExtent(features);     if (extent) {         extent = extent.expand(3.0);         map.setExtent(extent);     }      The extent of esri.graphicsExtent(features) matches the extent of the feature. Not sure why the map doesn't zoom on setExtent. **UPDATE** : I noticed that graphics layer is referenced to WGS84, while the map and graphic feature are State Plane. I'm just using the built-in `map.graphics` layer -- I assumed this would match the spatial reference of the map automatically. Still puzzled.


title: How to count number of feature classes in geodatabase using ArcPy?

text: I would like to count the number of feature classes in a geodatabase. Some of the classes are in feature datasets as well. I have tried get count management, but that returns the number of rows in each feature class in the GDB.


title: Build a website to track GPS signals

text: I need to build a solution that allows me to track cellphones/tablets and display on a visual map. I think there are two parts, one is understanding the general architecture. The other is more detail. I have seen a few websites that do what i want to do.For example,the Android Device Manager lets you see where your devices are on a map. I'd envisage the tablets to be Android based so even if there is a SAAS solution that I could plug into my web page would be perfect Basically I was thinking of installing a generic app(possibly android) on a tablet/cellphone then capturing the gps data somehow through a web interface and displaying the tablets onto a single map on my webpage. A SAAS solution would work well or alternatively a javascript package to read the captured GPS co ordinates and display on a map. Any help is appreciated. Regards


I develop an application and I want to erase overlay features in a feature class in the application.I use ArcEngine 10. I don't have another layer that have erase features then i cant use erase tool in overlay toolset How to erase using feature(s) within the same feature class? is there any tools or programmatically way to do this? Thanks for advance ![enter image description here](http://i.stack.imgur.com/FDTYl.png)


title: Determining Natural Origin for a GeoTIFF file (Lambert Conic)

text: I am trying to display an existing GeoTIFF file (from the FAA WAC set) in an application that requires the GeoTIFF NaturalOriginLongitude tag (FalconView 5.1.1). There are some older files in that set that display OK, and newer files that don't appear. I'm confident that it is that tag (from running the code in a debugger, and from trying "fake" longitude values in files). That is, I know how to add the tag (using listgeo and geotifcp) as described in an earlier question, and adding the tag makes the file appear. However the values that I'm putting in aren't providing an appropriately georeferenced display. I've tried zero (because that seemed like something easy) and I've tried using the same value as is in FalseOriginLongitude (which is what the older files seemed to contain). Setting NaturalOriginalLongitude to be the same as the FalseOriginLongitude value is definitely better than setting it to zero, but its a little off in Longitude (a few minutes) and a lot off in Latitude (a few degrees). I tried some likely-looking google hits (especially an EPSG Guidance Note 7 Part2), but I didn't find anything that gives me the natural origin longitude (lambda0) value that I need. Here is the output of listgeo for a specific file that does not display.               $ listgeo CG-18\ WAC\ 42\ North.tif      TIFFReadDirectory: Warning, CG-18 WAC 42 North.tif: unknown field with tag 42112 (0xa480) encountered.     TIFFReadDirectory: Warning, CG-18 WAC 42 North.tif: unknown field with tag 42113 (0xa481) encountered.     Geotiff_Information:        Version: 1        Key_Revision: 1.0        Tagged_Information:            ModelTiepointTag (2,3):              0                0                0                              -833364.681326397491974.541710076 0                            ModelPixelScaleTag (1,3):               126.87142        126.87142        0                            End_Of_Tags.        Keyed_Information:           GTModelTypeGeoKey (Short,1): ModelTypeProjected           GTRasterTypeGeoKey (Short,1): RasterPixelIsArea           GTCitationGeoKey (Ascii,21): "PCS Name = CG-18 WAC"           GeographicTypeGeoKey (Short,1): GCS_NAD83           GeogCitationGeoKey (Ascii,24): "GCS_North_American_1983"           GeogAngularUnitsGeoKey (Short,1): Angular_Degree           GeogSemiMajorAxisGeoKey (Double,1): 6378137                     GeogInvFlatteningGeoKey (Double,1): 298.257222101               ProjectedCSTypeGeoKey (Short,1): User-Defined           ProjectionGeoKey (Short,1): User-Defined           ProjCoordTransGeoKey (Short,1): CT_LambertConfConic_2SP           ProjLinearUnitsGeoKey (Short,1): Linear_Meter           ProjStdParallel1GeoKey (Double,1): 38.66667                    ProjStdParallel2GeoKey (Double,1): 33.33333                    ProjFalseOriginLongGeoKey (Double,1): -118.085                    ProjFalseOriginLatGeoKey (Double,1): 36                          ProjFalseOriginEastingGeoKey (Double,1): 0                           ProjFalseOriginNorthingGeoKey (Double,1): 0                           End_Of_Keys.        End_Of_Geotiff.          Projection Method: CT_LambertConfConic_2SP        ProjFalseOriginLatGeoKey: 36.000000 ( 36d 0' 0.00"N)        ProjFalseOriginLongGeoKey: -118.085000 (118d 5' 6.00"W)        ProjStdParallel1GeoKey: 38.666670 ( 38d40' 0.01"N)        ProjStdParallel2GeoKey: 33.333330 ( 33d19'59.99"N)        ProjFalseEastingGeoKey: 0.000000 m        ProjFalseNorthingGeoKey: 0.000000 m     GCS: 4269/NAD83     Datum: 6269/North American Datum 1983     Ellipsoid: 7019/GRS 1980 (6378137.00,6356752.31)     Prime Meridian: 8901/Greenwich (0.000000/  0d 0' 0.00"E)     Projection Linear Units: 9001/metre (1.000000m)          Corner Coordinates:     Upper Left    ( -833364.681,  491974.542)  (127d51'15.84"W, 40d 3'23.84"N)     Lower Left    ( -833364.681,  -31877.551)  (127d16'36.93"W, 35d21'26.95"N)     Upper Right   (  685286.216,  491974.542)  (110d 2'34.02"W, 40d10'41.37"N)     Lower Right   (  685286.216,  -31877.551)  (110d31' 8.78"W, 35d28'20.04"N)     Center        (  -74039.233,  230048.495)  (118d55'44.97"W, 38d 4'18.53"N)      There are also world files for each TIFF. Here is an example:               $ cat CG-18\ WAC\ 42\ North.tfw     126.8714200000     0.0000000000     0.0000000000     -126.8714200000     -833301.2456163971     491911.1060000764      There are also HTML files providing metadata, but it seems (to me) to be a subset of the GeoTIFF tags. Any suggestions for how to get a reasonable value for NaturalOriginLongitude from the information I have available?


title: How to fix "The entry point for the procedure sqlite3_open_v2 cannot be found in the dynamic link library sqlite3.dll"

text: ### Translation: Hello, I have a problem installing version 1.0.8. The following error message appears: > The entry point for the procedure `sqlite3_open_v2` cannot be found in the > dynamic link library `sqlite3.dll` _[The next sentence is ungrammatical in the original and therefore its translation is in doubt.]_ Previously I had version 1.7 Wroclaw, which I uninstalled before installing the new one. Thank you for your help. ### Original Bonjour, J'ai un problème pour installer la version 1.0.8, le message d'erreur suivant apparait: > Le point d'entrée de procédure `sqlite3_open_v2` est introuvable dans la > bibliothèque de liens dynamiques `sqlite3.dll` J'avais avant d'installer la version 1.7 Wroclaw que j'ai désinstaller pour installer la nouvelle. Merci pour votre aide.


title: online map of 242,888 locations: which service to use? (what if temporal heat map?)

text: I have a 161.8 MB large CSV file with geocoded locations (with latitude and longitude) for 242,888 lines. I want to have an embeddable online map of all this. I wonder which service could deal with this amount of data, and what you'd recommend. Relevant details: — if the data looks too big (162 MB), I can cut down on fields, some are extraneous. But others would be nice to link to the points on the map, even if they are not strictly there for the mapping. (So people could navigate back to the original data source, e.g.) — I have date field for a temporal map. It would be great to use this to let viewers filter which points to show, or even better, show animated versions if possible. — one field would be good to use for showing intensity, hot spots on the map (viz. concentration of procurement contracts in space and time) — the data is the universe of procurement contracts in Hungary since 1998. So good Unicode support is a must, as is global coverage (though I think both are pretty standard these days).


title: My KMZ works for older versions of GE, but not the newer ones

text: I have a KMZ file which is a time animation. It used to work in the older GE, such as version 5 (I think). I can load the KMZ file in the newer GE, but do not see the animation control bar to play it. Something needs to be changed in the KMZ file, but cannot figure out. Here is the beginning and end of the KML file:               <?xml version="1.0" encoding="UTF-8"?>     <kml xmlns="http://earth.google.com/kml/2.0">     <Folder>       <name maxLines="1">ADCIRC_pbl</name>       <Snippet maxLines="1"><![CDATA[]]>       </Snippet>       <Folder><name><![CDATA[ahframe_legend]]></name><visibility>1</visibility><ScreenOverlay><visibility>1</visibility><overlayXY x="1" y="1" xunits="fraction" yunits="fraction" /><screenXY x="1" y="1" xunits="fraction" yunits="fraction" /><Icon><href>ahframe_legend.png</href></Icon></ScreenOverlay></Folder>     <Folder>     <name maxLines="1">ahframe_1: Image</name>       <Snippet maxLines="1"><![CDATA[]]>       </Snippet>       <GroundOverlay>         <TimeSpan>             <begin>2005-08-28T23:00:00z</begin>             <end>2005-08-28T23:30:00z</end>         </TimeSpan>         <name maxLines="1">ahframe_1</name>         <Icon>           <href>ahframe_1.png</href>           <viewBoundScale>1.0</viewBoundScale>         </Icon>         <LatLonBox>           <north>30.46</north>           <south>29.44</south>            <east>-88.3531</east>           <west>-89.6952</west>            <rotation>0</rotation>           </LatLonBox>        </GroundOverlay>       </Folder>     <Folder>     <name maxLines="1">ahframe_2: Image</name>       <Snippet maxLines="1"><![CDATA[]]>       </Snippet>       <GroundOverlay>         <TimeSpan>             <begin>2005-08-28T23:30:00z</begin>             <end>2005-08-29T00:00:00z</end>         </TimeSpan>         <name maxLines="1">ahframe_2</name>         <Icon>           <href>ahframe_2.png</href>           <viewBoundScale>1.0</viewBoundScale>         </Icon>         <LatLonBox>           <north>30.46</north>           <south>29.44</south>            <east>-88.3531</east>           <west>-89.6952</west>            <rotation>0</rotation>           </LatLonBox>        </GroundOverlay>       </Folder>     <Folder>      The KML file ends like this:               <name maxLines="1">ahframe_77: Image</name>       <Snippet maxLines="1"><![CDATA[]]>       </Snippet>       <GroundOverlay>         <TimeSpan>             <begin>2005-08-30T13:00:00z</begin>             <end>2005-08-30T13:30:00z</end>         </TimeSpan>         <name maxLines="1">ahframe_77</name>         <Icon>           <href>ahframe_77.png</href>           <viewBoundScale>1.0</viewBoundScale>         </Icon>         <LatLonBox>           <north>30.46</north>           <south>29.44</south>            <east>-88.3531</east>           <west>-89.6952</west>            <rotation>0</rotation>           </LatLonBox>        </GroundOverlay>       </Folder>      </Folder>      </kml>


title: GIS projection for non-geographic data

text: I am working on two different types of data. One is related to mapping of Genome structure and another is creating a layout of the inventories inside Internet Data Center. What I am confused about is what kind of projection do we use (if any) is such conditions. The reason using GIS APIs here is very helpful is that these projects also require the same layering of the information. So how do I create or structure the data? Should it also be shapefiles? What should their projections be? Can I use a shapefile to be hosted in Geoserver without projections? Please let me know if you need any more info.


title: Openlayers text layer is not working

text: I am trying to load markers from a text file using openlayers text layer, but no marker is shown, text file and icons are in same directory. However I am able to show markers using marker layer. My code is here and text file is also attached, kindly help me out.               function init()     {         /*     Create a new Map with Basic controls which are passed as     arguments     */     var map= new OpenLayers.Map('map_element',     {         controls:[         new OpenLayers.Control.Navigation(),         new OpenLayers.Control.PanZoomBar(),         new OpenLayers.Control.LayerSwitcher(),         //new OpenLayers.Control.Attribution()         ],         numZoomLevels:19,         units:'m'     }     );           // default lan/lat projection in which points will be provided      var epsg4326 = new OpenLayers.Projection("EPSG:4326");         /*      Create Basic Layers to show on map     */     // openstreetmap layer     var osm_map_layer = new OpenLayers.Layer.OSM (     'Open Street Maps'     );     // google map layer     var google_map_layer = new OpenLayers.Layer.Google(     'Google Street Maps',     {sphericalMercator: true,     maxExtent: new OpenLayers.Bounds(     -20037508.34,-20037508.34, 20037508.34, 20037508.34)}     );          /*     vector layer to draw custom shapes on map     */     var vectors = new OpenLayers.Layer.Vector("Vectors");     /*     markers layer will be used to insert markers on map     */     var markers = new OpenLayers.Layer.Markers( "Markers");         /*      layers are added to map     */     pois = new OpenLayers.Layer.Text( "My points", { location:"./pois.txt", projection:     map.displayProjection});     map.addLayer(pois);     map.addLayers(     [     osm_map_layer,     //google_map_layer,     //vectors,     markers,     pois     ]     );       //          var lonlat = new OpenLayers.LonLat(48.9459301,9.6075669);     lonlat.transform(epsg4326, map.getProjectionObject());     map.setCenter(lonlat, 15);     markers.addMarker(new OpenLayers.Marker(lonlat));                    } // last closing bracket      and text file lat lon title description icon iconSize iconOffset   48.9459301 9.6075669 Title One Description one   Second line.      (click again to close) Ol_icon_blue_example.png 24,24 0,-24   48.9899851 9.5382032 Title Two Description two. Ol_icon_red_example.png 16,16 -8,-8


title: Geoserver create SqlView Layer with Geoserver Manager

text: I would like to register a new sql view layer to geoserver through the rest api. I am using geoserver-manager. The registration is ok, the layer is working fine, except wms getFeatureInfo, which returns: "Operation on mixed SRID geometries". If i manually set the sql view's geometry column's srid in the geoserver administration page, it works fine. But i cannot find any way to set it in the xml request. I'm sending the geometry info:               <geometry><name>geom</name><type>Point</type><srid>4326</srid></geometry>      And the geometry attribute:               <attribute><name>geom</name><minOccurs>0</minOccurs><maxOccurs>1</maxOccurs><nillable>true</nillable><binding>com.vividsolutions.jts.geom.Point</binding></attribute>      But i cannot set the attributes srid and the geometry info is not enough. The geometry attribute will be registered as srid -1. Any thoughts i can overcome this? Is there an api call to refresh the layer fields? thx.


title: QGIS attribute values are changing when I save and stop editing a layer

text: I'm working with QGIS 2.0.1 on Windows XP. I'm creating a shapefile of points with one attribute column for a 10 digit phone number. I've tried setting the type as both string and integer, but I'm getting the same problem. Once I save the shapefile and reload it all of my 10 digit phone numbers have been replaced with " -2147483648 ". This keeps happening even when I create a brand new shapefile and start from scratch. It doesn't happen to any of my other values, just my phone numbers.. I'm about at my wit's end. -Thanks


title: Why does QGIS's Clipper tell me "Computed -srcwin falls outside raster size" and how do I address it?

text: I'm trying to clip a heightmap to a small area that I'm interested in. When I attempt to use Raster > Extraction > Clipper I get an error dialog with the message ""Computed -srcwin falls outside raster size of..." with the dimensions. I tried to select a smaller area but continue to get the error message and am not sure what it's complaining about and how to resolve the issue so I can clip the heightmap. What am I doing wrong or need to adjust? ![enter image description here](http://i.stack.imgur.com/4pIaa.png)


title: How to apply a graduated renderer in pyQGIS?

text: Just trying to run the following code already given in the pyqgis cookbook to create graduated colors symbology (as in example with 2 classes) I define the source info for the shapefile but when I run execute the code within my test-plugin, it just adds the layer without classifiying. Just adds as usual, all polylines appear same. Any solution or working code/method? Thank you link to code in cookbook http://www.qgis.org/pyqgis- cookbook/vector.html#graduated-symbol-renderer               from qgis.core import  (QgsVectorLayer, QgsMapLayerRegistry, QgsGraduatedSymbolRendererV2, QgsSymbolV2,QgsRendererRangeV2)          myVectorLayer = QgsVectorLayer('C:/my_test_shape_file.shp', 'test_shp_file', 'ogr')          myTargetField = 'target_field'     myRangeList = []     myOpacity = 1          # Make our first symbol and range...     myMin = 0.0     myMax = 50.0     myLabel = 'Group 1'     myColour = QtGui.QColor('#ffee00')     mySymbol1 = QgsSymbolV2.defaultSymbol(myVectorLayer.geometryType())     mySymbol1.setColor(myColour)     mySymbol1.setAlpha(myOpacity)     myRange1 = QgsRendererRangeV2(myMin, myMax, mySymbol1, myLabel)     myRangeList.append(myRange1)          #now make another symbol and range...     myMin = 50.1     myMax = 100     myLabel = 'Group 2'     myColour = QtGui.QColor('#00eeff')     mySymbol2 = QgsSymbolV2.defaultSymbol(myVectorLayer.geometryType())     mySymbol2.setColor(myColour)     mySymbol2.setAlpha(myOpacity)     myRange2 = QgsRendererRangeV2( myMin, myMax, mySymbol2, myLabel)     myRangeList.append(myRange2)          myRenderer = QgsGraduatedSymbolRendererV2('', myRangeList)     myRenderer.setMode(QgsGraduatedSymbolRendererV2.EqualInterval)     myRenderer.setClassAttribute(myTargetField)     myVectorLayer.setRendererV2(myRenderer)     QgsMapLayerRegistry.instance().addMapLayer(myVectorLayer)       * * * Thank you so much Kelly. Your code works perfect. However I noticed that both the code (one I have typed in my first message and yours in the 2nd message) DO NOT work in QGIS v1.7.3 but WORKS in QGIS v1.8.0. I think this was a bug(?) that already resolved in v1.8.0 And one more question; Do you have any sample code for (natural breaks, Jenks) classification of the "numeric_attribute_field" based on given number of classes (i.e. no of classes will be a variable in the code, say "n", and I will pass that from the plugin GUI i.e. n = spinBox.value() ) Thanks.


title: Is it possible to install GeoServer on a web server?

text: Is it possible to install and use GeoServer with PostgreSQL on a website which is hosted by a web server provider? I am using a web server for my website. I want to add mapping feature on my website. I know PostgreSQL for storing my data like shp files (postgres db is provided by the webhoster) but I do not know how can I install GeoServer and use to fetch data from postgreSQL and display on my website using Openlayers / Leaflet.


title: Is anyone else working on PostGIS rasters for Ruby?

text: This is more of an answer looking for a question, really. If you're writing Ruby (and especially ActiveRecord) code and want to use those tasty pixels from PostGIS rasters, you've probably had to deal with eleventeen gazillion ST_Value() queries to get the values for the band, or squirt it into WKT and parse that, or to a TIFF on the disk. PostGIS 2.1.0 brings some new functions, in particular ST_DumpValues(), which allows you to get the whole 2D array of values for a raster band with one query. But AFAIK no one has an ActiveRecord adapter that allows that to be returned, much less a column of "raster" type. But there's a dodge: Using ActiveRecord's connection select_value method (good discussion here), you can code up a query as a string, and return the results -- that 2D array -- as a string. For a 2-row, 3-column raster, the string (if you specified a single band in your ST_DumpValues call) will look something like:               "{{1.0,2.0,3.0},{4.0,5.0,6.0}}"      A simple, fast pair of gsub! operations swaps the PostgreSQL "{" delimiters for Ruby's "[" ones, call eval() on the string, and Bob's your uncle: You now have the entire raster as a 2D array. It may not be elegant, but it certainly works! So my question: Is anyone doing something more sophisticated of which I should be aware? Here's the (admittedly simple-minded) Ruby code:               # Get and set raster band values out of PostGIS 2.1 rasters as Ruby arrays.     # This example assumes that the raster column is 'rast' and that the primary     # key for the table is 'rid'. Also, NULL pixels will probably tank the getter.     # Fixing those is left as an exercise :-)          class PostGisRaster < ActiveRecord::Base       self.primary_key = 'rid'       # Grab a band of a raster out of the database and convert it to a 2D array. Defaults to first band (1)       def self.to_a(tablename,id,band=1)         id = id.to_i         str = connection.select_value("select st_dumpvalues(rast,#{band}) from #{tablename} where #{primary_key}=#{id}")         str.gsub!('{','[')         str.gsub!('}',']')         eval(str)       end            # Update a raster band (default: first band) in the database with values from arr. Note that arr can be any       # subset of the size of the full raster. Corner is X,Y (or long,lat), with the upper left corner as [1,1].       def self.update_raster(tablename,id,arr,band=1,corner=[1,1])         update_str = <<-END         select st_setvalues(           (select rast from #{tablename} where rid=#{id}),           #{band},#{corner[0]},#{corner[1]},           array#{arr.to_s}::double precision[][]         )         END         update_str = "update #{tablename} set rast=(#{update_str}) where rid=#{id}"         puts "update_str: '#{update_str}'"         connection.update(update_str)       end     end


title: How to use POST method with Overpass QL and OpenLayers?

text: I am using OpenLayers and Overpass QL with GET method to retrieve data from the Overpass server. Some URLs tend to be very long and hence get rejected by the Overpass server. To get around this, I am thinking of using POST method to query the Overpass server using OpenLayers. Is there a way to use POST method with Overpass QL (not XML) in OpenLayers? And if yes, what (syntax) goes into the **url** property of OpenLayers.Layer.Vector?


title: Problems using the 'Split' tool with ArcGIS

text: I'm running into difficulty when using the `split` tool with ArcGIS 10. ![enter image description here](http://i.stack.imgur.com/a4J99.png) I am trying to split lines using a grid that I generated using the Fishnet command. Each gridcell has a unique text value as a name (FID converted to a string), which is my `Split Field`. However, I can't tell whether it is the size of my grid (300000 cells), or if it is due to no line being within a cell that causes the error. ![enter image description here](http://i.stack.imgur.com/CYgJE.png) I tried breaking my grid cell into smaller pieces of ~50000 cells, but I still got an error occuring. Then I tried selecting only cells that were intersected by the lines, and just splitting using those cells, and still got an error. The strange part about the error messages, were that the analysis would run for about 30 minutes or so, producing output, and then an error message would appear on screen, but no explanation - no error number or no further details. What I have no resorted to doing is splitting the grid into 500 smaller grids, and then splitting each subgrid on itself. Then I will loop through every cell, and clip the lines to each cell. The problem with this approach is that it is taking a long time (~20 hours so far). Are there any ways I can improve this process?


title: How to create Concentric Circles in QGIS?

text: I'm trying to build a planning tool for center pivot irrigation system. It involves selecting a field (polygon), a center point, and then an outer circle to set the outer limits of the system. From there you select the inner rings based on fixed length sections. to optimize the coverage area. Anyone have an idea where I can start?


title: Why isn't View>View Settings working in ArcGIS 10 ArcScene?

text: I am running ArcGIS 10 desktop verison with the 3d Analyst extension to view DEMs in ArcScene. I want to view the DEM from a specific xyz location, which should be achievable using the View Settings dialogue box in the View menu. When I click on View>View Settings nothing happens except the menu closes. I can't find any hint of why this isn't working in the ArcGIS help, forums or blogs. Does anyone know where I'm going wrong?