Unable to see the featureLayer but it is present on the map

I have created a `featureLayer` from `featureCollection` and I am able to add it to the map and when I execute, `map.graphicsLayerIds`, it also returns the id of the featureLayer, `["resultLayer"]`. But I am not able to see it on the map. I have also set the renderer but still it makes no difference. This is the application. It takes some time to execute (around a minute). This is the code I am using to create the `featureLayer`,               function addResults(results) {             //console.log(results);             var features = results.value.features;                  var layerDefinition = {                 "geometryType": "esriGeometryPolygon",                 "fields": [{                     "name": "FID",                     "alias": "ID",                     "type": "esriFieldTypeOID"                 }, {                     "name": "STATE_NAME",                     "alias": "State Name: ",                     "type": "esriFieldTypeString"                 }, {                     "name": "NAME",                     "alias": "County Name: ",                     "type": "esriFieldTypeString"                 }, {                     "name": "value",                     "type": "esriFieldTypeDouble",                     "alias": "Absolute Value"                 }, {                     "name": "RESULT2",                     "type": "esriFieldTypeDouble",                     "alias": "Percentage Change"                 }]     }             var featureCollection = {               "layerDefinition": layerDefinition,               "featureSet": {                 "features": features               }             };             var states_Symbol = new esri.symbol.SimpleFillSymbol(esri.symbol.SimpleFillSymbol.STYLE_SOLID, //EVENTS                         new esri.symbol.SimpleLineSymbol(esri.symbol.SimpleLineSymbol.STYLE_SOLID,  //EVENTS                         new dojo.Color([255,255,255,0.35]), 1),  //EVENTS                         new dojo.Color([125,125,125,0.35]));                  var popupTemplate = new esri.dijit.PopupTemplate({                     title: "{NAME}",                     fieldInfos: [{                         fieldName: "STATE_NAME",                         visible: true,                         label: "State Name: "                     },                     {                         fieldName: "value",                         visible: true,                         label: "Absolute Change: $"                     },                     {                         fieldName: "RESULT2",                         visible: true,                         label: "Percentage Change: "                     }                     ]                 });                  resultLyr = new esri.layers.FeatureLayer(featureCollection, {               id: 'resultLayer',               infoTemplate: popupTemplate             });             resultLyr.setRenderer(new esri.renderer.SimpleRenderer(states_Symbol));             map.addLayers([resultLyr]);                  dojo.connect(resultLyr,"onClick",function(evt){                map.infoWindow.setFeatures([evt.graphic]);             });                   }      Any idea what I am missing?


How do you create a button in ArcGIS and have it run code?

I'm brand new to ArcGIS and have only a year or two of programming experience. I know Java and a little C#. I have a Java app that processes geographical data. I need to have a button in ArcMap that, when you click it, sends map data to this Java app. This task was already implemented in VBA; it's my job to translate it to something else since VBA is no longer supported in ArcGIS. Any thoughts on how to accomplish this? I'm not sure where to begin. I'm guessing either using the embedded Python console or the ArcObjects SDK for Java or .NET. Just want to see what someone with experience would recommend. Thank you.


Where can I find high resolution data on residential land use?

I am currently searching for high resolution data on residential or commercial land use in the USA. I'd like to use this data to determine for a given location wether this location is in a house with people living in it or not. The data doesn't have to be perfect, as I only intend to do a case study on that. Therefore, it would be sufficient to find such data for a bigger city (e.g. Boston). Does anyone of you guys know if such data exists? As I am not planning to use the data commercially, I can't afford paying very much for it. With this restriction, I can't use business directories like InfoUSA or Dunn & Bradstreet, which are mention in this question. I considered using OSM, as it contains a node attribute `landuse` that can be used to mark a node as used for residential use. Unfortunately, the data in OSM currently doesn't seem to be very accurate in that respect. If I am wrong there, I would be very grateful for corrections and hints.


Why do ArcGIS Desktop 10 HTML popups give JavaScript/HTML Error?

Using ArcMap 10 Desktop HTML Popups. The issue I am having is that the html < select > elements aren't working - you can select it and use the arrow keys to navigate through(up/down) the list but you can't bring up a drop down with the visible list bit. When on desktop IE running in ie 7 or 8 mode selects don't work I get a popup notification but it wont allow me to unblock it. This a bug, is there a fix? can I upgrade the IE higher than 7/8? chromeframe? I just need to get the selects working.               < ? php          #$tmp = "< select id='something' >";          for(something){          $value = "value";          $name = "name";           #$tmp.="< option value="$value" >$name< / option >";                  }                 #$tmp.="< /select>";          echo $tmp;          ? >


How to modify a Marker object's icon after it's been drawn?

I'm displaying gps track data (e.g. points in time, like the path a vehicle covered) and I need to allow the user to "play back" the track, by having the markers change color one at a time in sequence. So imagine a sequence of yellow markers and then during playback I change the first one to a red marker, then the next one, etc. I've done a prototype using point features because it seemed easier at the time to implement this with Points because I could do the following (pseudo code):               p = getpointFromLayer();  // get the yellow point from the layer     p.style = myRedStyle; // change it's style     p.layer.drawFeature(p); // redraw the feature using the red style      Now I'd like to re-do this functionality using Markers because Markers seem to generally be more flexible for other functionality, but I don't know whether it's possible to somehow change a specific marker's appearance. Any ideas?


Is there a way to configure user logins to be case-insensitive in WMX (JTX) for ArcGIS 10?

I have some users who login on different workstations and the case of the username is different: John.Doe might be john.doe on another workstation. I would like WMX to treat this as the same user. I figure this might addresses by a windows configuration change, but I'm hoping I can leave the OS alone and just have WMX ignore case.


I am not able to display Ground Overlays in a KML file with OpenLayers. Does OpenLayers supports KML Ground Overlays? If not, how would I proceed in writing my own code or plugin for OpenLayers so that I can display Ground Overlays?


Calculate midpoint from a series of latitude and longitude coordinates

I've got a series of longitude and latitude coordinates which represent a building outline e.g.               -0.5485381346101759,53.2285150736142     -0.5482220594232723,53.22842450827133     -0.5482298619861881,53.22841205254449      ...(intermediate points not listed)...               -0.5483123769301657,53.22882101914848      How can I work out the midpoint? I've found tutorials that show how to do it if you've got three coordinates (e.g. http://mathforum.org/library/drmath/view/68373.html), but in many cases I've got more than three. Thank you


Calculate distance of points along a route QGIS

I have a line layer and a point layer, I need the **distance between the first and the last point measured along the line**. How can I do that using QGIS? Do I need to know linear reference? I don't have any knowledge on PostGIS.


QGIS 1.6.0 problem during printing & export to PDF

When trying to print or export my map to pdf something happens to my document. When trying to print, the outcome always results in the printer producing a solid grey rectangle, and when I try to export to pdf the file becomes damaged and will not open. I'm trying to print a 36x24 in map if the size makes any difference. Has anyone had a similar problem?


Learning resources for beginning differential topology for a programmer?

With this question I'm looking for suggestions about how and where to start learning differential topology starting from a good base of linear algebra and geometry, also consider that this know-how is supposed to be used in a programming context and for reading and writing algorithms.


Will hundreds of attributes cause performance problems in an ArcGIS Server layer?

This relates to ArcGIS Server 10.0 SP1, with data served from a file geodatabase polygon featureclass. 1) For performance reasons, is it a bad idea to have a featureclass with hundreds of attributes? 2) Is it a bad idea to create an ArcGIS Server map service (containing hundreds of attributes) from this featureclass? 3) Is there a performance implication if I build a featureLayer in the JS API, and only specify the attributes I currently need? eg, the map service might contain 500 attributes, but I define my feature layer using featureLayer.fields = [x, y, z] Thanks for any advice, and please let me know if you need any further details.


How to tell if a country has states and what the state codes are

With the help of this community I understood how to find country codes and query geoNames data based on that. What I am stuck on now is understanding how to break things down by state where needed. For example, take Canada. There are a few big Provinces like Ontario, Alberta. I just happen to sort of know that. But how can I get to know such information about other countries like Mexico or Russia, etc. Thanks, Alex


Mapserver - Hide Labels

I'm creating a map with some layers, each one has a label, the problem is that, Mapserver is displaying the labels of bottom layers, so it is mixing all the label. Is it possible to hide label that is under others layers?


How to create Where Clause from user input via Python list?

I am trying to perform a query based on values in a Python list. The where clause will then be used in an Update Cursor, so that the rows NOT equal to the values in the list will be deleted. This is what the list will look like after it is passed into the script from the user (using a multivalue parameter in the script tool), and the semicolon delimited string has been split. This is just an example of many choices the user can make. userSelectionList = ["40","77","95] selectionQRY = "RTE_1_NBR <> 'userSelectionList'" <-- I know this doesn't work Hard coded, it would look like: selectionQRY = "RTE_1_NBR <> '40' AND RTE_1_NBR <> '77' AND RTE_1_NBR <> '95'" But the input and number of selection depends on the user's selection. It should finally be something like this:               userSelectionList = ["40","77","95]     selectionQRY = "RTE_1_NBR <> 'userSelectionList'"          with arcpy.da.UpdateCursor(userI,["RTE_1_NBR"],selectionQRY) as userUC:         for row in userUC:             userUC.deleteRow()     del userUC      Is it possible to make the where clause from a user input/dynamic Python list? Thank you!


ogr2ogr giving "ERROR 1: Failed connect to 192.168.4.18:8181; Operation timed out" error?

This is my command, which was working OK yesterday:               ogr2ogr -f "GeoJSON" ../geojson/myfile.geojson ../data/gml/myfile.gml        -s_srs EPSG:27700 -t_srs EPSG:4326      The output is this:               ERROR 1: Failed connect to 192.168.4.18:8181; Operation timed out      I'm baffled, why would ogr2ogr need to connect anywhere?


Remove unneeded columns from geographic data file?

I have some GML data that has some unneeded columns in it. I'd like to remove these columns in order to reduce the file size. (To clarify: I mean that the file has many points in it, and each of these points has some metadata attached, and I don't need all of the metadata.) I'd also like to do this from the command line, as I have many files. I can convert the GML files into different formats using `ogr2ogr`, so my question is this: what's the easiest way to remove a column from a geographic data file, from the command line? I'm a Python programmer so my instinct would be to convert the file to GeoJSON, then use Python's JSON-parsing methods to read and rewrite the JSON file. But perhaps there's a more direct way, using a GDAL tool or similar.


Finding a common area by a set of ID variables?

I am trying to find a way to find a way to get matches of areas based on multiple sets of ID variables. I am using state-county and state-puma as an example below, but I am wishing to do this for more complex levels. Basically, these areas are divided up in different ways: one by county, and another by public use microdata area (PUMA). I am able to know what county and puma each area is in, but this only gives me parts of a county or parts of a puma. I want to find the aggregated area which represents only whole counties and whole pumas, which leads to me to believe this is a special type of clustering. For instance, in this table, there are four main types of matches. The first type represents a whole county and a whole puma, hence there is only one record for each. The second type represents five counties within one puma, so the puma is divided into five parts, hence five records. The third type is similar to the second type except one county is divided into three whole pumas, hence three records. The last type is a real pain, parts of counties are in some cases represented by parts of pumas, and vice versa.               ST_CO   ST_PUMA     39003   3902500          39039   3900100     39069   3900100     39125   3900100     39171   3900100          39017   3905401     39017   3905402     39017   3905403          39015   3905700     39025   3905600     39025   3905700      In essence, I am trying to get a result similar to the one shown below without doing it by hand or eyeballing it. For instance, Pair number 4 is linked together because county 39025 is made up of two different pumas (3905700 & 3905600), yet one of 3905700 is partially linked to county 39015 as well.               ST_CO   ST_PUMA PAIR     39003   3902500 1          39039   3900100 2     39069   3900100 2     39125   3900100 2     39171   3900100 2          39017   3905401 3     39017   3905402 3     39017   3905403 3          39015   3905700 4     39025   3905600 4     39025   3905700 4      Any tips or help would be appreciated as I am currently drawing a blank for ways to routinize this operation. Additional Note: Based on the advice to list what software to use I put down ArcGIS. However, I am open to trying many platforms that I am reasonably comfortable with (even though I am still a novice in some platforms). For instance, if there was a way to solve this problem in R, SAS, Stata, Excel, Python, SPSS, or QGIS, I am open to suggestions.


Extract multi values to point error

I am trying to extract values from my rasters : slope, aspect, landcover etc. for my animals group locations. I used Spatial analyst>Extract multi values to points but when I look in my attribute table after the process the results are the same for all the points and make no sence : e.g. 0 for a landcover, -3x10-38 for a slope etc. Do you have an idea what I am doing wrong? I am working with ArcGIS 10 Thank you


ArcPad 10: Accessing the Geography information of a feature through ArcPad Studio (VB Script)

I am trying to access the geography feature property of a feature being edited. I need to do it programatically (in order to auto populate a field using the geography) using VB Script. I cannot seem to figure out what the ArcPad object is that I want to access in the script. Any suggestions? Thanks EDIT How do I access the Points Collection of an added feature in ArcPad studio?


How to create histograms

I am very new at using GIS data and in particular QGIS (1.7.3) and I need some help to produce histograms regarding my data. I have a great number of data about an algae density collected in 8 sampling sites 4 times a year for 5 years. I would like to analyse them, underlining the differences in density, for each site, both among times and years but I don't know how to do this. I can create a layer for each considered year but I'm not able to image how to compare with histograms the density among years for each site. Can you please help me? I thank you in advance. Best regards Giovanni


Conversion problems from mappoint to screenpoint using ArcGis Javascript API

I am getting a very strange problem when converting from mappoint to screenpoint. I am using the following code:               var sp = esri.geometry.ScreenPoint(512, 151);         mp = map.toMap(screenPoint);      After converting from screenpoint(512, 151) to mappoint, I get the following map coordinates:               338638.7275175141, 4404774.250335829      Now when I perform the reverse conversion (i.e., from the same point above which is mappoint(338638.7275175141, 4404774.250335829) to screenPoint), I get a very different screenpoint (i.e. 31563417740, -72455974). I think I should get the same screenpoint(512, 151) coordinate. For this reverse conversion I am using the following code:               mp1 = new esri.geometry.Point(338638.7275175141, 4404774.250335829);     sp1 = map.toScreen(mp1);      With it, I get a very different screenpoint as I showed before and no the value (512,151). Any help will be appreciated,


Spreadsheet functions to query route planner for travel time/distance

I would like to achieve something whereby I have a spreadsheet such that the columns are:   * Column A - place name   * Column B - place name   * Column C - distance by road between places in columns A and B   * Column D - travel time by road between places in columns A and B I thought it might be possible using Google Docs' spreadsheet and its 'Google' functions, but I've not found any that might do the trick. In the end I could knock up an app to do it using the Google Maps API but would rather avoid it if I can.


How does Mapnik Map Similar Features to Images

# General Problem I want to be able to select all nodes inside of a bounding box that have similar properties. A good example would be: _Map features that represent vegetaion._ In other words: Basically everything a map renderer usually displays in green (trees/woods/forests/...). # Specific Problem OpenStreetMap data is my primary data source. OpenStreeMap manages a wiki page about vegetation which lists various tages OSM users usually use to tag vegetation. Using the Overpass API I was able to select most of the vegetation that appears green on the OSM web site but unfortunately not everything. The problems are inappropriately tagged map features. Since most users on OSM primarily look whether their added/edited nodes look correctly in OSM, my idea was to extract the specific mapping used for OSM from the Mapnik renderer source repository. There must be some information which maps OSM tags to the Mapnik model which eventually results in green areas. I had a quick look in `plugins/input/osm/` but was unable to find something. # Question Where is the bit of code and the respective data that decides that OSM tags like `wood=mixed` , `landuse=forest` and so forth are eventually displayed as green area (with little dark green trees)? Alternatively: Do you have any other ideas to extract all vegetation that the OSM web site renders greenish? # Disclaimer I'm not familiar with the OSM architecture nor the Mapnik source code, therefore my question might seem silly. So if you only know a little more than me, I'd be very happy if you could point me out the bits of code/data that'd help me :-). Thanks!


How to manually edit point coordinates in QGIS?

Is there a tool for shapefiles or a QGIS plugin available, which will enable me to edit the coords of a point manually ?               ("tool" -f1 "Tsunami Gauge"-move x1.abcdef -to x2.ghijkl)        I do have approximate points of which I the later receive highly preceice data and there I would lie to update my so far data, instead of creating new files.


Qgis is slow with complex 'postgis query' layer. How to speed up?

I have a postgis layer in Qgis that is a result of a complex query. Every screen update, Qgis runs the Query again and that's slow. Any suggestion on how to speed up? EDIT - Here is the query:                SELECT       parcela.id AS id,       st_union(areas.the_geom) AS geom,       (((((('p'::text || round(avg(to_number(dado.valor -> 'pequena'::text, '99'::text)) * 10::numeric, 2)) || ',m'::text) || round(avg(to_number(dado.valor -> 'media'::text, '99'::text)) * 10::numeric, 2)) || ',g'::text) || round(avg(to_number(dado.valor -> 'grande'::text, '99'::text)) * 10::numeric, 2)) || ',mo'::text) || round(avg(to_number(dado.valor -> 'morta'::text, '99'::text)) * 10::numeric, 2) AS "legenda",        round(avg(to_number(dado.valor -> 'pequena'::text, '99'::text)) * 10::numeric, 2) AS pequenas,        round(avg(to_number(dado.valor -> 'media'::text, '99'::text)) * 10::numeric, 2) AS medias,        round(avg(to_number(dado.valor -> 'grande'::text, '99'::text)) * 10::numeric, 2) AS grandes,        round(avg(to_number(dado.valor -> 'morta'::text, '99'::text)) / 10::numeric, 2) AS mortas            FROM petr.dado, petr.parcela, petr.parcela_rel, petr.areas       WHERE dado.parcela_id = parcela.id AND parcela_rel.area_id = areas.gid AND parcela_rel.parcela_id = parcela.id       GROUP BY parcela.id


How do I add a shapefile in ArcGIS via python scripting?

I am trying to automate various tasks in ArcGIS Desktop (using ArcMap generally) with Python, and I keep needing a way to add a shape file to the current map. (And then do stuff to it, but that's another story). The best I can do so far is to add a layer file to the current map, using the following ("addLayer" is a layer file object):               def AddLayerFromLayerFile(addLayer):       import arcpy       mxd = arcpy.mapping.MapDocument("CURRENT")       df = arcpy.mapping.ListDataFrames(mxd, "Layers")[0]       arcpy.mapping.AddLayer(df, addLayer, "AUTO_ARRANGE")       arcpy.RefreshActiveView()       arcpy.RefreshTOC()       del mxd, df, addLayer       However, my raw data is always going be shape files, so I need to be able to open them. (Equivantly: convert a shape file to a layer file wiothout opening it, but I'd prefer not to do that).


Is there an architecture for distributed geoprocessing?

Suppose I have 50 computers on my LAN. Each computer has a geodatabase for all parcel polygons in a particular state in the US. I'd like to write a geoprocessing task that finds all parcels valued over _x_ $/acre that are within _y_ feet of another parcel that is valued at less than _z_ $/acre. I'd like to formulate and run this query without knowing or caring that the data is distributed across 50 computers. Keep in mind the boundary conditions: I also want the query to return cases where expensive parcels in one state are near inexpensive parcels in another. Is there an architecture that supports this sort of distributed geoprocessing? **Update:** The architecture can be described abstractly, or as an implementation specific to Azure or Amazon Web Services. Or, preferably, as a typical office where computers sit idle at night with plentiful ArcGIS desktop licenses.


batch buffer with python in arcgis 10

Hi friends: i have some feature class(7 feature classes) and I want to batch buffer in ArcGIS 10 . I don't know how can I buffer at once. Please help me. sincerely


Free (or cheap) software for generating Geospatial PDF (GeoPDF)

I'm looking for a cheap software to generate Geospatial PDF (GeoPDF, ISO 32000 Adobe extensions or OGC recomendation). Currently, as far as I know, there are very expensive propietary solutions. See http://en.wikipedia.org/wiki/Geospatial_PDF for details. Adobe has retired Geospatial PDF support from Adobe Acrobat, and they propose to use expensive solutions form ESRI or TerraGo. Do you know any cheaper way to generate Geospatial PDF files?


How can I install gdal 1.9.0 on my linux ubuntu server 10.04

`apt-get install python-gdal` installs version 1.6.0 and this is already installed. I would like to install the 1.9.0 version. Do I have to de-install 1.6.0 first and then compile manually or is there an `apt-get` package that I have missed? thanks for any help, Rob


Getting MongoDB working with GeoServer 2

I'm trying to use MongoDB with GeoServer but it is actually not working at all. Someone once posted this: Getting MongoDB Working with GeoServer, but since he could solve his problem without any help there's not much info about how I can do it too. I got the MongoDB GeoTools plugin (what is the correct .jar file and on which folder should I put it?). I'm trying to make the "new data source" option appear like this guy did. If someone could help me in a "tutorial like" way I'ld be very happy :)


Searching an ArcGIS Server 10.1 search service with ArcGIS for Sharepoint

Could someone kindly tell me if it is possible to search an AGS 10.1 Search Service index using ArcGIS for Sharepoint? I am basically looking for a comparable search functionality to the Enterprise Search in ArcMap 10.1. Thanks


How can I merge geotiffs with a vrt file?

I have heard that I can use a vrt file to merge mulitple geotiffs into one? what is the syntax? gdal_tranlsate -of VRT file1.tif file2.tif file2.tif Something like this?


Calculate XYZ velocities between 2 ECEF points

> **Possible Duplicate:**   >  How do I calculate ECEF X,Y,Z velocities for GPS if I have no satellites If I have a starting location (Latitude/Longitude), a bearing (degrees), and a speed (m/s), how can I calculate ECEF X, Y, and Z velocity?


Calculating ECEF X, Y, and Z velocities between 2 points on Earth

> **Possible Duplicate:**   >  How do I calculate ECEF X,Y,Z velocities for GPS if I have no satellites This is my third attempt at explaining what I'm trying to do. If I have two ECEF X, Y, Z points (A and B), and the given time to travel from one point to the other is 1 second, how do I calculate the X, Y, and Z velocity from point A to point B. Taking the difference between point B's X, Y, and Z and A's X, Y, and Z result in a linear velocity, which I'm guessing doesn't happen on earth (Because it's a sphere). So, how would I calculate a non-linear velocity between point A and B's X, Y, and Zs? If my question is correct, I must be looking for the change in rotation between each X, Y, Z points / 1 second. So how would I determine the change in rotation between all three points? I hope this is clear enough. I'm attempting to simulate a GPS signal for which there are no satellites to triangulate from, but what my simulation knows is where to start, what direction to move to get to a end point, and how fast to travel that direction in m/s.


Why doesn't Gdaladdo make the background transparent?

I noticed that the backgrounds of my geotiffs were transparent up until I started added overviews with gdaladdo. Is there a way to combat this?`


what is the reason still I can't be success full, though i follow the dtail procedure?

I am a widowsXP user with QGIS 1.8.0. and I Downloaded the OSGEO installer from http://download.osgeo.org/osgeo4w/osgeo4w-setup.exe and following the detail procedure outlined on the page to install semi-automatic plunge. however, it respond "Unable to set up.ini.from http://download.osgeo.org/osgeo4. on the other hand, when I download semi- automatic plunge to instal manually from http://plugins.qgis.org/plugins/plugins.xml and extract in my directory i.e C:\Program Files\Quantum GIS Lisboa\apps\qgis\python\plugins, and activate it in manage plunge to use, still I can see response that it Couldn't load plugin SemiAutomaticClassificationPlugin due an error when calling its classFactory() method Traceback (most recent call last): File "C:/PROGRA~1/Quantum GIS Lisboa/apps/qgis/./python\qgis\utils.py", line 164, in startPlugin plugins[packageName] = package.classFactory(iface) File "C:/PROGRA~1/Quantum GIS Lisboa/apps/qgis/./python/plugins\SemiAutomaticClassificationPlugin__init__.py", line 65, in classFactory from semiautomaticclassificationplugin import SemiAutomaticClassificationPlugin File "C:/PROGRA~1/Quantum GIS Lisboa/apps/qgis/./python\qgis\utils.py", line 309, in _import mod = _builtin_import(name, globals, locals, fromlist, level) File "C:/PROGRA~1/Quantum GIS Lisboa/apps/qgis/./python/plugins\SemiAutomaticClassificationPlugin\semiautomaticclassificationplugin.py", line 51, in import sextante File "C:/PROGRA~1/Quantum GIS Lisboa/apps/qgis/./python\qgis\utils.py", line 309, in _import mod = _builtin_import(name, globals, locals, fromlist, level) ImportError: No module named sextante 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, f440b60 Python path: ['C:/PROGRA~1/Quantum GIS Lisboa/apps/qgis/./python', 'C:/Documents and Settings/JUCAVM/.qgis//python', 'C:/Documents and Settings/JUCAVM/.qgis//python/plugins', 'C:/PROGRA~1/Quantum GIS Lisboa/apps/qgis/./python/plugins', 'C:\Program Files\ArcGIS\bin', 'C:\PROGRA~1\Quantum GIS Lisboa\bin\python27.zip', 'C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27\DLLs', 'C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27\lib', 'C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27\lib\plat-win', 'C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27\lib\lib-tk', 'C:\PROGRA~1\Quantum GIS Lisboa\apps\qgis\bin', 'C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27', 'C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27\lib\site-packages', 'C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27\lib\site-packages\PIL', 'C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27\lib\site-packages\win32', 'C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27\lib\site-packages\win32\lib', 'C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27\lib\site-packages\Pythonwin', 'C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27\lib\site-packages\wx-2.8-msw- unicode', 'C:\PROGRA~1\Quantum GIS Lisboa\apps\qgis\python\plugins\fTools\tools'] what shall i do? what is my failure?


I have tiff file that has been modified and geotiff file with info, how do i convert my tiff to geotiff?

I have modified tiff file that i would like to turn to geotiff. My orginal geotiff has the following information that i got from gdalinfo: Driver: Gtiff/Geotiff Files: info.tiff Size is 1127, 831 (new file is 4096 3020 in size) Cordinate system is: * * * GEOGCS["WGS 84",                   DATUM[WGS_1984",               SPHEROID[WGS 84", 6378137,298.257223563,              AUTHORITY["EPSG","7030"]],              AUTHORITY["EPGS","6326"]],              PRIMEM["Greenwich",0],              UNIT["degree",0,0174532925199433],              AUTHORITY["EPSG","4326"]]      Orgin = (-17.385864000000002,28.226624000000001) Pixel Size = (0.000277826086957, -0.000277632972323) Metadata: AREA_OR_POINT=Area Image Structure Metadata: INTERLEAVE=BAND Corner Coordinates: Upper Left (-17.3858640, 28.2266240) Lower Left (-17.3858640, 27.9959110) Upper Right (-17.0727540, 28.2266240) Lower Right (-17.0727540, 27.9959110) Center (-17.2293090, 281112675 Band 1 Block=1127x3 type=int16, Colorinterp=gray NoData Value=32768 * * * Now i would like to place this new tiff that will be converted to same place with same information of the orginal data. I think i need to use gdal_translate, i have seen some examples but those are too confusing to understand and i need some wrap thing too i think. Would it be possible for someone to show me the lines i need to use with the provided information, so i could manage to do it by my self next time. I made the new image copy of the orginal, but with more resolution and slight tweaks, the new file is 4096 x 3020. Thank You Very Much!


How to Batch Import Collada (.dae) files to ESRI MultiPatch format?

I was wondering if anyone has come across a tool, or can point to code snippets, for batch importing Collada (.dae) files into a MultiPatch format. I have found a good "How to", but it is all manual. Also with this manual process the Collada file is not imported into is exact geospatial location. Any help would be appreciated. Here is the youtube link to the "How to" http://www.youtube.com/watch?v=udhm3yxVinc Thanks...


How to create a hectare raster from vector input data in QGIS

I have a vector file containing points with an attributed value (no. of employes at that location). I need to create a hectare raster with the sum of employes for each grid cell. I'm using QGIS and have so far played around with gdal_rasterize gdal_grid but with no success. I'm not even sure if I'm looking at the right tool. Can anyone give me a hint?


OpenLayers WMS Request speed against direct request

I'm using OpenLayers and I'm calling an Untiled 3rd party WMS service. I've run some tests to test performance\time differences between the following:   1. OpenLayers WMS call to the service (mostly slowest)   2. WMS URL directly into browser (sometimes slowest)   3. Calling 3rd party WMS service directly (quickest) I'm using the exact same image dimensions and bounding boxes each time for all calls. I'm seeing differences of 1 to 2 seconds between calls 1 and 3. I'm calling the WMS layer correctly using OpenLayers.Layer.WMS... I've done it before but the time differnce is odd. Can anyone enlighten me as to why there would be differences?


How to install Shapely, but got stuck into Python language

I have a MAC OSX 10.7.5 and am trying to learn how to use QGIS. Following a couple of tutorials I came to know I have to install shapely. I followed the instructions by the multi-cited http://tumblr.pauladamsmith.com, but I am afraid I did not understand neither how to start the installation. I downloaded all the required files (i.e. GDAL 1.10 and its python package but did not move them in any particular directory - i.e. the complete GDAL framework seems to be installed into Library and the GDAL python package is still in the Download folder). I did the following, thanks to the instructions of the link above:               MacBook-Pro-di-myname:Downloads myname$ cd /Users/myname/Downloads/GDAL-1.9.1      MacBook-Pro-di-myname:GDAL-1.9.1 myname$ CFLAGS=`/Library/Frameworks/GDAL.framework/Versions/1.10/Programs/gdal-config --cflags`     MacBook-Pro-di-myname:GDAL-1.9.1 myname$ LDFLAGS=`/Library/Frameworks/GDAL.framework/Versions/1.10/Programs/gdal-config --libs` python setup.py install      but then that's what I get:               running install     error: can't create or remove files in install directory          The following error occurred while trying to add or remove files in the     installation directory:          [Errno 13] Permission denied: '/Library/Python/2.7/site-packages/test-easy-install-3861.write-test'          The installation directory you specified (via --install-dir, --prefix, or the distutils default setting) was:          /Library/Python/2.7/site-packages/          Perhaps your account does not have write access to this directory?  If the     installation directory is a system-owned directory, you may need to sign in     as the administrator or "root" account.  If you do not have administrative     access to this machine, you may wish to choose a different installation     directory, preferably one that is listed in your PYTHONPATH environment     variable.          For information on other options, you may wish to consult the     documentation at:          http://peak.telecommunity.com/EasyInstall.html          Please make the appropriate changes for your system and try again.      Where is this error coming from? What I am getting wrong? Since this error impedes me to continue I cannot follow the other instructions at the link above, so I am pretty stuck with my work and cannot continue. I will appreciate any type of help: I am totally ignorant on Python language and programming as a whole, indeed it was an unexpected problem during my QGIS learning. Thanks a lot Fra


How to calculate heuristic shooting star algorithm?

I have a problem with the shooting star algorithm. It's my final project and I need some help. My question is how can I calculate the heuristic shooting star algorithm? I can use the **pgrouting shooting star** , but I don't understand how to calculate heuristic shooting star. I have searched the calculating heuristic shooting star in internet book but I can't find it now. * * * Can you tell me about shooting star algorithm, not just source code , Mr Mario?


Install PostGIS on Ubuntu 12.04

I am trying to install PostGIS on Ubuntu 12.04 and I think everything is installed correctly except for postgresql. I followed this post. I get the following error and nothing I do seems to satisfy the unmet dependency.               dpkg: error processing /var/cache/apt/archives/postgresql-9.1-postgis_2.0.1-2~precise3_amd64.deb (--unpack):      trying to overwrite '/usr/lib/postgresql/9.1/lib/rtpostgis-2.0.so', which is also in package postgresql-9.1-postgis2 2.0.0-1ubuntu1~ppa18~precise1     dpkg-deb: error: subprocess paste was killed by signal (Broken pipe)     Errors were encountered while processing:      /var/cache/apt/archives/postgresql-9.1-postgis_2.0.1-2~precise3_amd64.deb     E: Sub-process /usr/bin/dpkg returned an error code (1)      A little help for this ubuntu/gis novice.


how to get ascii grid from corine land use data?

I'm trying to get ascii grid with land use data (by one parameter - CODE_00). Afterwards I will import it into simulation tool. But I cannot figure out how to convert ascii grid from shapefile. Shapefile comes from CORINE (http://www.eea.europa.eu/themes/landuse/interactive/clc-download) or raster data from CORINE (http://www.eea.europa.eu/data-and-maps/data/corine-land- cover-2000-raster). What could be the best solution for such a task? This is how should look ascii grid(if I get lucky):               ncols         750     nrows         250     xllcorner     20     yllcorner     55     cellsize      500     NODATA_value  -9999     26 26 6 26 26 26......     25 8 11 23 22 .........


Unable to select Vector Layer

I am using the following example, to select a vector layer feature on click. The example demonstrates Highlighting and selecting the layer. I was able to successfully implement the Highlighting of the layer but I am not able to select the layer. This is my code               var highlightCtrl = new OpenLayers.Control.SelectFeature([vector_layer,QuebecSpecies], {                     hover: true,                     highlightOnly: true,                     renderIntent: "temporary",                     eventListeners: {                         beforefeaturehighlighted: report,                         featurehighlighted: report,                         featureunhighlighted: report                     }                 });                      var selectCtrl = new OpenLayers.Control.SelectFeature([vector_layer],                     {clickout: true}                 );          map.addControl(highlightCtrl);     map.addControl(selectCtrl);          highlightCtrl.activate();     selectCtrl.activate();      Where am I falling short?


Finding Satellite Imagery For Use Without Internet Connection

So, I am not a GIS expert. I just write software and I've done some work with ArcMap, ArcGIS Server, and FalconView. What I need is world wide satellite imagery I can use on a local machine. Street data would be nice but not a priority as I really need the terrain more. I am starting from the premise that cost isn't an issue. My mapping software will be installed on a special rugged lightweight laptop that will not have any internet connection. Basically, map data transfer will be by USB flash drive and/or DVD as needed before embarking (and possibly ethernet over USB vi jack or wifi but that can't be counted on). No matter when and where the map application is in use there won't be any internet (not even 3G/4G) at all. Think back woods, high mountains, post disaster etc. I've found the following possible providers that seem to sell image tiles and I'm wondering if there are any I might be missing?   * http://www.intermap.com/   * http://www.pictometry.com/   * http://www.terraserver.com/   * http://earthobservatory.nasa.gov/Features/BlueMarble/ Thanks for any tips or ideas! Matthew


which directory path does SEXTANTE work on?

I find recommendation that the SEXTANTE path shoud follow on http://fromgistors.blogspot.com/2013/06/configuration-of-semi-automatic.html saying "ATTENTION! The version 1.0.9 of SEXTANTE plugin has an issue that precludes some functionalities of the Semi-Automatic Classification Plugin. Please, install the version 1.0.8 of SEXTANTE plugin that is available at http://plugins.qgis.org/plugins/sextante/version/1.0.8/ and replace the SEXTANTE directory in the plugin directory of QGIS." however I get confuse to change (replace) directory. can one help me?


Is there a geo-fence map of Apple's or Google's reverse geocoding?

Is there a map that shows how lat/long maps to locality/county/sublocality/subadmin and possibly covered area?


Place image without using points in QGIS

I am trying to find a GIS software that does this, so I am trying QGIS but other suggestions are welcome. I have a basic shp boundary such as a rectangle. I want to be able to import an unreferenced image (such as a tif or jpg) and place it as a new layer without needing to reference points immediately. I want to be able to use handles or another tool to re-size or stretch it in comparison to my shp. I realize this file is not necessarily referenced in this state even though it is floating over a referenced layer. It does not necessarily have to be saved as a referenced file ever as long as I can see it as an overlay while I work. Basically because the image I am importing is not actual imagery and does not currently match anything on the map, I want to be able to place it in a way similar to placing and stretching an image in Microsoft Word. Here's a sample image I photoshoped of what this feature might look like: ![enter image description here](http://i.stack.imgur.com/Ek7bI.jpg) Is there a way to do this in QGIS, or suggestions of other software?


Generating a Web based map interface from ArcGIS export file

I have an XML export file from an ArcGIS geodatabase with the following fields:   1. OBJECT ID (-)   2. Shape (-)   3. Grid_Area_ha_1 (ha)   4. Shape_Length_1 (m)   5. Shape_Area_1 (m2)   6. Residue_Yield_t_a_Restrict_NEW (t/a)   7. Shape_Length (°)   8. Shape_Area ((°)2) Which is the quickest way to build a Web based map interface using this data file. Thanks


GML Simple Feature compliance test

I'm in the process of defining a GML schema from a database schema, the idea is for it to be GML simple feature compliant. However I can't find a way to confirm this compliance other than reading the spec and manually checking. I'd feel happier if there was an automated tool I could pass the schema to which could check all was well. The existance of this sort of tool is hinted at in the spec (Appendix A) but I can't find one. Does anyone know of such a tool?


Adding countries to database using osm2po one by one

Is it possible to import countries into postgresql using osm2po one by one, mean not whole europe at once but split and import it by countries to lower load, increase speed etc.


Another Student in Need of Career Advice

I would like to apologize in advance for asking yet another "what should I do with my life" type of question. I am currently a senior who is very close to finishing up my B.S. in Information Systems. Over the past 5 months I have become very much interested in GIS, after taking courses in both Geology and Oceanography. I have been spending some time independently studying/learning application systems such as ArcGIS and Google Earth. I have noticed that most GIS related jobs are reserved for Geography, Architecture, and Environmental Sciences grads, which is obviously understandable being that the underlying principles of GIS are rooted in those disciplines. My question is as follows: Do you feel that I would be considered for a position working in this field with an Information Systems degree, or should I move on to other aspirations? I am considering completing a 9 course GIS certificate program offered at one of the local universities upon completion of my degree, though a GIS related internship would be great right now. Thank you for reading, and any bit of advice would be greatly appreciated.


Error: Objects in this class cannot be edited outside an edit session

I'm simply trying to modify the Default Editing Widget ESRI sample to show my extent and service and am noticing the following error in Firebug when trying to add a point:   * {"error":{"code":400,"message":"Unable to complete operation.","details":["Unable to perform applyEdits operation."]}} I also get the following errors in my logs on ArcServer:   * a database error occurred   * objects in this class cannot be edited outside an edit session Any ideas on why my data is locking like this? Or how I can fix it?


3D polyline editing

I have a question similar to Poalo's question.   Polyline ZM: automatic vertex height errors correction.   It could be that I am using the wrong workflow.   I have a 3d polyline feature class in sde.   **EDIT:** I get updates (via a 2d shapefile) to that fc and move the polyline.   I start an edit, edit vertices and snap to the update shapefile polyline.   Now I want to "Update" the z elevations that I moved.   If I am not mistaken the only method to do that with an sde fc. Is to...   1. move features   2. interpolate to a new fc   3. truncate the sde table   4. load new fc with updated z.  Is this the "only" workflow available?   I see add z inormation, and feature to 3d by attribute. But that would assume I have a way to update the attribute with correct elevations. There are too many to do this manually.   Please let us know if there is a way to do this "in" the feature class.   **EDIT2:**   I would like to do this in a gui (not programmatic). I also do not use the vertex only option when interpolating.   The problem domain is that after editing sde I want to update the shapez values with a dem. I don't write code so if it can be done in a toolbox, calculator, or other interface I can accomplish without help with code. I am fairly aquainted with sql and vbscript for field calculation but only to tinker/load.


Create a circle or ellipse from points

Sorry if this is basic. I have an excel sheet with lat/long in it and a field with a radius, in meters, included. I would like to create a circle or ellipse feature class around these points. Example:               ObjID      Lat     Long       Radius(m)     0         26.514  -115.253      6     1         33.516  -111.235      54      I can create a point feature class out of this but i do not know how to create a circle or ellipse bases on the "Radius (m)" field. I am using Arc10.0 Thanks!!


OpenLayers is not adding marker to my map

I have an OpenLayers map and I create it using a JavaScript function. I then have a button which when pressed calls another JavaScript function which has code to add a marker to the map. For some reason the marker is not added. Code for OpenLayers map creation:               function init()         {             var map, baseLayer;             var options = {numZoomLevels: 5, isBaseLayer: true,};                  map = new OpenLayers.Map('map');                         baseLayer = new OpenLayers.Layer.Image("My map","map.png",              new OpenLayers.Bounds(-2000, -2000, 1000, 1000),             new OpenLayers.Size(1500,1500), options);             map.addLayer(baseLayer);             var center = new OpenLayers.LonLat(-100,-100);             map.setCenter(center ,0);               map.addControl(new OpenLayers.Control.LayerSwitcher());             map.addControl(new OpenLayers.Control.MousePosition());             var size = new OpenLayers.Size(21,25);             var offset = new OpenLayers.Pixel(-(size.w/2), -size.h);             var icon = new OpenLayers.Icon('http://www.openlayers.org/dev/img/marker.png',size,offset);             var markers = new OpenLayers.Layer.Markers( "Marker" );             map.addLayer(markers);     }         Code for adding marker:               function addMarker(x, y)         {                 markers.addMarker(new OpenLayers.Marker(new OpenLayers.LonLat(x, y),icon));         }      Any help is much appreciated :)


R: data.frame to raster object error

I have a data.frame consisting of a points with an associated continuos value and their X,Y coordinates. Altough it's not stored anywhere it's projection is a Lambert conformal conic. I think this is why when I try to make it a raster object:               rast <- raster(dataframe)      I get the following error: Error en .local(x, ...) : data are not on a regular grid How can I tackle this issue? Thank you.


How do I get the extent of a rotated data frame in map units from the Layout View?

I need to get the real extent of a rotated data frame in map units from the Layout View using ArcObjects. I've checked previous discussions on the subject: How do I obtain the envelope of a data frame in Layout? How can I transform layout page units to map units? but it seems there was no final answer on how to get corner coordinates for a rotated data frame. Also, the tool provided in the above post does not recognize rotated data frames. It is a bit confusing as even if you rotate the data frame in ArcGIS, the extent, indicated in the Data Frame Properties, remains unchanged. Is there a simple way to get the real extent directly from ArcObject, or do I need to calculate these coordinates using map scale, size of the data frame etc..? Any help would be appreciated.


How can i mark movement route to the map with QGIS?

I am beginner with QGIS and i am having problem with marking movement route to the map. So i would like to ask how can i mark movement route to the map with QGIS. I have a list of cities and i have to show how much people came from one city to another. For example: 5 people came Vantaa to Helsinki, 10 people came from Tampere to Helsinki, 50 people came from Joensuu to Helsinki etc, and i have to put all this kind of data to map. It should be shown with different thickness lines, where 5 is the thinnest and 50 is the fattest. I hope somebody can help me (some tutorial would be great :))


Does ArcGIS 10.1 have a auto-backup feature?

Did ESRI decide to dump the backup file when ArcMap crashes and didn't save the file? Did they change anything to that? Did they change that to MXD Doctor? Whenever I try to exit or save my project in ArcMap, it just crashes and I when I look for the backup file, I see nothing.


Why are ArcGIS Server SOC machines removed from host after Windows Update?

Whenever Windows service pack updated on our servers, SOC machines were automatically removed from HOSTS. Is this a known behaviour? Are there any ESRI tech support articles on this? Machines are Win2008 R2 with ArcGIS 10.


Can I supply my own REST API token to the ArcGIS Javascrpit API?

I am using the ArcGIS Javascript API to connect to the arcgis.com portal. To do so I am calling `portal.signIn()`. The issue is that I don't want the user to be responsible for performing the login, bo I want to generate the REST API token on my server and then use that token from the Javascript API. **Mock code**               portal = new esri.arcgis.Portal('http://www.arcgis.com');     dojo.connect(portal, 'onLoad', function() {         portal.id.token="<server-side writes token here>";         portal.signIn().then(function(loggedInUser){             console.log(loggedInUser);         });     });      EDIT: From what I can tell in this sample http://help.arcgis.com/en/webapi/javascript/arcgis/help/jssamples_start.htm#jssamples/widget_identitymanager_client_side.html I might be able to build my own credentials object to run esri.id.initialize() (which from what I understand is an implementation of esri.IdentityManager). I will try this and report back if it works.


Why won't QGIS let me write to a file?

I am writing a QGIS Python Plugin. I am trying to write data to a file using               with open("filename.txt", "w") as f:          f.write("file text.")     QMessageBox.information(None, "Test", "Written to file.")      It runs the above and displays the message but it does not write to the file. Why is it not writing to the file? No exceptions are raised. Thank-you in advance.


How to develop GUI in Python for geoprocessing?

I want to develop GUI in Python for geoprocessing, but I don't know what is the best GUI tool or library (Tkinter, PyQT, wxPython) for geoprocessing? Can I use PyQT for GUI programming in ArcGIS?


ArcMap UI with Python

I want to add a UI to my arc map using Python. Yes, Arc Map GP is not enough in my case... Thing is, i want to migrate all my VBA to Python. Here are what i already have :   1. Input forms (combo box,text box etc)   2. Data integration with SAP. The existing VBA can take feedback from SAP Web Service and then processed in Arc Map. So i need more than just GP (???). Right.. ? Please kindly give me suggestions.. Thanks


Is it possible to add arcgis-server services or ArcSDE datasets to Tilemill

Our web design department wants to use mapbox/tilemill to create webmaps for our public websites, but all our data is made available through ArcGIS-server or ArcSDE. Is there some way to use these directly in tilemill, maybe through OGR, or serverside script thats converts ArcGIS-services into geojson? I would prefer not set up a parallel storage in shapefile if I don't have to. Extra caveat: our designers are hardcore Mac-users.


OpenLayers IE performance issues

I have an Openlayers map set up, it is all working well except in IE7 (which most of our users have) where performance is slow. I load in several layers of vector features and then switch these layers on/off using a custom 'layer switcher' which is some buttons/icons on the webpage. In Chrome this works well, it's smooth and fast. IE7 is slow and there is usually a big delay when switching on a certain layer with 250 point features. The delay itself is ok, but it's easy for a user to think it hasn't worked and click the button again. I decided to temporarily replace the button image with a loading image, make layer visible, then change button back to its correct image:               function toggle(buttonClicked) {             button = document.getElementById(buttonClicked);             var layer = map.getLayersByName(buttonClicked);             if ( layer[0].visibility === false) {                 toggleButtonImg(buttonClicked, 'loading');                 layer[0].setVisibility(true);                 toggleButtonImg(buttonClicked, 'on');             } else {                 layer[0].setVisibility(false);                 toggleButtonImg(buttonClicked, 'off');             }          }      And the toggleButtonImg function:               function toggleButtonImg(buttonName, status){             button = document.getElementById(buttonName);                    if (status === "loading") {                 button.src = "images/loading.gif";                           }else {                 button.src = "images/" + buttonName + "_" + status + ".png";             }         }      And if it helps, the button code:               <input type="image" onClick="toggle('Postcode')" id="Postcode" src="images/postcode_off.png">      The issue is, when you click the button the image doesn't change, the browser 'locks up' as it is loading the layer, then the button changes to the 'on' version. (off = greyed out, on = solid). I've added some extra code to the toggleButtonImg() in the hope the image would update before loading the layer:               button.src = "images/custSymb/loading.gif";     alert("hey the button should now have changed");      but I click, image stays 'off', alert appears, I ok it, image is still 'off', layer loads and image changes to 'on'. Any ideas on how I can get around this issue and also any other tips to improve IE performance. e.g Will postcode_layer.removeAllFeatures(); be beneficial to perfomance when the postcode layer is not visible? (I realise I will have to reload it each time). Thanks


Differences in costs of open-source hosted GIS software?

I recently saw a listing of the available GIS products at thegismarketplace.com, which led to a couple of questions.   1. What is the difference between a paid version of GeoServer or MapServer and the free ones?   2. I also noticed that the site is charging for QGIS, which I don't understand either since QGIS is free. There is a related discussion here: comparison of open sources gis servers


Cannot reselect polygon after modify

I've got various vector layers on a OL map. On one layer I have a few draw controls. If I draw a polygon I can select it (the whole pop-up thing works). This also works for all the other features. Then I use the modify control to reshape the polygon. When I complete this action (with a clickout) I catch the afterfeaturemodified event and disable the modify control, re enable the selectcontrol for the map. Now the problem occurs. The clickHandler of the selectControl cannot get the modified feature from the event. I used drawlayer.layer.getfeaturefromevent(event). The weird thing is, this does work for all the other features that er still on the map. Doing things like:               drawlayer.layer.refresh();     or     drawlayer.layer.redraw();      doesn't help. I use the infoWindow for a form that actually saves the modified feature, so I need this to work. Also, I can't use the modify control more then once on the same feature. Which is also very annoying from a usuability point of view. If I look at the DOM, the drawLayer still has just one feature. I've got no idea how to solve this. Any help would be nice.


How to determine which WFS server software is being used?

I want to display WFS data from the following OGC compliant service: http://gis.arso.gov.si/wfs?STORE=ARSOGIS I haven't been successful yet, so now I'm wondering which WFS server software is being used by arso.gov.si. How can I find out what's running on their servers?


Simplify Raster Contours ArcGIS 10

Does anybody know a way to simplify a number of contour polylines derived from a raster dataset? I need to generate a map that shows contour polylines at relatively small scale. With the contour interval at a 1 meter interval, contours converge too closely, making the output lines difficult to distinguish at that scale. The spatial resolution of the original raster is 2 meters. I'm looking for something like the 'simplify line' tool. I know this tool exists in the 'Generalization' Toolset, but I haven't got an appropriate extension. I'd be grateful for your suggestions! Thanks in advance.


Create a One-to-Many Spatial View with ArcGIS 10.1

We have two datasets: parcels and owners. This is a many-to-many relationship, where a parcel can have multiple owners and there may be multiple properties on a single parcel (also with multiple owners: condos or strata property) We would like to publish this data to ArcGIS Server, and issue queries on it, so that we can get a single result for a single condo owner. Currently, when mapping the condo owners back to the parcel, we get all other condo owners, which gets confusing if the result was obtained by a text search. In ArcGIS 10 and earlier, it's possible to create a spatial view that generates a unique Object ID field for repeated geometries. This lets us publish a parcels layer with stacked parcels for owners. Is this possible with ArcGIS 10.1? The online documentation implies that ObjectID and Geometry must be obtained from the same feature class in a spatial view. The underlying database is SQL Server. We do not want to convert the data storage type.


Raster and shp overlap not working with user CRS mercator

I created a GeoTiff raster in python, that is projected with the following CRS: *"+proj=merc +lon_0=34.9586601257 +k=1 +x_0=-0 +y_0=1453756.01544 +a=6370997 +b=6370997 +units=m +no_defs"* When I load it in QGIS, I set project CRS from layer and enable on the fly CRS transformation. Then I open a shapefile of the country (Tanzania), but it doesn't overlap, even if I set the same CRS manually. I'm sure that the raster extent is correct. It looks that it is a problem from the shp. Here is a sceenshot: http://imageshack.us/f/193/wrongrastersinqgis.jpg/ It is working if I georeference my raster with decimal degrees WGS84... But not if I use the projection coordinates. That is strange. Screenshot: http://img841.imageshack.us/img841/9500/rastersinqgis.jpg I did the same with another raster I also created in python, that is projected with: *"+proj=lcc +lat_1=51.4259986877 +lat_2=51.4259986877 +lat_0=51.5159492493 +lon_0=2.60400009155 +x_0=163358.880084 +y_0=163358.804499 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs"* And with this raster, I can properly overlap the shapefile of the country! I don't understand why it is working with a CRS in lambert and not with a CRS in mercator. Any idea?


getting "TypeError :bounds is null" in firebug

i've got something error when I run my code in firebug.. here this code :               <script type="text/javascript" src="OpenLayers.js"></script>                  <?php                 mysql_connect("localhost", "root", "") or die("koneksi gagal");                 mysql_select_db("jogja") or die ("database tidak bisa dibuka");                          $query="SELECT * FROM jogja ORDER BY idkab";                     $result= mysql_query($query);                          $value=array();                     $kab=array();                          while($data=mysql_fetch_array($result))                     {                         $value[$data["idkab"]] = $data["value"];                         $kab[$data["idkab"]] = $data["kab"];                     }                          $bataskelas = array(0, 500, 1000, 2000);                     $warnakelas = array("#EEFFFF", "#6DD5FF", "#19A3FF", "#0048FF");         ?>                  <script type="text/javascript">                      var map;                 var jogja;                 var data = <?php echo json_encode($value); ?>;                 var kab = <?php echo json_encode($kab); ?>;                 var batas = <?php echo json_encode($bataskelas); ?>;                 var warna = <?php echo json_encode($warnakelas); ?>;                      function init(){                 map = new OpenLayers.Map('map',{                 controls: [                     new OpenLayers.Control.Navigation(),                     new OpenLayers.Control.PanZoom(),                     new OpenLayers.Control.ScaleLine(),                     new OpenLayers.Control.MousePosition(),                 ]});                      var context = {                     getColor: function(feature){                     var id = feature.attributes.IDKAB;                     var val=data[id];                     var color = "white";                          var i =0;                     while(i < batas.length-1){                         if (val < batas[i+1]){                             color = warna[i];                             break;                         }                         i++;                     }                         if (val !=null && color == "white"){                             color = warna[batas.length-1];                         }                         return color;                     },                          getLabel: function (feature){                         return wilayah[feature.attributes.IDKAB];                     }                 };                      var template = {                     strokeColor: "#808080",                     strokeWidth: 1,                     graphicZIndex: 1,                     fillColor: "${getColor}",                     label: "${getLabel}",                     labelAlign: "cm"                      };                      var style = new OpenLayers.Style(template, {context: context});                 var myStyle = new OpenLayers.StyleMap(style);                      jogja = new OpenLayers.Layer.GML ("GML", "gml/jogja.gml", {                             isBaseLayer: true, styleMap: myStyle});                      jogja.events.on({'loadend': onLoadEnd});                      map.addLayer(jogja);                 map.zoomToMaxExtent();                  }             function onLoadEnd(){                 var ext = jogja.getDataExtent();                 map.zoomToExtent(ext,false);             }                   </script>      and here screenshot error : ![enter image description here](http://i.stack.imgur.com/yPIhA.png) what should I do to avoiding this error?? because of this error the map don't show out.. -thanks for answering-


How can I download a feature class from an ArcGIS Online web map?

I'm new to ArcGIS Online (or ArcGIS.com as they're calling it now), and I'm confused as to how data is organized, and how it is accessible. I created a new map, chose Add -> Create Editable Layer, and set up an editable layer open to the public for people to submit input. Now I want to download this feature class to my computer. The only option I know of is to go to "My Content" then download the entire web map as a map package, but when I open it up, it only has the basemap and not the editable layer. **Is it possible to download this data as an individual feature class (shapefile), or at all?**


I want to georeference jpg map without having any digitizing data?

> **Possible Duplicate:**   >  I want to georeference jpg map without any digitizing data? I have only a jpg map but dont have any digital data for georeferencing that jpg.So is there any way for georeferencing my map (like googlemaps). I know only the lat/long of the map.i'm working with arcgis9.3? could you give me some reference as well.


Either GeoJson to TopoJson conversion problem, or projection-issue with D3.js

My 'handcrafted' cartogram (How to create contiguous area cartograms in QGIS or R?) is ready to embed in code - I want to present it online with D3.js. Therefore I produced a topoJSON file of it - but all I get are weird shapes all around the globe. Here is the (simplified) cartogram in geoJSON _(file:removed)_ - it works perfectly with QGis. Also I checked the geometry in QGis and it doesn't display an error [Vector > Geometry Tools > Check geometry validity]. I tried both: From QGis > geoJSON > topoJSON (with topojson), and QGis > ESRI shapefile > geoJSON (with ogr2ogr) > topoJSON (with topojson). Also I tried reprojection with ogr2ogr:               ogr2ogr "ESRI Shapefile" -t_srs EPSG:4326 outfile.shp infile.shp      And I did the whole thing with UTF-8 flag:               topojson --shapefile-encoding utf8 -o outfile.json infile.shp      Here is the (simplified) cartogram in topoJSON _(file:removed)_ , and it looks ok to me - even the json-validators are happy with it. But my D3.js code just results in weird shapes and lines.               <head>         <meta charset="utf-8">         <script src="http://d3js.org/d3.v3.min.js"></script>         <script src="http://d3js.org/topojson.v0.min.js"></script>         <title></title>         <style>              .unit {             fill: none;             stroke: #111;             stroke-linejoin: round;         }              </style>     </head>     <body>     <script>              var width = 600,             height = 600;              var projection = d3.geo.albers() // also: orthographic & mercator wont work             .rotate([0, 0])             .center([8.43, 46.8]) // Switzerland             .scale(100);              var path = d3.geo.path()             .projection(projection);              var svg = d3.select("body").append("svg")             .attr("width", width)             .attr("height", height);              d3.json("hexmap.json", function(error, hexm) {             var units = topojson.object(hexm, hexm.objects.hexmap).geometries[25];                  svg.append("path")                .datum(units)                .attr("class", "unit")                .attr("d", path);         });          </script>     </body>      I just don't know where I am going wrong... Is it the projection? The conversion? Are the files ok? **EDIT - SOLUTION**   1. use a transferable projection in the first place.   2. the correct ogr2ogr command:               ogr2ogr "GEOJSON" -s_srs infile.prj -t_srs EPSG:4326 outfile.geojson infile.shp      alternative:                    ogr2ogr "GEOJSON" -s_srs EPSG:21781 -t_srs EPSG:4326 outfile.geojson infile.shp


Where does the Style Manager store the files in and Why does the Style Manager does not show some edit color ramps?

![enter image description here](http://i.stack.imgur.com/ufMwq.jpg) Where are Color Ramp files saved? Every time I make edits, it does not show up on the Style Manager. Does it save it to the hard drive or in the Style Manager? I need more information. Thanks


Problems Refreshing View in Editor Extension code

I can't seem to find anyone else with this problem, but I'm hoping maybe I'm missing something obvious! I am writing an Editor Extension add-in using VB.net for ArcMap 10. The feature layer I'm rendering is from an SDE dataset. Once the user starts editing, I use a unique value renderer to display features, omitting some (the client wants features they have classified as 'inactive' to disappear - without being deleted - and my understanding is that a definition query won't work in an edit session, so I have modified the symbology instead). After I set the renderer, I want the display to refresh to show the new symbology. However, ActiveView.Refresh, ActiveView.PartialRefresh, and ScreenDisplay.Invalidate will not work. There is no error, but the screen does not refresh. I have stepped through using the debugger and confirmed that the program execution does reach those lines, but after 'executing' those lines nothing happens. The TOC updates properly, but the view does not. I can manually click the refresh button and it refreshes just fine, but I can't seem to make it work programmatically. I have tried brute forcing it by using ActiveView.Deactivate followed by ActiveView.Activate, and this does cause a screen refresh, but then when I try to draw the point, the mouse pointer and the marker symbol indicating the potential location of the new point are offset by about 4 inches. I'm completely at a loss here. Any ideas? I've included an exerpt from my code below (this is done nearly identically for 4 feature layers, so normally the select case would have 3 more entries, but I simplified it for this post). Thanks!!                   'define the colors to use in the renderers         Dim pOutlineColor As New ESRI.ArcGIS.Display.RgbColor         Dim pStructureColor As New ESRI.ArcGIS.Display.RgbColor              'make the structure outline color black         pOutlineColor.Red = 0         pOutlineColor.Blue = 0         pOutlineColor.Green = 0          'make the structure center color macaw green         pStructureColor.Red = 152         pStructureColor.Blue = 0         pStructureColor.Green = 230              'create the symbols         Dim pStructureSymbol As New ESRI.ArcGIS.Display.SimpleMarkerSymbol              'construct the structures symbol         pStructureSymbol.Style = ESRI.ArcGIS.Display.esriSimpleMarkerStyle.esriSMSCircle         pStructureSymbol.Size = 4         pStructureSymbol.Outline = True         pStructureSymbol.OutlineSize = 1         pStructureSymbol.OutlineColor = pOutlineColor         pStructureSymbol.Color = pStructureColor              'define renderers         Dim pStructureRenderer As New ESRI.ArcGIS.Carto.UniqueValueRenderer              'set the fields used in the renderers - we'll just have one, the 'Active' field         'set usedefaultsymbol to false so that inactive features           'will not draw (i.e., if Active <> 1 or Null)         'structures renderer         With pStructureRenderer             .FieldCount = 1             .Field(0) = "Active"             .UseDefaultSymbol = False             'add "1" to the value list for the renderer - i.e., display "1" values using pStructureSymbol             .AddValue("1", "Active", pStructureSymbol)             .Label("1") = "Active"             'add "<Null>" to the value list for the renderer using the same symbology as "1"             .AddReferenceValue("<Null>", "1")             'tell it we're using our own custom style             .ColorScheme = "Custom"             'tell it the field type for the 'Active' field is not a string (it's an integer)             .FieldType(0) = False         End With                           Dim pLayer As ESRI.ArcGIS.Carto.ILayer         Dim pFeatureLayer As ESRI.ArcGIS.Carto.IFeatureLayer2         Dim pGeoFeatureLayerStructures As ESRI.ArcGIS.Carto.IGeoFeatureLayer              'this makes the layer properties symbology tab show         'the correct interface.         Dim pPropertyPage As ESRI.ArcGIS.CartoUI.IRendererPropertyPage         pPropertyPage = New ESRI.ArcGIS.CartoUI.UniqueValuePropertyPage              Do Until pLayer Is Nothing             If TypeOf (pLayer) Is ESRI.ArcGIS.Carto.IFeatureLayer2 Then                 pFeatureLayer = pLayer                 If pFeatureLayer.FeatureClass.FindField("Active") > 0 Then 'only set the renderer if there is an 'active' field in the layer                     Select Case pFeatureLayer.FeatureClass.AliasName                         Case strStructuresLayerName                             pGeoFeatureLayerStructures = pFeatureLayer                             pGeoFeatureLayerStructures.Renderer = pStructureRenderer                             pGeoFeatureLayerStructures.DisplayField = "Active"                             pGeoFeatureLayerStructures.RendererPropertyPageClassID = pPropertyPage.ClassID                          End Select                 End If                 pFeatureLayer = Nothing             End If             pLayer = pEnumLayers.Next         Loop         My.Document.ActiveView.ContentsChanged()         My.Document.UpdateContents()                     My.Document.ActiveView.Refresh() 'this does nothing              'This also does nothing:         'My.Document.ActiveView.PartialRefresh(ESRI.ArcGIS.Carto.esriViewDrawPhase.esriViewAll, Nothing, pActiveView.Extent)              'This also does nothing:         'My.ThisApplication.Display.Invalidate(My.Document.ActiveView.Extent, False, ESRI.ArcGIS.Display.esriScreenCache.esriNoScreenCache)              'this makes the screen go blank:         'Dim pMaps As ESRI.ArcGIS.Carto.IMaps         'pMaps = My.Document.Maps         'My.Document.ActiveView = My.Document.PageLayout         'My.Document.ActiveView = pMaps.Item(0)         'My.Document.ActiveView.Refresh()              'this makes the cursor and the point to be drawn offset by about 4 inches:         'Dim pActiveView As ESRI.ArcGIS.Carto.IActiveView         'pActiveView = My.Document.FocusMap         'pActiveView.Deactivate()         'pActiveView.Activate(My.ThisApplication.Display.hWnd)


Building a disconnected web based mapping application

I have been tasked with building a disconnected, web-based mapping application. It must run in a browser, but it will run on a disconnected Windows based laptop. It needs to be able to do simple map displays as well as display of objects on map and historical displays of the route of an object. Does anyone have any suggestions about a particular software solution? I am guessing that I will have to install some sort of mapping server software that has all of the mapping information and then also build a web-based client that accesses the local mapping server software. Does that make sense? UPDATE: Here are more detail about the rough system requirements.   * The finished application will be web-based (runs in a browser)   * The entire system will be disconnected from the internet. So it can't rely on any internet-based mapping system like Google Maps.   * The main screen will display a map of the earth. It needs to show basic country boundaries but there is no requirement for any other information.   * The user should be able to zoom in and out. There will be some zoom limits although I don't know exactly what they are. The application is displaying flight paths for airplanes so hopefully that gives you some idea of the zoom needs.   * The default view will keep the icon for the airplane centered and the map will move beneath.   * The user should also be able to pan the map in any direction and then return to having the airplane icon centered.   * The system will need to be able to dynamically load some overlay maps that display some polygon type shapes.   * The main screen will also have text areas that display other types of information such as current airplane location, etc.   * Ideally the mapping software would be easily integrated using C# and ASP.NET These are the basic rough requirements. I have done a LOT of searching on the web for a solution. Most of the mapping software applications that I have found seem to be overkill. They are complex GIS systems that cost lots of money. I just need a simple mapping application like Google Maps that works offline and allows for the integration functions that I have outlined above. Does anyone have any recommendations? (I mistakenly started a new question that is basically the same thing as this. Here is a link to the second question: Looking for web-based mapping software recommendation)


arcpy spatial join input parameter problems

Using ArcGIS 10 SP4 with all extensions... I have a parent directory (J:\...\spatialjoin) containing two sub-directories of feature classes I want to spatially join. Directory 1 (target features) contains 14 point fc's of mosquito trap locations dispersed around a particular county. Each fc represents one week of trap samples. Fields included: trapname, lat, lon, collection date, # mosquitoes, species. Directory 2 (join features) contains 14 point fc's of rainfall amounts at 219 even-spaced points for a particular county. Each fc represents one week of rainfall. Fields included: point/stationID, lat, lon, weekly sum of rainfall at that point I want to spatially join one week of rainfall to one week of mosquito locations and repeat until all 14 pairs of files have been joined. Filename Example: target directory: traps_2008_week01_7292008.shp, traps_2008_week02_852008.shp, ..., traps_2008_week14_10282008.shp join directory: week_01.shp, week_02.shp, ..., week_14.shp The code below lists multiple workspaces and prints names features in each sub-directory. How can I modify it to input target and join features into the tool? Thanks               import arcpy, os     from arcpy import env          # Allow overwrites     arcpy.env.overwriteOutput = True          dir = r'J:\Spring2012\gis_project\spatialjoin'     arcpy.env.workspace = dir          outWorkspace = r"J:\Spring2012\gis_project\spatialjoinoutput"               try:              targetFeatures = r'J:\Spring2012\gis_project\spatialjoin\collectionweek2008'         joinFeatures = r'J:\Spring2012\gis_project\spatialjoin\weeklysum2008'         outfc = os.path.join(outWorkspace)                   def fcs_in_workspace(workspace):           arcpy.env.workspace = workspace           for fc in arcpy.ListFeatureClasses():               yield os.path.join(workspace, fc)           for ws in arcpy.ListWorkspaces():               for fc in fcs_in_workspace(os.path.join(workspace, ws)):                   yield fc              for fc in fcs_in_workspace(dir):             arcpy.SpatialJoin_analysis(targetFeatures+".shp", joinFeatures+".shp", outfc+".shp", "#", "#", "#", "CLOSEST", "2.392 Kilometers","#")          except:         # If an error occurred while running a tool, then print the messages.         print arcpy.GetMessages()


Topology problems with QGIS/GRASS

*edited file links below I'm working on several NGO projects at the moment, all of which are reliant on being able to use a district level map with some 490+ features, and I am running into some trouble with my geometries. I spent most of the weekend trying to figure out how to snap the features together in this file (https://dl.dropbox.com/u/49434182/492%20Kabupaten%20Google%20Mercator.zip), to get rid of overlaps and slivers. Running the fTools validation tool, it tells me that I have more than 200 problems with this file. So, being a sensible person, I wanted GRASS to do the work for me. However, no matter how I tweaked the parameters of v.in.ogr, I couldn't eliminate the gaps or slivers when I went to export it again. GRASS would generate one layer of (basically) the geometries I wanted, and two layers of slivers. I was able to get my hands on an official shapefile from the national statistics agency (found here: https://dl.dropbox.com/u/49434182/indokabkot_497_2010.zip), which although way more detailed, also has invalid/overlapping geometry according to fTools. I'm unsure of the right projection with this file, so that's also a bit of an issue. I'm using QGIS 1.8, and am at a bit of a loss as to where I've gone wrong. Can anyone provide pointers before my data is due for a presentation next week? Thanks!


Manually edit cell value in DEM in Arcmap10?

Is there any chance that I could manually pick a group of cell in DEM, then edit the value of them in Arcmap 10? I know about the Aris Grid Editor for Arcmap but it requires licence and it's too expensive..:( Thanks in advance for your answer!


WMS Requests - does 1 request = 1 response tile

I'm using Openlayer and GeoServer when I pan or zoom the map - OpenLayers kicks of a series of WMS requests. Does 1 WMS request = 1 tile in response. I'd have thought if you sent a WMS request with the map extents then we'd get multiple responses back to reduce the number of requests - is this possible?


Zoomin to Polygon

1) I have a list of records(states) in a table and an Hyperlink, when user clicks on Hyperlink(state hyperlink) I need to zoom in to repective polygon(state) on layer. How can I acheive this? 2)My layer inside the div is not central aligned. How can I make it central align? Srikanth


what dictates the drawing order of overlapping features with in the same layer/feature class in ArcGIS?

I know that in ArcGIS drawing order is controlled by layers, and if I want to make some features be drawn on top of other features I can use definition queries to separate them into different layers and position one on top of the other. But i was wondering (mostly out of curiosity) what dictates the drawing order of the 1000 overlapping polygons in my feature class, and if there's a way that I can control what features are displayed on top within one layer? Thanks


Isolate the largest patch in a raster in GRASS?

I have a raster with a series of patches (contiguous areas with same value). How do I isolate the largest patch? I know I could run r.stats, look through the long list of calculated areas, figure out which one is the largest and then use r.reclass to isolate it, but I'm sure there has to be an easier or more efficient way to do this.


Line defined by two geographical points crossing circle area with defined radius

I'm trying to calculate the amount of lines defined by two geographical points crossing a circle area defined by a center and radius. ![enter image description here](http://i.stack.imgur.com/Wb6es.jpg) I've tried to use algorithm discussed here, but in my case it just counts too many lines. Probably it's because of the fact that I'm trying to calculate it in Cartesian coordinates and points are defined in WGS 1984. Is there any other approach which could help my to calculate the lines crossing the circle?


matching up fields in different tables

I am going to attempt to explain my problem, and hopefully someone out there will have some insight/solutions/work arounds. I work for a University, and we have some utility/site people go out with their GPS and shoot points pretty often when there are new projects, construction, etc. The data dictionary that they have created has some field names that are longer than 10 characters, so when they export to shapefiles, the fields are being cut off to only contain the allowed limit. Once they bring in the shapefile into our existing base map, they copy and then paste the new items into the existing sde feature class, the data that is contained in the truncated fields won't copy. So, then they have to go in and redo some of the fields that didn't copy over. I know that fields have to match between tables in order to copy/paste data, but is there some kind of script, code, or work around that can match up these fields? I was hoping that there is some kind of field name match. for example, something that i could use to say that create_dat from the gps export file would equal create_date in the existing feature class. I hope that I made some sense in trying to explain this. Any advice or help would be greatly appreciated. I have attached an image of both tables and what data will and won't copy over. ![enter image description here](http://i.stack.imgur.com/5GOWV.jpg)


How to specify the pixel depth in model builder?

I have a simple model that calculates the NDVI from bands 3 and 4 of a raster dataset. The input bands are 8-bit unsigned, whereas the default output NDVI is signed 32 bit. How can I specify the pixel depth for the output NDVI as 8 bit unsigned? ![enter image description here](http://i.stack.imgur.com/7GDp3.jpg)


How to embed a PDF into an Excel spreadsheet?

This is what I want to do:   1. Create report using a polygon feature class to output areas of polygons using some logic and the ExportReport() function (This I have accomplished)   2. Create a PDF of the Polygons. Nothing fancy, I just use python to create a map, add a polygon layer to it and label using the OBJECTID and then use the arcpy.mapping.ExportToPDF() function to create a pdf (This I have accomplished)   3. Finally, I would like to embed the PDF created in Step 2 into into the spreadsheet created from the ExportReport() function in Step 1. I basically want to have this labeled PDF'd map embedded into the spreadsheet so that my end user can see which areas in the table are associated to which polygons in the PDF. I'm using ArcGIS 10.1 with Python.


How to add links to PDFs in a Shapefile?

I have a problem. I want an attribute in a shape file named "Link" to add. The works are all fine. But then I would like to add a link in my directory. This link accesses a PDF. Now, this PDF is not shown. Is there any other solution or there is nothing to be done? Thank you very much


Create a hyperlink to a pdf file in the attribute table

I want to add a hyperlink to a pdf file into the attribute table of a vector layer. I already read the documentation about evis and the answered questions about this here. It didn't seem difficult, but I still can't create a hyperlink...can someone give me a short step by step information about this? Unfortunately I have no access to youtube... I'm using Qgis (2.0.1) and evis is a plugin that should allow to create such hyperlinks via a event browser. I first added a column "Link" to the attribute table and put the link to the document in this column. this is the link to the documentation for evis: http://www.qgis.org/en/docs/user_manual/plugins/plugins_evis.html#figure- evis-2 Thanks a lot! G.H.


How to create stacked bar/column charts in QGIS?

I know that conventional bar charts (with the bars adjacent to each other) based on feature attributes can be plotted in QGIS. But can you please tell me if there's any way to plot the bars in the form of stacks, as in stacked bar/column charts, in QGIS? I tried to find an example in the QGIS manual and on the internet but wasn't able to get any information about it. Thank you.


Enter 4 Coordinates and get a Polygon QGIS

Is it possible to enter 4 Coordinates for getting a Polygon. I have to Clip lots of exact rectangles and i was just wondering wether you know a plugin. Thanks a lot danke! Alex


How to reproduce Stamen's Toner PostGIS setup?

I'm trying to render Stamen's toner map with Landez which uses Mapnik. The style I'm rendering is from Stamen's toner project. I'm using the one which only renders labels (no lines). I imported the portion of osm (not planet) using osm2pgsql and I've also ran the motorway script. I have not imported the natural earth, because I cannot download the data (https://github.com/Citytracking/toner/issues/53), but I doubt I need it for text labeling only. In the style its asking for map source from table (or perhaps view) such as "admin1_labels_50m_z5", but I have no idea the sources of these tables or views? How to create them?


QGIS adding points using decimal degrees

As duck blind permits are sold, we need to add the location using lat/long. we are working in ArcView now but would like to move to QGIS- Everything I was able to find seemed to deal with importing delimited text files. I was wondering if there was a tool like the "add a point using x,y" in QGIS like the one in ArcGIS. The user will need to add single point to the layer as the permit is issued, so I guess I'm looking for a way to edit a point file one record at a time using the coordinates.


Creating a point layer from lat/long coordinates

I have a spreadsheet that lists the lat/long coordinates of several locations I wish to map. I would like to produce a geocoded point layer in QGIS. How do I do this? Thanks!


How to put marker on vector layer with lat/lon on QGIS

I have a vector layer, and bunch of lat/lon values. I want to put markers on the vector layer for each lat/lon pair. I have never used QGis, would you help with it. thanks.


Issue openlayers and tilestream

Basically I created a MBtiles using Tilemill and it generates a map with six zooms. I implement it into tilestream. I added a layer on the openlayers module(Drupal) in XYZ:               site:8888/v2/geography-class/${z}/${x}/${y}.png       The problem is this URL sent by openlayers: `site:8888/v2/geography- class/-1/x/x.png` Each time, whatever the zoom I am on the map, the zoom value always stays at -1, the others values are okay. If I try to reach `site:8888/v2/geography-class/1/1/1.png` it works. What did I miss?


Change coordinates of a polygon

I wanted to know whether it is possible to change the coordinates of a Polygon by typing them in. I don't want to have to manually drag them into the right position.


What is the correct spelling: CSW or CS-W?

In different publications I read sometimes CSW or sometimes CS-W for the standard of transferring metadata-documents. What is the correct 'official' spelling?


How can I conduct Geographically Weighted Principal Component Analysis using ArcGIS, Python and SPSS/R?

I am after a description/methodology for conducting a Geographically Weighted Principal Components Analysis (GWPCA). I am happy using Python for any portion of this and I imagine SPSS or R being used to run the PCA on the geographically weighted variables. My dataset is composed of roughly 30 independent variables that are measured throughout ~550 census tracts (vector geometry). I know this is a loaded question. But, as I search and search, there does not seem to be any solutions out there. What I have come across are mathematical equations that explain the fundamental composition of GWPCA (and GWR). What I am after is more applied in a sense, that I am looking for what major steps I need to accomplish to get from raw data to the GWPCA results. * * * **I would like to expand on the first part with this edit due to the comments received below.** To address Paul... I am basing my interest in GWPCA off of the following paper: _Lloyd, C. D., (2010). Analysing population characteristics using geographically weighted principal components analysis: A case study of Northern Ireland in 2001. Computers, Environment and Urban Systems, 34(5), p.389-399._ For those who do not have access to the literature, I have attached screenshots of the particular sections which explain the mathematics below: ![Article](http://i.stack.imgur.com/WTvhw.png) And to address whuber... Without going into detail (confidentiality), we are attempting to reduce the 30 variables, which we believe are all very good indicators (albeit globally), to the set of components with eigenvalues greater than 1. By computing the geographically weighted components, we attempt to understand the local variances explained by these components. I think our primary goal will be to proof the concept of GWPCA, that is, show the spatially explicit nature of our data and that we cannot consider all independent variables to be explanatory on a global scale. Rather, the local scale (neighbourhoods) that each component will identify will aid us in understanding the multi-dimensional nature of our data (how variables can be combined with one another to explain a certain neighbourhoods in our study area). We hope to map the percentage of variance accounted for by each component (separately), to understand the extent of the neighbourhood explained by the component in question (aid us in understanding local spatiality of our components). Perhaps some other mapping examples but none come to mind at the moment. Additionally: The mathematics behind the GWPCA is beyond what I understand given my background in geographic analysis and social statistics. The application of the maths is most important, that is, what do I plug in to these variables/formulas.


ArcGIS API for JavaScript - Adding Point(s)

I need to add a single point, I am new to JS and ArcGIS API. I was looking for some code examples on how to create a point and/or a few other points (as needed) I can't use any additional pages/external sources, as in the code needs to be in my single page. I found examples like the button add a point, but I don't need a button I just need to place a point or points with their respective coordinates. Then style as necc. Can someone provide either code or a link so I can do this simple task? Thanks


Spatial regression

I have a question regarding performing spatial regression on data which residuals are spatially correlated.   1. If the residuals are spatially autocorrelated due to the presence of trend, how should I account for the existence of spatial autocorrelation?    2. And what about if there is no trend, but rather spatial autocorrelation in random variation (in residuals left after detrending)?    3. And what about if there are both trend and spatially autocorrelated random variation in my data?    4. Can I in both cases apply this formula to account for spatial autocorrelation no matter if it is the result of trend or random variation correlation: > >     y^*=y-ρ∑(w_ij*y_j) >     x^*=x-ρ∑_(w_ij*x_j) >   Thanks.


Leaflet wms offset on zoom out

I have Leaflet map with base layer (epsg4326) and wms layer of gps point (also epsg4326). I've noticed when I zoom out , the gps point offset to the left. Any idea why it happens ? Here you can see that when I zoom out the symbol moves left each zoom. ![enter image description here](http://i.stack.imgur.com/AikRV.png) ![enter image description here](http://i.stack.imgur.com/YkOXF.png) ![enter image description here](http://i.stack.imgur.com/AL1w3.png) This is the code:               function initLeaflet() {                 osm = L.tileLayer.wms(geoserver + "/OSM/wms", {             layers: 'OSM',             format: 'image/png',             transparent: true         });              coords = L.tileLayer.wms(geoserver + "/coords/wms", {             layers: 'Roham:DeviceInfo',             format: 'image/png',             transparent: true                 });              var baseLayers = {             "osm": osm         };              var overlays = {             "coords": coords         };              map = new L.Map('map2d', {             center: new L.LatLng(32.05522664582287, 35.05812072269409),             zoom: 9,             layers: [osm, coords],             zoomControl: true         });              L.control.layers(baseLayers, overlays).addTo(map);     }


postgis under win 7 "on demand" not as a service

how do I install postgis2.0 under win 7 to run just on demand. I dont want to run it permanently as a service. Similary as the (nonspatial) mysql can be run on demand via the xampp-tools (http://www.apachefriends.org/de/xampp.html)


Updated instructions for setting up TileCache on IIS 7.5?

This article seems great...but it is hard [for me] to follow when setting up for IIS 7.5: http://viswaug.wordpress.com/2008/02/03/setting-up-tilecache-on-iis/ Any updated instructions available with latest versions of TileCache and IIS?


Geowebcache REST 404's

I'm running Geoserver 2.1.3 (SVN 16668) and am trying to get REST working with GeoWebCache. I've got :8080/geoserver/rest working correctly, and am able to interact with that side of things. However, :8080/geowebcache doesn't exist. I'm told that the only two contexts that exist are /javadocs and /geoserver. I understand that there is such a thing as a REST api for geowebcache. How do I get to it? Is there some additional plugin that I need to download for GeoServer before it becomes accessible, or something of that sort? I haven't been able to find anything on the net about someone just missing the REST section of GeoWebCache. The rest of it works fine. Any ideas? Thanks!


Arcpy - How do I specify the precision and scale of a new field?

I'm wanting to create a field with a maximum of 1 decimal place and then calculate a value into it: (I know that all my values are < 10)               arcpy.AddField_management('C:/shapefile.shp', 'myField', "FLOAT", 2, 1, "", "", "", "NON_REQUIRED", "")     arcpy.CalculateField_management('C:/shapefile.shp', 'myField', "[MEAN]", "VB", "")      It seems that no matter what I try I can't get the number of decimal places I want - I get the full value from the field I'm calculating from : [0.888889] _and_ when I view the field properties I just see 0 for precision and 0 for scale. The data context is ESRI shapefile. And altering the numeric settings in the field properties is not an option. ![enter image description here](http://i.stack.imgur.com/YXaaq.png) Anyone have any ideas?


Add layer from a folder to Map using java script api

I know you can add layers to map through map service using javascript api. Can we directly add layers (placed in a folder on ftp) to the map using javascrip api. If so, can please explain it and provide me a sample. Thanks.


Are there any simple map data sources like world.dat

I recently came across this blog post www.gnuplotting.org/plotting-the-world/ that used a file that holds a basic map of the world. Is there a name for this kind of data. Also is there a repository of similar but higher resolution data. Direct link to data: www.gnuplotting.org/data/world.dat