title: Using Openlayers for GeoRSS from remote website

text: I am trying to get features from remote GeoRSS file with Openlayers... But I am stuck. Openlayers can not read remote file. I have searched on the net but I could not find any solution... Here is the code:               map = new OpenLayers.Map('map');          var wms = new OpenLayers.Layer.WMS("Osgeo WMS",         "http://vmap0.tiles.osgeo.org/wms/vmap0", {         layers: "basic"     });     map.addLayer(wms);          var earthquakes = new OpenLayers.Layer.Vector("T", {         strategies: [new OpenLayers.Strategy.Fixed()],         protocol: new OpenLayers.Protocol.HTTP({             url: "http://brainoff.com/test/rss.xml",             format: new OpenLayers.Format.GeoRSS()         }),     });     map.addLayer(earthquakes);          var select = new OpenLayers.Control.SelectFeature([earthquakes]);     map.addControl(select);     select.activate();          var draw = new OpenLayers.Control.DrawFeature(         earthquakes, OpenLayers.Handler.Point);     map.addControl(draw);          function toggle() {         if (document.getElementById("toggle-id").checked) {             draw.activate();         } else {             draw.deactivate();         }     }     document.getElementById("toggle-id").onclick = toggle;     toggle();     map.addControl(new OpenLayers.Control.LayerSwitcher());     map.zoomToMaxExtent();


title: ST_Box2D in PostGIS 2.0

text: I am getting an error trying to create a bounding box of a geometry using `ST_Box2D()`. Is the function deprecated? If so, what is the equivalent function?               SELECT ST_Box2D(ST_GeomFromText('LINESTRING(1 2, 3 4, 5 6)'));


title: Some highways segments ignored by osm2po?

text: Where i can read why some roads/segments are disappear from routing? We are trying to comment all "wtr.deny" rules in config file, but nothing changed. Examples added as links bellow. Screenshot1 Screenshot2


title: When should I merge two or more sub-basins into one?

text: I've been using Archydro to analyze one of the unknown hilly regions with many watersheds. Unfortunately, ArcHydro generates so many sub-basins that it looks really unrealistic and weird. However, I know that basins can be merged using the 'Basin processing toolbox'. But my question is: Does anyone know based on what criteria contiguous sub-basins are merged? Thanks.


title: Looking for examples on how GIS techniques facilitate utility surveying and management

text: Anyone get some ideas about the interface between underground utility and GIS? Utility surveying is quite a new topic to me and I want some examples on how GIS techniques facilitate utility surveying and management. Thx.


title: Why is my new Shapefile layer not showing up?

text: I am having trouble with a new shape layer not showing up on my map. The attribute table is showing up fine and all of the points are in there, however visually, there are no dots showing on the map. I have made sure that the new shape layer was at the very top of the stack. I tried placing the new shape layer on a blank new project and it worked perfectly as it should. Then I re- opened the original problem project, deleted all of the existing layers, saved it as temp with no data in it at all, and then closed and re-opened the project and re-loaded the problem shape layer into the new temp project. It did not show up again. I have gone back to my original data and re-created the csv that I initially used to create the shape file (did this a couple of times), in case there was some spurious content in there, however, this had no effect, and the problem remains. Using QGIS 1.8.0 and Windows 7 64-bit. However I installed QGIS onto a Linux machine and the identical problem occurred. Hope someone can help, or has an idea what it might be...


title: Can PGrouting be used to identify which path fits a known distance between a start and end node?

text: I am involved in a research project that uses historic transport data for England and Wales mostly in the nineteenth century. One particular problem concerns coaching data where all that is known is the point of departure (start node), destination (end node) and distance. We will soon have a GIS of turnpikes and other main roads. Could PgRouting or such like be used to identify which path best fits the distance between the start and end node? I am a newcomer to transport GIS. Best, Max Satchell Research Associate Department of Geography University of Cambridge


title: Issue on using Pyscripter with ArcGIS 10.2

text: Can someone please let me know how I can set the **Pyscripter 2.5.3.0 x64** works with Python 2.7.3 which comes with ArcGIS?   what I did was:   1- I installed the ArcGIS 10.2   2- I installed Pyscripter 2.5.3.0 x64 which forced my to install Python 3.3   3- I installed Python 3   now I cant use the Pyscripter for ArcGIS Pyrhon. I already add the arcpy in to Pyscripter IDE Options--> Special package as: os, wx, scipy, **arcpy** but still not working. I also check the "desktop10.2.pth" at c:\python27\ArcGIS10.2\Lib\site-packages which looks like:               C:\Program Files (x86)\ArcGIS\Desktop10.2\bin     C:\Program Files (x86)\ArcGIS\Desktop10.2\arcpy     C:\Program Files (x86)\ArcGIS\Desktop10.2\ArcToolbox\Scripts      Can you please let me know how i can fix this issue?


title: Query to find the Intersection coordinates of multiple Polyon

text: I am using the below Query to find out the intersection coordinates of 2 Polygon and its working fine,but how to use or modify the query to find the intersection area coordinates for multiple(or N number of) Polygons. SELECT ST_AsText( ST_INTERSECTION( ST_GeomFromText('POLYGON((1.25 9.5,1.375 18.25,7.25 17.875,7.375 9.5,1.25 9.5))'), ST_GeomFromText('POLYGON((1.25 9.5,1.375 18.25,7.25 17.875,7.375 9.5,1.25 9.5))') ) ); Plz let me know if any other alternative to this.


title: Clipping a SqlGeometry

text: I have been messing around with the SqlGeometry in the Microsoft.SqlServer.Types namespace in .NET, and have been rather impressed. Can anyone provide me with an example of how I could clip a geometry by another geometry? I know there are other libraries (ex. NetTopologySuite) that do this type of stuff in .NET, but I am trying to do this without any additional dependencies.


title: Dissolve not working in ArcGIS 10

text: The dissolve tool is not working for me when I use ArcGIS 10. It works perfectly when I use version 9.3. Any help?


title: Openlayers snapping the whole figure/feature

text: is it possible to snap (drag and snap) a feature to another feature in the same layer? I found OpenLayers.Control.Snapping but it's snapping only the points.


title: How to get Lon/Lat values from the postgis raster type column

text: I have a tiff raster in my PostGIS 2.0 database. I want to get the value lat/lon and extra data like temperature from the raster type column in postgis table. I want the table with columns lon,lat,temperature etc.. from raster type column(rast). For each corresponding lat and lon threre is a temperature value. The tiff file details by gdalinfo as:               Driver: GTiff/GeoTIFF     Files: atm_ncep_gfs_temp_20130422_003.tif     atm_ncep_gfs_temp_20130422_003.tif.aux.xml     Size is 720, 361     Coordinate System is:     GEOGCS["WGS 84",     DATUM["WGS_1984",     SPHEROID["WGS 84",6378137,298.257223563,              AUTHORITY["EPSG","7030"]],            AUTHORITY["EPSG","6326"]],         PRIMEM["Greenwich",0],     UNIT["degree",0.0174532925199433],     AUTHORITY["EPSG","4326"]]     Origin = (-0.250000000000000,90.250000000000000)     Pixel Size = (0.500000000000000,-0.500000000000000)     Metadata:     AREA_OR_POINT=Area     latitude#long_name=latitude     latitude#units=degrees_north     longitude#long_name=longitude     longitude#units=degrees_east     NC_GLOBAL#Conventions=COARDS     NC_GLOBAL#GRIB2_grid_template=0     NC_GLOBAL#History=created by wgrib2     time#_FillValue=9.999e+020     time#long_name=verification time generated by wgrib2 function verftime()     time#reference_date=2013.04.22 00:00:00 UTC     time#reference_time=1366588800     time#reference_time_description=forecast or accumulated, reference date is fixed     time#reference_time_type=3     time#time_step=0     time#time_step_setting=auto     time#units=seconds since 1970-01-01 00:00:00.0 0:00     TMP_2maboveground#_FillValue=9.9990003e+020     TMP_2maboveground#level=2 m above ground     TMP_2maboveground#long_name=Temperature     TMP_2maboveground#short_name=TMP_2maboveground     TMP_2maboveground#units=K     Image Structure Metadata:     INTERLEAVE=BAND     Corner Coordinates:     Upper Left  (  -0.2500000,  90.2500000) (  0d15' 0.00"W, 90d15' 0.00"N)     Lower Left  (  -0.2500000, -90.2500000) (  0d15' 0.00"W, 90d15' 0.00"S)     Upper Right (     359.750,      90.250) (359d45' 0.00"E, 90d15' 0.00"N)     Lower Right (     359.750,     -90.250) (359d45' 0.00"E, 90d15' 0.00"S)     Center      ( 179.7500000,   0.0000000) (179d45' 0.00"E,  0d 0' 0.01"N)     Band 1 Block=720x2 Type=Float32, ColorInterp=Gray      Min=211.200 Max=311.600      Minimum=211.200, Maximum=311.600, Mean=277.878, StdDev=21.574      NoData Value=9.9990002605540088e+020      Metadata:     _FillValue=9.9990003e+020      level=2 m above ground      long_name=Temperature     NETCDF_DIMENSION_time=1366599600     NETCDF_time_units=seconds since 1970-01-01 00:00:00.0 0:00     NETCDF_VARNAME=TMP_2maboveground     short_name=TMP_2maboveground     STATISTICS_MAXIMUM=311.60000610352     STATISTICS_MEAN=277.87782601544     STATISTICS_MINIMUM=211.19999694824     STATISTICS_STDDEV=21.573645315608     units=K      I have: postgresql-9.2 postgis-2.0 raster support raster imported into PostGIS with SRID 4326 How can i do this? Thanks in advance!


title: In Geoserver, how to give different colors of features within a single shapefile based on the value of a field?

text: In arcgis, this can be easily done by using the Symbology tab in the properties.I need to achieve a similar effect in Geoserver. How can I do it?


title: How to wire ESRI.ArcGIS.Desktop.AddIns.Extension to a custom AddIn Button click event?

text: I am developing custom ArcMap toolbar with various AddIn tools, button and extension. I see the sample how to access the ActiveView events in the extension Initialize() method.                private void Initialize()         {           // If the extension hasn't been started yet or it's been turned off, bail           if (s_extension == null || this.State != ExtensionState.Enabled)             return;                // Reset event handlers           IActiveViewEvents_Event avEvent = ArcMap.Document.FocusMap as IActiveViewEvents_Event;           avEvent.ItemAdded += AvEvent_ItemAdded;           avEvent.ItemDeleted += AvEvent_ItemAdded;           avEvent.SelectionChanged += UpdateSelCountDockWin;           avEvent.ContentsChanged += avEvent_ContentsChanged;                  ....                   // ? Possible to access AddIn button even here ?                   myCustomAddInToolbar.myCustomAddInButton.RaiseCustomEvent += HandleCustomEvent;       How can I wire my custom AddIn Button? I don't know how to access the AddIn button events from here. Is there any ArcMap(ArcObjects) Interface to access AddIn Toolbar/Button Events? Thanks


title: Freeware alternatives to eCognition?

text: Does anyone know any good freeware alternatives to eCognition? I need to do image segmentation and classification. I tryed the ecognition trial and it was really good. I'm doing the job with SPRING, but it has a lot of limitations... Is there any qgis plugin? Cheers


title: How to display previous and current year with dynamic text?

text: Is there a way to display previous and current year with dynamic text? If so, what is the formatting? I'm using ArcMap 10.0


title: How to Georeference Photos With No Spatial Reference?

text: I have about 100-200 field ground shot photos I would like to overlay in ArcGIS 10 that have no spatial information which I need to somehow develop using a photogrammetric formula. The camera was attached 2 meters above ground on a RTK GPS pole. One-to-One relationship between GPS points and photos. Accuracy is not a huge concern, but should at least be within 1 meter (which is approx. the field of view of each photo). The goal is to develop a way to auto-georeference each image using a script or model since it is possible to know the footprint size of each image in the coordinate system I'm using. **What I have to work with:** 100-200+ Ground-Shot Field Photos:               -camera height 1.228 meters (from tip of gps)          -camera angle ~90 degrees          -camera focal length 0.006 meters          -North orientation      GPS points shapefile (100-200+ records; 1 record=1 point=1 photo):               -XY location for each photo(approximately the centroid of photo)     -NAD83 UTM Zone 14N (meters)      **Basic procedures off the top of my head to accomplish without spending a ton of time georeferencing each individual photo:** *1) Determine the dimension (in meters) of one of the photos. I can calculate Photo Scale (PS = f/H) and Representative Fraction (RF = 1/(H/f)), but I'm needing to calculate length DE in order to find the equivalent corresponding length on ground AB (aka. the width (in meters) I will make for rectangular buffers in the next step).  ![photo scale](http://www.geog.ucsb.edu/~jeff/115a/lectures/scale_area/diagram.jpg) 2) Use those dimensions to create rectangular buffers around the GPS points 3) Spatial Join GPS points to the rectangular buffers (to get the picture name corresponding to each point) 4) Split rectangle buffer layer by attribute (using picture directory/name field) *5) Replace color fill of each rectangle buffer with the corresponding picture I'm having trouble with the photogrammetry formulas to use in #1. #5 is more of a hack to visualize the pictures to avoid the georeferencing situation. I don't expect anyone to 'solve my problem'. I like to post problems like these to a) collect my thoughts b) to possibly edify others. Your thoughts are appreciated. **UPDATED DETAILS:** -photos are ground shots -North-oriented -while this problem may seem impractical due to image size vs. scale, I have my reasons :D **Update: 6/20/2013** I calculated the geometric dimensions of each photo (1.52m x 1.14m). I then created rectangular buffers with those dimensions at each gps point centroid. I georeferenced one image by hand and it fits perfectly. **Current problem:** Using python, I want to scale down each photo to fit the same dimensions of its correpsonding polygon footprint.I looked at using warpfromfile tool in arcgis 10.1, but need to specify to and from ground control points for both the photo and the footprint which seems a bit tedious. Is there a tool that will auto-scale/fit a raster to the dimensions of a polygon?


title: Extract parcel layers from vector pdf

text: So, I have an interesting problem... I have a PDF file which has over 100 pages and each of these pages is displaying a map in vector format. By using some third party program I was able to find out that one of these layers is a parcel layer. Now, what I am trying to do here is extract just this layer for every page in this PDF document and then stitch each of the parcel layer pages back together. I am looking for some program or script which will automate some or all parts of this without screwing up the parcel line geometry. Additionally, I am assuming that I probable will need to manually align each of the extracted vector parcel pages which I guess is fine but what I think I need for this to work is some type of basemap background, I was thinking if one script can extract all the parcel vectors and then some other program can load all of them into some mapping program, me manually alligning each page and then finally being able to "join" the seperate pages into one large parcel vector map that the same (or yet another gis program) can append the correct lat/long data and turn it into a shapefile. Uggghhh, so... Does anyone know of a solution for this?


title: gdal_retile.py and gdal_translate Requirements

text: I basically need to use two GDAL utilities: gdal_retile.py and gdal_translate, on a Solaris 9 machine. I have tried compiling GDAL from sources, but this seems not to be working. So, since I have python on my server machine, can I just use the utilities without compiling gdal or these utilities depend on GDAL?


title: How to store an image in PostgreSQL and display it in an openlayers popup?

text: This is question regarding storing of image in PostgreSQL and display in openlayers popup through geoserver(using wms service). PostgreSQL is installed in my virtual machine. I tried to store image in my spatial table where I added new image column(bytea) then by using following query               update latlong set image = 'c:\\img.jpg' where gid = 1;      I could be able to store image but I think image not stored but its path. Now when I click on my overlays map then I can see other attributes plus image attribute also in popup but no image just some encrypted code which I did not be able to understand. So my first question is that if that's the correct way to store image so that I can get image attribute with image using wms service of openlayers. and second question if this is not right way then please give me some information how to accomplish that.I want to display image in popup with other attribute by using database which is PostgreSQL. Please if anyone can help I am desperate here...Thanks in advance. Best Regards, ADK


title: How to join edge of International Date Line in world file?

text: I am currently using world file (http://www.ngdc.noaa.gov/mgg/shorelines/gshhs.html) for international boundary information. Now I have difficulty with map projection that the map ends with International Date Line (IDL). Can anybody help me with method to join IDL together so that the projection of the map looks like google map where IDL is joined (like this-> http://gyazo.com/cd63a04362b844a728ad5b082a5bf532.png). Thank you for help.


title: How to get values from a raster attached to a line in ArcGIS 10.1

text: In MapInfo, using Vertical Mapper, I can do a line inspection and get the Max, Min, Mean etc values for a polyline from one or more open grids which is then added as an attribute to the polyline. I'm trying to do a similar process in ArcGIS 10.1 but have been unable to work out which tool I would need to run. Any suggestions, I have access to most extensions including Spatial Analyst and have ArcInfo level licencing.


title: Is it possible to always see all band values with Identify Features?

text: Does anyone know if it is possible to get the Identify Features tool in QGIS to show the band values for all the layers shown in the results window? To clarify, whilst I have several bands "showing" in the results window, the actual values for all but the topmost layer are folded away. Having to unfold the info each time I click on a new location is getting pretty tiresome.


title: Finding intermediate points and speed over a segment

text: I have a dataset which is in the form of ( Time , Car ID, Speed at that moment, Lat,Long) and I need to do some machine learning on the data to answer the following questions :   1. Best route from point A to point B    2. The speed at some point ( which can be on any segment) .  I am having a few problems because I have never dealt with this kind of data before ( I mostly worked with text before this ) :   1. The data is Hyderabad specific. I want to divide Hyderabad into segments but how do I go about dividing it based on points ? If I use 500m as minimum length and keep dividing , there can be many twists and turns which will go missing :(    2. If I get a random point , how do I decide which segment it falls on ? For example, I have two segments , each having start point ( lat,long) and end point ( lat',long') . If i get a third point , (l1,l2), how do I decide which segment it falls under ?  I actually have more questions, but if these get solved first, it would be awesome. Thanks a lot


title: PostGIS raster connection is not closed by MapServer

text: I'm using PostGIS raster in MapServer like this:               TYPE RASTER     DATA "PG:port=6432 dbname='mydb' user='postgres' password='1234' schema='public' table='truemarble' mode='2'"     PROCESSING "CLOSE_CONNECTION=ALWAYS"      My problem is that the database connection is not closed after a request. The PROCESSING directive should ensure the disconnection and this works fine with PostGIS vectors. Any ideas? I've tested with:   python_mapscript in a mod_python application   php_mapscript in a php-fpm application (Fastcgi) MapServer 6.2 and 6.4, GDAL 1.10, Ubuntu 12.04 and 14.04


title: Using postgis functions in ArcGIS Server 10.2 for Layer Definition Query

text: I asked the previous question Using PostGIS functions in ArcGIS Desktop for Layer Definition Query? And got an excellent answer. However this covered ArcGIS Desktop. In 10.0 it seemed that if I got this working for ArcGIS Desktop then it would work for ArcGIS Server. By this I mean, typing into the definition query dialog box would be the same as the following JavaScript               var layerDefs = [];     layerDefs[0] = "ST_DWithin(location, ST_GeomFromText('POINT(-97.52249000 35.47340000)', 4326), 0.1)";     dataLayer.setLayerDefinitions(layerDefs);      However when I use JavaScript for 10.2 in this manner it does nothing. The layer data stays exactly as before. There are also no error logs. What am I doing wrong here? Thanks


title: Changing $id in a shapefile in QGIS

text: I need to change the id (feature.id() in the python API or $id in the field calculator) of for features in shapefile layers in qgis, and I need to do it in Python. Is there a way? I tried something like what was suggested in the post How to change attributes with QGIS python? , but it did not work, because apparently id is not a field. In any case, is there a way of changing it? If not, is there a way to create a geographic feature and specify its id?


title: Find nearby cities by a driving distance radius

text: I'm studying GIS services to solve one specific need: Given x/y coordinates of one city, I want a JSON result that lists only cities (not facilities) that are nearby the specified city, within a driving distance radius (e.g. 100km). I came close to that with two: Geonames and ArcGIS, but I couldn't extract this kind of data on neither (maybe it's my lack of experience). Do you know a service that I can use for this specific need? (I don't want map rendering). If possible, can you give a query example using it?


title: Michigan Projection with GDAL or FWTools

text: Has anyone come across a specific projection for Michigan, specifically: NAD_1983_CORS96_StatePlane_Michigan_South_FIPS_2113_Ft_Intl When I look at spatialreference.org I do not find any instance of the projection, so I can locate the ESPG code. Thanks...


title: How to union polygons with common ID in QGIS?

text: I have a shapefile with lots of polygons and I need to union all polygons which have the same identifier-id into a new (multi-)polygon, ideally with the same identifier-id on another (new) layer. As there are lots of polygons I don't want to do this manually. Is there a simple way I can do this automatically in QGIS? I suppose This can be done with Spatialite, but that's only a guess. Please help a rookie. Any hints appreciated. Thank you!


title: How to calculate the bounding box by a given center and scale in Android?

text: Given these conditions:   * a scale like `1:50000`   * the center of the viewport is `100°E, 20°N`   * the size of the viewport is `400x600` How can I calculate the bounding box of the viewport? * * * **Update after atlefren's comment:** The geographic coordinate system of the map is EPSG:4490. We want to display them in a different projection like Mercator or latitude_longitude_projection(maybe this is the so-said un-projected). The size of the viewport is in pixels.


title: what is the alternative of setZoomLevel in vector based map

text: It is commonly used to generat map tiles at a certain scale and re-organize them for web map on pc. For example googlem map, and bing map. Then user will use the `setZoomLevel(int level)` or something else to set the zoom of the map. This is because the tiles is not continuous. Now I wonder what is the alternative method when use display the vector data ?


title: OpenLayers Problem | Unable to zoom Layer

text: I am trying basic examples of OpenLayers with my own Geoserver. When I tried to replace WMS with own local WMS path, I am unable to zoom layer after a specific scale.               var map;               function init(){         map = new OpenLayers.Map('map_element', {});                 var wms = new OpenLayers.Layer.WMS(                       'OpenLayers WMS',             'http:/localhost:7000/geoserver/test/wms',             {layers: 'test:builtup'},             {}         );                   map.addLayer(wms);         if(!map.getCenter()){             map.zoomToMaxExtent();         }     }      Please help me with possible error. ![This is the max scale](http://i.stack.imgur.com/LceUa.png)


title: What vector file type yields the greatest detail with the Heat Map plugin?

text: I'm using .shp files built from .csv files and having quirky results. Perhaps there is a better way to begin. The heat map appears to be accurately describing the data but calculates a huge area in relation to the data and so the cell size is way too large for smooth detailed map.


title: Combine hover info with onclick info in openlayers?

text: I am trying to do some complex thing. I wonna have:   1. popup with a name of a feature on mouse hover,   2. popup with more complex info ( name, foto, etc.) after mouseclick on the same feature - the same layer. It is possible to have those two things at the same time??


title: How to restrict zoom levels in the Google Maps API

text: Is it possible, in the Google Maps API v3, to restrict the zoom to specific levels (eg 16 and 7)? In my case I wish to show city points at a global scale, then show suburb polygons when zoomed in - but I only want those 2 levels of zoom rather than the full 20 levels. Thanks


title: How to develop a PHP GIS Applicationfor land parcel search?

text: I know this question has been asked numerous times but not from this angle. I am proficient in PHP and PHP Frameworks, my interest is now on GIS with the sole intent of creating a map solution where someone can query and know where a piece of land is located based on a UNIQUE Code. This can further indicate lease expiry and such. Basically a land sale/solution. My question is, how can I best implement this in GIS? Is google maps sufficient for this kind of system? For some reason I feel it takes something more in depth hence this. How do I go about achieving this goal? What are the prerequisites? I have zero experience in GIS, the most I have done is use google maps.


title: Discontinuous elevation contours at boundary of CanVec datasets

text: I am using Quantum GIS and I am working with the elevation contours from two adjacent CanVec datasets. The problem is that many of the contours' endpoints at the shared boundary do not coincide with their counterpart from the other dataset. I want to have these endpoints connected and then I want to merge the two datasets into one. Can anybody give any suggestions on how to do it?


title: Vector Legend item names in QGIS Map Composer

text: Is it possible to get the text for a legend item to go over multiple lines? I have some long item names and I have a set legend width that I have to use. Without being able to make the text go over multiple lines I am unable to get my legend to fit the page properly. Edit: I have added an image to better illustrate the problem. As you can see the text for area type one goes out of the box I have set up for a legend. Personally I think short names in legends work best, but project managers do not always agree with this. ![enter image description here](http://i.stack.imgur.com/BI67i.jpg)


title: select two markers & draw line between them in leaflet

text: I am very much new to leaflet. I have multiple markers plotted on my map in leaflet. now I have to draw line between two markers when I select them. Can any one help in doing this.


title: QGIS mixing up CRS

text: I´ve just started working with QuantumGIS for my thesis. I ordered shapefiles with the CRS: ETRS 89 UTM Zone 33N, but QGis layer properties now shows the layer has a different CRS (namely SWERE99 TM) and i did not change anything. I´m puzzled and can´t figure out what´s the problem. Thanks. Katja


title: Displaying network analysis service for Iran?

text: I am working on the ArcGIS API for JavaScript of ArcGIS for Server. I have a problem with network analysis services in Iran. There isn't any base map for routing in ArcGIS Online for Iran so how can write code for custom base map or dynamic layer? I create network data set but I don't know how can write code for it.


title: Converting GPS Positions to Google Maps

text: I have a set of GPS-coordinates in a database, which I want to compare to a coordinate supplied by the html5 location API. My coordinates look like this:                           <WGS84Coordinates>                     <Longitude>180325E</Longitude>                     <Latitude>591947N</Latitude>                 </WGS84Coordinates>      This is from Stockholm Central Station. If I enter these coordinates into google maps, I end up southwest of Stockholm in Google Maps. If I enter coordinates supplied by the html5 location API, I end up right next to Stockholm Central Station in Google Maps (which is my correct location). The same pattern applies to all the coordinates I supply form my GPS- coordinate-database. What is going on here? I've read that Google Maps uses WGS84, as does the coordinates from my database. I've also read that the html5 position API uses WGS84. Also, I've read something about Google Maps using Mercator Projection. I am unsure which projection the coordinates from the database has. Ideally, I'd want a converter to convert between the two formats. Any ideas?


title: Flow Direction Arrows in python

text: Any help would be appreciated. I am trying to write a python script to determine flow direction on a water pipe network. I have calculated the mid point of each section of pipe and created a point feature class. I now need to determine the flow direction. Originally used the ESRI cartographic line feature to do this but need the flow direction to be seperated from the polyline. When I drew the pipe network it was digitized in flow direction, so all i really need to do is access the from/to field of the polyline and populate the an angle field within the point feature class. ![adding flow direction arrows](http://i.stack.imgur.com/bY3D5.jpg)


title: Distribute ArcGIS services across multiple servers

text: What do you do when you have so many published services in ArcGIS that running even one instance of each service maxes out RAM/CPU? We've grown to the point where we have more published services than our current server can run at one time, and expect to continue to grow considerably over the next few years. Even at the minimum one instance per service, we're running out of RAM. Of course, doubling the RAM would get us by for probably another few months, but we will eventually need multiple servers, and I need to know how to work this out. My understanding is that when you add a new server to an existing site, all of your published services will have at least one instance per service on all servers. Is there a way of distributing these services across multiple servers without having multiple 'sites'? I konw that I could publish X amount of services to one server, and another batch to a separate server, but I'm assuming there's a better way to do this, and I'm hoping someone can shed some light on this. Thanks!!!


title: How do I get that polygon from Javascript ArcGIS API?

text: I would like to have a user draw a Polygon (or Line) onto a web map with the Javascript ArcGIS API drawing tools. For the life of me I can't figure out how to pull that polygon from the client and get it to my server. It baffles me. Help?


title: Controlling output linewidth linestyle using gdal_rasterize shapefile shp to Geotiff

text: I have successfully combined a shapefile (.shp) and a GeoTiff using gdal_rasterize. The shapefile displays as one pixel wide in the GeoTiff, but I need it to be wider and more visible, say, 5 pixels wide and red. Can anyone suggest an alternate way to achieve this? This needs to run as a script and not use desktop gis software.


title: How to split an oblong, irregular polygon layer into even pieces?

text: I have a polygon layer (landuse) that I want to split into 100 meter pieces. The polygons show landuse along a road. I want to split the polygons layer so that I can then calculate the percentage of every landuse per 100 meters. How do I make these splits? Any other ideas on how to make this calculation?


title: What are the default units of elevation/altitude/z in GeoJSON?

text: The GeoJSON spec says "The default CRS is a geographic coordinate reference system, using the WGS84 datum, and with longitude and latitude units of decimal degrees". It doesn't say seem to specify what the units of elevation/altitude/z are for points. Is there a default unit of altitude for GeoJSON? If not, are there common conventions, or a way to specify units?


title: Dividing the world space to geographic sub-space

text: I am trying to divide the world space to geographic sub-space from "general" to "specific" locations and without using latitude and longitude information. That is, I am trying to think about the world in someway like the following:               The most "general sub-space" > ... > The most "specific sub-space"      _Example_               Continent > State > Region > ... > City > Town > Address        1. What is (are) the most _used_ "division approach(es)"?   2. What is (are) the most _efficient_ "division approach(es)"?


title: Hyperlink relative paths to files in points layer in QGIS?

text: I would like to add hyperlinks to each point in a points layer in QGIS 1.7.0 running on Mac OS X 10.6.7, but I have several requirements: (1) Each hyperlink has to be a relative path to the linked file. (2) The linked files are mainly JPEG images, but may occasionally be HTML files. Either way, I want them to open upon click in the OS's default image/html viewer, regardless of what computer I open it on. (3) The linked filenames are already specified in a column in the point layer's attribute table. I want to use those filenames in they hyperlinks so I don't have to manually type in a path for each point (there are several hundred points, each pointing to a different file). How do I do this? Thanks.


title: Automate creating Rectified files from georeferenced images

text: A python 2.6.5 / ArcGIS 10 solution is preferred. UPDATE: I found a tool to do this "Create Ortho Corrected Raster Dataset" but the python script for it gives errors when I try.               import arcpy     import glob     import os          path = os.getcwd()     filename = path+"\Stage2_19"     #for pathname, directories, filenames in os.walk(path):     #    for filename in filenames:     #        if filename.lower().endswith(".TIF"):     arcpy.GetMessages()     arcpy.CreateOrthoCorrectedRasterDataset_management(filename+".TIF",\                                                        filename+".jpg",\                                                        "CONSTANT_ELEVATION", "0")          arcpy.GetMessages()     print "Created: " + filename + ".jpg"      Gives the following errors. > > > > Traceback (most recent call last): File > "Q:\Software\GIS\ESRI\ArcGIS10\CustomScripts\temp\Georef-Rectified.py", line > 16, in "CONSTANT_ELEVATION", "0") File "C:\Program > Files\ArcGIS\Desktop10.0\arcpy\arcpy\management.py", line 7815, in > CreateOrthoCorrectedRasterDataset raise e ExecuteError: ERROR 000425: Failed > to create ortho-corrected raster dataset Failed to execute > (CreateOrthoCorrectedRasterDataset). >>> >>> > ORIGINAL: I have georeferenced a few hundred survey plans using the Georeferencing toolbar in ArcGIS but did not Rectify the images (i.e create a new image) as I didn't want to create new files and take up more space. Now when I want to open them up in ERDAS or MapInfo they don't accept the prj created in the georeferencing process and require each using Rectify-->Save As Can I use python to automate the process - open the georeferenced Tiff files and save them as Rectified Tiff or JPG files with the same name. I noticed that there is an option for using gdal_translate but this gives me a syntax error and http://www.gdal.org/gdal_translate.html doesn't help to much (even when trying the example `gdal_translate -of GTiff -co "TILED=YES" utm.tif utm_tiled.tif`). Also not sure how it would handle the fact that the original file is only georeferenced and not rectified. Please help...thanks,


title: How can I transform GCS_WGS_1984 to GCS_North_American_1983?

text: I need to convert a GCS_WGS_1984 data set to NAD_1983_UTM_Zone_14N. What transformation should I use?


title: How to create network simulation / temporal models for videos?

text: I know of ArcGIS Tracking analyst and QGIS Time Manager. I have a transportation network with time based schedule for buses. I will create time aware shapefiles from this data. Is it possible to create a simulation that shows buses moving along a track (as per the time ) using some plugin? This is necessary as i want to convert it to a video. Essentially can you all list **all Add-ins / Plugins that deal with time based data** in popular GISystems.


title: Legend: Can I force long layer names to multiple lines based on a specified legend width?

text: Do to various limitations I almost always create static (graphic) legends but I do have the occasional map where impeccable cartography is secondary to functionality hence a dynamic legend is acceptable. I am Currently working on a map where I want the legend to be 4" wide with one legend column. I need to show very long layer names in the legend but **I don't know how to force labels to roll to multiple lines dynamically.** Thanks


title: HTTPS and working with non-secure services - ESRI jsAPI

text: **Situation:** I am currently working with an app that is required to be delivered over HTTPS because of an organizational security requirement. I currently do not have access to this servers administration, just the application code itself. **Problem:** The application relies on services from ArcGIS server (I have access to these servers). These services are only delivered via HTTP since we do not have the server configured for HTTPS. This causes issues with the app, since it is secure and requesting non-secure services. **Question** How do I get around this? I am guessing the best way would be to enable HTTPS on our ArcGIS Service server, but I would like to avoid doing that if at all possible. Can I use ESRI's JavaScript API Proxy rules to solve this (this would be optimal)?


title: Colorizing raster in QGIS 2.2 without using OSGeo4W command line?

text: Total beginner in QGIS. Following this guide right now: https://www.e-education.psu.edu/geog585/node/695 At step 4 I created the simple txt file which I now need to use to colorize the raster (slope.tif) I created before. Now it tells me to "open your OSGeo4W command line window" to execute some commands to do the colorizing. Where can I find this command window or a command line window for the GDAL Plugins? Do I have to install something separately? I wouldn´t mind to using command line windows. Or is there a way to do what is done in the tutorial just by using standalone QGIS Installation 2.2 (on Win 7 64bit system)? Perhaps by using the properties menu of the raster file? If yes please give some tipps for steps 4 to 8 of the guide. I guess the rest will be no problem.


title: In R, set NA cells in one raster where another raster has values

text: I have two rasters in R. I want to set values to NA in the first raster where the second raster has values. I think this should be simple, using the raster package, with two RasterLayer objects raster1 and raster2, both the same extent and snapped to each other. They are 29775x29930. I'm doing: newraster <\- raster1[is.na(raster2)] But this seems to take an unnecessary amount of memory, and keeps crashing R. My computer has 8GB of memory. Is there a less memory-intensive way to do this?


title: Can't get Tiger Geocoder to work in PostGIS

text: I've installed the PostGIS Tiger Geocoder according to these instructions and loaded in data for the state of New York. Unfortunately, when I try to geocode something...               ny_tiger=# SELECT g.rating, ST_X(g.geomout) As lon, ST_Y(g.geomout) As lat,      (addy).address As stno, (addy).streetname As street,      (addy).streettypeabbrev As styp, (addy).location As city, (addy).stateabbrev As st, (addy).zip      FROM geocode('350 Fifth Ave, New York, NY 10118') As g;                                                                      rating | lon | lat | stno | street | styp | city | st | zip      --------+-----+-----+------+--------+------+------+----+-----     (0 rows)      So that's not great. I noticed when loading things in I got a ton of errors like this:               Unable to open tl_2013_36123_addr.shp or tl_2013_36123_addr.SHP.     tl_2013_36123_addr.dbf: shape (.shp) or index files (.shx) can not be opened, will just import attribute data.      Is that the problem? How do I fix that?


title: Show graph in html popup for dynamic data of a postgresql table

text: I am trying to plot xy value from a postgresql database table. In the table there are multiple longitude and latitude corresponding to values. The table is like this, where all longitude and latitude of a is same same as b.                +----------+--------+-------+---------+--------+--------+---------+      | Time     | Long_a | Lat_a | value_a | long_b | lat_b  | value_b |      +----------+--------+-------+---------+--------+--------+---------+      | 12:10:04 |-94     | 30    | 1       | -95    | 34     | 9       |      +----------+--------+-------+---------+--------+--------+---------+      | 12:20:11 | -94    | 30    | 2       | -95    | 34     | 8       |      +----------+--------+-------+---------+--------+--------+---------+      | 12:30:23 | -94    | 30    | 3       | -95    | 34     | 7       |      +----------+--------+-------+---------+--------+--------+---------+      | 12:40:34 | -94    | 30    | 4       | -95    | 34     | 6       |      +----------+--------+-------+---------+--------+--------+---------+      | 12:50:44 | -94    | 30    | 5       | -95    | 34     | 4       |      +----------+--------+-------+---------+--------+--------+---------+      When I plot the data it plots fine in 2 different location point for a and b. But when I use the html popup and click on the icon of a or b it shows only one row not all 5 though they are in the same coordinate. My questions are,   1. How can I see all values of a same coordinate in the html popup ?   2. How can I create graph based on the time and value (Dynamic) of a and b points and show it in the html popup ? Edit: The image given below will make the question more specific, ![Required Output](http://i.stack.imgur.com/PcOz1.jpg)


title: Dissable layers in layermanager

text: I´m making an app using opengeo suite SDK and I added a LayerManager to show the layer list and the corresponding leyends. Layers load fine but now I want some of them off by default, added the "selected:false" option on those layers but it doesn´t work and the app loads with all layers enabled. Any idea how can I load the layers but disabled ?


title: Points to Line Feature in ArcGIS 9.3

text: I am writing to ask if there is the possibility to use extensions written for Arcview 3.2 in ArcGIS 9.3 With the older version of the gis package, i had an extension that allowed me to link together isolated points, in a large pool of coordinates, with the same attributes into a line. A typical example of my data is: -10 hours of navigation at sea made in 2 days resulting in, say, 1000 stored coordinates (500/day) -the daily navigation can be defined in different ways depending on the sea state (beaufort and duglas scales), weather condition, visibility and so on; my categories are, say, negative navigation "NN", positive navigation "PN" and off effort "X" (no navigation or transfer from a sector to another) -within the same day i can have all the three categories of navigation mentioned above, say in the following order: X, NN, PN, X, PN, X, NN, X -all the info about the navigation are stored in a single *.txt tab delimited file What I need, after importing and playing with the *.txt file in arcgis9.3, is a tool/extension that allows me to create a line, starting from the single points/coordinates, using those points/coordinates with the same attributes. say that i want to create a line that links those points that have the following characteristics: day, navigation status, beaufort scale, douglas scale, weather condition and visibility. if there is a series of 100 consecutive points with these characteristics they should be part of the same line. as soon as the next points, say, from 101 to 200 have at least one different value of the above attributes, they should be linked among them but not with the previous 100. basically, the final result should be a series of consecutive lines which i can know the attributes and the length. I could do all these things in Arcview 3.2 using a dedicated extension (i did not write it). Iow can i do the same things in arcvgis9.3 Many thanks for your help and advice.


title: How to use the URL for my Map Service in ArcGIS Server Manager?

text: I am having a Map Service from a ArcGIS Server Manager which I would like to add to the .cshtml page of my MVC project is this possible like the way shown in here by ESRI https://developers.arcgis.com/en/javascript/jssamples/index.html#latest_samples I am able to use the samples but how do I use my Map. Thank you in advance.


title: Accessing the esri js api over ssl

text: I changed the url i was using to fetch the esri js api to https://js.arcgis.com/3.7 and i can see in Fiddler that the call to that url is going out as ssl. But the file it pulls down (init.js) must be firing off other calls to fetch more components of the api and those are only going over http. Short of hosting the api myself is there a way to fetch the whole enchilada over ssl?


title: Generating a polygon tile index from multiple vector datasets

text: OK. well this is silly, in that there is no way it should be as difficult as it has been for me to find information regarding how to generate a polygon tile index from multiple vector (contour) tiles. With web publishing, tilecaches, etc. all the search results are not what im looking for. I have a folder full of about 50 contour shapefiles tiles that combined cover a US county. However there was no index distributed with the contour file collection, so I need to create one based on the bounding boxes, convex hulls, or some other function for each contour tile, and merge those together in a meaningful layer to be used to locate tiles necessary for a particular area... I'm failing to see a tool to do this in either QGIS or ESRI? I must be missing something so if someone could point me in the right direction, that would be great. Thanks


title: Qgis ,Area Management of Dealer Network - Analysis of Population versus dealer location and sales performance

text: I want to analyse a Dealers retail sales performance by location , measuring sales versus population in total and social class etc. I have the shape file for the population of Ireland which includes social class per area. Also I have the sales data per dealer per product category, and I know where they are located! So...now ,I would be very appreciative of any help where to start in Qgis, The way I envisage this is that I would have a map of Ireland with my dealers highlighted on the map with a specific radius , 30km for example , highlighted with a coloured circle. this radius would automatically estimate / calculate the population and it's demographics within the dealer area (30km radius) I would also like to be able to filter by product to try identify trends of a particular type of product sold per area. I Look forward to your help and hopefully patience as I am going this alone and am a beginner ,in the sense that I only installed Qgis yesterday. I installed the openlayer plugin and have my google street map of Ireland. Matt


title: I have a .sid image but can't figure out the CRS

text: There is a .map file:               OziExplorer Map Data File Version 2.2     Tiritir_2.sid     C:\OziExplorer\Maps\Islands\Tiritir_2.sid     1 ,Map Code,     Geodetic Datum 1949,WGS 84,   0.0000,   0.0000,WGS 84     Reserved 1     Reserved 2     Magnetic Variation,,,E     Map Projection,(NZG) New Zealand Grid,PolyCal,No,AutoCalOnly,No,BSBUseWPX,No      my best guess is NZGD49 North Island but when I use this projection the data from the GPS does not match. Any ideas?


title: Convert GPS coordinate to ArcGIS x,y point?

text: I'm trying to plot some `AGSSimpleMarkerSymbol's` in a `AGSGraphicsLayer`. The issue I'm running into is that the data I'm starting with are GPS coordinates (lat/lon). I'm using `AGSPoint` to plot the data, however, `AGSPoint` wants x and y points - not lat/lon. My question is, how do I convert my GPS coordinates to x and y points? Thanks in advance, Cole.


title: Loading GRASS gives error installing QGIS Chugiak on Windows 8.1?

text: I am experiencing a problem installing QGIS 2.4 Chugiak on Windows 8.1, resulting in a error loading Grass. This may be a problem only affecting certain installations, I'm not sure. I am a noob on QGIS and not so hot on IT. The installer QGIS-OSGeo4W-2.4.0-1-Setup-x86_64.exe does not report any errors, but when QGIS is launched a dialogue box reports that there is an error in finding the Grass installation and asking for the folder where this is located. Pointing to the folder results in it being unrecognised by the dialogue box, which asks the same questions again. This appears to have happened to others in the past (QGIS is not able to find GRASS. Choose the path of installation folder of GRASS) and occurs due to either of both of long path names not being recognised (or being correctly resolved to short names)and/or spaces in the file names.


title: How to automatically create all possible polyline routes connecting many point locations around a polygon obstruction?

text: Lets say I have 30 ports on the west of an island, and 30 on the east, and need to calculate all possible shipping distances between each and every port, with the routes having to steer around the geography of the island, which I have as a polygon. I could digitise by hand but this would take an immense amount of time, so ideally I would like to get ArcGIS to calculate these routes for me by connecting all the points with polylines, but with the polygon as a barrier the polylines have to avoid rather than simply connecting in a straight line over the land (polygon). As far as I can see in Network Analyst, if I had already had a network of potential routes, it would give me the shortest route option and take into account the polygon barrier. I could create a basic 'rough' network based on a vector grid, and get it to give me routes from this, but is there a more elegant solution? It doesn't seem like it should be an unusual task, but I can't find an existing solution? Any ideas? Thanks.


title: Geoserver, PostGIS and KML: Creating clickable nodes from OSM data

text: I have a PostGIS database that has OSM data imported into it using osm2pgsql and I'm using GeoServer to plot the points into a networked KML file for Google Earth. What I'd like to do is make the nodes that Geoserver pulls from PostGIS interactive. I'd like to be able to search for an element and have Google Earth query the PostGIS db for results. I'd also like the ability to click on a POI and have it pop up a dialog with all the attributes from PostGIS. So far I have a bunch of red points in Google Earth that the KML file added but they are not interactive. Any suggestions?


title: Elevation difference map using two DEM

text: I am studying glacier mass balance and have ASTER and SRTM DEMs. I want to make an elevation/volume change image such that it graphically represents changes in elevation difference in various parts of the glacier. If anyone can suggest an outline for how to perform it in ArcGIS 9.3, I would be highly beneficial for my project


title: Leaf Area Index (LAI) processing software

text: I took a bunch of hemispherical vegetation images with a fisheye lense in the field. Now I am looking for a software that can (semi-)automatically translate these images in LAI values. So far I struggled with CanEye but didn't get results. Anybody knows a good freeware solution?


title: Building a WebGIS with search functionality for a charity

text: I am playing with the idea of building a Basic Web based GIS system for a charity. The problem: I work for a charity that deals with people with Cleft Palates. They currently have a website, but need individuals to contact them to identify which branches of health teams / services / charities they have access to. I am wanting to develop an open sourced system that allows a user to **enter a Zip / Postcode, and for a list of the organizations that cater for there area to be produced.** Because this is a small charity finances are tight, so it has to be free / open source. Hence i am doing the leg work. I can digitise all the boundaries in question into .shp files. Also i can probably develop the front end. Thanks in Advance


title: Can not change layers on a WPF map

text: I am using WPF to display a map, which is created as MPK by ArcMap. I am trying to chang the DBF table and then update the map, but after I remove and then add new layers, the map won't change, the old layer is still there. What should I do to refresh the map? I think it may be because it is using cache or something. Thank you for any infomation.


title: How to install GRASS in QGIS Dufour?

text: How do I install GRASS to run in QGIS 2.0.1? GRASS is not listed in the Plugin Manager and attempting to access it via the Processing Toolbox gives me the following error: ![GRASS error](http://i.stack.imgur.com/xnC8J.jpg)


title: Graph of a specific area of USA

text: I am trying to get the network graph of a specific reason of USA. I have got a couple of options. One of them is using osm2po tool to extract the nodes and edge table from the osm map of a specific reason. I tried with this one and I could get the graph in the database table. But I think that there is property called geom_way, that defines the intermediate nodes of each edge. I mean suppose I have an edge like this               edgeid = 1;     source_node = 1;     target_node = 152     source_lat_lon = -93.27368,44.98204     dest_lat_lon = -93.27254,44.98124     geom_way = 0105000020E6100000010000000102000000040000006F641      Then if I decode the geom_way it will be the coordinates of the intermediate nodes between source_lat_lon and dest_lat_lon. Well I need to run the shortest path algorithm on this graph. So I was wondering if it could create some issues. Also I was wondering if there is any alternative tool or way to get the graph of a specific area of USA except osm2po.


Is there any way to add a google maps layer in a qgis web client project? I made a qgis project using a google maps layer from OpenLayers Plugin. In my desktop project, this layer works fine. But when I load the same project in the qgis-web-client all my layers are loaded, except the google maps layer.


title: Get coodinates from point shapefile

text: I converted a line shapefile to points as you can see in the attached image. Is there a way to know the coordinates of each point? Any help would be really appreciated! PD Im using QGIS 1.8.0 ![enter image description here](http://i.stack.imgur.com/9Hhkn.png)


title: Export selected converted geotagged photo points to a new feature class

text: I used the `GeoTagged Photos to Points` tool to convert a folder of geotagged photos to points, importing the photos into the file gdb as attachments. I now have a feature class of about 7 000 points. I have selected 100 of the points and now need to put this into a new feature class in a different gdb. However, the points lose their attachments when exported to a different gdb, and copying the original attachment table to the new gdb results in all 7 000 being copied. Trying to associated the match table with the new feature class using `Add Attachments` fails because the geotagged conversion tool automatically stores the image paths as relative not absolute. Is there a way for me to extract the selected photos along with their attached images to a new feature class?


title: How do I change the default location of a project in qgis-web-client?

text: we're using qgis-web-client to present our maps publicly over the net. I haven't worked with it a lot so I'm still playing around with it a bit. One thing I do not understand, yet: When qgis-web-client opens I just see a blank screen because the extents are not set correctly. But how do I change them? I also set `var useGetProjectSettings = true;` so actually qgis-web-client should zoom to the correct extents by default?! It's quite annoying to always click on the icon in the middle of the navigation arrows and then zoom to the correct location. Thank you for your help, Niklaas


title: How do I shorten the refresh rate of qgis_mapserv.fcgi?

text: Is there a way to shorten the refresh rate of qgis_mapserv.fcgi? After having uploaded a revision of a file these changes do not take effect in qgis-web- client. I guess it's because the computed tiles are cached by the qgis server. Is there any way to prevent this?


title: Detailed NTS Data for British Columbia

text: I'm looking to overlay a grid system on top of an open street maps map. The technology side isn't a big deal for me but I'm having some trouble finding the data I need. For Alberta and Saskatchewan I've been able to dig up appropriate DLS information to generate a grid down to the LSD level. I'm having more trouble with British Columbia. It is my understanding that BC makes use of the NTS. I have been able to find some data for it at http://ftp2.cits.rncan.gc.ca/pub/index/ but the files there seem to only map down to a 32km x 32km square. I assume there are at least two more detailed divisions below this level to come even close to 400m x 400m blocks which are denoted by LSD. Is there a publicly available source of data which might get me to the resolution I require? If you haven't noticed I have no idea what I'm doing with GIS so feel free to point out the numerous errors I've likely made in this question; I'm anxious to learn.


title: Configuration of QGIS Server with Apache Web Server on Windows 7 (64 bit)

text: I have installed QGIS Server(version 2.X) on my machine (Win. 7 64 bit OS & installed on C:\OSGeo4W64) while installing I found that for 64 bit Apache web server was not available hence its separately installed on C drive (C:\Apache24) Now tested Apache server installation and it works successfully in web browser for port 80. I am following this question & this link for installation reference so where should I put "OSGeoW64" folder in Apache web server ? to test `http://localhost/qgis/qgis_mapserv.fcgi.exe?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetCapabilities` this url ? Any help will be great !!! thanks in advance !!


title: Creating multiples buffers (or rings)

text: I would like know how could I creat and automate with a batch multiples buffers (or rings) between two cut buffers where is crossing roads. I creat rings using the attribute table but I can't automate and process and calculate the length of the roads inside these rings. Here I can show the photo of my result: ![enter image description here](http://i.stack.imgur.com/p0Lih.jpg) Best Regards!


title: Mann-Whitney U test for rasters

text: I'm doing a Mann-Whitney U test to compare two very large samples (sample size 1 = 13250; samlple size 2 = 38871) originating from a shapefile with 52121 features. I know t-tests are not recommended to compare rasters, because since rasters have so many values, they will almost surely detect a significant difference, no matter how small that difference might be (see this post in gis.stackexchange.com and point 3 of the first answer to this post). I ran the test on SPSS and got the following results:   * Test Statistics    * Mann-Whitney U 4094520,000   * Wilcoxon W 759591276,000   * Z -210,227   * Asymp. Sig. (2-tailed) ,000 While I did expect to find a difference between the groups, I don't know what to make of such large values. There is a paper which did the same thing I did. See the last sentence of the Material and Methods section. However, they obtained much smaller values http://icesjms.oxfordjournals.org/content/69/1/75/T5.expansion.html (they don't make clear if these are Mann-Whitney Us or z-scores). Granted they used a smaller sample (combined sample size = 5629), but the difference in the magnitude of the values still seem strange to me. **So, are my results simply the result of a very large sample, but still valid? Or should I use another test?**


title: How to get a signal when the geometry to a line in QGIS is changed

text: I am working on a plugin in QGIS and want to get a signal when the geometry to a feature(in this case a polyline) in my layer is changed. For example if I use the node tool to drag one of the nodes. In the qgis-api- documantion(http://qgis.osgeo.org/api/classQgsVectorLayer.html) i see that the line:   `void geometryChanged (QgsFeatureId fid, QgsGeometry &geom)` is listed under signals. So I used this command: `QObject.connect(self.iface.activeLayer(),SIGNAL("geometryChanged(QgsFeatureId, QgsGeometry)"),self.changedGeometry)` I have a method called changedGeometry:               def changedGeometry(self, intValue, qgsGeom):         QMessageBox.information(self.iface.mainWindow(),"Message", "Changed geometry. intValue: " + str(intValue) + "qgsGeom: " + str(qgsGeom))      When executed I get this error:               TypeError: C++ type 'QgsFeatureId' is not supported as a slot argument type      I have also tried `QObject.connect(self.iface.activeLayer(),SIGNAL("geometryChanged(int, QgsGeometry)"),self.changedGeometry)`. Then the program do not crash, but I don't get the message about changed geometry I am using QGIS version 1.8.0-Lisboa. Compiled against Qt 4.7.1 and these               import statements are in the top:                                                                  from PyQt4.QtCore import *         from PyQt4.QtGui import *        from qgis.core import *        from qgis.gui import *      Edit: I tried using qint64 instead QgsFeatureId, but it still don't work. `QObject.connect(self.iface.activeLayer(),SIGNAL("geometryChanged(qint64,QgsGeometry)"),self.changedGeometry)` Edit2: I also tried using the new style of connecting signals and slots. `self.iface.activeLayer().geometryChanged.connect(self.changedGeometry)` But it gave me this error message: TypeError: connect() failed between geometryChanged(qint64,QgsGeometry) and unislot()


title: Web browser capacity to display maximum records on Map using ArcGIS JSAPI 3.3/3.4

text: I have created ArcGIS JSAPI application in which I am displaying almost 50 thousand records at specific point. As usual, my browser hangs (until 30 thousand its showing but after 30 thousand it hangs) while displaying such large records so I want to know the browser capacity to display number of records. I am using ESRI 10.1 suit (ArcGIS Desktop and ArcGIS Server). I have changed default setting in ArcGIS Server to display records (default was 1000 and now its 80,000) In addition, I am searching on Google regarding the browser capacity to display number of records on map but it is not clearly mention plus ESRI support people also cannot given such article for ArcGIS JSAPI and I have search in the web help also unable to find any document. So is there any standard document or article that explains the maximum capacity to display number of records on map using ESRI JSAPI?? Tested on following Web browser: Firefox, Chrome, IE (very bad) and Safari (using latest version for all browsers) Thanks in advance : )


title: ArcGIS Server Security related justification

text: We just installed ArcGIS server(JAVA) and the network security has a few policy that requires us to provide justification on our setting. Does anyone have any reference document for the following questions? 1) Why must "arcgismanager" have access to the network folder? Basically they believe that noone should have access to the folder. Where can I find the document that states arcgismanager is required to the folder for configuration of arcgis server. 2) We have a batch job, and it requires ArcGISSOM, ArcGISSOC, ArcGISManager to have access to perform the batch job. Basically is the same as the above, we need to provide document that it is absolutely necessary for those 3 roles to have their ACL. where can I find the documents? 3) Remove Excessive Permission from CATALINA_HOME - conf, webapps, bin, logs Basically we set full control for all the above mention folder to all the 3 roles, they aren't happy with that and just want to bear minimum for the arcgis server(JAVA) to work. May I know if there is a set of mininum access right for those folders to the 3 roles? 4) 2.1.1 Upgrading of tomcat server in ArcGIS (ArcGIS JAVA). They want to know if it is possible to upgrade ArcGIS server internal tomcat servers -.- please kindly help me answer the above 4 questions with references to ArcGIS official documentation.


title: Point file to raster

text: I have a point file of trees in my study area, which includes the attributes of each tree. Each tree also has a crown radius associated with it. I need to summarize the attributes of these points into a 5-meter raster file. The 5-meter raster file placement is defined by a 5-meter fishnet. I would like a suggestion as to which method is the most statistically robust. One method I thought might work is: For a given 5 X 5 meter cell, take the area of each crown that occupies a cell and multiply it by the attribute of interest. Do this for each crown and sum these values and divide by the total area of crown occupying the given cell. i.e. Let the shape i be have the value of some attribute of the Tree … say Crown bulk density, (CBDi ) Then the value of the grid cell X = sum(CBDi * Ai) / sum (Ai) where Ai is the area of each crown within the cell. I don't know how to perform this method in ArcGIS though. Another simple method:   1. Do a point to raster and calculate the mean value and snap the raster to the defined raster.    2. Choose either bilinear or nearest neighbor for the resampling method in the environment settings. Thanks in advance for your help! Jim


title: psql:routing_dd.sql:28 ERROR

text: psql:routing_dd.sql:28: ERROR: could not load library "/usr/lib/postgresql/9.1/lib/librouting_dd.so": /usr/lib/postgresql/9.1/lib/librouting_dd.so: undefined symbol: __gmpq_equal Have no compile errors afaik or have seen. As soon as I'd like to import my DD extra functions it's not working properly: (dash) With Driving Distance psql -U postgres -f /usr/share/postlbs/routing_dd.sql routing psql -U postgres -f /usr/share/postlbs/routing_dd_wrappers.sql routing Versions: PostgreSQL 9.1 pgRouting - git'ed this morning (1.05?) CGAL 3.9 Any idea? Read that some similar isues have been fixed not long ago... Thanks! * * * Is not working! V1:   * Postgresql 9.1   * Postgis   * pgRouting freshly gitted   * all compiled from scratch \--> FAIL V2:   * Postgresql 9.1   * Postgis   * pgRouting   * all deps and libs downloaded from ubuntu repository \--> Still not working! Hopefully this will be fixed: https://github.com/pgRouting/pgrouting/issues/48#issuecomment-2789012 Thanks!


title: Getting Transparent Circles with Heatmap Plugin

text: With QGIS, I am trying to create a heatmap from a set of points I added as a layer using the "Add Delimited Text Layer" plugin. The points are a collection of complaints from around a city and I want to visualize portions of the city that have a higher frequency of complaints. Using the heatmap plugin, I choose the following options:   * Input Point Vector: my CSV layer of complaint points   * Output Raster: heat_map   * Output Format: GeoTIFF   * Radius: 1000 meters   * Decay Ratio: 0.1 I click OK and then choose NAD83, EPSG:4269 for coordinate system. Next, I get a large grey box with tiny transparent circles. I've tried switching the color map to pseudocolor and adding a transparent pixel list to make anything with a value of zero appear to be transparent. Then the whole map disappears. I'm not getting a heatmap at all and I've followed all of the popular guides I can find with Google. Advice would be greatly appreciated. Thank you!


title: How to show attribute table of object when user clicks on it?

text: sorry if I'm asking stupid questions but I really spent days trying to figure out how to add functionality to my map so that when user clicks on some map object he gets the attribute table on the side with the attribute data of that object. Is it even possible? I don't speak about popups, I speak about actual table. This is how my map looks now: ![map](http://i.stack.imgur.com/6Liqu.png) I would like that space on the east side of monitor populate with that attribute table when users click. I'm not saying that's the best solution but that's how I imagined that. But I'm open for suggestions. As you can see I use geoext, openlayers and extjs and honestly, I don't know lot about php or some other frameworks so I would be happy if the solution could be solved using geoext or smilair. I was looking around but nothing what I found is what I'm actually looking for. One more time, sorry if my question is unprofessional or something but I'm really kind of desperate :)


title: OpenLayers Plugin won't render maps correctly

text: In QGIS, I've created a map by adding a Tiger layer of my county, a KML of a water feature, and a geo-coded CSV of points of interest. I'm trying to add another layer using the OpenLayers Overview. When the Overview is open in my sidebar, it correctly identifies in its small window the location that I am viewing in the main window. However, when I try to add a layer, the I instead getting a repeating Google map that just appears as a band across the middle of the screen instead of a layer. Any ideas?


title: How to create schematic attribute as associated field using ArcObjects?

text: I'm trying to automate the complete creation of a schematic template with ArcObjects SDK for Java (10.2). I've managed to create the GDB, the schematic dataset and the schematic diagram class (template). Then, I created a schematic element class associated with a feature class taken from the same GDB, which works. I would like to share a field "NAME" between the feature class and the schematic element class, in order to label the nodes of a schematic diagram. So, I tried that:               ...      ISchematicElementClass nodeSchematicElementClass = schematicDataset.createSchematicElementClass("NODE", esriSchematicElementType.esriSchematicNodeType);      nodeSchematicElementClass.setSchematicDataSourceByRef(schematicDataSource);      nodeSchematicElementClass.alterAssociatedObjectClass(schematicDataSource, nodeFeatureClass.getObjectClassID());      nodeSchematicElementClass.setGeometryType(esriGeometryType.esriGeometryPoint);      nodeSchematicElementClass.setExternalQueryEvaluationMode(esriSchematicExternalQueryEvaluationMode.esriSchematicQueryBuildEvaluation);      nodeSchematicElementClass.setQueryString("SELECT * FROM NODE");      schematicDiagramClass.associateSchematicElementClass(nodeSchematicElementClass);           UID associatedFieldUID = new UID();      associatedFieldUID.setValue("{7DE3A19D-32D0-41CD-B896-37CA3AFBD88A}"); // taken from the file ...\Desktop10.2\bin\configuration\CLSID\esri.clsid.ecfg      ISchematicAttribute schematicAttribute = nodeSchematicElementClass.createSchematicAttribute("NAME", associatedFieldUID);      then I get the following ArcObjects error in the last line: > AutomationException: 0x80040a28 - Additional field can't be added : NAME in > 'Esri Schematics' I get the same error with any other field name. Does anybody have a solution? Am I doing something wrong with my attribute field association?


title: ESRI Javascript API Export to Shapefile and CSV

text: I am working on a Javascript Web Application using the ESRI Javascript API v3.3. My application allows the user to edit and query data and also provides a table for attributes. I would like the user to be able to export map data to a Shapefile or export the table data to a CSV. Is this possible with Javascript along? If not, my web application is running on Amazon's Elastic Beanstalk service with PHP. I've never written anything in PHP, but can this be done with PHP? How would I go about doing this in PHP?


title: Simplified flood "model" / estimation from elevation data

text: The goal is to present an estimation of potential flooded area. But still keeping it simple. I know this will not be accurate. The question is:   * will it be good enough for at gestimate / rough guideline? I am aware this is a crude method. I am using QGIS. ## Output I want to create flood area polygons of n m^3 of water. e.g. polygons showing the max. area that will be flooded if xx m^3 water is released. ## Input   * raster DEM   * m^3 of water to be released  ## Process   1. convert DEM raster layer to polygons based on elevation (in steps on 1m because of the DEM resolution).   2. set the area of the polygons as an attribute   3. create an volume attribute (area * 1m^3) to get an estimate of how much water the area will hold.   4. optionally rescaling the polygons to get estimations for volumes lesser than the polygon but greater than the "inner" smaller polygon (e.g. in fractions of 10).


title: Workflow for importing a flexviewer website into Flashbuilder 4

text: I have websites built using the flexviewer, but would like to import them into Flashbuilder 4 in order to enhance them. Is there a way to do this by importing the .xml or config.xml from the compiled flexviewer? Or is the workflow to build custom widgets in flashbuilder and deploy them to the flexviewer?


title: (When) should I use a GIS-enabled database?

text: I am a long-time programmer, new to GIS. I am trying to get a feel for the database aspects of GIS and I understand that there are some databases which are specialized for GIS use. Basically I am trying to understand whether to use such a database or to stick with the hugely popular, well established, widely supported (and free) MySql. The sort of applications I am liable to code would things like:   * fleet management (land or ocean based vehicles)   * employee tracking   * inventory control (at a granularity of a physical position accurate to with a metre or so)   * errm, that's about it, really I would normally expect to be tracking no more than 1 few hundred (maximum, possibly a few thousand) items. Sometimes the items will be located in the same (large) building, series of building, city, country or world-wide, depending on the application. I will occasionally be visually representing them either on a custom floor plan of a building or, more likely, in Google Earth or the like (more of which in another question). It seems to me that I am as well off with MySql and adding columns for lat/long or other positional data, but that might just be because I know MySql. Is there any reason why I should be looking at a more specialized database? Thanks very much in advance, and apologies for asking such a basic question. I _did_ google, but with no conclusive result, so thought that I would ask the experts.


title: GDAL bindings problem with ArcMap 10.1 Python tool

text: I'm really getting mad with this. ArcGIS sometimes can be so just...irritating! Forgive me if this question is by any means duplicating this one or similar ones, but I really can't find my way out of this problem. I have a x64 Windows Server 2012 machine (pretty cutting-edge, isn't it?) and a 32-bit ArcMap 10.1 installation on it, as well as a 64-bit ArcGIS Server installation. I came up with having two different Python 2.7 installations on my machine, both placed by ArcGIS:               C:\Python27\ArcGIS10.1\     #32 bit Python 2.7     C:\Python27\ArcGISx6410.1\  #64 bit Python 2.7      Now, what I want to do is to make an ArcMap tool which peforms calls to GDAL 1.9 Python binding library. My tool starts with a standard import of the gdal bindings:               import gdal     from gdalconst import *     [...]      I therefore installed (using gisinternals.com/sdk msi installers) BOTH 32 AND 64 bit GDAL core binaries and binding libraries, so I have:               #GDAL core binaries     C:\Program Files\GDAL       #64 bit     C:\Program Files (x86)\GDAL #32 bit          #GDAL bindings     C:\Python27\ArcGIS10.1\Lib\site-packages\osgeo     #32 bit     C:\Python27\ArcGISx6410.1\Lib\site-packages\osgeo  #64 bit          #PATH environment variable     echo %PATH%     C:\Program Files (x86)\GDAL;[...]      If I try to import the gdal bindings via both the Python interpreters, I get the followings:               C:\>C:\Python27\ArcGIS10.1\python.exe -m gdal     Traceback (most recent call last):       File "C:\Python27\ArcGIS10.1\lib\runpy.py", line 162, in _run_module_as_main         "__main__", fname, loader, pkg_name)       File "C:\Python27\ArcGIS10.1\lib\runpy.py", line 72, in _run_code         exec code in run_globals       File "C:\Python27\ArcGIS10.1\lib\site-packages\gdal.py", line 2, in <module>         from osgeo.gdal import deprecation_warn       File "C:\Python27\ArcGIS10.1\lib\site-packages\osgeo\__init__.py", line 29, in      <module>         _gdal = swig_import_helper()       File "C:\Python27\ArcGIS10.1\lib\site-packages\osgeo\__init__.py", line 25, in      swig_import_helper         _mod = imp.load_module('_gdal', fp, pathname, description)     ImportError: DLL load failed: %1 is not a valid Win32 application.          C:\>C:\Python27\ArcGISx6410.1\python.exe -m gdal     Traceback (most recent call last):       File "C:\Python27\ArcGISx6410.1\lib\runpy.py", line 162, in _run_module_as_mai     n         "__main__", fname, loader, pkg_name)       File "C:\Python27\ArcGISx6410.1\lib\runpy.py", line 72, in _run_code         exec code in run_globals       File "C:\Python27\ArcGISx6410.1\lib\site-packages\gdal.py", line 2, in <module     >         from osgeo.gdal import deprecation_warn       File "C:\Python27\ArcGISx6410.1\lib\site-packages\osgeo\__init__.py", line 21,      in <module>         _gdal = swig_import_helper()       File "C:\Python27\ArcGISx6410.1\lib\site-packages\osgeo\__init__.py", line 17,      in swig_import_helper         _mod = imp.load_module('_gdal', fp, pathname, description)     ImportError: DLL load failed: %1 is not a valid Win32 application.      and the very same error is given by the ArcMap Python shell when I try to run the tool:               Traceback (most recent call last):       File "C:\processors\URBAN_CLASS\URBAN_CLASS.py", line 2, in <module>         import GLCM_wrapper.GLCM_wrapper       File "C:\processors\URBAN_CLASS\GLCM_wrapper\GLCM_wrapper.py", line 2, in <module>         import gdal       File "C:\Python27\ArcGIS10.1\lib\site-packages\gdal.py", line 2, in <module>         from osgeo.gdal import deprecation_warn       File "C:\Python27\ArcGIS10.1\lib\site-packages\osgeo\__init__.py", line 29, in <module>         _gdal = swig_import_helper()       File "C:\Python27\ArcGIS10.1\lib\site-packages\osgeo\__init__.py", line 25, in swig_import_helper         _mod = imp.load_module('_gdal', fp, pathname, description)     ImportError: DLL load failed: %1 is not a valid Win32 application.      In addition, I tried to have the sys.path printed out from the ArcMap Python shell:               ['C:\\processors\\URBAN_CLASS',     u'c:\\program files (x86)\\arcgis\\desktop10.1\\arcpy',     'C:\\Windows\\SYSTEM32\\python27.zip',     'C:\\Python27\\ArcGIS10.1\\Lib',     'C:\\Python27\\ArcGIS10.1\\DLLs',     'C:\\Python27\\ArcGIS10.1\\Lib\\lib-tk',     'C:\\Windows\\system32',     'C:\\Program Files (x86)\\ArcGIS\\Desktop10.1\\bin',     'C:\\Python27\\ArcGIS10.1',     'C:\\Python27\\ArcGIS10.1\\lib\\site-packages',     'C:\\Program Files (x86)\\ArcGIS\\Desktop10.1\\arcpy',     'C:\\Program Files (x86)\\ArcGIS\\Desktop10.1\\ArcToolbox\\Scripts']      Any hint/idea on how to make my tool correctly import the gdal bindings? I'm completely stuck!


title: Counting occurrences across a record in ArcGIS 10.0

text: I'm a student at Millersville University, PA studying for a BA in Geography with a concentration in GIS. I've stumbled across this site for questions I've had over the past year or so with getting introduced to ArcGIS. So thank you all here for helping a student out :) I have come across a problem that I can't figure out and my textbook didn't offer any help so I'm hoping one of you can. I have a table which consists of 3 attribute joins, and I have a text field which has 2 unique values. I want to create a field in which I can summarize the number of times a particular text occours. For Example: Jx = Join #               ObjectID 1 |J1: A | J2: A | J3: B | CountFldA: 2     ObjectID 2 |J1: B | J2: B | J3: A | CountFldA: 1      I have a feeling that the answer lies somewhere in python but I have no formal experience in it. Thanks for any insight in advance! wxmeddler from Millersville.


title: Where to change GRASS GIS wxGUI language?

text: I use GRASS 6.4.2 on Russian version of Windows 7 x64. As I start wxGUI it is shown in Russian. But I don't like its translation and want to work with English interface. I have no clues how to solve it.


title: ArcGIS 9.3.1 Organizing point lables

text: I am creating a map with about 50 points and want to label the points numerically (i.e. 1, 2, 3, 4...) and provide a table with the corresponding point numbers in order with additional columns with info related to each point. My problem is I want the point numbering on the map to be grouped or organized so that the numbers are not randomly scattered all over the map (i.e. point 1 is near to point 2 etc) to assist in visual interpretation. My concern is if I just add a attribute field and asign numbers, on the map it will appear random. Any thoughts on how to address this? Thanks.


title: What is the ArcGIS equivalent to WEED tolerance?

text: Example:   User is digitising new (line, poly) features and accidentally adds additional vertices within very small distance of previous vertex, resulting in redundant vertices. I want to prevent this from happening within the current editing environment (ArcGIS 10.x). Back in workstation this could be accomplished using WEED tolerance parameter in the ARCEDIT editing environment. It seems there is no equivalent in ArcGIS. I've been offered a couple of answers to this question, both of which have been inadequate thus far:   **Answer - Change the XY resolution or XY tolerance of the feature class.**   _I absolutely do not want to do either of these as it creates geometry problems in the data. I don't want to change anything about the storage of the data itself, this would be like changing the coverage precision._ **Answer - Use the 'classic' snapping environment and snap to sketch vertices and set the tolerance to a minimum distance value.** _While this might eliminate some of these occurences, it does not actually **prevent** the user from creating redundant vertices. It also requires additional validation work afterwards to identify and fix possible errors (check geometry, topology etc.). I am trying to minimise errors at the editing stage while also reducing the amount of work to validate the data._   So the question boils down to "I used to be able to do this? Why can't I now?". The answer could well be "There is no equivalent, deal with it". I'm hoping this isn't the case though.


title: Is it possible to export GeoTIF current view to GeoTIFF with QGis?

text: ArcGis offers this, but I didn't found such a function in Qgis. Maybe I have simply not found the command.... Does someone know the command or is there an extenstion which offers such a functionality?


title: Setting default zoom level in Openlayers

text: GISers!!   With your help I was able get a shape file from the CAD file of the previous question. Now the next step is the display the shape file in Openlayers.... I wrote the script but the map being displayed in the browser is way too small... I need help setting the default zoom level... Can someone point me to an "simple" example that explains how to set the default zoom level?? The ones I've found a way to complicated.... Thanks   Chris


title: Has anyone use Network Links updates successfully with ArcGIS Explorer?

text: I've been trying to do updates of kml through a network link without success. Here is an example from the KML 2.1 tutorial. Load the following file (Point-load.kml) into AGX:               <?xml version="1.0" encoding="UTF-8"?>     <kml xmlns="http://earth.google.com/kml/2.1">     <NetworkLink>       <name>Loads Point.kml</name>       <Link>         <href>http://code.google.com/apis/kml/documentation/Point.kml</href>       </Link>     </NetworkLink>     </kml>      This should create two placemarks. Then load the following file (Update-load.kml):               <?xml version="1.0" encoding="UTF-8"?>     <kml xmlns="http://earth.google.com/kml/2.1">     <NetworkLink>       <name>Update</name>       <Link>         <href>http://code.google.com/apis/kml/documentation/NetworkLinkControl-Update.kml</href>       </Link>     </NetworkLink>     </kml>      This should update one of the Placemarks to a different name and does in Google Earth, but doesn't seem to work in ArcGIS Explorer. Cheers,   Alex


title: How to select all points within a circle in PostGIS?

text: I have set of points in a PostGIS Table (Type Geometry). Now, I need to filter points that lie in a circle with a center (C) and radius (R)! I expect C to be of type Point or Coordinates.. R to be in Meters!               SELECT * FROM my_points         WHERE ST_Within(point, ...)      However, I don't know how to complete it.. How can I exactly measure C and R?


title: Time Slider in CartoDB

text: I once again have a question regarding this map project I'm working on. I'm trying to add a time slider to my CartoDB map, and the slider appears, but instead of displaying the correct dates and the appropriate information, it's picking random dates. I've verified within CartoDB directly that the SQL query is correct. Below is my code - any thoughts? I feel it's probably something really simple that I'm missing, like a missed capital letter or something. Thanks!               <script>        function changeLegend(start) {         $('#legend').html(new Date(start))         }       function addTimeSlider(sublayer) {              var sql = cartodb.SQL({ user: 'ssci582-user12' })              // fetch time range              sql.execute('SELECT max(date), min(date) FROM History', function(data) {                var range = data.rows[0];                var max = new Date(range.max).getTime()                var min = new Date(range.min).getTime()                // update slider with range                  $("#slider").slider({                   range: true,                   min: min,                   max: max,                   values: [ min , max],                   change: function(event, ui) {                     // give feedback to the user on slide change                     changeLegend(ui.values[0], ui.values[1]);                   },                   stop: function( event, ui ) {                     // when user selects the dates, update the layer with the range                     var start = new Date(ui.values[0]).toISOString()                     var end = new Date(ui.values[1]).toISOString();                          // build sql                     sublayer.setSQL("SQL QUERY GOES HERE" [i've already verified it's correct]);                     }                 });                          changeLegend(min, max);              });          }          function main() {         var map;              cartodb.createVis('map', 'http://[username]/api/v2/viz/...json', {         shareable: true,         title: true,         description: false,         layer_selector: true,         center_lat: 35.79,          center_lon: -78.65,         zoom: 10     })          .done(function(vis, layers) {              layers[1].setInteraction(true);         layers[1].on('featureOver', function(e, pos, latlng, data) {             cartodb.log.log(e, pos, latlng, data);         });              var vislayers = vis.getLayers();         addTimeSlider(vislayers[1].getSubLayer(0));         map = vis.getNativeMap();              })              .error(function(err) {             console.log(err);         });          }      window.onload = main;


title: Display Lat/Lng coordinates from Oracle DB with ArcGIS HTML5 component

text: We have some experience in using ArcMap and ArcSDE, but this dates back to 2008, so we're not up to date with the recent developments in ESRI products. Now we're dealing with a database that has some Lat/Lng coordinates stored as plain `FLOAT` columns. The client wants an HTML5 map that displays some information based on these Lat/Lng coordinates: some point annotations, as well as some arrows between pairs of Lat/Lng coordinates. I proposed an architecture using Google Maps Javascript API v3 for the map component. I'd write a small web service that would deliver the Lat/Lng coordinates as `JSON`. In the browser, some small javascripts would display the Lat/Lng data in form of annotations and arrows. No big deal. However, due to high-level decisions by the customer's management, **the map component must be implemented using ESRI products,** so my idea with Google Maps and the simple JSON service won't work as intended. I see the following options:   1. **Basic scripting:** We're using roughly the same architecture as we would have used with Google Maps. So we need an HTML5 map component by ESRI that has a Javascript API, so we can display annotations and arrows via Javascript.   2. **WFS service:** We install/implement some kind of (OGC-compliant) WFS web service that reads the Lat/Lng values from the Oracle database and serves them as an OGC-compliant feature layer. Then, we're using an HTML5 map component by ESRI to consume that WFS service. No idea how I would implement the arrows between two Lat/Lng pairs.   3. **ESRI service:** Maybe ESRI provides a web service component with a generic data source, so we could wire that web service component to the Lat/Lng columns from the DB. This is basically the same architecture as in option 2, but with an ESRI service instead of some generic WFS component. Thanks for any hints. We're new to the HTML5 world of ESRI, so we're just looking for some **_what's-possible_** and **_how-to-get-started_** hints.


title: How can I insert a bounding box or polygon in a PostGIS geometry field?

text: I have a field in my table called 'geometry' that was created thus:               CREATE TABLE "public"."sites" (         "id" int8 DEFAULT nextval('sites_id_seq' :: regclass) NOT NULL,         #... A few MORE fields          "geometry" "public"."geometry",         CONSTRAINT "sites_pkey" PRIMARY KEY ("id")     ) WITH (OIDS = FALSE);     CREATE INDEX "index_sites_on_geometry" ON "public"."sites" USING gist (geometry);      I can create a point with the command:               UPDATE sites     SET geometry = ST_SetSRID (         ST_MakePoint (44.28, - 96.68, 495),         4326     ),     WHERE         ID = 12;      I would like to create a bounding box, but being new to the syntax have not been able to get the correct syntax. My best attempt:               UPDATE sites     SET geometry = ST_SetSRID (         ST_3DMakeBox (             ST_MakePoint (44.2, - 96.6, 400),             ST_MakePoint (44.3, - 97.0, 400)         ),         4326     )     WHERE         ID = 1123;      Gives the error: > Column has Z dimension but geometry does not I get this error even after trying   1. wrapping the `ST_SetSRID(...)` in a select statement as some online examples show (though the select statement returns a 2D result)   2. using `ST_Envelope(ST_Makeline` as suggested by @mattmakesmaps My apologies that I am new to this. I _have_ been looking through the documentation and googling `postgis insert bouding box` etc, but at my level, this isn't clear.


title: Script tool: Batch downgrade MXD version

text: I am trying to take the guts of this code (Batch convert ArcGIS 9x files to 10) and build it into a tool that we can more flexibly use and share around our organization. I think I'm almost there but when I run the tool as a script in Arc, nothing happens. Thoughts? Thanks for your help. I'm fairly new to Python but am pretty good on where the guts of it all is. mxdinloc is set to a _folder_ mxdoutvers is set to a _string - value list_ where the list choices are (10.1, 10.0, 9.3, 9.2) - no quotes EDIT: I ran through, got it working and wanted to share the final bit that I brought into my Script Toolbox:               #import modules     import arcpy     import os     import glob          #set parameters for inputs     mxdinloc = arcpy.GetParameterAsText(0)     mxdoutvers = arcpy.GetParameterAsText(1)          #set the values for input versions of ArcMap     if mxdoutvers == "10.1":         version = "10.1"         suffix = "_(10.1)"     elif mxdoutvers == "10.0":         version = "10.0"         suffix = "_(10.0)"     elif mxdoutvers == "9.3":         version = "9.3"         suffix = "_(9.3)"     elif mxdoutvers == "9.2":         version = "9.2"         suffix = "_(9.2)"          for pathname, directories, filenames in os.walk(mxdinloc):         for filename in filenames:             if filename.lower().endswith(".mxd"):                 print "\nUpdating: " + filename                 mxd = arcpy.mapping.MapDocument(mxdinloc + "\\" + filename)                 filename, ext = os.path.splitext(mxdinloc + "\\" + filename)                 print filename                 if mxd.dateSaved:                     print mxd.dateSaved                     mxd.saveACopy(filename + suffix + ext, mxdoutvers)                     currentMxd = filename + suffix + ext                     mapDoc = arcpy.mapping.MapDocument(currentMxd)                     print "Updated: " + filename + "\n"             else:                 print "\n" + filename + " is not a MXD and was not updated. \n"


title: Area surface doesn't match

text: I have three municipalities with a total surface of 178 sq km. I want to analyze land cover of those municipalities using both LANDSAT images with a pixel size of 30 x 30 meters and SPOT with a pixel size 10 x 10 meters. I have a raster of a whole county (> 1000 sq km) with land cover info derived both from LANDSAT and SPOT. When I clip the whole county raster with polygons of three municipalities, in case of clipping SPOT-based raster I get the correct surface of 178 sq km. That is the clipped raster area, the same as the surface of the three polygons/municipalities I am cutting with If I clip the LANDSAT raster, the surface is 25% smaller, according to the attribute table! I know there should be differences because the pixel size is not the same and on the edges of the polygon I am cutting with some small discrepancies will appear. But what happened with the 25% of the surface? I really have no clue what to do, I have done such things before, but never had such different results. Could it be something wrong with the data in the attribute table? (Even though I don't know what). Thank you in advance.


title: Import non spatial data into SQL Server using ogr2ogr

text: I am trying to import a CSV files in to SQL Server 2012 using ogr2ogr. The spatial data imports without any problems, but I cannot work out why the non spatial tables fail. I can import the data using other tools but it would be convenient if I could import everything using ogr2ogr so I can script it more easily. My VRT file looks like this...                   <OGRVRTDataSource>           <OGRVRTLayer name="test">             <SrcDataSource relativeToVRT="1">test.csv</SrcDataSource>             <FID>id</FID>             <Field name="id" src="Id" type="Integer"  />             <Field name="Test" src="Attrib" type="string" width="100" />          </OGRVRTLayer>         </OGRVRTDataSource>      And my command line:-                  ogr2ogr -append -f MSSQLSpatial "MSSQL:server=SQL01;database=TESTDB;trusted_connection=yes" test.vrt      The first error I get is:-               ERROR 1: Column ogr_geometry requested for geometry, but it does not exist.      After that there are errors trying to create the attribute field because the table does not exist. If I manually create the table in SQL Server with a single field named "id", then this command works (still with error messages).


title: Optimize Tilemill export of project with contour-lines (and fix out of memory problems)

text: I have a project in Tilemill (based on OSM-Bright) that covers Europe and also displays Contour-Lines. My problem is that the performance of Tilemill is pretty bad when displaying the contour lines and that the export (even only from a very small sample area) fails all the time with out-of-memory problems. My contours are in a postgres DB imported by `gdal_contour` and the table looks like this (default gdal_contour result, but I added the is_100m column and indexed it):               CREATE TABLE "public"."cont" (        "ogc_fid" INTEGER DEFAULT nextval('cont_ogc_fid_seq'::regclass) NOT NULL UNIQUE,        "wkb_geometry" "public"."geometry",        "id" NUMERIC( 8, 0 ),        "height" NUMERIC( 12, 3 ),        "is_100m" INTEGER DEFAULT '0',      PRIMARY KEY ( "ogc_fid" )      );     CREATE INDEX "cont_geom_idx" ON "public"."cont" USING gist( "wkb_geometry" );     CREATE INDEX "is_100m_idx" ON "public"."cont" USING btree( "is_100m" ASC NULLS LAST );      The `is_100m` column is used for displaying 100m lines. So in Tilemill I display them differently (a bit thicker and with labels). Is there any problem with my indexes? The machine I am running on has the following specs: OS: Debian Jessie/Sid in a VM (Host Kernel: 2.6.32-5 (Debian oldstable)) RAM: Host 100GB, the VM tilemill is running on: 59GB CPU: 2x Intel(R) Xeon(R) CPU X5650 (2.67GHz) -> 12 Cores. tilemill: v0.10.1-305-gb69b633 (github) nodejs: 0.10.29~dfsg-1 (Debian) postres: 9.3.4