Levelplot(RasterVis) two marginal plot scale?

![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()


Conversion table convert XY coordinates to GPS

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


Excellent app for OSM browsing on smartphone

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


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

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)


Deleting vector lines beyond image extents using QGIS?

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?


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

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)


How to find nearby objects using Easting and Northing range

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.


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

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!


CartoDB get rows in initialized map

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?


ArcPy - Add Layer from Cursor Selection

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()


refactoring a osmfilter to do a better work

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!?


ERDAS Imagine ecw compression settings

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


Reprojecting vector layer in OpenLayers

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!


Calculate distance and surface with mapbox api

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.


Geoserver layer rotation problem

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}                  }         );


Convert Sql Server Spatial data to GeoJSON & plot on map

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]         }     };


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

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.


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

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.


Conditional SQL Query

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?


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

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?


GEOTIFF not showing up in Tile Mill

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


Is it possible to rotate items in QGIS Print Composer?

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


Rasterfile clipping error: Cannot compute bounding box of cutline

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)


Search By Location

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.


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

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


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

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!


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

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">


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

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?


Read raster value at mouse position

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.


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

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?


Reading FileGDB (9.2-9.3.1) with OGR

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.


How to bypass errors in arcpy for/while loop?

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()


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

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?


ArcGis Raster Resampling Question

> **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.


Stop ArcMap from probing missing data sources on launch

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.


Getting WFS data from Geoserver into leaflet

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/*"


calculating home ranges with QGIS

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!


Info from multiple rows in single popup

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.


The arctic circle moves?

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.


Render GPX track data into video

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!


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

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.


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

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?


How to change attributes with QGIS python?

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).


Rotating the coordinates of a layer by 90 degrees

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.


Implement security to qGIS Server (for windows)

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 !!!


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

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"_


Determining whether point location is left or right of line?

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?


Comparing routes

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!


I want to show the information of several layers of a request made GeoJSON geoserver

Always enters the error function and shows me the alert 'what happen?'. This is the piece of code that causes problems:               map.addEventListener('click',onMapClick);     function onMapClick(e)     {         var latlngStr = '('+ e.latlng.lat.toFixed(3) + ", " + e.latlng.lng.toFixed(3)+')';         var BBOX = map.getBounds()._southWest.lng + ", " + map.getBounds()._southWest.lat + ", " + map.getBounds()._northEast.lng + ", " + map.getBounds()._northEast.lat;         var WIDTH = map.getSize().x;         var HEIGHT = map.getSize().y;            var X = map.layerPointToContainerPoint(e.layerPoint).x.toFixed(0);         var Y = map.layerPointToContainerPoint(e.layerPoint).y.toFixed(0);         URL = 'http://192.168.60.55:9090/geoserver/ems/wms?request=GetFeatureInfo&service=WMS&version=1.1.1&info_format=application/json&exceptions=application/vnd.ogc.se_xml&layers=ems:tblsevilla,ems:tblmanz,ems:tblcarre&styles=&srs=EPSG%3A4326&format=image%2Fpng&bbox='+ BBOX + '&width='+ WIDTH +'&height='+ HEIGHT + '&query_layers=ems:tblsevilla,ems:tblmanz,ems:tblcarre&feature_count=20&x='+ X +'&y='+ Y;                   function onEachFeature(feature, layer) {             if (feature.properties && feature.properties.PERIMETRO)             layer.bindPopup("Perimetro: "+feature.properties.PERIMETRO);             }         }              $.ajax({             type: 'POST',             url: URL,             dataType : 'json',             success: function (response) {                 alert("2");                 L.geoJson(response, {                     onEachFeature: onEachFeature                 }).addTo(map);             },             error: function (req, status, error) {                 alert('what happen?');             }         });     }      'm Leaflet and geoserver beginner, any help would be welcome. Thank you very much in advance to all.


Is there a simple interface from which I can extract average summer temperature and humidity for a list of latitude and longitudes?

I am interested in finding the summer average (either climatalogical mean or for a specific year) relative humidity and temperature for a set of globally distributed latitude and longitude coordinates. What I need are mean values for June, July, and August. I would prefer to use actual observed values, but a reanalysis (such as NARR or CRU-NCEP would suffice). Is there a web-based interface from which I can obtain this information? The closest that I have come to this   * Wolfram Alpha, but they only provide rainfall and temperature, not relative humidity   * worldclim.org looks promising, but requires me to download all of the data before using another program to extract information; it also does not have relative humidity.


How to batch a model created with "model builder" in arcmap 10?

I created a simple model that selects iteratively elements from a point dataset and calculates a cost distance raster for each of these points. It works fine when I run it for just one dataset, but when I try to batch it get the following error. ERROR 000865: Input raster or feature source data: I_atlanta_OBJECTID does not exist. ERROR 000581: Invalid parameters. Failed to execute (Cost Distance). Please any help will be greatly appreciate. Here a link to a snapshot of the model http://img88.imageshack.us/img88/7432/modelmultidist.png


Why are steps out of order when running model in Batch?

I have an issue using an ArcMap model in batch mode. The purpose of my model is to obtain the result of an OD cost matrix from a selected network and a set of origins and destinations. My goal is to be able to get these results without having to add the data to the display in ArcMap (time and resource consuming), but most of all, I want to be able to automatically process several OD cost matrix one after another. For these reasons, I rebuilt the OD cost matrix routine using ArcMap's Model Builder. ![enter image description here](http://i.stack.imgur.com/6gAj9.png) ![enter image description here](http://i.stack.imgur.com/MJ9KL.png) If I test the model on only only network, it works great. However, when I try to batch the process where each time I have a different network with different locations, I have a problem since ArcCatalog seems to mix the different steps. For example, if I batch with only two different networks... ![enter image description here](http://i.stack.imgur.com/riG94.png) Instead of having:               1. Create OD cost layer (1)     2. Load Origins 1 in network 1     3. Load Destinations 1 in network 1     4. Solve OD cost matrix (1)     5. Create OD cost layer (2)     6. Load Origins 2 in network 2     7. Load Destinations 2 in network 2     8. Solve OD cost matrix (2)      ArcCatalog does this:               1. Create OD cost layer (1)     2. Create OD cost layer (2)     3. Load origins 1 in network 2     4. Load origins 2 in network 2     5. Load destinations 1 in network 2     6. Load destinations 2 in network 2     7. Solve OD cost matrix (1)     8. Solve OD cost matrix (2)      I thought that when batching a model ArcCatalog would do the first task completely before even looking at the second. Why are the steps of my model mixed up? Thanks for your help!


Why does ModelBuilder fail to clear batch parameter inputs and only use first batch line's input paramters?

I have built a model to convert coverage files to shapefiles and also feature classes in a geodatabase. The model uses a sub model to run an iterator to loop through an initial input folder (Which is set to level 2-102A, etc) which has multiple sub levels of folders:   1. TRIM   2. 102A, 102B,102C, etc   3. 102001P, 102002P, etc.   4. 102001P_hydrology, 102001P_contours, etc.   5. With files of various extensions ( point, tic, arc). Then I switch to the main model and use the select data tool to group the various file extensions and collect values tools. At this point I merge the values and create two paths, one which uses the feature class to shp. and the other which uses the feature class to GDB. I have set the initial folder, wildcard and output Shp folders and GDBs as parameters and when I run the model 1 input at a time it works fine however there are 23 intial input folders each with 7 subfolders so the use of the batch option would really save time however I am running in to problems.   1. The model appears to be retaining some memory of past operations or trails, meaning I tried to batch the entire process using value lists for the parameters, which failed to produce the proper outputs so I returned them all to single value inputs, yet the model still opens with the value lists present.   2. When I ran the entire batch it completed the process however when I checked my output values they were all copies of the original line, just sent to the various output location I had specified. I think this has to do with the fact the initial merge output is set as an intermediate file in a scratch.gdb so it is getting overwritten each time the model is run. I believe inline variable substitution can fix this yet I am unsure of how to do this as the model fails: invalid output name, when I tried the parse tool for the initial input file to use the name as an inline variable in the merge output name.  I have attached the model for clarity and my attempt at the inline variable substitution for the merge output. ![main model](http://i.stack.imgur.com/Fdlaa.jpg) ![output name using inline varaible substitution](http://i.stack.imgur.com/99uKO.jpg) Initially, I also had another pathway which converted the individual files from the output of the collect values tool into feature classes in my geodatabase, but after a while of trying to batch the model our link submodels each selecting for different wildcard the model would no longer create individual outputs it just produced a merge, even though all preconditions were set, and they were separate paths. Eventually I had to pull this path out of the model and run in a separate model. Is there a way to wipe all the past memory without deleting the entire model and beginning anew?


How do I pass each parsed file name to the next tool in Model Builder?

I cannot get the Parse Path tool to pass each file name to the next tool. In the example shown, the Raster to Ascii tool will take the first entry in the 'filename' variable list as part of the output name. But when it runs the next iteration, the Raster to Ascii tool uses some generic file naming convention. 'filename' is set as a precondition for the R2A tool. Example: raster1.tif raster2.tif Run the Model in Batch mode, using tifs as parameters. 1st output = raster1_asc.txt 2nd output = RastTo_rast.txt <\--- this is wrong. ![Model Builder running as Batch](http://i.stack.imgur.com/ruQc7.jpg) Summary: Raster to Ascii tool sucessfully uses the filename variable as inline variable as part of the output file name. However, the the Raster to Ascii tool uses a generic naming convention for the following files (running Model as Batch). Does anyone have an insight as to why this is happening?


is there any free software to create a mapserver that can be installed directly in a subdomain without a server install?

i want to create a webmap page using my own shapefiles but all the mapserver softwares that i've seen need to be installed directly in the server. so, is there any software that allows me to create the mapserver just by putting the files in some folder and "execute" them (like php or js, for example)? thanks in advance =^)


Is there a way to decrease redraws in ARCGIS?

I'm using a huge database on arcgis desktop 9.3.1 and every time I do a pan/zoom I have to wait 0.5-1 minute for the whole database to redraw itself. Is there any way to decrease the amount of redraws ( like making the program redraw only once every five pan actions)?


ArcGIS 10 SP1 Custom Overlay Grid: need to add a labeling expression

I have a custom overlay grid whose labels are driven by an ugly attribute field which contains unnecessary leading zeroes and other labeling no-nos. How can I introduce a labeling expression to this grid? Thank you!


Throwing E_FAIL errors so edit operation is aborted

I need to abort an edit operation somewhere, determined by a certain condition. Heres the code I'm using:                   private const int E_FAIL = -2147467259;              private void _Events_OnDeleteFeature(IObject obj)         {             _logger.Info("Iniciando handler de OnDeleteFeature...");             try             {                 // do we need to abort it?                 // this method returns a COMException with code = E_FAIL;                  if (DoINeedToAbort(obj, geaEditOperation.DELETE))                     return;             }             catch (COMException comEx)             {                 throw comEx;             }         }      This is not showing me the appropriate message and it's not aborting the operation. BTW, this event is a IEditEvent, and it's being fired ok. The problem is to abort the operation, without calling _editor.AbortEditOperation(), which gives me an ugly message. This is where I looked it in ESRIs manual (I know that the following applies to class extensions) http://resources.esri.com/help/9.3/arcgisengine/ArcObjects/ao_start.htm#esriGeodatabase/IObjectClassEvents.htm Any tips?


sum of raster counts

I'm trying to get the sum of "counts" in a raster. I have the code below which can successfully read the counts of each row in raster attribute and stores the counts in a list provided below. > [763461] [407207] [232269] [137428] [83909] [55074] [33340] [24168] [16747] > [11488] [7010] [4323] [4021] [2552] [1497] [775] [72] However, I want to get the sum of the counts. I am providing the code and error below. I would really appreciate your help/ suggestion! code:               import arcpy     from arcpy import env     from arcpy.sa import *          #To overwrite output     arcpy.env.overwriteOutput = True          #Set environment settings     env.workspace = "C:/Subhasis/Test/7-17-14/test"          #set local variables-STI extraction     inRaster ="sti"     inSQLClause = "VALUE >= 9"          #checkout ArcGIS spatial analyst extension license     arcpy.CheckOutExtension("Spatial")          # Execute ExtractByAttributes     attExtract = ExtractByAttributes(inRaster, inSQLClause)           # Save the output      attExtract.save("C:/Subhasis/Test/7-17-14/test/STI-9")          inraster ="sti-9"     rows = arcpy.SearchCursor(inraster,"","","Count","")          for row in rows:         count = row.getValue("Count")         s1=[]         s1.append(count)         print sum(range(s1))      Error:               Traceback (most recent call last):       File "C:/Subhasis/Test/7-17-14/test/test.py", line 31, in <module>         print sum(range(s1))     TypeError: range() integer end argument expected, got list.


Where can I download ArcGIS Service Pack 2?

I understand SP2 is available for download. http://blogs.esri.com/support/blogs/supportcenter/archive/2011/04/20/arcpad-10-0-service- pack-2-released.aspx > We are pleased to announce that Service Pack 2 for ArcPad 10.0 is now > available to download at the Esri Customer Care website. Maybe I am missing something but when I log into the ESRI customer care website the Service Pack is not listed. Anyone?


Creating data driven pages via arcpy

We have a data driven pages enabled MXD. Refreshing this mxd in arcmap works and gives me 9000+ pages. Via python however, _mxd.dataDrivenPages.pageCount_ returns 0. Even after _mxd.dataDrivenPages.refresh()_. In ArcGIS resources, they say: > Four elements are required for the creation of Data Driven Pages. They are: > > Detail Data Frame, Index Layer, Name Field, Sort Field In ArcMap I use the FILENAME option in the Name Field. In python, I find that _mxd.dataDrivenPages.pageNameField_ is _None_. But I'm unable to set this property. Any ideas on what to look for to get my pages via python? Thanks! Code is:               import arcpy, oss          logger = oss.log()          mxd = arcpy.mapping.MapDocument(MxdFile)     ddp = mxd.dataDrivenPages     logger.info('Checking required params')     logger.info('dataFrame is : "%s"' % mxd.dataDrivenPages.dataFrame)     logger.info('currentpageid is : "%s"' % mxd.dataDrivenPages.currentPageID)     logger.info('index layer is "%s"' % dpp.indexlayer)     logger.info('pagenamefield : "%s"' % mxd.dataDrivenPages.pageNameField)          logger.info('Refreshing index layer')     ddp.refresh()     logger.info('ddp number after refresh : "%s"' % mxd.dataDrivenPages.pageCount)      This results in:               INFO - Checking required params     INFO - dataFrame is : "<geoprocessing Data Frame object object at 0x0000000010CC     6970>"     INFO - currentpageid is : "0"     INFO - index layer is "VISION.SWG_WPA_MAP_ATLAS_TB"     INFO - pagenamefield : "None"     INFO - Refreshing index layer     INFO - ddp number after refresh : "0"


Spatial interpolation with breaklines in R

I have a dataset consisting of irregularly-spaced points (x,y,value) that represent an area of ocean that contains some land. I need to interpolate these to a regular grid, ensuring that the interpolation routine does not consider data points that are on the opposite site of a land barrier (e.g. the other side of a peninsular). I understand that this means that I need to insert "breaklines" in the data, and I have polygons of the coastline that I can use for this. Do any R packages allow for spatial interpolation using breaklines? Ideally I think that natural neighbour interpolation would be the most appropriate for my dataset, but I suspect that the breaklines are more important than the interpolation method.


Contouring surface with barriers?

I'm using ArcGIS Desktop 9.3, and I want to contour a surface with barriers. Every time the result is a polyline contour shape file which covers only a part of the surface not the whole surface. However, when I use "Contour" without barriers it works fine and the resulting contour shapefile covers the whole surface. Can anybody help with this?


QGIS documentation in French?

I am working on a plugin in qgis and i have to make a document from A to Z about that but on French . and I am looking now for some documentation of python,QT , pyqt4, and QGIS in French.


Repaint the map after updating a feature in QGIS

I update a feature in a vector layer using the following Python code:               layers = QgsMapLayerRegistry.instance().mapLayers()          layer = layers.get('my_vector20140614125509945')          fid = 1     attrs = { 0 : "hello", 1 : 100 }     layer.dataProvider().changeAttributeValues({ fid : attrs })      The feature is updated but the map is not redrawn. I am using catagorized style. The value of 100 will be shown with a different color.


Writing values to feature class in Python script results in strange characters"?

The script below does a series of things, and towards the end it is supposed to write values pulled from a list to a new field, if the feature class name is found anywhere in a list element. It does write values but there are strange characters that appear, or it doesn't write the entire list element value. I can't tell what the problem is?               import os, shutil,arcpy,arcgisscripting     gp = arcgisscripting.create()     path = r'C:\Temp\FME_share'     copy_dir = r'C:\test'     gdb_dir = copy_dir+'\\'+'CAD_Export_to_GDB'     BB_gdb_path = gdb_dir+'\\'+'BoundBoxes.gdb'     file_ext_lst = ['.2dm','.2de','.3dm','.3de','.dgn']     ext_to_chng = ['.2dm','.2de','.3dm','.3de']     file_lst = []     TxtVal_lst = []     repr_lst = []     found = dict()     arcpy.env.overwriteOutput = 'TRUE'               for dirpath, dirnames, filenames in os.walk(path):       for filename in filenames:           fullPath = os.path.join(dirpath, filename)           name,ext = os.path.splitext(filename)           if ext not in file_ext_lst:               continue           if name not in found:               found[name] = fullPath          print 'creating unique list of file names...'     for name,path in found.iteritems():       file_lst.append(path)          print 'copying files to local drive...'     for f in file_lst:       file_nm = f.split('\\')[-1]       shutil.copyfile(f, copy_dir+'\\'+file_nm)          print 'changing file extensions of copied files to .dgn...'     for dirpath, dirnames, copiedfiles in os.walk(copy_dir):       for copyfile in copiedfiles:           c_name,c_ext = os.path.splitext(copyfile)           filename_new = copy_dir+'\\'+c_name+'.dgn'           if c_ext in ext_to_chng:               os.rename(copy_dir+'\\'+copyfile, filename_new)          print 'exporting dgn files to geodatabase and exploding complex features...'     for dirpath, dirnames, renamed_files in os.walk(copy_dir):       for rn_file in renamed_files:           rn_name, rn_ext = os.path.splitext(rn_file)           in_dgn = dirpath+'\\'+rn_file           out_gdb = gdb_dir+'\\'+rn_name+'.gdb'           if rn_ext == '.dgn':               gp.ImportCAD_conversion(in_dgn,out_gdb,'','Explode_Complex')          print 'renaming geodatabase output to include original dgn file names...'     for dirpath, dirnames, filenames in arcpy.da.Walk(gdb_dir):       for filename in filenames:           gdb_name = dirpath.split('\\')[3]           orig_dgn_nm = gdb_name.split('.')[0]           if not dirpath.endswith('.gdb'):               arcpy.Rename_management(dirpath+'\\'+filename,dirpath+'\\'+filename+'_'+orig_dgn_nm)          print 'joining dgn text attributes to point feature classes, then adding text values to a list for searching out coordinate system information...'     for dirpath, dirnames, filenames in arcpy.da.Walk(gdb_dir):       for filename in filenames:           fullPath = dirpath+'\\'+filename           if 'Point' in filename:               Point_FC = dirpath+'\\'+filename             arcpy.JoinField_management(Point_FC,'EntID',dirpath.replace('CADStaging','')+'TxtProp','EntID','TxtValue')             values = [row[0] for row in arcpy.da.SearchCursor(Point_FC,'TxtValue')]             unique_vals = set(values)             unique_vals_str = ' '.join(str(uv) for uv in unique_vals)                  if ('CCS 83' or 'CCS83') and 'ZONE 3' and 'METERS' in unique_vals_str:                 print 'Trying to assign coordinate system to',dirpath,'....'                 arcpy.DefineProjection_management(dirpath,26943)             elif ('CCS 83' or 'CCS83') and 'ZONE 3' and 'FEET' in unique_vals_str:                 print 'Trying to assign coordinate system to',dirpath,'....'                 arcpy.DefineProjection_management(dirpath,2227)             elif ('CCS27' or'CCS 27') and 'ZONE 3' in unique_vals_str:                 print 'Trying to assign coordinate system to',dirpath,'....'                 arcpy.DefineProjection_management(dirpath,26743)             elif ('CCS 83' or 'CCS83') and 'ZONE 2' and 'METERS' in unique_vals_str:                 print 'Trying to assign coordinate system to',dirpath,'....'                 arcpy.DefineProjection_management(dirpath,26942)             elif ('CCS 83' or 'CCS83') and 'ZONE 2' and 'FEET' in unique_vals_str:                 print 'Trying to assign coordinate system to',dirpath,'....'                 arcpy.DefineProjection_management(dirpath,2226)             elif ('CCS27' or 'CCS 27') and 'ZONE 2' in unique_vals_str:                 print 'Trying to assign coordinate system to',dirpath,'....'                 arcpy.DefineProjection_management(dirpath,26742)             else:                 print 'Couldn\'t find coordinate system'          print 'creating minimum bounding box for line feature classes...'     for dirpath, dirnames, filenames in arcpy.da.Walk(gdb_dir):       for filename in filenames:           if 'Line' in filename:               BB_in_fc = dirpath+'\\'+filename               BB_out_fc = BB_gdb_path+'\\'+'BoundBox'+filename.replace('Line','')               arcpy.MinimumBoundingGeometry_management(BB_in_fc,BB_out_fc,'RECTANGLE_BY_AREA','ALL')          print 'adding and calculating SHEET_NAME and SHEET_PATH fields to bounding box feature classes, then reprojecting to WGS84...'     for dirpath, dirnames, filenames in arcpy.da.Walk(BB_gdb_path):       for filename in filenames:           sheet = filename.replace('BoundBox_','')           BB_fullpath = dirpath+'\\'+filename           arcpy.AddField_management(BB_fullpath,'SHEET_NAME','STRING')         arcpy.CalculateField_management(BB_fullpath,'SHEET_NAME',"'"+sheet+"'","PYTHON_9.3","#")         arcpy.AddField_management(BB_fullpath,'SHEET_PATH','STRING',200)         for f in file_lst:             if sheet in f:                 f_str = str(f)                 print f_str                 arcpy.CalculateField_management(BB_fullpath,'SHEET_PATH','"'+f_str+'"',"PYTHON_9.3","#")         BB_desc = arcpy.Describe(BB_fullpath)         BB_CS = BB_desc.spatialReference         BB_repr_file = BB_gdb_path+'\\'+filename+'_repr'         BB_repr_CS = arcpy.CreateSpatialReference_management('WGS 1984')         if '1927' in BB_CS.name:             arcpy.Project_management(BB_fullpath,BB_repr_file,BB_repr_CS,'NAD_1927_To_WGS_1984_79_CONUS')         else:             arcpy.Project_management(BB_fullpath,BB_repr_file,BB_repr_CS,'NAD_1983_To_WGS_1984_1')         repr_lst.append(BB_repr_file)     merge_fc = BB_gdb_path+'\\'+'BoundBoxes_All'     print 'finally, merging all bounding boxes into one feature class!'     arcpy.Merge_management(repr_lst,merge_fc)


Where can i get Shape Files for AE - U.S. Armed Forces – Europe AP - U.S. Armed Forces – Pacific PR - Puerto Rico?

I need shape files for these : AE - U.S. Armed Forces – Europe AP - U.S. Armed Forces – Pacific PR - Puerto Rico Where can i find them ?


Need help with GeoExpress 9 by LizardTech/ArcGIS Server 10.1

This is a multi-part question and after much reading on different forums I still could use some help/advice. My setup: ArcGIS Desktop 10.1, Dedicated Server for GIS data with SQL Server 2008 R2 and ArcServer 10.1. My problem: I have switched from ArcServer 9.3 to 10.1 which no longer supports ECW format for free. We can't afford the ESRI Image Extension or license for ECW. I need a way to convert my 5 years of Aerial Imagery in ECW mosaics to another format that is compatible with ArcServer, not too large (less than 30GB for each year's imagery). These are County wide datasets so pretty large when not compressed (~250GB for each year of Imagery). I also have individual tiles for each in MrSid, JPEG and Tiff format. I have not found a solution using ArcGIS tools so that the datasets are small, fast over the network/ArcServer, and "seamless". I've recently tried a free trial of LizardTech's GeoExpress 9 to create one MrSid mosaic of all the smaller MrSid files. This would have worked great and is only 24GB with minimal loss in quality, the problem is that the entire mosaic is covered in small semi-transparent hash marks (X X X X) that look like this... Individual tiles: ![Individual Tiles](http://i.stack.imgur.com/y5Q0N.jpg) The GeoExpress 9 (Trial Version) Mosaic: ![GeoExpress 9 version](http://i.stack.imgur.com/UaEgs.jpg) Has anyone seen this before? This is MrSid Gen4 could Gen3 help with ESRI compatibility? Could this be due to only having a trial version of LizardTechs software? I have not received an answer from them yet. Or if someone has used this software maybe I just did something foolish and used a wrong setting? The best solution for me, since I also use the ESRI JavaScript API for our WebApps, would be to use ArcGIS software to create some sort of image or mosaic of each year's imagery. One more option could be an OpenSource solution that supports ECW format? I could bring that service into my Web Apps and use ArcServer for all other layers, but I have not found one yet. I have created raster catalogs, merged the tiles into one raster, etc. with ArcGIS but the original tiles boundaries show through due to slight color mismatching. I'd like my end result to be one "seamless" image. I'm new to imagery processing so any help is much appreciated. Thanks in advance!


How to create Series of Rasters that displays Standard Deviation from Mean?

Cell statistics question for you. I have 100 rasters of raw rainfall data in which the pixels each represent rainfall for that particular pixel on the landscape for a single year during the 100 year period. I have also created a raster that calculated the mean rainfall received in each pixel over those 100 years. From these raw data rasters, I want to create a new set of 100 rasters (one for each year) that displays the standard deviation from the mean of each pixel (of the entire 100 years) for each pixel. I have a feeling this will be likely solved with python coding, which I have a little bit of experience in, but I am unsure of what tool to use? It seems like Cell Statistics is on the right track, since it focuses on individual cells, but doesn't achieve what I am looking for. Any other ideas on how to tackle this?


I am trying to find PDF documentation for Maplex for ArcGIS

Does anyone know where to find the PDF documentation for Maplex for ArcGIS 9.3 ? There was a website from ESRI but it was broken link ? I need to review it because I am working on impaired streams and lakes and needs their labels to set it up properly. Trying to do it in Label Manager several times. Thanks,


How to choose workspace when uploading SLD to Geoserver with the REST API

I have 200+ SLD's I would like to upload to Geoserver using the REST API. I successfully uploaded an SLD using the gsconfig Python module. However, I would like to be able to set the workspace when I upload. Is this possible? Below is my script so far.               from geoserver.catalog import Catalog          cat = Catalog("http://chchchch:8080/geoserver/rest", "username", "pw")          with open("C:\Users\Documents\SLD\my_sld.sld") as f:         cat.create_style("my_sld", f.read())


Should I use a CMS for my web mapping application?

For my web mapping app, I am using GeoServer as server, PostGIS as spatial DB and OpenLayers on client side (as well as other JS libraries). The map is composed of base maps and vector layers. One of these vector layers is the "Public institutions" layer. For each category of institutions (Ministries, Tribunals...) there are persons _(I called them "contributors")_ responsible for editing data (Tel, fax, location on map...) about these institutions. The contributors are managed by an admin who specify the privileges of each one of them. Some modifications of the contributors need to be validated by the admin (for example deleting an institution)... So I thought maybe using a CMS would help. What do you think? There are some CMS like Drupal that integrates OpenLayers, but it seems the project is still in development link to project page. Any alternative or experience to share?


How to get estimates of areas of overlapping polygons in one shp to another?

Sorry for my ineptitude, but I have two vector shape files of Africa. One has polygons of countries, and one has polygons of tribes a long time ago. What I want to do is to find the intersection of these two polygons, and get some idea of how many tribes were contained in a country today, and what share of land in the country today that they took up. Hopefully that makes sense. So, I'd like to end up with some data with a row for a country, say, Nigeria, that has: NIGERIA | Tribe 1: 20% of land | Tribe 2: 40% of land | Tribe 3: 40% of land | Something along those lines. Is this possible in QGIS or ArcGIS? My preferred software is QGIS.


GDAL.RasterizeLayer not ported to the C# bindings?

I am trying to use the `OSGeo.GDAL.Gdal.RasterizeLayer(…)` from _GDAL 1.9.2_ and it requires some `SWIGTYPE` parameters:               OSGeo.GDAL.Gdal.RasterizeLayer(         OSGeo.GDAL.Dataset dataset,         int bands,         OSGeo.GDAL.SWIGTYPE_p_int band_list,         OSGeo.OGR.Layer layer,         OSGeo.GDAL.SWIGTYPE_p_void pfnTransformer,         OSGeo.GDAL.SWIGTYPE_p_void pTransformArg,         int burn_values,         OSGeo.GDAL.SWIGTYPE_p_double burn_values_list,         string[] options,         OSGeo.GDAL.Gdal.GDALProgressFuncDelegate callback,         string callback_data     );      Reading up on some old mails somewhere ( _I've lost the link_ ) regarding ReadRaster and WriteRaster in the C# bindings, **Tamas Szekeres** states that methods requesting `SWIGTYPE` parameters, cannot be used (yet). **Question:**   Does this mean that the RasterizeLayer method hasn’t been ported to the C# bindings or are there some clever way to utilise this method?


GPS problem in QGIS

I have a very strange problem that I cannot understand. I am trying to put together a demo-laptop with QGIS and gps. On other computers I have it working as expected, but on this particular box something is strange. What I get is a gps position about 680 meters south east of my position. The position is fixed in QGIS, it is not "walking around" as expected. I see the altitude is changing though, but not lon lat. I have tried both connecting directly through /dev/rfcomm0 and through gpsd, same result. In other clients everything looks ok, like Foxtrot GPS and gpsmon. In gpsmon, and Foxtrot gps I get another position than in QGIS and the position in the other clients is moving around. I have moved the laptop and tried at another location and I get the same behavior but with another position some hundred meters away. I did get desperate and re installed the laptop, but get the same behavior. I am totally confused. It cannot be a hardware problem I think since it is working on other clients in the same laptop. It shouldn't be a software problem either since I have re installed everything including the OS. As mentioned I get it working on other computers with the same software. My setup is: OS Ubuntu 11.04 QGIS compiled from trunk a few days ago I can think of two reasons, but don't know how to search: 1) Something in my hardware gives another configuration when I am compiling QGIS than on the other computors 2) I have some software-package present (that I apparently also re installed in that case) that is interfering. Any ideas is highly appreciated Thanks   Nicklas


Issues with arcpy.sa.Reclassify

I have a raster file with classifcation codes, and want to reclassify valid codes to 1, and all others to 'NoData'.               RMAP2 = arcpy.sa.RemapValue([[3.2, 1],                                  [4.1, 1],                                  [4.2, 1]])          LCA = arcpy.sa.Reclassify(LCA, "LCCODE_USE", RMAP2, "NODATA")      The output however only contains the values 0 and 127, with the 0 value coinciding with just one of the missing values. Any idea what the hell is going on?