title: What is correct entry to Input Location (GPFeatureRecordSetLayer)*?

text: I use: http://sampleserver1.arcgisonline.com/ArcGIS/rest/services/Network/ESRI_DriveTime_US/GPServer/CreateDriveTimePolygons/execute What is correct entry to **Input Location (GPFeatureRecordSetLayer)** ? I try:               {"type": "Point","coordinates": [-105.01621,39.57422]}      but I get : "Invalid value for parameter 'Input_Location'"


title: Displaying multiple records in CartoDB with common zip code

text: I am uploading multiple records and want to use five digit zip code as the georeference. I would like to display each of those records on a map. Some of the records have identical zip code fields, meaning that several records share a single zip code. In the map view, only one of those records, among the handful that share a single zip code, appears. Is it possible to show on the map all of the records that share that single zip code?


title: gdallocationinfo is not recognized

text: I tried to run **gdallocationinfo** from my FWTools shell and got a message saying: "gdallocationinfo is not recognized as an internal or external command, operable program or batch file". So I checked the installation folder and didn't see any file with this name. So I am guessing that something could be missing. I've also tried to search for the file on Google but without success. I am using FWTools 2.4.7 installed on Windows 7. Any idea on how to fix this? Thanks


title: Clip and copy files to same directory structure using OGR/GDAL with Python

text: We need to clip all our data (many TB) by a new boundary and save it while maintaining the same data structure. I want to use OGR/GDAL as the data is in shp/tab/kml/geoTif/ecw etc. I can find the file and create the data structure in python but need to use ogr/gdal to then clip the file and save it to the new directory in the same projection and file type as the input. I have written the process to do this using arcpy but that doesn't work for tab and ecw, so I would prefer to go opensource but am new to this... So the python is as follows --               import os, glob, sys, shutil, datetime, subprocess     top = r'C:\Program Files (x86)\MapInfo\Offline Dataset\MapInfo_Tables\Property Information' # Source directory for all files to check and process     #top = os.getcwd()     RootOutput = r'C:\junk' # Directory to copy/move files to     RootDirDelete='C:\\' #repatative name of root file path to remove from new directory name     clip_polygon='c:\junk\boundary.shp'     filecount=0     successcount=0     errorcount=0          print "Working in: "+ top          list =[]     f = open(RootOutput+'\\Success_LOG.txt', 'a')     f.write("Log of files Succesfully processed. RESULT of process run @:"+str(datetime.datetime.now())+"\n")     f.close()               for root, dirs, files in os.walk(top, topdown=False):         for directory in dirs:             currentPath=os.path.join(root,directory)             os.chdir(currentPath)             #arcpy.env.workspace = currentPath             print os.getcwd()             lstFCs = glob.glob('*.*')             print lstFCs             for fc in lstFCs:                 OutPutDir=RootOutput+"\\"+str(os.getcwd()).strip(RootDirDelete)                 filecount=filecount+1                 if OutPutDir.endswith(".gdb"):                     print "File in GDB...SKIPPING"                 else:                     try:                         os.makedirs(OutPutDir)                         print "Created: "+OutPutDir                     except:                         #pass                         print "Unexpected error:", sys.exc_info()[0]          ##                try:                     fcIn=currentPath+"\\"+fc                     fcOut=OutPutDir+"\\"+fc                     ext=fc[fc.find('.')+1:]                     if ext=='shp':                         subprocess.call(["ogr2ogr", "-f", "ESRI Shapefile", "-clipsrc", clip_polygon, fcOut, fcIn], shell=True)                         successcount=successcount+1                         print "Data clipped: "+fcIn+" to " +fcOut                     else:                         print 'extension not shp'      * * * I need   * ogr clip function to clip data by clip_polygon for vector, raster and elevation -in arcpy it is like #arcpy.Clip_analysis(outPath+'\'+outFileN+'.shp', clip_features, outPath+'\'+outFileN+AOI+'.shp', xy_tolerance) -in ogr2ogr as a process I guess the following would work (Error when using gdal/ogr and python to clip a "shp" file by a "shp") #subprocess.call(["ogr2ogr", "-f", "ESRI Shapefile", "-clipsrc", clipping_shp, output_shp, input_shp], shell=True)   * projection should be same as input (usually EPSG 28355 or 4283) UPDATE: I ran the above process and it gives the following errors, the subprocess when run independently reports a 1 but no file is written...               [Dbg]>>> fcIn     'C:\\Program Files (x86)\\MapInfo\\Offline Dataset\\MapInfo_Tables\\Property Information\\shp\\TRC_DCDB.shp'     [Dbg]>>> fcOut     'C:\\junk\\Program Files (x86)\\MapInfo\\Offline Dataset\\MapInfo_Tables\\Property Information\\shp\\TRC_DCDB.shp'     [Dbg]>>> clip_polygon     'c:\\junk\\boundary.shp'     [Dbg]>>> subprocess.call(["ogr2ogr", "-f", "ESRI Shapefile", "-clipsrc", clip_polygon, fcOut, fcIn], shell=True)     1     [Dbg]>>>      [Dbg]>>> subprocess.check_call(['gdalinfo', fcIn], shell=True)     Traceback (most recent call last):       File "<interactive input>", line 1, in <module>       File "C:\Python27\ArcGIS10.1\lib\subprocess.py", line 511, in check_call         raise CalledProcessError(retcode, cmd)     CalledProcessError: Command '['gdalinfo', 'C:\\Program Files (x86)\\MapInfo\\Offline Dataset\\MapInfo_Tables\\Property Information\\shp\\TRC_DCDB.shp']' returned non-zero exit status 1     [Dbg]>>>       UPDATED                       if ext=='shp':                 subprocess.call(["ogr2ogr", "-f", "ESRI Shapefile", "-clipsrc", clip_polygon, fcOut,  fcIn])      gives me the following error.               >>>      Working in: C:\Program Files (x86)\MapInfo\Offline Dataset\MapInfo_Tables\Property Information     C:\Program Files (x86)\MapInfo\Offline Dataset\MapInfo_Tables\Property Information\shp     ['TRC_DCDB.dbf', 'TRC_DCDB.prj', 'TRC_DCDB.sbn', 'TRC_DCDB.sbx', 'TRC_DCDB.shp', 'TRC_DCDB.shx', 'TRC_DCDB_Raw.dbf', 'TRC_DCDB_Raw.prj', 'TRC_DCDB_Raw.sbn', 'TRC_DCDB_Raw.sbx', 'TRC_DCDB_Raw.shp', 'TRC_DCDB_Raw.shx']     Unexpected error: <type 'exceptions.WindowsError'>     extension not shp     Unexpected error: <type 'exceptions.WindowsError'>     extension not shp     Unexpected error: <type 'exceptions.WindowsError'>     extension not shp     Unexpected error: <type 'exceptions.WindowsError'>     extension not shp     Unexpected error: <type 'exceptions.WindowsError'>     Traceback (most recent call last):       File "C:\junk\Clip_And_Copy_Files.py", line 45, in <module>         subprocess.call(["ogr2ogr", "-f", "ESRI Shapefile", "-clipsrc", clip_polygon, fcOut,  fcIn])       File "C:\Python27\ArcGIS10.1\lib\subprocess.py", line 493, in call         return Popen(*popenargs, **kwargs).wait()       File "C:\Python27\ArcGIS10.1\lib\subprocess.py", line 679, in __init__         errread, errwrite)       File "C:\Python27\ArcGIS10.1\lib\subprocess.py", line 893, in _execute_child         startupinfo)     WindowsError: [Error 2] The system cannot find the file specified     >>>       From FWTools ogr2ogr in Python I tried               >>> command = 'ogr2ogr -f "ESRI Shapefile" ' + fcIn + fcOut     >>> os.system     <built-in function system>     >>> os.system(command)     1     >>> os.getcwd()     'C:\\junk'     >>> fcOut     'C:\\junk\\Program Files (x86)\\MapInfo\\Offline Dataset\\MapInfo_Tables\\Property Information\\shp\\TRC_DCDB.shp'     >>> fcIn     'C:\\Program Files (x86)\\MapInfo\\Offline Dataset\\MapInfo_Tables\\Property Information\\shp\\TRC_DCDB.shp'     >>>       but again I get no output..


title: Refresh Map after Invoke URL

text: I am making a map on **Auto Desk Infrastructure Map Guide Studio**. Now I created a small ASP.NET website and hosted in local machine, Now on click of different buttons in flexible layout toolbar, I **invoke URL** for that asp.net page and it opens fine. I do some data querying on that page by selecting different options and execute query. The data changes in table in database and the effect is shown on the layer on map. Ok the issue is that I want to refresh the Auto desk map after invoke URL window is closed. And I want to do it by configuration in Map studio 2013. Does anyone have any idea how to refresh map automatically after **invoke URL? Or in other sense if the data in the data source of layer is changed.** **Addtional Info:** the map always refresh when I click the refresh button from top or if I drag it or zoom it..Otherwise new layer data stays same.


title: how to add dynamic tree and layer in geoext

text: i am using Geoext,openlayer and Geoserve for map. Now i want to create dynamic tree and layer using Geoext. Means any user can create a tree and add its layer (child) to that tree dynamically. I have seen many examples on geoext website tree panel like http://api.geoext.org/1.1/examples/tree.html.But this a static example. i want to that user can add any number of tree and it layer dynamiclly.It is possible by Geoext and openlay,If there is any other solution for it please suggest. How can i do it Please help me. Thanks in Advance


title: Unable to display legends of group layer in openlayers

text: i m using the example http://api.geoext.org/1.0/examples/tree.html but how can i show legends of tasmaia (group layer) in legend panel here is my working application link http://hsdi.in:8080/geoserver/www/webgis/index.html


title: How to rename ogr2ogr tables on export to Spatialite?

text: I'm exporting data from PostGIS to Spatialite using ogr2ogr. The export is working properly, except I need to be able to rename the tables that are being exported. Here's the command I'm using: `ogr2ogr -f "SQLite" testSQLite.sqlite -progress PG:"host=127.0.0.1 dbname=mydb user=admin password=myPassword active_schema=my_schema schemas=my_schema tables=building_shared,alarm_shared" -lco LAUNDER=yes -dsco SPATIALITE=yes -lco SPATIAL_INDEX=yes -overwrite` This exports a Spatialite database with two tables - **building_shared** and **alarm_shared**. I want the two tables to be named **building** and **alarm**. I tried adding the following SQL query to the end of the command, with no luck: `-sql "SELECT * FROM building_shared AS building"` How do I rename the tables on export?


title: AutoDesk Infrastructure Studio - How to enable POST for WFS services?

text: I'm attempting to get feature info with the following url: SERVICE=WFS&REQUEST=DescribeFeatureType&VERSION=1.1.0&TYPENAME=ns75382197:MAPSHEET_GRID With many, many typenames listed after this first one. Is there any way to enable the POST service analogous to this GET request?


title: What is the best way to join lots of small polygons to form a larger polygon?

text: I have the following layer using SRID 27700 in postgis: ![enter image description here](http://i.stack.imgur.com/5P53M.png) It's every administrative region in the UK, and (as you can see from the colour grouping) each of them has a text field specifying the county they lie in. What I'd like to do is to make larger county polygons from the smaller ones in a given county, so EG in the picture above all the teal colour polygons would form one large polygon from the single outer ring that contains all the polys in that colour, like wise all purple, brown, pink, grey etc should all form one polygon. I've already tried the following:               insert into parishesmerged (geometry)     select astext(multi(ST_Union(the_geom))) as the_geom from parishes     group by county_name      But it keeps generating broken geometries which I then have big problems processing further. If anyone has any further Ideas it would be appreciated. I'm essentially trying to make a simpler county level map with the major output areas in. Any solutions don't have to be in Postgis either, I have the full OS4Geo stack installed, the latest version of QGis and more utils than I can shake a stick at. The only things I don't have are the big boys like ArcGis (Although I may have an Old Mapinfo lying around somewhere) Cheers Shawty ------======== Update : 22/8/2012 ========------ Thanks everyone for all the answers so far, I'll go away and examine each one, see what happens then come back and document on my findings. For the record, the dataset I'm trying to create is to accompany a GIS book I'm on writing aimed at .NET programmers wishing to write GIS applications using .NET I know what I need to know about relatively straight forward stuff in Postgis such as Spatial Querying and Statistics gathering and if the publisher will allow me I do intend to mention gis.stackexchange and all those who frequent here :-) I'll be back soon.... Cheers Shawty ------======== Update : 25/8/2012 ========------ Well after wandering back off into cyber space and trying the suggestions below, the one that worked the best was 'Paul Ramseys' solution. Thanks Paul. I now have a nice simplified counties & boroughs file that's just simple enough for my book, but complex enough to allow me to demonstrate some interesting geo-spatial SQL. Even though Paul's solution ultimately was the one that worked for me, I also drew on the other answers for things like simplifying the polygon map and reducing the complexity further. On thing I did observe while doing this however, while ST_Collect is indeed faster than ST_Union, run for run it was also the one mostly responsible for broken geometries. My guess is the speed increase is obtained at the expense of less accuracy in the core function.


title: Feature Class - Adding domain to have dropdown box in Attribute Table

text: http://resources.arcgis.com/en/help/main/10.1/index.html#/as_a_stand_alone_feature_class/002200000005000000/ Line 23 - To associate a domain with this field, click the field next to Domain, click the drop-down arrow to see a list of the domains that apply to this field type, then click the domain. i am trying to get where, in the attribute table, I have a dropdown box. I follow the instructions in here and I am getting knowhere. What am I not doing, is there more to do then what this states? I am using ArcView 10.1 I go to the FEATURE CLASS PROPERTIES, under FIELDS, I have the field name "STATUS" and text. Then I go to SUBTYPES. The subtype field shows , what do I need to add to that field? I can't click anywhere in those boxes. I am just looking to have the dropdown boxes under status "producing", "leased" etc... Thanks, Z ![enter image description here](http://i.stack.imgur.com/6atQ5.jpg) ![enter image description here](http://i.stack.imgur.com/sUQ0q.png) ![enter image description here](http://i.stack.imgur.com/Xksuv.jpg) ![enter image description here](http://i.stack.imgur.com/7wTas.jpg)


title: filepath length errors in ArcMap

text: Is there a maximum directory or char length for file paths in ArcMap 10.0? I spent hours trying to trouble-shoot an untyped crash in both ArcMap and Catalog when trying to reclassify and subsequently delete a large raster. Moving the files to a shorter directory seems to have fixed the issue, but now I am curious if there is a certain filepath length you cannot exceed.


title: Type mismatch error in Calculate Field tool in ArcGIS 10.1

text: I'm trying to run a rather lengthy model (built in ArcGIS 9.3.1) in ArcGIS 10.1 and I'm recieving an error on a Calculate Field tool. The error is - `Type mismatch: 'iif'`. This error is refering to an expression in the tool: `iif(IsNull( [PID] ) ,NA,[PID])`. ![enter image description here](http://i.stack.imgur.com/jNIYi.jpg) I tried "if" instead of "iif" but then recieved a generic error. I never had this issue when running the model in 9.3.1. I'm not sure what to do here. Any suggestions?


title: Can not get text element to update?

text: I can not get text element to update? Text ELement has been added "Subdivision" Still not updating?..               class ButtonClass2(object):     """Implementation for FemaMapper_addin.button (Button)"""     def __init__(self):         self.enabled = True         self.checked = False     def onClick(self):         #set the map document reference         mxd = arcpy.mapping.MapDocument("CURRENT")         #set the dataframe reference using the map document and the first data frame in the list         df = arcpy.mapping.ListDataFrames(mxd)[0]         #set the layer reference using the map document reference and the second layer in the list of layers         layer = arcpy.mapping.ListLayers(mxd)[1]         #heading of the map         map_title = "Parcel Map"         #create a search cursor to get relevant information about a parcel         sc = arcpy.SearchCursor(layer)         #for loop to get the selected data         for row in sc:             #select the relevant fields to be put on the map             Sub = row.getValue("SUBDIVISIO")                   #text element of the map         elm = arcpy.mapping.ListLayoutElements(mxd, "TEXT_ELEMENT", "Heading")[0]         elmSub = arcpy.mapping.ListLayoutElements(mxd, "TEXT_ELEMENT", "Subdivision")[0]              #font of the heading         elm.text = '<FNT name ="Arial" size="32">' +  map_title + '</FNT>'         elmSub.text = '<FNT name ="Arial" size="14">' + 'Subdivision: ' + Sub + '</FNT>'              #Change the view to layout view         arcpy.mapping.MapDocument("CURRENT").activeView = "PAGE_LAYOUT"         #refresht the active view         arcpy.RefreshActiveView()         #delete search cursor and row variables         del sc         del row


title: How to configure the Ubuntu server after publishing the map using Geoserver (Open GeoSuite)

text: I've installed the **GeoServer** that comes with the **Open GeoSuite** on **Ubuntu**. I've created an HTML page (index.html) inside the /usr/share/opengeo-suite-data/geoserver_data/www/ for viewing maps using OpenLayers. index.html is working fine in local directory. To deploy my Web Map Server on the Web Server environment, how could I configure my Ubuntu Server? Many thanks in advanced.


title: OpenGeo Suite Dashboard Not Showing Up

text: I installed the OpenGeo Suite in my computer running Windows 7 Home Edition. When I installed it, I checked the box "launch OpenGeo Suite" and it ran. After installing several other software (XAMPP, Hamachi), the Suite Dashboard won't show up. As in there's no interface or taskbar icon. However, I checked the services tab in Task Manager, it was there and running. I tried changing the ports used by XAMPP, Skype, YM, and Hamachi, but the interface isn't there still. Any help is very much appreciated. Thank you.


title: PostGIS: Querying a record table and explanatory table

text: Basically I have two tables one containing a collection of records and one table that contains extra info about features in the first. i.e Records Table: GID SITE SPECIES GEOM and the second table has GID SITE SITE_desc Essentially the first table has many records for the same site while the second table only has one record (in one row) for each site with a description. What i'm wondering is if there is anyway i can query both tables (for a specific species for example) and have the field Site_desc be printed in a new column for each record where the sites are the same (SITE=SITE). Is this possible? Regards Paul


title: E00 file structure?

text: I am trying to write a java application that reads and visualizes data from E00 file. I cannot find any sources describing what is the structure of an E00 file and also what information does it contain exactly, line by line.


title: Select layer by attribute: return top 3 records

text: I have a feature class in a file geodatabase using ArcGIS 10.1. I want to create a model that selects the three features with the highest values in one field (Average Price). Basically I want my selection to include the three features with the highest average prices.I cant quite figure out how to accomplish this. Any suggestions would be greatly appreciated!


title: Placing Pins and following transformation on zoom map d3

text: So I am using this http:// bl.ocks. org/mbostock/4699541 framework on D3.js and I just figured out how to append circles to the map. However when you click on a state the circle moves out of place and I can't. Here is my code for appending circles:               d3.json("city_locations.json", function(new_us_places){     var my_cords = projection(new_us_places["CA"]["SAN DIEGO"]);     var my_cords1 = projection(new_us_places["CA"]["SAN FRANCISCO"]);     var my_cords2 = projection(new_us_places["CA"]["SAN FERNANDO"]);     var my_cords3 = projection(new_us_places["CA"]["LOS ANGELES"]);     var r = 5;     svg.append("circle")     .attr("cx", my_cords[0])     .attr("cy",  my_cords[1])     .attr("r", 5)     .attr("class", "HighlightCircle")     .on("click", function(d){alert("text")})     svg.append("circle")     .attr("cx", my_cords1[0])     .attr("cy",  my_cords1[1])     .attr("r", 5)     .attr("class", "HighlightCircle")     .style('fill','red')     .on("click", function(d){alert("text")})      Here is what it looks like now: http://i.imgur.com/02vqjA2.png This is what happens when you click a state to zoom: http:// i.imgur. com/rPd4HdP.png WHAT I want: (and hopefully able to keep that marker in place when click state to zoom): http://i.imgur.com/YUrZ0Kr.png edit: sorry for links with spaces, don't have enough rep for more than 2 links...


title: Using a point to cut one polyline

text: Given a polyline L and a point P which belongs to this polyline but that also is not necessarily one of the polyline's vertices, I'd like to measure the length from the start of the polyline L till P. I got P using an identify task with zero tolerance, so I can assume for right that the point belongs to the polyline. To reach my goals, I create a fake polyline L' that contains the points P and P', where P'=(Xp + Delta, Yp + Delta). Then I use L' to cut L and get two new polylines L1 and L2. Finally I use one geometry service to get the length I want. The code snippet to create the fake polyline, cut L and measure the length:               private void IdentifyCompleted(object sender, IdentifyEventArgs e)     {         if (e.IdentifyResults.Count <= 0)         {             return;         }              var geometryService = new GeometryService(                 @"http://myarcserver.com/ArcGIS/rest/services/Geometry/GeometryServer");         MapPoint pLine = new MapPoint();         // Some workaround here: creating and 'arbitrary' line starting with P         pLine.X = p.X + 0.1;         pLine.Y = p.X + 0.3;              PointCollection pointCollection = new PointCollection();         pointCollection .Add(p);         pointCollection .Add(pLine);              Polyline cutter = new Polyline();         cutter.Paths.Add(pointCollection);         cutter.SpatialReference = Mapa.SpatialReference;              // Transform the polilyne returned from the identify task in a list         List<Graphic> polylineL = new List<Graphic>();         foreach (var result in e.IdentifyResults)         {             polylineL.Add(result);         }              geometryService.CutCompleted += geometryService_CutCompleted;         geometryService..CutAsync(listaResultado, tesoura);     }          private void geometryService_CutCompleted(object sender, CutEventArgs e)     {         if (e.CutIndexes.Length == 2)         {             // Magically figured out what are the first geometry. It             // is out of the scope of this question             Polyline first = getFirstGeometry();             int theLengthIWant = first.Length();             // We got it ! ! !         }     }      I guess is too overload execute one identify task and then one cut operation. So my question is: is there a simpler to do this? I mean, is there a straight forward way to get the length from the start of one polyline till an arbitrary point on it? [EDITED] I found a lot of people using IHitTest to do what I want. Unfortunately, this interface belongs to ArcObjects API which is not available to me right now.


title: Updating records in a table using GetCount and an Update Cursor

text: What I have is a buffer which is then used to select a number of point features which intersect the polygon. I would like to take that number and add the value to a table. The row where the value will need to be located corresponds with the name of the original feature class the buffer was created from.               import arcpy     from arcpy import env     import os     from os.path import join                    arcpy.env.workspace = r'MyWorkspace'     updateTable = r'table'     shortField = "FieldName"     arcpy.MakeFeatureLayer_management("FieldName", "fieldname_layer")     fcs = arcpy.arcpy.ListFeatureClasses("*_Value", "POLYGON")     for fc in fcs:         select = arcpy.SelectLayerByLocation_management("fieldname_layer", "INTERSECT", short)         result = int(arcpy.GetCount_management(select).getOutput(0))         print short.replace('Buffer_','') + " " + str(result)      So far the best I can come up with is printing the file name less the Buffer_ portion which precedes the name I need. I cannot wrap my head around how to go from there to updating only relevant records using the update cursor. thank you for any insight you are able to provide.


title: GRASS v.out.dxf 3D into AutoCAD

text: So continuing with my last question (Create 3D dxf with OGR) I have decided to give the GRASS option a go. I have imported my polylines into GRASS using v.in.ogr and converted them to 3D using v.to.3d. I checked the lines using the GRASS info tool and the polylines all have the correct z values. I then exported it using v.out.dxf but this is where I am having troubles. The dxf contains the z values for the polylines: _Extract_               VERTEX     8     line     10     401120.708496     20     6875185.007069      30     539.000000     0       However when I open it in AutoCAD Civil 3D it doesn't recognize the z value and just makes all the lines 2D. I have checked them with the info tool in AutoCAD but the z values are all the same. Is the dxf created using v.out.dxf supported in AutoCAD in order to get correct 3D?


title: QGIS: using postgis can I define a layer using sql?

text: I have a vector layer (call it master) that is stored in postgis and consists of a number of line segments. I also have a collection of sub layers that consist of subsets of the master layer. I can set up a 'link' table master_id, layer_name the says which segments belong to which sub layer -- some on the edges will be common. Can I specify a layer using an join on the master and link tables? Alternatively can I use a view at the postgres level to do this. I have never used views or stored procedures although I know the general principles.


title: Is it possible to create a standalone python script to export "layout" to pdf?

text: I have a python script that successfully runs from ArcMap 10.2. I was wondering if it would be possible to rewrite it so that it is a standalone script that prompts the user for the shapefile and then runs some analyses (e.g. buffer, tabulate intersection, intersect, etc.), changes layer's symbology properties using a preexisting lyr file, zooms to its extent, turns on and off other layers in the mxd TOC template depending on a specific conditions, changes the map title, and finally exports resulting map as a pdf. The script is quite extensive and although much of the coding can be run as a standalone script, I am wondering if the section of the script dealing with manipulating the layout properties (legend, title) and the layers in the TOC as I want them to look so that I may create a map for exporting as a pdf can be done as a standalone script. It seems to me that some of these actions may only be possible inside an ArcMap session.


title: WMS GetMap MIME type error

text: After adding the following to the top of my index.html page                   <!DOCTYPE html>      I begin to receive the following error on my GetMap WMS request                   Resource interpreted as Image but transferred with MIME type application/vnd.ogc.se_xml      Any thoughts on why this would occur are appreciated. Here are the URL params of my WMS request:                   BBOX=-15130249.425402,2538320.1060898,-5869750.5745982,6461679.8939102         FORMAT=image/png         HEIGHT=802         LAYERS=EFETAC-NASA_1YrBaseline_previous1         MAXEXTENT=NaN,,,         PROJECTION=EPSG:900913         REQUEST=GetMap         SERVICE=WMS         SRS=EPSG:900913         STYLES=         TRANSPARENT=true         UNITS=m         VERSION=1.1.1         WIDTH=1893      Thanks, Derek


title: Error output with gdalwarp

text: Following is my raster file and mask shape file files I use qgis 1.8.0 to clip the raster with shape file by gdal tools but after I run warp in gdal tools it's just give me no data in the output GeoTiff?


title: Is there a way to refresh imported modules in a ArcGIS Python Toolbox?

text: When using a python toolbox (.pyt) in ArcMap, I'll typically follow a pattern where the .pyt file itself is simply a wrapper to collect input arguments and define the tools themselves. Supporting code is kept in separate unit-testable modules. Example:               import supporting_module          class MyTool(object):              ...              def execute(self, parameters, messages):             """The source code of the tool."""             some_input = parameters[0].valueAsText             some_output = parameters[1].valueAsText                  supporting_module.do_something(some_input, some_output)             arcpy.SetParameter(2, some_output)      This works out really well but I've run into one frustrating issue during development. Right-click -> Refresh on the .pyt in ArcMap only refreshes the .pyt file's code. It does not refresh the imported modules, so I have to close and re-open ArcMap whenever I change something there. Fortunately, since I'm testing the code independently, I don't have to do this a ton, but it's still a major hassle. Is there any way around this? Somewhat related - is there any way to fully refresh the Python console (i have a custom site-package that I must also close/reopen ArcMap to pull in changes from as well)? I'm using ArcMap 10.2.1.


title: Import SketchUp created KML contour lines model to QGIS

text: There's a neat trick to create contour lines based on the GE 3D terrain using Sketch-Up. The resulting KML however is a model ![enter image description here](http://i.stack.imgur.com/lyDD7.png) I have looked at ways of importing KML layers to QGIS, like it was suggested here, yet it doesn't accept this kind of kml file. So, is there a way to import such a thing into QGIS (or ArcMap for that matter)? and if not can such a model be split into its separate contour lines while retaining the z elevation attribute? Here's the model itself: ![enter image description here](http://i.stack.imgur.com/Rk7Yd.png)


title: Opinion of using Google Earth Pro vs ArcGIS Explorer

text: I would like to hear from anyone who uses either or both Google Earth Pro vs ArcGIS Explorer. The reason I am asking is because our agency has a very limited of budget to purchase license for each person. They are thinking about moving the license to centralize for everyone to use it. GEP cost agency money while ArcGIS Explorer is free and download. One of my co-workers wants to use GEP for his ecological site work and wants to use photos and share others. We have soil scientists here that will probably need to use GEP because they wanted to our soil line to overlay with the NAIP imagery that GEP uses it. One other possible is that they will use the soil line to edgematch with other soil lines from adjacent states. What is the pro and con for using either program ?


title: How to convert UTM coordinates into Decimal Degree coordinates?

text: I work on QGis 2.0.1 - Dufour on Macintosh OS X Mavericks. I have UTM coordinates (XXXX and YYYYY) in a table (with one column of X and an other one of Y). I need this coordinates in Decimal Degree. Could you tell me the detailed protocol to follow please ?


title: calculating a field in which null values may be present

text: I would like to calculate a field (short integer) using other fields (short integer) present in the table. In some cases these fields may have null values for one or more of the fields. For example, the new field name is 'total' and total is based upon adding field1 + field2 + field3. Some records (rows) have values for field1, field2, and field3, but some rows only have values for field1 and field2, then field3 is null. Do I need to go and replace all null values with a numeric value? Or, is there some workaround within the VB or python script? I'm using ArcMap 10.1 Thanks -al I've used a python example from below and have been unable to get it to work. Here is the code I'm using:               #calculate new field 'summed' in feature class 'test'     inFeatures='test'     fieldName01='summed'     exSummed = "stack(!field01!+!field02!+!field03!)"     codeBlock = "def stack(*args):         return sum(filter(None, args))"          arcpy.CalculateField_management(inFeatures, fieldName01, exSummed, "PYTHON", codeBlock)      where 'summed' is the field (short) I would like to populate and field01, field02, field03 are all short integers which may or may not have NULL's present. I receive an error when I run this... ![enter image description here](http://i.stack.imgur.com/xMFmI.png) I'm sure it's something simple I'm doing wrong. Thanks-al


title: Merge Multiple Tables into a New Table in PostGIS

text: I am looking to merge a number of individual tables into a new table in PostGIS. This is an easy task when working with Shapefile data, but I am unsure how to do this in PostGIS. Any help would be greatly appreciated. I think I use Append, but I am unsure even where to start. Thanks in advance.


title: Offset Google markers at the same address

text: What can I incorporate into my current code below to offset markers, as markers at the same address don't show. Here's example http://www.mediwales.com/mapping/test UPDATE: I've changed the way the map works by running it from a database.


title: is it possible to download planet files of OpenSeaMap data?

text: I want to visualize OpenSeaMap data on an Android device using Mapsforge. Since I need to do this in a disconnected manner (think boat out to sea) I need to generate a .map file using the Mapsforge Osmosis plugin, however to be able to do that I need a .osm file of the OpenSeaMap data. Can anyone point to where I might find one or some other way to proceed?


title: ArcGIS Model Builder perform split operation on each selected feature and delete identical on collection

text: I have created a model to select parcels and create segments for each parcel. If two parcels are sharing the same segment, delete duplicate. My model is attached. My question is: Before appending or deleting identicals, I need to perform the action on a whole. The model I created is executing for each feature. Please advice..! ![enter image description here](http://i.stack.imgur.com/1BCyH.png)


title: Python code does not work with RasterCalculator after exported from Model Builder

text: I have a model authored in Model Builder that uses "RasterCalculator" to process input raster datasets.                arcpy.gp.RasterCalculator_sa("OutRas = Con(\"%STATIC_forestland%\"  ==  0, -255, Con(\"%mapvictoria20121210%\" >= 115, -255, Con((\"%mapvictoria20121210%\" >100) & (\"%mapvictoria20121210%\" < 115), 100, Con(\"%mapvictoria20121210%\" <= -15, -255, Con((\"%mapvictoria20121210%\" < 0) & (\"%mapvictoria20121210%\" > -15), 0, \"%mapvictoria20121210%\")))))", mapvic1)      It works well as a model. However, when exported to a Python script and saved in the same Toolbox, the script does not work. The error message is as below               RuntimeError: ERROR 000732: Input Raster: Dataset %mapvictoria20121210% does not exist or is not supported      Failed to execute (RasterCalculator). I am sure the input raster are there. Can anyone give some help? Is it to do with raster naming issues? Thanks in advance. Cheers, Alex


title: opengeo-suite upgraded from 3.0 to 3.1 on Ubuntu and I can't start GeoServer

text: I'm having trouble starting GeoServer after updating all of the opengeo-suite applications via apt-get update, apt-get upgrade and apt-get dist-upgrade. GWC works, the dashboard is there as is geoexplorer. I can see geoserver in Tomcat's application manager, but I get a 404 tomcat error when I go to /geoserver. Log file is here: http://pastebin.com/TdnYRxAZ catalina.properties here: http://pastebin.com/acJCpAsa


title: Looking for a well used opensource .net WMS server that can show raster layers

text: We covered a large area with a plane and have all the images with georeferences avalible also. I would like to provide the images taken as a WMS and not knowing much about gis still, there are alot of options out there. I hope someone here can direct me to a good candidate to start developing against. I found something called MapArround and a .net example using it but that was only showing shape files. If anyone has an example in MapArround leveraging custom images that would be great. Otherwise, then any example of leverage our custom taken photos. We already have software that sticth images together and having overlap so basicly for every request we have the images needed to be provided. I took the time to look at the WMS specification and know we could implement the requests our self but if we can save time but just providing an extension into some already opensource wms solution that is preferable.


title: Minimum records in an excel for converting shape file

text: I'm Converting an excel with 3000 rows and 10 columns. My Problem is while exporting excel to a shapefile resulting that " **The Maximum Record length has been exceeded** ". Plz Suggest..


title: How to make a raster layer transparent via MapServer

text: I loaded a PostGIS raster Data through OpenLayer and Mapserver. The code inside the OpenLayer is as following :                localWMS = new OpenLayers.Layer.WMS( "Local WMS layers",                       "http://127.0.0.1:8080/cgi-bin/mapserv.exe?MAP=C:/MapServer/raster.map",                       {layers: 'polygon',transparent: true,alpha:true},                       { projection:new OpenLayers.Projection("EPSG:4326"),                         units:'m',                       isBaseLayer: false}                       );           map.addLayer(localWMS);      And Here is my Map File, I already set opacity to each class of the layer but nothing happens.               LAYER     DATA "PG:host=*** dbname=*** schema=public table=london where='rid=3' mode='1'"     NAME "postgis"     PROCESSING "SCALE=AUTO"     PROCESSING "NODATA=-9999"     STATUS DEFAULT     TYPE RASTER          UNITS METERS     CLASS       NAME "Minor Impact"       EXPRESSION ([pixel] >= 0 AND [pixel] < 1.29525335532 )       STYLE         ANGLE 360         OPACITY 50         COLOR 122 182 245         SYMBOL 0         WIDTH 1       END     END     CLASS       NAME "Moderate Impact"       EXPRESSION ([pixel] >= 1.29525335532 AND [pixel] < 2.59050671063 )       STYLE         ANGLE 360         OPACITY 50         COLOR 76 230 0         SYMBOL 0         WIDTH 1       END     END     CLASS       NAME "Severe Impact"       EXPRESSION ([pixel] >= 2.59050671063 AND [pixel] < 5.69911476339 )       STYLE         ANGLE 360         OPACITY 50         COLOR 252 228 91         SYMBOL 0         WIDTH 1       END     END     CLASS       NAME "Catastrophic"       EXPRESSION ([pixel] >=5.69911476339 )       STYLE         OPACITY 30         COLOR 232 16 20       END     END     END      The map looks like following :![enter image description here](http://i.stack.imgur.com/FSoOc.png) Then, I tried to set one of the opacity to less than 1, the red color just disappeared. ![enter image description here](http://i.stack.imgur.com/2Rp5K.png) But if I set the Opacity of the Layer to 50               LAYER     DATA "PG:host=*******************'"     NAME "postgis"     PROCESSING "SCALE=AUTO"     PROCESSING "NODATA=-9999"     STATUS DEFAULT     TYPE RASTER     OPACITY 50     UNITS METERS      Even more strange map appeared. ![enter image description here](http://i.stack.imgur.com/f2yXW.png) Does anyone know what's the problem about that? Or Does anyone know the proper way to set a transparent raster layer via MapServer ? Thanks in advance.


title: How to calculate the length of clipped lines using buffered points using ArcPy?

text: I need to create a buffer around each point of a points feature class. I would then want to use the point buffers to clip the streets feature class to calculate the length of the clipped streets using ArcPy. I need some guidance in summarizing the length of the streets regarding their types. The "summarizing" in the following script does not work, I receive the following message:               ERROR 000732: Input Table: Dataset roads_0 does not exist or is not supported     ERROR 000800: The value is not a member of SUM | MEAN | MIN | MAX | RANGE | STD |   COUNT | FIRST | LAST.     Failed to execute (Statistics).      My ArcPy script so far:               import arcpy          arcpy.env.workspace = "D:/Data/ToYahya"          # Perform the buffer which will create one circular buffer per site     arcpy.Buffer_analysis("test_point.shp", "Buffer_test.shp", "50 Meters", "FULL",   "ROUND", "NONE", "")          # Open a search cursor and use the geometry of each buffer polygon in turn to Clip roads     # Use the id of each buffer polygon as suffix on the out road shapefile     features = arcpy.SearchCursor("Buffer_test.shp")     for feature in features:           bufferPoly = feature.SHAPE           arcpy.Clip_analysis("road.shp", bufferPoly, "road_" + str(feature.FID))                arcpy.Statistics_analysis("road_" + str(feature.FID), "outStatsTable", "Shape_Leng", "type")     del feature,features


title: Different legend properties with ArcGIS API for JavaScript

text: I have a map template that I use to display a variety of feature layers and dynamic layers. I have found a snippet of code to add the legend to my map:               //add the legend     dojo.connect(myMap, 'onLayersAddResult', function (results) {         var layerInfo = dojo.map(results, function (layer, index) {             return { layer: layer.layer, title: layer.layer.name };         });              if (layerInfo.length > 0) {             var legendDijit = new esri.dijit.Legend({                 map: myMap,                 layerInfos: layerInfo},                 "legendDiv");                 legendDijit.startup();             }         });      The issue is, I can't figure out how to configure the names of the layers, labels, descriptions, etc. Most of the layers in my service are symbolized by Category, so, I want to be able to display those category names in the Legend. This is what I want to show..or something like it, but with a category or layer heading ![](http://i.stack.imgur.com/owQEM.png) This is what shows up. No heading, just the categorized layer without a title or anything. ![](http://i.stack.imgur.com/38qtO.png) Any help would be appreciated. Thanks


title: Script or model to create random points having multiple minimum distances within a polygon

text: Given: A single part polygon, total number of points to place inside that polygon, and table specifying the size of each subset (number of subsets n >= 1) of that total number of points and each subset's minimum spacing. Find: A (semi-)random distribution of points inside the polygon such that no point is closer than its minimum spacing to any other point in the polygon. For example, you are given a polygon and 100 points. Of the set of 100 points there are three disjoint subsets--A, B and C. Subset A contains 50 points having a minimum spacing of 7 m, subset B contains 25 points having a minimum spacing of 4 m, and subset C contains 25 points having a minimum spacing of 2 m. Distribute the points such that no point in subset A is within 7 m of any of the remaining 99 points, no point in subset B is within 4 m, no point in subset C is within 2 m. One possible method to solve this would be to start with subset A, distribute the points randomly using Create Random Points, create a 7 m buffer around the points, erase those buffers from the original polygon and use the erased polygon to randomly distribute the points in subset B. . .repeat through all subsets. I'm looking for a way to implement this (or a superior) solution using scripting or model builder. Ideally, the user would input the number of points in each subset, the minimum spacing for each subset, and the polygon feature class into which the points would be placed.


title: Buffer is giving me a different distance

text: I'm using ArcGIS 10.1. When I Buffered I changed the environments to compare the two different projections. WyLAM is a statewide Lambert Conformal Conic projection used for Wyoming. State plane in Wyoming is broken into four zones and looks like it buffering is more accurate in a smaller zone. I have a well pad that needs to be outside of a 4 mile buffer in Wyoming. I created one that is 4 miles in Wyoming East. Another person created one in wyLAM. I wanted to know why it was giving two distances.


Create Feature and store to different FGDB I'm creating a feature and I want to store it to a specific geodatabase. At the moment I'm using feature.Store() but it's not storing it to where I want. How do I specify a location? I want to store it to C:\College Stuff\Dissertation\2009 stuff\Geo\Geo\Full_3D.gdb I've been told that if I am looking to store in a particular database I will need to open the feature class from that database. I've looked at connecting to a geodatabase using this code                   Public Shared Function FileGdbWorkspaceFromPath(ByVal FGDB_WORKSPACE As String) As IWorkspace         Dim factoryType As Type = Type.GetTypeFromProgID("esriDataSourcesGDB.FileGDBWorkspaceFactory")         Dim workspaceFactory As IWorkspaceFactory = CType(Activator.CreateInstance(factoryType), IWorkspaceFactory)         Return workspaceFactory.OpenFromFile(FGDB_WORKSPACE, 0)     End Function      This is my CreateFeature sub                Public Sub CreateFeature(ByVal layername1 As String, ByVal x As Double, ByVal y As Double, ByVal FeatureClass1 As IFeatureClass)              Dim cPoint As IPoint = New ESRI.ArcGIS.Geometry.Point         ' To individually assign a value to the coordinates, use cPoint.X = x and cPoint.Y = y.         cPoint.PutCoords(x, y)                   ' Ensure the feature class contains points.         If FeatureClass1.ShapeType <> esriGeometryType.esriGeometryPoint Then             Return         End If              ' Build the feature.         Dim feature As IFeature = FeatureClass1.CreateFeature()         feature.Shape = cPoint                        ' Commit the new feature to the geodatabase.         feature.Store()          End Sub      To create the feature do I need to connect to the geodatabase in the CreateFeature sub? How do I do that? I'm modifying someone elses code and I am new to this so that's why I am very unsure Here is the GetFeatureClass function, this where the CreateFeature gets the feature from. GetFeatureClass gets a feature class from the active map. The feature class is a single shapefile, I want it to get a feature class that is in the active map but in a geodatabase               Public Function GetFeatureClass(ByVal layername1 As String) As IFeatureClass         Dim mxDocument As IMxDocument         mxDocument = m_application.Document         Dim map As IMap = mxDocument.FocusMap         Dim check As Boolean = False              Dim layerCounter As Integer         Dim FeatureClass1 As IFeatureClass = Nothing         Dim featureLayer As IFeatureLayer         Dim currentLayer As ILayer         For layerCounter = 0 To map.LayerCount - 1 Step layerCounter + 1             currentLayer = map.Layer(layerCounter)                  If (TypeOf currentLayer Is IFeatureLayer) And (currentLayer.Name = layername1) Then                           featureLayer = currentLayer                      FeatureClass1 = featureLayer.FeatureClass                      check = True                      layerCounter = map.LayerCount             End If         Next              If check = False Then             MsgBox("cannot find InputStops layer")             Return Nothing         End If              Return FeatureClass1     End Function


title: cross domain problem with :OpenLayers.Protocol.HTTP

text: this code works with local kml files, but not when the kml files are on my other server, I got 3 error :               OPTIONS http://192.168.0.150:9000/files/anonymous/AVGHDOP.kml 404 (Not Found) XMLHttpRequest.js:220     OPTIONS http://192.168.0.150:9000/files/anonymous/AVGHDOP.kml Origin http://localhost is not allowed by Access-Control-Allow-Origin. XMLHttpRequest.js:220     XMLHttpRequest cannot load http://192.168.0.150:9000/files/anonymous/AVGHDOP.kml. Origin http://localhost is not allowed by Access-Control-Allow-Origin.  so, I try to put headers in the OpenLayers.Protocol.HTTP. but it still is not working -_-. here is my code:               map = new OpenLayers.Map('map');     osm = new OpenLayers.Layer.OSM();     vectors = new OpenLayers.Layer.Vector("KML", {         strategies: [new OpenLayers.Strategy.Fixed()],         protocol: new OpenLayers.Protocol.HTTP({             url: url_file,//the file is on an other_server(a play server).             headers: {                 "Access-Control-Allow-Origin":"*",                 "Access-Control-Allow-Headers":"Access-Control-Allow-Origin"},             format: new OpenLayers.Format.KML({                 extractStyles: true,                 extractAttributes: true,                 maxDepth: 2             })         })     });     map.addLayers([osm, vectors]);     map.addControl(new OpenLayers.Control.MousePosition());     map.addControl(new OpenLayers.Control.EditingToolbar(osm));     vectors .events.register("loadend", vectors , function (e) {         map.zoomToExtent(vectors .getDataExtent());     });      I dont know if I can configure a proxy. So if I can solve this problem without it it would be great :)


title: Should I delete the existing replica when create a new replica with the same name?

text: I have successfully create a replica between 2 geodatabase and use a scheduler to synchronize the replica between the 2 geodatabase several months ago. But the scheduler didn't work because of changes in the original geodb. In this case I need to create the replica once again (according to user's request), my question is if I use the same replica name, is there a need of deleting the existing one? What will happen if I create the replica directly? Any potential dangers? Since this is a in operation system, I cannot afford of failure. EDIT: After my research from internet, I even cannot find a way to delete a replica.


title: QGis 2.4 - Python 2.7 PyQT4.Qtcore.QStringList Import failed

text: I receive an import error when trying to import QStringList from the Qgis Python console (or when loading a custom plugin).               from PyQt4.QtCore import QStringList     Traceback (most recent call last):      File "<input>", line 1, in <module>     ImportError: cannot import name QStringList      I receive no error when using the same import from an ipython console or executing a standalone python script. Python 2.7.5+ [GCC 4.8.1] Please advise


title: How to Create SLD Label with Non-Latin Languages

text: I create a OpenLayers WMS Layer with my map server url. I want to Style this with SLD_Body. I defined My SLD_Body and want to labeling it. but My Label Map is with non-Latin language and after executing the program , the labels is shown with Question mark(????). Is there any way to display labels in My language?? this is my code:                       var wmsm = new OpenLayers.Layer.WMS("wms", mapurl,           {               layers: Mylayers,               format: 'image/png',               "SLD_BODY": sld },               {"reproject":"true"}                                 );


title: How to set QGIS options in standalone pyqgis application

text: I have a standalone application which uses the QGS api to view and query map layers. One layer is a WMS-C layer which automatically caches tiles which is great. I want to also activate a setting which is present in QGIS but I can't work out how to implement it in my app, namely:   * Use render caching... (Settings -> Options -> Rendering) I found the QSettings key for this in QGIS by browsing the output of:               from PyQt4.QtCore import QSettings     s = QSettings().allKeys()     for k in s:         print k      The one I need to access/emulate/implement is:               Qgis/enable_render_caching      Where can I activate/implement these settings in my standalone app? Is this the setting that is read by `QgsMapCanvas.setCachingEnabled(bool)`?


title: Saving bookmarks in a database using javascript

text: I want to save a bookmark permanently in a database using javascript. I have an idea. "I have a service. so that i shall use that service to interact my code and the database." I am struck with this how to begin. I would like you to suggest whether this is possible. If so how to proceed with this?


title: Implementing Wrapper of PostgreSQL?

text: How can I implement the Wrapper of PostgreSQL? I found this useful looking wrapper from a FOSS4G routing with pgRouting Workshop but I do not know how to implement it!


title: What Z-factor to set in Curvature (spatial analyst)?

text: I'm making a slope curvature map, my DTM is 5*5 meters cells, in a meter grid projection. The elevation is in meters asl. When using the Curvature tool for different units of measure, the Z-factor has to be set. > If the x,y units and z units are in different units of measure, the z-factor > must be set to the appropriate factor, or the results will be incorrect. For > example, if your z units are feet and your x,y units are meters, you would > use a z-factor of 0.3048 to convert your z units from feet to meters (1 foot > = 0.3048 meter). Yet I'm not exactly sure if I need to set a Z-factor of 5, or use the default of 1, because in the "how curvature works" page, the illustration shows that the value for "L", is one raster cell (in my case, 5 meters) What is your opinion?


title: CartoDB - problem trying to add a geoserver WMS as a basemap and it says it doesn't contain supported projections (3857,900913)

text: I had suspected that the version of geoserver (2.0.2) did not support these 2 EPSG's . However, when I then look at the getcapabilities - I find that both EPSG's are present. http://er-info.dpi.vic.gov.au/services/wms?SERVICE=WMS&REQUEST=GetCapabilities any ideas?


title: Assign spatial reference for 1000 different ECW raster datasets using 1000 different *.prj files through python programming or .net

text: I am having 2000 ECW raster datasets without spatial reference and also I have 2000 *.prj files in another folder. I know only the way as right click on raster datasets-->properties-->Spatial References-->Import projection files to assign. It is too difficult to assign each and every file and I am unable to make ensure myself as I assigned correct projection files correspondingly or not. I came to know that if I am following my methodology cannot complete task at a given time. So help me through programming. Python or .dll add in tool. It is very urgent. Thanks in advance who are all participating in this topic.


title: Hitting Memory Limits of PostGIS ST_ConcaveHull / Incomplete Shrink Wrap

text: I am trying to use PostGIS 2.0 on Postgres 8.4 to do some mapping. I am working with the NOAA Shapefiles, specifically their US States & Territories file located here: http://www.nws.noaa.gov/geodata/catalog/national/html/us_state.htm The app I'm building needs to do polygon in polygon tests, and point in polygon tests (PiP). The shapefiles as-is are too detailed, and the queries are taking too long. I'm trying to get PostGIS to simplify the complex state shapes, I want to achieve a high quality shrink wrap of each state, so that rivers and bays are included in the states but small coastal islands or other tiny landmasses are ignored . My goal is to have two tables: both the states raw geometries table for high quality output situations, and also a separate states_simple table where I can use pre-compressed geometries to speed PiP tests. I have finally figured out how to properly implement ST_Concavehull, but I'm running into limits of either the server hardware, or the configuration of postgres. Here is an image to illustrate the issue I'm running into: ![ConcaveHull vs Raw Geometry](http://i.stack.imgur.com/8XWxh.png) Processing "square" states like Oklahoma using ConcaveHull process quickly and without any issue. States that have SOME coastal islands / complexity like Georgia process a little slower, but reasonably fast and come out in a way that will be useful for PiP test (i.e. tight shrink wrap): ![ConcaveHull Georgia](http://i.imgur.com/f9CnQKU.png) However, some states like New Jersey are exhibiting some of the same problems Florida has using ConcaveHull: ![ConcaveHull New Jersey](http://i.imgur.com/eM0KvzI.png) The top and bottom of New Jersey are EXACTLY what I want, but the large part of PA including Philladelphia, and the inclusion of Staten Island basically make this geometry useless for my purposes. So right now the query I'm using to retrieve these polygons is as follows (Example is FL):               SELECT ST_AsGeoJSON(ST_ConcaveHull(ST_Union(geom),.4)) as geom  FROM states WHERE id = 11;      **tl;dr:** If I turn down the tolerance from .4 to something like .6, the state comes out like a giant triangle. .4 seems to be the floor limit of ConcaveHull, anything lower than that and queries run for 3+ minutes without returning results. I've been following this article for tips and steps to increase performance of Postgres and PostGIS: http://workshops.opengeo.org/postgis-intro/tuning.html Changes I've made to postgresql.conf:   * shared_buffers = 256MB    * work_mem = 16MB    * wal_buffers = 1MB   * checkpoint_segments = 6    * random_page_cost = 2.0 The server all this is running on is a Turnkey Linux LAPP Server running on an Amazon EC2 Medium M1 Instance (64 bit, 1cpu, 3.75gb RAM) Am I going about this all wrong? Is there another freely available set of "smoother" CONUS states I can use? I'm theorizing that the relative size (read: huge) of the state geometries is what's gumming up PostGIS/Postgres... just not really sure how to diagnose or monitor that beyond running               ps aux | grep ^ postgres      and watching how long queries take to run. **My main questions:**   1. How can I monitor the resource usage of this process better?   2. Am I using the right PostGIS functions to accomplish this goal?   3. Is it a matter of turning down the tolerance of ConcaveHull, and increasing system limits? Please help! Thank you in advance for reading this far!


title: Displaying JS API Find Task Results in a Table without Dojo

text: I was wondering if anyone had examples of displaying ArcGIS Javascript API Find Task results from a map service in an html table without using a Dojo Grid. Any pointers would be greatly appreciated!


title: polygons won't buffer in ArcMap

text: I have been trying to buffer 15 polygons (-15m) in ArcMap. I did a batch buffer and 13 polygons worked fine, the other 2 have not buffered, an error occurs: WARNING 000117: Warning empty output generated. Executed (Buffer) successfully. End Time: Fri Jul 06 18:58:32 2012 (Elapsed Time: 1.00 seconds) All the input polygons are similar with full attribute tables so the output should not be empty for these 2! Thanks ally


title: FME - MapInfo to File Geodatabase Conversion Project

text: I'm about to start a big project converting 25,000 MapInfo files to File Geodatabase. I'm using FME to process the files. I want FME to search all files in a directory (and subdirectories). When FME encounters a TAB file which references a raster, I get the following error: > TAB files referencing raster data cannot be read through the MapInfo reader. > Instead, use an appropriate raster reader. Error fetching dynamic schema > definitions from reader `MAPINFO_1'. Is there some sort of format attribute for MapInfo files which will identify TAB files which contain geometry and which ones reference Rasters? I would also like to identify TAB files which reference Excel files or Databases. I could try skipping these by ignoring failed readers, ![enter image description here](http://i.stack.imgur.com/mjWbh.png) but I want to be able to list the TAB files that I'm skipping and store those in a spreadsheet or text file.


title: How is distance between zip codes calculated?

text: Specifically, where can I find a database of zip codes, their lat/lon centers, and intervening "obstacles" (for example, zip codes that are geographically close but separated by the Chesapeake Bay).


title: Coordinate points don't show on correct location in the open layers plugin Qgis

text: I know there are several questions concerning the open layers plugin and CRS's. I tried several things, but none of it solved my problem so far. For a rural area in Nepal a little healthcentre wants to know where their patients live. I am trying to map it for them. As a basemap I am currently using Open layers plugin Google Satellite view. Just to test if it is possible to show the coordinates of the houses in the basemap I tried to plot two points of known places with known coordinates (the temple and the school). Somehow when I import the points using Delimited Text layer they show a few kilometres southwest of the area. My project CRS is set to WGS84/Pseudo Mercator EPGS:3857 The CRS of the coordinates vector layer is set to WGS84 EPGS:4326 I hope someone can help me!


title: Arcmap: Create a projection to emphasize a given direction?

text: Suppose one had a dataset (centered on a given coordinate) of lines that were very close together in terms of longitude but also had important turns/bends far separated by latitude. In order to map this, I want to over-emphasize the x direction and de-emphasize the y-direction. Can one do this by creating a projection in ArcMap? If so, how?


title: What are the existing solutions to automate generalisation/simplification of geographical data?

text: I am looking for existing solutions (softwares, libraries, methods, etc.) to automate the generalisation/simplification of geographical data as presented there, and also there. How do you proceed when you need to simplify your geographical data? Do you know and use some tools related to this topic? What is your opinion on these tools? I already had a look at these ones:   * maptimize to simplify the representation of KML marker layers,   * Modelbuilder by ESRI   * Radius clarity by laserscan   * mapshaper to compute online filtering algorithms on geometries,   * Axpand by axes-systems The ICA commission on generalisation and multiple representation website is a huge information source on this topic. Many papers presents advanced methods and use cases related to this topic. EDIT: the **very** famous douglass & peucker filtering algorithm is often used for geometrical simplification and is present in almost all GIS softwares. However, its effect is to make geometries very bony/angular without really simplifying the geometry shape. For example, it is not very adapted for lines representing roads or rivers, that have to be smooth and not angular.


title: What is the best way to store Flickr's Shapefiles Public Dataset in SQL Server 2008?

text: I'm wanting to import specific locations from the Flickr Shapefile Dataset into SQL Server would like to know what is the best way to store the data in SQL Server 2008? Link the the 1.0 Dataset: http://code.flickr.com/blog/2009/05/21/flickr- shapefiles-public-dataset-10/


title: Gdal: How to reproject from equirectangular to mercator?

text: I work with raster ETOPO1, made a crop of my area of interest.               curl -L -C - 'http://www.ngdc.noaa.gov/mgg/global/relief/ETOPO1/data/ice_surface/grid_registered/georeferenced_tiff/ETOPO1_Ice_g_geotiff.zip' -o ./ETOPO1.zip     unzip -n ./ETOPO1.zip '*.tif'     # --- crop     gdal_translate -projwin -5 51 10 41 ETOPO1_Ice_g_geotiff.tif cropXL.tmp.tif     # --- resize for speed sake     gdalwarp -of GTiff -s_srs epsg:4326 -t_srs epsg:4326 -te -5 41 10 51 -ts 1200 0 cropXL.tmp.tif resized.tmp.tif     # --- reproject into mercator     gdalwarp -of GTiff -s_srs EPSG:4326 -t_srs EPSG:3857 resized.tmp.tif resized-mercator.tif      But this result into an incorrect file. **What does I do wrong ? How to reproject correctly from EPSG:4326 to EPSG:3857 ?** **Note:** command inspired by Re-project raster image from Mercator to Equirectangular * * * **Edit:**               yug@yug-PC:~/projects/WikiAtlas_scripts/back$ gdalinfo ./01_reliefs/crop_resized.tmp.tif -hist     Driver: GTiff/GeoTIFF     Files: ./01_reliefs/crop_resized.tmp.tif            ./01_reliefs/crop_resized.tmp.tif.aux.xml     Size is 1962, 1551     Coordinate System is:     PROJCS["WGS 84 / Pseudo-Mercator",         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"]],         PROJECTION["Mercator_1SP"],         PARAMETER["central_meridian",0],         PARAMETER["scale_factor",1],         PARAMETER["false_easting",0],         PARAMETER["false_northing",0],         UNIT["metre",1,             AUTHORITY["EPSG","9001"]],         EXTENSION["PROJ4","+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=@null +wktext  +no_defs"],         AUTHORITY["EPSG","3857"]]     Origin = (11365720.009993232786655,1654650.287315490888432)     Pixel Size = (323.425157454019200,-323.425157454019200)     Metadata:       AREA_OR_POINT=Area       NC_GLOBAL#Conventions=COARDS/CF-1.0       NC_GLOBAL#GMT_version=4.4.0       NC_GLOBAL#history=grdreformat ETOPO1_Ice_g_gdal.grd ETOPO1_Ice_g_gmt4.grd=ni       NC_GLOBAL#node_offset=0       NC_GLOBAL#title=ETOPO1_Ice_g_gmt4.grd       x#actual_range=-180, 180       x#long_name=Longitude       x#units=degrees       y#actual_range=-90, 90       y#long_name=Latitude       y#units=degrees       z#_FillValue=-2147483648       z#actual_range=-10898, 8271       z#long_name=z     Image Structure Metadata:       INTERLEAVE=BAND     Corner Coordinates:     Upper Left  (11365720.010, 1654650.287) (102d 6' 0.00"E, 14d42' 0.00"N)     Lower Left  (11365720.010, 1153017.868) (102d 6' 0.00"E, 10d18' 6.40"N)     Upper Right (12000280.169, 1654650.287) (107d48' 1.26"E, 14d42' 0.00"N)     Lower Right (12000280.169, 1153017.868) (107d48' 1.26"E, 10d18' 6.40"N)     Center      (11683000.089, 1403834.078) (104d57' 0.63"E, 12d30'36.90"N)     Band 1 Block=1962x2 Type=Int16, ColorInterp=Gray       Min=0.000 Max=0.000        Minimum=0.000, Maximum=0.000, Mean=0.000, StdDev=0.000       1 buckets from 0 to 0:       0        NoData Value=-2147483648       Metadata:         NETCDF_VARNAME=z         STATISTICS_MAXIMUM=0         STATISTICS_MEAN=0         STATISTICS_MINIMUM=0         STATISTICS_STDDEV=0


title: Check multiple feature classes have same spatial reference

text: As the title says I've several feature classes. I want to check whether they have same spatial reference using `arcpy`


title: Simplifying geometries (generalization)

text: What are the recommended ways of simplifying geometries? Keeping in mind projections, and simplifying state geometries for example. I've heard about converting to an 'equi-distant' projection that allows simplification without distortion, and then converting back to your chosen projection.


title: Sandbars feature extraction from satellite image

text: I want to know if we can analyse the permanancy of sandbars from a satellite image? Whether the formation of the sandbar is permanent or temporary.


title: How to use IObjectClassEvents?

text: I need to use IObjectClassEvents in my project. I don't know how to start, can any one give me a simple explanation, or an example in C#.


title: Where to get street data in an importable format?

text: I don't want to use WMS layers, but to actually store the data in my local server, where can I find street data from the world?


title: QgsLegend link error when building

text: I'm developping a software using QGIS. When building my software there is a LINK error.               1>mainwindow.obj : error LNK2019: symbole externe non résolu "public: __thiscall         QgsLegend::QgsLegend(class QgsMapCanvas *,class QWidget *,char const         *)"(??0QgsLegend@@QAE@PAVQgsMapCanvas@@PAVQWidget@@PBD@Z) référencé dans la         fonction "public: __thiscall MainWindow::MainWindow(class QWidget *,class         QFlags<enum Qt::WindowType>)"      (??0MainWindow@@QAE@PAVQWidget@@V?$QFlags@W4WindowType@Qt@@@@@Z)     1>C:\Documents and Settings\XXXXX\mes documents\visual studio 2010\Projects\XXXXX     \XXXXXX.exe : fatal error LNK1120: 1 externes non résolus      I'm using QGIS 1.7.2, QT 4.7.3, Visual Studio 2010 (with VC10), and there were no errors when I built the sources of QGIS. I've already put all the .libs i found in the link editor of VS10. In QGIS API, it is written that for the QgsLegend, I had to link the qgis_gui.lib, so had I. I don't know what I can do now. Can someone help me ? Or how to proceed in order to find the lib i missed. If you need more informations, just ask. I'm really blocked... Thank you in advance for the help !! Regards Jonathan


title: OpenLayers.Bounds transform not working as expected

text: I have problem with transforming bounds. It is not working for the first time, later it works just fine. At first i thought that something was still loading butsearching in firebug shows that nothing gets loaded and i think everything is provided for it to work: Here is the code:                       Proj4js.defs["EPSG:2180"] = "+proj=tmerc +lat_0=0 +lon_0=19 +k=0.9993 +x_0=500000 +y_0=-5300000 +ellps=GRS80 +units=m +no_defs";                  var proj2180 = new OpenLayers.Projection("EPSG:2180");             var proj4326 = new OpenLayers.Projection("EPSG:4326");                  var bounds = OpenLayers.Bounds.fromArray(record.get("llbbox"));             bounds.transform(proj4326, proj2180);             console.log(bounds);      This code is placed in function. When i run it for the first time i get following console output: 13.897781645001,48.960416159343,24.621296204369,54.83756633314 when second time i get: 143457.09623933677,137520.1675486639,876176.7241084952,786832.4073041743


title: How to get routes in the database that intersects circle?

text: I have routes (linestrings) recorded in the database. I want to create a query which returns the routes that have some commonality with a circle whose center and radius are passed as parameter. I tried this way:               SELECT *      FROM   routes     WHERE  ST_Touches(geom, (ST_Buffer(ST_GeomFromText('POINT(-12.988542214345332 -38.44563066959381)'), 100, 8)))       Is it correct?


title: ArcGIS Javascript API Editing - Create Feature from Geometry Object

text: I'm interested in passing a geometry object with certain attributes to an exposed feature class for committed edits. Essentially, I want to manually create a new feature without the template pane or editor widget. The geometry and attributes are both pre-defined, and the user just needs to select one pair from a group and push an "Add to Geodatabase" button. I'm reading through the documentation and I can't see a clear way to do this. I see that if I don't want the pane or widget, then I need to use a geometry service, but I don't see a "create feature" method, or something that accepts a geometry object or attributes.


title: Adding WMS and WFS layers from GEOSERVER in arcgis viewer for flex 3.0

text: I am interested in adding feature layers from geoserver in arcgis viewer for flex config file and consequently display them on the web map.IS that possible? If yes, could someone share how they did that and perform querying on the features (points and polygons)? Thank you very much. Simon.


title: Merging adjacent polygons in shapefile that has been split at tile boundaries?

text: I'm currently using building data from the Ordnance Survey Vectormap District and I've noticed that a lot of building polygons are split in 2 because of the way the OS split up the shapefile. This is what I'm seeing in QGIS – I've highlighted some of the buildings split at the tile boundary: ![enter image description here](http://i.stack.imgur.com/9FJZM.png) I can't manually merge the buildings as there are way too many of them that have been split at the boundaries (this is just a zoomed in example). I'm looking for an automated way to solve this. Unfortunately, each side of the split buildings has a different id so I can't easily dissolve them. How would you suggest I can automatically merge the building tiles? Thanks, Rob **Edit** I'm now using PostGIS to merge the polygons split at the tile boundaries. Here is the SQL statement that does this for me – it's many, many times faster than doing it with QGIS:               DROP TABLE merged;     CREATE TABLE merged AS     SELECT       -- Merge polygons that are within 0.01 metres of each other       (ST_Dump(ST_Union(ST_Buffer(the_geom, 0.01)))).geom AS the_geom     FROM unmerged;          -- Update the geometry_columns table     SELECT Populate_Geometry_Columns();


title: Delaunay triangulation doesn't work

text: I am trying to create a Delaunay triangulation of a set of 3692 points with QGIS 1.8.0. The resulting set of polygons appears unfinished, leaving a small part of the original map uncovered by triangles and a small set of points not used in any triangle. Is there a limit of the number of input points or is there some other constraint? Best regards from France


title: Why is not ST_GeomFromEWKB() working?

text: The following SQL is working fine from the terminal:               SELECT MIN(ST_Distance(       ST_Transform(ST_GeomFromText('POINT(-3.163916 55.926528)',4326),32630),       ST_Transform(ST_Boundary(ST_GeomFromEWKB(geometry)),32630)     )) FROM geo_table      BUT, when I try it through a python psycopg2 try-except it fails and postGIS returns a `WARNING: OGC WKB expected, EWKB provided - use GeometryFromEWKB() for this`. I am slightly confused here, any help would be much appreciated! I have more queries in my code that they use a stored geometry without the need of a wrapper like ST_GeomFromWKB. I don't understand why this happens. Cheers, A


title: OpenLayers does not render wms layer generated by MapServer

text:
I need to add a wms layer to my OpenLayers map and it should be generated by MapServer installed locally. **OpenLayers HTML file**               <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN" "http://www.w3.org/TR/html4/loose.dtd">               <meta http-equiv="Content-Type" content="text/html; charset=utf-8">     <title>OpenLayers Demo</title>

                **Map File**               # To enter comments into a mapfile add a # symbol at the beginning of the line               MAP     NAME parcels_map     STATUS ON     SIZE 420 300     EXTENT 489885 5450946 490904 5451637     UNITS METERS     SHAPEPATH "shapefiles"     IMAGECOLOR 255 255 255          WEB     TEMPLATE "template.html"     IMAGEPATH "C:\\xampp\\htdocs\\output\\"     IMAGEURL "/output/"     METADATA         "wms_title" "WMS Demo Server"         "wms_onlineresource" "http://localhost:8080/cgi-bin/mapserv.exe?map=C:\\xampp\\htdocs\     \mapfile.map"    "wms_srs" "EPSG:4326" "wms_feature_info_mime_type"  "text/html"     END      END          LAYER     NAME parcel_boundaries     DATA parcels     STATUS DEFAULT     TYPE POLYGON     CLASS     NAME "parcels_color"     COLOR 255 215 0     END     END          # end of map file     END      Obviously MapServer can generate the map using the given shape file. I have tested it by entering following URL to the browser.               http://localhost/cgi-bin/mapserv.exe?map=C:\\xampp\\htdocs\\mapfile.map But, when adding this URL as wms service URL, it just gives broken image tiles. Any help would be very much appreciated.


title: MapServer WMS feed and QGIS

text: I've managed to create a WMS feed from Mapserver showing a tileindex to some rasters and a shapefile and it appears in my browser OK. If I go to QGIS (1.8.0 Windows StandAlone) and I connect to my WMS feed I can see the layer names appear in the dialog box. However when I add the layers to the map I don't see any rasters or shp layers. I tried zooming to them etc but nothing is visible on the screen. Does anyone have any ideas why this might be?


title: GIS/Developer Skills

text: Currently, I have a strong knowledge ArcGIS/Python/JavaScript and basic RDBMS. I'm currently taking some database theory/sql classes and trying to figure out what course I should take next, I'm considering Java. I've been told ArcObjects and .NET is a thing of the past and to focus on Java, is this true? I see most GIS/Developer job postings requesting proficiency in ArcObjects .NET C#.


title: Calculating distance between borders when one polygon inside another polygon?

text: I am brand new to ArcGIS and having some trouble. I would like to calculate the distances between the border of one polygon to another polygon. However, these polygons overlap. Thus when I use the Near function, it results in zero because the polygons overlap at some points or one is within a larger one. I cannot use the Point Distance because I have shape files and so I have nothing to choose from in the drop down menu. Thank you for any help!!


title: How to install QGIS 2.0 on Mint 15?

text: Trying to use QGIS 2.0 on Mint 15 (based on Ubuntu 13) I installed, uninstalled and re-installed QGIS from various repositories (ubuntugis-unstable, then qgis.org/debian raring - not clear for me which one I should use). I don't seem to have any python support. Any help appreciated! If I try               apt-get install python-qgis      I get the message that               the packages libgdal1 libqgis2.0.1 qgis qgis-providers      are going to be uninstalled          Bump: Still no solution. I uninstalled QGIS and `apt-get autoremove`, kept only ppa:ubuntugis/ubuntugis-unstable. Now when I try               sudo apt-get install qgis python-qgis qgis-plugin-grass      the error message is               Les paquets suivants contiennent des dépendances non satisfaites :      python-qgis : Dépend: python-qgis-common (= 2.0.1-2+raring1) mais ne sera pas installé      qgis-plugin-grass : Dépend: grass642     E: Impossible de corriger les problèmes, des paquets défectueux sont en mode « garder en l'état ».      Can someone help?


title: Whats is Openlayers editor?

text: Staff has anyone ever used OpenLayers Editor? Openlayers researching on the internet I found this link: https://github.com/geops/ole It is an editor based on OpenLayers open source? How do I set it on my machine to test with PostGIS and GeoServer?


title: Options to Replace Interactive Mapping System

text: I have been tasked with replacing our organizations online interactive mapping system. It is your run of the mill County Government interactive map. Its bulky, unreliable, slow and out of date. It is currently running on `ArcIMS` and hasn't been supported in ages from ESRI. 44 layers are currently loaded onto the map. In general we are moving towards content specific applications and I have created several that divide those 44 layers into themed apps (planning & zoning, public safety, etc) I have done this using a combination of ArcGIS Online and ArcGIS JavaScript API. Those products work well for the smaller targeted applications, however, I can't seem to come up with a good solution for building the all-in-one interactive map replacement (it could be argued its not needed, but its a requirement). Is anybody aware of a good solution or template (preferably free) for this large scale application? I have recreated all of the map layers in services using ArcGIS for Server 10.1, I have marginally decent ArcGIS JavaScript API skills and an ArcGIS Online Organizational account.


title: Avoid text/shields from being too close using Mapnik xml

text: I'm having the problem of text/shields from different layers appear too close to each other on the map generated from Mapnik xml. The screenshot below shows the problem ![enter image description here](http://i.stack.imgur.com/Jeevr.png) I have three layers in this map.   1. A layer for shields, appeared as "rectangle 972" in the image   2. A layer for major roads, which has "SW 24th St" appear on the map   3. A layer for secondary roads The problem is the "rectangle 972" and "SW 24th St" are too close to each other, and the street names for the secondary roads are also concatenated with each other (I'm less concerned with this though). I couldn't find the right configuration option in Mapnik for this to look better. I have tried setting the minimum-distance in ShieldSymbolizer and TextSymbolizer, but I don't think that's the right way to do, and it didn't work of course.


title: Create isochrone from multiple start points

text: I want to create an isochrone. For this I could use the pgr_drivingDistance function with pgRouting (http://docs.pgrouting.org/dev/src/driving_distance/doc/index.html) but it only takes one start node. How could I compute isochrones efficiently from multiple start nodes?


title: Unable to create shapefile in QGIS

text: I am very new to QGIS and although I've been able to download the software and get aerial imagery using the plugins menu, I'm getting nowhere with creating a shapefile that I can use to draw a polygon on the aerial image. I've gone through the layer - create new shapefile procedure many times, but when prompted to save the shapefile on my C:drive, nothing happens. I've tried using "add vector layer" then searching for the file I thought I'd saved, but there is nothing there. Nor does the shapefile appear on the layers menu on the left of the project. Where am I going wrong please?


title: How to create a QGIS Installer for Windows with GUI customization?

text: I'm trying to build a custom QGIS installer for Windows with a stripped down GUI. I customized the GUI using the customization function of QGis and saved it as a .ini file. The question now is: How can I include this custom GUI in my installer? I found this very useful tutorial: http://linfiniti.com/2011/05/building-custom-qgis-installers-for-windows/ and in the comments it also says that it's possible to include it using the customization tool. However, it doesn't exactly say how and where to save the customization file in the installer package. Another possibility automatically load a custom GUI is explained here: http://linfiniti.com/2011/07/batch-customizing-qgis/ Thought about exporting the registry key customization.reg and to include it somewhere in the installer, but where? So that would mean I would have to export the Windows registry key customization.reg, save it somewhere in the installation package and write a batch file that imports this registry key after the installation is finish. Do you think that would be possible? So that's what I have found out so far. Does anybody have an idea what the easiest way would be to do this? Thanks! :)


title: Covisibility with QGIS/GRASS : r.los issue

text: I try to use r.los for a viewshed calculation. I specify the dem (Aster DEM clipped on a 2km*2km region). I specify x and y coordinates and run the command. Nothing append. I waited for 20 minutes without a result. The loading cursor stay at 0. No error message, no stop, nothing. I'm using r.los in QGIS/GRASS with an 1.8 version and under Windows 7. Do you kwno why ? Thanks


title: D3 for maps---at what stage to bring in data to the geo?

text: I'd like to map a world choropleth for display with D3, a la:   * http://bl.ocks.org/mbostock/4206573   * http://bl.ocks.org/jasondavies/4188334 I have a dataset I'd like to display that's keyed to ISO-alpha-3 keys. So...               danger.csv     iso,level     AFG,100     ALB,0     DZA,12      etc. Following the instructions on topojson, I know I can do...               wget "http://www.naturalearthdata.com/http//www.naturalearthdata.com/download/50m/cultural/ne_50m_admin_0_countries.zip"     unzip ne_50m_admin_0_countries.zip     ogr2ogr -f "GeoJSON" output_features.json ne_50m_admin_0_countries.shp -select iso_a3     topojson -o topo.json output_features.json --id-property iso_a3      to produce a worldmap json that is ID'd by ISO3. My question is: at what point in the workflow should I merge the data from danger.csv onto the geo data? I had previously worked with qGIS as a GUI, but where /should/ the merge happen? In the .shp? After the ogr2ogr? Dynamically in the browser after the topojson shrink (like here http://bl.ocks.org/mbostock/4060606 http://bl.ocks.org/mbostock/3306362)? I'm pretty good with python, but pretty new to javascript, and find myself copying and pasting Bostock examples more than actually being a generative coder there. (I also have a related, but more involved follow-up on Stackoverflow that maybe I should migrate here: http://stackoverflow.com/questions/18604877/how- to-do-time-data-in-d3-maps)


title: Convert PCIDSK file with seven layers to Geotiff without loss of layers

text: does anybody know a way to save a PCI Geomatica landsat file (PCIDSK, *.pix) as a Geotiff without loss of layers?? The Transferring layer function of Orthoengine does not work. With the export map function I can only save a three layer Geotiff. With Arcmap it is the same thing. What about R and QGis? Greetings, Stacky


title: Implementing Custom Selection Tool (Adding to current selection) with ArcObjects

text: I am trying to implement a custom selection tool in ArcMap 9.3.1 that is functionally similar to the default ArcMap selection tool with "Adding to current selection" enabled. ie. the user can click to incrementally select individual features one by one, while adding on to previously selected features. When the user clicks on an already selected feature, that feature should be unselected. How should I go about replicating and implementing this behavior for my custom tool? There are some selection-related code snippets lying around, but nothing really comes close. I was able to perform single selections by calling SelectByShape(geom, null, true), where geom = a small area around the clicked point. To support multiple selections at various locations, how then do I keep the previously selected features still selected as I select a new one, since SelectByShape(geom, null, true) only selects a single feature each time it's called. Clearly, I will need another approach, which probably involves me keeping a growing list of currently selected features, but I am not sure if there is a more straightforward solution. Has anyone attempted this before? Many thanks!


title: Is it possible to get the points used for the convex hull in osm2po?

text: I have a non-convex hull algorithm that I want to use for an app. With OSM2PO, it's possible to get a convexHull, but is there a way to get the points before the hull is constructed? More, it would be awesome for me if I can get only the points that are "path end", I mean points that are part of segments that are connected to no more than one other segment. (If you're intererrested by the non-convex hull algorithm, aka Alpha-shapes, I can share...) * * * Thanks for your quick answer, but I don't understand well your solution. I've seen this table, and it might be helpfull but I don't see how to query PGSQL (pass it millions indices to check ?) and my aim is to avoid PGSQL. I'dont know how you dit it but when you traverse the graph, there's probably something that say "Hey stop traversing this branch, you've reached the max cost or this branch has no more points". It should be possible to grab these "last" points, don't you think ? In order to be sure we understand well, here's a little example. Imagine this "graph" :                        7                         /                       6                      /   1---2---3---5---8----9---10              |              4  `router.getVisited()` will give me `1,2,3,4,5,6,7,8,9,10`. My dream : only get `1,4,7,10` Hope you can help me. Many thanks again. (sorry for answering my own question but too many chars to be posted as a comment...)


title: Difference between gdal versions in reading raster on two machines?

text: I am using raster data from Afripop.org and using the R package raster (raster_2.0-41) to do all my computation. Everything has been working out just fine until I tried sending some jobs to a remote linux machine. When I read in the raster on one machine the values are correct but on the other there are no NAs and the population (i.e. sum of all the values) WAY off. On the linux machine I get:               > civpop <- raster("data/apciv10v3.flt")          > civpop[1:10]     [1]   0  60  28 198   0  60  28 198   0  60          > civpop     class       : RasterLayer      dimensions  : 7656, 7327, 56095512  (nrow, ncol, ncell)     resolution  : 0.0008333, 0.0008333  (x, y)     extent      : -8.599802, -2.494213, 4.357029, 10.73677  (xmin, xmax, ymin, ymax)     coord. ref. : +proj=longlat +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +no_defs      data source : /home/andrewa/D4D/trunk/data/apciv10v3.flt      names       : apciv10v3      values      : 0, 255  (min, max)      One the other machine I (correctly) get:               R> civpop <- raster("data/apciv10v3.flt")     R> civpop[1:10]      [1] NA NA NA NA NA NA NA NA NA NA     R> civpop     class       : RasterLayer      dimensions  : 7656, 7327, 56095512  (nrow, ncol, ncell)     resolution  : 0.0008333, 0.0008333  (x, y)     extent      : -8.6, -2.494, 4.357, 10.74  (xmin, xmax, ymin, ymax)     coord. ref. : +proj=longlat +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +no_defs      data source : /Users/andrewa/Documents/D4D/trunk/data/apciv10v3.flt      names       : apciv10v3       If I use just rgdal I have the same issues so I guess it isn't an issue with raster. Any ideas? Incorrect Machine (using GDAL 1.8.0, released 2011/01/12):               > civpop.gdal = readGDAL("data/apciv10v3.flt")     data/apciv10v3.flt has GDAL driver EHdr      and has 7656 rows and 7327 columns     > civpop.gdal@data[1:10,1]     [1]   0  60  28 198   0  60  28 198   0  60      Correct machine (using GDAL 1.8.1, released 2011/07/09):               R> civpop.gdal = readGDAL("data/apciv10v3.flt")     data/apciv10v3.flt has GDAL driver EHdr      and has 7656 rows and 7327 columns     R> civpop.gdal@data[1:10,1]     [1] NA NA NA NA NA NA NA NA NA NA`


title: Erasing overlay geometry with in the feature class..!

text: What is the best practice for erasing overlay geometry? My scenario is: I have a polygon layer with 1000 features. i need to create a tool from model builder that select polygons with user defined attributes, and verify if there is any other polygons are overlapping. Then erase the overlay area from overlay polygons. EX: Polygon A is selected by attributes and polygon B is overlapping with A. Erase overlapping area from B. I am able to find the overlapping polygons, but when i use Erase tool it is creating new feature class but i need output should be modified with in the same feature class. This is nothing but **Clip in Edit toolbar** with " ** _Discard the area that intersects_** ".


title: Adding .e00 data into QGIS

text: I wonder if anyone can advise me on how to go about adding .e00 data into my QGIS project. I'd be very grateful for any advice, or a hint on where to find it. Thanks TVZ


title: developer problem with PyQGIS

text: I am a QGIS developer. We have been modifying the source code of QGIS for several months and have done good work in QGIS version 2.1. After that we tried to do the same modifications for QGIS version 2.5. A lot of code had to be modified, but we managed to build it. However, we encountered a problem in python, which is this:               Couldn't load PyQGIS.     Python support will be disabled.     Traceback (most recent call last):       File "", line 1, in        File "D:/qgis/build/x86/qgis2.5.0/./python\qgis\__init__.py", line 36, in          from qgis.core import QgsFeature, QgsGeometry       File "D:/qgis/build/x86/qgis2.5.0/./python\qgis\core\__init__.py", line 1, in          from qgis._core import *     ImportError: dynamic module does not define init function (init_core)     Python version:     2.7.4 (default, Apr  6 2013, 19:54:46) [MSC v.1500 32 bit (Intel)]     QGIS version:     2.5.0-Master 'Master', exported      I have looked up this error and found that it may be related to PYTHONPATH or PATH environment variables. I tried all possible values for both environmental variables without success. I am working on a 64-bit Windows 7, and the QGIS build is 32 bit.


title: How to install GRASS on Linux Mint?

text: I'm trying to figure out an error in GRASS network analysis, for several days now. I started to suspect that it has something to do with my OS, namely Windows 7. Hence I installed Linux-Mint 14 (Nadia) on my computer. I would like to run GRASS 6.4 (or latest) on Linux to test the function, yet can't figure out how to do so. Can someone help me with it? or direct me to some where on the web where it is explained? Please advice, Dor


title: compile QGIS-plugin for Windows: problem with gdal_i.lib version

text: While compiling my plugin for QGIS 2.2 there are no errors, but it ist still not loading in QGIS because of an dll error. (function not found) I think the problem is at the gdal version. The QGIS-About says it is compiled with GDAL 1.10.1 but my plugin is compiled with GDAL 1.11. So I tried to compile my plugin with referencing the gdal_i.lib form the installed QGIS named "C:\Program Files (x86)\QGIS Valmiera\lib\gdal_i.lib" but cmake says while configuring that it is gdal version 1.11. Why? I tried it with the lib fom QGIS Lisboa - same result. Although the libs have different sizes and dates. Where is the problem?


title: Permanent add Geoserver Server on Geoexplorer different URL

text: I have my Geoexplorer run on Open shift (asume http://geoexplorer-aaa- rhcloud.com) using this Git (https://github.com/jason-callaway/geoexplorer-on- openshift). My geoexplorer run perfectly, Problem is my geoserver are on another gear on openshift (asume http://geoserver-aaa-rhcloud.com/web). I have Github File for making Geoexplorer.war (https://github.com/boundlessgeo/suite/tree/r4.1/geoexplorer) and I know where should I put my new "Edited" Geoexplorer.war on 1st Git to push on open shift. My question is how to make my geoexplorer permanently add my geoserver url, at which part of the Git source I should add / edit so that I dont have to repeat step "Add New Server" each time I open / reload my geoexplorer site?


title: How to get Shapefile point coordinates into .xls?

text: How can I get just Coordinates (in a file) from point features (shapes) without big hussle. Something like Shape to .xls, with distinct coordinates? I want to use this coordinates immediately, so GPX is no solution. Thanks in advance for your help.


title: Line start and end nodes: determining what is a junction

text: I have a layer of points and polylines and want to determine which points are road junctions. I currently have all start and end nodes within 3m of my road network as a layer. However I cannot determine which ones are simple joins between sections or what is a junction, e.g: ![enter image description here](http://i.stack.imgur.com/TDBFi.png) In the image above, the highlighted point is not a junction, but it is where one section ends and another begins. Is there a way to distinguish between these and points like the others in the image without manual checking?