title: How do I upload a file geodatabase to fme server 2013?

text: I periodically load an SDE-geodatabase with about 1 GB of data that is delivered to me in a ArcGIS file geodatabase (.gdb) format. In order to run that job on FME server I need to upload the .gdb to the FME server. I tried to upload it in various ways but since the .gdb is a directory containing lots of files in Windows, the upload dialog box does not accept the directory, instead it wants at file. I also tried to zip the .gdb. and then I got a success but no data was read from the .gdb. This is a transcript from the log file: 2014-05-02 12:29:24| 1.9| 1.6|WARN |Generic Reader: Unable to guess a valid reader name from dataset 'D:\FME_TEMP\fmeservertst_10_0.gdb.zip_1399026562740_14984'. A NULL reader will be used instead. 2014-05-02 12:29:24| 1.9| 0.0|INFORM|Generic Reader: No format specific settings have been found for the NULL Reader -- default values as documented in the Reader/Writer manual will be used Can anyone put me on the right track? I use FME Server 2013 and ArcGIS 10.0


title: Merge uniquely named point shapefiles in 50,000 different folders into a single layer

text: I have point shapefiles in 50,000 different folders (each shapefile file name is unique) - about 33,400 shapefiles (not all folders will have a point shapefile). I like to open the point shapefile in each of the folder into QGIS Desktop 2.2.0 map frame and then hopefully run a merge on them to get a single point shapefile layer. I tried drag-dropping from Windows Explorer to QGIS, and QGIS froze initially and continue to add the point shapefiles for several hours and didn't seem to stop. I found the Load Them All plugin (http://goo.gl/HbJhhS) that seem to specify the feature type and also specify a string filter, this plugin also didn't seem to work... I was thinking that had I planned well, I could've automated and copied the point shapefiles into another single folder as soon as they were created and I could 've dragged-dropped into mapframe from this single folder. I don't know if this would've had worked. My analysis run took 18 days to run, so I don't want to re-run just to do that. I tried copy-pasting with Windows Explorer and it hung, may be because it was 100,000+ files... Would you know any trick to do this?


title: Changing the color ramp in arcgis javascript

text: Using the arcgis javascript api, I am displaying geotiffs. What I want to do is change the color ramp in these images. I've looked through the api and I have not seen anything. Does anyone know where I can get some examples of this?


title: Using conditional statements in Field Calculator

text: I have two fields, "P_Area" and "PC_Area", whose sizes I'd like to compare in a third field, "Areas_Equal". The "Areas_Equal" field should contain a "Yes" if the other two are equal and a "No" if they aren't. Could someone please tell me how I can do this with Field Calculator? In Python or VB. I've tried a both based on snippets I've found online, but I'm doing something wrong because I keep getting an error message.


title: Arc 10 Python .shp attributes - usings fields to populate other fields - if/then/else statement

text: Ive been looking for a solution and unable to find one so please help if able: I'm trying to avoid manually inputting Begin and End engineering station values already inputted into a seperate field. Thanks for any help on this. I'm dealing with(4) fields: Begin_Sta, End_Sta, Beg_End and Station. The Begin_Sta and End_Sta fields have engineering stations. The Beg_End field has either "BEGIN" or "END" inputted. The Station field is the one I want to populate using python in the field calculator. If Beg_End has "BEGIN" then the engineering station value in the Begin_Sta field is returned into the Station field, else the End_Sta value is used to populate the Station field. Is this possible using python in the Field Calculator?


title: Can I update features in an SDE database from python without using an ArcEditor License?

text: I have a non spatial data set that I would like to keep sync'ed in SDE as a spatial set. So I would like to be able to insert, update, and delete without using one of my two ArcEditor licenses since this will be running on a server. Can I do this? Clarification: I meant data set in the most general term, just a raw table in a DBMS that holds spatial data, but is not a part of any ESRI product. Full explanation: We use a very old MapObjects customization that was purchased from a 3rd party vendor. It stores binary shape information (conforming to ESRI shapefile standards thankfully) in a binary column in MSSQL server. Every time an edit is made to a record in that datbase a trigger records the edit. As of now, there is an application that runs on a server that checks out an ESRI ArcEditor license and processes the edits (with ArcObjects) that have been made and makes the appropriate change in an SDE database with feature datasets and feature collections. I would like to avoid this use of a license on an unattended server for what would seem to me to be a normal use of SDE and not something inherent that requires an ArcEditor instance.


title: Appropriate map projection for the UK?

text: As the title: what is the most appropriate projection to reduce distortion, if I'm just showing the UK on a map? Is WGS84 fine (since it's a relatively small area), or should I use some other kind of projection?


title: Problem in Adding line feature in Openlayers vector layer

text: I have added the WMS line feature in vector layer. But I am facing some display problem in it as shown below screenshot.Is there anyother way to do this? ![enter image description here](http://i.stack.imgur.com/yZit0.png) But there is no issue in point features. ![enter image description here](http://i.stack.imgur.com/s0i6E.png) Here is my code:                selectionLayer = new OpenLayers.Layer.Vector("Selection", { styleMap:                 new OpenLayers.Style(OpenLayers.Feature.Vector.style["select"]),         displayInLayerSwitcher: false,         isBaseLayer: false     });          osMap.addLayers([selectionLayer]);          var wmsFeatInfo = new OpenLayers.Control.WMSGetFeatureInfo({     url: 'http://xxx.xxx.x.xx:8080/geoserver/cite/wms',         title: 'identify features on click',                 queryVisible: true     });     wmsFeatInfo.infoFormat = 'application/vnd.ogc.gml';     wmsFeatInfo.events.register("getfeatureinfo", this, getWmsFeatAttributes);     osMap.addControl(wmsFeatInfo);     wmsFeatInfo.activate();          function getWmsFeatAttributes(e) {         selectFeatLst = {};         if (e.features.length > 0) {                    for (var ii = 0; ii < e.features.length; ii++) {                            var layer_ID = e.features[ii].fid.split('.');                            var featId = layer_ID[1];                            selectFeatLst[featId] = e.features[ii];             }                 }      }           function test(key){          selectionLayer.addFeatures([selectFeatLst[key]]); //Here line feature is adding       }      Please suggest me to fix this issue.


title: Finding a certain length of a line from a central point

text: I am not exactly sure on how to describe what I am trying to do so hopefully this picture will help. What I am trying to figure out: I have a point and from that point I have, walking paths from that point. Say I want to know how far a half mile is on the paths, Is there a a method or tool that will aid me in this? By just doing a 1/2 mile buffer on the point does not achieve what I want since the paths are not straight lines. I hope what I trying to do is clear enough![enter image description here](http://i.stack.imgur.com/4p3FR.png)


title: What are hardware requirements for a cloud hosted PostGIS & GeoServer?

text: I have been playing with a local installation of PostGIS & GeoServer. I would like to move my work online now to share it with several collaborators. The proof of the concept map will consist of simple WMS and WFS services mapped through OpenLayers and GeoExt. It will be based on a map of ~3000 administrative units. Initially it needs to serve several people only. From then we will be working on polishing final view and functionality of the application. There have been several suggestions made under answers to this question. For example, Nicklas is running his site on 512MB Linode (I think...). I was thinking about using Ubuntu Server on cloud VPS, most likely hosted on rackspace. My question: what would be the minimal server size I would need to use to start playing with this data online? There are few of them available. Once we decide to open the service to public rackspace promises seamless way of scaling up the service, so I would like to know what bare minimum RAM I need to start with. * * * _**Update 1_** : Similar discussion on GeoServer forum. _**Update 2_** : Useful tutorials on setting up stack on AWS micro instance using Ubuntu and CentOS.


title: What do we have to set in mapfile´s data parameter for SOS?

text: We are trying to run the MapServer SOS example at this link. However, we don´t know how to configure the data parameter in the mapfile. **What do we have to set in mapfile´s data parameter?** We need to relate the PostGis database with MapServer, and we don´t how. Any help with that? Thanks! * * * Here are some examples on the MapServer configuration with a SOS server. Thanks!


title: Picture Labels with leader line from point without annotation?

text: The closest thing I have is using a line call out in the text background, then placing the picture in the line call out border. But the picture tiles throughout. Is what I am trying to do even possible without adding graphics or annotation? Basically I want to do what Google earth does automatically with the placement of their labels? Any extensions besides Maplex? Any tools or solutions? Thanks


title: Using the C++ GDAL API how do I set the GDAL_DATA path?

text: If I want to change the directory of the GDAL_DATA path or point to a different version of a csv file for EPSG numbers how do I do that using the C++ API?


title: Design a routing system

text: I am trying to develop an application, firstly web, later mobile, that builds a route based on points of interest, basically it should show the route and the points of interest and when the user clicks it will display some kind of information but I'm really struggling starting to create this...I already have my POI spatial database and I haven't really tried anything, because I don't know where to start. Is there any example you can give me? Thanks in advance.


title: ArcPy Field Calculation with special characters

text: I'm trying to calculate a string field, using "µSv/h". Running this in the manual field calculator works just fine, but I've tried this in the Python window in ArcMap as well as the IDLE, but both fail. For the line below it says invalid syntax.               arcpy.CalculateField_management(tblTempAssayImportName,"Dose_units",""""µSv/h"""","PYTHON_9.3")      When I try to remove the additional three quotes, it comes up with a 000539 error. I've tried using `r` to escape, The line above comes straight from the snippet after running the field calculator manually. Not sure why this isn't working, but I'm assuming it's the special characters.


title: Python: how to speed up spatial search (nearest point)?

text: I'm working on a program that determines the closest location from a given point. The point cloud I'm testing against is very big (~800.000). However, this doesn't really explain why my implementation is so slow. This is my approach: First, I created a spatial index for the point shape `pntshp.ExecuteSQL('CREATE SPATIAL INDEX ON %s' % table_name)` I defined an array of buffer distances to narrow down the search radius. Which of course also means that I have to create a buffer for each point (which might be expensive). `BUFFER_DISTANCES = ( 0.001, 0.005, 0.01, 0.02, 0.05 ) # 100m, 500m, 1km, 2km, 5km` Then, the buffer is used as a spatial filter `node_lyr.SetSpatialFilter(buff)` If the filter returns `None` the buffer distance will be increased.               for buffer_d in BUFFER_DISTANCES:         buffr = get_buffer(xy_street,buffer_d)         ...      Then I am calculating the distance to the points returned by the spatial filter               p=ogr.Geometry(ogr.wkbPoint)     p.AddPoint(xy[0],xy[1])          for feat in node_lyr:         geom = feat.GetGeometryRef()         d = p.Distance(geom)         dist.append(d)             To get the closest point:               def get_closest_pnt(dist, node, how_close):         mrg = zip(dist,node)         mrg.sort(key=lambda t: t[0])         try:             return mrg[how_close]         except IndexError, ierr:             print '%s \ndist/node tuple contain %s' % (ierr,mrg)      It all works fine but is really slow. Creating a spatial index didn't show any effect, really. To calculate 100 points this implementations takes ~6,7 seconds. The program needs to be able to calculate the closest location for more than 2000 points as fast as possible. Any ideas on how to improve my approach? **EDIT** I tried different approaches to see where it gets me. I came across something very astonishing I want to share here. I implemented a simple lookup algorithm as described here, and one of the solutions that where suggested (the sorted set approach). The surprising fact is that performance is not only dependent on the implementation but even more so of the OSX. My original ogr/buffer algorithm turns out to be blazing fast on my OSX whereas it is painstaking slow on Linux (hence the question here). Here are my results (100 runs).               Method       |     OSX        |  Linux Ubuntu     ogr buffer   | 0:00:01.434389 | 0:01:08.384309          sub string   | 0:00:19.714432 | 0:00:10.048649          sorted set   | 0:00:01.239999 | 0:00:00.600773               Specs Mac OSX     Processor 4x2.5 GHz       Memory 8 GB 1600 MHz          Specs Dell Linux Ubuntu     Processor 8x3.4GHz     Memory 7.8 GB      If someone can explain why these differences occur, please don't hesitate.


title: The PrintingTools service produces garbled or pixelated legend patches in PDF

text: I'm building a web map with ArcGIS viewer for flex 3.6. I made some template for print task in different size A3, A4, A5. THe print task works well for A5 and A4 but it produces garbled or pixelated legend patch in output . I looked on ESRI support site and found this http://support.esri.com/de/knowledgebase/techarticles/detail/40538 I asked the server admin to perform the workaround for me but it still doesn't work. Is there any other way I can solve this problem??


title: Get the outer ring inner ring for a given polygon using ogr python binding

text: I am trying to read the geometry of some shapefiles, now I meet some problem when I try to identify the polygon with holes, since a polygon may(or not) have one hole(or multiple holes). This is the code I am using:               def _get_features(coord,datalayer):         bbox=_get_bbox(coord)         datalayer.SetSpatialFilter(bbox)              layer_def=datalayer.GetLayerDefn()         geom_type=layer_def.GetGeomType()              for feature in datalayer:             geometry = feature.geometry().Clone()             if not geometry.Intersect(bbox):                 continue                     geometry = geometry.Intersection(bbox)             fixed_geom=_get_coordinate(geometry)                  #maybe a polygon with hole             if isinstance(fixed_geom[0], list):                 if geom_type is wktPolygon:                     # now, how to get the rings? outer, inner, multiple inner?                              def _get_coordinate(geometry,result=None):         if result is None:             result=[]              if geometry.GetPointCount()>0:             for i in range(0,geometry.GetPointCount()):                 point=geometry.GetPoint(i)[:2]                 result.extend(map(_fix_coordinate,point))         elif geometry.GetGeometryCount()>0:             for i in range(0,geometry.GetGeometryCount()):                 geom=geometry.GetGeometryRef(i)                 _get_coordinate(geom,result)              return result      Any suggestion? BTW, is it possible a hole inside a hole of a polygon?


title: Change kml tag structure

text: I am working on google-maps kml. I have several shapefiles and using QGIS i have converted it to kml files. When i export the kml (save as kml in Qgis). It contains tag like   `<SimpleData name="Name">7</SimpleData>` but i want tag structure like   `<Data name="name"> <value>7</value> </Data>` I have export the shape as kml and update the kml change style add BalloonStyle for that i need to modify structure otherwise it will not display info in popup window.Here is the link for rules of tag structure to display popup in google map as per kml extended data      If we use `<SimpleData>` tag then popup will not display lables in google map      ![wrong popup](http://i.stack.imgur.com/bzNQi.png)   When i change the structure of tags to `<Data><Value>`it will works gr8 ,Expected popup window is ![expected popup](http://i.stack.imgur.com/y2uhW.png)   I just want to change the tag structure is there automated way to change the tag structure from kml file. like below `<simppledata>` to `<data><value>`               <SchemaData schemaUrl="#new_circle">                     <SimpleData name="Name">7</SimpleData>                     <SimpleData name="gid">1</SimpleData>                     <SimpleData name="active_sub">4596</SimpleData>                     <SimpleData name="net_additi">5181</SimpleData>                     <SimpleData name="gross_addi">3917</SimpleData>                     <SimpleData name="churn">4617</SimpleData>                 </SchemaData>`       to               <Data name="name">                         <value>7</value>                     </Data>                     <Data name="active_sub">                         <value>4596</value>                     </Data>                     <Data name="net_additi">                         <value>5181</value>                     </Data>                     <Data name="gross_additi">                         <value>3917</value>                     </Data>                     <Data name="churn">                         <value>4617</value>                     </Data>      here is the sample link in which i have commented the simpledata tag and add manually data and value tag.expected kml output To do it manually will take maximum time because i have thousands of records in other layers.   Please give me any suggestion. Thank you so much in advance.


title: Trying to build pyramids for large raster

text: I finally was able to make a slope with a rather large 5m DEM. But now I it usually wont finish drawing. When ArcMap asks to 'build pyramids' it always freezes sometimes instantly after accepting sometimes around 18%. Whether it is in Catalog or ArcMap 10.0 I cannot get it to finish. And of course even though it was frozen it did start making pyramids. So when I open the file later it shows up as a big black box. I tried building pyramids with 0 layers to remove them, but that doesn't seem to help. So the past couple times I've had to delete that copy and make another copy of the original. Is the file just too big? Has anyone else experienced this? I am using Arc 10 SP 5. Thanks for any insight.


title: OpenLayers World coverage

text: I gave a demo this afternoon and it couldn't have gone any worse :( I checked everything was working before I started (apart from the external internet connection).... My dependency on an external WMS layer of the world coastline cause my app to hang so I would like to serve out my own world coverage rather than rely on any external resource but where do I get the data and the sld from, I have a shapefile of the world coastline but it doesn't look half as good as the default openlayers basemap. My base layer is configured like so:               //create the basemap     basemap = new OpenLayers.Layer.WMS(          "Baselayer", "http://vmap0.tiles.osgeo.org/wms/vmap0?",          {layers: 'basic'},          {             singleTile: true,              ratio: 1,              isBaseLayer: true         }     );      I would like to be totally independent of external resources (vmap0.tiles.osgeo.org above) and rather serve this data out myself via geoserver (wms). Thanks for any help!


title: Worldwide population density data not grouped by country?

text: Are there any data sources available, preferably freely (but not necessarily), which show population density without regard to country or constituency boundaries of any sort? I'd like to analyze a map which shows the human concentration as accurately as possible, without "blurring" it by country or other similarly large areas, and instead using a grid of some sort.


title: Is ther open source tool for calculating Getis-O-ord (ARCGIS) and able to customize in own platform

text: i would like to know is there any open source tool for calculating Getis-O- ord(ArcGIS)? If so can i customize for own application or can i write my script within that application like qgis.


title: Convert ArcMap label expressions from VBScript to Javascript or Python

text: I have a large number of layers (hundreds) with label expressions written in VBScript. I am using ArcMap 9.3.1, and I'm now starting to export the layers as MPKs for use in our ArcGIS Runtime (WPF) application. I get an error that says that the layers can not be exported to MPK with Runtime support because it has label expressions written in VBScript. It will take to long, and too much manpower to manualy convert the label expressions from VBscript to Python or Javascript. I'm looking for a tool to do this automatically. **Another possible solution can be to enable the export of layers with VBScript expressions to MPK format. _This solution would be the best case scenario!_** Any help will be very much appriciated!


title: Save Visible WMS layer to file

text: Is there tool to save currently visible WMS layer to file in QGIS, or some place in QGIS where you can get needed xml to fetch it using gdal ?


title: How do I access a multivalue output parameter from a python script in modelbuilder in ArcGIS 10.0?

text: How do I access a multivalue output parameter from a python script in modelbuilder, in ArcGIS 10.0? That is, I have a python script that loops through several input layers (working through a script tool in a toolbox), performing a number of tasks and finally outputting a final layer. There is one output layer for each input layer. The number of input layers is typically 5, but not always. I have seen this link: How to use multivalue output parameter with ArcPy My entire script is too long to post, so I have posted a pared down version at the end of my post. You will see that my script currently creates output parameters (feature layers) as part of the loop, using an incrementing number variable as the output parameter number. The script works fine as long as I know how many input layers I will be looping through and I can set up the output parameters of the attached toolbox tool appropriately. But I need to have a tool(s) that don't know how many inputs there will be. I see from the referenced link how to set up a semicolon-delimited string as the output parameter, but I can't figure out how to do anything with it in modelbuilder. I also setup a boolean output parameter that simply shows that the script has run, in case the number of output layers is fewer than what I have set up for the tool--but I would like something better than that. I thought all the new iterators and other new modelbuilder-specific tools were supposed to be for things like this, but I can't get it to work. Of course, I could be setting up the tool parameters incorrectly, using the modelbuilder iterators incorrectly, I don't know. Does anybody know how to do this? Can you show me an example? The link I referred to doesn't actually show any examples of the modelbuilder as recipient of such an output parameter (the semicolon-delimited string) or describe what to do. Any help will be greatly appreciated!               # IMPORT MODULES       import arcpy, os, sys, string, traceback, datetime            # SCRIPT ARGUMENTS       # Select layer       inSelect = arcpy.GetParameterAsText(0)       # Input layers       inLayers = string.split(arcpy.GetParameterAsText(1), ";")       # Output workspace       outFolder = arcpy.GetParameterAsText(2)            # OTHER FRONT STUFF          try:         # SET VAR FOR INITIAL OUTPUT PARAMETER NUMBER         outputParam = 4         # SET VAR FOR COMPLETION OUTPUT PARAMETER         bComplete = False         # CREATE LIST FOR INTERMEDIATE DATA TO BE DELETED DURING FINALLY         intermediateData = []         # LOOP THROUGH INPUT LAYERS.         for lyr in inLayers:                  # SELECT FEATURES                  # PROJECT                  # BUFFER                  # MERGE                  # CREATE LAYER             # Create virtual layer             layerName = lyr + prjLayernameSuffix             virtualLayer = arcpy.MakeFeatureLayer_management(outFC, layerName)             # Create layer on disk - required for arcpy.mapping.Layer             layerFilename = layerName + ".lyr"             layerPath = os.path.join(outFolder, layerFilename)             arcpy.SaveToLayerFile_management(virtualLayer, layerPath, "RELATIVE")             # Find layer to symbolize from - some of this is left out from actual script             symbolLayerName = hcatype + "_sample_buffers.lyr"             symbolLayer = arcpy.mapping.Layer(os.path.join(tooldataPath, symbolLayerName))             # Create map layer from layer file - must use complete path             mapLayer = arcpy.mapping.Layer(layerPath)             # Symbolize map layer             arcpy.mapping.UpdateLayer(df, mapLayer, symbolLayer, True)             mapLayer.transparency = outputLayerTransparency             # Add map layer to current map             arcpy.mapping.AddLayer(df, mapLayer)             # Save a layer file with FC in case you want to re-use the layer             arcpy.SaveToLayerFile_management(mapLayer, layerPath, "RELATIVE")                  # CLEAR SELECTION                  # SET OUTPUT PARAMETER AND INCREMENT PARAMETER NUMBER             arcpy.SetParameter(outputParam, mapLayer)             outputParam += 1              # SET COMPLETION OUTPUT PARAMETER         bComplete = True         arcpy.SetParameter(3, bComplete)              # REFRESH         arcpy.RefreshActiveView()         arcpy.RefreshTOC()         mxd.save()          except:         # SET COMPLETION OUTPUT PARAMETER         arcpy.SetParameter(3, bComplete)         tb = sys.exc_info()[2]         tbinfo = traceback.format_tb(tb)[0]          pymsg = tbinfo + "\n" + str(sys.exc_type)+ ": " + str(sys.exc_value)          arcpy.AddError("Python Messages: " + pymsg + " GP Messages: " + arcpy.GetMessages(2))          finally:         # Delete vars and intermediate data         for dataObj in intermediateData:             if arcpy.Exists(dataObj):                 arcpy.Delete_management(dataObj)         del inSelect, inLayers, prjSpatRef, outFolder         # Other appropriate deletions are left out for clarity of example         del layerName, virtualLayer, layerFilename, layerPath, mapLayer


title: Using "Time" in PostgreSQL/PostGIS

text: I'm loading gtfs into a PostGIS database but I have a problem when loading my stop times. The data is in time format "00:00:00" and for 1:30 AM, my data is "25:30:00". I'd like to keep it that way because usually, night buses services start before midnight so I can keep track of the schedule by time, but PostgreSQL doesn't accept time after 24:00:00. I could modify my script (python) so when it comes to treating stop_times.txt, to modify the time to create something like 2530 so I'll have to deal with integer, but I'm not the only one that will have to deal with that db and my colleague would really prefer to deal with time instead of integers... Does anyone could think of a way to "bypass" this limitation? Thanks!!


title: Point Density and area

text: I would like to create circles around a point density layer that reflect the magnitude of the point density. I have created a point density estimation for a few towns. Each town varies in population, although the area does not change among them. I would like to create an "area of influence" around in town that reflects the population density. (Larger circles around higher populated towns.) Is there a way to do this? I am using ArcGIS 10.1


title: offline editing not working on qgis mobile

text: I would like to work with offline editing using Qgis desktop and Qgis Mobile, in the first essay i tried to use Qgis 2.2 and its offline editing plugin , it has worked fine then i moved the qgis project and its sqlite database to the mobile, when i open it with qgis mobile (1.9) it seems it doesnt read sqlite database and give me this error > There are still 1 unhandled layer that will be lost if you closed now i tried to open the sqlite database only,but it gives me error > Failure exploring tables from: "filepath" no such column type i tried to create offline project using Qgis destok 1.8 but it gives this error when i try to export. > Impossible to open spatialite database Any hints or tricks, to use Qgis desktop and Qgis mobile for offline editing ?


title: Control geoprocessing results related to default data values calculated for null inputs

text: I am looking for information on Default data values for the different data types in ArcGIS. We are running a series of Overlay Route Events (Linear Referencing) tool operations. If any of the values from the two input source tables are null then the values in result table are some default data value. For example for Date field types that have null values in the input tables the values that getting populated in the Overlay operation result are '12/30/1899 0:00:00'. This is happening in a file geodatabase at ArcMap 10.1 and Windows Server 2008. Does anyone have information on default data values for data types specifically in reference to geoprocessing? Any suggestions to get the null values to stay when the Overlay does it's calc? Field properties say allow null is yes. I received word from Esri that this is a known bug in the tool **NIM081760**. I'd be interested in hearing of any workarounds.


title: How to toggle ArcGIS webmap layers using ArcGIS Javascript API

text: I've created a webmap using ArcGIS Online that I want to render in a web app using the ArcGIS for Javascript API. In ArcGIS Online, I've saved the webmap with all of the layers off (visibility: false). I want set the visibility of one layer ("Layer X" below) to true using the API.                   arcgisUtils.createMap(mapid, "map").then(function(response){                   layers = response.itemInfo.itemData.operationalLayers;                   dojo.forEach(layers, function(layer){                       if (layer.title == 'Layer X') {                      layer.setVisibility(true)**;                  }                  console.log(layer.title + ", visibility: " + layer.visibility);                   });               });      **I've tried using .show() as well. Using the code snippet above, the visibility of "Layer X" doesn't change. The console log reports the layer title & visibility of all the layers except the one I'm trying to show. The response object returned from the .createMap method is a Deferred object. Is it possible to toggle the layers of a webmap using the ArcGIS for Javascript API, or is the response object immutable?


title: Postgis: Traslate geometry from EPSG:23030 to EPSG:21818, then to lat long

text: I'm from Colombia, South America. I'm a computer programming student and I have to develop a prototype for a GIS system for my small home town using Postgis >= 2.0 as database back-end. The local public works office gave me the old Arcgis shape files: one for urban roads, one for rural roads, one for mountain streams, one for lakes, and so on. All shape files contain the *.shp, *.shx and *.dbf files, except the *.proj file. This is the problem: The local goverment made these shape files long time ago using the EPSG:23030 reference system which locates in Western Europe (A very bad mistake with they still work and are not interested to change due to time and budget resources). However, Colombia is located within the EPSG:21818 Colombia Mainland. Using http //epsg.io/, curiously, I figure that coordinates x,y (in meters since EPSG:23030 is UTM) of the roads, streams and lakes correspond to locations in Nigeria, Africa. You can find a point located on the central park of my home town http //epsg.io/ using these coordinates: (-10153166.72, 1702714.02). So my question is: How can I use Postgis, for example, using the st_transform function, to translate the imported shape files data from EPSG:23030 to EPSG:21818? How can I "convert" the EPSG:21818 to lat, lon coordinates later? Is it possible to do all this with Postgis? Thank you for your replies. * * * UPDATE 2014-07-08 20:16 UTC-5 Thank you, mkennedy for your reply. I uploaded the shape files (no *.proj available), sql files (with EPSG:23030) and some pictures showing a sample point in UTM: 1106289.573,1131273.240 which I've taken as an example (it's an upper left corner on the central park). Since, the original shape files didn't have the projection file, I set it in Postgis with this command:               SELECT UpdateGeometrySRID('viasurb02','geom',23030);      When using QGis we can see some rural roads distort due to low bad quality when capturing the data, so I considered the urban roads layer (viasurb02.shp) for taking the sample point. Huh, That point is located in Benin, Africa. www.mediafire.com/?s1s6ae8bh8i3kev * * * UPDATE 2014-07-08 21:14 UTC-5 Thanks you Paul Ramsey for your replay. First, I've tried your suggestion about importing the shape files by setting the EPSG:23030 SRID:                shp2pgsql -s 23030 -W LATIN1 -D -I viasurb02.shp viasurb02_b | psql -d mycity -U me      Then updated the SRID of the table:               ALTER TABLE viasurb02      ALTER COLUMN geom      TYPE Geometry(MultiLinestring, 21818)      USING ST_Transform(geom, 21818);      However, the traslated geometry is impossible for EPSG:4326 because it is inf/inf ![The same sample point with EPSG:21818](http://i.stack.imgur.com/w0fwC.png) Then, I've tried your second suggestion assuming that the shape files already has the EPSG:21818:               shp2pgsql -s 23030 -W LATIN1 -D -I viasurb02.shp viasurb02_b | psql -d mycity -U me      And now, the same sample point is located in North Venezuela i58.tinypic.com/2jdny2x.jpg. I don't know if the local public works office created in the late 1980s, using a bad projection or using Arcgis as it just it were Autocad. Are there additional steps I can take to make another traslation geometry? Thank you for your reply.


title: Creating training areas for parametric classifcation

text: I am creating a signature file for a supervised classification for a subset of a LANDSAT image. I know what classes I would like to create, and have identified ideal areas for training. However, I am running into a problem in some of these areas in which the covariance matrix cannot be inverted, thus I am not able to run a parametric classification, such as Maximum Likelihood. I understand that the covariance matrix must be linear for this to happen and is a function of the pixels themselves and I cannot edit this. I am using layer- stacked LANDSAT 8 imagery with 14 bands (7 from spring imagery, 7 from winter imagery to detect leaf-on and leaf-off deciduous areas). Is there a way to determine areas that would lend themselves to fit into a linear covariance matrix? I ran across literature that suggests that an increase in bands makes it more difficult to achieve this linear satisfaction for covariance matrix inversion. I am using ERDAS 2014. It has the option to run a non-parametric (such as parallelepiped) along with a parametric (max likelihood), but I would like to stick to maximum likelihood. Would larger training areas fix this? Or more smaller areas merged together? It seems that when adding new training areas it is random as to whether the cov matrix will be invertible or not. Thanks in advance!


title: How to remove the description information when exporting to KML?

text: I have been working in ArcGIS to convert shapefiles into kml files using ArcGIS toolbox. The resulted kml/kmz contain some extra set of information in the description field that is used to show a pop-up screen when clicked on a record. This is really unwanted for me and I am required to delete this html links. I open each of the records and remove all html code from their Description field. Does anyone know whether there is an automatic way/tool to do the same, that can remove the code from all the records at a single cycle?


title: OpenLayers.Layer.MapServer not rendering but new OpenLayers.Layer.WMS does

text: I got the following code:               var layern2 = new OpenLayers.Layer.MapServer("Tech",                     "ServerPathTo_mapserv.exe",                     {                         map: "ServerPathTo_My.map"                     },                     {                         isBaseLayer: false,                         layer: "Layer1",                         version: "1.1.1",                         transparent: "true",                         visibility: true                     });      That does not render a non transparent layer without images (all titles are rendered with a missing image icon on the upper left corner)... But the following Code does render everything in the correct way:               var layern2 = new OpenLayers.Layer.WMS("Tech",                       "ServerPathTo_mapserv.exe?map=ServerPathTo_My.map&REQUEST=getmap&VERSION=1.1.1&LAYERS=Layer1&FORMAT=PNG", { layers: 'Layer1', TRANSPARENT: 'TRUE' });      Does anyone know what I'm doing wrong? Thanks in advance.


title: Scratch Environment

text: I have a GP script and I'm wondering if I should use Scratch GDB at the processes or just regular GDB, The problem is that this is a read only environment and my temp layers need additional processes and editing   1. What are the advantages using Scratch GDB for saving temporal data of geoprocessing script ?   2. Whats the difference between those options in terms of performance ?


title: Uninstall opengeo suite in Ubuntu

text: How can I properly uninstall the complete opengeo suite in my Ubuntu 12.04? I had installed opengeo suite but it seems there was some error in the PostGIS extension to the Postgres of the OpenGeo suite. So I need to now reinstall it. What is the best and safe way to uninstall it properly and reinstall a fresh copy?


title: Solving the Point in Polygon Problem Using Google Maps and Fusion Tables

text: The Point-In-Polygon problem is fairly well documented on this forum. My question is if there is a way to use Google's free tools, i.e. Google Maps, Fusion Tables, and possibly Google Refine, to solve this problem. I am ultimately attempting to create a heat map of neighborhood boundaries. In basic language here are the steps that I wish to perform: 1. Geocode a list of home addresses using Fusion Table's built-in feature. This yields a set of lat/long coordinates for each address. 2. Import a Shapefile of polygons into Google Fusion to store the boundaries (polygon coordinates). *I am able to get this far.. The next steps are where I am looking for some guidance. 3. Apply a point-in-polygon algorithm to identify which addresses fall within each boundary. What is the simplest way to do this considering that I am ultimately attempting to color-code each polygon depending on the number of points that lie within its boundary lines. Thanks,


title: v.dissolve renames the attribute table column and does not produce a dbf

text: When I run v.dissolve on my vector...               v.dissolve --o input=pass2@mapset1 layer=1 column=fttcp output=pass3      The dissolving appears to work but the column name in the attribute table has changed from 'fttcp' to 'cat'. Also, although the input has a dbf file, the output does not (so I assume it is storing the attribute table in a different format or something?). Is there a way to tell v.dissolve to retain the original column name or is there a way to rename it? Also, is there a reason v.dissolve doesnt write a dbf. Is there a way to get it to write the attribute table to a dbf or is there a way to convert the vector so it uses a dbf? I am using GRASS 6.4.2 on windows 7


title: What ArcGIS 10 geoprocessing tools support exporting tabular data to CSV?

text: I have found that most if not all geoprocessing tools that allow tabular output are limited to geodatabase tables and DBF files for supported output formats. Are there any that also support CSV output? Numerous tools supports CSV input, such as Table to Table, Copy Rows, etc. The only exception I have seen is a script tool in the Spatial Statistics toolbox called Export Feature Attribute to ASCII, but this only supports feature classes, not tables, and it chokes on non-ASCII characters:               <type 'exceptions.UnicodeEncodeError'>: 'ascii' codec can't encode character u'\xa0' in position 213: ordinal not in range(128)`      **Edit:** Thanks for the alternatives in the answers so far. I actually do use some of those those already myself, but just to be clear my question is specifically about the built-in geoprocessing tools in ArcGIS 10. Alternatives are welcome but I am more specifically wondering why the capability to export to CSV exists interactively in ArcMap but apparently not through the GP framework (without Python scripting). Incidentally, I trawled around the ArcGIS Ideas site a bit and did come upon a few existing ideas that may be of interest to some:   * Add Text Option to Table to Table   * Model Builder only Tool - export to csv   * Export a table directly to an MS Excel file (not specifically about CSV but has a lot of relevant discussion


title: How to calculate for every basin/sub-basin, drainage area at every elevation contour along stream?

text: With DEM as source, I created all necessary files like flow direction, accumulation,slope,streams,watersheds,pour point etc.... In my watershed(basin), I have many sub-watersheds(sub-basin) with streams and outlets at various elevations I need o do analysis on the following two basis.   1) I need to calculate for every basin/sub-basin, the drainage area at every elevation contour along the stream. Please ref figure for better understanding my view.   2) I the figure there are six elevation contours (2 black top-bottom, 4 red in between) in the basin/sub-basin. In the figure, outlet is at the bottom elevation. Now I want to change outlet to all above elevations along the stream. And for each elevation point outlet along the stream, I need to calculate drainage area and volume. I am trying various options with Arcgis,globalmapper,autocad,surfer,microstation and others.. Please give me some idea in doing this watershed analysis with less effort with any of the above said software/tools/extensions/toolbox with better accuracy. ![Analysis](http://i.stack.imgur.com/yCZuw.jpg) Thanks much for your time and support.


title: Importing Rasters to ArcSDE 9.3.1 using sderaster

text: OK, I am trying to import a raster to 11g using the sderaster command. It is a stand alone jpg with a world file associated, and pyramids. Am I right in assuming I cannot import a JPG, I can only convert to a jpg on the import process i.e. store it as jpg, but not import a jpg I am currently lost as importing the jpg file I get an "Error: Invalid input file specified" error. I know I can do this via sdeexport - sdeimport, but I believe the correct process should be the sderaster -o import process I am sure we were able to do this with ArcSDE 8.x Any ideas people? Thanks in advance.


title: ArcGIS JSAPI Toggle Identify (AMD style)

text: I am trying to convert the Toggle Identify example to AMD style. However, I am stuck with an empty map. Debugging in Firebug doesn't help too much. Could anyone help me to find where I am wrong? The code is in here: http://jsfiddle.net/neWnE/ I am trying to modify Blord Castillo's code from this thread: Identify with Layer Toggling Any support appreciated!!!


title: Help needed editing existing coordinate information with python

text: I'm hoping someone can help me with this. I would like to find a way to edit shapefile (polyline) coordinate data using python. Basically, I want to edit all the z-coordinates in my shapefile so that they are 1 metre lower than they previously were. I noticed that coordinate data can be manipulated en masse in the Esri Python tutorial video (https://www.youtube.com/watch?v=ge_U_w0PcFM, at 18:17 for anyone thats interested). They use a function within the field calculator to move all the x coordinates +5000m. I have reproduced the script below:               def shiftXCoord(Shape)       point = shape.getpart(0)       point.X  point.X + 5000       return point     shiftXCoord(!Shape!)      I am experiencing two main problems using this script as a solution to my problem. (1) The commentators do not go into any great detail about what particular parts of the script mean/are used for. I know the above is a function but I do not know what lines 2 and 3 represent. (2) Also, without the field calculator I am unsure as to how to link this to the shapefile I am working on. I would appreciate any help you guys are willing to provide on this. Either using the method described above or another method that I am not aware of. Also, if anyone knows of any links/literature that focuses on editing existing coordinate information I would appreciate a heads up. I am using ArcGIS 10.0 by the way, Thanks


title: Why does Exporting Document to JPEG using ArcObjects in Layout View result in bad picture quality?

text: I want to export the Active View with ArcObjects to JPEG. All Layers in the Document are raster layers (GeoTiff). I use this **code** which works quite well when in Data View, but when I perform the same export in Layout View, the picture quality is really bad. Has anyone an idea what the cause of this problem could be? If I do File -> Export Map picture quality is good in both Data and Layout View. So Esri seems to be doing it differently in ArcMap. EDIT: kenbuja had the answer. In case somebody needs the C# variant of the method he mentioned, have a look at the esri documentation


title: How current is GDAL in Java

text: I have used the ogr2ogr to import various formats to Postgres and have used GDAL utilities. Is there a complete Java port for these tools? Or are they simply wrappers calling JNI to C++ code? I have seen some Java stuff by googling but not sure how much are they current and kept upto date. Thanks, Ramesh


title: Fiona: won't install in Debian7.x-64bits/Python2.7

text: I need to install Fiona in my Debian7.x-64bits/Python2.7 machine. I already installed GDAL/OGR, pip, virtualenv, and six. I am using Python 2.7 so argparse is also installed. The installation protocol provided at the Fiona website is not completely clear to me:               $ mkdir fiona_env     $ virtualenv fiona_env     $ source fiona_env/bin/activate     (fiona_env)$ pip install Fiona      (Source: https://pypi.python.org/pypi/Fiona ) I followed the steps from both a non-root- and a root-terminal. In the first case some problems related with permissions arose. In the second the errors were related with an inability to find OGR directories. I am running the installing protocol from my home directory, e.g. me@mymachine:~($ or #). I am not sure if that is the problem. I already checked around for documentation in the Internet, but it seems that the installation protocol from Toblerity (referenced above) is the most explicit. I would appreciate the community help with the installation of Fiona in my rather standard Debian machine.


title: openlayersplugin qgis 2.0.1

text: I have installed qgis 2.0.1, great thing... but.. open layerplugin doesn´t work proper. When I try to load two different open layers I get the answer "could not draw google satellite because:" Also there are problems with loading one map. The correct display is only possible in the open layer projection EPSG:3852, when I try to change the project projection the display gets wrong. This is depending on the zoom factor of the display. At the same time the open layers collaps. I set the project properties on "on the fly"..so my latin is ending here. Many thanks for help! Andreas


title: Populating newly created fields

text: In ArcGIS, i created a new field in my shapefile table and would like to type names there, but cannot. what do I need to do in order to edit the field?


title: Rasterlite db with multi table

text: When working with Spatialite DBs, I can store multiple tables in single DB. But now I want to store multiple raster tables on a single Rasterlite DB. I tried QGIS for conversion but it converts only one table not more... (I am using Windows)


title: How to build query expression to display text from one-to-many relationship in ArcMap?

text: I have two tables in an enterprise geodatabase stored in PostgreSQL and accessed via ArcSDE: parcel and person 1 or + parcel can have 1 or + person (relationship is pll_person) ![enter image description here](http://i.stack.imgur.com/dGdh1.png) I can display all parcels with annotation "id" What I want is to display in each parcel the Name of the person(s) in the related table - see picture please : ![enter image description here](http://i.stack.imgur.com/yRYjj.png) **I want to display parcels proprietor name as annotation please any suggestion ?**


title: Change SDE Version in ArcGIS JavaScipt API depending on connected user

text: I am programming for an web editing application with ArcGIS API for JavaScript. I've published a version of ArcSDE with Feature Acces and I need to know how can I change in the web application from this version to another through the user that is connected to the application using the security manager of ArcGIS Server. I've seen that exist a parameter gdbVersion and SetGDBVersion of featurelayer, but I don't know how to use. I need to see some example of this because the information in the documentation is not enough. Can anyone help me? Any sample, please? Thanks


title: How to use label classes, expressions, and annotations?

text: I've always found that when it comes to python, vector analysis, QGIS gives ArcGIS a run for its money, hell some of the qgis plugins exceed ArcGIS's performance by a few hours; however, I always have found that compared to ArcGIS, QGIS's labeling system lacks the expressions, sql, and detailed representation. I've tried mapnik, and it's nice, but still it seems to pretty up the problem rather than fix it.   * How can I improve my own maps done on QGIS, especially with label classes, expressions, and annotations?   * What needs to be done to enhance QGIS's labeling system in terms of someone writing a plugin or something similar?


title: How to get a polyline from a polygon in ArcGis Server 9.3?

text: I'm want to get the lines of a polygon and then get the length of each line, I'm looking on how to do this, Can anyone show me some code? Does anyone have any pointers on how to do this? I'm using ArcGis Server 9.3.1 and Visual Studio .NET. Thanks.


title: Using ESRI.ArcGIS.DataManagementTools.RepairGeometry

text: I am trying to update some code from VB6 to VB.Net 2008 and have been generally successful. However, I came across an issue where, as it processes each polygon in the shapefile, it implements the following code. `Dim lngExtRing as Integer lngExtRing = pPolygon.ExteriorRingCount - 1   ReDim pExtRings(pPolygon.ExteriorRingCount - 1) pPolygon.QueryExteriorRingsEx(pPolygon.ExteriorRingCount, pExtRings(0)) pExtRings(lngExtRing)` Once that code runs and you hover over the pExtRings portion of pExtRings(lngExtRing) for some polygons (not all), the first popup lists a ComObject but if you expand the plus, it indicates "Nothing". I have run this on a shapefile that has 5 polygons and the two polygons that it fails on have negative areas. There are two ways that might resolve this, but I don't know how to implement either. One was to reverse the orientation if the area was negative, but I tried that on a different shapefile and it still failed on a polygon that did not have a negative area. The other, perhaps better, approach would be to implement the RepairGeometry class of ESRI.ArcGIS.DataManagementTools, but I have no clue on how to make it happen. If anyone has any insight or examples of how to use RepairGeometry, I would certainly appreciate being enlightened. Thanks Dale


title: Traceflowsolver

text: I am working with VBA in ArcGIS.I want to use "Traceflowsolver.findpath" but I recieved an error on "SegmentCost".I need to see an example about this subject .anyone can help me? Also my code is: " traceFlowSolver.FindPath esriFMDownstream, esriSPObjFnMinSum, junctionEIDs, edjeEIDs, count, Segment"


title: How to store historical versions of road (route) portions?

text: Part of my job is to analyze certain portions of road for our clients. A client will ask, I want to look at Road A from parcel XYZ to parcel RST and Road B from parcel RST to parcel OPQ. The way I do this currently is:   1. I select the road line layer around the parcels   2. Save the selection as a new shp layer   3. edit, trim, the layer as needed    4. produce a new shp that matches the client's request.    5. We analyze the hell out of the road The reason I do selection/save as is that, in future when road geometry changes we will still have the original requested road. The problem I find with this is that by the end of the year, i have 100+ mini layers of road segments. What are the issues in my process, and how can I improve on them?


title: Display photo in GeoExt popup

text: I have point feature with name PointPhoto. This feature have attribute:   * id    * area    * PhotoPath Now I am using GeoExt popup to display this attributes, but I cant find a way to display image from image path, this is the code:                <html>          <head>              <link rel="stylesheet" type="text/css" href="ext/resources/css/ext-all.css">              <script type="text/javascript" src="ext/adapter/ext/ext-base.js"></script>                  <script type="text/javascript" src="ext/ext-all-debug.js"></script>              <script src="../OpenLayers-2.9.1/OpenLayers.js"></script>              <script type="text/javascript" src="geoext/script/GeoExt.js"></script>                <link rel="stylesheet" type="text/css" href="http://localhost:8181/geoserver/openlayers/theme/default/style.css"/>             <script type="text/javascript" src="../proj4js/lib/proj4js-compressed.js"></script>             <link rel="stylesheet" type="text/css" href="geoext/resources/css/popup.css">              <script type="text/javascript">      var app, items = [], controls = [];      var store;     Ext.onReady(function() {                  app = new Ext.Viewport({                      layout: "border",                       items: items                  });              });              items.push({                  xtype: "gx_mappanel",                  ref: "mapPanel",                  region: "center",                  map: {                      numZoomLevels: 19,                      controls: controls                  },                  extent: OpenLayers.Bounds.fromArray([                          2189574.7042680467, 5418655.108865695,                         2199298.319819289, 5426457.253025988                         ]),                  layers: [new OpenLayers.Layer.WMS(                      "cite:PointPhoto",                      "http://192.168.1.2:8181/geoserver/cite/wms",                      {layers: "cite:PointPhoto"},                      {isBaseLayer: false}                  )] ,                 });                  controls.push(             new OpenLayers.Control.PanZoomBar({position: new OpenLayers.Pixel(2, 15)}),                     new OpenLayers.Control.Navigation()                );               controls.push(new OpenLayers.Control.WMSGetFeatureInfo({      autoActivate: true,      infoFormat: "application/vnd.ogc.gml",      maxFeatures: 3,      eventListeners: {              "getfeatureinfo": function(e) {                  var items = [];                  Ext.each(e.features, function(feature) {                      items.push({                          xtype: "propertygrid",                          title: feature.fid,                          source: feature.attributes                               });                  });                       new GeoExt.Popup({                      title: "Feature Info",                      width: 200,                      height: 200,                      layout: "accordion",                      map: app.mapPanel,                      location: e.xy,                      items: items[0]                  }).show();              }          }      }));                     </script>          </head>          <body>               </body>      </html>       ![enter image description here](http://i.stack.imgur.com/5ImLF.jpg) My question is how I can display image inside of this popup?


title: Arcpy PageLayout.dll Errors in Python

text: import os         import arcpy         arcpy.env.overwriteOutput=True         folderPath= raw_input ("Enter directory:  ")         filePath= raw_input ("Enter data repository:  ")         f = open((os.path.join (filePath, "MXDPrintout.txt")), "w")         print "Name of the file: ", f.name         os.chdir(folderPath)         for root, dirs, files in os.walk(folderPath):             for file in files:                 if file.endswith(".mxd"):                     results=os.path.join(root, file)                     myString= str(results)                     f.write(myString + os.linesep)             for filename in os.listdir (folderPath):                     fullPath = os.path.join(folderPath, filename)                     if os.path.isfile(fullPath):                         basename, extension = os.path.splitext(fullPath)                         if filename.lower().endswith(".mxd"):                             mxdName=filename                             mxd = arcpy.mapping.MapDocument(fullPath)                             mxd.relativePaths = True                             mxd.save()             del mxd         f.close ()      Scripting to parse a mature .gdb for all mxd's and set relative pathways to True before conducting data management. This script worked fine for a few runs but then failed giving this error message: "Access violation at address 65055D43 in module 'PageLayout.dll'. Read of address 00000000." Ive searched for a debug but cannot seem to track down the fix. The .dll is Esri code. Im running Arc 10.1 SP1. Any thoughts on how to debug or fix the .dll error?


title: Custom BatchProject Script

text: Here's my problem. My GPS collects data in WGS1984 and I have decided to start working in UTM instead of state plane. Needless to say, I would edit shapefiles, not paying attention to what coordinate system I was using and relied on the transformations. I have realized the error of my ways and want to correct my workflow. I would like to take my original files from my gps and do a batch project to UTM so I can correct the data. The directory structure of my gps software is such that each project has it's own folder with subdirectories of corrected/uncorrected data and the shapefiles. I would like to be able to load up the tool with the various files from different locations and save the UTM file to the same directory as the WGS1984 file, just with UTM at the end of the file name. What I have done is made a copy of the BatchProject script from ArcToolbox to modify, but I have little to no experience with Python. If this is not the way to approach this problem, please let me know. If I should be using just the regular project tool and put that inside of a loop of some sort, that's fine too. Any help would be greatly appreciated!


title: Using ECW with mapserver

text: I know that mapserver drop the support for ecw due to licensing issues. But if i understand well, as mapserver rely on GDAL, if GDAL support ECW, mapserver can serve ECW data. So I've compiled libecwj 3.3 , then compile GDAL. Now GDAL can read ECW. But in mapserver, ecw data don't work. The error is : > MsDrawMap(): Image handling error. Failed to draw layer named 'ortho'.; > MsDrawRasterLayerLow(): Unable to access file. Corrupt, empty or missing > file 'my_file.ecw' for layer 'ortho'. I tried to recompile mapserver. But with no effect. Does someone know if it's possible and where is the mistake ? My config is : - debian 7 - libecwj 3.3 - gdal 1.10.1 - mapserver 6.4.0 Thanks for your help !


title: google earth profile used in mashup or flex

text: I really like the way google earth profile works, 1. mouse over line, right click, select profile 2. slide the mouse over the created 3d profile an dsee where that is on the map. Has anyone an idea of whether this capability could be mashed outside of GE (I probably already know the answer)? Or if there is anythig esle as slick that could be integrated with flex or arcgis server? Thanks


title: Organizing Arcmap folders for portability

text: I am doing some basic constraint mapping in ArcMap 10.1 and need some pointers on organizing my project for portability. My project TOC in ArcMap contains shapefiles in the following structure: Project Details -Project Lands Features -POI's, roads, railways etc Constraints -Setbacks from features Typically the folder structure for our Arcmap and data looks like this: Arcmap file path: Client -> Project -> ArcMap -> Date -> mxd file Shape file paths: Client -> Project -> Data -> Data Type i.e. Project detail, feature, constraint -> Date -> shp file Often we will package up the Arcmap and shp files to send to the client. What is the easiest way to do this without breaking links in Arcmap? With the folder structure the way that it is, I am unable to only pull the Arcmap -> Date -> mxd file and Data -> Date -> shp files. I have played with exporting/importing to/from file geodatabases in ArcCatalog but don't think this is the solution other than it being just one file rather than a bunch of different shape files to send to the client. The other solution that I think should work would be making the paths relative but it seems whenever I change a folder name it breaks the data source and I end up having to reconnect everything which is error prone. Any tips would be much appreciated.


title: Adding feature classes from within multiple geodatabases within multiple folders

text: I have thousands of feature classes that I'm trying to add to an ArcMap 10.1 map document. I have a folder of U.S. DEM files, and within that folder is a folder for each U.S. state. Inside each state folder are multiple geodatabases (20-100, depending on the size of the state). Each feature class is then located within a feature dataset that is within a geodatabase. (C:\DEM Files\State_Folder\State_location.gdb\featuredataset\featureclass) I am simply trying to add all of the feature classes at once. I have successfully written code that will add feature classes one at a time. I have also attempted to list the feature classes but this was impossible because each feature class has the same name ("Elev_Contour"), and arcpy was unable to list them. Any help with this problem will be greatly appreciated! Thanks!


title: PostGIS views vs WFS Filtering operations on PostGIS table

text: In what use cases should I prefer writing PostGIS views for filtering datasets versus filtering with WFS filter encoding operators and vice versa? In both cases I would be serving the data via GeoServer through the PostGIS table or view. The main pro I can think of for the WFS filter approach is that the query can be encoded easily through the front end according to the filter specifications without needing to add extra logic to the backend code to handle the PostGIS queries. Also, how does the WFS filter fare performance wise against PostGIS?


title: Django PolygonField: how to display on map?

text: Extremely beginner question. I have a GeoDjango model with a MultiPolygonField. How do I retrieve the value of this field and display it as a polygon on a map?               class county(models.Model):         cd = models.CharField(max_length=10)         geom = models.MultiPolygonField(srid=27700)         objects = models.GeoManager()          def lsoa(request, code):         county = get_object_or_404(county.objects, cd=code)         return render_to_response('county.html', { 'county': county }, context_instance = RequestContext(request))      I know GeoJson is involved somewhere along the line, but I'm struggling with to get it out of Django, and with how to display it.


title: Can you create new annotation feature classes with Python in ArcGIS for Desktop?

text: I am working on a project that needs to create file geodatabases from scratch and I cannot find any information about how to make new feature classes as annotation with Python. Does anyone have any thoughts or suggestions on this? Is there another way to do this programmability outside of Python? This question was written for ArcGIS Desktop 9.3 but Answers using ArcGIS 10 for Desktop are also welcome.


title: How do I display geojson polygons on a map?

text: Here's a very basic question for you all. I've got 17,000 polygons in a GeoDjango database. How do I create a slippy web map that can display these polygons? In other words, do I:   * load a GMap, get the boundaries, each time the map moves, make an Ajax call to the Django db, load the polygons...   * make my own map tiles in advance (if so, how?)   * something else? I'd like users to be able to view the polygons on the map. I'd also like users to be able to search for a point, using GMaps geocoder or something similar, and load up the polygons within the borders of the map. What client-side technology should I be using? Thanks.


title: Openlayers proxy.php Curl settings for geoserver

text: We have to use a proxy to send a request to geoserver. I created a php file like this.               OpenLayers.ProxyHost = "proxy.php?url=";      is Openlayers settings.               localhost/myapp/proxy.php      is my proxy adres. Proxy content is like this.               <?php     $headers = array(         "Content-type: application/xml",         "Content-type: text/xml",     );           $ch = curl_init($_GET['url']);      curl_setopt($ch, CURLOPT_HEADER, 1);      curl_setopt($ch, CURLOPT_HTTPGET, 1);      curl_setopt($ch, CURLOPT_FOLLOWLOCATION, 1);      curl_setopt($ch, CURLOPT_RETURNTRANSFER, 1);      curl_setopt($ch, CURLOPT_HTTPHEADER, $headers);      $output = curl_exec($ch);      curl_close($ch);      echo $output;     ?>      I am using this proxy like this. `localhost/myapp/proxy.php?url=http://somedomain.com` and works. Web page is opening good. But when I try geoserver WMS get query, it gives error, NullPointerException. For example: `http://myip/geoserver/app/wms` is returning admin page, but `localhost/myapp/proxy.php?url=http://myip/geoserver/app/wms` is returning: **HTTP/1.1 200 OK Server: Apache-Coyote/1.1 Content-Type: application/xml Transfer-Encoding: chunked Date: Wed, 03 Sep 2014 14:18:11 GMT java.lang.NullPointerException null** `http://myip/geoserver/app/wms?request=GetCapabilities` is returning xml page, but               localhost/myapp/proxy.php?url=http://myip/geoserver/app/wms?request=GetCapabilities is returning plain text result.


title: Issues Merging in Model Builder after Running Script

text: I have multiple inputs running through various processes throughout my model. At the end of these processes they are all merged together and output into one final feature class. This was working fine until I decided to put a script tool at an earlier point in the model on one of my functions, essentially stopping the process if a certain value was false. The script tool works well, but now my final merged output is not performing correctly. I took the model apart in pieces and ran through to isolate the problem, and it would appear that when the script tool stops the process, even though the merge function has multiple correct inputs, it will not merge anything together. Any ideas on how to resolve this issue or a workaround? I'm using ArcGIS 10. Thanks.


title: Draw a skyplot with calculated ECEF values (in meters)

text: I have ECEF values that I calculated from ephemeris subframes 2/3. Using this information, I'd like to place these values on a skyplot based on my current location. I'm using a C# picture box for my skyplot image, and a custom control to show the dot and SV PRN. I am struggling in these areas:   1. How to merge the ECEF information with my current location to provide a useful x,y for my skyplot.   2. How to determine the bounds of my drawing. Currently i estimate that from edge to center of my drawing is 90 pixels. So i attempted to do some scaling 90 / 20,000,000 (max distance for satellite, guessed), then multiplied that by my ECEF.x, and ECEF.y values. I'd appreciate any help. Thanks, Jason


title: How to increment the data frame time by one day in ArcMap 10?

text: I often use ArcMap 10 to display dynamic text using the dataFrame time.               <dyn type="dataFrame" name="Layout" property="time" emptyStr="[off]"/>      The above might display something like               2011/08/16      However what I really want to be able to do is add 1 day to the dataFrame time and display it. This way I can display a date range               2011/08/16 - 2011/08/17      Thanks


title: How to use Route Task in Network analysis tools

text: I am using esri network analysis tool and I just want to draw the path between two clicked point on the map (by using route task) I use the example code block in the following link http://help.arcgis.com/en/webapi/silverlight/samples/start.htm#Routing which is included in the ESRI ARCGIS silverlight API sample page and I only change the route task and Dynamic layer URLs... In the end I always receive the error message as "Routing error: 'stops' parameter is invalid" Can anyone help me? Thanks a lot


title: Best ways to embed dynamic maps into PDF files

text: Working in consulting, I found that sending clients PDF file with my embedded maps is one of the best ways to achieve a wow-effect from them. The latest pdf file formats all support the layers feature (which arcmap makes use of); however they also are able to support javascript. Does anyone have ideas to use the power of javascript (or anything else) to make my maps within my pdf files more dynamic??? Which tools do you recommend? Any extensions to esri products? Thank you very much!


title: Geomedia resources

text: Are there any Intergraph Geomedia users visiting this site? As of today, there is not one question about Geomedia.   Do you know about any site where there is any discussion on Geomedia products (besides Intergraph Synergy site which is not very active lately) ? **Note:** would somebody, please, add tags Geomedia and/or Intergraph to this question because I don't have enough reputation to create new tags.


title: Add legend from raster catalog? ArcGis 10.1

text: Does anyone have a solution on how to add the legend from a raster catalog? It does not show up, either in version 10.1 or 10.2.1. It only shows the name "Legend".


title: 9.2 to 9.3.1 Desktop - How to troubleshoot slow data access with shps?

text: We have just upgraded from ArcGIS Desktop 9.2 to 9.3.1 **For a ~250MB .shp in 9.2 it took 1 minute to render, at 9.3.1 it now takes 50 minutes.**   * The format has to stay as shapefile.   * The data is stored on both NetApp and Dell Servers - problem consistent on both.   * Only an issue with GIS data, in particular shapefiles and only in 9.3.1.   * Have already have tested during the day and night and the times returned are consistent showing that network traffic is having no or very little effect on the poor performance.   * Also tested on servers that do not have virus protection to eliminate that from the equation.   * Already tried testing on several servers with differing switches (100Mbit and 1Gbit) but the results have been the same.   * Have tested on a newly built pc with a vanilla install of ArcGIS so it can’t be incorrectly uninstalled components or third party products causing the issue. I am putting together some further test-cases to carry out, to try and pinpoint the problem further before taking it to Esri:   * Run the tests at 9.2 with a large ESRI sample dataset. Run tests with it as a shp and as a file GDB.    * Re-test with 9.3.1 (+SPs) and again with v10   * Run tests in both ArcMap, ArcGIS Explorer and ArcReader (to see if it is specific to an application)   * Compare these results with how long it takes to open when stored locally? From what I know, there have not been any changes made to the underlying code and the way that ArcGIS connects to a network between 9.3.1 and 9.2. I also don't think that Esri changed anything in the way that ArcGIS reads shapefiles between these versions? Other than the above, the only other suggestion I can think of is to run WireShark to try and see how long it takes to send/receive packets from the client machine & the server. This will tell me two things: 1) How many packets are going both directions 2) how long each packet takes to be transmitted. My thinking is, if it is related to 9.3.1, then you would most likely see an increase in the number of packets. It could be a security configuration that at 9.3.1 we have encountered, but at 9.2 we did not. It would be best to capture a wireshark trace for both 9.2 and 9.3.1 accessing the same shape file on the same server. **Has anyone experienced anything similar to this?** **Does anyone have any other ideas on how to troubleshoot this?**


title: Geonames query for Alabama Cities does not bring two biggest cities

text: I have this query: http://api.geonames.org/search?featureCode=PPL&username=genadinik&country=US&style=full&adminCode1=AL&maxRows=100 Which brings back the top most populated cities in order by population. The biggest cities in Alabama are Birmingham and Montgomery. But they do not appear on this list, no matter how far down the list I extend the query. Is there a reason that this happens? Or some rule/design feature that I am not aware of that is causing this? Thanks!


title: Using unicode characters (arabic) in Select By Attribute of ArcGIS for Desktop?

text: I'm new to Python scripting with ArcGIS for Desktop and I have unicode characters (arabic) in my attribute table. I want to use these characters in my script but when I try to do that the script does not work. How can I use unicode characters in my script?               import arcpy     arcpy.env.workspace = "c:/test"     mxd=arcpy.mapping.MapDocument("CURRENT")     listlyr = arcpy.mapping.ListLayers(mxd,"","")     # hatta agar chand dataframe ham dashte bashim baz ham select ra anjam mdahad     for a in listlyr:         delim = arcpy.AddFieldDelimiters(a,"City")              arcpy.SelectLayerByAttribute_management (a, "NEW_SELECTION", delim+"='u'\xd9\x86\xd9\x88\xd8\xb4\xd9\x87\xd8\xb1''")      **u'\xd9\x86\xd9\x88\xd8\xb4\xd9\x87\xd8\xb1'** means 'نوشهر'


title: Problem with CadTools plugin in QGIS 2.2

text: I've just downloaded the CadTools plugin from QGIS Plugins page (can't do this through QGIS itself as I'm behind a proxy) and after installation I keep getting an error saying that it " _Couldn't load plugin cadtools due an error when calling its classFactory() method_ ". Looking further into this there appears to be a problem with a line in the cadtools.py file relating to _self.localPtah.pluginPath_. pluginPath seems to be set to look in the _i18n_ subfolder but I don't appear to have this subfolder in my installation. Has anybody else had a similar issue with CadTools or have I just missed a step with the installation??


title: What is the best way to convert WKT String to GML String

text: What is the best way to programatically convert WKT String to GML String? Every reasonable library (at best available in some Maven repository) can be used. Till now I came up with this solution using JTS. Its written in Java:               public String wktToGml(String wktGeometry) {        WKTReader wktR = new WKTReader();        Geometry geom = wktR.read(wktGeometry);             // write JTS to string        GMLWriter gmlW = new GMLWriter(true);        String gml = gmlW.write(geom);        return gml;     }      But I'm not sure if it's correct and most effective way. There are probably some problems with LF (linefeed) symbols. Probably there should be also solution using geotools library. I'm new to this field, I'm trying to gain some overview.


title: ArcGIS export layout with inserted EMF graphics

text: using ArcGIS 10 SP3 for some simples mapping tasks. But::: I have a layout with some inserted EMF graphics. These graphics are in fact formatted Excel- tables. The first rows in the tables are shaded with a light gray background. The EMFs are drawn correctly in the ArcMap layout, but dependent on the position on the page the shades disappear in the exported maps: if the graphic is positioned at the right page border the shades disappear, at the bottom or left side, the output will be correct. I have copied and pasted the same graphics a dozen times and placed it side by side. It seems, that all graphics which lay mainly outside the DIN A4 portrait format are not correctly exported, independent of the page settings (in my case landscape). In the meanwhile I have deleted ALL OTHER layout elements and the effect still is the same. Has anybody any idea? Is this one more of ArcGIS strange features?


title: Issue while running grass from subprocess

text: I have a map based web application using geodjango. I needed to display the crime hotspots within a selected polygon area. I found that v.kernel is very good for generating heatmaps/hotspots. So in my application I ran this command through python subprocess and created the raster map showing point density. I displayed this raster map over the polygon area in grass. But when I try to run the command through python using d.out.file command, error is raised.               command5 = "d.mon start=x1"      command6 = "d.mon select=x1"     command7 = "d.mon unlock=x1"     command8 = "r.colors byr"     command9 = "d.rast crimes_"+str(pname)     command10 = "d.vect "+str(pname)+" type=line,boundary"     command11 = "d.out.file output="+str(pname)+" format=png size=800,600"     command12 = "d.mon stop=x1"          saveCommands = [command4,command5,command6,command7]          for sc in saveCommands:         s = subprocess.Popen(sc,shell=True,env=environ)          s.wait()      Error is :               ERROR: Graph_Set: can't open Display     ERROR: Unable to get GIS_LOCK enviroment variable value     Problem selecting x1. Will try once more     ERROR: Unable to get GIS_LOCK enviroment variable value     d.mon select=x1     Exception AttributeError: AttributeError("'_DummyThread' object has no attribute '_Thread__block'",) in <module 'threading' from '/usr/lib/python2.7/threading.pyc'> ignored     ERROR: Unable to get GIS_LOCK enviroment variable value      All I need is to export the grass vector and raster layer to an image.


title: Cut a layer into many parts using a shapfile

text: I have a layer which is a map of a town, and a shapefile which shows the boundaries within that town. I wish to cut the layer up along each of the internal boundaries so that I end up with many layers, each containing the map of one small area as defined by shapefile. I know that I can use the clip tool to cut the map to the shape of the whole shape file, but I don't know how you use the boundaries contained within the shapefile to cut the map. I am using ArcGIS Desktop 10 SP5 with the ArcView License


title: How to alter this PostGIS ST_distance_sphere query to give the answer for all points in the table, not just one?

text: I have the following query which returns the expected result for a single point. I would like to alter this query to return the answer for every point in the table.               SELECT count(case when(                        ST_Distance_Sphere(                            (select the_geom                              from brc_property where property_r=20), the_geom) < 500000)                         then 1 else null end)      FROM brc_property


title: Geoserver cascading and fetching projection

text: I tried to use cascade of WMS from this portal in order to cache wms layers http://geoportal.dgu.hr/wms?request=GetCapabilities Native for this portal is EPSG:3765, but EPSG:31275 is also supported But while direct access to WMS returns correct result: http://geoportal.dgu.hr/wms?layers=DOF&FORMAT=image%2Fpng&SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&STYLES=&SRS=EPSG%3A31275&BBOX=5578043.9690343,5053103.9883878,5585211.9651635,5060271.9845171&WIDTH=256&HEIGHT=256 cascaded WMS with exactly the same parameters returns all-white image. And in documentation of Geoserver it says that while cascading if the requested projection exists on remote server it will serve it from remote server, otherwise it will reproject it. What could be going wrong, why isn't geoserver just returning correct image from remote service? Geoserver version used is 2.2.


title: Open source alternatives to Google's StreetView engine?

text: I am looking to create an open source javascript app using either Leaflet or OpenLayers. Is there any way to implement panoramic views (similar to Google's StreetView) into either of these, or any other open API? My panoramic views would be custom made from personal photos so I will not need to grab panoramic images from an existing service. I understand that it is possible to incorporate Google StreetView into non- Google apps, but is there anything out there that would give me similar functionality without using any Google API?


title: Is there any known reason for QGIS recurrently crashing in windows7?

text: I use (and been a big fan of) QGIS daily for 7 months now, and only recently - past couple weeks - been experiencing problems with QGIS recurrently crashing. This is even more serius cause it happens when i try to save the editions I'm doing. I'm using a windows7 machine and i'm performing the same task for 2 months now which consists on visually interpreting a satellite image to digitalize features to a shapefile. I haven't change anything and I don't understand why this occurs but it's very annoying to have QGIS crashing every half an hour... Any ideas? similar experiences? Thanks, Vasco Nunes PS: I'm using version 1.7.3 but this has been happening since I had 1.7.2


title: Error publishing feature service in GISserver 10.2

text: I installed ArcGIS 10.2 for Server just now. I started to create a map using the ArcMap application of ArcGIS for Desktop. I created My map with just: 1 shape file + 1 gdb file... there were not any error when I push the anlyze button(I mean I have checked the enterprise db for my gdb file to be published correctly). After all, when I want to publish my map as a feature service(I mean a service map with the 'feature access' capability enabled) I encounter an error like this(the feature service definition is created with no errors but just not published and running): ![error 001369 when publishing my feature service ](http://i.stack.imgur.com/gxKyC.jpg) I searched all the web but there where no solution for the error!!! actually there where no useful log in GP results window as you see below: ![GP results window](http://i.stack.imgur.com/4ITD8.jpg)


title: Sorting by geometry

text: I'm trying to implement sorting of features by set of fields. These fields can contain geometry one. For geometries I want to use adjustable XY- and Z-tolerance. My thought was to use ITableSort with custom ITableSortCallBack:               private class SortCallback : ITableSortCallBack     {         public int Compare(object value1, object value2, int FieldIndex, int fieldSortIndex)         {             // ... custom compare operations if value1/2 is IGeometry                  return 2;         }     }      But for SHAPE field value1/2 is always null. Documentation says: > For example, by default, a shape field or a raster field cannot be sorted > on. But I don't have default behavior because I use my custom comparer. However I can't do anything with shape because I always get null geometries. I use cursor to construct TableSortClass but it's guaranteed that rows contain values for SHAPE field. What am I doing wrong? Or maybe there is another easy way to implement sorting of features including the SHAPE field with adjustable tolerances?


title: How to find the center point of MULTILINESTRING?

text: How to find the center point of MULTILINESTRING, which lies on the multiline?


title: Permissive Open source planarize algorithm?

text: Is there a permissive open source implementation out there for a "planarize" algorithm (in the style of arcgis planarize - including user modifiable cluster tolerance)? NB this is not a duplicate of Is there an open source tool to planarize graphs? as alas I need a permissive license, LGPL or BSD, which GRASS doesn't have :(


title: How to set workspace from input Feature Class?

text: How to set workspace from input Feature Class ? is there an attribute of the source GDB for a feature class ? some times the FC is in a dataset and some time it isn't, How can i get the path of source GDB for setting workspace when the input is dynamic ?


title: Extract all pixels of given RGB value and render as MultiPolygons

text: Given a geotiff image with a colortable in it, what is the easiest way to extract each color as its own "layer or shapefile or multipolygon object"? Specifically, I am looking for a programmatic approach to solve this problem. (not a GUI app, QGis/ArcGis approach) Here's what gdalinfo shows:               $ gdalinfo GLOBCOVER_L4_200901_200912_V2.3.color.tif      Driver: GTiff/GeoTIFF     Files: GLOBCOVER_L4_200901_200912_V2.3.color.tif     Size is 129600, 55800     Coordinate System is:     GEOGCS["WGS 84",         DATUM["WGS_1984",             SPHEROID["WGS 84",6378137,298.257223563,                 AUTHORITY["EPSG","7030"]],             AUTHORITY["EPSG","6326"]],         PRIMEM["Greenwich",0],         UNIT["degree",0.0174532925199433],         AUTHORITY["EPSG","4326"]]     Origin = (-180.001388888888897,90.001388888888883)     Pixel Size = (0.002777777777778,-0.002777777777778)     Metadata:       Generate by=gdal_mean       Copyright=Copyright ©UCL Geomatics, BELGIUM 1999-2010       Authors=Sophie Bontemps <sophie.bontemps@uclouvain.be>,Eric Van Bogaert <eric.vanbogaert@uclouvain.be>,Pierre Defourny <pierre.defourny@uclouvain.be>       process begin time=2010-12-23T09:49:37       process finish time=2010-12-23T09:57:38       process files=CL5_GLOBCOVER-L5_CLASSIF_2009_V2.3.20101220.tif       AREA_OR_POINT=Area     Image Structure Metadata:       COMPRESSION=LZW       INTERLEAVE=BAND     Corner Coordinates:     Upper Left  (-180.0013889,  90.0013889) (180d 0' 5.00"W, 90d 0' 5.00"N)     Lower Left  (-180.0013889, -64.9986111) (180d 0' 5.00"W, 64d59'55.00"S)     Upper Right ( 179.9986111,  90.0013889) (179d59'55.00"E, 90d 0' 5.00"N)     Lower Right ( 179.9986111, -64.9986111) (179d59'55.00"E, 64d59'55.00"S)     Center      (  -0.0013889,  12.5013889) (  0d 0' 5.00"W, 12d30' 5.00"N)     Band 1 Block=129600x1 Type=Byte, ColorInterp=Palette       NoData Value=0       Color Table (RGB with 256 entries)         0: 0,0,0,255         1: 0,0,0,255         2: 0,0,0,255         3: 0,0,0,255         4: 0,0,0,255         5: 0,0,0,255         6: 0,0,0,255         7: 0,0,0,255         8: 0,0,0,255         9: 0,0,0,255        10: 255,255,100,255        11: 170,240,240,255        12: 170,240,240,255        13: 170,240,240,255        14: 255,255,100,255        15: 255,255,100,255        16: 255,255,100,255        17: 0,0,0,255        18: 0,0,0,255        19: 0,0,0,255        20: 220,240,100,255        21: 210,240,100,255        22: 210,240,100,255        23: 0,0,0,255        24: 0,0,0,255        25: 0,0,0,255        26: 0,0,0,255        27: 0,0,0,255        28: 0,0,0,255        29: 0,0,0,255        30: 205,205,102,255        31: 240,200,100,255        32: 205,205,102,255        33: 0,0,0,255        34: 0,0,0,255        35: 0,0,0,255        36: 0,0,0,255        37: 0,0,0,255        38: 0,0,0,255        39: 0,0,0,255        40: 0,100,0,255        41: 0,100,0,255        42: 0,100,0,255        43: 0,0,0,255        44: 0,0,0,255        45: 0,0,0,255        46: 0,0,0,255        47: 0,0,0,255        48: 0,0,0,255        49: 0,0,0,255        50: 0,160,0,255        51: 0,0,0,255        52: 0,0,0,255        53: 0,0,0,255        54: 0,0,0,255        55: 0,0,0,255        56: 0,0,0,255        57: 0,0,0,255        58: 0,0,0,255        59: 0,0,0,255        60: 170,200,0,255        61: 0,0,0,255        62: 0,0,0,255        63: 0,0,0,255        64: 0,0,0,255        65: 0,0,0,255        66: 0,0,0,255        67: 0,0,0,255        68: 0,0,0,255        69: 0,0,0,255        70: 0,60,0,255        71: 0,0,0,255        72: 0,0,0,255        73: 0,0,0,255        74: 0,0,0,255        75: 0,0,0,255        76: 0,0,0,255        77: 0,0,0,255        78: 0,0,0,255        79: 0,0,0,255        80: 70,75,0,255        81: 0,0,0,255        82: 0,0,0,255        83: 0,0,0,255        84: 0,0,0,255        85: 0,0,0,255        86: 0,0,0,255        87: 0,0,0,255        88: 0,0,0,255        89: 0,0,0,255        90: 40,100,0,255        91: 40,100,0,255        92: 40,100,0,255        93: 0,0,0,255


title: Automate Select by Location

text: I have 2 shape files (A and B) each with thousands of polygons. I want to select polygons from A which intersect a specific polygon from B, and save the selected polygons from A as a new shapefile. I can do this using the 'select by location' feature and 'create layer from selected features', however I want to repeat this for every single polygon in B. Is it possible to automate such a task? I do not know how to use the 'Python' programming, but may be able to have somebody do this for me if that is the best method? I would really appreciate any suggestions! Thank you, Zoe


title: Is there an open source directions webclient for pgRouting like the directions API from Google?

text: I am looking for an open source webclient with (Google) directions functionality (like draggable lines to create extra waypoints) based on an pgrouting database (with an osm2pgrouting topology). I noticed that there are some examples and demos on the pgrouting site (some are not working anymore), but they are all specific things people have done with pgrouting. They seem to have all created their own specific front end. I am looking for a more generic directions/routing client based on a pgrouting topology.


title: Georeference problem including EPSG 102067 projection

text: I'm Qgis user from Czech Republic. We have a spatial national coordination system (Krovak system) having positive X axes pointing to the south and positive Y axes to the west. Many geographical systems have problem to display such coordination system, so usually is used transformation which gives: X_transformed = - Y_Krovak Y_transormed = - X_Krovak The above described transformed coordination system is described in EPSG database as EPSG: 102067 To my sadness this projection was excluded from Qgis projection database starting Qgis 1.7.4. Many Czech wms servers work with this projection and I need their data. Many GIS projects are accomodated into this transformed coordination system. To bring EPSG 102067 projection into my QGis 1.8 I override its original projection database file which is distributed with projection database file of Qgis 1.7.3 (including EPSG 102067 projection). Name of this file is srs.db My problem is that I can not correctly geoference raster map into EPSG 102067 projection system. Above mentioned information possible could be part of the problem, or absolutely not. My attempts are resulting into georefferenced raster having height and width less than 1 meter. Georeference points and original bmp raster could be downloaded form http://www.1progeo.cz/download.php Login: qgis Password: qgis as a zip archive so to debug my mistake could be quite easy. Would you be so kind to help me to give me guidence how to georeference raster map correctly? Thank you very much. With best regards Jan


title: Select single feature from multiple types with WFS in OpenLayers

text: I have a WMS base and from that create a Vector Layer for selections. Then I create 3 Control.GetFeature() controls, each using WFS to select a different feature type upon clicking the map (the 'featureselect' event). Doing it this way, however, I end up with as many as 3 Features being selected, when I want a maximum of one. Now, ideally I'd like GeoServer to pick the best (closest) feature to the mouse click, but even using Javascript in the browser, how would I detect when all the features have come in via WFS so I can test the geometries? And what would be the best way to test them? (mostly multilines and multilinestrings) Here's the simplified code to give you the idea of what I've got now:               function init() {         var mapOptions = { /* yada yada */ }         var map = new OpenLayers.Map('map', mapOptions);              var wmsOptions = { /* yada yada */ }         var wmsLayers = 'parcelshape_lin,wmain_lin,wservice_lin,buildings_pol';         var wms = new OpenLayers.Layer.WMS(             "Testing",             baseUrl+"/geoserver/wms",             { layers: wmsLayers },             wmsOptions         );              var select = new OpenLayers.Layer.Vector("Selection", {             styleMap: new OpenLayers.Style(OpenLayers.Feature.Vector.style["select"])         });              map.addLayers([wms,select]);              var BuildingsControl = new OpenLayers.Control.GetFeature({             protocol: OpenLayers.Protocol.WFS.fromWMSLayer(wms, {                 url: baseUrl+"/geoserver/wfs",                 featureType: "buildings_pol"             }),             hover: false         });         BuildingsControl.events.register("featureselected", this, function(e) {             select.addFeatures([e.feature]);         });         BuildingsControl.events.register("featureunselected", this, function(e) {             select.removeFeatures([e.feature]);         });                   var WaterMainsControl = new OpenLayers.Control.GetFeature({             protocol: OpenLayers.Protocol.WFS.fromWMSLayer(wms, {                 url: baseUrl+"/geoserver/wfs",                 featureType: "wmain_lin"             }),             hover: false         });         WaterMainsControl.events.register("featureselected", this, function(e) {             select.addFeatures([e.feature]);         });         WaterMainsControl.events.register("featureunselected", this, function(e) {             select.removeFeatures([e.feature]);         });              map.addControl(BuildingsControlCtrl);         map.addControl(WaterMainsControl);              BuildingsControl.activate();         WaterMainsControl.activate();     }


title: Use arcgis with SQL Server 2008 as a backend

text: I'm wondering how best to use arcgis] and sql 2008 together. At the moment we run manual imports and exports between sql 2008 and arcgis 10 shape files. With sql 2008 having a geography data type I was expecting it to be simple to replace to the map files with a relational database, but for some reason it doesn't seem as simple as I expected. Has anyone got any ideas how I can use sql 2008 for data storage and do away with shape files? I've read about arcsde, but really don't understand the point of another layer. No idea what it's actually supposed to do. Thanks, Michael (I've also asked the same question here but think it was probably the wrong place)


title: create buffer around locator-derived point

text: I am trying to create a buffer around a point derived from a locator search, using this example: http://help.arcgis.com/en/webapi/javascript/arcgis/help/jssamples_start.htm ..but I am getting no results. I am not sure how to pass the correct information from the locator result graphic to the buffer functions (sorry if this has been answered before, i couldn't find a solution). Thanks, code is below..Jason                   ////////////// locator func ///////////////////////////////////            function showLocatorResults(candidates)        {         var candidate;         var symbol = new esri.symbol.SimpleMarkerSymbol();         symbol.setStyle(esri.symbol.SimpleMarkerSymbol.STYLE_SQUARE);         symbol.setColor(new dojo.Color([153, 0, 51, 0.75]));         var geom, graphic;              dojo.every(candidates, function (candidate) {                  if (candidate.score > 80)              {                // extent code removed for brevity                 var attributes = { address: candidate.address, score: candidate.score,    locatorName: candidate.attributes.Loc_name };                 geom = candidate.location;                 graphic = new esri.Graphic(geom, symbol, attributes);                      map.graphics.add(graphic);                 return false;              }         });              doBuffer(graphic);         }                    ////////// buffer funcs ////////////////////////////////////////////////          function doBuffer(g) {         //map.graphics.clear();              bufferParams = new esri.tasks.BufferParameters();         bufferParams.distances = [100];                  bufferParams.geometries =  [g.geometry] ; // problem here?              bufferParams.unit = esri.tasks.BufferParameters.UNIT_FOOT;         bufferParams.outSpatialReference = map.spatialReference;         geometrySvc.buffer (bufferParams, showBuffer);     }          function showBuffer(geometries) {            var symbol = new esri.symbol.SimpleFillSymbol(         esri.symbol.SimpleFillSymbol.STYLE_SOLID,         new esri.symbol.SimpleLineSymbol(           esri.symbol.SimpleLineSymbol.STYLE_SOLID,           new dojo.Color([0,0,255,0.65]), 2         ),         new dojo.Color([0,0,255,0.35])       );            dojo.forEach(geometries, function(geometry) {         var graphic = new esri.Graphic(geometry,symbol);         map.graphics.add(graphic);       });               }


title: How do I overrideMimeType in Openalyers.Request.Get

text: I can't find any info on this. When performing a GET request with the following format, I cant figure out how to add the request for "overrideMimeType("application/xml")":               request = OpenLayers.Request.GET({     url: "http://aaaaaaaa.com",     async: true,     callback: request_handler,     });      An Openlayers request like this creates the request object and sends it all in this single code. I know if I weren't using the OpenLayers API, I would do something like:               request.overrideMimeType("application/xml");      But I can't get this to work using using my format. Any ideas? Thanks!