title: Styling a raster on relief shading analysis

text: > **Possible Duplicate:**   >  How to create a style file for QGIS from a plain text file? I have a text file which contains the following data               eleva     R      G       B     0   153 255 153     100 179 255 179     200 237 250 230     300 255 242 242     400 255 242 237      and I need to convert the text file into a saved color style.How best can I achieve this in QGIS


title: What is the best way to automate the importing of symbology into feature classes in ArcMap 10?

text: I'm importing several shp files with the same schema using the FeatureClassToGeodatabase tool in Arcpy and during the process I would also like to apply some standard symbology to the feature class. I was thinking I could use the ApplySymbologyFromLayer but it always errors out when I try to do it with a feature class in a GDB. It works fine with shp files when using the ArcMap built-in tool. So something is different when using a feature class within a GDB for the ApplySymbologyFromLayer. So can I not keep a symbology stored with feature class when it's sitting in the Geodatabase? Meaning the symbology has to applied when I add the feature class to the data frame? Or can I use a Representation by using the Add Representation tool? When I use this it seems to always import the same symbology classification values from the original feature class the representation was created with. So I'm really at a lost here. Any direction would be appreciated.


title: How to create arcs between two points in Python?

text: I've been working on an python script automating the visualization of live point data. I have a table with two coordinate pairs. I'd like to connect these two point with an arc of a circle or similar (parabola??). I've been able to connect the two points with a straight line, but it gives me a boring visual. One of the major hurdles is my license level: ArcView. Anyone have any idea how I can generate a set of points that represent the path of a line between these two points? I can then use the points to polyline command. I know this is a tough one, but I know there are some amazing math whizzes in this community!


title: How to recategorize features in GRASS by geometry, like the ArcGIS Sort tool?

text: The ArcGIS Sort tool creates a new feature class with features sorted by geometry. For example, it can number features IDs from north to south, etc. Is there a way to renumber the `cat` field in a GRASS vector map based on geometry in a similar way?


title: Label Features with Arcpy

text: I had run the contour tool in my script to generate polylines.               Contour("rectExtract", "C:/fakepath/Class1.shp", 50,0)      Now I want to label the polylines with the contour-value. If I execute this script:               layer = arcpy.mapping.ListLayers(IMXD, "")[1]     layer.showLabels=True     arcpy.RefreshActiveView()      it will label the polylines with the ID-Value. How can I change this? Thanks for your help.


title: Does ArcGIS for Server support others using own credentials (Single Sign On) for Federated Services?

text: I have a number of secured ArcGIS 10.2 for Server (AGS) services which a number of other companies use as well. At the moment we have created domain accounts on our domain which they need to access the services. However, does ArcGIS 10.2 for Server allow me to setup my site as a federated resource partner so that the other companies can use their own credentials instead (Single Sign On). I know that arcgis.com and the Portal will allow you to do it but haven't found anywhere about doing it on your own AGS installation?


title: QGIS not interpreting Lambert Conformal Conic projection correctly

text: I'm a fairly new user of QGIS and I've run into an issue that I can't seem to solve on my own. I'm trying to work with the NRCan Atlas of Canada Waterbodies SHP file but I'm having an issues re-projecting other datasets to line up with this one. When I reviewed the .PRJ file packaged with the NRCan data, it indicates it is Lambert Conformal Conic projection with two standard parallels at 49°N and 77°N (see .PRJ text below); however, when I add the file to QGIS (v1.8.0-Lisboa) and view the layer properties, the CRS is displayed as follows: > USER:100001 - * Generated CRS (+proj=lcc +lat_1=49 +lat_0=49 +lon_0=-95 > +k_0=1 +x_0=0 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs) Note that the second standard parallel (which should be in there as '+lat_2=77') is missing... I'm hoping this explains why my other datasets aren't lining up when I project them, but I can't seem to get QGIS to recognize the second parallel. Suggestions? Or is this second parallel not an issue for re-projecting data? * * * ## Unmodified NRCan .PRJ file contents > > PROJCS["NAD_1983_Canada_Atlas_Lambert",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Lambert_Conformal_Conic"],PARAMETER["False_Easting",0.0],PARAMETER["False_Northing",0.0],PARAMETER["Central_Meridian",-95.0],PARAMETER["Standard_Parallel_1",49.0],PARAMETER["Standard_Parallel_2",77.0],PARAMETER["Latitude_Of_Origin",49.0],UNIT["Meter",1.0]]


title: Distance calculation PostGIS

text: I've got the following task to solve: I've got two tables; one contains buildings (houses) as points and the other contains train-stations (at_bhf_classified) as points. I need to find out the distance from each house to the next trainstation. Well, this is more or less easy. But now the issue: The trainstation are categorized. That means: Give the distance to the next station of category 1 if within a 25km radius. If no category 1 station is found within this distance, give me the distance to cat 2 station within the 25 km radius. So far I use the following statement to calculate the nearest category 1 station:               WITH dist_min       AS (SELECT h.h_key,                  Min(St_distance_sphere(h.geom, b.geom)) dist_min           FROM   house h,                  at_bhf_classified b           WHERE  St_distance_sphere(h.geom, b.geom) < 25000 and b.kat=1          GROUP  BY h.h_key)      SELECT h.h_key,         b.bhfname,        b.kat,        b.bhfnummer,        round(cast(d.dist_min as numeric),2)     FROM   dist_min d,         house h,         at_bhf_classified b      WHERE  St_distance_sphere(h.geom, b.geom) = dist_min         AND St_distance_sphere(h.geom, b.geom) < 25000     GROUP BY h.h_key, b.bhfname, b.kat, b.bhfnummer, d.dist_min     ORDER  BY h.h_key       Do you know how to include the second condition explained further up? Thanks for you help! Carsten


title: Identify Route Locations not working in ArcMap

text: When I click Identify Route Locations on a polyline m in ArcMap (10.2.1) it's not brining up the box with the value/measurement. This was working fine about 5 minutes ago, I've tried restarting ArcMap and re-opening the mxd but to no avail! Has anybody encountered this before or do they know of a fix? It flashes like it is working so I think it's something to do with displaying the result as opposed to finding the location.


title: How to define layer order in Openlayers?

text: I have a marker layer, two vector layers and two raster layers. The marker layer gets covered by the two raster layers. Is there a way to move the marker layer always on top? **Edited:** The order of Z-indexes works well in the vector layer, but when I add another WMS raster layer on the map with a lower Z-index, the raster layer still covers the vector layer. **Edited:** Could anyone point out the difference of the rendering order in terms of "LayerIndex" and "Z-Index" of a Layer? I think the layer index deal with the overlay order and the Zindex deals with the feature order in a vector layer, am I right?


title: ArcGIS for Android API Example/Tutorials

text: I know that the ArcGIS API for Android is still in beta (I've joined the beta program), but I was curious if anyone had produced or discovered any examples/tutorials/cookbooks for the API. The Help Resource is still severely limited. I have watched the ArcGIS API for Android Videos, but they do not provide any of the source code (contrary to what is stated in the videos themselves). I have also browsed through the API Reference but I've found the documentation to be lacking.


title: Why is my DEM not caching correctly?

text: Relating to a previous post, found an excellent source for data called Koordinates. I downloaded this 5m DEM of New Zealand, Wellington.   This is in a TIFF format and I have added to ArcMap 10, applied a colour ramp, some hillshading, and applied a standard deviation on the histogram: ![enter image description here](http://i.stack.imgur.com/D0JCu.png) I know that the MSD format does not support hillshading, so I published direct from the MXD. As a dynamic service, it draws fine (but slow as you would expect). So I cached the service using Mixed mode, and the compact storage type. Viewing the cache in both ArcGIS Desktop and REST, the resulting cache is a solid red: ![enter image description here](http://i.stack.imgur.com/KzRQn.png) At some of the closer levels, the expected symbolgy comes through in some tiles: ![enter image description here](http://i.stack.imgur.com/ErWLj.png) Have tried re-creating several times and still getting the same problem. **Would like some advice on troubleshooting this.**   Have not seen anything related in the Esri forums.


title: Is there any way to match images against google maps satellite view?

text: I've got an image that was snapped from google maps satellite view. I know the general location (within a county), but need to find the specific location of the top-down photo. Are there any tools that can search the google maps satellite images in this given area for a match on the image I have?


title: How to get WMS layer title instead of layer name from Geoserver?

text: I specified the raster file with "name" and "title" in Geoserver. For example, I have these in Geoserver:               name,  Title       AT:usa, "This is a map of USA"        Can I return the title ("This is a map of USA" ) as the layer name in Openlayers layer switcher when I call the WMS service for this layer by name? By default, it shows "usa" in the layer switcher.


title: How to set a custom center and zoom level for Zoom World icon in OpenLayers?

text: OpenLayers Zoom Pan Control has a Zoom World Icon. When it is clicked the map zooms all the way out. Is it possible to override this behavior and define a center and zoom level to which the map should be reset? If yes how? PS: My Base Map is Bing Map Thanks Shaunak


title: Set a theoretical (not empirical) variogram in ArcGIS Geostatistical Analyst

text: I'm trying to do some Unconditional Gaussian Simulations using the Gaussian Geostatistical Simulation tool (see here). However, I don't have any data that I am doing this from - just a theoretical parameterisation of a variogram model (lag, sill, nugget etc). The Gaussian Geostatistical Simulation tool requires "an input geostatistical layer resulting from a Simple Kriging model". Is there a way to generate one of these without any data underlying it, and therefore purely from a theoretical variogram model? It seems that the Gaussian Geostatistical Simulation tool will only be using the variogram parameterisation, not the data, so it seems like it should be possible - I just can't see how!


title: Setup a map server that runs from USB to whatever computer it is connected?

text: I would like to have a USB that I can bring with me to any computer. That USB should hold a map server, so that the client computer wont need internet and avoid using OSM or Google maps or any other online map service. I have another portable application that requests a WMS service to load maps, but i would like to direct that import to my local USB map/WMS server. Is PortableGIS the solution? I am using PortableGIS for QGIS and that works fine. However I cant find any documentation...


title: Load third-party APIs only when the corresponding map layer is requested?

text: OpenLayers allows multiple map layers to be hosted on a common platform. If the map application supports a wide variety of maps, say from Google, Yahoo and Microsoft, I find that upon initialization, the APIs from these map providers are loaded as well, regardless of whether I use them or not. I am not sure if the calls to these third party javascripts count towards the usage quota of these map service providers (seems to be). Besides, there is also an additional delay in loading of the web page. Having said that, I am looking for a way to put a placeholder in the layer switcher such that I can avoid declaring the third party map layer class and hence avoiding a call to the third party API script.


title: How to generate Buildings ID numbers inside parcels

text: In order to decrease manual data entry, I would like to automate a task to generate numbers for all buildings inside a parcel starting with 1. For example, I have 4000 parcels and each parcel has many buildings. An example of what I would like to achieve is shown in this diagram: ![enter image description here](http://i.stack.imgur.com/0z1LE.png)


title: Why can't I add GlobalIDs to the table and feature classes I created?

text: I'm creating a script to which make a blank geodatabase with a required dataset, feature classes, tables, subtypes, relationships. I've got everything to work expect for the damn GlobalIDs! I have no idea why it's not working.               #Creates GlobalID's for tables     arcpy.AddGlobalIDs_managment("Equipment", "Copper", "WorkOrder")      These are tables that sit directly inside my geodatabase. My workspace is set directly to the geodatabase by the user:               GDBpath = arcpy.GetParameterAsText(0)     arcpy.env.workspace = GDBpath      I keep getting this error:               <type 'exceptions.AttributeError'>: 'module' object has no attribute 'AddGlobalIDs_management'      I have no idea what this means or why I'm getting the error. Can I not add GlobalIDs to tables that only have fields and no records? Thanks for any help.


title: move 1 vertex at a time rather than the entire shape

text: Somewhere I have changed settings that allow me to move 1 vertex at a time while editing a shape. It now moves the entire shape or scales the shape! Any suggestions on how to fix this?


title: osm2po service with custom router

text: I want to make custom router which while calculating route costs, would take into account additional features. Here is the use case: Lets imagine paid motorway. Lets suppose there are 2 toll points. When router calculates route, it should, depending on user settings, find another way avoiding toll points, or inform user that his route crosses toll points here and here.. My first idea was to extend EdgeRouter2 class and override calcEdgeCost function, where I could check for given edge id in database and see if it contains toll points and act accordingly, but I want to be able to update my map and it looks like binding toll points table with edges cannot be done reliably neither by segment id nor osm_id (Is OSM_ID also part of osm2po graph?). Now I think it would be better to store toll points using lat,lon coordinates, calculate where those points are road network, and split edges at the given points. Does this new approach make sense? What osm2po's classes should I use to accomplish this? I have a gut feeling that in this case using DefaultRouter2 instead of EdgeRouter2 would be better as it looks like it's more vertex centric - it doesn't have calcEdgeCost but it has (private) calcVertexCost function..


title: Google fusion tables: Fill in lon/lat columns via geocoded address

text: I have a list of 100 companies that I want to geocode. I have all of their addresses but none of their lon/lat coordinates. I have three columns related to location:   1. Address   2. Longitude   3. Latitude Is it possible that when file->geocode is clicked (relating to the address column) that it fills in the longitude and latitude columns with the coordinates? I can then read these columns from there on in avoiding the need to geocode again. Update On further research it's seems I have to store the geocoded result in a database... not fusion tables as it doesn't seem possible with how I want it to work.


title: How can I do this geometric-operation in QGIS?

text: My problem: I have 3 polygons (a, b, c) on layer 1 and two other polygons (1 and 2) on layer 2. Now I want to "crop and melt" (forgive me I don't know the terminus technicus) my three polygons (a,b,c) with the polygons (1 and 2) that (on a new layer) I have two new areas: area I (marked in red in my picture) and area II (marked in blue in my picture). Each of my new areas has two parts. ![enter image description here](http://i.stack.imgur.com/2l4MG.png) Thank you very much!


title: How do I select only those features, that are connected?

text: I have a polygon layer of water reservoirs and a polyline layer of canals. In my data, some canals are originating from reservoirs, and leading to other canals. On the other hand there are some canals that are isolated or not connected to anything. (How this came about is story for another day) I need to delete the isolated canals and those canal networks(i.e. connected set of canals) that are not connected to a reservoir. I assume that i will have to use some kind of network tracing, but I am not sure where to start. I have ArcGIS desktop, along with Network Analyst to solve this task. How would you go about this task?


title: QGIS Python error on OS X while trying to simplify shapefile geometries

text: I am trying to reduce the file size of the shapefile by simplifying geometry in QGiS on OS X. But I get the following error...               An error has occured while executing Python code:          Traceback (most recent call last):       File "/Applications/QGIS.app/Contents/Resources/python/plugins/fTools/tools/doSimplify.py", line 255, in run         self.runSimplify()       File "/Applications/QGIS.app/Contents/Resources/python/plugins/fTools/tools/doSimplify.py", line 280, in runSimplify         shapeFileWriter = QgsVectorFileWriter( self.outputFileName, self.outputEncoding, shapeFields, wkbType, crs )     TypeError: arguments did not match any overloaded call:       QgsVectorFileWriter(QString, QString, dict-of-int-QgsField, QGis.WkbType, QgsCoordinateReferenceSystem, QString driverName="ESRI Shapefile", QStringList datasourceOptions=QStringList(), QStringList layerOptions=QStringList(), QString newFilename=None): argument 2 has unexpected type 'NoneType'       QgsVectorFileWriter(QgsVectorFileWriter): argument 1 has unexpected type 'QString'          Python version:     2.7.1 (r271:86832, Jun 16 2011, 16:59:05)      [GCC 4.2.1 (Based on Apple Inc. build 5658) (LLVM build 2335.15.00)]               QGIS version:     1.8.0-Lisboa Lisboa,           Python path: ['/Applications/QGIS.app/Contents/MacOS/../Resources/python', '/Users/admin/.qgis//python', '/Users/admin/.qgis//python/plugins', '/Applications/QGIS.app/Contents/MacOS/../Resources/python/plugins', '/Library/Frameworks/SQLite3.framework/Versions/3/Python/2.7', '/Library/Python/2.7/site-packages/numpy-override', '/Library/Frameworks/GDAL.framework/Versions/1.9/Python/2.7/site-packages', '/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python27.zip', '/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7', '/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/plat-darwin', '/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/plat-mac', '/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/plat-mac/lib-scriptpackages', '/System/Library/Frameworks/Python.framework/Versions/2.7/Extras/lib/python', '/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/lib-tk', '/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/lib-old', '/System/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/lib-dynload', '/System/Library/Frameworks/Python.framework/Versions/2.7/Extras/lib/python/PyObjC', '/Library/Python/2.7/site-packages', '/Applications/QGIS.app/Contents/Resources/python/plugins/fTools/tools']      Any idea what the issue could be ? Gaurav


title: How to convert MapInfo TAB and/or MIF to SHP using free tools or ArcGIS for Desktop core product?

text: How do I convert MapInfo TAB and/or MIF files to SHP using only free tools or ArcGIS itself? I've not had much luck trying to do this in the past and I don't want to have to shell out for the ArcGIS Data Interoperability extension which seems to be the 'official' way to go about this.


title: Can shapefiles be converted to .TAB files?

text: I read this response to this question, "Yes, if you have access to MapInfo, you can directly read Shapefiles (as .shp is an "open" format) and save it as TAB files!" Once it is saved as a .TAB, are you then free to modify the data in MapInfo?


title: How to change language direction from left to right in arcgis viewer for flex in layer list widget?

text: I have created my web map app using ArcGIS Viewer for Flex and want to change the direction of layer names in the layer list to the Arabic language direction (i.e. right to left; as below) ![my web mapp app](http://imageshack.com/a/img822/3144/huod.png)


title: Why does ArcGIS Publisher appear to mirror ecw Lidar files on southern axis of file?

text: ArcGIS Publisher appears to be mirroring the ecw Lidar files on the southern axis of the file. Has anyone experience this? Is there a setting in ArcGIS Publisher that I need to turn on or off? The way i noticed is by comparing two different imagery data sets.


title: ArcSDE Geodatabase MS SQL Server restore of .bak file

text: If we have ArcSDE 9.3.1 running on MSSQL Server 2008 R2 Enterprise and we do a full backup of a geodatabase. Then I try to restore the .bak file to a MSSQL Server 2008 R2 Express on a machine that doesn't have SDE installed, does this create a proper geodatabase? What about the sde user in the previous geodatabase..does that need to be created prior? We want to do some testing of converting an APDM geodatabase from Oracle to MSSQL Server but would like to do a lot of testing locally on a MSSQL Server Express.


title: Copying shapefiles within windows explorer (or - extension)

text: Is there a user-friendly way to copy shapefiles to a different file location? ArcCatalog is not fast enough for this purpose, creating new map and navigating to this new workspace takes a lot of time.. Does an extension for windows explorer (version 7) exist (maybe specifying an advanced filter) ?


title: How to draw rectangle to clip layers using python addin in ArcGIS 10.1?

text: I want to clip rasters in data view using interactive rectangle (to define Area of Interest) and add them to current dataframe. Following is a piece of code which is not returning anything. Any suggestions to correct it would be appreciative:               import arcpy     import pythonaddins          def __init__(self):         self.enabled = True         self.cursor = 1         self.shape = 'Rectangle'          def onRectangle(self, rectangle_geometry):         """Occurs when the rectangle is drawn and the mouse button is released.         The rectangle is a extent object."""              extent = rectangle_geometry              ras1 = arcpy.Clip_management(r"C:/temp/ras",                                     extent.XMin extent.YMin extent.XMax extent.YMax,                                     ras1, "#", "#", "NONE")         arcpy.RefreshActiveView()         return ras1


title: How to auto-create a schema.ini file for a .csv?

text: I have been given some data in a csv file, which I need to import into my ArcSDE geodatabase. ArcGIS is having problems recognizing the field types, and is not importing some data. It is recognizing some fields as Integer, when it should have been Text. I know that we are supposed to create a schema.ini file, which will contain the appropriate data type, like in this question. My csv file has over 50 fields, so I was wondering if there was a way to auto- create this schema.ini file, and then edit it to the correct data types, instead of creating this file by hand, from scratch.


title: Python Mapping in Matplotlib Cartopy Color One Country

text: I have plotted a map of the world using matplotlib Cartopy. Now I want to select a specific country in the map in this case the United States and change the color. I think this is possible but not sure? I also, don't know how to filter the data for a 'country name' or other data that might be contained in the file.               import matplotlib.pyplot as plt     import cartopy     ax = plt.axes(projection=cartopy.crs.PlateCarree())     ax.add_feature(cartopy.feature.LAND)     ax.add_feature(cartopy.feature.OCEAN)     ax.add_feature(cartopy.feature.COASTLINE)     ax.add_feature(cartopy.feature.BORDERS, linestyle='-', alpha=.5)     ax.add_feature(cartopy.feature.LAKES, alpha=0.95)     #ax.add_feature(cartopy.feature.RIVERS)     ax.set_extent([-150, 60, -25, 60])      ![WorldMap](http://i.stack.imgur.com/6XSE8.png) All I want is two colors one for the United States and one color for all other countries.   Thanks for the help since I'm new to mapping via python.


title: What is the highest resolution of population density data that I can obtain for USA freely?

text: Ok so as far as my understanding goes the most high resolution dataset available that factors in many other things (other than Census data ) at a 1X1 km resolution is Landscan data created by Oak ridge National Laboratory (ORNL). But, these are not available freely. What is the closest alternative you know of? I am planning to use US Census Tiger data for 2010 and the lowest resolution it provides is of Census Tracts which are variable in their size (see here). Can anyone suggest me any other alternative (transformed) datasets ?


title: The topology cannot be validated because it is checked out

text: when I try to validate an existing topologiy in my personal GDB I got this error "the topology cannot be validated because it is checked out". How to solve it? thanks! It is a personal GDB. I'm working in a big MGCP project... I recieved the gdb already with the topology, so perhaps come from Create Replica


Is it possible to make a DEM from contours in QGIS? I have a shp file (contours) and their heights, can I use QGIS for produce a DEM without using GRASS? I find this creating a dem from contours but the answer is for Arcgis.


title: How to generate DEM from contours using QGIS?

text: Is it possible to generate digital elevation model using contour lines in QGIS 2.4 version?


title: How can I get clean road attributes from OSM data?

text: I would like to use osm road data for a project where I need to broadly distinguish between primary, secondary, and tertiary roads. OSM data seems just fine to work with Arcgis' Network Analyst. When I download data for WeoGeoMarket, selecting only highways, I face a lot of different classes, sometimes mis-coded (such as 'teriary' instead of 'tertiary'). However, if I visually check on OSM interactive map, I can see that the legend is not messy at all. My question is therefore the following: is it possible to get OSM road data with clean feature attributes? Thanks for your help!


title: Data visualization and maintenance tool

text: ## Background Looking to leverage an existing GIS system for exploring organizational data. ## Architecture The following figure represents a high-level overview of the system's desired features: ![](http://i.imgur.com/Y01ph.png) The most basic usage would be as follows:   1. The user visits a web site.   2. The system presents a map (having regions, cities, and buildings).   3. The user drills-down on the map to a particular building.   4. The system provides a basic CRUD interface.   5. The user can view and modify information about personnel (e.g., their assigned teams), equipment (e.g., network appliances), applications, and the building itself (e.g., contact and phone numbers). Ideally, all the components should be open-source (or otherwise free). ## Problem This must be a small project that needs a quick (but functional) prototype, mostly to confirm whether or not such a system would be useful in the long term. ## Questions   * What software components would you use to quickly develop a working prototype?   * What open-source solutions already exist, if any? ## Ideas Here is what I am thinking:   * PostGIS \- Define the regions, cities, and sites   * Google Maps \- Display an interactive, clickable map   * geoJSON \- Protocol between PostGIS and Google Maps   * Seam \- CRUD interface ## Custom Development For example, this would entail:   1. Installation and configuration      1. Configure SSH for remote logins     2. Subversion (or git)     3. PostgreSQL     4. PostGIS     5. Java     6. Tomcat     7. Seam     8. JasperReports   2. Enter GIS information into PostGIS   3. Aggregate data sources into PostgreSQL database   4. Develop starting page for map interface   5. Develop clickable Google Maps interface   6. Develop summary reports   7. Develop CRUD interface using Seam for data maintenance Surely something like this already exists? Thank you!


title: How to use Zoom Slider like the older version in Arcgis Javacript API 3.3?

text: I'm using Arcgis Javacript API 3.3. In this version zoom slider look like this. ![enter image description here](http://i.stack.imgur.com/YXZph.png) And here, this is the older version. ![enter image description here](http://i.stack.imgur.com/bbNF9.png) Is there any way to make it like the older version?


title: How to Deal with Overlapping Coordinates in Statistical Analysis

text: I decided to make an archaeological map about certain tombs distributed in the area of Crete and make some statistical analyses using the spatial statistical tools of GIS. ![enter image description here](http://i.stack.imgur.com/PHYhH.jpg) Specifically, I'll be using the nearest neighbor analysis to find out whether they are clustered and a Kolmogorov-Smirnov test to examine their relation to their elevation. I was faced with a dilemma - for several tombs the coordinates were exactly the same. Now this could have happened because they were constructed very close to each other and the GPS did actually capture the same numbers or the archaeologist who did the survey thought he should just take one point for the whole cemetery. How can I deal with this problem? Would it be wrong if I performed the test using more than once the same coordinates? The maximum is three tombs having the same coordinates... Thank you in advance for your answers :)


title: gdal_calc.py outputs huge files

text: I think this is a general enough question to apply to more than just gdal_calc, but perhaps not. When I run               gdal_calc.py -A map.tif --outfile=deforestation_00-10.tif --NoDataValue=0 --calc="A >= 7"      on a 38000 x 38000 62.5MB tif, the output is 1.37GB (and still 38000 x 38000). I feel like i'm probably missing something big here, i.e. how tif data is stored. [the same thing happens when I run it using the raster calculator in qGIS] Thanks in advance. Bonus points--anyone know how to include and/or logic in the --calc="..." in order to evaluate the following: (map@1 >= 7 AND map@1 <= 9) OR ((map@1 >= 10 AND map@1 <= 12) * 2)


title: Breaking a range of values into buckets for a gradient

text: If I have a range of values and want to make a choropleth map with a gradient, what algorithms are out there to determine the value ranges for the different colors in the gradient? For example, if I have a dataset with the values 1,10,40,50,2000,123,23,68,39,21,10,800,600,200 and want to color 14 states based on those values, what is the best way to determine the ranges? Something like [1-50]= light red, [51-100] red, [101-2000] dark red


title: Where does the copied data get placed when publishing a service to ArcGIS Server

text: I use ArcGIS for Server 10.1. What is not clear to me is, where does server save data when I upload it to it? I mean the case, when data in the service is not referenced, but it is inside the service definition file. Is it saved in some folder? Or is it saved in some geodatabase included in the server? Where? Thank you.


title: OpenLayers/GeoServer/Oracle combo input/output scenario

text: Im completely new to GIS. I must finish already started project that uses Oracle as DB for geoobjects, Geoserver, OpenLayers. Task is to draw something on a map (OpenLayers must handle that drawing i assume...) and then save it to Oracle database(as a layer i think). Ive looked through some manuals on geoserver, but maybe you can help me to quickly start in this "simple" input/output scenario?   Thank you.


title: How will I get route along 3 points?

text: Can any one help me out to find route for three points. the points are start- via(Mid)-End. I want route to get create from start point to via(mid) point and via point to end point. Can any one can provide any help to me regarding this issue. using this link:[http://workshop.pgrouting.org/index.html. How it will be done Sql query to retive same from php file in link above


title: How to style features created by DrawFeature control?

text: I've been following this tutorial: http://workshop.pgrouting.org/chapters/geoext_client.html#select-the-start- and-final-destination It contains an Openlayers.Control.DrawFeatures control defined in the following code sample. You can also see the lines where the author comments _"if we want to apply a special style to the start point we should do this here"_. Problem is: I don't know how to apply a style in this setting and cannot find any examples using DrawFeatures control in this way. How can I have the start point use a different style than the end point using this DrawFeatures control?               DrawPoints = OpenLayers.Class(OpenLayers.Control.DrawFeature, {          // this control is active by default     autoActivate: true,          initialize: function(layer, options) {         // only points can be drawn         var handler = OpenLayers.Handler.Point;         OpenLayers.Control.DrawFeature.prototype.initialize.apply(                 this, [layer, handler, options]             );     },          drawFeature: function(geometry) {         OpenLayers.Control.DrawFeature.prototype.drawFeature.apply(                 this, arguments              );         if (this.layer.features.length == 1) {             // we just draw the startpoint             // note: if we want to apply a special style to the              //       start point we should do this here         } else if (this.layer.features.length == 2) {             // we just draw the finalpoint             // note: if we want to apply a special style to the              //       final point we should do this here                  // we have all what we need; we can deactivate ourself.             this.deactivate();                     }     }     });


title: How to find the angle between intersecting features in two featureclasses?

text: I have two intersecting line featureclasses. I want to find the angle at each point of intersection using ArcGIS 10 and Python. Can anyone help?


title: Find projection from dataset

text: I have a dataset that I have pulled from a MaxSea 3d database. I am trying to figure out what projection is used, so I can transform and use the data. I have tried following this guide http://support.esri.com/en/knowledgebase/techarticles/detail/24893, for a small sample of data (300k points). Firstly, the fields are named lat and lon, indicating a Geographic Coordinate System (datum), but the data given looks more like it is in a Projected Coordinate System, becuase the numbers are in 8-9 digits range. The data on file is collected around the Faroe Islands, and if the data is in Geographic Coordinates, it should be in or around 62°N and 7°W ± a few degrees, so I can conclude that the data is not in any format where the leading digits are degrees. (not decimal degrees or DDMMSS.ss, as suggested here). That leaves me with either a proprietary, unknown projection (possible!), or it is in a Projected Coordinate system. I keep on following the Esri guide, to figure out what Projection is used, and this is where I need help, because the Esri guides only covers the U.S. > If the data is in the United States and shows an extent in which the > coordinates to the left of the decimal are 6, 7 or 8 digits, the data is > probably projected to the State Plane or UTM coordinate systems. ![screenshot](http://i.stack.imgur.com/aHMBT.png) Here I see that there are 9 digits (!) left to the decimal point, and when I try to add some reference data, I find that my data is outside the extent of the earth. What would you try, I you were asked to figure out the projection of the data? Here is a sample of the data (50k points). The fields are: id, level, lat, lon, weightdepth_cmap, valuedepth_cmap, weightdepth, valuedepth. The fields of interest are `lat`, `lon`.


title: Can QGIS projects be opened in ArcGIS for Desktop?

text: I am about to start a project with a number of people for a nonprofit organization. I am wanting to create the files in QGIS, but some of the people that I am working with are city officials and might feel more comfortable using ArcGIS for Desktop as opposed to QGIS. If I were to create project files in QGIS could others open up them up in ArcGIS for Desktop? Or is there a way for me to export it to a file that could be opened in ArcGIS for Desktop? Or maybe there is file conversion plug-in that could be used?


title: PostGIS: Find points that lie within a set of coordinates

text: I have a PostGIS db, and I want to find the points that lie within some region (a bounding box). I have two sets of coordinates, from which I can get all four rectangle points that form the box (if needed). My data column in question is named 'point' and it is also of type point.   * Is there anyway to specify four sets of coordinates (lat/long) and get all the points that lie within the box.   * Or specify two points and let the DB work out the rectangle's corners, and return the points within Just in case I am not being clear as to what I want to achieve. The equivalent with 'vanilla' sql if I had a lat and long field instead of a point would be:               SELECT * FROM myTable where lat> xMin AND lat < xMax AND long > yMin and long < yMax      **UPDATED EDIT:** I am trying underdark's solution. At first I didn't have the ST_MakePoint constuctor (now I do) and I still get a very similar error (just on a different character).               SELECT * FROM myTable WHERE ST_Within(ST_MakePoint(point),GeometryFromText('POLYGON((75 20,80 30,90 22,85 10,75 20))',4326))      and I am getting this error:               ERROR:  function st_makepoint(point) does not exist     LINE 1: SELECT * FROM triples WHERE ST_Within(ST_MakePoint(point),Ge...                                               ^     HINT:  No function matches the given name and argument types. You might need to add explicit type casts.               ********** Error **********          ERROR: function st_makepoint(point) does not exist     SQL state: 42883     Hint: No function matches the given name and argument types. You might need to add explicit type casts.     Character: 39      **EDIT:** In the short term I can solve it with:               SELECT * FROM triples WHERE box '((point1),(point2))' @> point      But I will have to work out why none of the PostGIS functions are not working for me.


title: How to generate a world file with rotation

text: I`m writing code to geo reference some images that I get with a rc plane. I have the central coordinate of the image because its atached to the center of the plane. But I have problems when I open the images + jpw in quantum gis, their center doesn't match a point shape that I made with the coordinates that I assume as the center of the plane. My doubt is, the 2 last lines of the world file represent the coordinate of the center of the upper left pixel of the image, but if the plane is heading south when the image has taken this is will be the lower right pixel coordinate? What I'm doing right now is moving the four coners of my image to 0,0, rotating it by my heading, and moving it back, am I doing something wrong? The formulas that i`m using to calculate the first four lines of the world file are those:               A = Cos((PI / 180) * Heading) * XPixSize;     D = Sin((PI / 180) * Heading) * -YPixSize;     B = Sin((PI / 180) * Heading) * -XPixSize;     E = Cos((PI / 180) * p.Heading) * -YPixSize;      ![enter image description here](http://i.stack.imgur.com/Z33zL.png) **EDIT:** I'm changing the aproach, now I'm rotating the image with a program to orientate it to north, but I'm still having problems to georefece it, I have the centroid coordinate in Lat, Lon, both sizes of the original image in meters, and the size in pixel of the new, and the old images, but i'm still having problems to create the worldfile. Just one detail: when I rotate the image, the program creates a image bigger than the original, to let the rotated image fit in. My question is now similar to this one: georeferencing aerial photos with only known centroid


title: How to use user made png/jpeg icons in QGIS 2.0?

text: I can only add svg icons in qgis 2.0. How can i use png or jpeg or gif images as icons in qgis2.0?


title: Two Controls at the same time

text: I have a problem about modifying control for multiple layers . on the map there are tow layers i would like to modify. in addition for the second layer, i have a select control to show a popup i tried:               layer1_control=new OpenLayers.Control.ModifyFeature(layer1, {         id: "ModifyFeatureControl1",         mode: OpenLayers.Control.ModifyFeature.RESIZE | OpenLayers.Control.ModifyFeature.DRAG});          layer2_control= new OpenLayers.Control.ModifyFeature(layer2, {         id: "ModifyFeatureControl2",         mode: OpenLayers.Control.ModifyFeature.DRAG});          layer2.events.on({         featureselected: function(e) {             createPopup(e.feature);         },         featureunselected: function (e) {         }     });      the problem is these controls are not working in same time. only the last added layer control!


title: How to modify the edit widget in ESRI flex viewer 3.0?

text: I am running ArcGIS Server 10.0 (sp2) on Windows server 2008 64-bit with IIS 7. I am using the ESRI flex viewer version 3.0. I have created a feature layer and published it as a rest service and included a template and use the edit widget with this feature service to allow specific users to modify a single attribute on road features on the feature layer. The Edit widget, as is, will allow you to select multiple features (in this case, roads), but you cannot update all selected features at once, you have to do it one at a time. I want to add a button to the window that appears (the attributes window) that will allow you to apply the change to all selected features. I have flex builder 4.6 and have created some custom widgets, but this change is a little beyond my skillset currently. I've looked in the editwidget.mxml and editwidgetattributeinspectorskin.mxml, but don't see where actual edits are made in either of these files. I'm thinking the actual codeing for applying edits to the selected features is in a different section of the source, but don't know where to look. Any help is greatly appreciated.


title: How to align two rasters of exact same cell size, extent in ArcGIS 10?

text: I have two raster datasets. One is already georeferenced. The second one is not georeferenced but should have the exact same extent and cell size as the first one. It is supposed to exactly overlap the first one. Is there an easy way to have the second raster adopt the same georeferencing information as the first one? This is all in ArcGIS 10. Thanks! **EDIT** I should add that the second raster was derived from the first one. I exported the first one and ran some analyses on it using a C++ tool I wrote, now I want to get the output back into ArcGIS. So the two layers will definitely overlap, it's just a matter of giving the second raster the right cell size/extent.


title: Why does this simple Python OGR code create an empty polygon?

text: Please excuse my ignorance...I am new to Python OGR. I am trying to create a simple square polygon and fail...this produces an empty polygon. Where did I go wrong?               from osgeo import ogr     driver = ogr.GetDriverByName('ESRI Shapefile')     datasource = driver.CreateDataSource('c:/temp/testPolygon.shp')     layer = datasource.CreateLayer('layerName',geom_type=ogr.wkbPolygon)          #create ring and polygon object:     myRing = ogr.Geometry(type=ogr.wkbLinearRing)     myRing.AddPoint(0.0, 0.0)     myRing.AddPoint(0.0, 10.0)     myRing.AddPoint(10.0, 10.0)     myRing.AddPoint(0.0, 10.0)     myRing.AddPoint(0.0, 0.0)#close ring     myPoly = ogr.Geometry(type=ogr.wkbPolygon)     myPoly.AddGeometry     #this produces myPoly.GetArea() of zero!          #create feature polygon:     feature = ogr.Feature( layer.GetLayerDefn() )     feature.SetGeometry(myPoly)     layer.CreateFeature(feature)      Even the following returns an area of zero: from osgeo import ogr               #create polygon object:     myRing = ogr.Geometry(type=ogr.wkbLinearRing)     myRing.AddPoint(0.0, 0.0)#LowerLeft     myRing.AddPoint(0.0, 10.0)#UpperLeft     myRing.AddPoint(10.0, 10.0)#UpperRight     myRing.AddPoint(10.0, 0.0)#Lower Right     myRing.AddPoint(0.0, 0.0)#close ring     myPoly = ogr.Geometry(type=ogr.wkbPolygon)     myPoly.AddGeometry     print (myPoly.GetArea() ) #returns zero, not 100...


title: How to transfer data and replica information from a file geodatabase to ArcSDE?

text: I've got an arcgis 10.1 file geodatabase that contains a number of featureclasses and a number of one-way child replicas. Is there any tool that I can use to transfer the data and replica information verbatim into an ArcSDE geodatabase? I'm giving the replica file geodatabase to an external person and I want him to easily be able to move the data and replica information to an ArcSDE database which he will then version and publish as a geodata service that I'll replicate to. The Extract Data Wizard in the distributed data toolbar does half the job by transferring the featureclasses verbatim (including GlobalIDs). But replica information does not appear to be transferred. Create Replica doesn't appear to help in this situation. Any help appreciated. Regards, Dion.


title: Slow performance when selecting features

text: I am trying to use the following code to select the features that have the attribute values that I pass in. The performance is really slow when it comes to a lot of values e.g >16k values to match features. Is there something I'm missing or a better way of doing this? // Assuming values.Any() is true:                    var featureSelection = (layer as IFeatureSelection);          var whereClause = new StringBuilder();     var count = 0;     while (true)     {         var values = valuesToSelect.Skip(count*1000).Take(1000).ToList();         if (!values.Any()) break;         var thousandsClause = new StringBuilder();         foreach (var value in values)         {             thousandsClause.Append(",");             thousandsClause.Append(value);         }         thousandsClause.Remove(0, 1);         thousandsClause.Insert(0, fieldName + " IN (");         thousandsClause.Append(")");         whereClause.Append(" OR " + thousandsClause);              count++;     }     whereClause.Remove(0, 4);     IQueryFilter queryFilter = new QueryFilter();     queryFilter.WhereClause = whereClause.ToString();          featureSelection.SelectFeatures(queryFilter, esriSelectionResultEnum.esriSelectionResultNew, false);     featureSelection.SelectionChanged();           Sample WHERE IN clauses: 761 values to select 122,072 values to select Edit-- So far with recommendations from some of the comments, I've improved the code and the results are faster but still not desirable. Here is the code I've begun trying using the same WHERE IN clause:               IQueryFilter queryFilter = new QueryFilter();     queryFilter.WhereClause = whereClause.ToString();     queryFilter.SubFields = "OBJECTID";          var featureLayer = layer as IFeatureLayer2;     if (featureLayer == null) continue;          IFeatureClass featureClass = featureLayer.FeatureClass;     var dataset = (IDataset) featureClass;     IWorkspace workSpace = dataset.Workspace;     ISelectionSet selectionSet = featureClass     .Select(queryFilter, esriSelectionType.esriSelectionTypeHybrid,esriSelectionOption.esriSelectionOptionNormal,workSpace);          var featureSelection = layer as IFeatureSelection;     if (featureSelection != null)     {         featureSelection.SelectionSet = selectionSet;         featureSelection.SelectionChanged();     }


title: A way to log the output from the Msys terminal in OSGEO4W?

text: Is there a way to log all the text output from commands in the Msys terminal to a text file when using OSGEO4W in Windows? I'm using GDAL commands and specifically wish to record any error messages that come up. Any help much appreciated!


title: Is there a way to copy/paste data & limit the attribution to match the attribution in the Target Data Set?

text: In an effort to avoid digitizing, I want to copy specific watercbodies from a provincial database (shapefile format) and paste it to my taget layer, which are waterbodies significant to my project. I don't want all the columns of attribution to paste over to my target data set (ie - I don't need the "GUT_NUMBER" column in my DBase). Is my best approach to simply delete the attribute columns I don't want when I'm done, or is there a better/more efficient way to do this?


title: Zooming Google Maps to fit the extents of a Fusion Table layer

text: How can I zoom a Google Maps control to fit the contents of a Fusion Table layer? The layer only contains pushpin / marker point data. One method is to perform a REST query on the Fusion Table and to calculate the bounds, and zoom to these (i.e. as illustrated here: http://gmaps- samples.googlecode.com/svn/trunk/fusiontables/fitbounds.html ). However this seems very inefficient - especially for larger datasets, because all of the table data has to be sent to the client and compared. This question appears to ask the same query, but the answer does not apply. The answer simply queries the data for the current map extents. Nice and efficient, but the opposite of what I want. I want to zoom to display all available data.


title: how do I position a popup relative to a point in Openlayers?

text:
I have a fixed size FramedCloud Popup that I want to position relative to the point calling it, always above the point and with its right side 50px to the right of the point. I am looking at the API and see positionBlocks and blocks, but I do not understand how to use them and if that's what I need. I'm really looking for. I tried setting popup.relativePosition, hoping to find an offset parameter somewhere, but relative position doesn't seem to be having any effect.               var vector = new OpenLayers.Layer.Vector("Points",{         eventListeners:{             'featureselected':function(evt){                 var feature = evt.feature;                 var popup = new OpenLayers.Popup.FramedCloud("popup",                     OpenLayers.LonLat.fromString(feature.geometry.toShortString()),                     new OpenLayers.Size(275,71),                     "

" + feature.attributes.title +"

",                     null,                     true                 );                 popup.imageSrc = 'img/popup.png';                 popup.autoSize = false;                 popup.imageSize = new OpenLayers.Size(275,71);                 popup.relativePosition = "bl"                 popup.fixedRelativePosition = true;                 //popup.positionBlocks = {}                 feature.popup = popup;                 map.addPopup(popup);             },             'featureunselected':function(evt){                 var feature = evt.feature;                 map.removePopup(feature.popup);                 feature.popup.destroy();                 feature.popup = null;             }         }     });


title: how to get Coodinates after Resizing a feature Openlayers

text:
I'm trying to display the coordinates of a Feature or vector layer after I use the modify radio button to resize my vector, but it does not work, nor on my drag feature. I'm using `afterFeatureModified()`, maybe this is wrong? It only works with my `onFeatureAdded()` when I finish drawing a circle. Below is my code, Any assistance is much appreciated.                                                                                <link href="css/style.css" rel="stylesheet" type="text/css" />

             <input type="radio" name="type" value="none" id="noneToggle"                    onclick="toggleControl(this);" checked="checked" />                  navigate

             <input type="radio" name="type" value="regular" id="regularToggle" onClick="toggleControl(this);" />             draw regular polygon

             <input type="radio" name="type" value="modify" id="modifyToggle"                    onclick="update();" />                  modify feature

                     <input id="rotate" type="hidden"                             name="rotate" onChange="update()"  checked="checked"/>

                     <input id="resize" type="hidden"                            name="resize" onChange="update()"  checked="checked"/>                   <input id="keepAspectRatio" type="hidden"                            name="keepAspectRatio" onChange="update()" checked="checked" />

                     <input id="drag" type="hidden"                            name="drag" onChange="update()"  checked="checked"/>


title: Is it possible to get the coordinates of a user clicked point in a current MXD with ArcGIS 10.0 ArcPy?

text: I would like to know if it is possible to get the coordinates from a user clicked point in the current map document - arcpy.mapping.MapDocument("CURRENT")?


title: Merge polygons in postgis

text: I have a hand-digitized coverage stored in postgis. There, I have a few shiver-polygons I would like to get rid of by including them into one of the neighbouring polygons. are there any commands that can be used to do that? Something like update coverage set the_geom=st_union(...) ... Of cource, I want to use the attributes from the large polygon - I must be able t decide which way the union goes.


title: Printing with the ArcGIS Javascript API?

text: Any examples of how to print a map using the JSAPI? All from the client?


title: Why is this ExtJS combo box not filtering?

text: I am trying to place a combo-box in the top of my web application which will search an internal WFS from Geoserver. I am using the combo-box so I can have the "type ahead" functionality which I think will help then narrow down the site they are looking for. Then once I solve this issue I will have a listener which will get the info and zoom to the location on selection. I am working of a test at then moment to see if I can implement it in an existing site, but I am running into some trouble. The combo-box doesn't filter the data store, it just returns a full list of the datastore regardless of what i type (see picture at bottom of post). The code I have used is below, unfortunately I don't have an external facing site to show. I originally was hoping to set my WFS output format to JSON and just load that into a JSON data store but I kept getting errors that I was unable to decipher. Thanks Ando               var queryLayer = "dam";     var protocol= OpenLayers.Protocol.WFS({         url: "/geoserver/ows",         featureType: queryLayer,         featureNS: "www.awebsite.com.au", //Just for SE proper featureNS website used in application         version: "1.0.0",         srsName: "EPSG:4283"     });          var assetStore = new GeoExt.data.FeatureStore({         fields: [             {name: "name", type: "string"}             //{name: "type", type: "string"}         ],         proxy: new GeoExt.data.ProtocolProxy({                 protocol: protocol                           })          });          new Ext.Viewport({         layout: "border",         items: [{             region: "center",             id: "mappanel",             //title: "Map",             xtype: "gx_mappanel",             map: map,             tbar:['->',{                     xtype: 'combo',                     fieldLabel: 'Search for an asset:',                     store: assetStore,                     displayField: 'name',                     forceSelection: true,                     triggerAction: 'all',                     minChars: 1,                     displayField: 'name'                                                                         }],             layers: [layer],             extent: extent,             split: true         }]     });      ![enter image description here](http://i.stack.imgur.com/qxKeQ.png) _Edit_ This is the query that is getting sent to the server, its not passing the name from the combo box to the query so I assume that why everything is being returned. Is there a configOption I am missing here??               <wfs:GetFeature xmlns:wfs="http://www.opengis.net/wfs" service="WFS" version="1.0.0" xsi:schemaLocation="http://www.opengis.net/wfs http://schemas.opengis.net/wfs/1.0.0/WFS-transaction.xsd" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"><wfs:Query typeName="feature:dam" xmlns:feature="www.mywebsite.com.au"/></wfs:GetFeature>


title: Many of the Esri samples fail the WC3 HTML validator. Is this a problem?

text: I'm running into problems with some ArcGIS Server JavaScript code (page elements are positioned different in IE than FF/Chrome), so I checked it against the HTML validators at http://validator.w3.org and http://html5.validator.nu Many errors were found, so I returned to the official Esri JS samples to ensure I was using valid code. Surprisingly, I found that many of the samples also failed the validators. Examples of errors include:               Error: Attribute label not allowed on element button at this point.     <button id="dropdownButton" iconClass="bingIcon" label="Base" data-dojo-type="dijit.form.DropDownButton">      from the Basemap sample. Is it a problem that the official Esri samples fail the validators? Is this related to the use of DOJO, and therefore not an actual problem?


title: Examples which use map reduce to process GIS data

text: Have anyone a view examples how to use mapreduce effective with GIS data?


title: Best program for crowd-sourced polygons on Gmaps

text: Looking for a way to crowd-source polygon creation on GMaps (in a wiki environment) that doesn't involved ESRI's expensive server software. Requirements:   \- a bunch of different wiki pages can be updated with a map showing a single polygon. \- all the polygons over the entire site are combined into one master map \- preferably everything in KML so Fusion Tables can be involved for styling, but not necessary Any ideas?


title: Do OpenLayers maps always show/flash an empty screen while zooming?

text: When double-clicking on an OpenLayers map, the screen is momentarily blank while the next zoom level draws. Is this avoidable? It's not ideal, given that many other mapping APIs show a smooth zooming effect. Thanks


title: Updates to ArcSDE via Oracle trigger

text: Is it possible and/or feasible for an Oracle trigger to periodically update ArcSDE via a REST geoprocessing service? Does anyone have any experience with such a thing? The reason for this is that I have an Oracle database that is separate from the SDE database, and there are number of business rules that we want to run first, before inserting the geometry into SDE. So the basic workflow would be:   1\. Point features via ArcGIS Server via ArcSDE loaded onto web map   2\. User moves point or edits attributes   3\. Edits to points first get written to Oracle, then, via trigger to ArcSDE


title: join grid or polygon layers in R

text: I have two grids. Cells in each grid have an attribute. Size of grid cells differs between the two grids and the two grids do not overlay each other perfectly. I would like to create a new layer, or group of polygons, or grid, that results from joining the two existing grids. Each cell in the new grid should retain the attribute from each of the original grids. Area not occurring in both original grids can be discarded from the new grid. Such discarded areas would only have one attribute (and one missing observation for the second attribute). How can I do this is R? So far I have tried `over` in the `sp` package and `unionSpatialPolygons`. I have also considered the `PBSmapping` package, which allows for creating intersections of polygons. This post describes how to join polygons within a layer, but perhaps not how to join polygons in different layers: joining polygons I have written code to compute the (x,y) coordinates of the four corners of each cell that results from joining the two original grids. That code is at the very bottom of this post. Perhaps I can create a new polygon layer from the resulting data frame. However, I suspect that process has already been implemented in a package. Below is code to create the two original grids, assign an attribute to each cell in each grid and plot the two grids. There will be 16 cells in the new grid. Those cells generally will not be squares, but rather rectangles of various shapes and sizes. Thank you for any advice.               library(rgdal)     library(maptools)     library(raster)     library(rgeos)     library(gridExtra)          setwd('c:/users/mmiller21/simple R programs')          set.seed(1234)          # create grid 1          grd1 <- GridTopology(c(2,2), c(1,1), c(4,4))     polys1 <- as(grd1, "SpatialPolygons")          # assign projection to grid          proj4string(polys1) = CRS("+proj=longlat +datum=NAD83 +no_defs +ellps=GRS80 +towgs84=0,0,0")          # create fake atttribute data for each grid cell          poly.data1 = data.frame(f1=runif(length(row.names(polys1)), 1, 10))     row.names(poly.data1) <- paste0('g', 1:length(row.names(polys1)))          # convert grid to a SpatialPolygonsDataFrame:          poly.df1 = SpatialPolygonsDataFrame(polys1, poly.data1)          # create grid 2          grd2 <- GridTopology(c(3.8,3.8), c(2,2), c(2,2))     polys2 <- as(grd2, "SpatialPolygons")          # assign projection to grid          proj4string(polys2) = CRS("+proj=longlat +datum=NAD83 +no_defs +ellps=GRS80 +towgs84=0,0,0")          # create fake atttribute data for each grid cell          poly.data2 = data.frame(f2=runif(length(row.names(polys2)), 100, 200))     row.names(poly.data2) <- paste0('g', 1:length(row.names(polys2)))          # convert grid to a SpatialPolygonsDataFrame:          poly.df2 = SpatialPolygonsDataFrame(polys2, poly.data2)          jpeg(filename = "joined.grid.jpeg")          plot(poly.df2, col='blue', xlim=c(1,8), ylim=c(1,8))     text(coordinates(poly.df2), labels=row.names(poly.df2))     axis(1, lty = 0, lwd = 0)     axis(2, lty = 0, lwd = 0)     plot(poly.df1, add=TRUE)     text(coordinates(poly.df1), labels=row.names(poly.df1))          dev.off()      ![enter image description here](http://i.stack.imgur.com/pxOye.jpg) Here is code to compute the (x,y) coordinates of each cell in the new grid:               x1 <- seq(1.5,5.5,1)     y1 <- seq(1.5,5.5,1)          x2 <- seq(2.8,6.8,2)     y2 <- seq(2.8,6.8,2)          x <- c(x1,x2)     x <- sort(x)     x     # [1] 1.0 2.0 2.8 3.0 3.8 4.0 4.8 5.0 5.8 6.8          new.x <- x[x>= min(x2) & x <= max(x1)]          y <- c(y1,y2)     y <- sort(y)     # [1] 1.0 2.0 2.8 3.0 3.8 4.0 4.8 5.0 5.8 6.8          new.y <- y[y>= min(y2) & y <= max(y1)]          z <- merge(new.x, new.y)          plot(z$x, z$y)


title: How to handle line data in Mapinfo?

text: I have laoded my **lines data** in MS-SQL in one table. I have **points data** in DBMS as well in separate table . When I open both table the map window would show lines and points as I want but now when I join more points the data is updated to **points table**. And since point table is table to store points only and line table is to store line data only. How is it possible to solve this issue so that when I connect points Line table is updated? > I am using > >   1. Microsoft SQL 2008 >   2. Connect the dots Tool to connects points ( actually I need line data in > line table) >   3. Mapinfo 10.0 >


title: Is there an open source GIS to view and edit CityGML models?

text: I'm looking for an open source GIS that supports viewing and editing of CityGML models. So far, I've found a list of open source applications that support CityGML, but they all seem to be stand-alone applications: http://www.citygmlwiki.org/index.php/Open_Source. I would also be interested in your experience with CityGML in a GIS (also commercial) context.


title: How to use the ISpatialFilter interface with ArcGIS 10.1?

text: From what I've heard till now: 1) ArcGIS from version 10.1 and forth will only support http connections: local connections will not be supported any longer. 2) The code               ESRI.ArcGIS.Carto.IMapServer iMapServer = (ESRI.ArcGIS.Carto.IMapServer)ctx.ServerObject;      or the code               ESRI.ArcGIS.ADF.Web.DataSources.ArcGISServer.MapResourceLocal mfunct = (ESRI.ArcGIS.ADF.Web.DataSources.ArcGISServer.MapResourceLocal)gisResource;                     ESRI.ArcGIS.Carto.IMapServer iMapServer = (ESRI.ArcGIS.Carto.IMapServer)mfunct.MapServer;      can't both work with an internet service connection because they require a local connection. 3) The ESRI.ArcGIS.ADF.Web.SpatialFilter object cannot be set to identify what kind of query I'd like to perform. I mean, with the ESRI.ArcGIS.Geodatabase.ISpatialFilter interface you can do something like               spatialFilter.SpatialRel = ESRI.ArcGIS.Geodatabase.esriSpatialRelEnum.esriSpatialRelIntersects;      to specify what you want as result. SpatialFilter hasn't this feature. I need the ISpatialFilter spatial relationship feature. 4) If you put it all together, it figures out that you need an IMapServer- implementing object to perform the ISpatialFilter query               ESRI.ArcGIS.Geodatabase.IFIDSet fidSet = iMapServer.QueryFeatureIDs(mapDescription.Name, layerID, queryFilter);      So if you read back points 1 and 2... what am I supposed to do? I've also asked this on the official forums but still no answers. Please help because I'm getting crazy on this.


title: Applying osm2pgsql .diff patch in windows 7

text: When attempting to run osm2pgsql in command prompt, I ran into the password prompt as described here. This requires installing a patch. I'm unfamiliar on how to install a .diff patch on Windows 7. Do I need to install GnuWin as shown here? Thanks for all the help!       **EDIT:** This is the same question as this one on Stack Overflow.


title: Why do some maps show Sandy Island?

text: I need to ask you guys after reading this article. Any speculations on why some maps showed it? What can be the cause of inaccuracies in (web)Maps?


title: From pgRouting shortest_path back to osm node id

text: I'm successfully performing path finding queries like `shortest_path` on data imported with `osm2pgrouting`. The tables I have inside my pgrouting DB are the following:                              List of relations      Schema |        Name         |   Type   | Owner      --------+---------------------+----------+-------      public | classes             | table    | user      public | geography_columns   | view     | user      public | geometry_columns    | table    | user      public | nodes               | table    | user      public | spatial_ref_sys     | table    | user      public | types               | table    | user      public | vertices_tmp        | table    | user      public | vertices_tmp_id_seq | sequence | user      public | ways                | table    | user      I get the way of a given path between a source and dest node this way:               SELECT * FROM ways         JOIN          (SELECT * FROM             shortest_path('SELECT gid AS id, source::int4 AS source, target::int4 AS target, length::double precision AS cost FROM ways',3,7,false,false))        AS route ON ways.gid=route.edge_id;      ways's table has the following structure:                Table "public.ways"         Column    |       Type       | Modifiers      --------------+------------------+-----------      gid          | integer          | not null      class_id     | integer          |       length       | double precision |       name         | character(200)   |       x1           | double precision |       y1           | double precision |       x2           | double precision |       y2           | double precision |       reverse_cost | double precision |       rule         | text             |       to_cost      | double precision |       the_geom     | geometry         |       source       | integer          |       target       | integer          |       Now `gid`, `source` and `target` values aren't original osm node's ids (the node's ids contained in the `.osm` file from which the table has been filled). I suppose that `pgRouting` create custom ids in order to index ways and edges. From what I can understand the original ids of nodes contained in the .`osm` files are the ones listed into nodes table. I would like, given a row of ways table, recover the original osm id of source and target. How can I do that?is it possible?


title: Are there any common rules for WMS clients to comsume the WMS in tiles?

text: Looking at OpenLayers and other clients. They requests small 256x256 tiles of the map. I am wondering if there is a common ruleset for these tiles to relax the server computation on raster layers. Lets say a client need to show the map at BBOX, do it just generate a BBOX that fit exactually what it wants or do it ajust the box to fit some cache able bounding box that servers may implement? My consern is that creating a WMS server wont scale well because clients can ask it to generate images as it fits. If the server site could just serve some cached tiles then it would scale much better. I have seen there is something called WMTS, but our clients is asking for a WMS, so we cant just give them a WMTS.


title: ArcGIS 9.3.1 User Privileges

text: I know that you can grant privileges for users using `ArcCatalog` for the whole `Dataset` as letting them select, update, delete and insert feature classes. But my question is how could I grant the privileges for each feature class in the `Dataset` for users. Example :               Users     Feature Class     _____     _____________      1             Ftr 1      2             Ftr 2      User 1 can read and edit on Ftr1 but can only read Ftr2, and User 2 can read and edit on Ftr2 but can only read Ftr1. While Ftr1 and Ftr2 are within the same `Dataset`


title: ArcGIS key numbering labels for every feature in one layer

text: I would like to use key numbering when labeling specific layer. But the tricky thing is I want to label every feature in this way.   By default, key numbering is last option in Strategy order (source). And there's no possibility to change it.      Is there any workaround to force ArcMap to use key numbering for every feature?


title: How to get the absolute raster row / column of rastercell when using ArcObject Raster Cursor?

text: When using the ArcObject RasterCursor (see snippet below), how do I get the absolute row / column of the current pixel value?               IRasterCursor rasterCursor = raster.CreateCursorEx(null);          do     {         IPixelBlock3 pixelblock = rasterCursor.PixelBlock as IPixelBlock3;              Array pixels = (Array)pixelblock.PixelData[rasterBandIndex];              for (int i = 0; i < pixelblock.Width; i++)         {             for (int j = 0; j < pixelblock.Height; j++)             {                 // How to get absolute row / column of this pixelvalue                 dynamic pixelValue = pixels.GetValue(i, j);             }         }     }     while (rasterCursor.Next());


title: Bing Open layer shows ok in QGIS and Composer Preview, but is misaligned after exporting to PDF

text: I have QGIS 2.2 running on a Mac with iOS Maverick. I have prepared a map with layers with different CRS, including a Bing Open Layer as base layer. The project has the 'On the fly coordinate correction' option activated, and it works ok, although sometimes things go out of position, but nothing a refresh from a pan doesn't solve. Paraphrasing part of question #51986: "Sometimes when changing the scale (scrolling in and out) the [Bing] layer updates (...) wrong in (...) position compared to the polygon layer. But as soon as I "pan" the map everything reconstructs right. Now the problem I have concerns the print composer." "The result of any kind of export - pdf, or image or svg - of my map (vektor layer + [Bing] layer in the background) differs from what is shown in the preview of the print composer: the vector layer and the google layer positions are mismatching (...)" Could this be related to On the fly correction? Any ideas? Thanks! ![Screen capture from QGIS 2.2, with layers correctly aligned](http://i.stack.imgur.com/Iv4rO.png) ![Screen capture from PDF created with QGIS Composer, showing the misaligned Bing layer](http://i.stack.imgur.com/zfuqY.png)


title: Is there a purge function to remove layers from a qgs file which are not present in the Layers pane?

text: I have some QGIS projects that I've been using for long time, adding and removing layers, updating reference to deleted or moved files etc. Opening with a text editor the qgs file I see that it contains numerous references to files, and this references are not showed on the Layer tree pane. I have here different cases: references to existing file, references to missing files, references to database connections. I got aware f this when I opened my project on another machine, and QGIS reported a missing database connection to vector tables not present in the layer tree view. Manually deleting references from the qgs file solved the problem, but I'm managing huge projects and it's near to impossible to manually check the project consistency. I've tried the "Remove empty layers" plugin, but the ones I'm trying to remove are not empty layers. Is there any way to 'purge' qgis file?


title: Merge two geodatabases with identical structure

text: I'm using an ArcGIS addon which stores all created data (shapes and tables) in a gdb. The addon is made for single user use but to safe time I want to split the work with someone else. After that I want to merge both geodatabases into one. Is this possible? Line segments get IDs which are referenced in the tables. So by merging, identical IDs would have to be changed - in all tables in the same way. Maybe by adding to each ID in one of geodatabases the same number, like ID=ID+1000. But I don't know how to do that in ArcGIS either. Also it seems that the merge function in the toolbox only accepts one file add a time. So I would have to merge each table seperately because I don't want to change the structure. The number of tables has to stay the same, otherwise the addon probably won't work anymore. Is there an easy way to do this. Is this possible at all? Robert


title: Error in viewshed (arcgis)

text: I have a point data as observer point with field attributes for Spot, offsetA, offsetB, Radius1, Radius2, Vert1, Vert2, Azimuth1 and Azimuth 2. For the raster file- a dsm with elevation data. They have the same projection and coordinate system But when I run the viewshed analysis I get this error: ERROR 010067: Error in executing grid expression. No valid Observation points present. Failed to execute (Viewshed). Any suggestions how to solve this?


title: How to query and display layer using ArcGIS Python script tool?

text: I have database which contains point shapefiles and raster files for 3 species. My aim is to develop a script tool that should have a search tab to enter the species name. also, a check box to select what layers to display. how to develop this. give me hint to develop python script for querying and display a layer. also, i should add the script tool as a button in the toolbar.


title: How to place map legends and titles in Geoserver?

text: I am new to geoserver and tried to learn from basics & familiar with GIS concepts. I published a shape file in geoserver, Now I want to add a title & legend to the published file. Please advise me the solution.


title: Point of Intersection for a Ray and Earth's Surface

text: Say I have a ray vector originating from somewhere out in space. An example might be a ray of light from the sun. How can I calculate the point of intersection (if it exists) between the ray and the Earth's surface? I'm using Cartesian coordinates (ECEF) and I'd like the Earth's ellipsoid Geometry to factor into the calculation.


title: Is Django's LayerMapping capable of putting a static value into all features?

text: I'm using GeoDjango's LayerMapping utility to try and import a set of features. The catch is that aside from the attributes in the source layer, I want all the destination features to have the same value of another attribute. Specifically:               > sourceLayer = theDataSource[5]     > sourceLayer.fields     ['DataSources_ID', 'Source', 'Notes', 'Shape']     > from django.contrib.gis.utils import mapping     > fieldMapping = mapping(theDataSource, layer_key=5)     {'datasources_id': 'DataSources_ID', 'source': 'Source', 'notes': 'Notes', 'shape': 'MULTIPOLYGON' }      The destination model has another field, 'owningmap', which is an integer field. I would like to load data from the OGR source via...               > lm = LayerMapping(DestinationModel, theDataSource, fieldMapping, layer=5)     > lm.save()      ... but I would like every feature to have the same integer value put in to the 'owningmap' field. As though I could somehow add a "static" mapping to fieldMapping...               > fieldMapping['owningmap'] = 2      Is this possible? Otherwise my last resort is to load via LayerMapping, then calculate the value for all the just-loaded features in a separate database transaction... Thanks!


title: How to get X Y coordinates and cell value of each pixel in a raster using Python?

text: I’m really new to Python and I would like to know whether there is a quick method to get cell values of a raster pixel by pixel and the coordinates (map X Y coordinate of centre of each pixel) using Python in ArcGIS 10? To describe this further, I need to get the map X, map Y and cell value of the first pixel and assign those three values in to three variables and repeat this step for rest of the other pixels (loop through the entire raster). Thanks. * * * I think I need to describe my question more. Problem is, I need to get the X Y location of a pixel of the first raster and get cell values of several other rasters corresponding to that X Y location. This process should be loop through every pixel of the first raster without creating any intermediate point shapefile as it is going to be really really time consuming since I have to handle a raster with nearly 8 billion pixels. Also, I need to do this usin Python in ArcGIS 10. @JamesS: Thank you very much for your suggestion. Yes this would work for one raster but I need to collect the cell values for several other rasters as well. The problem is, after getting the X and Y coordinate of the first pixel of the first raster, I need to get the cell value of the second raster corresponding to that X, Y location of the first raster, then third raster and so on. So, I think when looping through the first raster, getting X and Y location of a pixel and getting cell values of the other raster corresponding to that location should be done simultaneously but I’m not sure. This can be done by converting the first raster in to a point shapefile and performing Extract multivalues to point function in ArcGIS 10 but I’m unable to follow that method because I don’t want to create a shapfile as it is going to be really slow and I HAVE TO FIND A SOLUTION USING PYTHON, not with any existing tool in ArcGIS. @hmfly: Thanks, Yes this method (RastertoNumpyarray) will work if I can get the coordinate of a known row and column value of the array. @whuber: I don’t want to perform any calculations, all I need to do is write X Y coordinates and cell values in to a text file and that’s all


title: WFS getFeature Request Blank XML for some layers

text: I have following url request to geoserver. the first request results blank page but the second page results in number of feature count. Why it there blank xml for some layer while other gives number of features.                http://192.168.70.65:7070/geoserver/wfs?request=GetFeature&typename=GisData:sites&version=1.1.1&resultType=results          http://192.168.70.65:7070/geoserver/wfs?request=GetFeature&typename=NCellGisData:sites&version=1.1.1&resultType=hits


title: Best Platform to develop QGIS Python Plugins

text: I am currently developing QGIS Python Plugins in a Windows 7 64-bit environment. Having said this I am stumbling across various issues which from what I have been able to gather have a less tendency of happening in the other platforms (i.e. Linux and OSX Mac). I would like to ask the community to answer the following question: **What is the best platform to develop QGIS Python Plugins?** Thanks


title: Select by Attribute on Date field using wildcards?

text: I have a date field within a file geodatabase (10.1). An example of a value in this field: **19/11/2012 08:49:56** I want to select all rows where the time is between 9pm to 11pm AND 6am to 9am. I don't care what day, month, year it is. So I want to make use of wildcards to ignore these values. What would be the correct statement to put in the Select By Attributes to achieve this?


title: Is it possible to add leaders to point symbols if they overlap?

text: I've added symbols to one of my point shapefiles, but for several datapoints this has resulted in symbols stacked on top of each other because there are multiple sample types for one sample location (i.e. water, sediment, and benthic samples collected at one station). The result is a cluttered mess! Can I bring a few of my symbols away from the location and have leaders point to the sampling station? Thanks! Update: Thank you for the info, guys! I'm working with ArcView so I can't use the "Disperse Markers" tool. If I download QGIS and deal with my markers in that environment, will my changes "stick" if I bring the layer back into ArcView? (I haven't used QGIS before; is this even possible?) Update October 19: Thank you for the posts. The Maplex extention suggestion doesn't seem to apply to my situation as I am attempting to remove overlap from symbols, not text. I am tring to determine if I can bring some of my symbols (I have 20 different symbols in the layer in question) out on leaders to aviod clutter. Is this even possible? I don't mind putting in a extra time to "fiddle" with my symbols. Thanks!


title: Identifying topological relationships using PostGIS

text: I'm trying to perform a basic typology analysis using PostGIS. My objective is to find any polygons that touch other polygons. To do this, I thought the `ST_GetFaceEdges` would work (reference). I want to check every polygon in my database and list all the other polygons that touch it. In the image below I would expect the result to say that two of the polygons (which are buildings) touch one building, and the results of the other 4 to say that they touch 0 polygons. ![enter image description here](http://i.stack.imgur.com/Y0VYn.jpg) However, I'm having some difficulty understanding what to do. When I tried to copy the example there were a few parts of it that I did not understand.               -- Returns the sequence, edge id     -- , and geometry of the edges that bound face 1     -- If you just need geom and seq, can use ST_GetFaceGeometry     SELECT t.seq, t.edge, geom     FROM topology.ST_GetFaceEdges('tt',1) As t(seq,edge)         INNER JOIN tt.edge AS e ON abs(t.edge) = e.edge_id;      I'm not sure whether `topology` is the name of a table, column or a part of the function. I assumed that it was the table, but I'm not sure.               -- try out on sample of data     -- Table is called 'TestArea', column used is 'fid_1'     SELECT t.seq, t.edge, geom     FROM  TestArea.ST_GetFaceEdges('fid_1', 1) As t(seq,edge)         INNER JOIN tt.edge AS e ON abs(t.edge) = e.edge_id;      I'm also not sure what the function of the inner join is - does this join the result to the original object?


title: Should a GDAL Plugin be the same version as GDAL aready built into QGIS

text: Should a GDAL Plugin be the same version as the GDAL that's aready built into the QGIS version installed? Or,can another different version of GDAL be built and it installed as a Plugin for QGIS, even though the version of GDAL is different? Is it reading QGIS's GDAL installation or a standalone GDAL installation for the plugin? I'm going to try and create a Ubuntu QGIS Lisboa plugin with File GDB support. And, I'd like to do the same for Windows 7 QGIS Lisboa. It's going to be installed using a standard osgeo4w (rather than advanced) install, but I'd like to build a plugin for it. The reason I'd like to use standard install is that I'd like to have all the standard Python, PostGIS and Spatialite packages installed as well, and not sure how to do that with "Advanced", and still get GDB support. If I can't get QGIS to use a new GDAL plugin featuring GDB support, the I'd like to atleast be able to add another standalone version of GDAL, and use the command line ogr2ogr to convert File GDB to shapefiles and dbfs (for the additional domain tables). Thanks in advance for your advice.


title: File Geodatabase XY Resolution and XY Tolerance

text: Here - http://gisintelligence.wordpress.com/2011/04/19/limiting-xyz-geometry- drift/ I found the description of XY Resolution and XY Tolerance in "Arcgis World". I have a question - my cadastral data in shapefile format have accuracy 1 cm. So, if I want to work with my data in Arcgis File Geodatabase I should use XY Tolerance 0,01 m and XY resolution 0,001 m ( http://goo.gl/uqDrG   * is that ok? I mean - that values are proper for topology matters and so on?  Thanks for answer


title: Mapnik TextSymbolizer value as integer

text: I'm working with TextSymbolizer in Mapnik. I have a shape file with a attribute which is a **real**. I wish to print the value of this attribute as **integer** and I don't find how I can do it in Mapnik. I must not change the source shape file. Can someone advise of the recommended way to do this?


title: Diferent Line color and Diferent Markers at the same time

text: i'm about to develop a GIS Track App, but until now, i'm testing with some codes, that i've been creating, but for my project we are going to use XML as data records, today, i am able to create a lines and Markers, but i want to achieve the creation of Line and Markers depending each vehicle, e.g: Vehicle 1:( Line Color=Blue and Marker=Lotus.ico) Vehicle 2:(Line Color=Red and Marker= Viper.ico) here is my XML code: IdDispositivo= IdCar , IconoDispositivo=Icon and LineaColor=LineColor of every Car               <?xml version="1.0"?>     <t_QryHistorial>     <Historial>         <IdDispositivo>1</IdDispositivo>         <IconoDispositivo>resources/imagenes/lotus.ico</IconoDispositivo>         <Fecha>2013-03-31</Fecha>         <Hora>09:41:00</Hora>         <Longitud>-101.36206</Longitud>         <Latitud>19.76807</Latitud>         <LineaColor> #FF0000 </LineaColor>      </Historial>     <Historial>         <IdDispositivo>1</IdDispositivo>         <IconoDispositivo>resources/imagenes/lotus.ico</IconoDispositivo>         <Fecha>2013-03-31</Fecha>         <Hora>10:41:00</Hora>         <Longitud>-101.31812</Longitud>         <Latitud>19.73511</Latitud>         <LineaColor> #FF0000 </LineaColor>     </Historial>     <Historial>         <IdDispositivo>1</IdDispositivo>         <IconoDispositivo>resources/imagenes/lotus.ico</IconoDispositivo>         <Fecha>2013-03-31</Fecha>         <Hora>10:45:00</Hora>         <Longitud>-101.26318</Longitud>         <Latitud>19.74610</Latitud>         <LineaColor> #FF0000 </LineaColor>     </Historial>     <Historial>         <IdDispositivo>1</IdDispositivo>         <IconoDispositivo>resources/imagenes/lotus.ico</IconoDispositivo>         <Fecha>2013-03-31</Fecha>         <Hora>10:50:00</Hora>         <Longitud>-101.25220</Longitud>         <Latitud>19.79004</Latitud>         <LineaColor> #FF0000 </LineaColor>     </Historial>     <Historial>         <IdDispositivo>1</IdDispositivo>         <IconoDispositivo>resources/imagenes/lotus.ico</IconoDispositivo>         <Fecha>2013-03-31</Fecha>         <Hora>11:00:00</Hora>         <Longitud>-101.21924</Longitud>         <Latitud>19.80124</Latitud>         <LineaColor> #FF0000 </LineaColor>     </Historial>     <Historial>         <IdDispositivo>1</IdDispositivo>         <IconoDispositivo>resources/imagenes/lotus.ico</IconoDispositivo>         <Fecha>2013-03-31</Fecha>         <Hora>11:20:00</Hora>         <Longitud>-101.15331</Longitud>         <Latitud>19.88569</Latitud>         <LineaColor> #FF0000 </LineaColor>     </Historial>     <Historial>         <IdDispositivo>1</IdDispositivo>         <IconoDispositivo>resources/imagenes/lotus.ico</IconoDispositivo>         <Fecha>2013-03-31</Fecha>         <Hora>11:40:00</Hora>         <Longitud>-101.17331</Longitud>         <Latitud>19.88669</Latitud>         <LineaColor> #FF0000 </LineaColor>     </Historial>     <Historial>         <IdDispositivo>1</IdDispositivo>         <IconoDispositivo>resources/imagenes/lotus.ico</IconoDispositivo>         <Fecha>2013-03-31</Fecha>         <Hora>12:40:00</Hora>         <Longitud>-101.18331</Longitud>         <Latitud>19.90669</Latitud>         <LineaColor> #FF0000 </LineaColor>     </Historial>     <Historial>         <IdDispositivo>1</IdDispositivo>         <IconoDispositivo>resources/imagenes/lotus.ico</IconoDispositivo>         <Fecha>2013-03-31</Fecha>         <Hora>12:42:00</Hora>         <Longitud>-102.17331</Longitud>         <Latitud>20.88669</Latitud>         <LineaColor> #FF0000 </LineaColor>     </Historial>     <Historial>         <IdDispositivo>2</IdDispositivo>         <IconoDispositivo>resources/imagenes/viper.ico</IconoDispositivo>         <Fecha>2013-03-31</Fecha>         <Hora>12:42:00</Hora>         <Longitud>-102.17331</Longitud>         <Latitud>20.88669</Latitud>         <LineaColor> #3333FF </LineaColor>     </Historial>     <Historial>         <IdDispositivo>2</IdDispositivo>         <IconoDispositivo>resources/imagenes/viper.ico</IconoDispositivo>         <Fecha>2013-03-31</Fecha>         <Hora>12:42:00</Hora>         <Longitud>-102.17331</Longitud>         <Latitud>20.88669</Latitud>         <LineaColor> #3333FF </LineaColor>     </Historial>     <Historial>         <IdDispositivo>2</IdDispositivo>         <IconoDispositivo>resources/imagenes/viper.ico</IconoDispositivo>         <Fecha>2013-03-31</Fecha>         <Hora>12:42:00</Hora>         <Longitud>-102.17331</Longitud>         <Latitud>20.88669</Latitud>         <LineaColor> #3333FF </LineaColor>     </Historial>     <Historial>         <IdDispositivo>2</IdDispositivo>         <IconoDispositivo>resources/imagenes/viper.ico</IconoDispositivo>         <Fecha>2013-03-31</Fecha>         <Hora>12:42:00</Hora>         <Longitud>-102.17331</Longitud>         <Latitud>20.88669</Latitud>         <LineaColor> #3333FF </LineaColor>     </Historial>      The XML This is a history of each route of each vehicle here is the image, what i want to do it: ![enter image description here](http://i.stack.imgur.com/lmgHZ.jpg) I'm Working with: OpenLayers and ExtJs 4 (GeoExt 2) Do you have some example how to achieve what i wants? Thanks to all. * * * Thanks for your answer, i'm working with Openlayers and Extjs 4, and here is my image, about what i want: imgur.com/ovVU3SA Finally I choose to do what I need using the KML format, and I can already get two separate lines, but I can not even get the icon, line color and line width, any idea that I can be I doing wrong? Here is my code:               var layers = [];                var mapPanel=this.getRefgridMapaJson();                  epsg4326 =  new OpenLayers.Projection("EPSG:4326");                 projectTo = new OpenLayers.Projection("EPSG:900913");                 var kmlStyleMap = new OpenLayers.StyleMap({'default':{                         strokeColor: "${color}",                          strokeOpacity: 0.5,                         strokeWidth: "${width}"                     }});                          var vecLayer = new OpenLayers.Layer.Vector("vector",{                              strategies: [new OpenLayers.Strategy.Fixed()],                              protocol: new OpenLayers.Protocol.HTTP({                                     url: "data/rutas.kml",                                     format: new OpenLayers.Format.KML({                                         formatOptions: {                                           extractStyles: true,                                            extractAttributes: true,                                           maxDepth: 2,                                           internalProjection: projectTo,                                           externalProjection: epsg4326                                         }                                     })                                 }),                     //styleMap: kmlStyleMap                                                 });                layers.push(vecLayer);                mapPanel.map.addLayers(layers);            If i enabled the line: //styleMap: kmlStyleMap , my map does not show anything. Here is my KML code:                <?xml version="1.0" encoding="UTF-8"?>      <kml xmlns="http://earth.google.com/kml/2.0">      <Style id="route_1">      <IconStyle>       <Icon>       <href>resources/imagenes/lotus.ico</href>        </Icon>       </IconStyle>       <LineStyle>         <color>FF0000</color>         <width>5</width>        </LineStyle>      </Style>       <Placemark id="route_1">        <name>West Campus</name>        <description>         <![CDATA[Prueba ruta 1.]]>         </description>        <styleUrl>#route_1</styleUrl>        <LineString>         <coordinates>         -101.36206,19.76807          -101.31812,19.73511         -101.26318,19.74610         -101.25220,19.79004         -101.21924,19.80124         -101.15331,19.88569         -101.17331,19.88669         -101.18331,19.90669         -102.17331,20.88669         </coordinates>           </LineString>          </Placemark>       <Style id="route_2">          <IconStyle>         <Icon>          <href>resources/imagenes/viper.ico</href>         </Icon>        </IconStyle>        <LineStyle>        <color>3333FF</color>         <width>5</width>        </LineStyle>        </Style>        <Placemark id="route_2">        <name>East Campus</name>        <description>           <![CDATA[Prueba Ruta 2]]>         </description>         <styleUrl>#route_2</styleUrl>        <LineString>         <coordinates>         -103.17331,20.88669         -103.19331,21.88669         -103.20331,22.88669         -103.121331,23.88669               </coordinates>        </LineString>        </Placemark>      Thanks again.


title: Converting altitude in feet to pixels on a mercator map image

text: I know the question might seem farfetched given that the map is 2d. But I essentially want to project items-aircraft-above the map and I would like their altitude to be vaguely 'realistic' in terms of map size etc. The whole thing is 3D so it is possible to see the altitude. I have the altitude in feet. Is it simply a case of world diameter / image width = number of feet represented in a pixel. altitude in feet / number of feet represented in a pixe I get the feeling that would be too simple. Thanks


title: Set a user scale to use QGIS for surfaces documentation

text: I would like to use QGIS for the documentation of mural paintings (so I don't need georeferecing) I would like to set my own scale from a raster file (measure a know distance and use it as reference). I'm reading the manual but I don't find how to do that... Thanks a lot!


title: How to import values from feature to variable list in ModelBuilder?

text: I'm trying to create a model to find POI by selection by atrribute. I have feature layer (POI with three fields Kat_1, Kat_2, Kat_3 ) and 3 parameters ( Variables: Kat_1,Kat_2,Kat_3). How do I import value list from feature to variables?


title: ArcGIS Kernel Density with polyline, search radius / bandwidth calculation

text: I'm conducting a density analysis of a series of polylines and I'm hung up on determining an appropriate search radius. The default calculation performed by arcgis is to take the shorter dimension of the input feature extent divided by 30. However, I would prefer to use the formula set forth by Bailey and Gatrell in Interactive Spatial Data Analysis (1995). They suggest that the bandwidth for a kernel analysis be set to 0.68n^(-0.2) where n=the number of features being analyzed. This calculation is straightforward for points as each point is the same size. However, the polyline dataset I am working with has lines of varying length. Simply taking the number of lines as n doesn't seem appropriate. I have considered converting the polylines to points but before I begin with this approach I am curious if there might be other valid approaches? Thanks in advance for any help you can provide.


title: Where do I find TerraGO GeoPDF hands-on tutorial online?

text: Where can I find a TerraGO GeoPDF hands-on online tutorial? I'd like some examples I can use on my computer to learning how to use it. My agency has their toolbar on our computer system.