Empty raster returned after GRASS i.atcorr with Landsat8 data

I've searched at length on previous forums but no solution seems to work for me. I'm trying to run the i.atcorr module on GRASS (version 6.4.3 on a Macbook Pro, Mavericks OSX) with Landsat8 imagery. Parameters didn't seem to be loading from the provided metadata file (MTL_txt), so I went back to the GRASS manuals and tried to input values by hand, here is what I used: 8 = geometrical conditions=Landsat ETM+ 05 04 03.81 98.117 23.131 = month day hh.ddd longitude latitude ("hh.ddd" is in GMT decimal hours) 1 = atmospheric mode=tropical 0 = aerosols model=none (ASSUMING NO AEROSOLS) 0 = visibility [km] 0.732 = aerosol optical depth 550 nm -1.408 = mean target elevation above sea level [km] -1000 = sensor on board a satellite 0 = THIS IS WHAT I'M WORRIED ABOUT: There is no pre-defined coding for Landsat 8 satellites and I am unsure how to define my own spectral conditions?? I checked I was in the correct region (i.e. used g.region) and I've tried using both radiance and reflectance at sensor values (computed using i.landsat.toar, which seems to have worked nicely) but the rasters returned are empty :( **Here is what I'm running (in verbose)** i.atcorr -a --overwrite --verbose iimg=Part1toar_B1@Charlie icnd=/Users/charlottewhitham/Documents/PhD/My research/my GIS/grassdata/NGH/Charlie/.tmp/Charlie.local/819.2 oimg=P1toar_atcorrB1 * **_*_ ** _*_ ** _*_ ** _*_ ** _*_ ***** 6s version 4.2b * * *   * geometrical conditions identity *   * \------------------------------- *   * etm+ observation * * *   * month: 5 day: 4 *   * solar zenith angle: 23.55 deg solar azimuthal angle: 103.24 deg *   * view zenith angle: 0.00 deg view azimuthal angle: 0.00 deg *   * scattering angle: 156.45 deg azimuthal angle difference: 103.24 deg * * *   * atmospheric model description *   * \----------------------------- *   * atmospheric model identity : *   * tropical (uh2o=4.12g/cm2,uo3 =.247cm-atm) * * *   * aerosols type identity : *   * no aerosols computed * * *   * spectral condition *   * \------------------ *   * constant *   * value of filter function : *   * wl inf=27640754752056197120.000 mic wl sup= 0.000 mic * * *   * target type *   * \----------- *   * homogeneous ground *   * constant reflectance over the spectra 0.000 * * *   * target elevation description *   * \---------------------------- *   * ground pressure [mb] 862.22 *   * ground altitude [km] 1.408 *   * gaseous content at target level: *   * uh2o= 2.056 g/cm2 uo3= 0.244 cm-atm * * *   * atmospheric correction activated *   * \-------------------------------- * Atmospheric correction... Atmospheric correction complete. (Fri Dec 6 12:16:10 2013) Command finished (22 sec)  **And here is a description of the output:** | Layer: P1toar_atcorrB1@Charlie Date: Fri Dec 6 12:12:58 2013 | | Mapset: Charlie Login of Creator: charlottewhitha | | Location: NGH | | DataBase: /Users/charlottewhitham/Documents/PhD/My research/my GIS/grassda | | Title: ( P1toar_atcorrB1 ) | | Timestamp: none | |----------------------------------------------------------------------------| | | | Type of Map: raster Number of Categories: 0 | | Data Type: CELL | | Rows: 7381 | | Columns: 7551 | | Total Cells: 55733931 | | Projection: UTM (zone 47) | | N: 2666715 S: 2445285 Res: 30 | | E: 522315 W: 295785 Res: 30 | | Range of data: min = NULL max = NULL | | | | Data Description: | | generated by i.atcorr | | | | Comments: | | i.atcorr -a iimg="Part1toar_B1@Charlie" iscl=0,255 icnd="/Users/char\ | | lottewhitham/Documents/PhD/My research/my GIS/grassdata/NGH/Charlie/\ | | icnd_lsat1.rtf" oimg="P1toar_atcorrB1" oscl=0,255 As you can see the range of data is NULL to NULL... not good. I would be very grateful if somebody could please help me with this - I'm really unsure what to try next! Thank you very much in advance and please let me know if I need to provide more info, Cheers, Charlie


Creating outline shp - georeferencing existing map/pictures

So I need to create shapefiles (that need to be adapted and transfered for other programs) with the district outline of a couple of Caribbean countries. Is georeferencing (using the _Georeferencing toolbar_ ) of pictures/maps online a good solution? Is there any other method I can approach, possibly more professional? What other measure should I take in order to make the shapefiles created can be used after (without basemap)? (Only thing off is Coordinate system, which in my case is WGS 1984) Realise this is a newbie question again, sorry, but help will be very much appreciated. Thanks


How to handle the New ESRI Ocean Basemap

Some of you may have seen this recent (4/28/14) ESRI blogpost - The New Ocean Basemap. ESRI, in their infinite wisdom, has split up their Ocean Basemap into 2 separate services. One with just the map tiles. And one with just the labels. And they will 'retire' the old basemap in 2015. My first thought: Cool, now I'll be able to drop the labels from their basemap if I want. My second thought: But what about those times when I don't want to drop the labels...? With basemaps being mutually exclusive, am I really going to need to build special circumstances into my mapping applications (to say whenever the ESRI Ocean Basemap is active, make sure this other ESRI Ocean Label overlay/layer is also active...etc). What a PITA. Are they really making it this difficult to display their Ocean Basemap **with** the Labels??? Have I mis-read their blogpost?   Am I missing the better/obvious solution? BTW I've created an openlayers example in a jsfiddle so you can see the 3 ESRI offerings in place (Original Ocean Basemap, New Ocean Basemap, and the New Ocean Labels).


Add x y data from existing attribute table

In Quantum GIS, I already imported an excel or dbf table. Now I need to tell Quantum GIS which columns are Latitude and Longitude. How do I do that? I know there is Add Delimited Text Layer which allows that, but it seems Unreasonable to me to convert XLS or DBF file to text file before importing! There must be some direct way of telling QGIS that _these_ parameters of the imported attribute table are actualy x y data.


How to restrict level crossings in route?

I want to evaluate if it is useful to bypass a level crossing through a nearby tunnel or overpass when the LC is closed or if it is worth waiting. I want to use the OSM data for realistic scenarios and calculate the most time efficient route for closed and opened LC. Is it possible to solve this problem with pgRouting? Or can you give me some hints where to start?


map stops at map.graphics.add(graphic);

function getLocRes(results) {     //map.graphics.clear();     var symbol = new esri.symbol.SimpleFillSymbol(esri.symbol.SimpleFillSymbol.STYLE_SOLID, new esri.symbol.SimpleLineSymbol(esri.symbol.SimpleLineSymbol.STYLE_SOLID, new dojo.Color([80, 0, 0]), 2), new dojo.Color([80, 0, 0, 0.5]));     var combinedExtent = new esri.geometry.Extent();     console.log("er",results.features.length);     for (var i = 0;i < results.features.length; i++) {console.log("i=",i);         var polyGraphic = new esri.Graphic(results.features[i], symbol);         console.log("polyGraphic",polyGraphic);         polyGraphic.setSymbol(symbol);         //map.graphics.add(polyGraphic);     }     combinedExtent = esri.graphicsExtent(results.features);     map.setExtent(combinedExtent, true);      } I made a function that shows multiple results in polygon. when I test it with out any map.whatever sort of working however, neither of map.graphics.clear() or map.grahpics.add() make the map work am I missing something here?


How to publish from ArcMap to GeoServer (not using GeoCat) - GeoCat Alternative?

How can I publish a map from ArcMap 10 to OpenGeo's GeoServer without using GeoCat Bridge? Are there any alternative solutions to accomplishing this? Thanks!


what criteria should I consider in buying a pocket GPS

It's my first time making this kind of choices, and I really don't know how to go with a choice or another. I basically know that it should have good **_durability_** , a good **_battery life_** and good **_coverage_**. What are other **_practical_** considerations I should be aware of ? And what about **_accuracy_** ? I actually need to have a pocked GPS to pick points from the electrical network (the nodes are the electrical posts), all to create a cartography of 1/25000 scale. I should cover 500 points or so. Thanks. P.S : this question is better being a wiki, I think.


OpenGeo Suite 3.1 doesn't work anymore when setting up Client SDK 3.0.2

We've installed OpenGeo Suite 3.1 on Ubuntu 12.04 LTS and everything works fine. When trying to set up the OpenGeo Client SDK 3.0.2 (following http://suite.opengeo.org/opengeo-docs/apps/clientsdk.html), however, the OpenGeo Suite doesn't work anymore. After inserting the SDK in PATH and adding ANT_HOME and JAVA_HOME in /etc/environment, tomcat6 starts up, but doesn't serve anything in the browser anymore. The severe error in the log looks as follows:               24 Jul 15:35:24 WARN [config.CustomEditorConfigurer] - Passing PropertyEditor instances into CustomEditorConfigurer is deprecated: use PropertyEditorRegistrars or PropertyEditor class names instead. Offending key [org.geotools.util.Version; offending editor instance: org.geoserver.platform.util.VersionPropertyEditor@3f705316          24.07.2013 15:35:24 org.geoserver.platform.GeoServerExtensions checkContext     SCHWERWIEGEND: Extension lookup occured, but ApplicationContext is unset.          24 Jul 15:35:24 ERROR [geoserver.global] -     - GEOSERVER_DATA_DIR: /usr/share/opengeo-suite-data/geoserver_data      Any help would be greatly appreciated. Are there conflicts between the versions? A new SDK doesn't seem to exist yet. EDIT: The error from the log doesn't seem to have to do anything with our problem, since it is also logged in another working geoserver environment.


display a choropleth map using D3 and Leaftlet

I am trying to get a leaflet map to display a choropleth map using d3. I have written this code using samples I have found however the map just displays as the black fill, not the color classes I defined. This is the code I am using. Anyone have an idea about what I'm doing wrong?               <!DOCTYPE html>               <meta charset="utf-8" />          <meta name="viewport" content="width=device-width, initial-scale=1.0">          <link rel="stylesheet" href="http://cdn.leafletjs.com/leaflet-0.5/leaflet.css" />                    <div id="map" style="width: 900px; height: 500px"


Find locations of raster maximum within polygons: Basin Delineation

I'm trying to identify the Location (lat/lon) and upstream area of every raster point that flows into the Arctic Ocean (there are a lot so it has to be fully automated). I'm starting with the following rasters:   * Flow Direction   * Accumlated Flow Area My approach so far has been to use the above rasters to create a delineated basin shapefile. From there I was hoping I could use the shapefile to look at the Flow Accumulation Raster and identify the location of the maximum value in each basin. In terms of output, I just need a table that shows [Basin#, Longitude, Latitude, UpstreamArea] for every basin. Any ideas or tips would be apreciated. [Note that I've already looked here but I don't know how to restrict the raster to only the area bounded by the polygons in the shapefile] * _Using: ArcGIS10, Python, etc._ Thanks, Joe


Is there are QGIS plugin to allow the 3d visualisation of geological borehole data similar to the functionality of Target for ArcGIS?

I'm very new to QGIS, I have started using it as my new company is not willing to pay for ArcGIS and Target at this stage so I'm hoping QGIS can substitute not just in the interim but for good. What I am looking for is a way to load in my borehole database and create 3d visualizations of the data, make cross sections and if possible model seams. Is there any capability within the QGIS world to undertake this kind of work?


Displaying information on a map using an InfoWindow

I am wanting to display information in an InfoWindow when the user clicks on a point on the map. I have implemented the InfoWindow in XAML, but I do not understand the code behind it. Any information will be appreciated. Here is my (what I feel is relevant) XAML:                <esri:InfoWindow x:Name="MyInfoWindow"                                  Padding="2"                                  CornerRadius="20"                                  Background="{StaticResource PanelGradient}"                                  Map="{Binding MyMap}"                                  ContentTemplate="{StaticResource MyFeatureLayerInfoWindowTemplate}"                              MouseLeftButtonUp="MyInfoWindow_MouseLeftButtonUp" />      If you need anymore information, let me know. I have the information being displayed correctly at the moment in a MapTip. Here is the XAML for that:               <esri:FeatureLayer.MapTip>                             <Border CornerRadius="10" BorderBrush="Black" BorderThickness="2" Margin="0,0,15,15" Background="{StaticResource PanelGradient}">                                 <StackPanel Margin="7">                                     <TextBlock Text="{Binding [WVAL_ID], StringFormat='ID: \{0\}'}" FontWeight="Bold" Foreground="White"  />                                     <TextBlock Text="{Binding [SIZE], StringFormat='Size: \{0\}'}" FontWeight="Bold" Foreground="White" />                                     <TextBlock Text="{Binding [WVAL_USE], StringFormat='Use: \{0\}'}" FontWeight="Bold" Foreground="White" />                                     <TextBlock Text="{Binding [WVAL_TYPE], StringFormat='Type: \{0\}'}" FontWeight="Bold" Foreground="White"  />                                     <TextBlock Text="{Binding [OPENS], StringFormat='Opens: \{0\}'}" FontWeight="Bold" Foreground="White"  />                                     <TextBlock Text="{Binding [TURNS], StringFormat='Turns: \{0\}'}" FontWeight="Bold" Foreground="White"  />                                     <StackPanel Orientation="Horizontal" >                                          </StackPanel>                                 </StackPanel>                             </Border>                         </esri:FeatureLayer.MapTip>


Using ArcGIS ModelBuilder with a "While" Iterator

I am struggling to get the ArcGIS ModelBuilder "While" Iterator to work as I would expect it. I have a hunch it's user error, but even reading the help hasn't seemed to get it working for me. I have a simple model that calls a Python script that simple does a `time.sleep(15)` and then sends back the conditional parameter of "TRUE" or "FALSE" using `arcpy.SetParameterAsText(1, "FALSE")` or similar for True statements in an if/else statement. I've also tried every combination of TRUE, True, and true, as outputs from my Python script without success. ![Model with WHILE Iterator](http://i.stack.imgur.com/G4R7u.jpg) I can see in the results that my conditions are being set correctly:               Executing (While): IterateCondition FALSE TRUE     Start Time: Tue Nov 20 16:10:06 2012     Succeeded at Tue Nov 20 16:10:06 2012 (Elapsed Time: 0.00 seconds)      and               Executing (While): IterateCondition TRUE TRUE     Start Time: Tue Nov 20 16:26:41 2012     Succeeded at Tue Nov 20 16:26:41 2012 (Elapsed Time: 0.00 seconds)      but the WHILE Iterator continues to fire the sleep script, even if I change the output parameter of my script. Any ideas on why I cannot seem to get the sleep script to stop firing? What should I be passing from my SetParametersAsText() in order to get the While Iterator to stop firing the Sleep script?


Is it possible to export cross sections from QGIS to HEC-RAS

Is it possible to create cross sections along a river in QGIS and then export the data to HEC-RAS? Thanks


Create Array From Selected Features Based on the Order in Which They Were Selected

**I need to create an array of features from a selection based on the order in which they were selected.** As in: I select feature 5,2,3,4. I need to have an array so I can apply attributes (with python) to to these features based on order in which they were selected. This is too apply part numbers to pre-mapped survey plans. To add a part to each is time consuming when their are 50 parts. This would allow me to select them in the order of part number, then apply the correct part with a button I have already created with python. I've been trying for a long time to figure this out and I have run out of solutions. I am pretty well versed in Python and with that ArcPy. Seem's like this should be easier then it is. Thanks for any help.


Link attribute value to Fanout in FME?

I would like to use attribute value from Workbench as an attribute that says Yes or No to the Fanout options in Workbench settings. So apply value from WB to the advanced fanout settings of WB. Any idea? Someone, please?


If statement in ModelBuilder to check if field exists

I'm trying to create the following in ModelBuilder:               Set input feature class       if field1 exists         execute the processes in branch1       else         add field1         execute the processes in branch2         merge branch      While I have no problem with the processes in the two branches, the `IF` check at the start doesn't seem to be working the way I want it to work. I am using this code in the `Calculate Value` code block               def ifFieldExists(Input_FC, SUMMARY):       import arcpy       fields = arcpy.ListFields(Input_FC)       for f in fields:         if f ==" SUMMARY":           return "true"         else:           return "false"      and expression               ifFieldExists(r"%Input_FC%", "%SUMMARY%")      This is the precondition for `branch1`, so if it returns `true`, the field exists so `branch1` executes. The other branch's precondition returns the reverse. When I run it on an input that does contain the field, the precondition returns `false` instead of `true`, so the wrong branch executes.


How to convert fake donut holes in polygons to real ones?

I have this problem where I have to import shapefiles into our system which must be geometrically valid, but the software that everyone uses to produce them seems to output files where every polygon that should contain a hole is instead made up from a single outer ring which goes around the boundary and then goes inward along some arbitrary path, traces around what should be the inner ring, and then doubles back on itself back to the outside. What I seek is a way to resolve these by somehow dissolving the vertices along the arbitrary path in and back out of the hole, to make a seperate inner and outer ring which will form a valid polygon. I have tried the trial version of ArcGIS and it seems to do this automagically when you snap the vertices on the double-back part together, but my business doesn't have a copy of ArcGIS and I might have a hard time convincing my boss to buy it, so I was hoping QGIS would allow me to do such a thing? Screenshot: ![description](http://i.imgur.com/zGwtG.png) I couldn't figure out how to get any WKT out of this sorry, but from playing around with the vertices I'm quite sure that it's one single ring that doubles back on itself for a bit.


How to compute a map-to-screen transformation?

Every projection calculation I have come across dealing with the Lambert Conic projection converts the given Lat & Long to meters. So now the magical question: **how on earth do you take the meters and convert to X, Y in pixels?** I have a map of the US. The following is the projection info:               ullat=53.3638     ullon=-130.964     lrlat=18.2977     lrlon=-69.1058          sp1=0     sp2=39     rlon=-95.5     rlat=36.5     width=640     height=480      I converted this http://svn.onelab.eu/deprecated.PLCWWW/imports/import-001-head-04-24/plot- latlong/plot-latlong over to php, but get returned meters. I am working on an open source project. Most of the users will not be able to install Proj4, so that is not an option (even there the return is in meters). Anybody?


Understanding null shape type

I have some shapefiles that I am reading and then loading into tables in the database. The shapefiles with no zoom (eg. "au_postcode") have no shapes with type null. When I tried to load a shapefile with zoom (eg. "au_postcode_100"), I saw some null shapes. It would be great if someone can please explain to me how the null shapes work and why would a shape file be null. Its kinda confusing because to me a shape is at least a point and if there is a physical location, there would be a shape for it. I am new to GIS and if something in this question does not make sense or I've used the wrong terminology, please let me know. Thanks


How to create a polygon from other polygon boundaries?

I want to write an algorithm that creates a polygon inside a hole made by a set of other polygons. It's something like a flood fill and it's like the "boundary" tool in Autocad. The image shows an example: the algoritm should fill the enclosed white space in the center (after a click on it). Any theory, pseudocode or examples would help a lot. ![enter image description here](http://i.stack.imgur.com/IIc4P.jpg)


Looking for a sample dataset of international addresses

I need a test dataset of 100 addresses in each country's local langauge. This is approximately 10,000 addresses. The don't need to be geocoded or resolved to address components. For example, 1 microsoft way seattle wash xxxxx, or 5814 pontiac dr san jose, ca 95123 Почтовый адрес: 125009, Москва, Театральная пл., 1. are comlete addresses with enough data. Similar ones for addresses in each world langauge would be fine. Can I get such a non geocoded, random data set for low or no cost? This is to test a google geocoded address cache (as per the API recommendations)


Link non-spatial SQL Server query to a Geodatabase?

Is it possible to build an object (of whatever kind) in ArcEngine code that links the geodatabase data with a non-spatial SQL Server query? How would you go about doing this?


DIY easy and ~accurate small-scale landscape survey

## Summary Is there any equipment that I can rent 'inexpensively' (less than $400) that would allow 1-2 people to accurately (+/- 6") capture the topography of hilly terrain of perhaps 250' on a side in an afternoon? ## Details My wife and I are about to embark upon a redesign of the landscape around our home in the mountains. Our lot started out on a generally 1 in 3 slope, and this has (of course) been changed in many interesting ways during building. There are small hills, patios, semi-flattened areas followed by large falloffs. As part of our designing, we want to capture exactly what the landscape looks like now, digitally. (Then we'll sculpt in the computer what we want, taking into account existing boundaries and surrounding contours, and see roughly how much dirt we need or will need to remove.) We performed a ghetto survey this weekend with strings and poles, roughly sampling a grid of elevations every 10'. The results we got produced a reasonable-looking representation of the terrain. However, the 10' sampling was a bit too coarse for my tastes and—more importantly—we were only able to sample about 1/4 of the area we need before our measuring poles weren't tall enough or we ran out of reference objects to attach to. So: what hardware exists that I can rent to do this task for me? Ideally I'd like to take a pole, walk all over the property and every few feet stick it on the ground, press a button, hear it go "beep!" and have it record an XYZ sample that is accurate within +/- 6". After doing this for a couple hours I'd plug it into the computer and download some DEM, CSV, or other accessible data of all the points in a non-proprietary format. It looks like Trimble offers some neat gadgets in this arena, and I see that some companies rent this sort of gear out, but not being a landscape survey engineer professional I cannot figure out exactly which, if any, products meet my needs. Can someone who knows this product space please help suggest what I should be looking for? Thanks! _Note: I originally asked this question on Home Improvement, butas suggested there and approved here on meta I'm cross-posting it to this site._


Using Python to combine middle parts of two fields into third field

Below is how I used the Update Cursor in my final script. As you can see getValue was used a few times, perhaps there is a cleaner way to do this? But my concatenated field wouldn't populate properly otherwise.                   sourceRows = arcpy.UpdateCursor(inFeatures)         source = sourceRows.next()               while source:             fieldString1 = source.getValue(fieldName1)             fieldString2 = source.getValue(fieldName2)             concatString= (str(source.getValue(fieldName2))[1:7] + str(source.getValue(fieldName1))[1:7])             source.setValue(fieldName3, concatString)             sourceRows.updateRow(source)             source = sourceRows.next()         del source, sourceRows      I have a script that creates three fields in a feature class. The script adds the X and Y coordinates to the first two fields. I would like data found in the middle parts of the first two fields to populate the third field. Does anyone know how to find the middle part of a field using Python? I see where VBScript seems to have a MID function but I haven't found anything for Python yet. Below is what the third field (fieldName3) should look like.               fieldName1     fieldName2   fieldName3     12345.321       36789.321   67892345     24680.5425      98765.432   87654680      Do you see how I want data to be combined in fieldName3? Below is the script as I have it now.               import arcpy     from arcpy import env     arcpy.env.overwriteOutput = True          env.workspace = "D:/gisfiles/test.gdb/"          inFeatures = "fc"     fieldName1 = "xCoord"     fieldName2 = "yCoord"     fieldName3 = "Combined"     fieldPrecision = 18     fieldScale = 11          expression1 = "float(!SHAPE.CENTROID!.split()[0])"     expression2 = "float(!SHAPE.CENTROID!.split()[1])"          arcpy.AddField_management(inFeatures, fieldName1, "DOUBLE", fieldPrecision, fieldScale)     arcpy.AddField_management(inFeatures, fieldName2, "DOUBLE", fieldPrecision, fieldScale)     arcpy.AddField_management(inFeatures, fieldName3, "DOUBLE", fieldPrecision, fieldScale)          arcpy.CalculateField_management(inFeatures, fieldName1, expression1,                                 "PYTHON")     arcpy.CalculateField_management(inFeatures, fieldName2, expression2,                                 "PYTHON")          sourceRows = arcpy.UpdateCursor(inFeatures)     source = sourceRows.next()          while source:             **#Instead of 12345 I need an expression to combine fieldName1 and fieldName2**             source.setValue(fieldName3, "12345")             sourceRows.updateRow(source)             source = sourceRows.next()     del source, sourceRows


How to make beautiful maps in R?

There are quite some nice possibilities to analyse spatial data in R and in the context of my current project I would like to use R more frequently to do that. Until now I plot my maps with ggplot2 package that brings a lot of practical tools to plot and explore data. Still my maps are not nearly as good looking as the ones I used to plot in ArcGIS. So for publishing purposes I wonder if there are any good tutorials, books, practical tips, packages etc. that I could use to experiment a little bit and make my maps more good looking. Any advice would be warmly welcome :)


How to avoid corrupt geometries when translating between TopoJSON and GeoJSON?

I'm trying to display some geojson vector features on a leaflet map. The size of my GeoJSON is too big (25.4mb of data, 8.47mb gzipped), so I'm trying to use a TopoJSON format to reduce the size without loosing detail. I converted my geojson data with topojson tool and just basic params: topojson --id-property okato_code -p name -o adm.topojson adm.json That significantly reduced the file size (down to 1.54mb of data, 201.3kb gzipped). But when I translate it back ( using topojson.object() ) to GeoJSON on the client, it looks overgeneralized and like it has some incorrect topologies when zooming in. I thought there was no generalization happening during translation to/from topoJSON, then why are the geometries different? How to make them identical? My goal is just to reduce the size of the file without any simplification of geometries. Here's the map - http://basemap.ru/json/ You have to wait for the big GeoJSON to download, it's included for comparison, but after it loads, you can see a distinct difference on high zoom levels. GeoJSON layer is orange and TopoJSON layer is white. I've included the layer switcher to turn topojson and geojson on and off. The data itself is here: http://basemap.ru/json/json/adm.geojson http://basemap.ru/json/json/adm.topojson The script with translation is http://basemap.ru/json/mapscript.js If it makes sense, my map is in Asia Lambert Conformal Conic projection (http://spatialreference.org/ref/esri/102012/)


How do I fix this QGIS sextante error? NameError: name 'sextante' is not defined

A colleague of mine posted this the other day. We haven't got any nibbles yet and I wanted to elaborate on the issue. We are having trouble running the r.mapcalculator function and receive the error below ![enter image description here](http://i.stack.imgur.com/GQ0v9.png) I have gone into the file "HistoryDialog.py" and changed a reference to Sextante on line 101 to 'sextante' (original text shown below). I think I am clutching at straws here, because after this change I receive another error altogether. ![enter image description here](http://i.stack.imgur.com/LaPnY.png) A good googling returned nothing related to the 'sextante is not defined error'. Surely someone else has come across this before. Also for us it is happening on Ubuntu and Windows. We installed the latest dev version through Osgeo4w installer but no luck!!


How to calculate altitude with given lat,long position?

I'm new to GIS. I want to calculate the elevation of any point on earth giving lat, Long (for example like google elevation Api, or the earth tools: http://www.earthtools.org/webservices.htm # cheigit) based the ASTER GDEM MODEL V2 I went to download data from there, and gives me geotif files. Then what should I do? Thanks


What is difference between Geo BI (Geo-business Intelligence) and GIS (Geographic Information System)?

I tried finding information with difference between Geo Business Intelligence and Geo IS but unfortunately no clearly answer being found. What is different between Geo BI and GIS and its advantage/disadvantage?


LinearRing getVertices implementation

I'm fairly new to OpenLayers so please excuse me if this is an ignorant question. I'm curious why the LinearRing getVertices implementation               getVertices: function(nodes) {         return (nodes === true) ? [] : this.components.slice(0, this.components.length-1);      },      does not return the last point since I was under the impression a LinearRing was supposed to represent a closed ring. I was using getVertices to serialize to another format which requires closure and found this interesting. Can someone enlighten me?


How to create shaded relief in ArcGIS?

How can I create shaded relief from contour lines (.shp format)? I see Hillshade in ArcToolbox, but it asks for a raster file. I have ArcGIS 10.0.


GIS research areas in Telecommunications?

I am doing MS in GIS and Remote Sensing; i have a telecommunications Engineering Bachelors degree. Can any one please guide me what sort of GIS research can i do for telecommuncations industry???


Modelling two-way single-track/one-lane width road in ArcGIS Network Analyst?

I don't know how make one-lane roads(only room for one car at a time) for one car waiting for another car to pass by in ArcGIS. Yesterday, I asked about the wait time when the two cars meet in the road, one car should wait for the other car to pass that as How to calculate wait time on a one-lane road with two-way travel? But I realized that I can't even make a one-lane road... How can I make this road in ArcGIS Network analyst (especially, VRP (Vehicle Routing Problem))?


Why does QGIS not detect the CSV/DBF column headers?

I am trying to import a csv or dbf into qGIS so I can perform a join on a shapefile. However, qGIS is not detecting the first row as header names. I know I've done similar joins before and there's never been a problem detecting column headers. Any help would be much appreciated. Thanks.


Displaying 3D subsoil models in ArcScene

I was wondering if anybody had any ideas on how to display a 3D subsoil model in arcscene? The data we will have will be exported from rockworks in point form - and then we need to convert it to - TINs (? i assume). I have never done any extensive work in arcscene before and i couldn't find anything too helpful on the internet!


How to determine which EditEvent is fired BEFORE the feature is created?

I have written an AddIn-Tool and want to handle the OnSketchFinished-Event. The Way               ((IEditEvents_event)_myEditSketch).OnSketchFinished += OnFinishedSketch;      doesn't really work for me, cause I need the get the event before the feature is created. I rather need something like the OnFinishSketch that was fired by the IEditTask. Does anyone have an idea how to solve this problem? I'm very thankful for any suggestions.


Can a raster map in GRASS be saved in JPEG/BMP format?

I'm trying to save a raster map in GRASS in jpeg/bmp format without using the GUI for d.out.file/d.save. Is there a script to do it? I tried r.out.gdal but it produced a map that is black and white.


How can I import ALOHA .pas plume file into QGIS?

I am trying to find a way to import ALOHA's plume model (.PAS) into QGIS? I have not been able to find anything showing how to accomplish this outside of using ESRI-based tools.


Why is st_intersects giving bad results?

I am using PostgreSQL version 9.2.2 with PostGIS 2.0.2 (or just 2.0, not sure), and need to do simple queries to find points that intersect with a polygon. The following test case returns only partial results for an unknown reason. If the minimum longitude value of the polygon is changed from 0 to 1 to 2 to 5 to 20 (leaving the latitudes limits as they are), different numbers of the points are returned, even though all of the points fall into the polygon for all these values.               CREATE TABLE testing.tsttab     (       lat real,       lon real,       tstpoint geography(Point,4326)     );          INSERT INTO testing.tsttab (lat, lon, tstpoint) values (36.1, 22.1, ST_GeographyFromText('POINT(22.1 36.1)'));     INSERT INTO testing.tsttab (lat, lon, tstpoint) values (36.2, 22.2, ST_GeographyFromText('POINT(22.2 36.2)'));     INSERT INTO testing.tsttab (lat, lon, tstpoint) values (36.3, 22.3, ST_GeographyFromText('POINT(22.3 36.3)'));     INSERT INTO testing.tsttab (lat, lon, tstpoint) values (36.4, 22.4, ST_GeographyFromText('POINT(22.4 36.4)'));     INSERT INTO testing.tsttab (lat, lon, tstpoint) values (36.5, 22.5, ST_GeographyFromText('POINT(22.5 36.5)'));     INSERT INTO testing.tsttab (lat, lon, tstpoint) values (36.6, 23.1, ST_GeographyFromText('POINT(23.1 36.6)'));     INSERT INTO testing.tsttab (lat, lon, tstpoint) values (36.7, 23.2, ST_GeographyFromText('POINT(23.2 36.7)'));          select * from testing.tsttab     where st_intersects(     st_geogfromtext('POLYGON((2 35, 2 37, 40 37, 40 35, 2 35))'),tstpoint);      If I convert to GEOMETRY rather than GEOGRAPHY, it performs correctly. Am I mis-using something?


ArcObjects - Create and Draw a Transparent Polygon

I'm trying to create a transparent Polygon in ArcMap; however, it's never transparent. Here's the code I have.               // Create polygon     IPolygon polygon = new PolygonClass();     // create & add points,......          // Create lineSymbol     ISimpleLineSymbol lineSymbol = new SimpleLineSymbolClass();     // set color & style,....          // Create fillSymbol     ISimpleFillSymbol fillSymbol = new SimpleFillSymbolClass();     fillSymbol.Outline = lineSymbol;     fillSymbol.Color = //set to color red//;     fillSymbol.Color.Transparent = 0;     fillSymbol.Style = esriSimpleFillStyle.esriSFSSolid;          // Draw Shape     myMap.DrawShape(polygon, fillSymbol);      RESULT ==> A solid red polygon is drawn. The Transparent setting is ignored. Thoughts?


NASA World Wind Pick From A Group Of Icons

I am using an IconLayer in NASA World Wind (Java SDK) to show positioned items on a map which works fine except for when several are in a group and start to overlay one another. The display of the icons in a group like that is not so much of a problem as I need to allow the end user to easily select the correct icon in that "pile" easily. A workaround would be to allow for a right-click on an icon in that group would generate a context menu that shows a list of the items in that "pile" so they could select the one they want. Creating a context menu on right-click is easy enough but I do not see how I could create a buffer from that point to do a simple geometry relationship like "contains" in NASA World Wind. I have searched the API a bit and have not seen anything that looks like it has potential so I was hoping that someone here might now of either how to do this in NASA World Wind or a simple way of including or achieving this functionality.


Double tap on the map not work in the chrome of Nexus 7

I am stuck with a problem which is the double tap not working in the chrome of Nexus7. Say you have code like below. please help to review it .               <!DOCTYPE html>                   <title></title>         <meta http-equiv="X-UA-Compatible" content="IE=7, IE=9,IE=10" />         <meta name="viewport" content="width=device-width, initial-scale=1.0, maximum-scale=1.0, user-scalable=0">         <link href="http://js.arcgis.com/3.9/js/dojo/dojo/resources/dojo.css" rel="stylesheet">         <link href="http://js.arcgis.com/3.9/js/dojo/dijit/themes/claro/claro.css" rel="stylesheet">         <link href="http://js.arcgis.com/3.9/js/esri/css/esri.css" rel="stylesheet">

                     Could someone please help me out ? **Updated** I test above script in below devices. Please review it.   1. Surfari in IPad2 (Double tap works fine.).   2. IE11 in Surface (Double tap not works) (and Why the feature esri-touch test is false? and the feature esri-pointer test is true.).   3. Chrome in Window7 computer(Mouse double click works).   4. Chrome in Nexus7 (Double tap not works).


Where can I download the last default.style for osm2pgsql? I get `http://svn.openstreetmap.org/applications/utils/export/osm2pgsql/default.style - this directory does not exist` when I try to download it.


problems of display modis ndvi with TIMESAT

I use MODIS MRT tool mosicia tow tiles MODIS NDVI data, then use TIMESAT image view to display the mosaic data, encounting file size waring : ![enter image description here](http://i.stack.imgur.com/Fp4W3.png) and after drawing data, the sequences of two tiles is wrong, original data display sequence is right to left, but display result is left to right. ![enter image description here](http://i.stack.imgur.com/jlMsH.png) but when i use ARCGIS or envi to display the mosaic data, the display sequence is ok. ![enter image description here](http://i.stack.imgur.com/s8FiS.png) so i want to know what is wrong with my mosaic data?


Leaflet: Is there a way to use click event just like mouseout?

I made a map with leaflet by following the official Leaflet tutorial on interactive choropleth maps. I would like to make the map compatible with smartphones, so I added a click event, which acts like hovering on a feature.               function onEachFeature(feature, layer) {     layer.on({         mouseover: highlightFeature,         mouseout: resetHighlight,         click: highlightFeature     });     }      The consequence: When a user is clicking, the `hightlightFeatures` remain highlighted. Is there a way to prevent this? In other words: Is there an event on the lines of `mouseout`, in order to activate `resetHighlight` by clicking another feature?


How to find all intersections and count lines arriving at each intersection?

I'm having a hard time getting off the ground with this, and am trying to move away from ArcGIS / arcpy. I have a shapefile with a network of polylines. I want to find/count all 2-way intersections, 3-way intersections, 4-way, etc. This would involve, I think:   * splitting all lines at intersections   * deleting identical segments (with some spatial tolerance)   * create points at the end of each new line segment   * buffer the line ends, and count the number of points near/at each line end I would like to do this in a python script. Is this overall plan easily implementable in pyQGIS? An outline of the calls I would be using would be very helpful.


How to count the number of fields in a table using arcpy?

arcpy contains a ListFields method, which "Lists the fields in a feature class, shapefile, or table in a specified dataset". The resulting object has a `count` method, but I can't see how to return the field count from this (shouldn't this return an integer?). Aside from iterating through the fields and incrementing a counter, how can I return the field count?               import arcpy     fields = arcpy.ListFields(<table>)     count = fields.count          [Dbg]>>> count     <built-in method count of list object at 0x0EDEABE8>


Images exported from ArcGIS Online as File Geodatabase inaccessible to ArcGIS for Desktop?

I am a new user of ArcGIS Online. I have a system that includes features and images of the objects and I use Esri ArcGIS for Desktop to perform data collecting. Before exporting the whole database as FGDB I use attachment management to download the images as FGDB zip file. How can I reach the images in the exported FGDB file using ArcGIS for Desktop?


node.js GeoLib fails on polygon that straddles the anti-meridian

geolib.isPointInside does not appear to work correctly when the polygon longitudes straddle the anti-meridian. I maintain a system with lots of polygons that do straddle the antimeridian (I'm in Australia :-) so it's important to me. see the test code below that illustrates the problem when it's run in node:               var geolib = require('geolib');          // -------------     // Test 1 This succeeds     // -------------     var polygon = [         {latitude: -46, longitude: 1},         {latitude: -46, longitude: 2},         {latitude: -47, longitude: 2},         {latitude: -47, longitude: 1}     ];     console.log('Test 1: Polygon:',polygon);     console.log('should be true', geolib.isPointInside({ latitude:-46.5, longitude:1.5}, polygon))     console.log('should be false', geolib.isPointInside({ latitude:-47.9, longitude:1.2}, polygon))     console.log('should be true', geolib.isPointInside({ latitude:-46.8, longitude:1.1}, polygon))     console.log('------------');          // -------------     // Test 2 This succeeds     // -------------     var polygon = [         {latitude: 10, longitude: -10},         {latitude: 10, longitude: 10},         {latitude: 20, longitude: 10},         {latitude: 20, longitude: -10}     ];     console.log('Test 2: Polygon:',polygon);     console.log('should be true', geolib.isPointInside({ latitude:15, longitude:-5}, polygon))     console.log('should be false', geolib.isPointInside({ latitude:15, longitude:-15}, polygon))     console.log('should be true', geolib.isPointInside({ latitude:15, longitude:-5}, polygon))          // ----------     // Test 3 : AntiMeridian 1: This fails     // ----------     var polygon = [         {latitude: 10, longitude: 175},         {latitude: 10, longitude: -175},         {latitude: 20, longitude: -175},         {latitude: 20, longitude: 175}     ];     console.log('Test 3: Antimeridian 1: Polygon:',polygon);     console.log('should be true', geolib.isPointInside({ latitude:15, longitude:177}, polygon))     console.log('should be false', geolib.isPointInside({ latitude:15, longitude:-170}, polygon))     console.log('should be true', geolib.isPointInside({ latitude:15, longitude:-176}, polygon))     console.log('------------');          // ----------     // Test 4: Antimeridian 2: This fails     // ----------     var polygon = [         {latitude: 10, longitude: -175},         {latitude: 10, longitude: 175},         {latitude: 20, longitude: 175},         {latitude: 20, longitude: -175}     ];     console.log('Test 4: Antimeridian 2: Polygon:',polygon);     console.log('should be true', geolib.isPointInside({ latitude:15, longitude:177}, polygon))     console.log('should be false', geolib.isPointInside({ latitude:15, longitude:-170}, polygon))     console.log('should be true', geolib.isPointInside({ latitude:15, longitude:-176}, polygon))     console.log('------------');


Hotkey in ArcGIS for Desktop to delete layer from Table of Contents?

I'm using ArcGIS 10.2 for Desktop and wish to remove/delete a layer from the Table of Contents. I cannot find such hotkey. I have tried the `DELETE` button but that does not do it. Also such hotkey is missing in the ArcGIS Help.


Is it possible to import an autocad feature to a feature layer without exporting to .shp/.dwg

Does anyone know if it is possible to import a feature (point, line or polygon) drawn in Autocad directly to a feature layer in ArcMap, so in essence a "feature sucker", rather then exporting the feature to a shp or dwg then importing into ArcMap then copying the feature into the correct feature class. EDIT: The datatype that I want to target is a feature class on ArcSDE (9.2). Versions are ArcGIS 9.3.1 and AutoCAD 2011. Thanks.


Google Maps/Google Earth - KML Overlays and Custom Orientation

I am am trying to get the floor plans for our various buildings imported into 4 separate KML files, one for each floor level so that I can toggle them with javascript. I'm able to overlay our floor plan images and properly scale them in Google earth, but once exported the KML doesn't display properly. Please take note that whatever the orientation of the Google map is, I can't get my overlays to display properly. On the actual API based map, I have it set to be oriented with West at the top instead of North, which seemingly throws us another curve ball, as when I attempt to load the KML per the API's docs and refresh the API based map, I end up with nothing but "There are no images for this zoom level." My KML:               <?xml version="1.0" encoding="UTF-8"?>     <kml xmlns="http://www.opengis.net/kml/2.2" xmlns:gx="http://www.google.com/kml/ext/2.2"     xmlns:kml="http://www.opengis.net/kml/2.2" xmlns:atom="http://www.w3.org/2005/Atom">     <GroundOverlay>     <name>Untitled Image Overlay</name>     <Icon>         <href>http://img694.imageshack.us/img694/2600/walther.png</href>         <viewBoundScale>0.75</viewBoundScale>     </Icon>     <LatLonBox>         <north>37.21924494476235</north>         <south>37.21876115002708</south>         <east>-93.2637540015874</east>         <west>-93.26445522810076</west>         <rotation>87.56327787656277</rotation>     </LatLonBox>      My API v3 Javascript:               <script type="text/javascript">         function initialize() {             var options = {                 center: new google.maps.LatLng(37.22090,-93.26371),                 zoom: 18,                 minZoom: 18,                 maxZoom: 21,                 heading:-90,                 mapTypeId: google.maps.MapTypeId.SATELLITE             };             var map = new google.maps.Map(document.getElementById("map_canvas"), options);                  //here's where it breaks!!!             //it just won't load the kml...or it is and just breaks beyond reason!             var ctaLayer = new google.maps.KmlLayer('test.kml');             ctaLayer.setMap(map);         }       </script>      Here's what it looks like in Google earth (correct): http://img811.imageshack.us/img811/6348/gearth.jpg Here's what it looks like in Google maps (for testing only, using it to load the KML...until I also figure out how to get my API based map to even display the layer without breaking): http://img838.imageshack.us/img838/2237/gmapsr.jpg Does anyone here have any idea what I need to do to get the display working properly, with regular google maps and more importantly my API based google map?


Plotting a circle on a Web map to view the resulting coverage over a central location

Not exactly a GIS question but I hope no one will mind. Does anyone know of a website (something equivalent to Google Map) where you can indicate a location, give a radius and you can view the resulting coverage over that central location ? Thanks.


3D Analyst missing from Tool Boxes in ArcGIS 9.3

I have a missing 3D Analyst from ArcGIS 9.3 (ArcMap - ArcInfo). I need it for a task to be done. How can i install it? Does anybody know it? Should i donwload the extension or is there any other solution?


Looking for an Android Object-relational mapping programming layer, PostGIS oriented

I'm working on a Android app which should manage data on a remote PostGIS database. Before diving into actual programming, I'm searching for ORMs for Android. I've found a couple, namely androrm and ormlite. However, these don't seem to give me an interface to query objects by the geometry column. For example, find all objects within a radius from the current location... I'd need to write SQL for this... What are my best options here? Is there a way to completely avoid PostGIS SQL? Or not really? Thanks in advance. (EDIT, 19/08/11: changed the title)


How to get an array of simple lines from one multipoint linestring in PostGIS?

How can I extract some LINESTRING consisting of 3 or more POINTs to several LINESTRINGs each of 2 POINTs in PostGIS?


Explanations for Various GIS Job Titles

I've been out of school and working in the GIS field (mostly the field of a desk) for over five years now. The one consistency I've found when working for various companies and applying for various jobs is that there is no consistency in the way our job titles correlate to the actual tasks we accomplish on a daily basis. Does anyone have a good list of the general tasks / job description for some of the GIS job titles:   * GIS Technican   * GIS Analyst   * GIS Administrator   * GIS Developer   * GIS Applications Specialist   * GIS Professional   * GIS Manager   * GIS / IT Technican   * GIS Technical Architect I know each title varies by region and company, but a general list would be a great resource. If I have left any out, please feel free to add.


Easiest way to create historical street maps using QGIS

I want to make some historical street maps using QGIS. I understand how to upload a historical OS map from DigiMaps to QGIS, but I want to actually 'draw' the streets myself in order to have some detailed maps. The area I am dealing with was subject to a slum clearance in the 1960s and I am looking to map it how it was in the late-nineteenth and early-twentieth centuries. Could anyone offer advice? Thanks in advance!


Error: already exists.(stg_e_filealreadyexists) when using IMapDocument.Open method?

I got the errror: already exists.(stg_e_filealreadyexists) when using the method IMapDocument.Open() as below codes:               if (File.Exists (strDirectory + "\\temp.mxd"))         File.Delete(strDirectory + "\\temp.mxd");     System.IO.FileStream oFileStream = new System.IO.FileStream(strDirectory + "\\temp.mxd", FileMode.Create);     oFileStream.Write(MyImage, 0, MyImage.Length);     IMapDocument mapDoc = new MapDocumentClass();     mapDoc.Open(strDirectory + "\\temp.mxd"); // Throw exception      Do you know why and how to fix this error ? Thanks and regards, Tai


how to transfer from traditional GIS software developing to internet area

I am a GIS developer, formerly doing some secondary development with GIS software. I now want to obtain some internet-based GIS area in my career. Somebody advised me to first read some networking/linux books. When prepared, then tranfer to some company to provide experience and so on. My question is. Is there any opensource project belonging to this type: server side provide some service, for example, provide some navigation service(or other GIS area service). I want to learn something about these software platforms and how to realize the functionality provided by server to client, Scheme about the whole thing. Ideall6 such type of thing is using in c++/c. Easier for me to learn. My former area and expertise has been using desktop software. Does anyone have any advices? thx.


Is it possible to draw a point interactively and execute an intersection operation (drill down) to identify spatial relationship with other layers?

Is it possible to draw a point interactively and execute an intersection operation (drill down) to identify spatial relationship with other layers using GeoServer? I am a newbie with GeoServer…


Citing background layer from ArcGIS Online that is used in ArcMap

How do you cite a background map layer that is used in ArcMap that is retrieved from ArcGIS Online? For example the background layer is suppose to be credited from the U.S. Geological Survey and in the description it was digitized from the U.S. Geological Survey Professional Paper 1183.


Looking for a list of tools/add-ins for ArcGIS Explorer

I'm hoping to compile a list of add-ins & resources for ArcGIS Explorer. I intend to install it at home, and apparently, there are a number of server feeds out there that can significantly augment the capabilities of this **free** ESRI software (IMS/WMS servers, toolbox items, etc, etc). I thought many of you might be interested in having access to such a compilation as well, especially since it's a free software. I'm quite anxious to see what's out there!! Try your best to give at least a brief description of the add-in/server feed you post in your answers. Thanks everyone


How to sample a 3D raster along a 3D path (or 4D raster along 4D path)?

I have a 3d raster dataset of atmospheric variables (lat x lon x altitude). I would like to sample this raster dataset along a 3d path. Any idea how I can achieve this, preferable with qgis or postgis or rasdaman? **Edit:** Actually, my data is 4d (time x lat x lon x altitude), but having a 3d solution would already be great.


QGIS MergeShapes Error

I am using the MergeShapes plugin with QGIS 1.8.0 on a Windows 7 computer. The issue I am having is when I perform a merge from a directory based on one geometry type, and then switch the geometry type to process a second time (ie: process all polygons in c:\Temp and then process lines in the same directory). This is a limitation when you need to preform multiple merges based on geometry. Any help would be appreciated. Below is the error from MergeShapes. An error has occured while executing Python code: Traceback (most recent call last): File "C:/Users/ryangarnett/.qgis//python/plugins\mergeshapes\mergeshapesdialog.py", line 141, in accept self.progressFiles.setRange( 0, self.inputFiles.count() ) AttributeError: 'NoneType' object has no attribute 'count' Python version: 2.7.2 (default, Jun 12 2011, 15:08:59) [MSC v.1500 32 bit (Intel)] QGIS version: 1.8.0-Lisboa Lisboa, 6416f38 Python path: ['C:/PROGRA~2/Quantum GIS Lisboa/apps/qgis/./python', 'C:/Users/ryangarnett/.qgis//python', 'C:/Users/ryangarnett/.qgis//python/plugins', 'C:/PROGRA~2/Quantum GIS Lisboa/apps/qgis/./python/plugins', 'C:\PROGRA~2\Quantum GIS Lisboa\bin\python27.zip', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\DLLs', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\plat-win', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\lib-tk', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\qgis\bin', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages\PIL', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages\win32', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages\win32\lib', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages\Pythonwin', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages\wx-2.8-msw- unicode', 'C:\Users\ryangarnett\\.qgis\python\plugins\qgis2kml', 'C:\Users\ryangarnett\\.qgis\python\plugins\qgis2kml', '/usr/share/qgis/python', '~/.qgis/python', '.', 'C:\PROGRA~2\Quantum GIS Lisboa\apps\qgis\python\plugins\fTools\tools', '~/.qgis/python', '.']


Proxy Issues for WFS and GeoServer - Running on Local Machine

I see that a number of similar questions have been asked, but I cannot find the answer to mine so I will be as specific as I can. I am trying to create an OpenLayers map that can run without an internet connection. I have created my own OpenStreetMap tiles and am storing them locally on my hard drive. I have also downloaded local copies of the OpenLayers.js and OpenStreetMap.js files, etc. The basic map functionality is all working. I am now trying to plot a number of latitude/longitude points on the map. These points are stored in a PostGIS database, and I have loaded them into a layer in GeoServer. I am attempting to create a WFS layer in my OpenLayers map that will plot the data from the GeoServer layer, but as of right now the points are just not showing up on the map. I understand from previous questions that I probably have to deal with Access-Control-Allow-Origin. I have installed Apache Server, downloaded and modified the proxy.cgi file to allow 'localhost' and 'localhost:8080', and placed the proxy.cgi file in the cgi-bin directory. My question is: Do I need to use Apache Server if I am trying to run this all locally? I have not yet been able to get this to work and I imagine it is because I don't know what I am doing with Apache. I appreciate any suggestions you may have. Thanks, -Rob _**_**EDIT ** _*_** Thanks for the suggestions! iant, I have successfully moved everything into the GeoServer www directory, and I can open my map by double-clicking on the .html file. Strangely enough, I get a 404 NOT_FOUND error when I try to use a url like the one you posted. Unfortunately, my WFS layer still won't show up. In GeoServer, I have the following settings: \- Workspace Name: MyWorkspace \- Workspace URI: myworkspaceuri \- Layer Name: testpoints \- Native SRS: EPSG:4326 I am including my code below:               // Start position for the map (hardcoded here for simplicity)         var lat=42.3583;         var lon=-71.0603;         var zoom=10;              var map; //complex object of type OpenLayers.Map              //Initialise the 'map' object         function init() {                  // allow testing of specific renderers via "?renderer=Canvas", etc             var renderer = OpenLayers.Util.getParameters(window.location.href).renderer;             renderer = (renderer) ? [renderer] : OpenLayers.Layer.Vector.prototype.renderers;                  map = new OpenLayers.Map ("map", {                 controls:[                     new OpenLayers.Control.Navigation(),                     new OpenLayers.Control.PanZoomBar(),                     new OpenLayers.Control.Permalink(),                     new OpenLayers.Control.ScaleLine({geodesic: true}),                     new OpenLayers.Control.Permalink('permalink'),                     new OpenLayers.Control.MousePosition(),                                         new OpenLayers.Control.Attribution()],                      maxExtent: new OpenLayers.Bounds(-8182919, 5035362.5, -7784059.5, 5304535),                 maxResolution: 700,                 numZoomLevels: 18,                 units: 'm',                 projection: new OpenLayers.Projection("EPSG:900913"),                 displayProjection: new OpenLayers.Projection("EPSG:4326"),                      layers: [new OpenLayers.Layer.OSM(                             "Local Tiles",                              "openlayers/tiles/${z}/${x}/${y}.png",                              {numZoomLevels: 18,                              alpha: true,                              isBaseLayer: true,                             //tileOptions: {crossOriginKeyword: null}                             }                         ),                              new OpenLayers.Layer.Vector("Points", {                             strategies: [new OpenLayers.Strategy.BBOX()],                             protocol: new OpenLayers.Protocol.WFS({                                 version: "1.1.0",                                 url: "localhost:8080/geoserver/wfs",                                 srsName: "EPSG:4326",                                 featureNS: "localhost:8080/workspaceuri",                                 featurePrefix: "MyWorkspace",                                 featureType: "testpoints"                             }),                             renderers: renderer                          })                          ]                       } );                           var switcherControl = new OpenLayers.Control.LayerSwitcher();                 map.addControl(switcherControl);                 switcherControl.maximizeControl();                  if( ! map.getCenter() ){                 var lonLat = new OpenLayers.LonLat(lon, lat).transform(new OpenLayers.Projection("EPSG:4326"), map.getProjectionObject());                 map.setCenter (lonLat, zoom);             }         }


How to enter metes and bounds descriptions in ArcMap with only a Basic license?

I am using ArcGIS for Desktop 10.1 Basic without extensions. How can I convert a metes and bounds description such as below? For example, 47 degrees, is there a way you can convert the degrees minutes seconds while using CONTROL G (Direction / Length)? I switched in the Editing Options the direction type to Quadrant Bearing and Direction Units to DMS.               SW, SOUTH 0 DEGREES 02' WEST, 110.51 FEET;      SOUTH 47 DEGREES 52' WEST, 246 FEET;      SOUTH 74 DEGREES 41' WEST, 167.25 FEET;      NORTH 75 DEGREES 19' WEST, 274.85 FEET;      NORTH 41 DEGREES 49' WEST, 339.45 FEET;      NORTH 89 DEGREES 14' WEST, 42.60 FEET;      AND EAST 876.75 FEET TO BEGINNING, CONTAINING APPROXIMATELY 4 ACRES,LESS 4 ACRES IN THE NE CORNER OF SAID SW/4


How to draw a boundary feature of an area?

Which shape file is more appropriate for drawing a boundary of an area. I mean would I take polyline feature or polygon feature for drawing a boundary? I want to digitize a map of metropolitan area. First I have made a manual map and now I am trying to digitize.In this map there are moujas, thanas, rivers, roads, different important location such a market place, police station. So in this case, there is overlapping of borders. Because a thana border can be a district border at the same time. So In this case which feature I chose.


What is a simple way to create a popup table on map

Using arcgis javascript, is there a way to create a popup table after i perform some processing? Basically my user will select an extend and all address that fall within the extend will be shown.


Select empty Postgis (Postgresql) geometry column

I try to select rows in Postgis that have an empty geometry column. This works but gives me the opposite:               SELECT * FROM foo where geom <> '';      And this fails:               SELECT * FROM foo where geom = '';          >>>ERROR: parse error - invalid geometry     >>>SQL state: XX000     >>>Hint: You must specify a valid OGC WKT geometry type such as POINT, LINESTRING or POLYGON      And this returns nothing:               SELECT * FROM foo where ST_IsEmpty(geom);      Does anyone know the correct syntax?


How to remove 'Edit Feature' option in 'View Map' task of ArcGIS Mobile?

We have a Windows Mobile application which is based on ArcGIS mobile's out of the box tasks (version 10.0). We are using 'View Map' task directly. These are the options available to 10.0 ![enter image description here](http://i.stack.imgur.com/9XSfh.png) While upgrading to version 10.1 we have this additional 'Edit Features' which is not desirable for our workflow. ![enter image description here](http://i.stack.imgur.com/nQQya.png) So I want to get rid of this Edit Feature option. Is there a configuration file where I can get rid of it? so technically I am asking to control the behavior of this : {"assemblyQualifiedName":"ESRI.ArcGIS.Mobile.Client.Tasks.ViewMap.ViewMapTask, ESRI.ArcGIS.Mobile.Client","name":"View Map","description":"Browse the map display"}, \- thanks


In cartodb, how can I insert a new row in a table via API?

We've discovered cartodb and we want to use it in our project, but in the API documentation we've only seen the sql.execute() method (docs) in order to insert a now row in a table. Is this the way to insert points dynamically in our map?


Geoserver sql view - parametric errors

I can not get it to work. I started with complex sql and ended with 'hello word' example, but I get errors: My sqls are:   * SELECT t.id, t.user_id, t.coordinates, t.size, t.price FROM my_table t WHERE t.id=%id_param% Error: Failed to create SQL view: ERROR: syntax error at or near ")" Position: 100   * SELECT t.id, t.user_id, t.coordinates, t.size, t.price FROM my_table t WHERE t.id='%id_param%' Error: Failed to create SQL view: ERROR: invalid input syntax for integer: "" Position: 100 I am not sure, what I am missing. I am using approx 2 months old OpenGeo suite (4.0.2). > About GeoServer General information about GeoServer Build Information > Version 2.4-SNAPSHOT Git Revision f295d37d70c3b1f0e51942e8b031253e4d75459 > Build Date 20-Jan-2014 22:29 GeoTools Version 10-SNAPSHOT (rev > 3bb3df4f410750baf9ebc48c9c475a70642aadb6) GeoWebCache Version 1.5.0 (rev > 1.5.x/f981cafd349886bb2a1ef45067f26b1aa8e4aff4)


BigTIFF support in libraries?

What is the level of support for BigTIFF in image processing libraries? GDAL supports this functionality so long as you're building against libtiff4(beta), and I would assume most open source C libraries would as well, since I assume most of them build against either libtiff or GDAL directly. Does ESRI software support bigtiff? Does Java software support bigtiff (presumably through Java Advanced Imaging)?


East North Format

I have a database with a list of address and their location as an East and North coordinate. An example is: East=431833 North=336684 Postcode=DE22 4LF When I try to convert this address using this tool: http://www.engineeringtoolbox.com/utm-latitude-longitude-d_1370.html The result in longitude and latitude in Google gives me a location in Africa, not England where the address actually is. Can anyone tell me which format the East and North coordinates are in and how to convert them?


How to change encoding of layer?

After loading a point layer out of a text file (plug in add delimited text layer) saved in ISO 8859-15 (Latin 9) out of Access, I did save the layer under Qgis 1.8 with the ISO 8859-15 encoding. When going into the properties of the layer, it stays in UTF-8 (layer properties - general - Provider specific options - Encoding). No matter what encoding I set there, after closing and reopening the properties window it still stays in UTF-8. Encoding seems to be still correct within the desktop Qgis, but when opening the project in the Qgis-web-client, not all special characters are correctly displayed. Does anybody know a remedy?


Get Lat & Lng from GeoDjango OSMGeoAdmin

Django OSMGeoAdmin with PointField uses OpenLayers and it's hard to use. I want to add 2 new input elements which are lat and lng. Using jQuery, when values of the input elements are change, I want to change point on the OpenLayer. How can I do that? No complete code or instructions are necessary, main idea is enough. Thanks


What was the rationale for the non-standard UTM zones near Norway?

UTM zones 30V, 31V, 32V, 31X, 33X, 35X, and 37X do not line up with the rest of the UTM zone grid. ![UTM zone map](http://i.stack.imgur.com/GOwfn.jpg) Why was this allowed? According to Wikipedia, UTM was developed by the U.S. Army Corps of Engineers in the 1940s. What prompted them to have special zones for Norway and not for any other country?


With CTM how can I get the latitude and longitude bounding box of a geoPDF?

Embedded in a geoPDF is info on the CTM - Coordinate Transformation Matrix. I am unfamiliar with this CTM. I know how to use python-numpy. How (or is it even possible) for me to use the CTM information to get the latitude and longitude bounding box of the image? Is more information need such as the pixel bounding box of the image? Is there a simple matrix calculation some one can explain to me that I can do in numpy?


Multiple layer line selection based on multiple layer point, then calculate the length of lines

I need to measure water network from a source to multiple meters. Source point is connected to water main line, then water main line is connected to service line, service line is connected to meters. i have 4 shapefiles.......... source point, meters, water main line and service line. I have to calculate distance from source to each meter individually. can anyone help me in VBA (ArcGIS 9.2 or 9.3).


TerraGo geopdf vs. Esri geopdf

I am looking for feedback from people who are familiar with both TerraGo and Esri's respective geopdf functions. I understand TerraGo offers a broader spectrum of functionality, but wonder if TerraGo's enhancements have been worth anyone's while. Any input is welcome; thank you!


How can I export county level data in a .shp file to a text (csv) file in ArcMap?

I have a shape file created as per the answer to my previous question, which includes county level means for each county in the US. How can I now export this data in a .csv file, e.g. with just the state, county, and county_mean fields?


Any issues in converting a AutoCAD 2010 DWG to FeatureClass?

I have a small .NET script that uses ArcObjects such as CADtoFeatureClass to convert DWG files to ESRI FeatureClasses. This standalone program has been working great with DWGs created in previous versions of AutoCAD but I've been told that my code is now not working on DWGs from AutoCAD 2010. As a temporary solution I can have the users save their DWG files in AutoCAD to previous versions but to permanently fix the problem I need to find out more details. So does anyone know or has anyone had a similar issue in trying to convert DWG files to FeatureClasses using these types of ArcObjects on AutoDesk 2010 DWG files? Thanks for any advice. best regards,


Why do my out fields change what map points are displayed on my map?

I have a feature layer within my map control. When I do not list out fields on the feature layer control within the LayerCollection, I see the points, and they appear as expected as I navigate around the map. BUT, when I add a single outfield that exists, the map points no longer appear! I know the layer can display them, because when I filter the layer for a specific value within the outfield attribute, The point appears. when it appears though, navigating around the map will change if the point is displayed. for example, if I zoom in, where the point no longer is in the map, and pan to its position, it no longer is there. When I zoom out or in, it appears, but if the extent does not display the point when zooming, it disappears. What the heck is going on?


Refer to zoom level in TileMill expressions

Is it possible to do something like this:               #ways[highway="trunk"] {       line-color:#c99;       line-width:@zoom/10;     }      Obviously that example doesn't work. Is there another way to do arithmetic expressions involving the zoom level?


How to implement WCS in geoserver?

I have to add raster data in geoserver using Web Coverage Service. I have raster data in geotiff format. The problem is that the geotiff image I have added to the geoserver is using WMS not the WCS facility. Please guide me how to add geotiff image using WCS service. The 3 services of WCS are working properly but when I check the source code in javascript file, it shows the following code for the source of geotiff image ..../geoserver/webmap/wms?LAYERS=webmap%3Ayam2&STYLES=&FORMAT=image%2Fjpeg&SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&SRS=EPSG%3A32643&BBOX=913482.6930837,6340905.95,746408.1930837,6421838.25&WIDTH=542&HEIGHT=330, And in this code the image is using WMS service. Also, it uses GetMap request of WMS. I can't understand how it supports all 3 services of WCS perfectly even though uploaded geotiff image is using WMS service.


GeoTool CSV to Shape Conversion Error

In the Feature Tutorial Example               GeometryFactory geometryFactory = JTSFactoryFinder.getGeometryFactory();      generates error: > method getGeometryFactory in class > org.geotools.geometry.jts.JTSFactoryFinder cannot be applied to given types; > required: org.geotools.factory.Hints found: no arguments reason: actual and > formal argument lists differ in length Is this the right way?               GeometryFactory geometryFactory = JTSFactoryFinder.getGeometryFactory(null);


Select Feature shows outline not filled

public static void SelectBy(this IActiveView view, ILayer2 targetLayer,          string filterClause, esriSelectionResultEnum selectionResult ) {              IFeatureSelection featureSelection = targetLayer as IFeatureSelection;         if ( featureSelection == null ) return;              IQueryFilter queryFilter = new QueryFilterClass();         queryFilter.WhereClause = filterClause;              IColor selectColor = featureSelection.SelectionColor;         selectColor = Local.ColorToIColor( System.Drawing.Color.OrangeRed );         featureSelection.SelectionColor = selectColor;         view.PartialRefresh( esriViewDrawPhase.esriViewGeoSelection, view.FocusMap, null );              featureSelection.SelectFeatures( queryFilter, selectionResult, false );         ISelectionSet selectionSet = featureSelection.SelectionSet;         view.PartialRefresh( esriViewDrawPhase.esriViewGeoSelection, view.FocusMap, null );     }      If I use this code to make a selection calling it with:               view.SelectBy( layer, "[KeyField] = 10000.0",            esriSelectionResultEnum.esriSelectionResultNew);        My feature polygons are always FILLED. If I set the featureSelection.Symbol field to a SimpleLineSymbol (and set SetSelectionSymbol=true) it shows nothing. How do I do a feature selection that is NOT filled (as happens when I use IMap.SelectByShape()?)


Is there any online tool for converting coordinates between different EPSG?

Every time I want to do conversion, I need to get into the console to run proj.4. Is there any online tools for converting coordinates between different EPSG?


Modify the Zoom Buttons in ArcGIS.com JavaScript Application

I'm writing a JavaScript application to view web maps from ArcGIS.com. It is similar to the sample found on Esri's Samples page: http://developers.arcgis.com/javascript/samples/ags_createwebmapid/ I've run the Chrome Developer Tools on that page and it looks like the Zoom buttons in the upper left are autogenerated by the referenced JS libraries. Right now I just want to add a "Full Extent" button, but I may make additional changes later on. Is there a straightforward way to modify those buttons in the JavaScript code? Thanks.


No tables listed in connection to MSSQL 2012 in QGIS

I've (finally) made a connection to my MSSQL 2012 database, but no tables are listed. Here's the tree view of the connection, and the connection properties: ![enter image description here](http://i.stack.imgur.com/wh53Q.png) Some notes and questions:   * I've read that geometry_columns table is NOT required   * Does QGIS indeed work with MSSQL 2012, or is only 2008 supported?   * This database contains both spatial and non-spatial tables   * Using QGIS 2.2 Thanks!


Using QGIS, how can I create attributes within a polygon layer based from a selection in that layer?

I'd like to add a column, or attribute rather, to selected rows within a polygon layer. I would simply like to make a binary system, where those that are selected receive a "1,' let's say, and all others are represented by a "0." How can I do this?


Determine sun elevation without header file?

I have some AWiFS and LISS-III data that I'm trying to convert from at-sensor spectral radiance to top-of-atmosphere reflectance. The problem is, the data didn't come with any header file. All I have is an XML document but it doesn't provide any sun-earth geometry. I could use something like: http://aa.usno.navy.mil/data/docs/AltAz.php knowing the acquisition date, but I can't determine the time of acquisition. I'm using ERDAS Imagine 2013 and the data-sets were obtained from Bhuvan Gateway to Indian Earth Observation at http://bhuvan-noeda.nrsc.gov.in/download/download/download.php Any help or advice would be greatly appreciated. Thanks


Multilateration Code

I am currently writing some code to do multilateration using TDOA and currently struggling and wondering if anyone has any reference code that can be used as a test case with some sample data that I have! At the moment I am unsure if it is the code I have written or whether it is the geometry of the points that is off. One example group of packets that I have received is               station1 = [50.919149572907784,-1.005043823556842,90.729];     station2 = [50.8304,-1.21956,0];     station3 = [51.0569,-1.26224,100];          stations = [station2,station1,station3];     timestamps = [28704.175275109, 28704.175334203, 28704.175334875];      My code currently locates this to 50.910447339754654,-1.2347406519401736 which I am 99% sure is wrong! Does anyone have any code that can verify this at all? One bit of code I am testing against and get very similar results is https://github.com/bistromath/gr-air-modes/blob/master/python/mlat.py


Short-term utility of ArcGIS/GIS certifications?

(Correct me if this belongs in meta) For undergraduate students considering a career using GIS but who will not be graduating with a degree in it, is it a worthwhile investment of time to pursue a GIS certificate/certification? If so, what criteria are important to consider when selecting a test/program?


repeat kml on wrapped dateline layer

I'm using openlayers3 with a based map and a kml over it as a vector layer, i would like to see kml panning outside the dateline. I tried the code below but i see kml only one time. I used the code below               var vector = new ol.layer.Vector({     isBaseLayer       : true,         wrapDateLine      : true,     sphericalMercator : true,     wrapX             : true,     source: new ol.source.KML({      projection: 'EPSG:3857',      url: 'data/kml/2012-02-10.kml',     })     });          var map = new ol.Map({      layers: [raster, vector],      target: document.getElementById('map'),       view: new ol.View2D({       center: [876970.8463461736, 5859807.853963373],       zoom: 10      })     });      Do you have any suggestion for this issue? thanks


feature label not work in xhtml file

Good Morning: I have a problem, I develop an jsf application with openlayers api, this example http://openlayers.org/dev/examples/vector-features-with- text.html work fine with html page, but when I change to xhtml page the label text doesn´t showing. I hope anyone can help me, Thanks for your time and help. Here is my code.               <?xml version="1.0" encoding="UTF-8"?>     <!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">               <meta http-equiv="Content-Type" content="text/html; charset=utf-8"></meta>     <meta name="viewport" content="width=device-width, initial-scale=1.0, maximum-scale=1.0, user-scalable=0"></meta>     <meta name="apple-mobile-web-app-capable" content="yes"></meta>     <title>OpenLayers Labeled Features Example</title>         <link rel="stylesheet" href="../openLayers/theme/default/style.css"         type="text/css"></link>     <link rel="stylesheet" href="../openLayers/style.css" type="text/css"></link>
OpenLayers Labeled features example

         vector, feature, labeling, symbolizer, light

         Label vector features with a text symbolizer.

This example shows drawing simple vector features with a label.


What's the difference between "Copy Features" and "Feature Class to Feature class"

I have a feature class within a feature dataset in a FGDB, and I need a copy of it in another FGDB. In ArcGIS 10.1, I've used the "Copy Features" tool, and it looks to be the same (according to the "Feature Compare" tool). Using the "Feature Class to Feature Class" tool, I receive the following messages from "Feature Compare":               Field: Field OBJECTID domain fixed is different (Base: false, Test: true).     Field: Field Shape domain fixed is different (Base: false, Test: true).      Does anyone have any insight on the differences between "Copy Features" and "Feature Class to Feature Class"? What do the messages from Feature Compare mean?


Featureclass search failing in C++. But working in C#

I need to perform search on a FeatureClass using ISpatialFilter for a given Point (x,y). I wrote a program in C# to do that and its retrieving the features as expected. But the similar program in C++ is failing for some feature classes and working for some. I am wondering where it went wrong. The below c++ code works fine for featureclass _cbsa_. But doesnt work for _BlockGroups_.                   //Create Point         IPointPtr ipPoint = NULL;         ipPoint.CreateInstance(CLSID_Point);         ipPoint->PutCoords( dX, xY);              //Create Envelope         IEnvelopePtr envelope(CLSID_Envelope);         hr = ipPoint->get_Envelope(&envelope);              //Open Feature Class         IWorkspaceFactoryPtr pWorkspacefactory(CLSID_FileGDBWorkspaceFactory);         IFeatureWorkspacePtr featureWorkspace;         IWorkspacePtr workspace;         hr = pWorkspacefactory->OpenFromFile(CComBSTR(gdbpath,0,&workspace);         featureWorkspace = workspace;         IFeatureClassPtr featureClass;         hr = featureWorkspace->OpenFeatureClass(CComBSTR("CensusBGroups"),&featureClass);              // Set ISpatialFilter Properties         BSTR shapeFieldName;         hr = featureClass->get_ShapeFieldName(&shapeFieldName);              ISpatialFilterPtr spatialFilter(CLSID_SpatialFilter);         spatialFilter->putref_Geometry(envelope);         spatialFilter->put_SpatialRel(esriSpatialRelIntersects);         spatialFilter->put_GeometryField(shapeFieldName);              // Search         IFeatureCursorPtr featureCursor;         featureClass->Search(spatialFilter, false, &featureCursor); // featurecursor doesnt contain any featueres         IFeaturePtr feature;              // Get First Feature         hr = featureCursor->NextFeature(&feature); //feature is null         IFieldsPtr totalfields;         hr = feature->get_Fields(&totalfields);


Changing projection in KML from Google Earth to Maps

I have created a bunch of map tiles for Google Maps mistakenly using the Google Earth projection. Is there a tool available to batch warp the KML files to the correct projection or will I need to generate the tiles again from scratch?


List items are replaced with new reference when adding

Based on the query i am getting some **features** in **pFeatureCursor**. I am adding that **features** into a **List(IFeature)**. When iam iterating the **pFeatureCursor** and adding the values into **pFeatureList** all values in the list are replacing with the new reference **pFeature** Please refer my code below               //List is declared globally     List<IFeature> pFeatureList = new List<IFeature>();          //In a method i am trying to add features in a list     //based on queryfilter     pFeatureCursor = pFeatureClass.Search(pQueryFilter, true);     IFeature pFeature = pFeatureCursor.NextFeature();     while (pFeature != null)     {       pFeatureList.Add(pFeature);       //lstQuerValue is the listbox       lstQueryValue.Items.Add(pFeature.get_Value(1).ToString());            //When below line is execute whatever the recent pFeature        //value is replacing all the items and duplicating in pFeatureList       pFeature = pFeatureCursor.NextFeature();      }      Please give some idea how to avoid it.


How stable is MapWindow on a 64 bit machine

One of my favorite ArcGIS extensions (N_SPECT) is no longer going to be supportd under ArcGIS 10. It will be totally rewritten as a MapWindow plug-in, so I need to get up-to-speed with that in case a project comes along where I need to use the tool. I have a 64 bit Windows 7 machine, but the MapWindow website lists the 64 bit version as Experimental. Has anyone here used the 64 bit version? Have you found it to be less stable than the 32 bit version? * * * To clarify - I'm looking at the desktop application and use of associated plug-ins. I'm an end-user, and will not be working on or developing anything with the MapWindow Active X Control. The N-SPECT tool was process intensive as an ArcGIS extension, and I assume it will be in MapWindow, too. That's why I'm looking at the 64 bit version. I suppose all of this may be moot depending on how NOAA develops their new N-SPECT tool?


Can I use LASTools to classify point data?

I am presently using LAStools to analyse LiDAR data. Can LAStools classify the points into various features such as vegetation, roads, buildings etc? If so, how is it done?


How to edit Python Add-In and get all changes to work?

I have created an Add-In using Python. However, when I edit the Python script for that Add-In and re-create the '.esriaddin' file, I don't get updated results when using that Add-In. Do I need to re-create the Add-In each time when I edit my Python script to cater to changes?


shapefile batch clip and merge combined

I have two shapefiles, one is polyline, the other polygon. For each polygon I would like to clip the polyline vectors inside. Using arcpy I already found part of a possible solution: http://gis.stackexchange.com/a/19934/6458               import arcpy          arcpy.env.workspace = r'D:\Projects\GDBs\slowbutter.gdb\IPAS'     rows = arcpy.SearchCursor('HspAOI')     for row in rows:         feat = row.Shape         arcpy.Clip_analysis('HspWBD_HU12', feat, 'HspWBD_HU12_' + str(row.getValue('NAME')), '')      This however would create a new shapefile per polygon. I have 8000 polygons and would have to merge these using arcpy.Merge_management(). Is there maybe a way to not have to merge 8000 separate shapefiles? I work with arcpy and ogr so the latter might also be an option?


Dissolved polygon shows 'ghost' of previous 'child' boundaries

I dissolved the polygons specifically to get a unified outline without any interior boundaries showing. Now I've got a partially broken interior line within the polygon. This is for cartographic display reasons; I know that I can display the polygon without a boundary, but I'm already using that technique to display another layer in the map. I used the dissolve tool in ArcMap 10 (ArcInfo) SP4. **Why does this tool leave artifacts of previous polygon boundaries??** ![Polygon artifacts](http://i.stack.imgur.com/hfvmO.png)


Alternate of Near Distance Tool of ArcGIS 10.1

I have basic license level of ArcGIS 10.1, and near distance tool is not available with it. Can anyone please suggest of another similar purpose tool, under basic license, to get distance of one polygon from another polygon or any workaround? Thanks.


How to generate Digital Surface Model (DSM) from aerial image?

Can I use Erdas LPS 2011 for automated DSM generation with low resolution (5m to 15m) by one aerial image? I don't have a stereo pair images for the city, thanks.


Atmospheric Calibration for Landsat 7 ETM+ Gap Filling

This is for thesis work involving mapping of above-ground mangrove forest biomass using Landsat 7 ETM+ imagery. I am assessing correlations between biomass with vegetative indices, e.g. NDVI, and then mapping biomass. I am using a 04/02/2013 Landsat 7 ETM+ image with Scan Line Corrector off, so I need to gap fill. I am using a 08/03/2013 image for gap fill. Both images have been atmospherically calibrated before gap fill using DOS_improve. However, the reflectance of the two images still contrast to each other. I do not know why the reflectance is contrasting, or how to better match the reflectance of the 08/03/2013 image to the 04/02/2013 image.


Why doesn't SQL sum(dollar) work in SQLITE?

I've got a sqlite table with real numbers, two decimal place precision, in a mydollar_column. (There are no null values, but some of the values are negative numbers) For some reason, using sql on those values, always returns 0.00000. I've tried this in Qspatialite (with QGIS 1.8) and Spatialite_GUI. Here's an SQL example: ` SELECT sum(mydollar_column) FROM myTable ` If I use the sum function on an integer column, it works. For example, I used `sum(pkuid)`on the primary key, and the sum() works as expected. As a test, when I export the table to csv file, the values of my dollar_column are real numbers with two decimal place precision, for example "1234.56" But Qspatialite and Spatialite GUI always shows the values as, for example "$1,234.56" I've also used the total(mydollar_column) function with the same unsuccessful results.   1. Does anyone know why I'm unable to calculate my dollar_column?    2. How does Qspatialite and Spatialite_GUI interpret the values of mydollar_column as dollars, when sqlite is typeless columns?


Google maps android vector rendering : how do they make it that fast?

Most of the question is in the title. If you are on android (I believe its the same on IOS) and use maps you will notice the quality of the maps app and how the vector maps are drawn in a very fluid manner. From what I understand, they are delivering vector tiles wich are encoded using protocol buffers which is pretty easily achievable on my side. However I was wondering how they managed to achieve real-time rendering of vector tiles in such a fluid manner. I am currently using mapsforge to render vector based maps but the rendering times of tiles are pretty far from what I can see in google maps. Any ideas on how I can approach such awesomness ? What kind of voodoo tricks are they using ? Thanks for any insight !


When will QGIS version 2 be available?

Just wondering, how close is the public release of version 2? I do realise that I can download 1.9, but was wondering about the public version.


How to choose between Markers and Features for drawing Openlayers objects?

Yesterday I was trying to create first mock-ups of my GPS tracker interface, and was exploring the available OpenLayers objects for drawing on a map. I began with simple task - placing a marker on a map, on a Markers layer. I added the appropriate image to the marker, learned how to move it on a map, and added a Popup to the onclick event that opens additional info about that marker. Works great. Since I wanted this marker to represent the vehicle, the next step was going to be to populate the content of the marker's Popup using AJAX. When the onclick event occurs I'll send the marker's ID to my python script, fetch the latest data from my postGIS database and update the marker and it's popup appropriately. But then I saw in the OpenLayers documentation that markers don't have ID attribute, and that I'll probably have to end up using vector or feature object. OK, both of them have id attribute and createMarker methods, so I'll easily transfer my Markers into Features/Vectors. After that I saw an example on the web where someone just initializes id attribute to the marker like this:               marker = new OpenLayers.Marker(new OpenLayers.LonLat(90,10),icon.clone());     marker.id = 353;      And it works fine, I can now add various attributes to my markers and use the as needed. For example:               marker.id = 353;     marker.popups_id = 382;     marker.data = "Marker for a vehicle";      So, after I've done what I wanted in the start I started wondering - what are the main reasons for these 3 different Classes if I can take Marker and fill it with various data to make it look like a feature. Am I doing it wrong, is there a better approach? It seems to me that I could create my map objects using Markers only, but I don't want to regret that later if I hit it's limits. Are there any size/performance limits and recommendations when to use Markers and when Features or is it all based on testing and refactoring as needed? Just to add, my application will be a simple GPS tracker. So I'll be using OpenLayers to display objects, move them, animate them, show their properties and similar tasks. I'll have to support having minimal 300 active objects on a map in the same moment, and refresh their coordinates at least every 15seconds. TIA for all your advices that could possibly lead me in the right direction of creating the effective interface for my application.


Where to find the variables of QgsPolygon and QgsPolyline?

now I want to use QgsPolygon (and Qgspolyline,too) class in python I have a example by using QgsPoint but I don't know how to use in QgsPolygon following is my Qgspoint example                   feature = QgsFeature()         qpoint = QgsPoint( float(x), float(y) ) #get qpoint's attitube here         feature.setGeometry(QgsGeometry.fromPoint(qpoint))#set point's geometry here         if numpy.isnan(val):             pass #no assign value means assign NULL value         else:             attr_val = round(val,4)             feature.addAttribute(0,QVariant(attr_val))         feature_list.append(feature)#append value here      I want to change QgsPolygon(or QgsPolyline) instead of QgsPoint but I couldn't find the variables of QgsPolygon(or QgsPolyline) where to find the variables of QgsPolygon and QgsPolyline? 2013/06/27 edit: I have another problem here ,                 for s in dataset:             x=s[0]             y=s[1]                   feature = QgsFeature()             qpoint = QgsPoint( float(x) , float(y) )             polygon=[]             polygon.append(qpoint)                  feature.setGeometry(QgsGeometry.fromPolygon([polygon]))      but it's still get                       TypeError: QgsFeature.setGeometry(): argument 1 has unexpeced type 'NoneType'      how to fix it? Have I mistake the meaning of your answer?


How to fix strange rendering behaviour of OSM Data in Mapnik?

My current setup so far is a running Apache webserver with tilecache.cgi script and a mapnik rendering engine connecting to my postgis database containing my OSM data. The problem is very easy to describe. All roadnames, roadnumbers, descriptions are to close to each other. I set up tilecache with MetaTile=true, TileSize=8x8 and it doesn't change. What am I doing wrong? Thanks! ![Screeny](http://i.stack.imgur.com/xbXyP.png)


What would be a good symbol for the map legend itself?

I am reformatting a mobile web map. Many of our users are equipped with mobile phones rather than large tablets, so screen real estate is at a premium. I've decided to remove the legend from default view for these devices. Users will have to click on or hover over a small control button to view the legend. So naturally, I now need a symbol that conveys "legend" to place on the control button. Leading candidates:   * A key: intuitive, but it says "key" more than "legend". A difference many users might not care about, but I do...   * A question mark: it's engaging and implies it can answer questions about the map.   * ~~Will Smith~~   * The word "legend". This is a last resort. The other control buttons have icons, so this would kind of throw a wrench in the works, but it definitely conveys the concept of a legend. I'm not completely satisfied with any of these, which is why I've come to you. What says legend without saying legend?


Mobile, Offline viewer for geographic data

Currently, field crews consumer geographic data via laptop (no internet, no gps). ArcReader. I would like to get them shapefiles (or some derivative of a shapefile) on a mobile platform (smartphone) - because of its usability (built in GPS) and to get them acclimatized to using (trusting) a smart phone, among many other reasons. It needs to be offline (we have limited cell service). It needs to contain some basemap information (which I can providen if needed; imagery, road centerlines, etc...). No editing, just viewing. What are my options here? Is there an app I can download that provides this functionality? Very surprised ESRI doesn't have an offline viewer for data on smartphones or tablets.


How to access a table in a personal geodatabase using C#

I've a personal geodatabase which has some shape files included. To get a "complex" relation I've created a single table which will be have 1:1 or 1:N references which the different shape files (atm. 3). Now I want to know how I can:   1. access the table in the personal geodatabase,   2. add new entries to this table,   3. remove/modify existing entries using C#. Does anybody done this before and can help me to solve this question? Thanks a lot!


How to configure web map data layer that can accept csv data uploaded by user?

I made this map: http://tinyurl.com/a368268 on arcgis online. How do I get the same map but instead you can upload a new csv every time? Instead of having to first add it as data layer on arcgis online and then configure pop ups after doing that. I basically want what I made but be able to upload a new csv with new data but in the same format. The fields I'm using currently are: ImageName,Latitude,Longitude,PicURL


Python package for Google Maps API v3?

I am looking for a Python package that can generate the JavaScript to display Google Maps API v3. I would like to integrate such package into a django website. I know pymaps exists but it's for Google Maps API v2 only. Do any other solutions exist?


pgRouting - Route between multiple destinations - in order

I want to calculate the route between multiple points, respecting the order of the list. For example, if I have the following list:   1. Lat: -40, Long: -30   2. Lat: -41, Long: -31   3. Lat: -40.5, Long: -30.5 I want the route to go from point 1 to 2, and then from point 2 to 3.   * I **don't** want to calculate the best path from 1 to 2 and 3 (Travelling Sales Person) because that would give me a path 1-3-2, as the 1st point is closer to the 3rd point than the 2nd one.   * I **don't** want to calculate the best path from 1 to 2, and from 1 to 3 (Mutliple destination Shortest Path Dijkstra), because that would give me a path from 1 to 2 and from 1 to 3. **I want** it to calculate the best path from 1 to 2, and then from 2 to 3. All with just 1 query, using pgrouting. It would be something like a TSRP but preserving the order of the points. Google already does that (https://developers.google.com/maps/documentation/javascript/directions#Waypoints) and if you want, they can optimize the order of the waypoints (using TSRP)