title: Is it possible to figure out which "Path"(PointCollection) is clicked on within a PolyLine?

text: I am using ArcGIS WPF's `Polyline` to plot several lines (`PointCollection`s). Adding a Path(`PointCollection`) looks something like this.               PointCollection collection = new PointCollection();     PolyLine myLine = new PolyLine();          collection.Add(firstMapPoint);     collection.Add(secondMapPoint);          myLine.Paths.Add(collection);      `PolyLine` is a Geometry in which I can set to a `Graphic` and plot on a map. Once that `Graphic` is clicked it of course is clicking the entire graphic. I am visually showing several Paths(`PointCollection`s). Is it possible to figure out which Path is being clicked on from the user? **Note** I have a solution now where I just have 1 `PointCollection` per `PolyLine`. So I am adding one `PolyLine`/`Graphic` for each line. I would rather not do this if I can avoid it.


title: Reverse geocoding polygon co-ordinates?

text: I have a large dataset of polygons for property sections that I need to reverse geocode. It was recommended that I use a GIS like Quantam to do this however from my research I am not 100% sure if a GIS will support doing this. The technique that I was going to use until I found GIS was to get the center co-ordinates of the polygon and then pass these co-ordinates into a reverse geo-coding api like google or openstreetmap. However I am not sure if this is very accurate. The polygon data is in the projection format WGS 84 (EPSG:4326 Lat/Long) Here is an example of what one of the polygons data looks like: > MULTIPOLYGON (((174.76394635 -41.34094115,174.763921416699986 > -41.340960566699998,174.763611683300013 > -41.340978616699999,174.763597816700013 -41.3409202,174.763793716700008 > -41.340691133299998,174.7638963 -41.34071535,174.763958516700001 > -41.34073005,174.76394635 -41.34094115))) What is the best way to reverse geocode this polygon data? How can I do it on a large scale?


title: CartoDB Convert from the_geom to google's lat long

text: I have a CartoDB table that when I query it via the rest api I get something like this:               {       "name": "sas",       "street": "la plata 2390",       "cartodb_id": 1,       "created_at": "2013-06-26T00:46:21.412Z",       "updated_at": "2013-06-26T00:58:26.799Z",       "the_geom": "0101000020E6100000FDA4DAA7E32D4DC0E3FC4D28445041C0",       "the_geom_webmercator": "0101000020110F0000064377E72DC858C1BB2873CAD2614FC1",       "cartodb_georef_status": true     },      I want to conver the the_geom_webmercator to a lat, long pair of numbers, to use with leaflet like this:               L.marker([-34.617917, -58.433257]).addTo(map)       .bindPopup('A pretty CSS3 popup. <br> Easily customizable.')       .openPopup();      Which would be the best way to achieve it? Create a view or query on CartoDB? (I tried with ST_Transform, but I don't know what parameter to use) Use a specific leaflet option to pass the webmercator value? Use some js library to do the conversion? -- EDIT I've found the st_asText(the_geom) that returns me something like: POINT(-58.35851 -34.62708) and also st_x(xxx) and st_y(xxx) functions which returns each coordinate. so I create a the_geom_lat and the_geom_long columns and executed the following query               update test set     the_geom_lat = st_y(the_geom),     the_geom_long = st_x(the_geom)      I still would like to know if there's any js function that could achieve the same, so I wouldn't have to create a cartoDB view


title: Open layers with Arcgis Rest api identify layer feature

text: We are using ArcGIS Rest services on Openlayers. The layers are placed on the map correctly But how to identify a layer feature on mouse over or mouse click and get the details of the particular layer feature...? Is there any function in Openlayers library to do it..?


title: Equal Area Projection for Central Ethiopia (37N)

text: Can anybody please tell me what is the best equal area projection for central highlands of Ethiopia. I am assuming Albers Equal Area Projection. Any suggestions?


title: How to identify image footprints of an aerial photograph

text: I need a tool that will take aerial images with Lat, Long, and Alt in their EXIF data and produce foot prints of the images. We need to verify that the images cover the photographed area with proper overlap.


title: Why is ArcGIS Server so slow in creating mapbook PNG

text: We have 3 MXDs for creating 3 different mapbooks. Our customer want a daily export to PNG. For this, I use a Python-Script, which works fine for me on my Desktop. But I need it on ArcGIS Server. There it works also, but is really slow. In fact, it takes up to 2.5 times on Server than on Desktop. Two of the MXDs have to export 173 PNG, the third exports only 32. Times on Desktop: One of the 'big' MXD (with 173 PNG) needs ca. 5 hours ( ~106 seconds per PNG), but on the Server it needs ca. 12 hours ( ~255 seconds per PNG). I tried a SDE as datasource (which is the way i normally want to go) and a fGDB. Also Labelling on and off. But there where no specific difference in processing speed. Server is an ArcGIS Server 10.2.1. Can anyone give tips for a better processing time or have the same issue? And please, know comments like these mapbooks are too big for a daily export ;) Another question, if someone have experience with that: How much could multithreading give here?


title: ArcGIS 10.1 Aggregate Points cannot create output

text: I'm trying to use the tool Aggregate Points with output created using in_memory but the tool keeps failing with: "ERROR 000210 Cannot create output \Output_Tbl" The tool works fine if I use a GDB feature class as output. Does anyone know how to use output for in_memory with this tool? ![enter image description here](http://i.stack.imgur.com/1kktK.png)


title: v.buffer angular buffer doesn't take all the polygons into account

text: I am trying to create an angular buffer around building footprints with the v.buffer command (QGIS 2.2). There is no error message, but the algorithm doesn't take into account all the buildings. This happens only with the flag "make outside corners straight". The "normal" buffer works perfectly for the same dataset. I've tried with v.buffer.distance and v.buffer.column, the result is the same. Do you have any ideas? Can it be that it worked in a previous version of QGIS? Any feedback is highly appreciated!


title: How to show coincident point with different symbology in webgis?

text: our team have develop an webgis using ArcGIS silverlight API 3.1 .we have some points in database.These points are save in a table in sql server 2008R2 not in Geodatabase.The points have same longitudes and latitudes. I want to show these points with different symbols( colors or sizes) . I surf in the internet and found clustering . But the problem is that i have coincident points and clustering is use when the latitude and longitude of point are different. How can i show my points in webgis ? Thanks for advance


title: Is there a Small Windows install of QGIS?

text: I am looking to find a faster alternative to ArcReader. I think QGIS might be faster but I was wondering if there was a smaller install without all the GRASS etc for windows. I just want to get a fast map out to my users in the field. The ability to Zoom, Pan, Identify and perhaps search for an address would be nice.


title: MXD Performance - number of layers vs number of datasets

text: My mxd currently has 5 buildings layers accessing 5 different buildings datasets (one for each of 5 counties). Most often there would only be a single county's buildings in the current view, if that matters. Assuming there is nothing else causing speed issues, **would my mxd likely draw faster (or worse) if I combined the 5 datasets into a single buildings dataset?** And if it does, what would be providing the speed increase? Having less layers to display, or having to access fewer datasets? **EDIT:** I appreciate the input on asking whether the data is local or on a network, what data type it is, etc. However, I am more interested in **whether there is a significant cost associated with the number of layers in the mxd** (such as the map having to iterate through each layer when drawing), **whether there is a significant cost associated with the number of separate datasets accessed** , and whether one cost is typically more costly than the other. I'm pretty sure that the cost of accessing multiple shapefiles is greater than the cost of accessing the same data as file gdb feature classes, but how do those compare to the cost of rendering layers? I was hoping that someone knowledgeable of the inner workings of ArcMap/ArcEngine would know how these two variables effect map drawing performance. *All of the data is stored in file geodatabases on the local machine. Currently the building datasets are in 5 separate file gdbs, but I plan to move them all to the same gdb soon.


title: Checkout Spatial Analyst extention, ArcObjects code

text: **Update** After investigating more closely the CODE 2 sample below (at the end of this post), it does check out spatial analyst (sa) extension, however it does not reflect it in the Tools>Extensions menu option as being checked. However you can still use the sa tools. --------------Original question--------------- I would like to checkout spatial analyst extension for ArcMap 9.3.1 via this ArcObjects code below, but when I run the code (CODE1) I get these two message box statements and the extension does not checkout under the Tools>Extensions dialog. > "You license has already been initialized please check you implementation" AND > "SpatialAnalyst license has been checked out for use" I'm not sure what exactly the first message box means? Any suggestions would be greatly appreciated. I also found this code (CODE2) to turn on spatial analyst extension, but it does not work either. Maybe I'm missing something obvious with this one.? Thank you! CODE 1 * * *               Private Sub UIButtonControl6_Click()            Dim licenseStatus As esriLicenseStatus            licenseStatus = CheckOutLicense(esriLicenseProductCodeEngine)            If (licenseStatus = esriLicenseNotLicensed) Then              licenseStatus = CheckOutLicense(esriLicenseProductCodeArcView)            If (licenseStatus = esriLicenseNotLicensed) Then              licenseStatus = CheckOutLicense(esriLicenseProductCodeArcEditor)                If (licenseStatus = esriLicenseNotLicensed) Then             licenseStatus = CheckOutLicense(esriLicenseProductCodeArcInfo)           End If         End If       End If            'Take a look at the licenseStatus to see if it failed       'Not licensed            If (licenseStatus = esriLicenseNotLicensed) Then              MsgBox "You are not licensed to run this product"            'The licenses needed are currently in use           ElseIf (licenseStatus = esriLicenseUnavailable) Then              MsgBox "There are insufient licenses to run"           'The licenses unexpected license failure            ElseIf (licenseStatus = esriLicenseFailure) Then              MsgBox "Unexpected license failure please contact you administrator'"            'Already initialized (Initialization can only occur once)            ElseIf (licenseStatus = esriLicenseAlreadyInitialized) Then              MsgBox "You license has already been initialized please check you implementation"            'Everything was checkedout successfully            ElseIf (licenseStatus = esriLicenseCheckedOut) Then              MsgBox "Licenses checked out successfully"            End If              Dim bAlreadyCheckedOut As Boolean            bAlreadyCheckedOut = m_pAoInitialize.IsExtensionCheckedOut(esriLicenseExtensionCodeSpatialAnalyst)            If (Not bAlreadyCheckedOut) Then              licenseStatus = m_pAoInitialize.CheckOutExtension (esriLicenseExtensionCodeSpatialAnalyst)              If (licenseStatus = esriLicenseUnavailable) Then                MsgBox "All SpatialAnalyst licenses are currently in use"              ElseIf (licenseStatus = esriLicenseCheckedOut) Then                licenseStatus = esriLicenseCheckedOut                MsgBox "SpatialAnalyst license has been checked out for use"              Else           MsgBox "Unexpected licensing failure contact you administor"         End If       End If            End Sub          Private Function CheckOutLicense(productCode As esriLicenseProductCode) As esriLicenseStatus            Dim licenseStatus As esriLicenseStatus       Set m_pAoInitialize = New AoInitialize       CheckOutLicense = esriLicenseUnavailable            'Check the productCode            licenseStatus = m_pAoInitialize.IsProductCodeAvailable(productCode)       If (licenseStatus = esriLicenseAvailable) Then              'Check the extensionCode              licenseStatus = m_pAoInitialize.IsExtensionCodeAvailable(productCode, esriLicenseExtensionCodeSpatialAnalyst)         If (licenseStatus = esriLicenseAvailable) Then           'Initialize the license           licenseStatus = m_pAoInitialize.Initialize(productCode)         End If       End If            CheckOutLicense = licenseStatus          End Function      CODE 2 * * *                 Dim gp As Object                      gp = CreateObject("esriGeoprocessing.GPDispatch.1")                      If gp.CheckExtension("spatial") = "NotLicensed" Then                          MessageBox.Show("NotLicensed")                      ElseIf gp.CheckExtension("spatial") = "Unavailable" Then                          MessageBox.Show("Unavailable")                      ElseIf gp.CheckExtension("spatial") = "NotInitialized" Then                          MessageBox.Show("NotInitialized")                      ElseIf gp.CheckExtension("spatial") = "Available" Then                          gp.CheckOutExtension("spatial")                          MessageBox.Show("Checkout")                      End If


title: Create a report wizard fails in ArcGIS 10 (win7 x64) - but works in ArcGIS 10 (winxp 32bit)

text: The report creation wizard will not launch for me when I try accessing it from both the View menu and from an attribute table. I am able to open a saved report template though. This seems like it is a Windows 7 issue since the wizard does launch under Windows XP. Or it could be a 64bit vs 32bit issue. I see this has been reported elsewhere (look in the comments) I actually ended up doing the task using Mail Merge in Word 2007. It was so easy to format it together from an exported attribute table - I'd definitely do this again, even if I could get the report wizard to work.


title: Looking for a basic PostGIS viewer to use as basis for own application

text: I'm looking for a rather basic PostGIS viewer. It should be possible to pan, zoom and retrieve coordinate information on mouse click. No need to support coordinate transformations. It should introduce as little further dependencies as possible. Preferably written in C++.


title: How to copy a geometry from one feature to another without attributes?

text: With QGIS it's easy to add a new part to an existing feature/geometry but I can't see how to copy/paste this part from another existing feature/geometry. Is it possible? Thanks


title: Are there any repositories to find custom map tiles for an OpenLayers implementation?

text: I can't seem to find an WMS that looks the way I want. I've also checked all of the regular OSM places too and haven't quite found anything.   Can anyone recommend any places to find custom map tiles that I could use in my OpenLayers implementation? Thanks


title: What is the best way to visualize subsurface drillhole and borehole geology data within ArcGIS?

text: I've heard some folks in mining are using Target for ArcGIS extension by Geosoft to   visualize and analyze subsurface drillhole and borehole geology data within ArcGIS. Does anyone have experience with this tool? What are your impressions? What is the learning curve like? Is there any other comparable extension for ArcGIS with similar functionality?


title: empty network dataset creation

text: I'm trying to create a new network dataset via ArcCatalog in a SDE database. I've imported a street featureclass into my dataset, go to create the network dataset with all the default options. After I hit go, it spins for a minute then comes back saying the creation is complete, but that "the new network dataset has been created no features were found within the extent of the network dataset." What causes this and how do I resolve? Thanks.


title: Creating geodatabase from CSV using python

text: I have a file with weather data (daily) for multiple weather stations for multiple years. I am looking for some python script to automate the database creation process and its interpolation on a daily basis. The weather data looks like this:               YEAR    MONTH   DAY     Date        STID    TMAX    TMIN    TAVG     2000    1       1       1/1/2000    WENT    72.63   44.02   58.27     2000    1       1       1/1/2000    DKOM    70.66   40.77   55.76     2000    1       1       1/1/2000    OLLU    73.15   33.47   51.67     2000    1       1       1/1/2000    LKNP    -999    -999    -999     2000    1       1       1/1/2000    OWNB    70.88   32.01   52.2     2000    1       1       1/1/2000    ANTL    71.85   31.13   51.66      There is a separate shapefile that is associate with the location of the weather stations. How do I create database using python script so that it can be used later for spatial interpolation? Thanks in advance.


title: How to determine if a polygon is self-intersecting using either ArcPad 8.0 or ArcMap 9.3.1?

text: I am try to automate a sync replica mxd a nightly basis. The issue I am running into is that in some the of the AXFs (that I get from the users and check-in) have polygons that are self-intersecting and is causing problems with the sync replica. The sync errors out and the polygons need to be fixed manually. Is there any way to detect if a polygon is self-intersecting in ArcPad 8.0, using VBA or some other technique, to stop the creation of the polygon? If not, is there some tool in ArcMap 9.3.1 that can be run again the mxd to detemine the self-intersecting polygons and fix them? Any ideas will be greatly appreciated. Thanks. Tim


title: How to split Polygons by line grid?

text: I want to split a polygon by a intersecting line into different parts, exactly like in this post(Qgis - Is there a way to split a polygon layer with a line layer?). I have a vector grid to cut the polygons, so I used the intersect tool(vector grid with polygons) first, to get the lines from the vector grid intersecting the polygons I want to cut. Afterwards I followed the procedure like in the post mentioned above (union and then polygonize), but it does not provide a feasible solution, because some layers are cut as I want it and others are not. Does anyone has any workaround about it? I read here (Make a polygon from intersecting lines in QGIS, using Polygonize or similar? ) that there have been some problems with polygonize, maybe this is the reason? FYI: I am running Windows 8, 64 bit, QGIS Desktop 2.4, and installed it by the downloaded file. Thanks for your answers!


title: Add Stacked Column Chart in InfoWindow

text: I want to add a a dojo stacked column chart to my InfoWindow (not pop up) that displays when a user clicks on a county in my map. The chart is based on land use information in each county layer from 3 different years. I'm using this ESRI tutorial to customize my javaScript code: http://developers.arcgis.com/javascript/samples/widget_formatInfoWindow/ Here is a link to my sample fiddle: http://jsfiddle.net/jr_js_dev/HzUC7/1/ **Issues:**   1. My county Feature layer shows sometimes and doesn't show other times. Does anyone know why? When I enter in the sample layer from ESRI, the ESRI layer shows(sampleserver1.arcgisonline.com/ArcGIS/rest/services/Demographics/ESRI_Census_USA/MapServer/3"). But when I switch back to my layer, it doesn't show. (I checked, and the REST service for the layer I want works.)   2. When my layer does show, I am able to click on a county, and an InfoWindow shows up with data based on the county layer and content set in the InfoTemplate. However, the stacked chart doesn't show up in the InfoWindow in the "chartNode" div I created for it. Anyone know how to get the chart to show in the InfoWindow? **Note** : I created a test stacked column chart with sample data that displays how I want my chart to look in the InfoWindow as a example.


title: Calculating fields using a list of features

text: I've written some basic code which creates a new .dbf and adds two fields to it: 1 for today's date, another which should contain the names of shapefiles located in another folder. What I'd like to do is then write the names of the shapefiles to the newly created .dbf.               env.workspace = r'shapefile folder'          AT = r'location of new .dbf'     appendField = "ProjName"     dateField = "DateProc"     arcpy.AddField_management(AT, appendField, "TEXT")     arcpy.AddField_management(AT, dateField, "DATE")          shpList = arcpy.ListFeatureClasses()     for shp in shpList:         print shp         appendExpression = '"'+str(shp).replace('.shp','')+ '"'              arcpy.CalculateField_management(AT, "ProjName",appendExpression, "PYTHON")         arcpy.CalculateField_management(AT, dateField, 'time.strftime("%d/%m/%Y)', "PYTHON" )      I believe my issue is a result of having the table and feature classes in two separate workspaces. I am using 10.1. Thanks for any tips which may get this working.


title: Is it possible to symbolise via a related table?

text: I´m using ArcGIS for Desktop with an Advanced license level. I´ve got a feature class with several points in it; each point contains information for location and Name. E.g. Location [Berlin] , Name [Deponiestrasse 1] Each Point describes a disposal site. Now I want to connect the point information with a Table, which contains several information’s for each point. The Table contains information for Name e.g. [Deponiestrasse 1] and several types of disposal. This types can be [concrete] ; [gravel] ; [brick] ; [toxic waste] and so on, up to 10 for one point. E.g. point 1: Name [Deponiestrasse 1] and types are [concrete] and [gravel]. My Idea to connect them, was to use the RELATE function. But now my Problems are starting. I want to draw around each point a simple “ symbol box”, a small box with different colors e.g. concrete-blue, gravel-red, brick-green……. and so on. But via using the SYMBOLOGY function I can´t symbolise by using a related table. So my questions are: Is there any way to do this in this way? Or is there any possibility to write many types like [concrete] ; [gravel] ; [brick] into one attribute field of one point? Maybe you have another idea how to solve my problem. I hope you understand my problem. I used google etc. but couldn´t find a solution.


title: Identifying geographies that match certain climatic profiles through GIS

text: A certain post-haste project has arisen, and I would like to put this question to the collective mind of GIS.SE for ideas and references. The long and short of this research is trying to find roads that lie within geographic areas that fit a given climatic profile. In this case, I need to find areas with relatively cold falls, winters, and springs that get snow. These places may not have a lot of snowfall on an annual basis, but the snow _persists_ for longer periods due to the colder weather patterns. Are there any websites that already collect similar data in a format ready for GIS analysis? If not, what should I look for when gathering data? There may not be a singly perfect answer for this question, but I mostly am seeking direction to refine my search appropriately. Time is short. Thank you!


title: List index out of range error

text: I've been writing a script that clips large datasets based on a selected feature. It works fine when I run it in a test workspace where only related layers are active and running. However, when I run the script in my organisation's main workspace where alot of unrelated layers are active I get the error "Runtime error : list index out of range". This has baffled me as I have the same layers active in this workspace as I have in the test workspace yet I am getting this error, could it be something to do with the dataframe? or something else? Anyway, here is the script (apologies about the length):               import arcpy     from arcpy import env     from arcpy import mapping          env.workspace = "C:\MyArcGIS\Clip_tool\uFMfSW\uFMfSW_Banded_Vector_SE00.gdb"          arcpy.env.overwriteOutput = True          edrn_node = "C:\\MyArcGIS\\EDRN_Shp\\EDRN_NODE.shp"     ID = "EDRN_043066"     query = '"LCC_DRN_ID"' + " = '" + str(ID) + "'"     polygon_extent102 = "C:\\MyArcGIS\\Clip_tool\\poly_extent102"     depth_1in100 = "uFMfSW_SE00_DEPTH_1in100_BV"     depth1in100 = "C:\\MyArcGIS\\Clip_tool\\uFMfSW\\uFMfSW_Banded_Vector_SE00.gdb\\dep1in100"     depth_1in30 = "uFMfSW_SE00_DEPTH_1in30_BV"     depth1in30 = "C:\\MyArcGIS\\Clip_tool\\uFMfSW\\uFMfSW_Banded_Vector_SE00.gdb\\dep1in30"     depth_1in1000 = "uFMfSW_SE00_DEPTH_1in1000_BV"     depth1in1000 = "C:\\MyArcGIS\\Clip_tool\\uFMfSW\\uFMfSW_Banded_Vector_SE00.gdb\\dep1in1000"     extent_1in30 = "uFMfSW_SE00_EXTENT_1in30_BV"     extent1in30 = "C:\\MyArcGIS\\Clip_tool\\uFMfSW\\uFMfSW_Banded_Vector_SE00.gdb\\ext1in30"     extent_1in100 = "uFMfSW_SE00_EXTENT_1in100_BV"     extent1in100 = "C:\\MyArcGIS\\Clip_tool\\uFMfSW\\uFMfSW_Banded_Vector_SE00.gdb\\ext1in100"     extent_1in1000 = "uFMfSW_SE00_EXTENT_1in100_BV"     extent1in1000 = "C:\\MyArcGIS\\Clip_tool\\uFMfSW\\uFMfSW_Banded_Vector_SE00.gdb\\ext1in100"          mxd = arcpy.mapping.MapDocument("CURRENT")          arcpy.SelectLayerByAttribute_management ("edrn_node", "NEW_SELECTION", query)          df = arcpy.mapping.ListDataFrames(mxd, "Layers") [0]     df.zoomToSelectedFeatures()     df.scale = 2000          dataframe = mapping.ListDataFrames(mxd, "*")[0]     frameExtent = dataframe.extent     XMAX = frameExtent.XMax     XMIN = frameExtent.XMin     YMAX = frameExtent.YMax     YMIN = frameExtent.YMin     pnt1 = arcpy.Point(XMIN, YMIN)     pnt2 = arcpy.Point(XMIN, YMAX)     pnt3 = arcpy.Point(XMAX, YMAX)     pnt4 = arcpy.Point(XMAX, YMIN)     array = arcpy.Array()     array.add(pnt1)     array.add(pnt2)     array.add(pnt3)     array.add(pnt4)     array.add(pnt1)     polygon = arcpy.Polygon(array)          arcpy.CopyFeatures_management(polygon, "polygon_Extent102")          arcpy.Clip_analysis(depth_1in30, "polygon_Extent102", "Depth1in30")     arcpy.Clip_analysis(depth_1in100, "polygon_Extent102", "Depth1in100")     arcpy.Clip_analysis(depth_1in1000, "polygon_Extent102", "Depth1in1000")     arcpy.Clip_analysis(extent_1in30, "polygon_Extent102", "extent1in30")     arcpy.Clip_analysis(extent_1in100, "polygon_Extent102", "extent1in100")     arcpy.Clip_analysis(extent_1in1000, "polygon_Extent102", "extent1in1000")          arcpy.Delete_management("polygon_Extent102")       I've also noticed when I convert the above to a script tool that I have to reconnect the geodatabase files before I get the expected results.


title: Error getting text attributes from shapefile with python / pyshp

text: I have a shapefile, containing only points in XY-coordinates and a string attribute for each point. I would like to read the coordinates and the attribute for postprocessing in python, but can't get it to work. I have tried the following.               import shapefile     r = shapefile.Reader("subregion_59117")     print r.fields     shapes = r.shapes()     print len(shapes)     print shapes[0].points[0]      gives               [('DeletionFlag', 'C', 1, 0), ['LON', 'N', 16, 8], ['LAT', 'N', 16, 8],     ['SUBREGION', 'C', 25, 47]]          59117          [176.0, 63.0]      So I can read the point coordinates, but               records = r.records()     sr = r.shapeRecords()     print len(records), len(sr)      gives just               0 0      meaning, these lists don't contain any elements. So I don't know how to get my hands on the string attributes of the points? However, if I open the shapefile in QGIS, and save as geojson (a JSON textfile, thus human readable), it contains the features as a list of elements such as:               { "type": "Feature", "id": 0, "properties": { "LON": 176.0, "LAT": 63.0, "SUBREGION": "Eastern Europe" }, "geometry": { "type": "Point", "coordinates": [ 176.0, 63.0 ] } }      So apparently my shapefile contains all the information, I just don't know how to read it with pyshp. (So yes, I could always convert my shapefile to geojson with QGIS, and postprocess that with python, but this feels overtly complicated.) ~~My shapefile is here.~~ **Edit:** My shapefile was generated by Saga GIS, and as explained by Tuure Laurinolli's answer, Saga GIS seems to have a bug so the attributes in my shapefile were marked with an invalid flag, and also pyshp seems to be very picky in this, interpreting everything but the correct flag as "deleted" and refusing to read it. Converting the broken shapefile to shapefile by e.g.                ogr2ogr file_new.shp file_broken.shp      fixes the problem.


title: Added PNG dissapears when I export/print map

text: Ive added a PNG image of a chart to my map and when i try to export as any image file or print the map the PNG image is not there. Im using arcmap 9.3. in layout view im adding image which ive tried as a jpg and png. it appears on the layout itself but once i print or export the map the image is no longer visible on the printout or export. Also, the image does not appear on the print preview. Ive added image files to many of my maps before so im not sure what would be causing this


title: Splitting lines automatically at specific distances using ArcGIS Desktop?

text: I have a set of polylines which I need to split at varying distances from their starting points. In general, I'm looking for a something similar to the SPLIT function within the editor tool. There it is possible to specify a certain line length where the split should be applied. For instance: I have a line A of 3200m length, that I need to split into two separate geometries, one from 0 til 1700m, and one from 1700m til 3200m. Does anyone have accomplished this using Python? Or this there no way around ArcObjects? Or might there be a different solution I'm missing? I'm using ArcGIS Desktop 9.3.1. on WinXP with Python 2.7.2.


title: Split a line from point but with a certain distance

text: I am using ArcGIS 10.2. I need some suggestion about this problem, maybe there is a certain tool, but I think no. I think some python code... but I do not know from where to start :( I have two shapefile layers: some polylines with points intersected. So, what I want to do **automatically** is to SPLIT those lines from a fixed measure (for example 10 meters) before and after this Point placed on the polyline. The Points are on the line. For example, if I have a 100 meters length polyline, and the point is situated in the middle, I want it split on the 40th and 60th meter.


title: Get the LatLon of middle points of a road segment?

text: Suppose I have the LatLon of the two end points as well as a small number of middle points of a road segment. What is the best way to calculate the LatLon of a middle point _P_ of the road segment given the road distance between _P_ and one end point of the road segment?


title: Convert shapefiles from WGS84 to Miller Projection

text: I have some shapefiles of the world from Natural Earth at **WGS84**. I would like to convert them to the **Miller Projection**. I have **ogr2ogr** , but I'm not sure what the parameters for this conversion should be, or if this conversion can be done with it.


title: Create Feature Class From XY Table Using Python

text: I have some txt files that are basically attribute tables (object ID, latitude, longitude, area). In ArcMaps, I have been converting these files into feature classes and saving them in a geodatabase using the following steps:   * Right click file Create Feature Class > From XY Table Click on the Coordinate System of Input Coordinates button and then select coordinate system Save as File and Personal Geodatabase feature classes and output to selected geodatabase Since I have a large number of files, I'd like to script this using Python. I don't have a problem iterating throuygh the file folder and getting a list of all files, but which function from the ARC Python library should I use? The is Create Feature Class tool the best one I've been able to find so far, but it created an empty feature class. i have tried opening the Geoprocessing > results tab while creating the feature class manually, so that I can copy as script, but for these actions I am not seeing any entries in the result pane.


title: How to use python to set a variable to a value from a vector attribute table for GRASS

text: I'm adding to a python script for GRASS v.transects.py. I want to set a variable to equal a value from a field in the vector's attribute table using a for loop to cycle through each row. This is one of my first GRASS scripts, I'm much more comfortable with arc. For arcpy I would use a search cursor and getValue:               vector = "river.shp"     lines = arcpy.SearchCursor(vector)     for line in lines:          attributeValue = str(line.getValue("FIELD"))      The existing script has the right for loop already (for line in vector), its just a matter of figuring out the last line, how to have it read the field in the attribute table of the line segment it is currently cycling through.


title: Overlaid raster gets shifted around when printing to pdf or when resizing graphics device with mouse

text: I've spent the past couple of days learning how to use the `maps`, `raster`, and `rgdal` packages to create nice images, but I'm running into trouble when trying to make pdf's for actually presenting. I'm overlaying two rasters -- one has my data and the other is a mask. The problem is that the mask moves around depending on the size of the output pdf, AND also with simple resizings of the graphcis device. I've made a little minimal example here:               library(raster)     library(maps)     library(mapdata)     library(rgdal)     lat = seq(from=36,to=40,by = .01)     lon = seq(from=-121,to=-117,by = .01)     d = expand.grid(lon,lat)     colnames(d) = c('lon','lat')     d$z = sin(d$lat)-cos(d$lon)+rnorm(nrow(d))     head(d)     coordinates(d)=~lon+lat     gridded(d)= TRUE     r = raster(d)     dd = expand.grid(lon[lon<120],lat[lat<37])     colnames(dd) = c('lon','lat')     dd$z = 1     coordinates(dd)=~lon+lat     gridded(dd)= TRUE     rr = raster(dd)     map('usa',xlim=c(-123,-116),ylim=c(35,41))     plot(r,add=TRUE)     plot(rr,add=TRUE,col="black",legend = F)     dev.copy2pdf(file = "example.pdf",height=8,width=8)      Here is the output when first rendered on my screen: ![enter image description here](http://i.stack.imgur.com/r63PQ.png) And here is the output when I resize the window! ![enter image description here](http://i.stack.imgur.com/qsd69.png) and the example PDF is different again! (Unfortunately SE doesn't seem to support pdf images, but if you run the code you can see it on your own machine) Any idea how to work around this behavior? Thanks!


title: Which interfaces/classes do I use to access the highlighted rows in an ArcMap attribute table?

text: I've tried to use various interfaces to access SelectionSet of table windows/tables but all I get back is the selected records (cyan). I want the highlighted rows/features (yellow). I would like to 1) programmatically specify the highlighted row/feature and 2) programmatically access the highlighted rows(s). I saw some old examples on the old Esri forums where the yellow feature was manually being drawn on the map, but surely there's a way to access in ArcGIS 10?


title: (Smooth) Isoline Extraction from TIN - Algorithms

text: Can anyone point me to a clear document or algorithm for extracting smooth isolines from a Triangular Irregular Network? I've seen some examples using the jagged linear methods, but can't seem to find good documentation on how to do the process using a smoother method like Natural Neighbor / Radial Basis / Quadric / etc... Thanks!


title: ArcGIS Surface Length Tool equivalent in FME

text: I am trying to calculate the 3D length of a line. I currently have a rather complex and cumbersome model built in ArcGIS Model Builder that uses a tool called "Surface Length". I was wondering if there is an equivalent tool in FME? Ultimately I want to calculate the 3D length of a shapefile. I have DEM data for the same area. Thanks!


title: If a shapefile with geometry esriGeometryPoint is valid is it guaranteed to have a non null featureclass?

text: I have been given a point shapefile but it's feature class is null. I haven't seen this before so I'm wondering if the data is bad or if it's something I need to handle with code.               //I was using this to determine if my layer was a point layer, but it's FeatureClass is null     public void isPointsLayer(ILayer layer)     {     IFeatureLayer fl = layer as IFeatureLayer;     if (fl == null && fl.FeatureClass == null) return false;      esriGeometryType type = fl.FeatureClass.ShapeType;     if (type == esriGeometryType.esriGeometryPoint) return true;     }      * * * edit: I guess ILayer.Valid returned false, thanks for the answers.


title: Very Inaccurate Geocodes in Canada (Quebec)

text: _Edit_ : Found this page detailing the standard geocoding abbreviation types which may be of use when trying to match foreign street types. It appears that we'll have to create a much more robust version. The problem is as listed in the title: Our current geocoding attempts (using a proprietary list of addressed and a proprietary streets network) in Quebec are coming up with abysmal match rates (40-50% overall). We have tried looking at problems with characters sets and postal code inconsistencies, but nothing seems to be nailing it down. Presuming these two factors are **not** the issue, what suggestions might be offered as possible pitfalls when geocoding in this area? Important notes: 1. Geocoding with ArcGIS 10  2. Using FSALDUs as post code in "ABC XYZ" format.  Here are some example addresses.                   ARC_Street          ARC_ZIP     CITY                           STATE     1134 BRIDGE ST W        A7R 3L1     NAPANEE                         ON     75 RUE PRINCIPALE       G0A 1A0     SAINTE-CHRISTINE-D'AUVERGNE     QC     43 RUE PRINCIPALE       G0A 1A0     SAINTE-CHRISTINE-D'AUVERGNE     QC     325 RANG SAINT-JACQUES  G0A 1A0     SAINTE-CHRISTINE-D'AUVERGNE     QC     325 SAINT-JACQUES RANG  G0A 1A0     SAINTE-CHRISTINE-D'AUVERGNE     QC     125 RUE PRINCIPALE      G0A 1A0     SAINTE-CHRISTINE-D'AUVERGNE     QC     168 RANG SAINTE-ANNE N  G0A 1A0     SAINTE-CHRISTINE-D'AUVERGNE     QC     168 RANG SAINTE-ANNE    G0A 1A0     SAINTE-CHRISTINE-D'AUVERGNE     QC


title: How to convert field values from lower case to upper case in ArcGIS 10.1?

text: Using the Field Calculator in ArcGIS 10.1, I am able convert the first letter in a value from lower case to upper case. How do you convert the entire value to upper case, though? I can't seem to locate a VBScript for it. For example, I want to convert _small_ to _SMALL_.


title: How to automatically add vertices to lines every x meters?

text: I've made a vector grid using lines; the grid coverage is 100 km x 100 km and the lines are at vertical and horizontal intervals of 1000 metres. Each line in the grid is defined by two vertices, one at each end, and each line has an ordinate attribute for labelling. My wish is to to add vertices to the lines, so that each line will have vertices at intervals of (say) 500 metres. Does anyone know how to do this? The background to this is that the grid will be displayed on OpenLayers and given the strangeness of Google Mercator the lines need to be pegged to the ground with vertices at frequent intervals. One other obvious use for the ability to add vertices to grid lines would also be to make nicely curved lat/lon graticules for projected maps.


title: Get Annotations selection Arcobjects

text: i have annotations stored in mxd document, when i select some annotation , i tried to get the selection through graphics selection, so i got always the rectangle of the dataframe (because i'm working on a axPageLayoutControl), and when i use features selection i get nothing; so none of them return the annotation selection. How could i get the annotation selection in Arcobjects ? i'm using Arcobjects 10.1


title: Parse error while inserting values into polygon feature

text: I created a polygon as follows               create table boundary (id int, name text);      table created;               select AddGeometryColumn ('',boundary', 'latlon' 4326,'POLYGON', 2);      geometry column added; then i'm trying to insert values as follows               INSERT INTO boundary (id , name, latlon) VALUES ('1', 'river', GeometryFromTerxt('POLYGON(82.240 16.988,82.243 16.991,82.240 16.988)', 4326));      finally i got an error as follow > <--parse error at position (position number) within geometry.


title: How do I automatically fill a table column?

text: I would like to solve basic problem. I need to automatically fill the table. How is in Qgis possible to fill it like : if column SPECIES contain string DOG -> write to column CODE integer 1? Thanks for any advice.


title: Which spatial weighting function does ArcGIS' GWR use?

text: I am running Geographically Weighted Regression and originally was using ArcGIS to carry out the analysis until I read/heard mixed reviews about it. I have now tried a different programme to run it (SAM) and get different results. Does anyone know what the spatial weighting function is in ArcGIS? I think it is Gaussian but am not 100% sure. I cannot select this with Arc but there are different options to chose from with SAM. Thanks!


title: How to use SQLite db with QGIS plugin?

text: Hi I would like to use SQLite for storing some information in my plugin for QGIS. I have a problem: if I put the db let say test.db among the resources as file, after compilied the resorce.qrc file, QGIS find the db, but if I insert some row in the db, qgis does not find any more I have to compile again the file, I tried to recompile the file every time I commit something in the db, it seems to work, but is a dirty solution and I do not think I can use it on every system different from the mine, is there a better solution?


title: Online map with rivers

text: Is there any internet map source (zoomable, tiled, i.e. "google like") with good display of rivers? I.e. well labelled even on larger scales? So that I'm for example able to track Donau river (with label) on scale like this.   * Basic google map is just terrible, rivers are not well displayed, often lost on larger scales, and in 90% cases not labelled.   * OpenStreetMap - similar, furthermore they denote roads as blue (!!!) which is even more confusing! (really bad practice)   * OpenCycleMap - even worse... Thanks in advance!


title: Albers Conic versus State Plane Zone in terms of accuracy

text: I'm trying to determine if an Albers Conic (Alaska Albers) is suitable as a projection for a corridor 800 miles long by 50 miles wide, that goes from the north coast of Alaska to the south coast. This corridor spans 2 state plane zones (west to east). I need to support survey grade accuracy in my GIS for the corridor, and Albers won't cut it I think. Should I create a new coordinate system composed of a merge of the 2 state plane zones to address the accuracy need? thanks,


title: Basic question about downloading SRTM 3 arc resolution data

text: I want to be able to download high resolution 3 arc degree resolution (90 meters) from NASA's earthexplorer.gov site. However it is a bit confusing for a beginner like myself. I want to be able to do this for my country i.e. India. I go here http://earthexplorer.usgs.gov/ and by pressing shift I zoom into the area of interest by creating a rectangle. Then it displays 4 coordinates and then I ask for datasets. I get so many options and each are of very small sizes. Am I supposed to download all of it or just any one of it ? If I am supposed to download all of it is there a way to do it all in one download ? I presume then I will have to use progamming code to concatenate these separate .hgt files together.


ossim-extract-vertices.exe (in commandline) not working on WIN7 op. I would like extract the image bounds from the pan-sharpened Landsat 8 image. 3 years ago it worked :-( . Have you any solution? BR Janos


title: Creating Home ranges (MCP) of marine mammals excluding solid land

text: I use ArcGIS 9.3 and the Hawth's tools for Home range analyses. I try to create MCP of marine mammals which all works fine. However, the MCP include some land as well. How can I create a MCP so that it only considers the water area? Thanks!


title: How to display places names when using Quantum GIS and shapefiles?

text: I am new to the world of geo systems. What I would like to do is displaying some information on a map for given regions. I downloaded shape files from cloud made web site and load layers using Quantum GIS. Now I would like to display the places names but I do not know how to do. Are these information already present in the file how can I do to display them ? If so, I guess It should be possible to do joins with the data I want to display. Can you please help me on this subject ?


title: GeoTIFF file creation from .Tiff file

text: im not sure if this question is suitable for GIS I'm working on a software that displays a map, like this ![enter image description here][1] i also have some export functions which exports the the about image as `.bmp/.Tiff/.Jpg/.PNG/.Wmf/.SVG` This is .TIFF export ![enter image description here][2] Now i have to export it as `GeoTiff` also The software i am working on is in `Delphi 7`. **information** i have come across this steps to creating a geotiff but since i am new to GeoTIFF , i am stuck after the 1st point. **question** so please tell me is it possible to convert the `.TIFF` file to `.geotiff`? **what i have tried** remotesensing geotiff how to write geotiff files in delphi possibly-compressed but i cant figure out how to create the meta data for the tiff to create the geoTiff


title: What exactly does the ArcGIS Server service do? SQL Server service?

text: What is the purpose of the ArcGIS Server service? It's not used at all to do a direct connect to an SDE database, right? That would just use the SQL Server service? While we are on the ArcGIS Server topic, what exactly is the windows user account that ArcGIS Server creates (default 'arcgis') used for? What is its purpose? If I create an SDE database and create a user for it, when that user connects to the database using a direct connect from his laptop- what is going on? What windows user account on the server is being used, or is that user just directly connecting to the SQL Server service over the network? Any reading material on this subject so I can better understand what's going on? Thanks for any help.


title: How to modifiy how selection highlights are displayed?

text: I have got a layer of points that are showed with a litlle symbol. When I select one, a red square comes around it. The problem is that this points are inside other layers (polygons) and some of them are red too. I would need to change the color of this square and the width, so that it becomes easily visible. I could not find anything related to this in the documentation, neither in the asked questions. Can you please help me for that? Thank you for your answer Best regards Thank you for your answer. I already tried this, but it doesn't affect the red square that remains exactly the same. My selection colour is setted to Blue under the Project's properties and the square is still red. And one thing that could be really usefull for me, would be to have the width of the square's lines increased. Do you have an idea for that ? Thank you


title: How to build a model creating shp polygonal objects using for that the parameters of xls file?

text: Could somebody advise which instrument I can use in ModelBuilder (ArcGIS10) for build a model which will create shp polygonal objects from a table (created previously in Excel) with necessary parameters (that include coordinates of latitude and longitude and distances). Some details.   I need to create many polygonal objects by 4 steps:   1. Create 2 circles with the same center (with coordinates X1, Y2 – they will be in the table) and radiuses Rmax and Rmin (they will be in the table);    2. Create a new polygonal object (“bagel”) = Circle (max) – Circle(min);    3. Create a polygon triangle by coordinates of 3 their vertexes (all coordinates will be in the table);    4. Find a resulting polygonal object (as a shapefile) that will be equal a space of intersection of “bagel” and that triangle.


title: Combining raster images?

text: I've decided to work on GIS implementation in transportation. I've these small tiles of images, with world files. I want to combine them to create a base map. I don't know coding/programming, so please suggest some other other method to combine these tiles.


title: Failed to publish globe service on ArcGIS Server. Error 001369

text: I'm trying to publish a globe service on our ArcGIS Server and in ArcGlobe I get this error message > " _Failed to publish service._ " ---> Went to Results: [h=1]001369 :  > Failed to create the service.[/h] [h=2]Description[/h]The creation of the > service failed. This is a generic error. Many different conditions could > cause this failure. Check the server logs for the specific issue. [h=2]Solution[/h]Check the server logs for the specific issue. --> Went to ArcGIS Server Logs: > SEVERE 22.08.2013, 10:15:20 Failed to construct instance of service > 'glob_4.GlobeServer'. Server SEVERE 22.08.2013, 10:15:20 Failed to > initialize server object 'glob_4': 0x80030003: null Server SEVERE > 22.08.2013, 10:15:20 (null) glob_4.GlobeServer Unfortunately I really don't know how to solve that. The service is published (I can see it on the server) but apparently it cannot start. I also tried saving a service definition file and then publishing it on the server-> same problem. So thanks for any advice. PS: We use ArcGIS Server 10.1 on Linux


title: How to generate features to test different implementations of LinearRing?

text: I have some code that throws the following exception:               Points of LinearRing do not form a closed linestring      and I want to check if my code is calling this correctly by supplying some features that either violate or do not violate the constraints of LinearRing.validateConstruction. Is there an existing Java utility that can be used to generate appropriate feature data for this?


title: Agent Analyst extension for ArcGIS

text: I am trying to create agent based model for simulating movement of occupants in different rooms of different floors of a building. I am using Agent Analyst tutorial downloaded from ESRI website and PhD dissertation of Elizabeth Groff as a guide. Since the building has three different floors, I have created a 3D network dataset and created a toolboox which can find shortest route, closest facility and OD matrix from room centre to all the emergency exits. I have created the following according to the two guides above.   1. csv file containing the activity space (i.e. the door of the room occupied, nearest exits of the floor, bends on the staircase and exit of the ground floor).   2. Point shapefile of path that the agents have to follow.   3. csv file containing the list of all neighbouring nodes.   4. Point shapefile of occupants. There are more than one occupants in few rooms. I intend to show the movement of these occupants as vector agents in Agent Analyst.   5. Point shapefile of all emergency exits. I also want it to be include in the model. So I want to:   1. Show the movement of these agents from their present location to the nearest exit on ground floor.    2. The occupants are of different sex and age, thus I want the agents to move at different speed from one node to another and while they are moving downstairs. How can I do that? Lastly, can anyone tell me if movement of agents can be shown in ArcScene because I tried to show the movement, it does not. I tried the same in ArcMap and it works. I haven't tried the same in ArcGlobe. I am using ArcGIS 10 and Agent Analyst extension for ArcGIS. Thanks in advance. ......


title: Which Graduated Symbology to use in QGIS

text: I am trying to map demographics data for a county in the UK based upon a combination of Census data. I want to compare the the averages for the county alone to the averages of the UK as a whole and display it using a sort of traffic lights symbology (red, orange, yellow, green, where red is way above average and green is way below average). **My problem is I'm not sure about what type of graduated symbology to use. Can someone please explain what Graduated Symbology is best to use in my situation?** I used Natural Breaks (Jenks) at first and was very happy with that but then I read on the ESRI website: " _Natural Breaks (Jenks) are data-specific classifications and not useful for comparing multiple maps built from different underlying information_." This makes me think using Jenks and then comparing the county averages to the national averages or comparing different categories (health data versus accessibility data) isn't appropriate. So then I thought it would be best to use Equal Interval or Quantile symbology or define my own bands (split the results in four so I can assign colours to each 25%). Any advice/help/information would be great. ![enter image description here](http://i.stack.imgur.com/MUmDP.jpg)


title: Creating management zones in NDVI ortho maps with contour like lines

text: This question is hard to explain so if you have a possible answer I would love to hear it. I work in the precision agriculture business which heavily relies on NDVI calculations but need to create like zones of crop health for field surveyors and field management purposes. The end users for maps I create are not GIS savvy so presentation is everything. Raster Classification, zonal fills, and reclassification's are not what I want because they are vague and distort the data. On the attached map there are geo-coordinate points over areas with possible health strains and areas that appear to be very healthy. **How do I create zone lines, almost like contour lines, to batch similar points into like areas?**   For example, in the Northwest portion of the map, the three points are over similar areas. How do I create a boundary lines showing different zones?   I'm looking for 5 areas:   i) worst;   ii) bad;   iii) OK;   iv) good;   v) best,   separated by contour like lines creating zones. I don't want to disrupt are high resolution photography. ![enter image description here](http://i.stack.imgur.com/gP0sg.png)


title: Can QGIS be installed on a Server

text: My question is whether QGIS can be downloaded and installed on a server for multi-user access. Rather than installing the software locally, can QGIS be stored on a server and accessed by many?


title: Qgis install can't find Gdal framework (Mac OS X 10.7)

text: After having installed Mac OS X 10.7 I can't install Quantum Gis as it maintains that I have not installe the GDAL Framework, which I have done (several times). Have anyone experienced the same problem (and solved it)?


title: What is the most common way of displaying geodata from Postgis on Leaflet?

text: For example, in my database I've some points with geographic coordinates. I want to display them on Leaflet JS layer. Which instruments should I use? Maybe there is some API...


title: Where can I get a multipoint shapefile?

text: I'm trying to find a multipoint shapefile to do some testing for an application I'm working on, but I cannot seem to find one (and do not have ArcGIS to create one). Can anyone suggest a place to find a sample multipoint shapefile?


title: What's the usual term for "ground plan based"

text: I don't know if "ground plan based" is the right word in the following context: A computer program reads a ground plan and produces a result ground plan. The result groundplan's deviation from the input ground plan must be within a prescribed tolerance to be called "ground plan based". Possible words for this situation might be: ground plan based, topographic planimetry, planimetric reconstruction But I'm not sure which is the common word? Which word would you use to describe that both of the ground plans match nearly 100%?


title: Anyone have success with OpenLayers and ArcGISCache?

text: * OPENLAYERS Has anyone had success using OpenLayers and the ArcGISCache layer beyond the examples? I'd like to mix tiles with **_different resolutions_** and although they seem to scale correctly they do not translate correctly. Code to translate the tiles appears to be missing with the exception of a % translation necessary after the scale:               transformDiv: function(scale) {              // scale the layer div              this.div.style.width = 100 * scale + '%';         this.div.style.height = 100 * scale + '%';              // and translate the layer div as necessary              var size = this.map.getSize();         var lcX = parseInt(this.map.layerContainerDiv.style.left, 10);         var lcY = parseInt(this.map.layerContainerDiv.style.top, 10);         var x = (lcX - (size.w / 2.0)) * (scale - 1);         var y = (lcY - (size.h / 2.0)) * (scale - 1);              this.div.style.left = x + '%';         this.div.style.top = y + '%';     },        * ESRI It appears that ArcGIS JSAPI has an even greater limitation...all the tiles must have exactly the same resolution. Can someone from ESRI discuss/confirm this limitation?


title: How to project geographic coordinates (lat/long) to a cartesian system (albers equal-area)?

text: I am using this matlab package http://www2.ocgy.ubc.ca/~rich/map.html#examples to project the latitudes,longitudes of my data to a specific projection. For this particular case I am using albers equal-area. Now when I do the projection, I am using the following code               m_proj('albers equal-area');     [x y] = m_ll2xy(longitude,latitudearray);      My question is when I specify the projection method, do I need to specify the geographical bounds. I mean do I need to pass the latitude, longitude bounds to the function m_proj I have got the following results. I am not sure if these points are correct.               lat                   lon                  x                      y     24.590960444058823    -119.6310020520101   -0.0778982228028     -0.3893279      Can anyone provide some suggestions?


title: Geoserver with ECW image file

text: I am trying to show ECW ortophoto as base layer in openlayers.   So i found this thread on GIS SE:                   Geoserver 2.1.1 with ECW (v3.3) support          1. Download Geoserver 2.1.1 from http://downloads.sourceforge.net/geoserver/geoserver-2.1.1-bin.zip 1a. extract to c:\bin\geoserver          2. Download imageio 1.1.1 installer from http://java.net/projects/imageio-ext/downloads/download/Releases/ImageIO-Ext/1.1.x/1.1.1/windowsInstaller/windows32-imageio-ext-installer-gdal-mrsid-ecw-1.1.1.zip 2b. Install it to c:\imageio 2c. Remove c:\bin\geoserver\webapps\geoserver\WEB-INF\lib\imageio-ext-1.0.8 2d. Copy c:\imageio\lib*.* to c:\bin\geoserver\webapps\geoserver\WEB-INF\lib (imageio installer will also copy gdal dll, ecw sdk dll to your java\bin directory)          3.Add Environment variables for GDAL_DATA c:\imageio\gdaldata GDAL_DRIVER_PATH c:\imageio\gdalplugin          4.Start geoserver      I did all this steps except number 3, what its means to Add Environment variables?   How I do that?


title: How to convert WKB to radians or utm in postgresql?

text: I am using PostgreSQL with postgis, the geometries are in UTM and stored as WKB with several shapes types (point, polygon, linestring...). Is there a way to convert the WKB to a list of points in radians or utm without the string parsing of the shape type and each point? Thanks.


title: OpenLayers - feaures from WMSGetFeatureInfo are not shown on layer (but they exists)

text: I have a project with google as base layer and geoserver serving WMS layers (structures - multipolygons) that are in EPSG:2039 The layers are showing fine on google base layer (I guess OL re-project them). Problem is , once I do WMSGetFeatureInfo and retrieve all the features (as in this example http://dev.openlayers.org/releases/OpenLayers-2.13.1/examples/getfeatureinfo- control.html) I don't see them on the map , the layer does have the features because I used firebug to query it (highlightLayer.features.length) but they coords are not in EPSG:900913 and I have no idea how to make them show. This is my code :                   new OpenLayers.Control.WMSGetFeatureInfo({                 url: geoserver,                 infoFormat: "application/vnd.ogc.gml",                 name: "infoControl",                 maxFeatures: 5,                 queryVisible: true,                 eventListeners: {                                getfeatureinfo: function (e) {                         if (e.features.length > 0 )                         {                                             highlightLayer.removeAllFeatures();                                                 highlightLayer.addFeatures(e.features);                             highlightLayer.redraw();                         }     });      I've tried reprojecting the feature before adding it to the layer by iterating e.features:                      Ext.each(e.features, function(feature) {                                       highlightLayer.addFeature(feature.geometry.transform(epsg2039,epsg900913));                                         });      but without any success. Is there any solution I can use in code ? (I don't want to use different base layer not reprojecting my geoserver data to 900913). Thank you !


title: How to visualize azimuthal data with uncertainties?

text: I am trying to make a figure showing azimuthal data with a different range of uncertainties at each point. This oldschool figure from a 1991 paper captures the "bowtie plot" idea that I'm aiming for: ![From Hillhouse and Wells, 1991. "Magnetic Fabric, Flow Directions, and Source Area of the Lower Miocene Peach Springs Tuff in Arizona, California, and Nevada"](http://i.stack.imgur.com/Q12Ns.png) Any suggestions on how I could make a similar figure? I'm a relative newbie at GIS stuff, but I do have access to ArcGIS through my university. My Arc experience has been limited to making geologic maps, so I haven't had to do anything too exotic. I've poked around in the symbology options in Arc and QGIS but haven't seen any settings that I thought would do the job. Note that this isn't just a matter of rotating bowtie-shaped symbols by azimuth; the angular range of each "bowtie" needs to be different. I'd rate my Python skills as 'strong intermediate' and my R skills as 'low intermediate', so I'm not averse to hacking something together with `matplotlib` and `mpl_toolkits.basemap` or similar libraries if necessary. But I thought I'd solicit advice on here before going down that road in case there's an easier solution from GIS-land that I just haven't heard about.


title: Help needed for georeferencing in QGIS

text: I am new to Quantum GIS. I just happened to browse to a tutorial on georeferencing in QGIS and I found it easy and interesting. I also calibrated the map and saved the points as per WSG84 standard. Now after I closed the software and when I restarted I am unable to see the coordinates when I load the GeoTIFF image. Can you please guide me where I am going wrong? Why am I not able to load the coordinates? How can this problem be solved?


title: can't add JPG files into QGIS?

text: I can't add JPG files into QuantumGIS. I get an error message "C:\Spatial_Days\Data\Working\AP_Coast_Export.jpg is not a supported raster data source". I am using QGIS 1.7.4 Wroclaw with GDAL 1.8.1 installed from OSGEO4W. Any ideas? Are som JPG libraries in the wrong place? Regards Richard


title: How should I get a Raster out of ArcMap and into Google Earth (.kmz) with transparency

text: RasterLayerExportClass will only create a .png with a single value transparency, however there is a bug in Google earth which prevents this from working. The other export options don't support transparency at all? (jpg, tiff) I do notice that the 3D analyst toolbox performs this correctly with the "Layer to KML" function, this I have no source for.


title: Identify primary routes in Postgis

text: I have a table of linestrings that make up a route network. I want to identify the primary routes and separate out the interconnections. Problem is that the primary routes are split into separate linestrings and have no ID to identify them separately. The diagram below illustrates the problem. There are three major routes parallel to each other and two intersections. I want to drop out the intersections and create three single linestrings to represent the three major routes. Then the intersections would be retrieved as separate geometries. I had thought that I could look for vertex's with 3 vertices connections at both ends and drop those out, but this can happen on a major route line. Also the highlighted interconnection is made up of two linestrings, so would not fit into this rule. ![major route problem](http://i.stack.imgur.com/qbMFv.png)


title: ArcGIS - Table columns lose formatting when I make a new copy of shapefile (Decimals, Thousands separator, Highlighting)

text: When I make a copy of a shapefile or I join two shapefiles and then create a subsequent shapefile, my number type columns (in the attribute table) will default back to numbers with many decimals and no thousands separator. This after I've cleaned up every column the way I like it. I also lose the highlighting of any columns I have a highlight on. I'm just looking for a way to avoid having to reformat every column if I happen to make a copy of the shapefile. That or is there a way to set a default of 2 decimal places and a thousands separator for the attribute table? I'm trying my best to get away from shapefiles, but we need it to transfer between multiple programs. Any ideas would be appreciated.


title: Hardware and software for capturing bathymetric lake elevations?

text: What hardware and software would you use for capturing bathymetric elevations (depth) in inland waterbodies (lakes and rivers)? Your device and software budget is $1,000. The budget is a soft target, to get setup, and doesn't include actual data capture. A 17' open boat with outboard motor and a ruggedized laptop running Windows XP Pro and Arcmap (Arcview license) is at your disposal, if you need them. The data captured should be suitable for building a raster elevation model and bathymetric contours. Depths are not expected to exceed 200m. Desired accuracy should be comparable to or better than a standard consumer handheld GPS (Garmin GPSMap, etc.). _(Actually the budget is a sham. If you've got a good solution for $5k or $50k I'd like to know about it too.)_


title: Finding lat/lon of corners of geotiff file

text: I am a GIS newbie. I just recently managed to turn a GeoTIFF file into a bmp file using gdal_translate. Now I need to find the latitude/longitude coordinates of the corners. I would prefer to use more of the gdal code or utilities, but that is not essential. Does anyone have a good way to do this?


title: Where can I find equipment that will allow me to recreate a form of Google Streetview?

text: I am looking for more affordable equipment that would allow me to recreate the concept of Google's Streetview, as well as software that would make a comparable UI map (without it "belonging" to Google, ie. not having everything available via Google Earth). What I would like to create are scientific maps along rivers, streams, and wetlands, with sampling data (temp., PH, TDS, dissolved oxygen...) available to the user along the path in various increments. Essentially I am looking for camera equipment that can take 360° shots, that can record the images with X,Y values that can be input to a digital map and brought up on screen the way Street View looks. I know this is more of an digital imaging question, but it will be set up through a GIS system for user interface. Thanks!


title: what is the suitable projection coordinate system for Iran

text: I am doing analysis for province wise Iran. I am using WGS1984 as my geographical coordinate system, and I do not know which is suitable Projection coordinate system for entire Iran. I feel the UTM will not be suitable because the country shares several zones in UTM system.


title: How to plot Storm Prediction Center (SPC) Shapefile using python

text: I am trying to make a plot of the SPC's thunderstorm outlook using Python. I have downloaded the shapefile from this website http://www.spc.noaa.gov/products/outlook/day1otlk.html. Then I use this code to try to read the data.                   from mpl_toolkits.basemap import Basemap         #draw basemap         m=Basemap(projection='mill',llcrnrlat=20,urcrnrlat=50,llcrnrlon=-130,\           urcrnrlon=-60,resolution='c')         m.drawcoastlines()         m.drawcountries()         m.drawstates()         m.drawmapboundary(fill_color='#FFFFFF')         shp_info = m.readshapefile('..\day1otlk_20140709_2000_cat','thunder')      But then I get ValueError: Shapefile must have lat/lon vertices - it appears this one has vertices in map projection coordinates. Convert the shapefile to geographic coordinates using the shpproj utility from the shapelib tools. I am pretty new to this but I understand I need to convert this shapefile to geographic coordinates but the error message is too vague and I have no idea how. Any help would be appreciated. Thanks.


title: How to calculate Network Service Areas in QGIS?

text: Is there any way in QGIS to apply service area buffers e.g if I have a point on a map I wish to show the service coverage of 1 mile, 2 mile, 3 mile, breaks using the road network?


title: How are bound variables specified in an ArcGIS REST SOE?

text: In the documentation for using handler functions with ArcGIS Server Rest SOE's, they don't explain how the keys in the boundVariables collection are determined. When I search through the entire solution for "customLayersID" I don't see anything. How would I know to use "customLayersID" as a key?               private byte[] CustomLayer(NameValueCollection boundVariables, string outputFormat,         string requestProperties, out string responseProperties)     {         responseProperties = null;              //The LayerID.         int layerID = Convert.ToInt32(boundVariables["customLayersID"]);              //Execute.         CustomLayerInfo layerInfo = GetLayerInfo(layerID);              string json = layerInfo.ToJsonObject().ToJson();              return Encoding.UTF8.GetBytes(json);     }


title: How to merge/union non contiguous polygons using Oracle Spatial functions

text: Is there any way of merging non contiguous polygons using the SDO_AGGR_UNION or a similar dissolve function? When i run the below query in oracle sql developer i managed to convert geometry type from polygon(2003) to multipolygon(2007), but in a gis system some polygons are not included on preview. `select shape from ( select nationalcode as nt,SDO_AGGR_UNION(SDOAGGRTYPE(shape, 0.5)) as shape from auadmunits where nationalcode = '1225' group by nationalcode) ;` I know that union is for adjacent shapes, but i am wondering if anyone else has faced this issue and has found a approach, not only handling oracle spatial. Thank you.


title: Guideline/scheme for naming properties in GeoJSON

text: Are there any guidelines/requirements for naming properties in GeoJSON? Should I add a vendor-prefix for my properties which only make in use with a specific application? I already read the specification however I didn't found any further information on this matter. **Examples:**               {         "foobar-id": 21     }      or               {         "com.foobar": {             "id": 21         }     }


title: Test for polygon orientation

text: I'm want to merge multiple geometries with Geotools (starting off with one geometry) but the result is a polygon that has a clockwise orientation. This gives a invalid geometry. My code for the merge:                   GeometryCollection geometryCollection =         (GeometryCollection) factory.createGeometryCollection(list.toArray(new Geometry[list.size()]));         Geometry all =  geometryCollection.union();      Why would something give an invalid geometry? And what can I use to test for orientation? Thx


title: How to switch to edit mode vector layers, in pyqt application based on qgis?

text: I do QGIS based application, I need to implement the ability to modify the geometry of vector layers in manual mode (editing nodes), how to do it?


title: OpenLayers legend graphics not using proxy

text: I'm using the OpenLayers/geoserver stack and my WMS requests are configured to go via a proxy using the following:               OpenLayers.ProxyHost="../gis/geoserver.seam?url=";      So far so good, all future WMS request are now channelled through this proxy script and WMS images are generated/returned as expected. My problem arises when viewing the legend, the request to get the layer image (WMS request) ignores the OpenLayers proxy definition and tries connecting direct (client browser->geoserver server which is not externally accessible - hence the proxy). Does anyone have any idea how I force the legend image requests to use the proxy? Thank you!


title: Adding attribute list to combobox in a python plugin

text: I am developing a simple plugin at the moment which contains 3 combo boxes:   * comboBox_1   * comboBox_2   * comboBox_3 Combobox 1 & 2 lists the vector layers available, however I would like the 3rd combobox to list the attribute column that both layers have in common (from the attribute table).                   # run method that performs all the real work     def run(self):         # show the dialog         self.dlg.show()         for layer in self.iface.legendInterface().layers():             if layer.type() == QgsMapLayer.VectorLayer:                 self.dlg.comboBox.addItem(layer.name())                 self.dlg.comboBox_2.addItem(layer.name())          Any tips/advice would be most welcome as I am new to this and followed various tutorials and the python cookbook.


title: QGIS - Atlas. Outputs missing

text: This really teases me. Somehow the Atlas-function in Print Composer are missing some of my polygons from my layer, when doing a print. I´ve looked at them once again,- they are there, and the geometry is right, but not all them are printed. How to solve this?


title: Updating Arc Toolbox paramaters with validation script

text: I'm simply trying to disable parameter 1 if parameter 0 if entered. I modified the following, but it does not work, do I need to return anything?               def updateParameters(self):       """Modify the values and properties of parameters before internal       validation is performed.  This method is called whenever a parmater       has been changed."""            if self.params[0].value == False:           self.params[1].enabled = False       else:           self.params[1].enabled = True            return [self.params[0],self.params[1]]


title: Proj4 Task not working - Northing Easting to Lat Lng

text: I usually post in StackOverFlow - if this is the incorrect area to post - I apologise. I have trying to convert a lot of northings / eastings to lat lng I have code below working: In ruby using this gem https://github.com/dark- panda/ffi-proj4                class PostCodeNode < Node       require 'ffi-proj4'        attr_accessor :postcode, :sector, :district            def self.from_northerns_and_eastings(northing, easting, postcode)          srcPoint = Proj4::Point.new(easting, northing)          srcProj = Proj4::Projection.new('+proj=tmerc +lat_0=49 +lon_0=-2 +k=0.9996012717 +x_0=400000 +y_0=-100000 +ellps=airy +datum=OSGB36 +units=m +no_defs')          dstProj = Proj4::Projection.new('+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs')          dstPoint = srcProj.transform(dstProj, easting, northing, 0)          lat, lon = to_degrees(dstPoint.x), to_degrees(dstPoint.y)          new lat, lon, postcode       end             def initialize(lat, lon, postcode)          @postcode = postcode.gsub(/\s+/, '')          @sector   = if @postcode.length == 7 then @postcode[0..4] else @postcode[0..3] end          @district = @postcode[0..2]          super(lat, lon)        end       end     end      However, the lat / lon that I am out putting end up in the middle of the sea! So I know I am doing something very wrong!                    dstPoint = srcProj.transform(dstProj, easting, northing, 0)      This line in particular is where the "magic" happens! In the test library on the git repo the test specifies some unknown numbers to me, I am not sure if I have my input data correct! https://github.com/dark-panda/ffi- proj4/blob/master/test/transformation_tests.rb               def setup         @proj_wgs84 = Proj4::Projection.new(["init=epsg:4326"])         @proj_gk    = Proj4::Projection.new(["init=epsg:31467"])         @proj_merc  = Proj4::Projection.new(["proj=merc"])         @lon =  8.4293092923         @lat = 48.9896114523         @rw = 3458305         @hw = 5428192         @zw = -5.1790915237       end            def test_gk_to_wgs84         to = @proj_gk.transform(@proj_wgs84, @rw, @hw, @zw)              assert_in_delta @lon, to.x * Proj4::RAD_TO_DEG, TOLERANCE         assert_in_delta @lat, to.y * Proj4::RAD_TO_DEG, TOLERANCE         assert_in_delta 0, to.z, TOLERANCE       end            def test_wgs84_to_gk         point = @proj_wgs84.transform(           @proj_gk,           @lon * Proj4::DEG_TO_RAD,           @lat * Proj4::DEG_TO_RAD,           0         )         assert_equal @rw, point.x.round         assert_equal @hw, point.y.round         assert_in_delta @zw, point.z, TOLERANCE       end      There is no documentation with this ruby gem so I am doing a little bit of guess work, but I think I am nearly there. Can someone see my error, so I can return the right lon/lat. I am writing this in Jruby so any better Java Libraries could also be used :) Thanks!


title: Why does IDatasetContainer2.CreateDataset fail with ArcObjects at ArcGIS 10.1?

text: I have a small application used to check and generate network, created with the help of the following tutorial: http://help.arcgis.com/en/sdk/10.0/a...00000w7000000/ It worked just fine with ArcGIS 10, but since I had to upgrade to 10.1, the application fails on this line:               INetworkDataset networkDataset = (INetworkDataset)datasetContainer2.CreateDataset(deDataset);      The error code is "System.Runtime.InteropServices.COMException was unhandled Message=Exception from HRESULT: 0x80042262", which really isn't helpful. Is there something new I could have missed, or some help on what could cause such an error ?


title: Select features that do NOT intersect in ArcGIS (10.0)

text: I used the "create spatially balanced points" in Geostatistical Analyst in ArcGIS 10 to generate a large set of points with equal inclusion probability but a "spatially balanced" dispersal. I then buffered these using different buffer settings, with the idea of explicitly testing the influence of sampling unit size on population estimates (in my case, tree cover) using a fixed area sampled. This latter detail is key - I sum of the area for the sum of the selection of buffered points to be the same between comparisons. A toy example: what's the difference sampling 10 polygons 1 ha each in area (10 ha total) versus 100 polygons, each measuring 0.1 ha (again, 10 ha total). The point layer created by the tool qualitatively has the desired dispersion, but when you start buffering points, some start to overlap, reducing the total area (plus violating any possibility of independence between samples). What I'm looking to for is an efficient way to find features that DO NOT intersect. I saw a post asking a similar question for PostGIS, but was looking for a way to approach the question using ArcGIS. Any help is greatly appreciated!!!


title: How to limit input of Python script parameter to polyline only when using Layer file?

text: When the parameter is feature class I can filter the input to polyline only. How to do this when the input is **layer file** at CURRENT MXD? I'm looking for something that will operate before the script starts running. I'm using Layer file so the user can choose the input from the data frame at the current MXD.


title: PostGIS raster functions aspect and slope slow

text: I am trying to compute the raster and slope values from a DGM table (8184 x 5159) imported from tif file tiled in 500x500 tiles. In my query in which I merge the DGM table with a vector (point type) table consisting of 3000 (or 15000) points this query takes extremely long. Explain analyze tells me for a single point to compute aspect and slope it takes more than 3s. I set on both tables an index (gist (st_convexhull(rast)) on the DGM and a gist(geom) on the point table). Any idea how to speed up this query? Markus explain analyze select st_value(st_aspect(dem.rast,1,'32BF'),sat.geom,true) as aspect, st_value(st_slope(dem.rast,1,'32BF'),sat.geom,true) as slope from safer.st_dgm_tiled dem, safer.aoi_bolz_rsat_a sat where ST_Intersects(dem.rast,sat.geom) limit 1; ## And this is the query plan Limit (cost=0.00..1101.35 rows=1 width=64) (actual time=3332.439..33051.935 rows=10 loops=1) -> Nested Loop (cost=0.00..1101.35 rows=1 width=64) (actual time=3332.434..33051.913 rows=10 loops=1) Join Filter: _st_intersects(dem.rast, sat.geom, NULL::integer) -> Seq Scan on st_dgm_tiled dem (cost=0.00..3.87 rows=187 width=32) (actual time=0.103..0.119 rows=59 loops=1) -> Index Scan using aoi_bolz_rsat_a_gist on aoi_bolz_rsat_a sat (cost=0.00..3.35 rows=1 width=32) (actual time=0.006..0.007 rows=0 loops=59) Index Cond: ((dem.rast)::geometry && geom) Total runtime: 33065.589 ms (7 rows)


title: Feature attachments in Geoserver: Possible?

text: Is is possible to have a feature attachment functionality with Geoserver services, similar to what ArcGIS server offers in its feature service? If Geoserver does not have a built-in functionality to support add/remove/download of feature attachments, could anyone suggest a way to do so? Additional info: Feature attachments are the files that you want to associate with some features. For instance, geo-tagged images could be attached to corresponding point features, or some textual reports could be attached to land parcels. This functionality is a built-in for some formats including local and ArcSDE geodatabases. That is, the attachments reside in the geodatabase, and the relation is not provided by hyperlinks. This approach could provide better consistency between features and their attachments. ESRI products provide means for add/remove/query/download of feature attachments. For instance, your web-app can interact with feature attachments through ArcGIS server's REST interface.


title: How to create a circle vector layer with 12 sectors with Python/pyqgis?

text: From a given coordinate and radius I need to create a circle with 12 sectors. The sectors should start from north (0 degrees) in 30 degree steps clockwise. I could also imagine to create the sectors by triangles, if that would be easier.


title: Using the esri analysisTool Erase class in C#

text: I am trying to use eSRI.ArcGIS.AnalysisTools.Erase class but I keep getting blank layers. I can get the layer that can be used for the input and the layer for the erase. I don't really understand the output value though. Can that be a layer? Is it supposed to be a featureclass? If it is a layer or featureclass is it modified when Erase is called? Documentation on this feature is rather limited and doesn't give much of any explanation of how this class works. Is there a place where I could get better examples of how to use the Erase class in C#?               Erase erase = new Erase(inputLayer, intersectionsLayer, erasedLayer);      Could I do the following in order to modify the input layer:               Erase erase = new Erase(inputLayer, intersectionsLayer, inputLayer);      input layer is a polygon layer, roads or rivers on a map. intersectionsLayer is the intersection of the roads and rivers on the map erasedLayer is supposed to be the output.


title: how can I find sDNA toolbox for arcgis?

text: I am working on spatial structure of the city and security. I need the sDNA Toolbox to add in gis toolboxes. please help.It is a real emergency.


title: Where can I find historical political affiliation data for the United States?

text: I am collaborating on a research paper where one of our variables is political affiliation across the United States. I will be creating a series of historical maps. I'm looking for the percentage of the population that is Republican, Democrat, and Independent in each zip code in the US. If zip code isn't available, county would work too, or statewide if nothing else. I've searched online and have only found 2008 and 2012 election data. Does anyone know where I could find this historical political affiliation data? I need data from 1980 to present day.


title: Small issue with field calculator

text: What am I doing wrong in this field calculator logic? Pre-logic               def re_score(elev_t,f2, contour):      f2=contour      if elev_t.endswith('.5'):       f2 = "Midi"      return f2      Expression for field "Contourtyp"                re_score(!Elev_t!,f2,!Contourtyp!)      I just want to rename the contour type where the contours are at half meter intervals. The 'Elev_T' is a text column with all the elevations (converted to text as I couldn't think of a way to do it with numerical elevation data which has .25, .5, 75 and 1m intervals). I want all .5 called 'midi' and .25/.75 remain as 'minor' and 1m remains as 'basic' > Executing (Calculate Field (3)): CalculateField 328000_8123000 Contourtyp " > re_score(!Elev_t!,f2,!Contourtyp!)" PYTHON_9.3 "def re_score(elev_t,f2, > contour):\n f2=contour\n if elev_t.endswith('.5'):\n f2 = "Midi"\n return > f2" Start Time: Thu Apr 12 17:10:40 2012 ERROR 000539: : unexpected indent > (, line 1) Failed to execute (Calculate Field (3)). Failed at Thu Apr 12 > 17:10:40 2012 (Elapsed Time: 0.00 seconds) I have tried different indentation and no luck... Thanks


title: Dataset of unknown islands?

text: Recently I've reviewed ASTER GDEM v2 dataset and found many tiny (at least > 1ha) islets that's missing in other dataset such as WVS. I've change the altitude to 10m and they're still there. Anyone/Organization ever attempted to documented/confirmed those? Or those are just defects? I can't see them from Google Maps as the area isn't part of interest (no data). Note: - ASTER itself limit their scope (1 square degrees) to area with at least a significant island. So there are possibility of unknown islets everywhere else. - The area that I reviewed is not area of polar ice extent.


title: System Requirements for 3d GIS & 3d Modelling

text: I'm looking for system (hardwares) required for view, edit & analyse 3d GIS & modelling data. And what is the decent monitor size should be ? single/Dual screens ? Modern systems comes with 8GB Ram, 2GB Graphic Card, 1TB space, 3 button mouse. But do we need anything else for a workstation ? which has decent capability for 3d GIS & 3d modelling data ? And do you reckon 2 different set of systems for viewing and editing capable workstation ? your help is really appreciated. Thanks Bas


title: Can GDAL Retile run multi-threaded?

text: gdal_retile is currently using only one CPU of the available eight and it is taking a while (48 hours already) to retile 5TB of data.


title: How to get WKT or GEOJson by country or city name with Overpass API?

text: I'm trying to understand Overpass API, I'm experimenting with it, but cannot get what I expect. Overpass API doesn't return me neither polygon shapes, nor WKT, nor GEOJson, nor even a highlight of the area I'm querying for. My goal is to get either WKT or GEOJson of current OSM relations (exactly polygons (a set of lng/lat points)) for countries and cities. For example, I want to get a polygon for France. I make these queries in http://overpass- turbo.eu/:               relation         ["name"="France"]         ["coverage"="polygon"]         ["type"="boundary:administrative"];     out body;      The response I get, no any polygon lat/lng data:               <?xml version="1.0" encoding="UTF-8"?>     <osm version="0.6" generator="Overpass API">     <note>The data included in this document is from www.openstreetmap.org. The data is made available under ODbL.</note>     <meta osm_base="2013-11-08T11:43:01Z"/>               </osm>      The second version of my request, I provide only a name:               relation         ["name"="France"];     out body;      Now I get dilated information about the country, but no polygons either:               <?xml version="1.0" encoding="UTF-8"?>     <osm version="0.6" generator="Overpass API">     <note>The data included in this document is from www.openstreetmap.org. The data is made available under ODbL.</note>     <meta osm_base="2013-11-08T11:37:02Z"/>            <relation id="11980">         <member type="node" ref="17807753" role="admin_centre"/>         <member type="relation" ref="1362232" role="outer"/>         <member type="relation" ref="1263863" role="outer"/>         <member type="relation" ref="1263868" role="outer"/>         <member type="relation" ref="1628087" role="outer"/>         <tag k="ISO3166-1" v="FR"/>         <tag k="TMC:cid_58:tabcd_1:Class" v="Area"/>         <tag k="TMC:cid_58:tabcd_1:LCLversion" v="8.00"/>         <tag k="TMC:cid_58:tabcd_1:LocationCode" v="2"/>         <tag k="admin_level" v="2"/>         <tag k="currency" v="EUR"/>         <tag k="flag" v="http://upload.wikimedia.org/wikipedia/commons/c/c3/Flag_of_France.svg"/>         <tag k="land_area" v="administrative"/>          ...... and so on ......      How to get a set of lng/lat (WKT or GEOJson) of the territory (countries and cities only) by their name?


title: ogr2ogr DXF to PostgreSQL some LINESTRING coordinates mirrored to the negative longitude

text: I have been using **ogr2ogr GDAL tool** to migrate geospatial data from DXF file to PotgresQL database without much problems. But now there is a special problem I can not fix. I have a **DXF file with LINESTRINGs** that when I import them **to the PostgreSQL database** some of the linestrings are imported in the correct geometric coordinates but **some other LINESTRINGs are imported with the negative longitude** , so that LINESTRINGs are mirrored to the west. When I open the DXF file with my QUANTUM GIS viewer I can see all the LINESTRINGs together with the correct coordinates, but when I export them to the postgres database and I print them in a map (by a web application) or later I export then I see some of the LINESTRINGs mirrored "away home" from the original coordinates. That is my ogr2ogr commandline:               ogr2ogr -select "" -f PostgreSQL "PG:host=192.168.*.* user=**** password=**** dbname=****" C:\PATH\FILE.dxf -s_srs "EPSG:23031" -a_srs "EPSG:23031" -lco GEOMETRY_NAME=the_geom -nlt LINESTRING      There is no error output from the ogr2ogr tool, so it is difficult to get the reason of that problematic. Any help or recommendation? Thank you very much. **EDIT:** I have tried without much effect the ogr2ogr opions: -lco GEOM_TYPE=geometry -lco DIM=2 --config PG_USE_COPY YES -lco PRECISION=YES -clipsrc xmin ymin xmax ymax


title: How to convert between KML and Esri's shapefile (SHP) format?

text: My users are sending me point data that were digitized using GoogleEarth. How can I convert their KML to a shapefile?


title: Good KML conversion program

text: I recently used the GPSlogger app and created KML files from my routes. I haven't been able to successfully convert them to shape files. Does anyone know a good downloadable program that converts file types?


title: How can I convert Shapefiles (SHP) to KML?

text: > **Possible Duplicate:**   >  How can I convert KML to ESRI's shapefile format? I have SHP files, I want KML files, what tools can convert?


title: qgis build under ubuntu 12.04, but all the python plugins are not there

text: All the .so plugin show correctly, only the python plugin do not. I am confused and do not know how to fix this. I think there must be some python problem . The build process is successful.


title: How to avoid or deal with "unable to save" errors within ModleBuilder ArcGIS (10.1)?

text: Do you have references to best practice, or examples where you have recovered Models or can you share your expertise on how to avoid the "Unable to save model " error within Model Builder. Alternatively can you advise on situations to avoid which cause this error to occur. There is very little help information on this topic, so any shared advice would be much appreciated. There is a linked question here Do you use a proven workflow to edit and document ArcGIS ModelBuilder models? - that shares advice on how to track the building of models, so you always have a clean copy of the model to go back to. Please comment there if you have comments about the process of tracking edits / changes to Models. If you have specific expert advice about avoiding this specific error, or dealing with the error once it arises then I'd be very grateful.


title: EPSG code as dynamic text in ArcGIS

text: Is it possible to get ArcGIS to show the EPSG code for a coordinate system as dynamic text? The default text is quite verbose (see below), and I think just showing the "Coordinate System" is probably enough. However, it would be quite useful if I could also show "EPSG:32632" in addition to this.               Coordinate System: WGS 1984 UTM Zone 32N     Projection: Transverse Mercator     Datum: WGS 1984     False Easting: 500,000.0000     False Northing: 0.0000     Central Meridian: 9.0000     Scale Factor: 0.9996     Latitude Of Origin: 0.0000     Units: Meter      The documentation for ArcGIS 10.1 doesn't mention it, but I'm hoping it's just been missed out, or there is a work-around.


title: Postgis array query

text: I have a Postgis database with land cover vector data with the following layout:               select * from CLCtable limit 0;      outputs:               gid | code_06 | id | remark | area_ha | shape_leng | shape_area | the_geom      My problem is the following: I have a 2D point array and I want to get the 'code_06' value for each point. I'm currently doing a query inside a for loop that goes through all the points. I end up with a query for each point which takes too much time if I have for instance a 100x100 2D point array. The query I'm doing is something like this:               select a.code_06 from CLCtable as a      where ST_Contains(a.the_geom, ST_GeomFromText('POINT(1234 4567)', 3035));      Is there some way to this more efficiently?


title: How do I use a DHM in tiff format for distance and slope measurements?

text: I am a student and I have to work out a project about soil loss on a plot. I need basic measurements of this plot such as length, width, area and slope for my erosion calculations. I have got my hands on a 4 GB large DHM measured by LIDAR in a TIFF file format. My question is now, how do I preform these calculations on this file format? I guess I need to convert it so I can use it in IDRISI or ArcMap?


title: Prevent ArcGIS saving shapefiles into workspace

text: How can I stop ArcGIS creating output shapefiles as results of various tools? I have got a script returning a certain value. However in order to get this value a have to use several functions such as `arcpy.CreateThiessenPolygons_analysis` or `arcpy.Clip_analysis`. These functions (tools) return maps and save them into my workspace. I would like to work with them in the script but I dont want them to be saved into the workspace since I wont need them anymore. Thank you


title: Finding the latitude/longitude that is in the middle of 2 others?

text: Ok, so I will have a system/application that will always at the most have 2 points (latitude/longitude) on a map per street. So if the user wants to add another point on that street I want my app to calculate if the new point is outside the bounds of the previous two points. So basically it will be like a max min type situation, and I want to know if this new point is a new max or a new min or if it should be thrown out, because I only track the max and mins. What mathematical magic hoop am I going to have to jump through to figure that out? or should I go to math.stackexchange.com?


title: How to do parallel queries in postgis

text: Is there any way to do queries in postgis using multiple cores in a transparent way, other than splitting my queries among several client connections ?


title: Create a View with a new geometry field in PostGIS

text: I have a **table** in my data base. This table was created with differents fields as 'x','y' and 'z'. Now I need to make a **view** to generate a geometry column, but I found problems, because the created geometry field is not correctly named. I need to force my column "geom geometry" to "geom geometry(PointZ, 25830)". I did this: st_force_3dz(st_setsrid(st_makepoint(table.x, table.y, table.z), 25830)) Any solutions? Thank you.