title: Concatenate multiple string / text fields some of which may have NULL values ArcGIS

text: I am trying to repeat the behavior of this code from the ESRI website. I want to replicate the output example in `ALL_Type` field, so that it concatenates a field but it deals (ignores) any that have NULL values. Unfortunately this example no longer works at 10.2.2 because all example of NULL return a "None" text. Instead I just want a blank to show. ![enter image description here](http://i.stack.imgur.com/VjHgM.jpg) ![enter image description here](http://i.stack.imgur.com/hZZRY.jpg) And here is the code.                  "*args" allows this routine to accept any number of field values.     # the values are passed as a Python tuple, essentially a     # non-editable list     #     def concat(*args):            # Initialize the return value to an empty string,       # then set the separator character       #       retval = ""       sep = "_"            # For each value passed in...       #       for t in args:         # Convert to a string (this catches any numbers),         # then remove leading and trailing blanks         #         s = str(t).strip()              # Add the field value to the return value, using the separator         # defined above         #         if s <> '':           retval += sep + s            # Strip of any leading separators before returning the value       #       return retval.lstrip(sep)


title: How to draw a rectangle graphic element (given feature and extents in meters) with ArcObjects?

text: I am trying to render a rectangle graphics box onto the map, given the center location of an polygonal feature and the extent dimensions (length and width) in meters. The feature has the spatial reference of its shape set to WGS84 though. I am able to get the location of the feature (an IPoint) by taking its centroid. Is there then a straightforward way to accomplish the drawing of such a box with meter-dimensions around that point? Looking at the API, IEnvelope seems to be what I need, but will I have issues with unit-conversions (since I need the extent dimensions to be in meters)? And how should I go about rendering the IEnvelope? Any help is appreciated, thanks!


title: Need help for getting a high resolution satellite image for digitising features in ArcGIS?

text: I want to digitise some places of my city in ArcGIS. I found ArcGIS 'World Imagery' layer which I can use it as a basemap for digitisation. But the problem is it works only while the internet is connected. Is there anyway to save the cache of the layer at a particular zoom level so that I can work on it offline even after I restart the computer. If no, please tell me any other way- may be with google or bing imagery, I need to digitise important features of the city. So i require resolution at house level.


title: How to detect(condition) if Query Layer Connection is failing?

text: I have a job (arcpy script) which relies on a query layer but it fails sometimes. It has to do with the query layer connection (the point where job fails randomly)- which is totally under the control of ORACLE team. I want to capture this to discuss with ORACLE team. What should be the right layer property to capture this?


title: How to use ArcGIS 10.1 to make Python scripts for ArcGIS 9.3?

text: Recently I have been making scripts to run geoprocessing tools for a client who is running ArcGIS 9.3. However I only have access to ArcGIS 10.1. As a result I cannot use model builder in 10.1, and have to script everything with arcgisscripting (`import arcgisscripting`). I cannot find many arcgisscripting examples beyond the 9.3 help files. My questions are: is there a way to reverse convert 10.1 (arcpy/modelbuilder) scripts to 9.3? Is there a place with good arcgisscripting .py examples? Has anyone else done this?


title: How to draw parallel or curved labels for OSM data in QGIS?

text: I'm trying to recreate a Google style map using OpenStreetMap data and QGIS. I've been following underdarks many wonderful guides but am getting stuck at the labels. The OSM data was taken from geofabrik and imported into a postgreSQL/postgis database using osm2pgsql. When I try to label the roads it seems to place them as points and does not give the option for parallel or curved labeling. You can see the result below. ![enter image description here](http://i.stack.imgur.com/VQUaz.jpg) When I click on the advanced tab for labeling all I see are the data defined placement options and no check box for curved or parallel labeling. I'm new to QGIS so I hope I'm missing something really simple.


title: Hydrograph in GRASS

text: How to make hydrograph for a basin in GRASS? While surfing the GRASS home page, I came to know that some time back in GRASS 4.x a command _r.water.fee_ exist. I understand, it is missing in GRASS 6.4. Is there any equivalent command for this? Any alternate? Thanks in advance....


title: Change Attirbute Name (Fiona)

text: I have a shapefile that has 2 attributes named "X" and "Y". I want to change this to "Lon" and "Lat". My (limited) experience tells me that I should read in the original shapefile, preserve the original driver, crs, make a new schema with "Lon" and "Lat" defined, and then write a new/new named shapefile. Is there a quicker/easier way with fewer steps?               import fiona          with fiona.open(filein, 'r') as source:         source_driver=source.driver         source_crs=source.crs         source_schema=source.schema         print(source_schema)          new_schema= {'geometry': 'Point',              'properties': { 'Lon': 'float:16.6',                               'Lat': 'float:16.6',                               'SP': 'float:16.6',                               'X_meter': 'float:19.11',                               'Y_meter': 'float:19.11',                               'Depth': 'float:19.8',                               'Line_Name': 'str:50',                               'SurveyID': 'str:50'}}          with fiona.open(fileout, 'w',             driver=source_driver,             crs=source_crs,             schema=new_schema) as c:         print(len(c))          c.closed


title: How to rasterize a point layer? Cell value should be 1 if feature attribute is 1

text: I’m a fairly experienced GIS user, but I’m new to QGIS. I’m using version 2.2, and I’m trying to convert a point data layer into a raster data layer using the Rasterize tool. The point data layer has an integer field with a value of 1 for all records, and I want to make a raster data layer which has a value of one for each cell which has at least one point in it (some cells have more than one point in them). When I use the Rasterize tool, it creates the raster data layer fine, but all cells with have points in them appear to have a value of 0.995968. I’ve tried editing the gdal code (even making the resulting raster data layer an integer layer), but I still cannot get it to produce a raster with a value of 1 rather than just a number very close to one. Does anyone have any suggestions of where I might be going wrong and/or how to solve this issue? On a related note, is it possible to use the Rasterize tool to work out the average or total values for all the points which fall in each grid cell? If so, how (I can’t find anything in the gdal instructions themselves which would suggest this is possible, but I might be missing something). If this isn’t possible to do with the Rasterize tool, are there other tools in QGIS which can do this? Thanks for any help anyone can give on this. All the best, Colin


title: pass an array of geometry to postgresql function

text: I am trying to pass an array of geometry to a postgis function using JDBC with JAVA:                           CallableStatement proc = conn.prepareCall("{ ? = call somefunc(?) }");                 proc.registerOutParameter(1, 7);                 PGgeometry[] arr = new PGgeometry[2];                 arr[0] = a;                 arr[1] = b;                 proc.setArray(2, conn.createArrayOf("geometry", arr));                 proc.execute();      a,b are valid PGgeometry types. I am getting this exception:               Exception in thread "main" java.lang.AbstractMethodError: org.postgresql.jdbc3g.Jdbc3gConnection.createArrayOf(Ljava/lang/String;[Ljava/lang/Object;)Ljava/sql/Array;      Is this even possible to do? I found no good answer for that. Thanks!


title: Carto CSS variables

text: Does anyone know whether it's possible to use a variable within an attribute selector. For example I define my variable at the top of the file:               @lowZoomMax: 9;      And then reference it as follows:               #settlement {      [zoom > @lowZoomMax] {     ...     ...     }     }      This doesn't appear to work in the map I currently have - I get `Invalid Code: @lowZoomMax` (in the second bit of code, not the definition).


title: How can I get my raster analysis to produce a larger extent in ArcGIS

text: Ok so the situation is as follows. I have an vector shapefile and within that area there are three points. I converted the points into raster and now I want to do a Euclidean distance analysis in ArcMap. My window is clipped to the extend of the vector shapefile so that the Euclidean distance analysis will only cover the inside of my area. Instead getting a nice output the output only covers a certain part of the vector shapefile instead of the full extent. so this is the map with the point. ![enter image description here](http://i.stack.imgur.com/0dxr5.png) and this is afer the analysis ![enter image description here](http://i.stack.imgur.com/7DI1T.png) Does anybody know why this is and how I can fix it?


title: How to reproject EPSG:102218 with QGIS?

text: I have a raster file. ArcMap says it is in `NAD_1983_USFS_R6_Albers` (EPSG:102218). QGIS says it is in `USER:100000 - * Generated CRS (+proj=aea +lat_1=43 +lat_2=48 +lat_0=34 +lon_0=-120 +x_0=600000 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs)`. I assume QGIS doesn't know that coordinate system. So I can also not re- project it with QGIS, right? Is there a work around to still reproject it with QGIS? I don't want to re-project it with ArcMap.


title: Materialized View Does Not Appear In Geometry_Columns View

text: I have defined a materialized view in a PostGIS 2.1-enabled database running PostgreSQL 9.3.1:               CREATE MATERIALIZED VIEW canvec_data.contours_m AS       SELECT         contours_original.ogc_fid,          contours_original.elevation::integer AS elevation_orig,          CASE             WHEN "substring"(contours_original.code::text, 1, 3) = '257'::text THEN round(contours_original.elevation::double precision * 0.3048::double precision)             ELSE contours_original.elevation::double precision         END AS elevation_m,          contours_original.wkb_geometry::geometry(Linestring,3578)      FROM canvec_data.contours_original      I had, perhaps erroneously, expected that this resultant entity (a table?) would be listed in the `public.geometry_columns` view, and thus be properly interpreted by GIS viewer software such as QGIS. The materialized view does not appear in the `public.geometry_columns` view. What am I doing wrong? Thanks!


I am migrating my Application built in Version 9.3.1 to version 10 In my code I am getting the exception: `Retrieving the COM class factory for component with CLSID {4147A362-EFC4-11D3- A0A4-00C04F6BC626} failed due to the following error: 80040111` at this place: `IExportOperation pExportOp = new ExportOperationClass();` I have already bound my application with bindmanager code: `ESRI.ArcGIS.RuntimeManager.Bind(ESRI.ArcGIS.ProductCode.Engine);` Every thing works fine like creating IWorkspaceFactory InitiallizeLicense, the exception arises as this statment is executed: `IExportOperation pExportOp = new ExportOperationClass();` Is there any solution to avoid this exception? Please help. Thanks


title: Unicode Characters in GRASS ps.map instructions file

text: I am trying to put text on a ps.map generated map and the text contains the degrees symbol. I can insert the symbol into the instructions file, and it looks fine, but when I generate the map, the result is Â° Here's the instruction: text 99% 63% 60 °N fontsize 16 ref left end Does anyone have advice on how to get unicode symbols to display correctly with ps.map?


title: why are my buffers not lining up with their original points despite setting the layer crs as all the same crs when I reopen my project

text: I am having my entire project in pseduo mercator (google earth raster image) and when I create buffers for points and then open back up my project all the buffers are offset from their points. Update: pseudo-mercator-error So my buffers I got to line up with their original centers, but other points that I meant to place within proximity were offset. I think it has something to do with meters and relative buffer distance being the reason why the buffers were fixed, but not some of the points.


title: OpenLayers Icon e PopUp

text: Ciao, sto creando una pagina in php con una mappa ... in OpenLayers Examples ho trovato due interessanti esempi:   * esempio che mostra come usare diversi tipi di poupup (http://openlayers.org/dev/examples/popupMatrix.html)    * esempio per inserire icone personalizzate (http://openlayers.org/dev/examples/markers.html)  Il mio problema è che non riesco ad integrare i due esempi... dal primo:               echo"              ll = new OpenLayers.LonLat($lng,$lat).transform(                  new OpenLayers.Projection('EPSG:4326'), new OpenLayers.Projection('EPSG:900913')              );              popupClass = OpenLayers.Popup.Anchored;              popupContentHTML = '$T';              addMarker(ll, popupClass, popupContentHTML, false, true);       dal secondo:               echo"              ll = new OpenLayers.LonLat($lng,$lat).transform(                  new OpenLayers.Projection('EPSG:4326'), new OpenLayers.Projection('EPSG:900913')              );              var size = new OpenLayers.Size(50,50);              var offset = new OpenLayers.Pixel(-(size.w/2), -size.h);              var icon = new OpenLayers.Icon('{$row['enclosure']}',size,offset);              markers.addMarker(                  new OpenLayers.Marker(                      (ll                      ),icon                  )              );      ";       Ho notato che nel primo addMarker richiama una funzione, che sembra simile a markers.addMarker ... sinceramente non so dove mettere ", icon" nel primo o nella funzione che richiama? EDIT - Translation to English (Tradotto in inglese) Hello, I am creating a php page with a map ... in OpenLayers Examples I found two interesting examples: • Example showing how to use different types of poupup (http://openlayers.org/dev/examples/popupMatrix.html) • example to insert custom icons (http://openlayers.org/dev/examples/markers.html) My problem is that I can not integrate the two examples ... the first:               echo"                  ll = new OpenLayers.LonLat($lng,$lat).transform(                      new OpenLayers.Projection('EPSG:4326'), new OpenLayers.Projection('EPSG:900913')                  );                  popupClass = OpenLayers.Popup.Anchored;                  popupContentHTML = '$T';                  addMarker(ll, popupClass, popupContentHTML, false, true);       The Second:               echo"              ll = new OpenLayers.LonLat($lng,$lat).transform(                  new OpenLayers.Projection('EPSG:4326'), new OpenLayers.Projection('EPSG:900913')              );              var size = new OpenLayers.Size(50,50);              var offset = new OpenLayers.Pixel(-(size.w/2), -size.h);              var icon = new OpenLayers.Icon('{$row['enclosure']}',size,offset);              markers.addMarker(                  new OpenLayers.Marker(                      (ll                      ),icon                  )              );      ";       I noticed that in the first addMarker calls a function, which looks similar to markers.addMarker ... I honestly do not know where to put "icon" in the first or in the function calling?


title: Python Script to convert a mass of KMLs to a single fGDB

text: I have a mass of KML files, approx 350+ that I need in a single fGDB. I lifted a string of code from ESRI help to attempt to do this - I got it working to the point of creating all the individual gGDBs (one per KML) but its failing with an invalid character in a variable somewhere. Below is the code, and then the error that it throws out. I've added some print commands to show whats in the variables right before it fails. Help? Also, this is my first time trying to use python for actual work, so I apologize in advance for n00bish mistakes.               # Name: BatchKML_to_GDB.py     # Description: Converts a directory of KMLs and copies the output into a single fGDB.     #              A 2 step process: first convert the KML files, and then copy the featureclases          # Import system models     import arcpy, os          # Set workspace (where all the KMLs are)      arcpy.env.workspace="C:\\Users\\mperez\\Desktop\\IT-Paths"          # Set local variables and location for the consolidated file geodatabase     outLocation = "V:\\2155\\active\\215510456\\GDB"     MasterGDB = r'AllKMLLayers.gdb'     MasterGDBLocation = os.path.join(outLocation, MasterGDB)          # Create the master FileGeodatabase      arcpy.CreateFileGDB_management(outLocation, MasterGDB)          # Convert all KMZ and KML files found in the current workspace      for kmz in arcpy.ListFiles('*.KM*'):       print "CONVERTING: " + os.path.join(arcpy.env.workspace,kmz)       arcpy.KMLToLayer_conversion(kmz, outLocation)               # Change the workspace to fGDB location     arcpy.env.workspace = outLocation          # Loop through all the FileGeodatabases within the workspace     wks = arcpy.ListWorkspaces('*', 'FileGDB')     # Skip the Master GDB     wks.remove(MasterGDBLocation)          for fgdb in wks:              # Change the workspace to the current FileGeodatabase       arcpy.env.workspace = fgdb                # For every Featureclass inside, copy it to the Master and use the name from the original fGDB       featureClasses = arcpy.ListFeatureClasses('*', '', 'Placemarks')       for fc in featureClasses:         print r"COPYING: " + fc + r" FROM: " + fgdb          fcCopy = fgdb + os.sep + 'Placemarks' + os.sep + fc         print r"fcCopy = " + fcCopy         print r"MasterGDBLocation = " + MasterGDBLocation         print r"fgdb = " + fgdb         arcpy.FeatureClassToFeatureClass_conversion(fcCopy, MasterGDBLocation, fgdb[fgdb.rfind(os.sep)+1:-4])               # Clean up     del kmz, wks, fc, featureClasses, fgdb      Output from shell showing error:               COPYING: Points FROM: V:\2155\active\215510456\GDB\canal-move-2-golden-gate-blvd.gdb     fcCopy = V:\2155\active\215510456\GDB\canal-move-2-golden-gate-blvd.gdb\Placemarks\Points     MasterGDBLocation = V:\2155\active\215510456\GDB\AllKMLLayers.gdb     fgdb = V:\2155\active\215510456\GDB\canal-move-2-golden-gate-blvd.gdb          Traceback (most recent call last):       File "C:\Users\mperez\Desktop\batch_kml_to_gdb.py", line 46, in <module>         arcpy.FeatureClassToFeatureClass_conversion(fcCopy, MasterGDBLocation, fgdb[fgdb.rfind(os.sep)+1:-4])       File "C:\Program Files (x86)\ArcGIS\Desktop10.1\arcpy\arcpy\conversion.py", line 1547, in FeatureClassToFeatureClass         raise e     ExecuteError: ERROR 000354: The name contains invalid characters     Failed to execute (FeatureClassToFeatureClass).


title: Watershed Delineation

text: I am trying to delineate just the catchment area of some lakes, and am wondering the correct way to do this. I already have a filled DEM, along with a flow direction and a flow accumulation raster. Is there a way to accomplish this task without the use of pour points? Or do I have to use them. If that is the case, can a pour point be located in the middle of a lake (assuming that the middle of a lake is the highest area of flow accumulation). Whenever I do try to run this data, whether with or without a pour point, I get something that does not look completely right. The catchment area does not encompass the entire lake, and the actual watershed itself does not accurately represent the flow accumulation arrows. Am I doing something wrong, or is there a simple cause to my problem. Maybe if someone could take a quick second to step through the watershed toolset, that would be helpful. Am using ArcGIS 10 SP1, with Spatial Analyst, ArcView license.


title: Combining several rasters before r.surf.contour?

text: I have been trying to generate an elevation model from contour lines using first _v.to.rast.attr_ and then _r.surf.contour_. Generally, it works pretty well, execpt for one issue. My data contains a lot of sea without contours, and _r.surf.contour_ spends a lot of time in the sea region where there are no pixels in the raster generated from _r.to.rast.attr_. Actually, I'm running one now on a reasobly large area, and after 24 hours, it is at 70%. The solution for me would be simple. I have sea polygons (and even lake polygons with elevation). I could generate a separatr raster for that using _v.to.rast.attr_. But then I need to combine the two rasters into a single raster before running _r.surf.contour_. How do I do that?


title: Change color of shape in Layer File

text: I currently have a layer file - in this case one of the United States. This layer contains a shape for each state so if I were to open the TOC Control and expand the "United States" I would see Alabama, Alaska, etc etc. Currently, if I wanted to I could change the color of all the states through the following code:               IRgbColor rgb = new RgbColorClass();     ColorDialog myColorDlg = new ColorDialog(); ///This is the windows color picker          DialogResult res = myColorDlg.ShowDialog();          if(res == DialogResult.OK)     {        IFeatureLayer featureLayer = layer as IFeatureLayer;        IFeatureClass featureClass = featureLayer.FeatureClass;            rgb.Red = myColorDlg.Red;       rgb.Blue = myColorDlg.Blue;       rgb.Green = myColorDlg.Green;            IGeoFeatureLayer geoFeatureLayer = layer as IGeoFeatureLayer;            ILineSymbol lineSymbol = new SimpleLineSymbolClass();       lineSymbol.Width = 1.0;       lineSymbol.Color = (IColor)rgb;            ISimpleFillSymbol fillSymbol = new SimpleFillSymbolClass();       fillSymbol.Outline = lineSymbol;       fillSymbol.Color = (IColor)rgb;            ISimpleRenderer simpleRenderer = new SimpleRendererClass();       simpleRenderer.Symbol = (ISymbol)fillSymbol;       geoFeatureRenderer.Renderer =  (IFeatureRenderer)simpleRenderer;            myAxMapControl.Refresh();     }      The question now is, how can I change each state separately. I'm aware that if I use the built in TOC control I should be able to do this with ILegendClass; however, I wish to use a different TOC Control.


title: Programmatically uncheck layer in AxTOCControl

text: I am running my map in dynamic display mode, and I want to programmatically uncheck a layer that is not a IDynamicLayer that is visible in the AxTocControl -- its just a regular layer file. I thought that `layer.Visible = false;` with a map refresh would work, but bafflingly it does not. I'm using ArcEngine 10 VS2010 & C#


title: How to add data to relationship table in ArcGIS 10.1 database design?

text: I have created 2 tables HORSE and VENDOR and a relationship table VENDOR_HORSE   With the foreign key in Horse being the Primary key in VENDOR (VendorID)   and the foreign key in VENDOR being the primary key in HORSE (HorseID).   It is a many to many relationship and ArcGIS automatically creates the link table -.   * yet I had to create the Relationship table which is simple and has extra fields (besides the horseID and vendorID) as Notes (text), Service (coded domain),  and Date (most important to solve many to many matter and have a relational database.   I cannot add data to the VENDOR_HORSE table.   I click on start an edit session, just as I do to add rows to the HORSE and VENDOR tables, with no problems, yet I cannot add rows to the Relationship table.   ESRI help gave instructions for adding rows yet they do not work.   That is,   1. I selected a row in HORSE,    2. I selected a row in VENDOR,    3. then double clicked on OBJECTID to have the table fields appear in the Attribute window and NOTHING. What am I missing? Help! Thanks Bette


title: Finding Normal to a given point in ECEF coordinates

text: I have a point of reference called Ground Reference Point (GRP) given in Earth Centred Earth Fixed (ECEF) coordinates as given below: (1.4892e+6, 6.1102e+6, 1.0637e+6) (all distances in meters) I want to find a normal to the earth ellipsoid at the above mentioned GRP again in ECEF coordinates. What is the process to calculate it?


title: What is a Multilinestring?

text: What is multiline string in gis? I have found when I use pgrouting certain edges have this geometric property mentioned as multiline string. I am not sure what's the purpose of this geometry and how it actually defines an edge. Any insights?


title: Correct way to create a DEM from point data in csv files

text: I know this might seem a somewhat general question, but I have to start somewhere. I have a lot of csv files which contain point data containing height, aspect and slope measurements at a resolution of 12,5 meters. I need to create a dem from this data and have ArcGIS 10 with spatial analyst and 3d Analyst. I did this before with 9.3 but I am not sure if my methodology is correct.   1. Create point shapefiles from the csv file   2. interpolate to raster (here I never know whether to use spline, kriging or some other method)? Last time I did something wrong and rounded all the values up to 1m even though the data was to cm accuracy so this time I want to do it right! If anyone has any tips or advice, please write, Thanks, Rob


title: Having trouble getting specific items from a list. What am I doing wrong?

text: I am a total beginner when it comes to python. I received some help already getting this code started but I hit a wall. I have a directory with 200 folders, each with a geodatabase, feature dataset, and point line and polygon feature classes. Each feature class has a field "Label" in the attribute table. I need to be able to run a search through each feature classes' attribute table and return a specific string value. For example, I need to find all feature classes where diabase is in the label field. The output would give me only the feature classes that contain diabase in the label field. I created a temp feature layer from each feature and did a select by attribute to see if get count returns anything. Here's my code:               import arcpy     arcpy.env.workspace = "X:"     ws = arcpy.ListWorkspaces("*","Folder")     arcpy.env.overwriteOutput = True     tempLayer = r"in_memory\templayer"     field = "Label"     value = "diabase"     where = "%s = '%s'" % (field,value)     dhList = []     for w in ws:          arcpy.env.workspace = w          gdb = arcpy.ListWorkspaces("*", "Access")          for fc in gdb:              arcpy.env.workspace = fc              fcl = arcpy.ListDatasets("*", "Feature")              for fcc in fcl:                  arcpy.env.workspace = fcc                  fccl = arcpy.ListFeatureClasses("*", "All")                  for stil in fccl:                      arcpy.MakeFeatureLayer_management(stil, tempLayer, where)                      arcpy.SelectLayerByAttribute_management(tempLayer,"SUBSET_SELECTION")                      count = str(arcpy.GetCount_management(tempLayer).getOutput)                      print count                      if count > 0:                         dhList.append(stil)     for item in dhList:          print item      The output is giving me a list of all of the feature classes in the directory instead of only the feature classes with diabase. What am I doing wrong? Thanks for any and all help in advance.


title: ArcSDE for ArcGIS Desktop

text: I have ArcGIS desktop (info) license. I want to know if I can configure SDE, where in I can work with versioning of geodatabases.


title: How to move Layer to back?

text: from my Plugin I've loaded a rasterlayer using the following code:               fileName = "/path/to/raster/file.tif"     fileInfo = QFileInfo(fileName)     baseName = fileInfo.baseName()     rlayer = QgsRasterLayer(fileName, baseName)     if rlayer.isValid():         QgsMapLayerRegistry.instance().addMapLayer(rlayer)         ...      Now I wonder how I can move the freshly loaded layer to the back, that it lies behind all previously loaded layers. Thanks and regards


title: How to align an HDF layer

text: Forgive me, but I am completely new to GIS. I am using QGIS. I have a raster layer (HDF format) that does not align with other layers. All layers are projected WGS84 (EPSG:4326) and still the problem persists. I know that the HDF layer is defined by rows and columns that need to be converted to longitudes and latitudes but I don't know how to do that.   Also, the HDF file, for some reason, can't be exported as CSV file. ![enter image description here](http://i.stack.imgur.com/gTm0h.jpg)


title: How can I add all the 4 corners as GCP in an image that I want to reproject?

text: I am working on a Modis - oceancolor image and I want to reproject it to wgs84. I have the upper left/right and lower left/right values from the metadata of the image and I want to add them all in the command line of `gdal_translate` or `gdal_warp`. But in these two commands I can only add the upper left and lower right lat/long (x,y). Do you have any idea how I can add all the corners of the image in the script?


title: Does 'r.patch' act only on "loaded" rasters?

text: Can I run the GRASS command 'r.patch' on raster files which are not "loaded" into the current GRASS session? That is, rasters which are simply in my folder but not in a GRASS location & mapset? Thanks.


title: How to reproject a raster file in QGIS with datum transformation?

text: I want to reproject a tiff raster file from projection A (ETRS89/UTM Zone 32) to projection B (DHDN/Gauss-Kruger zone 2). The two projections have a different datum so I have to use a datum transformation (NTv2 BETa2007 in my case). I want to create a new raster file with projection B because I want to use it in another software. So reprojection on-the-fly is not what I need. I have tried Raster > Transform (or similar - I do not know the English menu) but I can not see how to define a datum transformation here. How can I do a permanent reprojection of a raster file in Qgis with a datum transformation?


title: What is a good way to retrieve geographic data along with imgs and mp3 files?

text: I'm working in the Django backend of an audio guide for tourism with mobile devices. Basically, the client (Phonegap app for Android and iOS) can follow a route over a map (Leaflet client) and see images or listen audio files in certain points of the route. Now, I want to provide offline service. So, the client can choose a route and _download it_ , meaning _download all the POIs, including images an audio files_. Of course, I can send the coordinates of the POIs in a json, or a gpx file. But what concerns me is the audio files and the images. Some random thoughts:   * The devices can store the coordinates in SQLite, and the audio files and images in the SD card   * Store the audio files and images as BLOBs in SQLite is not an option, because some files can be really big (like 100MB)   * My first approach is to create a JSON with the coordinate points and the name of the files related with them, and compress it along with the mp3 and imgs. Then, send that zipped file to the clients. Sounds... gross. Is there any recommended way to do this? Send geolocated audio and images, basically. Many thanks in advance


title: GRASS and QGIS- " not recognized...format" but it worked before

text: Prologue: I used my Ubuntu file viewing GUI to move data files around and "clean house" a bit. I started a new GRASS session and tried to add NED elevation rasters which I have successfully added before. I got the error message "`/home/greg/GRASS/Loc- Oakland/NED_elev_rasters/CopyNED1-9Oaklnd/10959312\10959312\w001001.adf' not recognised as a supported file format." Since the software recognized it before, why the change? Is it because I moved the files? QGIS gave me the same sort of error message, although I had successfully loaded these files before. Thank you.


title: What is the correct approach to handle mouse events for a tool that is called by a button?

text: I have a dockable window with a 2 buttons that call 2 different tools (Draw Marker and Select Features), What i want to happen is when the user clicks on the relevant button on the dockable window and then clicks on the map, a onMouseUp event is called and it runs code i need to run. Where does the onMouseUp event need to be placed so this can be done? Thanks


Well it might be a simple question. Given a set of 2D or 3D points, the question is: **How to find the center of geometry of an object?** According to the following figure, the center of geometry differs from the center of mass if it is calculated in the simplest form i.e., homogenous density of mass. The problem appears, indeed, in the computation of those. Commonly, one approach is to average **X** coordinates and **Y** coordinates separately i.e., find an average position to the given points (here in 2D). This may be used as centroid for the set of points representing an object. As shown, _due to the extra vertex along the bottom edge,_ for a simple rectangle the resulting centroid is **(0.5,0.4)** while the correct answer is **(0.5,0.5)**.   Note that the example given is too simple. The problem of interest however is for complex shapes in 2D and objects in 3D for which only coordinates of vertices are available.   BTW, an efficient computational way is of interest. **Update 1:**   Just to mention that I have checked some web links such as Wikipedia's however my current problem is that there are group of 2D and 3D points wish to find a point as representative for those. Thus centroid became of interest. The points are given without any topological information. You may consider them as point cloud. The demonstration here provided to make it clear that the commonly known averaging of coordinates (see for example this) may be incorrect as shown in the example. Hope the point is clear now. ![enter image description here](http://i.stack.imgur.com/hOxTV.png) **Update 2:**   Here are some implementations for comparison:   * **aa** = accepted answer below   * chull = convex-hull of points i.e., the golden polygon   * cent = centroid proposed in Wikipedia and discussed in **aa** as the polygon centroid   * centl = centroid of polyline as explained in **aa** Visually, `centl` looks better representative for the geometry given compared to `cent`. Two others look promising here but usually they are too biased if the dispersion of the points were inhomogeneous as it is a usual case.   And also consider that although convex-hull makes problem reasonably simpler however it may generate too long and too short edges without any symmetrical positioning in the space, that is, awareness is necessary if you do simple averaging (i.e., without weighting) for both cases: entire points (green) or convex-hull polygon vertices (blue). ![enter image description here](http://i.stack.imgur.com/QO3Vg.png) **Update 3:**   One application could be found in this post.


title: problem to Disable the security in server.dat file

text: I am trying to remove the security from web services.While seaching , I found that we have to set               <SecurityEnabled>True</SecurityEnabled>                            To          <SecurityEnabled>False</SecurityEnabled>      Now after replacing this `true` with `false` , I saved the file. Now when I restart my computer and reopen this file, It is again set to `<SecurityEnabled>True</SecurityEnabled>` .I am unable to understand why this is happening?


title: What is the highest obtainable resolution of FREE elevation data in the UK?

text: I am seeking closure with this question guys. I am aware of various sources of free elevation data ~30+ m, and I have been searching for higher resolution FREE data for a while now. This already exists somewhat (Where to acquire elevation data for UK?), but may have changed since a year. I'm afraid I may know the answer already, but to put an end to my misery and false hope: **What is the highest resolution, freely available (outside of academic institutions), elevation data covering the ENTIRE UK? - is 10 m available?** Thanks.


title: How to convert FROM Navit to PostGis or SHP?

text: I have a Navit BIN file which I want to use in PostGis or QGis, so I need to import it to PostGis or convert it to a shapefile. I have found a lot of information on how to convert data TO the Navit format, but nothing on how to do the reverse. Is this possible? Is there an exporter? Many thanks


title: What's the best way to display thousands of polygons in Openlayers?

text: I am trying to create an Openlayers map that will display thousands of polygons. And when the user clicks on each, I want some of the metadata to display. Working with this large of a dataset seems like it will be challenging. But I am trying to figure out if there is a standard way to go about this. At this point, I am think I should either   1. display my polygons with WMS and use getfeatureinfo   2. Display with WFS, Load with Strategy.BBOX, and use GetFeature to get meatadata. Is there a logical way I am missing to load all of this data into a map and have a good speed. I was looking at this page and trying to figure out what they did: http://protectedplanet.net/ Thanks for any help!


title: GeoTools WFS Feature retrieval

text: I am trying to use geotools to retrieve the list of features served by a WFS. I am following the tutorial provided at http://docs.geotools.org/latest/userguide/library/data/wfs.html. The code that I have is as follows:               // Step 1 - connection parameters     //     Map connectionParameters = new HashMap();     connectionParameters.put("WFSDataStoreFactory:GET_CAPABILITIES_URL", capabilitiesUrl.toString() );          // Step 2 - connection     DataStore data = DataStoreFinder.getDataStore( connectionParameters );          // Step 3 - discouvery     String typeNames[] = data.getTypeNames();     String typeName = typeNames[0];     SimpleFeatureType schema = data.getSchema( typeName );     System.out.println( "Schema Attributes:"+schema.getAttributeCount() );          // Step 4 - target     SimpleFeatureSource source = data.getFeatureSource( typeName );     System.out.println( "Metadata Bounds:"+ source.getBounds() );          SimpleFeatureCollection features = source.getFeatures();     SimpleFeatureIterator iterator = features.features();     while( iterator.hasNext() ){         Feature feature = (Feature) iterator.next();         System.out.println(feature.getName());     }      However, my DataStore object is always null. I have tried this using geotools version 9.1, version 10-SNAPSHOT and version 11-SNAPSHOT. I am using Maven to manage the dependancies. When I use version 11-SNAPSHOT, I keep getting a DataStoreFinder NoClassDefFound on running the above code. When using version 11-snapshot, I cannot see any of the geotools jars in the maven local repository which is probably why this error occurs. However, I have not yet managed to resolve it. I know that geotools can be used to retrieve features as shown above and would appreciate any help in pointing out what I could possibly be doing wrong. Thanks! UPDATE: I have included the pom files below. The first POM results in the following error on running mvn clean install:               Exception in thread "main" java.lang.NoClassDefFoundError: org/geotools/data/DataStoreFinder         at com.envitia.geopackage.geopackagecreator.dataharvestor.WFSDataHarvestor.getFeatureGeometries(WFSDataHarvestor.java:295)         at com.envitia.geopackage.geopackagecreator.dataharvestor.WFSDataHarvestor.getFeatureType(WFSDataHarvestor.java:274)         at com.envitia.geopackage.geopackagecreator.dataharvestor.WFSDataHarvestor.readData(WFSDataHarvestor.java:176)         at com.envitia.geopackage.geopackagecreator.App.main(App.java:22)      POM 1:               <project xmlns="http://maven.apache.org/POM/4.0.0" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"         xsi:schemaLocation="http://maven.apache.org/POM/4.0.0 http://maven.apache.org/xsd/maven-4.0.0.xsd">         <modelVersion>4.0.0</modelVersion>              <groupId>com.sample.harvestor</groupId>         <artifactId>Creator</artifactId>         <version>1.0-SNAPSHOT</version>         <packaging>jar</packaging>              <name>Creator</name>         <url>http://maven.apache.org</url>              <properties>             <project.build.sourceEncoding>UTF-8</project.build.sourceEncoding>             <geotools.version>11-SNAPSHOT</geotools.version>         </properties>              <dependencies>             <dependency>                 <groupId>junit</groupId>                 <artifactId>junit</artifactId>                 <version>3.8.1</version>                 <scope>test</scope>             </dependency>             <dependency>                 <groupId>org.geotools</groupId>                 <artifactId>gt-main</artifactId>                 <version>${geotools.version}</version>                 <type>jar</type>             </dependency>             <dependency>                 <groupId>org.geotools</groupId>                 <artifactId>gt-opengis</artifactId>                 <version>${geotools.version}</version>                 <type>jar</type>             </dependency>             <dependency>                 <groupId>org.geotools</groupId>                 <artifactId>gt-api</artifactId>                 <version>${geotools.version}</version>                 <type>jar</type>             </dependency>         </dependencies>         <build>             <plugins>               <plugin>                 <artifactId>maven-assembly-plugin</artifactId>                 <configuration>                   <descriptorRefs>                     <descriptorRef>jar-with-dependencies</descriptorRef>                   </descriptorRefs>                   <archive>                     <manifest>                       <mainClass>com.sample.harvestor.creator.App</mainClass>                     </manifest>                   </archive>                 </configuration>               </plugin>             </plugins>         </build>              <repositories>             <repository>                 <id>maven2-repository.dev.java.net</id>                 <name>Java.net repository</name>                 <url>http://download.java.net/maven/2</url>             </repository>             <repository>                 <id>osgeo</id>                 <name>Open Source Geospatial Foundation Repository</name>                 <url>http://download.osgeo.org/webdav/geotools/</url>             </repository>             <repository> <!--Add the snapshot repository here-->                 <snapshots>                     <enabled>true</enabled>                 </snapshots>                 <id>opengeo</id>                 <name>OpenGeo Maven Repository</name>                 <url>http://repo.opengeo.org</url>             </repository>         </repositories>     </project>      The next POM does not return any errors on mvn install, but on running the application, I get the following error:               Exception in thread "main" java.lang.NullPointerException         at com.envitia.geopackage.geopackagecreator.dataharvestor.WFSDataHarvestor.getFeatureGeometries(WFSDataHarvestor.java:298)         at com.envitia.geopackage.geopackagecreator.dataharvestor.WFSDataHarvestor.getFeatureType(WFSDataHarvestor.java:274)         at com.envitia.geopackage.geopackagecreator.dataharvestor.WFSDataHarvestor.readData(WFSDataHarvestor.java:176)         at com.envitia.geopackage.geopackagecreator.App.main(App.java:22)      This NullPointer is due to the data object being null as mentioned earlier. POM2:               <project xmlns="http://maven.apache.org/POM/4.0.0" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"         xsi:schemaLocation="http://maven.apache.org/POM/4.0.0 http://maven.apache.org/xsd/maven-4.0.0.xsd">         <modelVersion>4.0.0</modelVersion>              <groupId>com.sample.harvestor</groupId>         <artifactId>Harvestor</artifactId>         <version>1.0-SNAPSHOT</version>         <packaging>jar</packaging>              <name>Harvestor</name>         <url>http://maven.apache.org</url>              <properties>             <project.build.sourceEncoding>UTF-8</project.build.sourceEncoding>             <geotools.version>10-SNAPSHOT</geotools.version>         </properties>              <dependencies>             <dependency>                 <groupId>junit</groupId>                 <artifactId>junit</artifactId>                 <version>3.8.1</version>                 <scope>test</scope>             </dependency>             <dependency>                 <groupId>org.geotools</groupId>                 <artifactId>gt-main</artifactId>                 <version>${geotools.version}</version>                 <type>jar</type>             </dependency>             <dependency>                 <groupId>org.geotools</groupId>                 <artifactId>gt-opengis</artifactId>                 <version>${geotools.version}</version>                 <type>jar</type>             </dependency>             <dependency>                 <groupId>org.geotools</groupId>                 <artifactId>gt-api</artifactId>                 <version>${geotools.version}</version>                 <type>jar</type>             </dependency>         </dependencies>         <build>             <plugins>               <plugin>                 <artifactId>maven-assembly-plugin</artifactId>                 <configuration>                   <descriptorRefs>                     <descriptorRef>jar-with-dependencies</descriptorRef>                   </descriptorRefs>                   <archive>                     <manifest>                       <mainClass>com.sample.temp.harvestor.App</mainClass>                     </manifest>                   </archive>                 </configuration>               </plugin>             </plugins>         </build>              <repositories>             <repository>                 <id>maven2-repository.dev.java.net</id>                 <name>Java.net repository</name>                 <url>http://download.java.net/maven/2</url>             </repository>             <repository>                 <id>osgeo</id>                 <name>Open Source Geospatial Foundation Repository</name>                 <url>http://download.osgeo.org/webdav/geotools/</url>             </repository>             <repository> <!--Add the snapshot repository here-->                 <snapshots>                     <enabled>true</enabled>                 </snapshots>                 <id>opengeo</id>                 <name>OpenGeo Maven Repository</name>                 <url>http://repo.opengeo.org</url>             </repository>         </repositories>     </project>


title: Adding ArcMap standard toolbar to the form in C# .NET WinForms project

text: I have Windows Forms project with ArcObjects framework. I have already added axMapControl1 to the form. ![enter image description here](http://i.stack.imgur.com/zbdMe.png) The question is How can I add a standard toolbar to the form and link it with axMapControl1 ? Like this: ![enter image description here](http://i.stack.imgur.com/oyUaK.png) Is there any standard toolbar control in the ArcObjects or should make it by myself?


title: ArcGIS tool for indicating changes in land use

text: I have a couple of shapefiles in which there is some information about land use (polygons) such as water, forests, soils, buildings etc. Each of the shapefiles represents a certain year of the same area of interest. I want to indicate the area which was the biggest change over the years. Can you suggest me a procedure to do it in ArcGIS 10 ?


title: Retrieve WMS service from GeoWebCache with OpenLayers

text: Ok guys, here we go ! For my needs, I basically added a **standalone version** of **GeoWebCache**. After some configuration, the **GeoWebCache** retrieve the data from **GeoServer** using Capabilities. Now, in my **OpenLayers** code, I would like to question **GeoWebCache** to retrieve my data, and not **GeoServer** (if it is possible). For the moment, I have something like this (working) :                   var fond_carto = new OpenLayers.Layer.WMS(             "Fond cartographique 2",              "http://179.172.277.77:8080/geoserver/wms",             {                 layers: 'departement',                 TILED: 'true'             }      I tried stuff like that (not working) :                   var fond_carto = new OpenLayers.Layer.WMS(             "Fond cartographique 2",              "http://179.172.277.77:8080/geowebcache/sevice/wms",             {                 layers: 'departement',                 TILED: 'true'             }       Is there a way to accomplish that ? After some reading, I assume I'm a bit lost with this standalone version of **GeoWebCache**. It is written it works as a proxy between client ( **OpenLayers** ) and server ( **GeoServer** ). But I read that the TILED: 'true' option means that **GeoServer** will send the layers for tiling to **GeoWebCache**. As I desactivated the integrated version of **GeoWebCache** in **GeoServer** , and as my WMS request are still tiled, I assume this is true. But then it's not really a proxy. Could it works as follow: _openlayers => standalone geowebcache => geoserver => standalone geowebcache => openlayers_ Or it has to work that way : _openlayers => geoserver => standalone geowebcache => geoserver => openlayers_ **Thanks for any tips or any clarifications :-) !** # EDIT After some more configuration in geowebcache.xml, I still have the error but I can add some details to help. This is the URL of the request that is not displaying (pink tiles) :               http://179.172.277.77:8088/geowebcache/sevice/wms?LAYERS=departement&TILED=true&SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&STYLES=&FORMAT=image%2Fjpeg&SRS=EPSG%3A4326&BBOX=2.8125,48.515625,3.515625,49.21875&WIDTH=256&HEIGHT=256:8080/geowebcache/sevice/wms?LAYERS=departement&TILED=true&SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&STYLES=&FORMAT=image%2Fjpeg&SRS=EPSG%3A4326&BBOX=2.8125,48.515625,3.515625,49.21875&WIDTH=256&HEIGHT=256 :      and when I try to copy/paste this url, I got a GWC error :               400: Thread 35 Unknown layer departement. Check the logfiles, it may not have loaded properly.      I don't get why I got this error. Here's my geowebcache.xml description of the layer :               <wmsLayer>       <name>departement</name>       <mimeFormats>         <string>image/gif</string>         <string>image/jpeg</string>         <string>image/png</string>         <string>image/png8</string>       </mimeFormats>       <gridSubsets>         <gridSubset>           <gridSetName>EPSG:4326</gridSetName>         </gridSubset>       </gridSubsets>       <wmsUrl>         <string>http://179.172.277.77:8080/geoserver/wms</string>       </wmsUrl>       <wmsLayers>aaz:departement</wmsLayers>     </wmsLayer>


title: How to get Database name for layer in ArcObjects?

text: I am writing an ArcGIS ArcMap Add-In and need to determine what database a given layer from. How can I determine this information? I attempted to use the FeatureLayer to get a FeatureClass, and then get the FeatureDataset from that. From FeatureDataset I can get the name and parse out the first component that seems to be the DB name. However this approach seems really hacky and does not work if the layer is not part of a FeatureDataset. When I look at the layer's properties, it clearly shows the Database (see the highlighted line in the screenshot) so I think I can get the DB from a layer but I don't know what interface to use. ![enter image description here](http://i.stack.imgur.com/bpHL1.gif)


title: ArcEngine: Drawing graphics returned from BackgroundWorker

text: My application is drawing graphics objects on the map(as simple as adding IElement to GraphicsContainer). When I create graphics elements and add them to GraphicsContainer in BackgroundWorker DoWork method, ActiveView Refresh method does not show any added graphics. Are there any triks to take into account when drawing on the map in different thread? (if at all possible) Any help will be greatly appreciated.


title: Measure area of a polygon in specific units using ArcPy

text: I am trying to find the area of a particular polygon (`OriginalPolygon`, below), and I can do that... But, I'd like it in acres regardless of the original units of the polygon.(Data is coming into this script in a range of projections, but we're always interested in "acreage.") If I was calculating geometry to store in a field, then I could use `feat.area@acres` \-- but I don't need this stored as an attribute, I just need it as a temporary variable to tell the user. This question implies that unless I'm actually calculating geometry, the units of the shape area measurement are tied to the projection and therefore fixed. So: Is there any way to calculate geometry in such a way that the value is only a numerical output that I can display as the script runs, or am I stuck with making an attribute field, calculating geometry, reading the field value to tell the user, and then deleting the field?               rows = arcpy.SearchCursor(OriginalPolygon)     shapeName = arcpy.Describe(OriginalPolygon).shapeFieldName     for row in rows:         feat = row.getValue(shapeName)         PolygonArea = feat.area


title: How to locate a point on google maps/earth using cadastral maps in Kenya

text: I'm an architect, so not too well versed with all this datum/transformation talk. My little problem: i want to be able to enter a position using UTM or deg/sec/min coordinates into google earth and fly there. the cadastral maps we have come with eastings and northings and a general reference that reads: coordinates referred to Equator and 37* E (37deg E), all in meters station:A +X (Easting):116379.70 -Y(Northing) : 6309.10 any help out there??


title: Unmodified OpenGeo Suite Geoexplorer proxy fails with HTTP 500, regardless of servlet

text: We're running an Opengeo Suite (3.0.1 currently) and an unknown amount of time ago the proxy redirects to other data stores GeoExplorer does* started failing with HTTP 500 (Internal server error). *               var app = new GeoExplorer.Composer({             authStatus: {{status}},             proxy: "../proxy/?url=",      applied to every store with an absolute URL(doesn't matter whether valid or not), like the default               suite: {url: "http://v2.suite.opengeo.org/geoserver/ows      so that the first call becomes `GET /geoexplorer/proxy/?url=http://v2.suite.opengeo.org/geoserver/ows?_dc=4734342332&SERVICE=WMS`. In this environment a newly installed and unmodified Suite does this, and an unmodified GeoExplorer webapp does the same thing in Tomcat. The specific error, same from jetty as Tomcat, is " _Error: No valid JSGI response: [object Object] (ringo/jsgi/connector.js#88)_ ". The environment is a Windows Server 2008. In other environments (on other machines) the Suite has no such issues, so I'm wondering whether it has to do with some network setting or other in Windows? Further, I'm not completely sure what this proxy is for and have been unable to find any documentation for it. Any hints would be very welcome :) Thanks for any hints. Edit. Further research has revealed that this error message (HTTP500, ringo/jsgi/connector.js:88) appears if the machine hosting the servlet has no internet connection. It has been observed on a different machine, running Win 7. I wonder if the servlets in the affected Win 2008 server environment are somehow prevented from reaching the internet, and if this can be detected in any logs / audit files of Win 2008 server R2?


title: How do I convert MULTIPOLYS into LINESTRINGS in PostGIS

text: I have a table with a geom column as a geometry(MULTIPOLYGON) I would like to create an additional column of type geometry(LINESTRING) How can I do this using a PostGIS Query? I need to be able to use the linestring in Geoserver to style labels along the line.


title: Area calculation of an object without coordinates

text: I have a photo of an object with a scale and I want to calculate its area. I guess I have to create a polygon of the area to be calculated to use the 'field calculator'. But how do I tell the program what a centimeter (cm) is?


title: Getting compilation error on esriLicenseProductCode.esriLicenseProductCodeArcView

text: I am getting compilation error on esriLicenseProductCode.esriLicenseProductCodeArcView   1. I have installed Arcgis 10.1 and authorized ArcEditor. I am not able to authorized ArcView.   2. I have not installed Runtime Engine. Am I missing anything any suggestions.


title: Represent full globe in Oracle Spatial

text: In SQL Server 2012 I can do this:               -- Construct a new FullGlobe object (a WGS84 ellipsoid)     DECLARE @Earth geography = 'FULLGLOBE'          -- Calculate the area of the earth     SELECT PlanetArea = @Earth.STArea()      Is there a similar construct in Oracle Spatial? Alternatively what would be the best way to represent the entire earth in an 8307/4326 SRID? Is there anything in the WKT standard that speaks to this?


title: How to install Quantum GIS Lisboa and scipy

text: I'm trying to install some plugins of the QGIS but an error occurs. The program asks the scipy program. I installed Python, the Numpy and scipy. Also installed the program R. However, QGIS continues asking for scipy. Could you help me?


title: Convert Gauss-Kruger to UTM

text: How to covert Gauss-Kruger coordinates to UTM coordinates with GeoTools java library? As I have figured out the location (zone, country) is important. I need the calculation for Slovenia.


title: Join is deleting fields in arcpy

text: I'm working on a script to automate the conversion of CAD .dwgs of rooms to feature classes, I'm nearing completion and am now working on the second portion of the script (I've seperated it into two modules as the first part takes about an hour to run and the second part should execute MUCH faster). Part of the script involves adding two fields to one layer, and then joining it to another layer. The problem is, when I execute "arcpy.AddJoin_management" the fields I added aren't in the attribute table. Whats puzzeling is the fields show up when I crate a shapefile from the layer when I perform "arcpy.CopyFeatures_management" before I execute the join.               import arcpy     from arcpy import env     import os     import math          # Set the MTolerance to 0.02     arcpy.env.MTolerance = 0.02                         # Set overwrite option          arcpy.env.overwriteOutput = True          #Set Workspace          arcpy.env.workspace = "in_memory"          #input dataset from part 2     converted_CAD = r"C:\Users\JOE.AWESOME\test_script_output\test.shp"          #create a minimum boundary geometry     arcpy.MinimumBoundingGeometry_management (converted_CAD, 'bounding_geom', 'RECTANGLE_BY_AREA','', '', 'MBG_FIELDS')                         #create layers for an attribute join     arcpy.MakeFeatureLayer_management (converted_CAD, 'CAD_layer')     arcpy.MakeFeatureLayer_management ('bounding_geom', 'bounding_geom_layer')          #Add a field to the table for the difference between the Minimum Bounding Geometry and the CAD_layer area          arcpy.AddField_management ('CAD_layer', 'area_diffe', 'FLOAT')          arcpy.AddField_management ('CAD_layer', 'area_sqM', 'FLOAT')                    #join the two layers for analysis          arcpy.AddJoin_management ('CAD_layer', 'FID', 'bounding_geom_layer', 'ORIG_FID')          #calculate the area of the polygons in the CAD_layer          arcpy.CalculateField_management ('CAD_layer', 'area_sqM', '!shape.area@squaremeters!', 'PYTHON_9.3')          arcpy.CopyFeatures_management ('CAD_layer', r"C:\Users\JOE.AWESOME\test_script_output\join_test.shp")      I've tried changing them to a FLOAT from a DOUBLE, but that didn't change anything. I'm stuck, maybe it's something little. In the words of Michael Bolton, "I always miss some mundane detail."


title: Exporting selected features to a new shape file

text: I have a street layer and I only need to export a subset of features to a new shape file where the "TestDone" attribute field is not '0' (Zero). My current code ( which is pasted below ) exports the total street layer. I can use a query filter to get the non zero values out and assign it to an ICursor. But then not sure how to conenter code herevert it to a IfeatureClass. ( Not really sure whether the ICursor and then trying to convert to IFeatureClass is the correct way to go even ) How can I just get features that have a non zero value in the "TestDone" attribute field and then convert that to a feature layer as the second argument in the ExportFeatureClass method ( please see below code) needs a IFeatureClass. Any help with some code guidance will be greatly appreciated. _**_*CODE _*_ ****               private void btnExport_Click(object sender, EventArgs e)     {     SaveFileDialog saveFileDialog1 = new SaveFileDialog();     saveFileDialog1.Filter = "Shape file|*.shp";     saveFileDialog1.DefaultExt = "shp";     saveFileDialog1.AddExtension = true;     saveFileDialog1.InitialDirectory = @"C:/";     saveFileDialog1.Title = "Check data";     saveFileDialog1.ShowDialog();     if (saveFileDialog1.FileName != "")     {     string strFile = System .IO.Path.GetFileName ( saveFileDialog1.FileName) ;     IMap pMap = default(IMap);     IWorkspaceFactory pWorkspaceFactory = default(IWorkspaceFactory);     IWorkspace pWorkspace = default(IWorkspace);     IFileName pFileName = default(IFileName);     IFeatureLayer2 pFeatureLayer = default(IFeatureLayer2);     IFeatureClass pFeatureClass = default(IFeatureClass);     IPropertySet pPropertySet = default(IPropertySet);     IGeometryDef pGeometryDef = default(IGeometryDef);     IGeometryDefEdit pGeometryDefEdit = default(IGeometryDefEdit);     string deletefile = null;          pMap = clsFunctions.pMxDocument.FocusMap;     pFileName = new FileNameClass();     pFileName.Path = System .IO.Path.GetDirectoryName (saveFileDialog1.FileName);     pPropertySet = new PropertySet();     pPropertySet.SetProperty("DATABASE", pFileName.Path);     pWorkspaceFactory = new ESRI.ArcGIS.DataSourcesFile.ShapefileWorkspaceFact ory();     pWorkspace = pWorkspaceFactory.Open(pPropertySet, 0);          IGeographicCoordinateSystem pGCS = default(IGeographicCoordinateSystem);     SpatialReferenceEnvironment pSpatialReferenceEnv = default(SpatialReferenceEnvironment);     ISpatialReference2 pSpatialReference = default(ISpatialReference2);     pSpatialReferenceEnv = new SpatialReferenceEnvironment();     pGCS = pSpatialReferenceEnv.CreateGeographicCoordinateSys tem(Convert.ToInt16(esriSRGeoCSType.esriSRGeoCS_WG S1984));     pSpatialReference = (ISpatialReference2)pGCS;          ILayer pLayer;     int intI = 0;              for (intI = 0; intI <= pMap.LayerCount - 1; intI++)         {         pLayer = pMap.get_Layer(intI);         if (pLayer.Name == cmbStreets.Text)         {         pFeatureLayer = (IFeatureLayer2)pLayer;         }         }          pFeatureClass = pFeatureLayer.FeatureClass;     pGeometryDef = new GeometryDef();     pGeometryDefEdit = (IGeometryDefEdit)pGeometryDef;     pGeometryDefEdit.GeometryType_2 = pFeatureClass.ShapeType;     pGeometryDefEdit.SpatialReference_2 = pSpatialReference;          ExportFeatureClass(strFile, pFeatureClass, pWorkspace, pSpatialReference, pGeometryDef);     }                    public void ExportFeatureClass(string strLayerName, IFeatureClass pFeatureClass, IWorkspace pOutWorkspace, ISpatialReference pSpatialReference, IGeometryDef pGeometryDef)     {     IFeatureClassName pInFeatureClassName = default(IFeatureClassName);     IDataset pDataset = default(IDataset);     IWorkspaceName pOutWorkspaceName = default(IWorkspaceName);     IFeatureClassName pOutFeatureClassName = default(IFeatureClassName);     IDatasetName pDatasetName = default(IDatasetName);     IExportOperation pExportOperation = default(IExportOperation);     IWorkspace pInWorkspace = default(IWorkspace);          pDataset = (IDataset)pFeatureClass;     pInFeatureClassName = (IFeatureClassName)pDataset.FullName;     pInWorkspace = pDataset.Workspace;     pDataset = (IDataset)pOutWorkspace;     pOutWorkspaceName = (IWorkspaceName)pDataset.FullName;          pOutFeatureClassName = new FeatureClassNameClass();     pDatasetName = (IDatasetName)pOutFeatureClassName;     pDatasetName.Name = strLayerName;     pDatasetName.WorkspaceName = pOutWorkspaceName;          pExportOperation = new ExportOperation();     pExportOperation.ExportFeatureClass((IDatasetName) pInFeatureClassName, null, null, pGeometryDef, pOutFeatureClassName, 0);          //Cleanup     pInFeatureClassName = null;     pDataset = null;     pOutWorkspaceName = null;     pOutFeatureClassName = null;     pDatasetName = null;     pExportOperation = null;     pInWorkspace = null;     }


title: Why does arcpy append from file geodatabase to SDE fail?

text: I'm trying to append a feature class from a file geodatabase to an Arcsde featureclass using this line of code: `arcpy.Append_management("parcel_copy", sde_path+"parcel" ,"NO_TEST")` This line works under Python console in ArcMap , but when I execute it under Pyscripter IDE, all the code works fine except the append statement, and it does not return an error. * * * Thanks for you responses, i have tried to isolate the `append` statements in a little script, and it worked, but i steel don't understand why it has not worked. i wonder if i have not deleted a featurelayer refering to a featureclass in SDE, that locks this featureclass, so i cannot append to it. any explanation is welcome


title: Looking for gis server for pictometry that can do bing-like "birds eye" tilted view and rotation

text: I need an open source gis server that can present these pictometry images, I guess they are also called 2.5D. The camera on a plane takes 5 shots, one orthogonal, other 4 in all directions. North, west... 2.5D images are the main focus here, user would have to be able to rotate the tilted view in 4 directions. Additionally I would like to have those standard gis server features like vector layers, queries, labels, redlining... So far I have decent experience with mapguide, open source and the the autodesk one, and I can't recall those features. No rotation or tilting at all, I might be wrong. I have very basic xp with geoserver and mapserver. I'm looking for advice in which direction to go with research. I was thinking geo-server is the best bet, but what interface... Also, I cant really say that I've been working with any kind of 3d desktop gis. I would take some advice on that too. Any suggestions? Thx.


title: GRASS r.out.mat memory issues

text: we have built a model with pythn/GRASS that calls the r.out.mat module. The area is quite large 25,000x25,000 cells. The machine we are running it on is Windows 7 machine with 16GB of RAM. After a while the script aborts complaining about an out of memory issues. Now, I would not have expected to have such issues with so much RAM even though the manual page for r.out.mat does say that the module loads the whole map into RAM so it requires a huge amount of RAM. Any idea on how we can optimise the script avoid such problems? If this is a limitation of the module, I was wondering whether r.out.mat could be used to read the map area in patches, dump out the results in single files and then some other module patches all the r.out.mat results together. Any suggestions on how we can overcome this issue? Thanks Antonio


title: Arcpad 10...using sql with VBScript

text: I am trying to get a count of related records (Tableform) in order to validate Editform entry (using VBS). Any help getting myCount? I Know there are some syntax issues here, but is this even possible in Arcpad axf? Seems simple enough, no brainer using VBA but something like this:               Sub Validate()         myCount = (Select count(*) from tblRelatedRecs where ID = Pcontrol.value Group by ID            If myCount > 5 then             Application.MessageBox ("Huma, you erred!")             ThisEvent.Result = False         End If     End Sub


title: Zonal statistics returns wrong values in QGIS

text: I used zonal statistics for a simple task: to transmit raster pixel values to vectormap which contains squares as big as the pixel size is. The problem is that the values of the outcome are not correct, for example: pixel value on the rastermap is 11, but when I look the same square on the vectormap, it is 76. What can cause the problem? Thank you for answering!


title: Problems with other applications after installing QGIS 2.0.1

text: Since the official QGIS installation 2.0.1 (Virtual Server) various core applications from other departments in the municipality do not work always properly. Does anyone know the problem? Thank you for your help Best regards


title: Viewer Function

text: Is there a Viewer function/display in QGis as it is in ArcGIS. Or is there an other possibility to limit the access (just visualisation, no analysis). Thank you for your answer.


title: Snippets applicable for ArcObjects 10.2.1

text: I am trying to create ArcGIS desktop add-in using Visual Studio Professional 2012. But when I use ArcGIS Snippets Finder to insert code, I found the applicable versions as follows: ' Applicable ArcGIS Product Versions: ' 9.2 ' 9.3 ' 9.3.1 ' 10.0 But I am using ArcGIS 10.2.1, and the ArcObjects obtained from ArcGIS 10.2.1 Desktop installation CD. I want to know where can get the snippets applicable for ArcGIS 10.2.1 (or above). And what version of Visual Studio can support ArcObjects for 10.2.1?


title: Socket Server to receive GPS location data

text: I am working on a web application which receives data sent from trackers via TCP Packets containing location information. Up until this point we have been working with dummy data within the mysql database but we are now ready to start building the server which will receive data from the remote devices. As we are just exploring our options at the moment, I thought I would see does anyone have any suggestions on the best route to take when setting up a server to receive TCP packets of data and then store them in a remote database? Our general requirements are as follows:   * Hosted on Amazon EC2   * Allow multiple concurrent connections to receive data (potentially in the thousands)   * Store data in remote MySQL database Any guidance would be greatly appreciated!


title: ImageLab Program Problem

text: I am new to GeoTools. I was trying the samples in http://docs.geotools.org/latest/userguide/tutorial/raster/image.html.   I am not able to run the ImageLab sample code. I know it would be some lack in understanding or a simple mistake. I am getting the following exception: Exception in thread "main" java.lang.UnsupportedOperationException: Trying to get a reader from an unknown format. at org.geotools.coverage.grid.io.UnknownFormat.getReader(UnknownFormat.java:62) at org.geotools.demo.ImageLab.displayLayers(ImageLab.java:93) at org.geotools.demo.ImageLab.getLayersAndDisplay(ImageLab.java:77) at org.geotools.demo.ImageLab.main(ImageLab.java:51) What actually needs to be done? Also, is it possible to use normal jpeg image or any other way around? In the GeoTools website the mentioned bluemarble.tif is not available and I have been using normal jpg image. In another website while seraching for more info on this, I read that it uses a world file. What is a world file and how it is related to this and how do I use it? Please help, I have been stuck with this for a long time. Thanks in advance.


title: Naming layers in DXF using OGR2OGR

text: We have been using ogr2ogr to export data from PostGIS into DXF. Using the default behaviour of the export, all data are put into layer '0'. We would put data into a layer name from the PG table name. It is not really clear if the correct solution is manipulating the header.dxf under the GDAL library or if this can be achieved by using some of the command line options from ogr2ogr. Thanks for any input. Rgs. Ulrik


title: Postgres Log Error

text: I am using Geoserver + GeoWebCache + Postgres for an application. On the application side there is no bugs but on postgres I am getting this log. Please explain and help me to solve this > 2014-03-14 10:14:50 IST ERROR: relation "public.gt_pk_metadata" does not > exist at character 15 2014-03-14 10:14:50 IST STATEMENT: SELECT * FROM > public.GT_PK_METADATA 2014-03-14 10:14:50 IST ERROR: relation > "public.gt_pk_metadata" does not exist at character 15 2014-03-14 10:14:50 > IST STATEMENT: SELECT * FROM public.GT_PK_METADATA 2014-03-14 10:14:50 IST > ERROR: relation "public.gt_pk_metadata" does not exist at character 15 > 2014-03-14 10:14:50 IST STATEMENT: SELECT * FROM public.GT_PK_METADATA > 2014-03-14 10:14:50 IST ERROR: relation "public.gt_pk_metadata" does not > exist at character 15 2014-03-14 10:14:50 IST STATEMENT: SELECT * FROM > public.GT_PK_METADATA 2014-03-14 10:14:50 IST ERROR: relation > "public.gt_pk_metadata" does not exist at character 15 2014-03-14 10:14:50 > IST STATEMENT: SELECT * FROM public.GT_PK_METADATA 2014-03-14 10:14:50 IST > ERROR: relation "public.gt_pk_metadata" does not exist at character 15 > 2014-03-14 10:14:50 IST STATEMENT: SELECT * FROM public.GT_PK_METADATA > 2014-03-14 10:14:50 IST ERROR: relation "public.gt_pk_metadata" does not > exist at character 15 2014-03-14 10:14:50 IST STATEMENT: SELECT * FROM > public.GT_PK_METADATA 2014-03-14 10:14:50 IST ERROR: relation > "public.gt_pk_metadata" does not exist at character 15 2014-03-14 10:14:50 > IST STATEMENT: SELECT * FROM public.GT_PK_METADATA 2014-03-14 10:14:50 IST > ERROR: relation "public.gt_pk_metadata" does not exist at character 15 > 2014-03-14 10:14:50 IST STATEMENT: SELECT * FROM public.GT_PK_METADATA > 2014-03-14 10:14:50 IST ERROR: relation "public.gt_pk_metadata" does not > exist at character 15 2014-03-14 10:14:50 IST STATEMENT: SELECT * FROM > public.GT_PK_METADATA 2014-03-14 10:14:50 IST ERROR: relation > "public.gt_pk_metadata" does not exist at character 15 2014-03-14 10:14:50 > IST STATEMENT: SELECT * FROM public.GT_PK_METADATA 2014-03-14 10:14:50 IST > ERROR: relation "public.gt_pk_metadata" does not exist at character 15 > 2014-03-14 10:14:50 IST STATEMENT: SELECT * FROM public.GT_PK_METADATA > 2014-03-14 10:14:50 IST ERROR: relation "public.gt_pk_metadata" does not > exist at character 15 2014-03-14 10:14:50 IST STATEMENT: SELECT * FROM > public.GT_PK_METADATA 2014-03-14 10:14:50 IST ERROR: relation > "public.gt_pk_metadata" does not exist at character 15 2014-03-14 10:14:50 > IST STATEMENT: SELECT * FROM public.GT_PK_METADATA 2014-03-14 10:14:50 IST > ERROR: relation "public.gt_pk_metadata" does not exist at character 15 > 2014-03-14 10:14:50 IST STATEMENT: SELECT * FROM public.GT_PK_METADATA > 2014-03-14 10:14:50 IST ERROR: relation "public.gt_pk_metadata" does not > exist at character 15 2014-03-14 10:14:50 IST STATEMENT: SELECT * FROM > public.GT_PK_METADATA


title: Does the ArcGIS JavaScript API let me build / package my application using Dojo?

text: I'm currently getting to grips with the new (3.1) ArcGIS JavaScript API and how it bundles Dojo and relies on Dojo's Asynchronous Module Definition (AMD) approach. I'm familiar with AMD and recently completed a project using Require, Backbone, and OpenLayers in combination - overall a very pleasant experience. However I'm worried by   * how tightly-coupled Dojo is with the ArcGIS JavaScript API,   * that trying to use Require instead of Dojo AMD throws an exception, and   * that Dojo AMD doesn't seem to have the same user-base and momentum as Require, making it difficult to research / resolve problems without hitting the forums My main concern is that I will write an ArcGIS JavaScript API application and not be able to build / package it into a single / small number of scripts. This will leave me with an application making 100+ HTTP requests at startup for all of the modules it depends on. I've seen this article covering 3.0 but it doesn't mention the build process and I believe 3.0 is quite different. I've also seen this piece on building / packaging an application that uses the ArcGIS JavaScript API but the comments suggest this is pre-Dojo 1.7 AMD. Have ESRI really only half implemented the AMD approach, and in a way that will actually make my application slower? Why can I not find ESRI documentation on building / packaging a JavaScript API-based application? Any thoughts or experience here would be much appreciated. **EDIT** Based on the answers so far it seems that no, ESRI hasn't provided any information on building your Dojo application after inflating the number of scripts involved. It seems the only useful information available is from non-official sources and slightly outdated (enough to cause problems). This leads to me wonder why more people aren't vocal about the API's new improvements / deficiencies   * am I being too critical / demanding in expecting some more information here?   * is no one using the latest version of the API?   * should I use an older version of the API that doesn't depend on AMD?


title: Is it possible to set the output of a spatial etl tool in model builder?

text: I am using the interoperability extension and building spatial ETL tools, but I can't seem to specify the output location of the ETL tool so it's impossible to use this dynamically. Has anyone figured out how to do this?


title: Keyboard shortcut for GPS data capture in QGIS

text: I would like to setup keyboard shortcuts for GPS data capture with QGIS. Using GPS Information dialog box, I could use Mouse-Left-Click to "Add feature" and "Add track point", but I would like to assign keyboard shortcuts for them so that I can use a programmable Mouse. Settings --> Configure Shortcuts lists only few commands Is there a way I can call the API behind that and assign to a shortcut key? I appreciate your time. Thanks, Thiru


title: What multiplier would you use for calculating estimated population?

text: Given a number of address points in a given area, what multiplier would you use for calculating estimated population? I have seen 2.5 in a few places, but I was wondering what the generally accepted multiplier would be.


title: ArcGIS/arcpy: Can a variable be set to reference an in_memory workspace object?

text: This is an offshoot question related to a question I posted earlier. I have the following arcpy code that inputs a Python list of polygon geometry objects to the Union_analysis tool, which then writes the result to the in_memory workspace. The result is then input to the Dissolve_management tool, which in turn writes the result to memory. This all occurs near the end of a function, and the final goal here is to return the single geometry from the Dissolve to the calling function, where it can be passed along. Hence my question: **How can (can I?) I set the Dissolve output to a variable so that I can return it as a geometry object?** The first thing that will be done with it is to read its geometry into a point array to use as the inner ring for another polygon.                   if len(lstHSGeoms) > 1:             #Union geometries into one geometry             arcpy.Union_analysis(lstHSGeoms, "in_memory\\lstHSUnionGeom", "ONLY_FID")             arcpy.Dissolve_management("in_memory\\lstHSUnionGeom","in_memory\\lstHSOneGeom", None, None, "SINGLE_PART")             hsOneGeom = lstHSOneGeom      The last line is the problem here. No surprise, arcpy complains that the global name `lstHSOneGeom` is not defined. Putting in the whole in_memory path as the assignment value, i.e. `"in_memory\\lstHSOneGeom"` just assigns the simple string to the `hsOneGeom` variable. There is no reference to the in_memory object. I originally coded this block as the ArcGIS help files would have you do, which is to create empty arcpy Geometry objects and to set them as the outputs of the tools. When then used as inputs, they were rejected as invalid. At least I could keep a reference to the geometry objects in a variable, but the tools just wouldn't recognize them as valid input. So, I went with in_memory workspace update that you see in the snippet above. (Thanks to @dmahr for this solution in the previous question referenced above.) I now don't get errors from the tool, but I do get an error from the last line. Any suggestions? BTW, I should mention the vitals: ArcGIS/ArcView 10.0, Windows 7, 8 GB RAM.


title: Disable Certain MapTips in Graphics Layer for ArcGIS Silverlight

text: I am using the ESRI Silverlight API to allow users to use the Identify functionality on a map. The graphics layer on this map is set to use Map Tips, and a pin will be applied to the graphics layer when a user clicks on the map. However, in specific cases, I want the map tips to be disabled. I am using the MouseEnter event on the graphics layer to catch each graphic as the user rolls over them. Once inside, I can determine whether or not I want the Map Tip to show up based on the graphics' attributes. Is there a way for me to disable the Map Tips inside of the MouseEnter event? I don't want to disable them for the entire graphics layer. Thanks, Seth


title: Getting points within a range from a postgres database?

text: ok so this is query i'm making which returns no result:               SELECT          building.gid     FROM          building     WHERE         ST_Dwithin(building.geom, ST_GeomFromText('POINT(long lat)',4326), 400)      If I run this in QGIS, I get a syntax error, but if I run this on pgAdmin, i get no error, but simply 0 rows.... my long and lat are in degrees. what could possibly go wrong?


title: How to convert ESRI SHP to KML with correct coordinate system?

text: I am working on this SHP file, containing geographical district of Hong Kong. I would like to export this to KML and export it to Google Map for overlay. After I converted it to KML in QGIS (just save as kml format) and opened it in Google Map, it seems the coordinate is out of the entire map world region. Is there any conversion I need to do from ESRI SHP coordinate system to Google Map compatible coordinate system? Thank you for any pointer!


title: How to use Google Maps API V3 to diplay circle with Openlayers

text: I found bellow code to display Circle by click event with Radius 60.0m.(google maps API v3). http://stackoverflow.com/questions/15111403/google-maps- api-v3-change-circle-radius-by-click-event I want to display a circle (within a limit of 3000m, 5000m ) and the popup window to appear attribute details of buffer zone. (school name and the school type & etc.,). Please help me to properly place a this JavaScript code into my HTML page.


title: How to realize two overlaying layers? One with Popups, the other one with Info-div!

text: I'm trying to display two overlaying layers (geojson). Layer1 should give me information by hovering over it (works fine based on the Leaflet Choropleth Map Example. Layer2 is situated within Layer1 and should give me information by clicking which does not work. It seems as if Layer1 is in front and disables the interaction with the other one. I already tried the `bringToFront()` and `bringToBack()` methods to handle the problem. To checkout what I mean look here: http://jsfiddle.net/uFPwB/2/


title: Are there any resources or standards when writing a GIS data specification?

text: I have been searching for sources of GIS data specifications to assist with documenting a number of datasets that I'm currently working on. I'm used to looking at Ordnance Survey specifications, such as the one located here: http://www.ordnancesurvey.co.uk/docs/technical-specifications/os-mastermap- water-layer-technical-specification.pdf. I'm guessing that a good specification will be one which is independent of applications and formats, and describes the data in the detail which is necessary for its purpose and intended audience. Are there any resources to assist with the creation of GIS data specifications? Is there an industry standard way to do them so that they can be easily understood by both users and GIS developers? I'd also be interested to know what kind of documentation people do generally create for datasets which are to be widely used within their organisations. Many thanks.


title: Where is the network message log tab in QGIS on Windows7?

text: I want to find out what web requests is QGIS exactly sending to my web server when it connects to a WMS and WFS layer. In my Linux installation of QGIS 2.0 I can see three tabs in the Log Messages window: General, Plugins, Network. In my Windows installation of QGIS 2.0 I can only see two tabs in the Log Messages window: General, Plugins. How do I activate the "Network" log messages tab in QGIS 2.0 on Windows 7?


title: Is it possible to use .bundle and .bundlex files in QGIS?

text: I am working with a colleague that has given me UK OS data and Aerial image files to work with in .bundle and .bundlex form. This isn't a huge amount of data as we are primarily working to the 1:50,000 scale to show what is going on in a forest environment. I need to open this and so far am getting nowhere. Initial searches reveal nothing on QGIS. I did find something about opening it with ESRI GIS Desktop, but I am not on that system. If this can't be done with QGIS, I have an ArcGIS 3.2, on an old machine at home. Is that a possibility? Although I think I will run into speed and memory problems on that four year old laptop. All the files are on my Macbook running QGIS 1.8. (QGIS 2.0 not supported with OS 10.9.1) OK... Nothing like Saturday night for chasing down the meaning of arcane GIS terms. The compact cache groups many tiles together in one large file called a bundle. A bundle contains up to 16,384 tiles. The result is a cache with dozens or hundreds of files, instead of thousands or millions. If you look at a compact cache on disk, you can see the bundle files with the .bundle extension. You'll also see some corresponding index files with the extension .bundlx. It seems all ESRI GIS product read .bundle and .bundlex files? http://resources.arcgis.com/EN/HELP/MAIN/10.1/index.html#//0154000002Z3000000.HTM * * * FINAL UPDATE: Curious to note that having posted this here, and talked to a half dozen others I get no answers. The person who gave me the files is hitting a wall, doesn't know what it is. It works fine on her system, but we can not get anything out that will open in QGIS. We started to export AI files but as the files are scalable we need to printout dozens and dozens of sections on her system and she has no time to do that. (What happened to MrSid? it was much easier.) Upshot question here.... Is this -really- such an oddity that nobody knows what it is or how to come at it? The slightest guidance would be most helpful. I am now pressing up against a deadline. \--> Meanwhile I've unpacked a Six year old PC and spent two days cleaning it out and downloaeding the 60 day eval version of Esri, to have a go with that this afternoon. I don't have high hopes about this either to be honest. But we shall see. Nothing is easy at the moment, but if it were easy...


title: Looking for examples of bayesian modelling and smoothing in a GIS

text: I am looking to assemble clear examples of Bayesian modelling / smoothing in a GIS. I have yet to see GIS being used to conduct Bayesian style analytical models and would like to see what is currently out there. I am interested in both the application of Bayesian modelling and its current state of inclusion in traditional and/or non-proprietary (opensource) software. At the moment I am very interested in Bayesian modelling for transportation planning, specifically transit ridership and modal split/share estimations.


title: Striping Artifact in Aster DEM Mosaic

text: I'm working with SRTM and ASTER data for a university project. My goal is to delineate watersheds and evaluation of these with ArcGis 10 and ArcHydro. But a striping error occurs with the flow direction or hillshade as seen in the attached picture. I have set advanced arcmap setting to Cubic and also projected the data with cubic resampling. ![](http://i.imgur.com/3UVnrEC.jpg) This Strange gdaldem hillshade artifacts is a smiliar question, but for QGis :( Can I get rid of this? And if not, does this alter the following steps after flow direction?


title: Is there a way to keep labels and symbols from displaying inside an extent indicator?

text: I've been working on a utilities map and have been asked to turn off the labels and symbols within the extent indicator on the larger scale map and show them only in the detailed map. Is there a way to do this without converting labels to annotation and making a copy of and editing each individual shapefile specifically for this map?


title: MapInfo Professional 11.5 sample data

text: would you be able to point me in the direction of finding a copy of the sample data that comes on the CD's for learning MapInfo Professional 11.5?


title: Two boundary sets intersect - some are genuine diffferences, others slivers resulting from poor digitising. How to differentiate?

text: **Problem:** I have two boundary sets. Intersecting reveals that there are genuine differences between them - as expected. However, there are also a host of slivers due to their being digitised by different agencies. These mismatches may have statistical implications. I have tried several approaches to clearly differentiate one mismatch type from the other. **Current thinking:**   1. I get the MBR of each polygon in boundary set 1 and get the ratio of the area of the MBR to the area of the actual polygon. I do the same for boundary set 2.   2. I get the MBR of each intersection (these have been disaggregated to singlepart) and get the ratio of the area of each intersection's MBR to the area of the actual intersection.   3. The distribution of the ratio for slivers is likely to be materially different to the distribution of the ratio for genuine differences.  **QUESTIONS** Do you agree that this seems like a sensible approach? How easy is it to obtain the MBR (minXmaxXminYmaxY) for a polygon in ArcGIS 10.1 Advanced? Either point-and-click or Python suggestions welcome. If my thinking is flawed, please point out how, and suggest an alternative.


title: Atmospheric correction for Landsat 8 in ERDAS?

text: I am attempting to run some classification on Landsat 8 imagery, however, I need to perform atmospheric correction on the images in order to reveal true spectral signatures. I am using ERDAS, and the autonomous atmospheric correction tool does not support the Landsat 8 sensor. Unfortunately, I do not have access to ENVI which supports landsat 8, or a licence to ATCOR, which seems to be a ERDAS add on for atmospheric correction. Any Ideas? Thanks


title: How to repair Geometric Network in Versioned ArcSDE Geodatabase?

text: Our Versioned ArcSDE workflow has a feature dataset with Geometric Network in the parent (Default), which users derive in their own versions. It is now observed that in some versions, the geometric network is no longer valid and does not abide connectivity rules.   * How can we repair connectivity/enable Geometric Network again in these versions?    * Or should we do this on the default dataset, rather than in user versions?    * What are the tools which may help us to retrieve the Geometric Network ?


title: How to map and color localities (provinces) with an external CSV using Quantum GIS

text: I use a Swiss shape file map (called PLZO_PLZ) available on http://www.cadastre.ch/internet/cadastre/fr/home/products/plz/data.html (ZIP file Shape LV03) PLZO_PLZ contains a column called PLZ which correspond to a swiss locality. For instance number PLZ 1202 is in Geneva. By the way, PLZ can also be called NPA (same meaning). My aim is to get a colored map according to the NPA frequency coming from an external CSV file: http://pastebin.com/raw.php?i=iAjV88KQ For instance NPA not appearing in the CSV should be not colored (white) and NPA coming few time should appear in light blue and the more frequent the NPA is the more dark blue should be the map. According to my CSV NPA 1203 should be the darkest. Is it possible to obtain what I am looking for with Quantum GIS? (I tried to use "Joins" function in proprieties but it's not working) ![enter image description here](http://i.stack.imgur.com/1ygmV.png) **Edit:** I try to join like this : ![enter image description here](http://i.stack.imgur.com/HP9dV.png)


title: Practical Building Extraction Algorithm Implementation

text: I've been reading many papers over 3D building extraction using LiDAR data and Aerial images. I see that there are many different algorithms for doing such things, such as this and this. Those papers, as far as I can see, describe mostly abstract concepts and math. I am a software developer, so how could I practically implement said algorithms into a piece of software? Also, what algorithm is best for extracting building data with minimal vegetation interference?


title: Where to obtain free 15m ASTER image for New Mexico?

text: I am looking for a free Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) image that will let me to perform NDVI for the pilot crop health/soil moisture that can go as far back ten years or so. It needs to have:   * 15 m (meters)   * temporal   * Ten years back or so   * Perform NDVI This will be in a small area somewhere near Albuquerque, NM.


title: Mimimum distance between points (Postgres/PostGIS)

text: I have a table (tab_Cities with cities which own more than 50000 inhabitants. The table also contains lat, long and the_geom as columns. I can already see the points via QGIS. Now I want to write a postgresql query which gives me the minimum distance of the points but I only now how to get the distance of two points. But now I have around 60 points and I want to have the smallest distance of ALL points. How can I do this? I found some similar questions in other threads but they were using two tables. I am using one table only. Thank you for help! :-).


title: PyQgis question: How to change or set active layer?

text: Developing a qgis plugin, I found that if no layer is selected in the Layers panel I can iterate over all layers and do anything to them, but "activate" them. Could someone give me a hint please? Thank you!


title: Read layer symbology class color with ArcPy

text: I am using ArcPy within ArcGIS 10.1 and attempting to read symbology from a *.lyr file created by a raster dataset within ArcMap. My goal is to be able read out the RGB (color) values assigned to the class breaks on the dataset. It’s seems easy enough to get the lyr file into an object with:                   lyrFile = arcpy.mapping.Layer("C:/.../ReadLyrFileSymbology/conus_foresttype6.img.lyr")      However, after that I am unable to find the method or property to get symbology color (RGB) on a layer. Any pointers would be greatly appreciated.


title: Get Cell Value Webservice for ArcGIS Raster (LiDAR) without using ArcGIS

text: I want to expose the getCellValue function on Lidar Data, as a web service. But I don't have access to ArcGIS (server or desktop). What options do I have to fetch cell value, given lat lon, from a raster in ArcGIS format? It would be great if someone can point me to some open source code/library that I can use to create a web service. I am open to any other option except ArcGIS. I know that Geoserver has a raster rendering feature: can it also find the cell value for a lat-long input?


title: Convert Lambert Conformal Conic to WGS84

text: I am trying to convert my "Lambert Conformal Conic map with 2 standard parallels" to the basic/standard WGS84 projection that Google Maps use. I have 30 of these maps which I need to convert. They are all in `.tiff` format with a corresponding `.tfw` file each. What my final goal is:   1. Convert Maps to WGS84   2. Merge all 30 maps   3. Change them to tiles For now I am only struggling with the conversion of projection. Does anybody now of any ways/programs/techniques or anything (preferably free) to accomplish this? _I have heard/read about GDAL etc. but I do run Win Vista, not MAC._


title: Geowebcache "seed" doesn't seem to work

text: I'm using Geowebcache seed command (via web interface at `http://localhost:8080/geoserver/gwc/rest/seed/workspace:layerName`) and it looks like everything runs smoothly but actually no tiles are saved, at least where I'd expect them to be saved `(/data_dir/gwc)`. I have GeoServer installed on top of Tomcat 7 as WAR. I also checked the logs (enabling `GEOTOOLS_DEVELOPER` mode) and it doesn't give any errors. Actually if enabling `VERBOSE` logging mode I also get something like...               ...         2014-07-07 09:38:08,766 DEBUG [wms.map] - Writing png image ...     2014-07-07 09:38:08,771 DEBUG [wms.map] - Writing png image ... done!     2014-07-07 09:38:08,772 DEBUG [wms.map] - Writing png image ...     2014-07-07 09:38:08,776 DEBUG [wms.map] - Writing png image ... done!     2014-07-07 09:38:08,777 INFO [seed.SeedTask] - GWC Seeder Thread-1 completed (re)seeding layer workspace:layerName after 448 tiles and 45 seconds.     ...      ...and it look like everything went fine. But in the end I have no tiles. By the way, the tiles I'm trying to save to EPSG:900913 projection coming from EPSG:3045 PostGIS store (already done this operation successfully on another VM having GeoServer installed on top of Tomcat6 via `apt-get` instead than via WAR). Any suggestion on how to inspect/solve the issue? **Update:** Following @EvilGenius advice I was able to find the missing tiles (in fact they quite strangely were served by the Geowebcache to the frontend) in `/tmp/tomcat7-tomcat7-tmp/geowebcache/workspace_layerName`. Now the question becomes: how comes they are not copied/moved to the final location (`data_dir/gwc` folder)? again, how can I inspect that?


title: Unable to locate PostGIS in OpenGeo Suite using command line

text: I am working on this new ubuntu 12.04 server and I need to work on the Postgres database which is PostGIS enabled. The problem is I am not sure how to locate the database using command line. The server already has OpenGeo suite installed on it and so it definitely has PostGIS. It is already installed on the server. I just need to locate it. How do I locate it?


title: Annoying little line in color ramp legends - how do I remove it?

text: This was never a problem with ArcGIS 9 but it appears that in 10.0 Esri has added small horizontal lines as the default when including color ramps in legends. I want to get rid of these annoying lines, but I haven't figured out how to. Anyone else run into this problem?