title: Error when creating insert cursor in arcpy

text: I've written several scripts that copy data from MS Access db into a new, empty ESRI Geodatabase. I've pretty much copied the code and changed it for each table -- the first two scripts have worked as expected, but the third is **throwing an error on the line used to create the Insert Cursor.** All of the variables are storing correctly at the top of the script, but the Insert Cursor is never created. Error Message: > Runtime error Traceback (most recent call last): > > File "", line 49, in File "c:\program files > (x86)\arcgis\desktop10.1\arcpy\arcpy__init__.py", line 1005, in InsertCursor > return gp.insertCursor(dataset, spatial_reference) > > File "c:\program files > (x86)\arcgis\desktop10.1\arcpy\arcpy\geoprocessing_base.py", line 383, in > insertCursor self._gp.InsertCursor(*gp_fixargs(args, True))) > > RuntimeError: ERROR 999999: Error executing function. Code:               #Import modules, environment settings     import arcpy     arcpy.env.overwriteOutput = True          # Database connection variables     geodatabase = 'P:/GIS_Development/ArcGIS Database Design/Landowner.gdb/Parcels'     accessDatabase = 'Y:/Database/LANDOWNER/NewLandowner.mdb/TblParcels'          #Create cursor for ACCESS Database and New Landowner GDB     accessCursor = arcpy.SearchCursor(accessDatabase)          # Create a variable that stores the value for column in a given row in the access table          for accessRow in accessCursor:              parcelID   = accessRow.getValue('Parcel_ID')         flag       = accessRow.getValue('Flag')         agreeID    = accessRow.getValue('Agreement_ID')         agreeStat  = accessRow.getValue('AgreementStatus')         ownerType  = accessRow.getValue('Owner_Type')         add1       = accessRow.getValue('Parcel_Address1')         city       = accessRow.getValue('Parcel_City')         county     = accessRow.getValue('Parcel_County')         state      = accessRow.getValue('Parcel_State')         zipCode    = accessRow.getValue('Parcel_Zip')         acreage    = accessRow.getValue('Acreage')         landClass  = accessRow.getValue('Land_Class')         mapGrid    = accessRow.getValue('Map_Grid')         parcelNum  = accessRow.getValue('Parcel_Num')         parcelCom  = accessRow.getValue('Parcel_Comments')         enterRes   = accessRow.getValue('Entry_Restrict')         enterCom   = accessRow.getValue('Entry_Comments')         contacted  = accessRow.getValue('Contact_Made')         appStatus  = accessRow.getValue('ApprovalStatus')         owner      = accessRow.getValue('Owner')              # "Paste" the stored value from the access table into the new geodatabase         gdbCursor = arcpy.InsertCursor(geodatabase)         row = gdbCursor.newRow()              if state:             state = state.upper()              #Fill out fields with variables that are storing MS Access data         row.Parcel_ID            = parcelID         row.Flag                 = flag         row.Agreement_ID         = agreeID         row.Agreement_Status     = agreeStat         row.Owner_Type           = ownerType         row.Address_Line_1       = add1         row.City                 = city         row.County               = county         row.State                = state         row.Zip_Code             = zipCode         row.Acreage              = acreage         row.Land_Classification  = landClass         row.Map_Grid             = mapGrid         row.Parcel_Number        = parcelNum         row.Parcel_Comments      = parcelCom         row.Entry_Restrictions   = enterRes         row.Entry_Comments       = enterCom         row.Contact_Made         = contacted         row.ApprovalStatus       = appStatus              # Insert new row into Geodatabase         gdbCursor.insertRow(row)


title: What is the easiest way to calculate borehole deviation?

text: ![enter image description here](http://i.stack.imgur.com/Qneo3.jpg) Does anyone know if there was a program or specific method (calculations) in order to create a deviation value based on querying specific geological and drilling variables?


title: Geometry has incorrect latitude-longitudes

text: * `SELECT ST_SRID(geom) AS s FROM zones GROUP BY s;`   returns a single record: `4326`   * `SELECT Find_SRID('public', 'zones', 'geom');`   returns `4326`   * `SELECT AsText(geom) FROM zones`   returns geometry like: `MULTIPOLYGON(((495468.88595 5464400.7263,495468.89 5464401.38005,...)))` Shouldn't my geometry been standard lat-lngs? Something like:               MULTIPOLYGON(((-122.968967906831 49.31253667013,-122.970356138605 49.3125609975662...)))      I tried executing both of these queries, but neither worked: `SELECT UpdateGeometrySRID('zones','geom',4326);`   `UPDATE zones SET geom = ST_Transform(geom);` How can I transform the geometry into my desired format?


title: Need to learn more about Precision and Accuracy

text: I am currently taking Introduction to Image Interpretation online and there is a section in it that talks about Precision and Accuracy. I wanted to know if someone out there has any good information that I can use to learn more about it.


title: Using Spatial Analyst Sample tool

text: I have a bunch of rasters and a point file. I need the information from ll the raster files to be added to my points. My usual workflow is to use the Spatial Analyst sample tool. It works fine and I can do a spatial join to get the information added to the point file. However, furstratingly the sample tool never seems to offer useful field headings and with a lot of rasters it can become confusing as to which value relates to which raster. Does anyone know of a script or workaround that can produce a table with field headings?


title: ArcGIS forces me to close the program using select layer by attribute

text: I am using ArcGIS 10.2 At the moment, I am trying to select from a layer of different polygons specific polygons from a specific area but when I am using the tool select Layer by attribute, ArcGIS get stack and It forces me to close the program, what it could be the problem?


title: Could not commit changes to layer X

text: I am editing a layer's (Z) attribute table. I made some columns, used the field calculator to fill them up, made some more, went back and deleted some and then tried to save the edits and got the error message: "Could not commit changes to layer Z. Errors: ERROR: 3 attributes not deleted. Success 2 attributes added Success: 10 attribute values changed" Then I closed the error box and the layer remains in edit mode and one can't save it. I tried to save the Z layer and got the same "Could not ..." but also with "OGR error deleting field 3: Invalid field address" I saved the project and then reloaded the project to find that all the columns in layer Z had been wiped out except an ID layer. The layer Z is ruined. v 1.8.0. What is going on? - thanks.


title: How to update a SLD file in GeoServer?

text: I'm working on my own geoportal. I have prepared a php page which users can describe their own SLD configuration. The SLD file will be stored in GeoServer style folder. But when I tried to view maps I couldn't see the changes on map. But if I go to the styles editor and click to the Submit button in my new SLD files, I can see the new SLD configuration on my map. I think Submit button is making some another configuration that ı don't know. Are there any way refreshing the style edits?


title: Heatmap algorithm to visualise point diversity

text: Can anyone suggest an algorithm to generate a heatmap for visualising point diversity? An example application would be for mapping areas of high species diversity. For some species, every single plant has been mapped, resulting in a high point count, but with very little meaning in terms of the diversity of the area. Other areas genuinely have high diversity. Consider the following input data:               x    y      cat     0.8  8.1    B     1.1  8.9    A     1.6  7.7    C     2.2  8.2    D     7.5  0.9    A     7.5  1.2    A     8.1  1.5    A     8.7  0.3    A     1.9  2.1    B     4.5  7.0    C     3.8  4.0    D     6.6  4.8    A     6.2  2.4    B     2.2  9.1    B     1.7  4.7    C     7.5  7.3    D     9.2  1.2    A      and resulting map: ![grid plot](http://i.stack.imgur.com/8kp2h.png) In the upper left quadrant, there is a highly diverse patch, while in the lower right quadrant, there is an area with high point concentration, but low diversity. Two ways to visualise the diversity could be to use a traditional heatmap, or to count the number of categories represented in each polygon. As the following images show, these approaches have limited use, since the heatmap shows the greatest intensity in the lower right, while the binning approach would look exactly the same if there were only one category (this could be addressed by increasing the size of the polygon bins, but then the result becomes unnecessarily granular). ![heatmap](http://i.stack.imgur.com/xYFtG.png) ![enter image description here](http://i.stack.imgur.com/k3r9h.png) One approach I thought of to do this would be to prime a traditional heatmap algorithm by the number of points of different categories within a defined radius, and then using that count as the weight for the point when generating the heatmap. However, I think this might be prone to unwanted artefacts, like mutual reinforcement leading to very sharp results. Also, closely mapped points of the same type would continue to show up as high concentrations, just not to the same extent. Another approach (probably better but more computationally expensive) would be:   1. Calculate the total number of categories in the dataset   2. For each pixel in the output image:      * For each category:        * calculate the distance to the nearest representative point (r) [probably limiting by some radius beyond which influence is negligible]       * add a weighting proportional to 1/r2 Are there already algorithms that I'm not aware of to do this, or other ways to visualise diversity? **Edit** Following Tomislav Muic's suggestion, I have calculated the heatmaps for each category, and normalised them using the following formula (QGIS raster calculator):               ((heatmap_A@1 >= 1) + (heatmap_A@1 < 1) * heatmap_A@1) +     ((heatmap_B@1 >= 1) + (heatmap_B@1 < 1) * heatmap_B@1) +     ((heatmap_C@1 >= 1) + (heatmap_C@1 < 1) * heatmap_C@1) +     ((heatmap_D@1 >= 1) + (heatmap_D@1 < 1) * heatmap_D@1)      with the following result (comments under his answer): ![normalised sum](http://i.stack.imgur.com/iXlVH.png)


title: Arcobjects in Python, inherit from Abstract Class

text: As for question, I am trying to adjust Grid Map labels via Arcobjects with Python. For now I do have:               from comtypes.client import CreateObject, GetModule     import arcpy          def CType(obj, interface):        """Casts obj to interface and returns comtypes POINTER or None"""        try:            newobj = obj.QueryInterface(interface)            return newobj        except:            return None          def NewObj(MyClass, MyInterface):        """Creates a new comtypes POINTER object where\n\        MyClass is the class to be instantiated,\n\        MyInterface is the interface to be assigned"""        from comtypes.client import CreateObject        try:            ptr = CreateObject(MyClass, interface=MyInterface)            return ptr        except:            return None          esriCarto = GetModule(r"C:\Program Files (x86)\ArcGIS\Desktop10.0\com\esriCarto.olb")     esriCartoUI = GetModule(r"C:\Program Files (x86)\ArcGIS\Desktop10.0\com\esriCartoUI.olb")     esriMapUI = GetModule(r"C:\Program Files (x86)\ArcGIS\Desktop10.0\com\esriArcMapUI.olb")     esriFrame = GetModule(r"C:\Program Files (x86)\ArcGIS\Desktop10.0\com\esriFramework.olb")          arcpy.SetProduct('Arcinfo')          pApp = NewObj(esriFrame.AppROT, esriFrame.IAppROT).Item(0)     pDoc = pApp.Document     pMxDoc = CType(pDoc, esriMapUI.IMxDocument)     pLayout = pMxDoc.PageLayout     pGraphContLayout = CType(pLayout, esriCarto.IGraphicsContainer)     iFrameElement = pGraphContLayout.FindFrame(pMxDoc.ActiveView.FocusMap)      As far as I understand, iFrameElement is an interface of an abstract class from which I need to inherit attributes (pointer) to MapFrame object. How do I do that? How do it get to object with IMapGrids interface. Any suggestions?


title: Rotate Selected Features ArcGIS 10 - VBA

text: I am attempting to loop through polygon features and rotate pairs of features by a specified angle. Pseudo logic: 1) Open up a standalone table in the map document and retrieve a unique identifier that will used to select a pair of polygons from a feature class. This code works just fine. Code retrieves the identifying value for selecting the polygons and the rotation angle. 2) Create a selection using the unique identifier. In this application only two records will be returned. 3) Rotate the selected polygons around the SelectionAnchor 4) Advance to the next standalone table record, repeat making selection and rotation Every works fine up to the point that the rotation routine attempts to apply QuerryCoords on the selection anchor. The SelectionAnchor is seen as empty geometry. If I manually select the polygon set, and run the RotateSelectedFeatures routine, it produces the desired outcome. Code for selection features:               Sub RotateSelectedFeatures()         Dim sHRMashUP As String         Dim intLayerCnt As Integer         Dim pMXDoc As IMxDocument         Dim pMap As IMap         Dim pFeatureLayer As IFeatureLayer         Dim pDataset As IDataset         Dim pUID As New UID         Dim pEditor As IEditor         Dim pRotatePoint As IPoint         Dim pEnumFeature As IEnumFeature         Dim pFeature As IFeature         Dim pRotateSet As ISet         Dim pFeatureEdit As IFeatureEdit         Dim pFeatureSelection As IFeatureSelection         Dim pQueryFilter As IQueryFilter         Dim sSQL As String         Dim pActiveView As IActiveView         Dim pInvalidArea As IInvalidArea         Dim pStartPoint As IPoint         Dim origX As Double         Dim origY As Double              Set pMXDoc = ThisDocument         Set pMap = pMXDoc.FocusMap         Set pActiveView = pMap              pUID.Value = "esricore.Editor"              Set pEditor = Application.FindExtensionByCLSID(pUID)              'Name of layer with home range polygons         sHRMashUP = "HomeRangePairs_Mashup_2"              'Loop through the table of contents to see which layer is the target home range poygons identified above         For intLayerCnt = 0 To pMXDoc.FocusMap.LayerCount - 1             If pMXDoc.FocusMap.Layer(intLayerCnt).Name = sHRMashUP Then                 Set pFeatureLayer = pMXDoc.FocusMap.Layer(intLayerCnt)             End If         Next intLayerCnt                   If pEditor.EditState = esriStateEditing Then             pEditor.StopEditing False         End If              If pEditor.EditState = esriStateNotEditing Then             Set pDataset = pFeatureLayer.FeatureClass             pEditor.StartEditing pDataset.Workspace                  ' Select features using the X,Y coordinates read in from text file             Set pFeatureSelection = pFeatureLayer 'QI                  'Create the query filter for selecting on the field with the coma delimited xy location read from the text file             Set pQueryFilter = New QueryFilter                  sSQL = "XYCoords = " & Chr(39) & g_sXYCoords & Chr(39)                  pQueryFilter.WhereClause = sSQL                  ' Invalidate only the selection cache             'Flag the original selection             pActiveView.PartialRefresh esriViewGeoSelection, Nothing, Nothing                  'Perform the selection             pFeatureSelection.SelectFeatures pQueryFilter, esriSelectionResultNew, False                  'Make sure something has been selected             'Debug.Print g_sXYCoords             If pEditor.SelectionCount = 0 Then                 pEditor.StopEditing False                 Exit Sub             Else                 Call RotateSelectedFeatures             End If                  'Flag the new selection             pActiveView.PartialRefresh esriViewGeoSelection, Nothing, Nothing         End If     End Sub     Public Sub RotateSelectedFeatures()       Dim pEditor As IEditor       Dim pID As New UID       Dim pEnumFeature As IEnumFeature       Dim pFeature As IFeature       Dim mySet As ISet       Dim pSpatialReference As ISpatialReference       Dim pStartPoint As IPoint       Dim pEndPoint As IPoint       Dim pFeatureEdit As IFeatureEdit       Dim pLine As ILine       Dim origX As Double       Dim origY As Double       Dim Count As Integer            'Get a reference to the editor extension       pID = "esriEditor.Editor"       Set pEditor = Application.FindExtensionByCLSID(pID)       'Make sure something has been selected       If pEditor.SelectionCount = 0 Then         MsgBox "Nothing Selected"         Exit Sub       End If       'Create an edit operation enabling undo/redo       pEditor.StopOperation "Rotate Selection"       pEditor.StartOperation            'Get the selection from the Editor and add all the features to a new Set       Set pEnumFeature = pEditor.EditSelection       Set mySet = New esriSystem.Set       pEnumFeature.Reset       For Count = 0 To pEditor.SelectionCount - 1         Set pFeature = pEnumFeature.Next         mySet.Add pFeature       Next Count            'Reset the Set       mySet.Reset            'MoveSet requires a line to specify the new location       'Use the selection anchor as a starting point for the line       Set pStartPoint = pEditor.SelectionAnchor.Point       pStartPoint.QueryCoords origX, origY            ' Rotate the set       Do         Set pFeatureEdit = mySet.Next         If Not pFeatureEdit Is Nothing Then           pFeatureEdit.RotateSet mySet, pStartPoint, g_dblAngle         Else           Exit Do         End If       Loop            'Stop the Edit Operation       pEditor.StopOperation "Rotate Selection"  'String to appear in Edit menu            'Refresh the Map to redisplay the moved features       pEditor.Display.Invalidate Nothing, True, esriAllScreenCaches          End Sub      The Where clause uses two variables set in the calling routine If a selection is produced, this routine calls the RotateSelectedFeatures routine which acts on the selected set. Bails on the pStartPoint.QueryCoords origX, origY line due to empty geometry.


title: Is there a way to copy the mxd file with all of the files in that in ArcGIS?

text: Obviously the mxd doesn't contain the gdb files, and the linked raster files. So when I copy a project, I need to search manually the 49 raster maps from the other 4000. Is there an easier way to copy this?


title: Help with ListFields and CalculateField_management syntax

text: I need some assistance with the syntax for a CalculateField_management calculation as I have still yet to master all the python syntax rules. What I am trying to do is find the max value in a variable number of fields in order to populate another field. I am using the ListFields function to discover the desired fields to choose from, but getting that list into the formula is giving me some difficulty.               import arcpy, os, string     LAYERS = arcpy.GetParameterAsText(0)     SLOSHFILEDS = [f.name for f in arcpy.ListFields(LAYERS,"","DOUBLE")     arcpy.CalculateField_management (LAYERS, "MAXSURGE", max(SLOSHFILEDS))      I have tried any number of different string combinations for the max() calc to no avail (not that this particular variation shows that). Adding/Changing the following to the script doesn't give me the syntax error that I would recieve with the above, but it does give me a "The calculate value is invalid for the row with ObjectID = 0..." x 18,526 (or however many rows are in my table) and then does nothing to the table except populate my MAXSURGE field with 0's.               SLOSHFILEDS = arcpy.ListFields(LAYERS,"","DOUBLE")     fieldNameList = []     for field in SLOSHFILEDS:             if not field.required:                 fieldNameList.append(field.name)     arcpy.CalculateField_management (LAYERS, "MAXSURGE", max(fieldNameList))      Hard coding the field names into the formula works great, but of course, I will not always have the same number of fields or same field names to work with. Any help would be appreciated. -John


title: locator service to use with ArcGIS Javascript API v2.1

text: I am hoping to find a locator service that will locate by address or intersection, and that will work with version 2.1 of the ArcGIS Javascript API. I tried this locator given in the examples, which works well with 3.1, but it doesn't seem to work with 2.1 : http://help.arcgis.com/en/webapi/jav...r_address.html Unfortunately I think I need to stick with v2.1. SO far I haven't had any luck. Thanks in advance, Jason


title: How to Build a graph from a map

text: I spent days looking in all questions here and trying tutorials. I need to build a graph for a region (It will be my input in simulation) I downloaded Postgresql and Postgis and osm2pgsql. Now i can load .osm files into a database. I want to have a complete database with roads ids so i can generate a tbale of edges and nodes then build my graph nodes = intersections and edges = roads for a city using python. Please any detailed help with steps.I saw a lot of keywords (Postgresql/postgis, Networkx, QGIS ..) but i couldn't manage to assemble a good tutorial I am a very beginner.


title: Why the result of a merge of multiple raster is so big?

text: I try to merge 14 geotiff like this : ![enter image description here](http://i.stack.imgur.com/0UsMo.jpg) Each geotiff is about 50Mb. I need a geotiff at the output My workflow :               gdalbuildvrt -input_file_list list.txt test.vrt       (where my list contains the name of the tifs) Then :               gdal_translate -of Gtiff test.vrt test.tif     Input file size is 79841, 59955      It works, but the result is a geotiff of 13,3 Gb ! For 14 files, each 50 Mb, I attempted a geotiff of 700 Mb, not 13 Gb. I know that gdal does not compress by default, so I tried this command :               gdal_translate -of Gtiff -co COMPRESS=JPEG test.vrt test_compressed.tif      But the "merge" of the file is too big for JPEG compression :               Input file size is 79841, 59955     0ERROR 1: JPEGPreEncode:Strip/tile too large for JPEG     ERROR 1: WriteEncodedTile/Strip() failed.     ERROR 1: JPEGPreEncode:Strip/tile too large for JPEG     ERROR 1: WriteEncodedTile/Strip() failed.     ERROR 1: An error occured while writing a dirty block     ...      So I tried another workflow, and converted all my tifs in jpeg (14 Mb each), build a vrt file and translated it with LZW compression. But the output geotiff is about 5 Gb. Could you tell me what is the best practice to do the job, and if it is possible to obtain one geotiff of 14*50Mb ? I didn't try it, but I thought about merge this tifs in photoshop, then re- georeference with upper left/lower right coordinates. With this workflow I think I will have 14*50 Mb, but I'm not sure. And I want to learn gdal best practises so I didn't try it for the moment Sorry for my english Thank you EDIT : > _coming to bites: if input is tif with 8bit and export is 32bit by default > you will get serious trouble. so make sure you keep your byte definition as > it is. And remember: the full tif will prob. have 20x 50mb as a tiff is > always rectangular_ If I understand, the number I pointed in green on this screenshot have to be the same on the left and on the right ? ![bits](http://i.stack.imgur.com/Rx52O.jpg) > _Your output image will have more pixels than the sum of your input images, > but this does not explain the large difference. I suggest that you look at > the characteristics of your images based on gdalinfo in order to see what > compression is used and check that the extents are correct._ The 14 tifs of 50 Mb were originally 14 tifs of 700 Mb that I processed with gdal_translate with -co COMPRESS=JPEG. I compressed the raster in order to reduce the number of Mb, but maybe it was not a good idea ? This screenshot represents the 2 gdal info of the same geotiff ( 01.tif) , at the left of the screenshot is is the gdalinfo of the non-compressed Gtiff of 700 Mb, end a the right the same Gtiff with COMPRESS=JPEG, so 50 Mb, with the diff in green: ![non-compress and compress gdalinfo of the 01.tif file](http://i.stack.imgur.com/TusKn.jpg) According to me, the extents are correct because in qgis it matches with other source of data and satellite imagery. > *assuming your input images have the same size, it makes 20000 * 12000 > pixels per input images, which is large for an image of 50 Mb, maybe you are > crossing the extent of your coordinate system when you create the mosaic.* I'm not sure to understand what you mean by "crossing the extent". But I tried to open my 5 Gb LZW in QGIS, and the extent is good, because it matches with other source of data. Your answer makes me realize that the Gtiff have not the same size, do you think it could be the cause of the increasing of the size when merged ? Because gdal prefers file of the same size. I made a gdalinfo on each Gtiff to get its size, there is a very small difference between the size of Gtiffs :               02.tif Size is 19956, 11981     03.tif Size is 19959, 11993     04.tif Size is 19961, 11992     05.tif Size is 19958, 11993     06.tif Size is 19958, 11990     07.tif Size is 19956, 11984     08.tif Size is 19956, 11993     09.tif Size is 19958, 11993     10.tif Size is 19958, 11989     11.tif Size is 19958, 11985     12.tif Size is 19958, 11993     13.tif Size is 19959, 11993     14.tif Size is 19960, 11994      > _Then you should look at the pixel depth of your images : if your input were > in Bytes, >then you should keep bytes. gdal_translate -of Gtiff -ot Byte -co > COMPRESS=LZW test.vrt test.tif_ I tried this command but the gdal told me that the size of the tiff was exceeded.               Input file size is 79841, 59955     0...10...20...30...40...50..ERROR 1: TIFFAppendToStrip:Maximum TIFF file size ex     ceeded. Use BIGTIFF=YES creation option.     ERROR 1: WriteEncodedTile/Strip() failed.      But if I have to create a big tiff, it does not resolve my problem because it is more than 4 Gb. Is the pixel depth important in my case ? ( HD photo of maps, then georeferenced, not DEM) > _Remark 1 : Converting your images to jpeg before building a vrt doesn't > help and you might loose data._ It is not serious if I loose a little bit information. I prefere not to loose it, of course, but if I have to it's not a problem. I was convinced that the output would be lighter if I worked with jpeg, but as a conclusion, it's not true when the output is Gtiff. So this is not a good solution. I give up this solution. > _> Remark 2 : using a vrt is helpful: are you sure that you need GTiff ?_ Yes, I need a Gtiff, because I have to import it in a mobile application which need geotiff input to work ( I think the app can take geospatial pdf input too, but I never work with it, and I want to understand my problem with gdal, because it is not the first time I have it). EDIT 2 : Sorry for the late. I tried -co tiled=yes -co bigtiff=yes -co compress=jpeg -co photometric=ycbcr and i tried -co TILED=yes -co BLOCKXSIZE=512 -co BLOCKYSIZE=512 These 2 commands work well, I have a size of ~700 Mb. It is exactly what I expected. Thank you very much Now I have another problem : it can't be opened quickly by QGIS. I have to wait 15 minutes ( but I quit before QGIS open the tif successfully). I don't know why. And in my android app, it does not work (maybe cause of "tiled=yes"). I have to read some doc on my own, thank you very much for your help


title: DeleteField_management gives ERROR 000308 Invalid field type?

text: I am attempting to delete a list of fields in a feature layer but am getting an error: Execute error: Failed to execute. Parameters are not valid. ERROR 000308 Invalid Field type Here is a code snippet:               arcpy.MakeFeatureLayer_management(roads, tmprds)          arcpy.DeleteField_management(tmprds, ["FIELD_NAME"])      I have listed the fields in my feature layer using:               fclass = tmprds     fieldnames = [f.name for f in arcpy.ListFields(fclass)]     print (fieldnames)      and I get a list of layers formatted like:               [u'FIELD1', u'FIELD2, u'FIELD3']      I have tried using both IN_MEMORY and real locations which don't seem to make a difference. Where am I going wrong? Thanks


title: How to sort Google DynamicMapsEngineLayers?

text: I tried to add several Google Maps Engine layers using the following code:               gmeLayer = new google.maps.visualization.DynamicMapsEngineLayer({         accessToken: token,         layerId: layerId     });     gmeLayer.setMap(map);      I would like to display the ploygon layer on the bottom and point layer on the top. How can I set the order for each layer? Any help is appreciated!


title: How do you display each polygon segment while digitizing a polygon with the Polygon construction tool?

text: I am almost positive that I've digitized the polygons and that I have seen each segment shown as I dropped in vertices during the sketch. I figure that I inadvertently configured something and now I can't get it back. Once I finish the sketch (by double-click or F2), the polygon redraws and displays as expected. The symbol has an outline and there is only one feature template. Here's the polygon, mid-sketch: ![enter image description here](http://i.stack.imgur.com/ejOqd.png) I'm expecting lines for the first 2 segments, but not getting them. Here's the polygon after I've finished the sketch: ![enter image description here](http://i.stack.imgur.com/x3rxZ.png) Here is my configuration while digitizing the polygon: ![enter image description here](http://i.stack.imgur.com/qxD7W.png) And this is the construction tool I'm using:   ![enter image description here](http://i.stack.imgur.com/gUnWs.png)


title: Compute the latitude of a center of mercator projection

text: Say I'm looking at a portion of a map that's displayed using Mercator projection with the view bounded at 40°N and 20°N. How do I compute the latitude of the point that's in the centre of that view? I figured it's not 30° and can compute it by firstly projecting the bounding latitudes onto a square using Mercator, computing the middle point, and then projecting back onto a sphere. But is there a simpler way of achieving the same without having to project to and from Mercator?


title: How can envelope angle be calculated?

text: Azure SQL has a hemisphere limitation where it errors on any geography having an envelope angle greater than 90 degrees. My coworker wants to set up some detection in the web server to catch such instances before they hit the database. How can envelope angle be calculated from a series of coordinates that define a polygon or multipolygon?


title: Levelplot(RasterVis) two marginal plot scale?

text: ![enter image description here](http://i.stack.imgur.com/aWmZ2.png) How can I add/display scale to two marginal plot?               RAD1998.all <- stack(list.ras)          list.files <- mixedsort(list.files(paste(getwd(), "/1998bil/", sep = ""), full.names = F, pattern = ".asc.bil"))      #list.files          outlist <- substr(list.files, 31, 43)     #outlist          png(paste("1998",outlist[[1]],".png",sep=""))     levelplot(RAD1998.all[[1]])     dev.off()


title: Conversion table convert XY coordinates to GPS

text: I have converted INline and Xline from Seismic lines to XY coordinates: How do I get to GPS? I am looking for a simple conversion table for field crews


title: Excellent app for OSM browsing on smartphone

text: I am looking for a decent app to browse OpenStreetMap on my Android smartphone in the offline mode. I have already bought OsmAnd+ and the Skobbler's OSM Navigation and Maps, and while their UIs are exactly at the opposite ends of the scale, neither offers some simple and important functionality. OSM Navigation and Maps interface is beautiful (a map is a visual thing, so eye-pleasing is paramount for efficient use) and simple in comparison to what OsmAnd+ does. It makes the maps fast to render. Navigation is good. What is lacking dramatically is all the additional data present in OSM (POIs). OsmAnd+ knows a lot more, but creates screens like this: ![enter image description here](http://i.stack.imgur.com/q8htd.png) Each tiny street is drawn with two lines, and all remote shops and train stations are marked. At the same time, getting the information from the program is almost impossible. The street names are shown randomly. I cannot click a street and get information about it. I cannot ask the program to hide most of the things and keep just what I want. So my question is: is an app I am searching for available? I want to be able to   * see only relevant information (aka layers)   * click an object to see its data


title: Setting the correct zoom level from Sector bounds for flat map using WorldWind SDK

text: This is not coming out correct for me. Setting the zoom level based on some other code I got on here using the EarthFlat model, which does not set the correct zoom for the placemarks within a sector bounds.               private void invoke() {         SwingUtilities.invokeLater(new Runnable() {                  public void run() {                      ArrayList<Position> pos = new ArrayList<Position>();                 pos.add(Position.fromDegrees(32, -117));                 pos.add(Position.fromDegrees(48, -72));                 pos.add(Position.fromDegrees(40, -70));                      final Sector bounds = Sector.boundingSector(pos);                      setWwPanel(new WWPanel());                 final WWPanel wwPanel = getWwPanel();                 wwPanel.setPreferredSize(new Dimension((int) getSize().getWidth(), (int) getSize().getHeight()));                 wwPanel.getWwd().addGLEventListener(new GLEventListener() {                          public void display(GLAutoDrawable e) {                         //System.out.println("Map displayed");                     }                          public void dispose(GLAutoDrawable e) {                         System.out.println("Map disposed");                     }                          public void init(GLAutoDrawable e) {                         System.out.println("Map GL initialized");                     }                          public void reshape(GLAutoDrawable e, int arg1,                             int arg2, int arg3, int arg4) {                         System.out.println("Map reshaped");                     }                      });                      LayerList layerList = new LayerList();                 BMNGWMSLayer bmng = new BMNGWMSLayer();                 layerList.add(bmng);                 ViewControlsLayer vcl = new ViewControlsLayer();                 vcl.setShowPitchControls(false);                 vcl.setShowVeControls(false);                 vcl.setShowHeadingControls(false);                 layerList.add(vcl);                 Model flatModel = new BasicModel(new EarthFlat(), layerList);                 ((EarthFlat) flatModel.getGlobe()).setProjection(FlatGlobe.PROJECTION_MERCATOR);                 wwPanel.getWwd().setModel(flatModel);                 wwPanel.getWwd().addSelectListener(new ViewControlsSelectListener(wwPanel.getWwd(), vcl));                      RenderableLayer rl = new RenderableLayer();                 for(Position curPos : pos) {                     PointPlacemark pm = new PointPlacemark(curPos);                     pm.setAltitudeMode(WorldWind.RELATIVE_TO_GROUND);                     pm.setLabelText("Test");                     PointPlacemarkAttributes pma = new PointPlacemarkAttributes();                          pma.setImageColor(Color.GREEN);                     pma.setLabelFont(Font.decode("Verdana-Bold-11"));                     pma.setLabelMaterial(Material.GREEN);                     pm.setAttributes(pma);                     pm.setValue(AVKey.DISPLAY_NAME, "Text displayed on mouse over");                          rl.addRenderable(pm);                 }                      layerList.add(rl);                      BasicOrbitView fov = (BasicOrbitView) wwPanel.getWwd().getView();                 fov.setCenterPosition(Position.fromDegrees(bounds.getCentroid().getLatitude().degrees, bounds.getCentroid().getLongitude().degrees));          // this is what i am using to set zoom ??????                 double radius = wwPanel.getWwd().getModel().getExtent().getRadius();                 double viewRange = Math.sqrt(Math.pow(bounds.getMaxLatitude().degrees-bounds.getMinLatitude().degrees, 2)+Math.pow(bounds.getMaxLongitude().degrees-bounds.getMinLongitude().degrees, 2));                 double altitude = radius * Math.sin(Math.toRadians(viewRange * 0.5));                      fov.setZoom(altitude);                      add(getWwPanel(), BorderLayout.CENTER);             }              });     }      It sets the zoom fine, but the placemarks are off the view bounds. Expected result: ![enter image description here](http://i.stack.imgur.com/12iRr.png) Undesired/actual result: ![enter image description here](http://i.stack.imgur.com/GNr12.png)


title: Deleting vector lines beyond image extents using QGIS?

text: I am using QGIS 2.2 I have one vector line layer (roads) that stretches beyond my base image (see below). ![enter image description here](http://i.stack.imgur.com/dNUkM.jpg) I wish to have the roads "clipped/trimmed or deleted" right at the edge of the image's extent (see below): ![enter image description here](http://i.stack.imgur.com/u3PiE.jpg) How can I achieve this?


title: LIDAR Dataset displaying 500,000 meters to the East of correct position

text: First post on here.. I've been hung up for a couple days. Why does my LIDAR dataset show up ~500,000meters to the east of the correct position even though a spatial reference/coordinate system is provided? Some background:   1. The LIDAR data were downloaded from the USGS website as .zip files   2. Unzipped, the .las files have an associated .xml file with correct coordinates for tile position (UL Corner, UR Corner etc.)    3. Projected coordinate system (State Plane NAD 1983 Florida East Meters) is listed in .xml   4. The FL County layer in the image below is in the correct location. How can I permanently bind the projection/spatial information to the .las file (not need the separate manually created .prj file)? I'm trying to use the LAStools (from lastools.org). When I use the 'LASclip', the clipping polygon is in the correct location but arcmap thinks the .las file is 500,000 meters east still, so nothing (or everything) gets clipped. Thanks in advance! ![text in .xml](http://i.stack.imgur.com/Ql12Z.png) ![enter image description here](http://i.stack.imgur.com/NRXjl.png)


title: How to find nearby objects using Easting and Northing range

text: I have a bunch of objects which have associated Eastings and Northings in format OSGB36 for example:               265397  869562     547297  179002      If I have a starting Easting and Northing, and wanted to find related objects within X km is it possible to compute a Start/End Northing and Start/End Easting range that I can search for programatically? I am using C#. I.e. searching for objects using pseudo-code :               Easting >= 265390 && Easting <= 265400 &&     Northing >= 869560 && Northing <= 869564      would this find nearby objects? What resolution do these coordinates operate to? I can operate in miles or km, just need to know the scale - if such a thing is possible. Sorry this is new to me.


title: Hiding selected features in ArcGIS for Desktop using ArcPy or ModelBuilder?

text: I'm using data driven pages to create a series of maps of point features. Some of my features overlap spatially and I only want to show some of them on each map. Once I select the features which I want to show, is there any way that I can hide the rest using ArcPy or ModelBuilder? My whole setup is automatic, I don't want to do this manually!


title: CartoDB get rows in initialized map

text: Is possible to get the data in this example? I mean something like:               map = cartodb.createVis ...     map.data.each ...      ? Or, where in documentation can I find something about it? or, which keyword should I look for?


title: ArcPy - Add Layer from Cursor Selection

text: I am trying to add layers from the fields of a reference layer using select by location and a cursor. Also, want to prompt the user with a message box asking if they want to add the selected by location. I am able to print string value from the reference layer when the user responds with "Yes" with the message box, but get an error when passing the string value as a path using "arcpy.mapping.addLayer".               def onMouseDownMap(self, x, y, button, shift):              mxd = arcpy.mapping.MapDocument("CURRENT")         df = arcpy.mapping.ListDataFrames(mxd,"*")[0]         #layer path         lyr1 = arcpy.mapping.Layer("C:\ArcGis\Aerial_Imagery.lyr")         #Add Layer         arcpy.mapping.AddLayer(df, lyr1, "TOP")                   point = arcpy.Point(x, y)         ptGeometry = arcpy.PointGeometry(point)         #Select layer by points/mouse down on map         arcpy.SelectLayerByLocation_management("Aerial_Imagery.lyr", "INTERSECT", ptGeometry)         #cursor to display selected layers on field name "Layer_Path"         cursor = arcpy.SearchCursor("Aerial_Imagery.lyr")         field = "Layer_Path"         for row in cursor:             selpath = (row.getValue(field))             message = "Would you like to add this Imagery:" + selpath             Addbutton = pythonaddins.MessageBox(message, "Imagery Finder", 4)             if Addbutton == "Yes":                 addlyr = arcpy.mapping.Layer(selpath)                 arcpy.mapping.AddLayer(df, addlyr, "TOP")#this breaks                      print selpath #this works                  else:                 pythonaddins.MessageBox("this is the answer to NO", "no answer", 0)             #Remove Layer         for lyr in arcpy.mapping.ListLayers(mxd, "Aerial_Imagery.lyr"):             #print "Remove Layer", lyr             arcpy.mapping.RemoveLayer(df, lyr)              #Refresh view and TOC         arcpy.RefreshActiveView()         arcpy.RefreshTOC()


title: refactoring a osmfilter to do a better work

text: hello dear gis experts new to the gis-technique: i am workin with the following tools: this               http://wiki.openstreetmap.org/wiki/DE:Osmconvert      and that               http://wiki.openstreetmap.org/wiki/DE:Osmfilter      i get the following example: wget download.geofabrik.de/openstreetmap/europe/germany.osm.pbf osmconvert germany.osm.pbf --all-to-nodes -o=germany.o5m osmfilter germany.o5m --keep="amenity=shop or shop=*" -o=shops.o5m osmconvert shops.o5m --csv="@id @lon @lat shop name" --csv-headline -o=shops.csv that runs fine - now i want to refine: to get such a result- [made out of xslt]               <node id="52768810" lat="48.2044749" lon="11.3249434" version="7" changeset="9490517" user="wheelmap_visitor" uid="290680" timestamp="2013-10-07T20:24:46Z">     <tag k="addr:city" v="Olching" />     <tag k="addr:country" v="DE" />     <tag k="addr:housenumber" v="72" />     <tag k="addr:postcode" v="82140" />     <tag k="addr:street" v="Hauptstraße" />     <tag k="amenity" v="restaurant" />     <tag k="cuisine" v="mexican" />     <tag k="email" v="info@cantina-olching.de" />     <tag k="name" v="La Cantina" />     <tag k="opening_hours" v="Mo-Su 17:00-01:00" />     <tag k="phone" v="+49 (8142) 444393" />     <tag k="website" v="http://www.cantina-olching.com/" />     <tag k="wheelchair" v="no" />      doable!?


title: ERDAS Imagine ecw compression settings

text: We have some jp2 or img satellite images that we compress to .ecw using ERDAS Imagine. We have realized that with some images we are losing too much information while doing that. What i wanted to ask is about preferences on ERDAS Imagine. Target ECW Compression (Greyscale): Value Target ECW Compression (RGB): Value Target ECW Compression (Multispectral): Value ![enter image description here](http://i.stack.imgur.com/ziVSR.png) Do I need to adjust them each time I compress an image or is there an optimal setting that i can use? (We are using mostly 16 bit satellite images as input) Or where else could the problem be that the ecw files are so foggy and unsharp compared to the input files. img: ![enter image description here](http://i.stack.imgur.com/4Y6IV.jpg) ecw: ![enter image description here](http://i.stack.imgur.com/FnrXf.jpg) Thanks in advance, Yasmin


title: Reprojecting vector layer in OpenLayers

text:
I have a simple OpenLayers map, where I would like to display the OSM map (the base map, in EPSG:900913) with some points on top of it (a vector layer, in EPSG:3763). I use proj4js. I do convert a point from EPSG:3763 to EPSG:900913, to set the map center. I also added the MousePosition control, and I can see the coordinates in the displayProjection system. so, the conversion runs fine between the two EPSGs. I created a new vector layer, with the EPSG:3763 projection. I add a point to that layer pragmatically, in the layer projection EPSG:3763. The point in the layer goes to Africa, into the ocean, around 0,0. My questions is: Does OpenLayers can reproject the entire layer? If it does, what am I doing wrong? From the documentation, it not not clear if OL does the reprojection or not of the layer. In one example, the reprojection is done by MapServer. On the other, the reprojection is done when reading the geojson. An alternative solution, is to convert the point coordinates to EPSG:900913 before adding them to the vector layer. But I really would like to see how to reproject an entire vector layer, in a different coordinate system. My code is:                                     <title>Multiple projections</title>             <meta http-equiv="Content-Type" content="text/html; charset=UTF-8" />                                       <link rel="stylesheet" type="text/css" href="http://dev.openlayers.org/releases/OpenLayers-2.13.1/theme/default/style.css" />

                    Happy mapping!


title: Calculate distance and surface with mapbox api

text: Is there a way to draw lines and polygons (with geojson) and calculate distance and polygon surface with the mapbox api or leaflet? I wish to have the same tools as in geojson.io with the extra ability to calculate line length and polygon surface. For lines I will use the Leaflet method _distanceTo_ but for the polygon surface I have no clue.


title: Geoserver layer rotation problem

text: I have a mapping system setup with PostGIS, Geoserver 2.4.1 and OpenLayers 2.12 I am trying to use the geoserver vendor specific parameter "angle" to rotate a layer, and it seems to work fine. However, If I zoom in and then pan the map, after I depress the mouse button the layer suddenly moves to the opposite direction of the pan movement. The amount of pixels the layer moves is proportional to the angle I rotated it, the higher the angle, the larger the unwanted layer movement. Below is the setting up of the layer. Anyone had this problem before? Thanks, Luis Lima                  // setup single tiled layer        untiled = new OpenLayers.Layer.WMS(                  "project:vw_1_shapes_hole_2 - Untiled", "http://192.168.2.60:8080/geoserver/ezcourse/wms",                  {                      LAYERS: 'ezcourse:vw_1_shapes_hole_2',                      STYLES: '',                  angle: '120',                        format: format                  },                  {                      singleTile: true,                        isBaseLayer: true,                                      yx : {'EPSG:4326' : true}                  }         );


title: Convert Sql Server Spatial data to GeoJSON & plot on map

text: I have a sql server table with Geo spatial data. An example of row would be (with Column names):               type: streetline     code: 231001     geog: 0xE6100000011 ........     Centroid 0xE61000000C.......     geom: 0xE6100000011 ........      I am using C# with MVC3 to draw the above shapes on leaflet map. I am fetching the above sql data into a datatable. Using the sql function "ToString()" on column "geom" I get the following:               "LINESTRING (-1.131510412 52.65531, -1.13286 52.65559)",     "POLYGON ((-1.1116360 52.6409953, -1.1116683 52.6413, -1.11146723 52.641317, -1.11133263 52.6413572, -1.1113059))",      The question is how do I convert the above to GeoJSON to be able to plot on the leaflet map. Here is an example of the expected output (GeoJSON) I want:               var geojsonFeature = {         "type": "Feature",         "properties": {             "name": "Coors Field",             "amenity": "Baseball Stadium",             "popupContent": "This is where the Rockies play!"         },         "geometry": {             "type": "Point",             "coordinates": [-104.99404, 39.75621]         }     };


title: Is there a way to remove this thick gray border from my raster symbology in the legend?

text: I have this problem with different raster layers, in this instance it's land use data. In my layout, and in a rendered PDF of my map, the legend entry for the raster layer has an ugly, thick gray border around my raster color patches. Here's how they look when rendered to PDF: ![enter image description here](http://i.stack.imgur.com/I0Wok.png) Sometimes it's even worse, and ArcGIS appears to randomly choose a thick border for some colors and a thinner, lighter one for others, making my legend very ugly. Here's another example in a finished PDF: ![enter image description here](http://i.stack.imgur.com/3uMsV.png) Here's what the same legend looks like in the ArcGIS layout. The border width looks the same, nothing like the PDF. The color differences are not as pronounced as in the PDF but you can see they're there: ![enter image description here](http://i.imgur.com/8vwZIGz.png) Is there any way to edit or remove those gray borders? I've searched a bit online and the most recent thread I've found is this one from the ArcGIS forums, with a couple workarounds that I may resort to, but I'm hoping for a simpler solution. FYI, my rasters' symbology is set to use Unique Values, but using Colormap gives the same results.


title: Strategies for managing (and ultimately merging) duplicate business locations in a spatial database

text: I am working on a project right now that involves using multiple (very large) "location directories" to identify with the highest possible inclusiveness, a diverse set of business locations. These locations are often geotagged close to one another but typically not directly on top of one another. For example, if a walmart was in one database at the NW corner of its parking lot, in another database it might be in the SE corner. In this example the two points might be as far as 2000ft apart. I'm using the extreme examples to create some basic buffer rules (after consolidating brandnames to identify same-brand store locations) for allowing my program to identify duplicate locations. But there are also some business locations that have closed. Some are geocoded directly on top of the "center of town" (i.e. they had an address of "Cityname, State" and got geocoded incorrectly in the middle of town). There are a lot of other edge cases and I guess what I'm asking is, is there some python code package or qgis plugin or ruby gem (haha, just kidding) or some other larger abstract object oriented design pattern (way of thinking) that somebody learned as a Computer Science major to handle (entity validation & consolidation) problems like this? I'm perfectly will and able to cook up my own solutions but don't want to miss out on sophisticated techniques and/or waste my time reinventing potentially less-effective, more time-consuming solutions for managing entity databases. **Am perfectly willing to modify, close this question or migrate it to another forum if a mod or mods thinks this something is not appropriate here. Please comment before taking any administration or negative voting action.** The technology I'm using currently is a ruby on rails database running a postgis database. I extensively utilize activerecord to manage my data. I also use qgis plugins a lot and am building a growing capability set with python though still very much a beginner.


title: Conditional SQL Query

text: I've been running a code that searches a desired field in a feature class and performs an specified SQL query based on whether the field is string or double. The variables are just temporary fillers to test the code, I'll be shifting them to GetParameterasText() as soon as I know that all of the code works. I need the code to work so that it can be turned into a ArcGIS toolbox script.               import arcpy     from arcpy import env     env.overwriteOutput = True          #Create working variables     fc = "C:/gisclass/assignment5/Results/CITIES.shp" #vector feature class     affectedField = "CAPITAL"     oldValue = -1     newValue = 2     integer = 0          fieldlist = arcpy.ListFields(fc)               for field in fieldlist:         if field.type == "String":             queryString = """"%s" LIKE '%s'"""%(affectedField, oldValue) or """%s" = '%s'"""%(affectedField, oldValue)         elif field.type == "Double":             queryString = """"%d" >= '%d'"""%(oldValue, integer) or """"%d" <= '%d'"""%(oldValue, integer)          with arcpy.da.UpdateCursor(fc, (affectedField,), queryString) as cursor:         for row in cursor:             row[0] = newValue             cursor.updateRow(row)          mess = arcpy.GetMessageCount()     print arcpy.GetMessage(0)     print arcpy.GetMessage(mess-1)      I have gotten the string where clause to work, but the integer clause is still producing the same result - Invalid SQL query. Any idea on what the solution to the problem is?


title: How do you insert a cartoDB map using Cartodb.js into a wordpress site?

text:
I've been looking at http://developers.cartodb.com/documentation/cartodb- js.html so that I can put a customized map into a wordpress site. I've been copying and pasting the following code into the wordpress text editor. However, it is not displaying a map on my website.

           How do I display a map in wordpress, using cartodb.js?


title: GEOTIFF not showing up in Tile Mill

text: All: Forgive any lapses in terminology here, although English is my native language, GIS most definitely is not, and I am _totally_ lost. Basically, I'm trying to convert a GEOTiff to a mbtile. I downloaded a sectional from the FAA website, and its SRS is this: +proj=lcc +lat_1=46.66666666666666 +lat_2=41.33333333333334 +lat_0=42.3 +lon_0=-97 +x_0=0 +y_0=0 +datum=NAD83 +units=m +no_defs I added it as a layer in Tile Mill, and it didn't show up. I tried playing with style.mss, and then read that "since Mapnik is currently unable to reproject raster data sources, to load them in TileMill you must ensure they are in Web Mercator projection." First, does this apply to my situation? Is a GEOTiff considered a raster data source? If so, when I go to Raster > Projections > Warp, what do I use for the Target SRS? Thanks for any insights, Michael


title: Is it possible to rotate items in QGIS Print Composer?

text: Is it possible to rotate items in Quantum GIS Print Composer? I just didn't found and before submit feature demand just double-check this. Pekka


title: Rasterfile clipping error: Cannot compute bounding box of cutline

text: I have a raster file and want to clip it, a vector file is supposed to be the mask layer. I use the most obvious method (Raster->Extraction->Clipper). No matter what I do, I get this error (see pic): > **Cannot compute bounding box of cutline** I tried it with changing the raster format, changing the projection etc. Of course I made research in older posts/Google too. I have no Idea how to solve this problem. If I clip by extend it does work. Weird thing is, I remember doing this three months ago, same method and I did not get this error. I also would (instead of a solution for this error) welcome an alternative way to perform this action. I just need to clip a rasterfile as described. Can someone help? Thanks. ![error](http://i.stack.imgur.com/qWCmF.jpg)


title: Search By Location

text: Using Arcmap 9.3 I was curious if anybody has a script or vba tool that will select items from a layer that are not within a selected polygon or not touching the selected polygon. Thanks for your input.


title: Altering ArcGIS RoadClass Driving Directions for Walking Network (Network Analyst)

text: I'm using ArcGIS 10 Network Analyst Extension to try and incorporate their driving directions functionality into a walking network. I am wondering if there is any way to alter the RoadClass Values to be more suitable for a Walking Network. For example: **RoadClass Value** : 3 **Type of Road Class** : Ramps **Example** : Take ramp and go on US-59 N Is there any way that I can modify or add a RoadClass value to be "Stairs" and instead of saying "Take ramp..." it says "Traverse stairs..." or "Elevators" saying "Use elevator"? Thanks, eric


title: Adding metadata to a geographic gridded data-set and projecting it

text: I asked a question before that nobody answered...but I kind of figured it out myself. But I still remain with some things that I just don't seem to get right. Can anybody please find the time to help me? It really seems like something trivial... Thanks! - I want to transform a geographicly gridded latlon data-set into a polar stereographic projection with gdalwarp. For this I need to add the correct metadata to my gridded data-set, and I can't seem to get this right. It is regularly spaced with 0.04 degrees from 0:360 longitude and -86:-60 Latitude (Antarctica). Only thing I have is the grid, so it comes without a header. I however have added the current header to the dataset:               ncols        9001     nrows        651     xllcorner    0     yllcorner    -86.000000000000     dx           0.04     dy           0.04      The first problem I get is that when I gdalinfo it I get:               Corner Coordinates:     Upper Left  (   0.0000000, -59.9600000)      Lower Left  (   0.0000000, -86.0000000)      Upper Right (     360.040,     -59.960)      Lower Right (     360.040,     -86.000)      Center      ( 180.0200000, -72.9800000)       Which seems wrong (as it is 0.04 off). Possibly I have to use a gridsize of 9000/650 instead of 9001/651, but also this is confusing to me as I then omit a line. Do I set the header correctly? Other than defining the grid and the coordinates like this, what do I need to do with the projection data? Need I add -s_srs WGS84 (this is what I do) or something like -s_srs EPSG:3031 (South Polar Stereographic)? This gives me the following:               Coordinate System is:     GEOGCS["WGS 84",     DATUM["WGS_1984",         SPHEROID["WGS 84",6378137,298.257223563,             AUTHORITY["EPSG","7030"]],         AUTHORITY["EPSG","6326"]],     PRIMEM["Greenwich",0],     UNIT["degree",0.0174532925199433],     AUTHORITY["EPSG","4326"]]     Origin = (0.000000000000000,-59.960000000000001)     Pixel Size = (0.040000000000000,-0.040000000000000)      I want to finally project it onto the following system:               PROJCS["unnamed",     GEOGCS["WGS 84",         DATUM["WGS_1984",             SPHEROID["WGS 84",6378137,298.257223563,                 AUTHORITY["EPSG","7030"]],             AUTHORITY["EPSG","6326"]],         PRIMEM["Greenwich",0],         UNIT["degree",0.0174532925199433],         AUTHORITY["EPSG","4326"]],     PROJECTION["Polar_Stereographic"],     PARAMETER["latitude_of_origin",-71],     PARAMETER["central_meridian",0],     PARAMETER["scale_factor",1],     PARAMETER["false_easting",0],     PARAMETER["false_northing",0],     UNIT["metre",1,         AUTHORITY["EPSG","9001"]]]      Origin = (-2800500.000000000000000,2800500.000000000000000)     Pixel Size = (1000.000000000000000,-1000.000000000000000)     Metadata:     AREA_OR_POINT=Area     Image Structure Metadata:      INTERLEAVE=BAND     Corner Coordinates:     Upper Left  (-2800500.000, 2800500.000) ( 45d 0' 0.00"W, 54d40' 3.36"S)     Lower Left  (-2800500.000,-2800500.000) (135d 0' 0.00"W, 54d40' 3.36"S)     Upper Right ( 2800500.000, 2800500.000) ( 45d 0' 0.00"E, 54d40' 3.36"S)     Lower Right ( 2800500.000,-2800500.000) (135d 0' 0.00"E, 54d40' 3.36"S)     Center      (   0.0000000,   0.0000000) (  0d 0' 0.01"E, 90d 0' 0.00"S)      Now there are 2 things occuring with my results when I warp it like this:               gdalwarp -s_srs WGS84 -r cubic data1.txt data2.tif      (i use either an existing data-set (I want to remap so...) or set the -te and -tr parameters, and get the same results....)   1. I get a warped file which fails to do the entire file (there is probably something going wrong at the 0/360 longitude transition. As expected it omits the value that falls at -360.04   2. I get only half of the file!! Only the right side of the projection is done... So my questions are:   1. How do i correctly set the correct projection parameters and how do I deal with polar stereographic data? Why can't I properly set the latitude and longitude range?   2. What in earth is going wrong when it only projects the right side? Thanks for the help!


title: How to configure IIS site as application for ArcGIS Server 10.2?

text: I am having several errors displayed after installing the webadaptor onto my site. They all contain this error message "This error can be caused by a virtual directory not being configured as an application in IIS" all issues point to my web.config file. If it's of any help the first two of four errors are caused by these two lines so it leads me to believe that it is being caused by the aspnet_client. I have already made this into an application so I'm a bit lost at this point.               <authentication mode="Windows" />     <membership defaultProvider="AspNetSqlMembershipProvider">


title: Literature on the appoximation of trajectory points with (smooth) curves

text: We are gathering trajectories of GPS points of moving objects (cars, bikes, etc.) and would like to compute approximate curves (that may self-intersect) that would represent those clouds of points. What I am looking is a literature lists (preferably textbooks if not, papers) that I need to familiarize myself in order to solve this problem. Any comments?


title: Read raster value at mouse position

text: I have raster layer imported in QGIS. How can I read value of the variable represented by the raster at mouse position? I just have a new raster layer and I want to read values at random places to have a basic idea about the raster. I prefer not to create point layer for this purpose, but if it is needed, please tell me how to do it that way.


title: How to get sum of table column values into Python variable?

text: Im creating a script to calculate some new fields in Arcgis10.1 using the arcpy library. I need to sum the total of a table column in order to calculate a new field. I wanted to ask if it is possible to store all the values in an array so i can get the total sum and use the variable in the new field calculation?


title: Reading FileGDB (9.2-9.3.1) with OGR

text: I have been tasked with processing USGS's National Hydrography Dataset. However, it appears that the only available format is FileGDB (version 9.2-9.3.1). Looking at OGRs support for FileGDB, it looks like they only support those created by ArcGIS 10 and above. Does anybody have any ideas how I would be able to view and export this data to a more usable format? If anybody is curious about the data, it can be directly downloaded through here.


title: How to bypass errors in arcpy for/while loop?

text: I have a handy script tool that loops through a workspace and renames and copies shapefiles to a feature dataset. However, if there is a corrupted shapefile somewhere in the workspace the script fails and stops processing. How do you handle errors such as this? Is there a way to print the error file and continue processing the next shapefile in the for loop to completion?               import arcpy     from arcpy import env          # Allow overwriting of output       env.overwriteOutput = True          # Parameters       env.workspace = arcpy.GetParameterAsText(0)      state = arcpy.GetParameterAsText(1)     gdb = arcpy.GetParameterAsText(2)          # Get a list of shapefiles in folder       fcs = arcpy.ListFeatureClasses()           # Find the total count of shapefiles in list       fcCount = len(fcs)           # Set the progressor      arcpy.SetProgressor("step", "Copying shapefiles to geodatabase...", 0,fcCount, 1)           # For each shapefile, copy to a file geodatabase          try:         for shp in fcs:                        # Define name for the output points              fc = str(state + shp[0:9])                  # Update the progressor label for current shapefile               arcpy.SetProgressorLabel("Loading " + shp + "...")                   # Copy the data               arcpy.CopyFeatures_management(shp, str(gdb + "\\" + fc))                   # Update the progressor position               arcpy.SetProgressorPosition()          except Exception as e:         print "An error has occurred"         print e          arcpy.ResetProgressor()


title: Calculation of average point distances from one another within a polygon grid

text: I am working with processed LiDAR data of forested areas in ArcMap 10, with various area based statistics extracted for each 30x30m polygon grid-cell in a shapefile. There are approximately 20,000 of these cells. I also have a layer of points relating to individual tree locations for the whole site. What I need to do is calculate the average distance between the tree-points and their standard deviations within each grid-cell polygon. And then write those two values to the corresponding grid-cell. I have used the GME toolbox (http://www.spatialecology.com/gme/) to derive this data for areas corresponding to my field site locations; however this is simply not feasible over such a large area. Does anyone have any ideas on how I can solve this issue?


title: ArcGis Raster Resampling Question

text: > **Possible Duplicate:**   >  How do you decide what interpolation method to use for resampling raster > data? I want to interpolate a certiain target variable over a large territory. I trained a random forest using a sample of points on that territory where the target is known (I have many regressor layers / variables). The information of these regressors is available over the whole territory (thus the capability to interpolate the target). The thing is that the regressor layers are available in different resolutions, I used resample to force them to the same resolution (1km^2). Overall the model works rather well (This on a 20% test set). And the output map looks good. But I wanted to know if the different techiniques available for doing resampling have anything to do with the quality of the final output map? What are the pros and cons of the different resampling techiniques available in ArcGis (Nearest neighbor assignment, Bilinear interpolation, Cubic convolution). I'm only referring to continuous variables. From ArcGis Help: Bilinear interpolation or cubic convolution should not be used on categorical data since the categories will not be maintained in the output raster dataset. However, all three techniques can be applied to continuous data, with nearest neighbor producing a blocky output, bilinear interpolation producing smoother results, and cubic convolution producing the sharpest.


title: Stop ArcMap from probing missing data sources on launch

text: When migrating a map from one environment to another (say, from an existing client environment to our development environment) we're often given MXD documents that represent the map. These documents are often connected to SDE data sources. Generally it is not too much effort to load the data into our own database, launch the MXD, wait for it to finish not finding any data, and re-path the layers. However, if the map document has dozens or even hundreds of layers, it can be a long wait. Is there any way to stop ArcMap from even trying to load those layers? A flag somewhere that says "I know the data are missing, just skip that step"? This applies to both ArcGIS 9.3.1 and 10.


title: Getting WFS data from Geoserver into leaflet

text: I'm trying to use a ajax call to load WFS data from geoserver into a GeoJson layer in leaflet. When I do the following message in the console:               XMLHttpRequest cannot load      http://localhost:8080/geoserver/cite/ows?service=WFS&version=1.1.0&request=…rmat=text/javascript&outputFormat=json&format_options=callback:loadGeoJson.      Origin http://localhost is not allowed by Access-Control-Allow-Origin.       I know that if you are trying to load wfs into OpenLayers you need a proxy.cgi as per this thread: XMLHttpRequest cannot load http://localhost: /geoserver/wfs. Origin http://localhost is not allowed by Access-Control- Allow-Origin That script seems to relate only to OpenLayers. My question is can I use this for leaflet and if not how to I get the wfs data? There must be a way round this problem but I can't see anything in the main documentation and I've searched almost exhausted my searching elsewhere. * * * I have added my script for the leaflet page...               var map = L.map('map').setView([50.370, -4.14], 14);               var basemap = L.tileLayer('http://{s}.tile.cloudmade.com/a8ed148d84c2438ba62c2742cac051ae/33356/256/{z}/{x}/{y}.png', {         maxZoom: 18,         attribution: 'Map data © OpenStreetMap contributors, CC-BY-SA, Imagery Â© CloudMade'     }).addTo(map);          var myLayer = new L.geoJson(null,{     onEachFeature: function      (feature, layer)      {layer.bindPopup(feature.properties.comment);}      }).addTo(map);               function loadGeoJson(data) {myLayer.addData(data);      }          var geoJsonUrl = "http://localhost:8080/geoserver/cite/ows?service=WFS&version=1.1.0&request=GetFeature&typeName=cite:masts&srsName=EPSG:4326&format=text/javascript&outputFormat=json&format_options=callback:loadGeoJson";     $.ajax({              url: geoJsonUrl,              dataType: 'jsonp'      });       and the httpd.conf...               #     # This is the main Apache HTTP server configuration file.  It contains the     # configuration directives that give the server its instructions.     # See <URL:http://httpd.apache.org/docs/2.2> for detailed information.     # In particular, see      # <URL:http://httpd.apache.org/docs/2.2/mod/directives.html>     # for a discussion of each configuration directive.     #     # Do NOT simply read the instructions in here without understanding     # what they do.  They're here only as hints or reminders.  If you are unsure     # consult the online docs. You have been warned.       #     # Configuration and logfile names: If the filenames you specify for many     # of the server's control files begin with "/" (or "drive:/" for Win32), the     # server will use that explicit path.  If the filenames do *not* begin     # with "/", the value of ServerRoot is prepended -- so "logs/foo.log"     # with ServerRoot set to "C:/Program Files/Apache Software Foundation/Apache2.2" will be interpreted by the     # server as "C:/Program Files/Apache Software Foundation/Apache2.2/logs/foo.log".     #     # NOTE: Where filenames are specified, you must use forward slashes     # instead of backslashes (e.g., "c:/apache" instead of "c:\apache").     # If a drive letter is omitted, the drive on which Apache.exe is located     # will be used by default.  It is recommended that you always supply     # an explicit drive letter in absolute paths to avoid confusion.          #     # ServerRoot: The top of the directory tree under which the server's     # configuration, error, and log files are kept.     #     # Do not add a slash at the end of the directory path.  If you point     # ServerRoot at a non-local disk, be sure to point the LockFile directive     # at a local disk.  If you wish to share the same ServerRoot for multiple     # httpd daemons, you will need to change at least LockFile and PidFile.     #     ServerRoot "c:/wamp/bin/apache/apache2.2.22"          #     # Listen: Allows you to bind Apache to specific IP addresses and/or     # ports, instead of the default. See also the <VirtualHost>     # directive.     #     # Change this to Listen on specific IP addresses as shown below to      # prevent Apache from glomming onto all bound IP addresses.     #     #Listen 12.34.56.78:80     Listen 80                    #     # Dynamic Shared Object (DSO) Support     #     # To be able to use the functionality of a module which was built as a DSO you     # have to place corresponding `LoadModule' lines at this location so the     # directives contained in it are actually available _before_ they are used.     # Statically compiled modules (those listed by `httpd -l') do not need     # to be loaded here.     #     # Example:     # LoadModule foo_module modules/mod_foo.so     #     LoadModule actions_module modules/mod_actions.so     LoadModule alias_module modules/mod_alias.so     LoadModule asis_module modules/mod_asis.so     LoadModule auth_basic_module modules/mod_auth_basic.so     #LoadModule auth_digest_module modules/mod_auth_digest.so     #LoadModule authn_alias_module modules/mod_authn_alias.so     #LoadModule authn_anon_module modules/mod_authn_anon.so     #LoadModule authn_dbd_module modules/mod_authn_dbd.so     #LoadModule authn_dbm_module modules/mod_authn_dbm.so     LoadModule authn_default_module modules/mod_authn_default.so     LoadModule authn_file_module modules/mod_authn_file.so     #LoadModule authnz_ldap_module modules/mod_authnz_ldap.so     #LoadModule authz_dbm_module modules/mod_authz_dbm.so     LoadModule authz_default_module modules/mod_authz_default.so     LoadModule authz_groupfile_module modules/mod_authz_groupfile.so     LoadModule authz_host_module modules/mod_authz_host.so     #LoadModule authz_owner_module modules/mod_authz_owner.so     LoadModule authz_user_module modules/mod_authz_user.so     LoadModule autoindex_module modules/mod_autoindex.so     #LoadModule cache_module modules/mod_cache.so     #LoadModule cern_meta_module modules/mod_cern_meta.so     LoadModule cgi_module modules/mod_cgi.so     #LoadModule charset_lite_module modules/mod_charset_lite.so     #LoadModule dav_module modules/mod_dav.so     #LoadModule dav_fs_module modules/mod_dav_fs.so     #LoadModule dav_lock_module modules/mod_dav_lock.so     #LoadModule dbd_module modules/mod_dbd.so     LoadModule deflate_module modules/mod_deflate.so     LoadModule dir_module modules/mod_dir.so     #LoadModule disk_cache_module modules/mod_disk_cache.so     #LoadModule dumpio_module modules/mod_dumpio.so     LoadModule env_module modules/mod_env.so     #LoadModule expires_module modules/mod_expires.so     #LoadModule ext_filter_module modules/mod_ext_filter.so     #LoadModule file_cache_module modules/mod_file_cache.so     #LoadModule filter_module modules/mod_filter.so     LoadModule headers_module modules/mod_headers.so     #LoadModule ident_module modules/mod_ident.so     #LoadModule imagemap_module modules/mod_imagemap.so     LoadModule include_module modules/mod_include.so     #LoadModule info_module modules/mod_info.so     LoadModule isapi_module modules/mod_isapi.so     #LoadModule ldap_module modules/mod_ldap.so     #LoadModule logio_module modules/mod_logio.so     LoadModule log_config_module modules/mod_log_config.so     #LoadModule log_forensic_module modules/mod_log_forensic.so     #LoadModule mem_cache_module modules/mod_mem_cache.so     LoadModule mime_module modules/mod_mime.so     #LoadModule mime_magic_module modules/mod_mime_magic.so     LoadModule negotiation_module modules/mod_negotiation.so     LoadModule proxy_module modules/mod_proxy.so     LoadModule proxy_ajp_module modules/mod_proxy_ajp.so     LoadModule proxy_balancer_module modules/mod_proxy_balancer.so     LoadModule proxy_connect_module modules/mod_proxy_connect.so     #LoadModule proxy_ftp_module modules/mod_proxy_ftp.so     LoadModule proxy_http_module modules/mod_proxy_http.so     #LoadModule rewrite_module modules/mod_rewrite.so     LoadModule setenvif_module modules/mod_setenvif.so     #LoadModule speling_module modules/mod_speling.so     #LoadModule ssl_module modules/mod_ssl.so     #LoadModule status_module modules/mod_status.so     #LoadModule substitute_module modules/mod_substitute.so     #LoadModule unique_id_module modules/mod_unique_id.so     #LoadModule userdir_module modules/mod_userdir.so     #LoadModule usertrack_module modules/mod_usertrack.so     #LoadModule version_module modules/mod_version.so     #LoadModule vhost_alias_module modules/mod_vhost_alias.so     LoadModule php5_module "c:/wamp/bin/php/php5.4.3/php5apache2_2.dll"          <IfModule !mpm_netware_module>     <IfModule !mpm_winnt_module>     #     # If you wish httpd to run as a different user or group, you must run     # httpd as root initially and it will switch.       #     # User/Group: The name (or #number) of the user/group to run httpd as.     # It is usually good practice to create a dedicated user and group for     # running httpd, as with most system services.     #     User daemon     Group daemon          </IfModule>     </IfModule>          # 'Main' server configuration     #     # The directives in this section set up the values used by the 'main'     # server, which responds to any requests that aren't handled by a     # <VirtualHost> definition.  These values also provide defaults for     # any <VirtualHost> containers you may define later in the file.     #     # All of these directives may appear inside <VirtualHost> containers,     # in which case these default settings will be overridden for the     # virtual host being defined.     #          #     # ServerAdmin: Your address, where problems with the server should be     # e-mailed.  This address appears on some server-generated pages, such     # as error documents.  e.g. admin@your-domain.com     #     ServerAdmin admin@localhost          #     # ServerName gives the name and port that the server uses to identify itself.     # This can often be determined automatically, but we recommend you specify     # it explicitly to prevent problems during startup.     #     # If your host doesn't have a registered DNS name, enter its IP address here.     #     ServerName localhost:80          #     # DocumentRoot: The directory out of which you will serve your     # documents. By default, all requests are taken from this directory, but     # symbolic links and aliases may be used to point to other locations.     #     DocumentRoot "c:/wamp/www/"          #     # Each directory to which Apache has access can be configured with respect     # to which services and features are allowed and/or disabled in that     # directory (and its subdirectories).      #     # First, we configure the "default" to be a very restrictive set of      # features.       #     <Directory />         Options FollowSymLinks         AllowOverride None         Order deny,allow         Deny from all     </Directory>          ProxyRequests Off     ProxyPreserveHost On          <Proxy *>         Order Deny,Allow         Deny from all         Allow from 192.168.11     </Proxy>          ProxyPass /geoserver http://localhost:8080/geoserver/web     ProxyPassReverse /geoserver http://localhost:8080/geoserver/web     <Location /geoserver>         Order allow,deny         Allow from all     </Location>          #     # Note that from this point forward you must specifically allow     # particular features to be enabled - so if something's not working as     # you might expect, make sure that you have specifically enabled it     # below.     #          #     # This should be changed to whatever you set DocumentRoot to.     #     <Directory "c:/wamp/www/">         #         # Possible values for the Options directive are "None", "All",         # or any combination of:         #   Indexes Includes FollowSymLinks SymLinksifOwnerMatch ExecCGI MultiViews         #         # Note that "MultiViews" must be named *explicitly* --- "Options All"         # doesn't give it to you.         #         # The Options directive is both complicated and important.  Please see         # http://httpd.apache.org/docs/2.2/mod/core.html#options         # for more information.         #         Options -Indexes FollowSymLinks              #         # AllowOverride controls what directives may be placed in .htaccess files.         # It can be "All", "None", or any combination of the keywords:         #   Options FileInfo AuthConfig Limit         #         AllowOverride all              #         # Controls who can get stuff from this server.         #          #   onlineoffline tag - don't remove         AllowOverride All         Order allow,deny         Allow from all          </Directory>          #     # DirectoryIndex: sets the file that Apache will serve if a directory     # is requested.     #     <IfModule dir_module>         DirectoryIndex index.php index.php3 index.html index.htm     </IfModule>          #     # The following lines prevent .htaccess and .htpasswd files from being      # viewed by Web clients.      #     <FilesMatch "^\.ht">         Order allow,deny         Deny from all         Satisfy All     </FilesMatch>          #     # ErrorLog: The location of the error log file.     # If you do not specify an ErrorLog directive within a <VirtualHost>     # container, error messages relating to that virtual host will be     # logged here.  If you *do* define an error logfile for a <VirtualHost>     # container, that host's errors will be logged there and not here.     #     ErrorLog "c:/wamp/logs/apache_error.log"          #     # LogLevel: Control the number of messages logged to the error_log.     # Possible values include: debug, info, notice, warn, error, crit,     # alert, emerg.     #     LogLevel warn          <IfModule log_config_module>         #         # The following directives define some format nicknames for use with         # a CustomLog directive (see below).         #         LogFormat "%h %l %u %t \"%r\" %>s %b \"%{Referer}i\" \"%{User-Agent}i\"" combined         LogFormat "%h %l %u %t \"%r\" %>s %b" common              <IfModule logio_module>           # You need to enable mod_logio.c to use %I and %O           LogFormat "%h %l %u %t \"%r\" %>s %b \"%{Referer}i\" \"%{User-Agent}i\" %I %O" combinedio         </IfModule>              #         # The location and format of the access logfile (Common Logfile Format).         # If you do not define any access logfiles within a <VirtualHost>         # container, they will be logged here.  Contrariwise, if you *do*         # define per-<VirtualHost> access logfiles, transactions will be         # logged therein and *not* in this file.         #         CustomLog "c:/wamp/logs/access.log" common              #         # If you prefer a logfile with access, agent, and referer information         # (Combined Logfile Format) you can use the following directive.         #         #CustomLog "logs/access.log" combined     </IfModule>          <IfModule alias_module>         #         # Redirect: Allows you to tell clients about documents that used to          # exist in your server's namespace, but do not anymore. The client          # will make a new request for the document at its new location.         # Example:         # Redirect permanent /foo http://localhost/bar              #         # Alias: Maps web paths into filesystem paths and is used to         # access content that does not live under the DocumentRoot.         # Example:         # Alias /webpath /full/filesystem/path         #         # If you include a trailing / on /webpath then the server will         # require it to be present in the URL.  You will also likely         # need to provide a <Directory> section to allow access to         # the filesystem path.              #         # ScriptAlias: This controls which directories contain server scripts.          # ScriptAliases are essentially the same as Aliases, except that         # documents in the target directory are treated as applications and         # run by the server when requested rather than as documents sent to the         # client.  The same rules about trailing "/" apply to ScriptAlias         # directives as to Alias.         #         ScriptAlias /cgi-bin/ "cgi-bin/"          </IfModule>          <IfModule cgid_module>         #         # ScriptSock: On threaded servers, designate the path to the UNIX         # socket used to communicate with the CGI daemon of mod_cgid.         #         #Scriptsock logs/cgisock     </IfModule>          #     # "C:/Program Files/Apache Software Foundation/Apache2.2/cgi-bin" should be changed to whatever your ScriptAliased     # CGI directory exists, if you have that configured.     #     <Directory "cgi-bin">         AllowOverride None         Options None         Order allow,deny         Allow from all     </Directory>          #     # DefaultType: the default MIME type the server will use for a document     # if it cannot otherwise determine one, such as from filename extensions.     # If your server contains mostly text or HTML documents, "text/plain" is     # a good value.  If most of your content is binary, such as applications     # or images, you may want to use "application/octet-stream" instead to     # keep browsers from trying to display binary files as though they are     # text.     #     DefaultType text/plain          <IfModule mime_module>         #         # TypesConfig points to the file containing the list of mappings from         # filename extension to MIME-type.         #         TypesConfig conf/mime.types              #         # AddType allows you to add to or override the MIME configuration         # file specified in TypesConfig for specific file types.         #         #AddType application/x-gzip .tgz         #         # AddEncoding allows you to have certain browsers uncompress         # information on the fly. Note: Not all browsers support this.         #         #AddEncoding x-compress .Z         #AddEncoding x-gzip .gz .tgz         #         # If the AddEncoding directives above are commented-out, then you         # probably should define those extensions to indicate media types:         #         AddType application/x-compress .Z         AddType application/x-gzip .gz .tgz         AddType application/x-httpd-php .php         AddType application/x-httpd-php .php3              #         # AddHandler allows you to map certain file extensions to "handlers":         # actions unrelated to filetype. These can be either built into the server         # or added with the Action directive (see below)         #         # To use CGI scripts outside of ScriptAliased directories:         # (You will also need to add "ExecCGI" to the "Options" directive.)         #         #AddHandler cgi-script .cgi              # For type maps (negotiated resources):         #AddHandler type-map var              #         # Filters allow you to process content before it is sent to the client.         #         # To parse .shtml files for server-side includes (SSI):         # (You will also need to add "Includes" to the "Options" directive.)         #         #AddType text/html .shtml         #AddOutputFilter INCLUDES .shtml     </IfModule>          #     # The mod_mime_magic module allows the server to use various hints from the     # contents of the file itself to determine its type.  The MIMEMagicFile     # directive tells the module where the hint definitions are located.     #     #MIMEMagicFile conf/magic          #     # Customizable error responses come in three flavors:     # 1) plain text 2) local redirects 3) external redirects     #     # Some examples:     #ErrorDocument 500 "The server made a boo boo."     #ErrorDocument 404 /missing.html     #ErrorDocument 404 "/cgi-bin/missing_handler.pl"     #ErrorDocument 402 http://localhost/subscription_info.html     #          #     # EnableMMAP and EnableSendfile: On systems that support it,      # memory-mapping or the sendfile syscall is used to deliver     # files.  This usually improves server performance, but must     # be turned off when serving from networked-mounted      # filesystems or if support for these functions is otherwise     # broken on your system.     #     #EnableMMAP off     #EnableSendfile off          # Supplemental configuration     #     # The configuration files in the conf/extra/ directory can be      # included to add extra features or to modify the default configuration of      # the server, or you may simply copy their contents here and change as      # necessary.          # Server-pool management (MPM specific)     #Include conf/extra/httpd-mpm.conf          # Multi-language error messages     #Include conf/extra/httpd-multilang-errordoc.conf          # Fancy directory listings     Include conf/extra/httpd-autoindex.conf          # Language settings     #Include conf/extra/httpd-languages.conf          # User home directories     #Include conf/extra/httpd-userdir.conf          # Real-time info on requests and configuration     #Include conf/extra/httpd-info.conf          # Virtual hosts     #Include conf/extra/httpd-vhosts.conf          # Local access to the Apache HTTP Server Manual     #Include conf/extra/httpd-manual.conf          # Distributed authoring and versioning (WebDAV)     #Include conf/extra/httpd-dav.conf          # Various default settings     #Include conf/extra/httpd-default.conf          # Secure (SSL/TLS) connections     #Include conf/extra/httpd-ssl.conf     #     # Note: The following must must be present to support     #       starting without SSL on platforms with no /dev/random equivalent     #       but a statically compiled-in mod_ssl.     #     <IfModule ssl_module>     SSLRandomSeed startup builtin     SSLRandomSeed connect builtin     </IfModule>               #     # Configure the proxy module of Apache     #     ProxyRequests On          <Proxy *>     Order deny,allow     Allow from all     </Proxy>          ProxyPass /geoserver http://localhost:8080/geoserver     ProxyPreserveHost On     ProxyStatus On              Include "c:/wamp/alias/*"


title: calculating home ranges with QGIS

text: i am facing some problems when trying to calculate animals home ranges with QGIS 1.8.0. One the one hand, when trying to calculate MCP with the Home range plugin, i get error messages. on the other hand, I am trying to download the AniMove SEXTANTE plugin but I cannot find it in the QGIS Ptyhon plugin installer! I would aprechiate detailed instructions how to calculate home ranges! Thank you very much!


title: Info from multiple rows in single popup

text: If you have a table like this: `Name | Item | State ---------|--------|----------- Alice | Apple | Alabama Ben | Banana | Alabama Caitlin | Cherry | California ` Assuming you have created polygons based on the states, how could you create an 'info window' which combines the details of all rows in the same state? So for example, the popup for Alabama might be something like:               State: Alabama     People: 2     Items: Apple, Banana      As far as I can tell, the default behaviour for Cartodb would be to place one polygon on top of the other, and for the info window to display the info from the top most polygon only.


title: The arctic circle moves?

text: I heard today that the Arctic Circle, along with it's siblings the Antarctic Circle and the Tropics of Cancer and Capricorn, are different from year to year. Is that true? Can you please explain how/why this is so? I'd always thought they were defined by particular latitudinal parallels. e.g. 66°N or S etc.


title: Render GPX track data into video

text: I have two geeky passions while driving: recording GPX data to upload to OSM as well as recording on board videos like this one. Since these kind of videos are sort of boring for others to watch I was wondering if maybe there's an application out there to convert the track data into a video which I could then overlay over my original recording showing real time stuff like time of day, speed, elevation, etc. Basically a HUD. If there was a way to show part of the map I'm driving through that would be dreamy. If you know of such a piece of software, please help. Thanks in advance!


title: What tools do you recommend for this kind of web map?

text: I have one tourist map (ski resort) with buildings and roads in .jpg format. It's problem because north is down and south is up (it is that kind of drawn map with drawn mountain up on map where is actually south). I need to put it on the website and also put some pins on it. What tools and ways do you recommend me? It needs to be done fast so please give me some simple solutions if possible. Thank you.


title: query to get predicted location after traveling on pgrouting route for some hours, postgis 2, postgresql 9.3

text: I am using Postgis 2.0 on Postgresql 9.3 on Ubuntu 14.04 x64. I am trying to find an efficient way to fulfill this routing request. So I am able to get the shortest route using the `source` and `target` locations. However once I have that route, I would like to know where a car or person will be `x` hours after starting from the `source` location. I know that the cost function has the travel time in hours, so I should be able to find the longitude and latitude location after traveling for some number of hours. However, I am not familiar enough with pgrouting to know the most efficient way to run this kind of query. Any suggestions?


title: How to change attributes with QGIS python?

text: I have this piece of code which I would like to use for inserting values into column. It does not throw an error, the changes are not written to the source (shapefile).               from qgis.utils import iface     from PyQt4.QtCore import *          layers = iface.legendInterface().layers()          for layer in layers:         name = layer.name()         if name.endswith('x'):             provider = layer.dataProvider()                features = provider.getFeatures()             for feature in features:                 feature.setAttribute('attr', 'a')      I understood there are two ways of editing datasources: directly with dataProvider or with edit buffer. If I got it right, it should'n be necessary to call commit action when using dataProvider. Any help really appreciated, thank you (using QGIS 2.2 Valmiera).


title: Rotating the coordinates of a layer by 90 degrees

text: I have got a dwg file with the plans for a building, with rooms and everything. (There's been some samples taken inside the building which I want to display.) Now, since it's of a building, the person making this dwg file didn't care much for coordinate systems, and it's all in some (supposedly) arbitrary system centered around the building. This is not an issue. The problem is that this file (which I have converted to shapefiles) does not honor North being upwards, and that is what I now need to fix. Basically I want all my shapefiles to be rotated 90 degrees, so that North is upwards. It doesn't really help to rotate my data frame, since that would only change North inside my document and not for my shapefiles' coordinates. Can I change the content of the prj-file in some way, or how could I achieve this? Using ArcMap 10.2 basic, Python 2.7, FME 2014, Win 7.


title: Implement security to qGIS Server (for windows)

text: I have few map services running on QGIS Server(v 2.2.X on Windows 7 32 bit) and anyone can connect or access those services (WMS & WFS). Is there any way to implement security in QGIS Server? So that I can set, assign username & password to those map services. I researched on my side and following few links in which provides few steps such as apply password to specific folder, creating .htaccess file, also researched official site of QGIS Server but cannot find any valid document. Any help will be great!!! Thanks in advance !!!


title: Command Line alternative to using QGIS GUI to save a Raster as a Rendered Image

text: I have a raster image in PNG format which I can open in _QGIS but not in Windows Picture Viewer_. To do that i figured out that I have to save my raster as a **_'Rendered Image'_** from the **_'Save As'_** (option by right- clicking on it in QGIS GUI). So great then I could now open my rendered PNG in Windows Picture Viewer. _"Now I need a little help with as to how to do this (i.e. save the Raster as Rendered Image) without the QGIS GUI and via a Command Line"_


title: Determining whether point location is left or right of line?

text: I have one line layer and one point layer. Points are located both sides of line. I want to know whether points are located on the right side or the left of the line. Is it possible to determine this using ArcGIS for Desktop or QGIS?


title: Comparing routes

text: I am looking for a way to compare routes of 2 or more persons to be taken by car. Person 1 has a starting point A1 and a destination point A2.   Person 2 has a starting point B1 and a destination point B2.   [...] What I would like to accomplish:   1. Check, if A1 and B1 are close to each other.   2. Calculate walking and/or driving distance from A1 to B1.   3. Save the individual routes of Person 1 and Person 2 in a way, that allows for performant comparison to other routes / route points.   4. Identify, if the routes to be taken of Person 1 and Person 2 are similar enough, to have at least same route parts. Are there already (open source) examples of code around that I haven't found yet? Or can you give me advice on how to approach this problem? Any help is highly appreciated. Thanks!