title: GPS Coordinates to WKID 102100

text: I'm currently creating a project for ios and I'm using the latest ArcGIS IOS api. I have a mapView that is using the spatial reference skid 102100 and i am receiving lat long coordinates from an external GPS device (NMEA 0183 messages). I believe NMEA 0183 uses WGS84. I want to display the location from the GPS to my mapView. I tried creating a AGSPoint based on the GPS lat long with WGS84 spatial reference and projecting it to my mapView's spatial reference but the point seems to be at the middle of my map, which i assume is wrong because im in Canada. Can anybody help me pls? I think the reason could be that my mapView has a unit of measure of Meters and the GPS is sending dddmm.mmmmmmm data (d=degrees, m=minutes). here is my projection code:               AGSPoint *tempPoint = [[AGSPoint alloc] initWithX:longitudeValue y:latitudeValue spatialReference:[AGSSpatialReference wgs84SpatialReference]];         AGSPoint *projectedPoint = (AGSPoint *)[self.geometryEngine projectGeometry:temp2 toSpatialReference:self.mapView.spatialReference];      Do still need to convert the lat/long to meters before i project it? if yes, how? Thanks.


title: How to restore a back up file into pgadmin III

text: My backup file containing about 50-60 tables , but I am able to just load only 6 tables for a new db and 31 tables for postgis dbase . And this was the error when I restore a backup file.               pg_restore: [archiver (db)] Error while PROCESSING TOC:     pg_restore: [archiver (db)] Error from TOC entry 1393; 0 0 SHELL TYPE box2d postgres     pg_restore: [archiver (db)] could not execute query: ERROR:  type "box2d" already exists         Command was: CREATE TYPE box2d;      I should I resolve this problem. I am using PostgreSQL 9.1 and please help me how can I load all tables without any errors.my backup file contains 90 MB.


title: "could not get any sources: no image returned from source WMS" error in mapproxy demo

text: I am trying to connect my mapserver to MapProxy. First part of the map file is shown below               MAP     FONTSET "fonts.lst"     IMAGETYPE png     MAXSIZE 4000     SIZE 800 800     EXTENT -20000000 -20000000 20000000 20000000     UNITS meters     IMAGECOLOR 153 179 204     WEB        METADATA           "ows_enable_request" "*"           "wms_srs" "EPSG:4326 EPSG:3857 EPSG:2154 EPSG:310642901 EPSG:4171 EPSG:310024802 EPSG:310915814 EPSG:310486805 EPSG:310702807 EPSG:310700806 EPSG:310547809 EPSG:310706808 EPSG:310642810 EPSG:310642801 EPSG:310642812 EPSG:310032811 EPSG:310642813 EPSG:2986"           "labelcache_map_edge_buffer" "-10"           "wms_title" "osm france"        END        IMAGEPATH "/tmp/ms_tmp/"        IMAGEURL "/ms_tmp/"     END     DEBUG 1     CONFIG "MS_ERRORFILE" "stderr"     CONFIG "PROJ_LIB" "/home/mwh/osm-demo/basemaps"     PROJECTION        "init=epsg:3857"     END     LAYER        TYPE POLYGON        STATUS ON        GROUP "default"        DEBUG 1        NAME "land0"        PROCESSING "APPROXIMATION_SCALE=full"        PROJECTION            "init=epsg:4326"        END        MINSCALEDENOM 332808204        MAXSCALEDENOM 99999999999        DATA "data/TM_WORLD_BORDERS-0.3.shp"        CLASS            STYLE                COLOR 242 239 233            END        END     END     LAYER        GROUP "default"        STATUS ON        MAXSCALEDENOM 99999999999        MINSCALEDENOM 332808204        TYPE LINE        NAME "borders0"        PROJECTION          "init=epsg:4326"        END        DATA "data/boundaries.shp"        CLASS           STYLE              WIDTH 0.5              COLOR "#CDCBC6"           END        END     END      I followed the Tutorial given in MapProxy website to start my proxy by starting the proxy. The tutorial example works fine. But when I use my own .map file it shows an error " **could not get any sources: no image returned from source WMS** " My mapproxy.yaml is shown below               services:       demo:       wms:         md:           title: MapProxy WMS Proxy           abstract: This is the fantastic MapProxy.           online_resource: http://mapproxy.org/           contact:             person: Your Name Here             position: Technical Director             organization:              address: Fakestreet 123             city: Somewhere             postcode: 12345             country: Germany             phone: +49(0)000-000000-0             fax: +49(0)000-000000-0             email: info@omniscale.de           access_constraints:             This service is intended for private and             evaluation use only. The data is licensed             as Creative Commons Attribution-Share Alike 2.0             (http://creativecommons.org/licenses/by-sa/2.0/)           fees: 'None'          sources:       test_wms:         type: wms         req:           url: http://osm.omniscale.net/proxy/service?           layers: osm            mapserver:         type: wms         supported_srs: ['epsg:3857', 'epsg:4326']         req:           url: http://localhost/cgi-bin/mapserv?           layers: default           map: /home/mwh/osm-demo/basemap/osm-google.map           transparent: true         coverage:           bbox: -20000000.0,-20000000.0,20000000.0,20000000.0           bbox_srs: 'epsg:3857'          layers:       - name: mapserver         title: Mapserver Test Layer         sources: [mapserver]      Currently I am stuck in this problem please help me to find some answer. Thank you.


title: Get Conflict Value from ConflictWindow

text: When I reconcile a version in Arc Map, if it has any conflicts it will stop the reconcile operation and open a window. The window shows the conflicts of that particular version with red dots. My task is to read those values and generate a report which contains only conflict fields.


title: How do I create a line between two points with a common field?

text: Does anyone know how I could create a line between two points that have a common code? For example: I have the start and end point of a pipe as a lat/long for a network of pipes, each with a different ID. I want to create a line feature for each of those with the pipe ID which I can then calculate the length from. Anyone know of any arcscript tools in ESRI software or QGIS tools I could use to accomplish this without coding? Thanks!!


title: How to configure GeoWebCache that it can generate tile files conform Slippy Map format

text: The tile files generated by GeoWebCache using gridset 900913 are like: \EPSG_90013_6\02_02\0032_42.png, etc. How can I configure GeoWebCache that it is generating the files conform slippy map format? That is \6\32\21.png, etc.


title: My map is little bit off

text: I am pretty new in GIS things and i have been learning how to use QGIS. Most of the maps are fine, but yesterday I went to GIS center at my University and I got the maps of my country that I will be using for my study. But when I add them in QGIS they are not on the right spot. I checked the coordinate system and it is WGS 84. when I add the species, they are on the right place and map is little bit diagonally towards the right upper corner. I checked the species with the other map that I downloaded from internet and they are on the right place. Also the species and the other map that I downloaded from internet are WGS 84. Any suggestions how can I correct this?


title: Best practice for aerial photo stitching and georeferencing

text: I need to georeference 53 air photos from 1981 (1:25,000 scale) and was wondering what would be the most efficient, defensible, and accurate option. For context, the air photos are following a highway taken along a valley bottom and therefore the relief would not be significant as would be the case if mountains were involved. Essentially, I want to quantify the change in surface area around this highway from 1981 photos to current day; how do I process the former images? Typically, one would georeference each air photo using current satellite imagery. However, in my case, I feel this would be difficult because the area of interest has changed through time, and therefore I fear I will not have enough reference points from current satellite imagery to georeference each air photo from 1981. Therefore, I thought I could first stitch the 53 air photos together using control points from the adjacent air photo taken the same day; essentially creating a mosaic from 1981. Then I could georeference the entire mosaic using current satellite imagery; I think on a larger scale I could find enough reference points to accurately geofrerence the mosaic. Would this process cause too much local distortion that I would no longer be able to accurately quantify surface area? Lastly, I'm curious if orthorectifying the image would significantly improve the surface estimation? I'm not looking for complete precision, but would like it to be accurate enough for a publishable manuscript.


title: How do I change graduated color class values with python?

text: I am trying to use python to adjust my graduated color symbols' class breaks and labels. The layer is a query layer tied to an SQL DB and I want to be able to execute my python (after the database has been updated) to automatically author a map and have the layer/legend reflect the changes made in the database. I've tried this script but it doesn't work and I get an error I don't recognize.               import arcpy     mxd = arcpy.mapping.MapDocument("current")     df = arcpy.mapping.ListDataFrames(mxd, "Layers")[0]     lyr = arcpy.mapping.ListLayers(mxd, "Unsafe_Driving", df)[0]     if lyr.symbologyType == "GRADUATED_COLORS":         lyr.symbology.valueField = "CSA.DBO.%Unsafe_Driving_1.VARIANCE"         lyr.symbology.classBreakValues = [-.01, .01, .0866]         lyr.symbology.classBreakLabels = ["-8.66% to -1.00%", "-.99% to 1.00%","1.01% to 8.66%"]           arcpy.RefreshActiveView()      > Runtime error Traceback (most recent call last): File "", line 7, in File > "c:\program files (x86)\arcgis\desktop10.1\arcpy\arcpy\arcobjects_base.py", > line 87, in _set return setattr(self._arc_object, attr_name, cval(val)) > RuntimeError The goal here is to have 3 classes:   1. Lowest value to -1.00%    2. -.99% to 1.00%   3. 1.01% to the highest value.  Since the lowest and highest values can change, do I need to use a query to dynamically insert those into the lyr.symbology.classBreakValues? Any suggestions or code examples would help tremendously.


title: Return features based on users credentials Geoserver WFS

text: We are using LDAP based authentication with geoserver which is connected to a PostGIS database. Users insert rows (features) into the table (layer) using WFS-T after authentication. Now when users query for data from this layer (also after authentication), we would like to return data that is applicable to this user based on their credential (i.e) only view data that was inserted by that user or group. We are storing the username in the table. We get the username because this is defined as non-null column in our table and users are required to provide it as a part of their WFS-T insert. I read about service level security and layer level security, but what we need seems more like row level or feature level security. We can use JDBC based authentication instead of LDAP if needed, but do not know if this will solve the problem. I am pretty comfortable with Java and some what knowledgeable on spring security and hence can deal with customizations if needed. Please advice.


title: going from hillshade to percent shaded (gdaldem)

text: I assumed hillshades output percent shaded (scaled to 255) at the given azimuth and altitude but this is incorrect. A flat raster at 80 degrees altitude has a maximum value of 251 instead of 255. How can I go from a hillshade to percent shaded and can I simply divide by the cosine of the 10 degrees in this case to get the correct value or is it more complicated? Here is some `R` code with `gdaldem` that may help you see what I'm seeing:               library(raster)     r <- raster(ncol=100,nrow=100)     values(r) <- rep(1,ncell(r))     writeRaster(r,"flatelevation.tif")     system("gdaldem hillshade -az 180 -alt 80 flatelevation.tif hillshade_flat.tif")     hs <- raster("hillshade_flat.tif")     summary(values(hs))     summary(values(hs)/cos((90.0-80.0)*pi/180.0))      I will try to parse the source code but it's been ten years since I've used C++... In the meantime if someone with general knowledge about hillshades can point me the right direction it would be greatly appreciated. ADDED: So I thought I could get a "percent of the tile that is shaded" from running a hillshade. Upon further reading a hillshade outputs the illuminance of each tile. So, the cosine is only part of the equation, I would need to pull in the slope and aspect as well for additional trigonometric readjustments, and maybe more. Turns out there are tools that are functionally equivalent for that I need called "shadow masks" that only output 1 or 0, indicating whether the tile is shaded. This nail biting article referenced in `gdal`'s source code helped enlighten me, in addition to this question on cloud shading, this question on solar maps, and this question on r.sun where a contributor indicates that r.sunmask is slow. A couple available tools: GRASS r.sunmask RSGISLib Elevation Module: shadowmask USGS ISIS Shadow Of course I need the speed of `gdal`'s hillshading so I'll keep this question up until I or others come up with a solution that can process a 1000 square mile area in a second.


title: I would like to reverse Geocode from the infowindow

text: I am just learning the ropes with CartoDB and am curious to know if this is possible. I have a CartoDB map with points and I am entering data into this map through the CartoDB editor. I would like to be able to Google reverse geocode the street address of my point automatically. Ideally I would create a link on the infowindow labeled 'Find address' this would then run a script which would find the address of that point and INSERT it into that points specified column on my table. I would like this to run entirely from the CartoDB Editor as this is where I update my data from. Does anyone know if this is possible and if so where is the best place to start? If this won't work are there any other reverse Geocode solutions out there?


title: How to solve jpeg2000, jp2. crash & errors on arcgis 9.3.1 windows 7 64-bit

text: I have several problems at work. One is that jpeg2000 crashes ArcMap but mostly ArcCatalog when trying to read them. This is only occurs on machines with Windows 7 64-bit (we don't have Windows 7 32-bit). What can cause this? I have deleted templates, re-installed, everything in my power. But there is obviously some compatibility issue, because it works great on Windows XP.


title: Buttons (Ok, Cancel) disappear in "Apri Modulo" from Attribute Table in QGIS 2.0.1

text: QGIS 2.0 has changed the visualisation of a form opened from the attribute table: instead of opening the form in a separate window, the form is opened in the same one. This is not a problem, **but all the buttons disappear...** I created my own form with Qt4 and all the logic with PyQt and the majority of the constraints are associated to the Ok button: now my logic is no longer useful. It is possible to re-enable these buttons? Thank you, Elisabetta **- Update -** I'm trying to do this, for the moment I have no update... I copy here part of my code, maybe someone can help me.               from PyQt4.QtCore import *       from PyQt4.QtGui import *       from definizioni import *          myDialog = None     denominazione = None     denominazioneCheck = None          def formOpen(dialog,layerid,featureid):       global myDialog       myDialog = dialog            global denominazione, denominazioneCheck       denominazione = myDialog.findChild(QLineEdit,"denominazione")       denominazione.textChanged.connect(denominazione_onTextChanged)       if (denominazione.text() == "NULL"):         denominazione.setText("")       else:         temp = denominazione.text()         denominazione.setText("null")         denominazione.setText(temp)            buttonBox = myDialog.findChild(QDialogButtonBox,"buttonBox")       # Disconnect the signal that QGIS has wired up for the dialog to the button box.       buttonBox.accepted.disconnect(myDialog.accept)       # Wire up our own signals.       buttonBox.accepted.connect(validate)       buttonBox.rejected.connect(myDialog.reject)       self.connect(buttonBox, SIGNAL("accepted()"), self.accept)          def denominazione_onTextChanged(text):       global denominazioneCheck            if not validateStringNotNull(text):         denominazione.setStyleSheet("background-color: rgba(255, 107, 107, 150);")         denominazioneCheck = -2       else:         if not validateAlpha(text):           denominazione.setStyleSheet("background-color: rgba(255, 107, 107, 150);")           denominazioneCheck = -1         else:           denominazione.setStyleSheet("")           denominazioneCheck = 1          # Valido la form     def validate():       if denominazioneCheck == -2:         messageBoxErrore(myDialog,"La denominazione di un comune non puo' essere vuota")       else:         if denominazioneCheck == -1:           messageBoxErrore(myDialog,"La denominazione di un comune non puo' contenere numeri")         else:           # Return the form as accepted to QGIS.           myDialog.accept()      I need that the field of the form respect some constraints, as you can see. **In case that constraints are not fulfilled I don't want that the record is stored in the database (I don't want accept() signal is called!).**


title: OpenLayers programmatically select base layer

text: I have several base layers added to a map in an OpenLayers based application. There is only one that can be displayed at a time, as they are all base layers. When the user changes the selected base layer, I store the name of the selected layer as a HTML5 local storage item (an improved cookie). When re- opening the application, I want to programmatically set the same base layer be displayed. It works fine with `layer.visibility = true` if `isBaseLayer` is `false`, but not in the case when `isBaseLayer` is set to `true`. Is there anything else I could try for this?


title: ArcMap Text Annotations Display in Proportion to Map Scale

text: My question is probably of an 'easy fix' however I've been struggling with it for a while! I'm trying to set up an annotation labels for several feature classess in a map and would like them to display proportionally - relatively to map scale. I use ArcGIS v 10.1 and Maplex label engine. So - as an example, I am working on a regional map at a scale of 1:700K. My annotation on a FC (within GDB) is set to 12pt size. When I zoom in to a larger scale of 1:350K my annotations appear exactly the same size and are consequently tooo big! I understand one can work with map document annotation or geodatabse annotation. I've seen some posts about setting a data frame reference scale to zero (eg http://forums.arcgis.com/threads/74004-Annotation-labels-at- different-scales), which should make mxd annotations display proportionally to map scale. However cannot seem to add custom reference scale of '0' to the customized list. I'd greatly appreciate your comments on best practices or any other advices you may have! Many thanks, Magda


title: Convert polyline to presence absence raster in arcpy

text: I have a polyline that I am converting to a raster in arcpy. This is occuring within a loop and at each iteration I change the processing environment to a different area of interest. In this case it is a specific polygon (like below)               whereClause = '"FID" = ' + str(poly.FID) # Create a clause to select only the current record     arcpy.MakeFeatureLayer_management(shapefile, "currentPolygon", whereClause)             arcpy.env.extent = "currentPolygon"      At each iteration I would like to create a raster where cells containing the polyline are given a value of 1 and cells without polyline are given a 0 within the bounds of the env.extent. What is the key to creating the 0 values? Does it come from the arcpy.env.nodata parameter? I can make the conversion in my loop but cannot quite get the binary format I am envisioning.               arcpy.FeatureToRaster_conversion("roadsLayer", "VALUE", "roadsGrid.tif",)


title: ArcObjects 10.1 IDatabaseCompact lock file

text: I'm trying to automate one of task by creating code that will execute the IDatabaseCompact.Compact method to a GDB file. Everything seems to work fine but there is always a .lock file generated after the Compact method is called. I am planning to do something else after the Compact process but it errors out because of the lock file. Is there any way to close/delete the lock file? BTW, this is how my process works: first, i bind and initialize my system using: ` ESRI.ArcGIS.RuntimeManager.Bind(ESRI.ArcGIS.ProductCode.EngineOrDesktop) _ao = New ESRI.ArcGIS.esriSystem.AoInitialize() _ao.Initialize(ESRI.ArcGIS.esriSystem.esriLicenseProductCode.esriLicenseProductCodeBasic) ` I then call my compact method:               Public Sub CompressFileGDB(gdbFilePath As String)         Dim dbCompact As IDatabaseCompact         Dim ws As IWorkspace = FileGdbWorkspaceFromPath(gdbFilePath)              If (TypeOf ws Is IDatabaseCompact) Then             dbCompact = CType(ws, IDatabaseCompact)             If (dbCompact.CanCompact) Then                 Dim didFinish As Boolean = True                 Try                     dbCompact.Compact()                 Catch ex As Exception                     didFinish = False                 Finally                     dbCompact = Nothing                     ws = Nothing                     _ao.Shutdown()                          RaiseEvent OnProcessDone(didFinish)                 End Try             End If         End If     End Sub      The compact method works, there is a significant change in the GDB file size. I also tried to set my attributes to Nothing after the process but it doesn't seem to help that much.


title: Merging h5 files from VIIRS to gtiff

text: I am trying to mosaic all the data from VIIRS of a day into a single gtiff. I've already done it with MODIS data using GDAL. The problem it seems is that the coordinates of each h5 files doesn't seem to be in the same file as the data. Which gives this as an error when I try a gdalwarp: > ERROR 1: Unable to compute a transformation between pixel/line and > georeferenced coordinates for > > > HDF5:SVI01_npp_d20140812_t0007456_e0013260_b14451_c20140812061326517212_noaa_ops.h5://All_Data/VIIRS-I1-SDR_All/Radiance. > > There is no affine transformation and no GCPs. Do you have an idea what could be done to get the localisation data and to merge them together ?


title: How to overlay a county boundary on a raster file?

text: I want to overlay California county boundary top of raster layer. Here is is my code:               data(coastsCoarse)     data(countriesLow)     library(gtools)     library(maps)     library(maptools)     library(sp)          # Raster file     RAD2012 <- raster(paste(getwd(), "/2012bil/2012ASC5min_ppt_spas1306_0001_20121222_0005_UTC.asc.bil", sep = ""))          list.ras <- mixedsort(list.files(paste(getwd(), "/2012bil/", sep = ""), full.names = T, pattern = ".asc.bil"))           # stack raster file     RAD2012.all <- stack(list.ras)          # png file tile(time stamp)     list.files <- mixedsort(list.files(paste(getwd(), "/2012bil/", sep = ""), full.names = F, pattern = ".asc.bil"))           outlist <- substr(list.files, 31, 43)          # plotting part     # I am using countriesLow     # I want to improve base layer from countriesLow to county boundary      ![enter image description here](http://i.stack.imgur.com/LM6Ty.png)               for(i in seq(from=0, to=length(list.ras)-1, by=1)){           png(paste(outlist[[i+1]],".png",sep=""))         par(mfrow=c(1,1),oma=c(2,2,2,2))         for(j in 1:1){               plot(RAD2012.all[[i+j]])                   plot(countriesLow, add = T)                   #plot(countriesHigh, add = T)                   title(outlist[[i+j]])         }           dev.off()     }


title: Dealing with overlapping features in OpenLayers

text: If users create a polygon feature that completely overlaps another one, like this … ![](http://i.stack.imgur.com/He9pi.png) … then there's no way you can select the inner feature, i.e. the one that's overlapped by the other. I'm currently using a simple `SelectFeature` like this:               var select = new OpenLayers.Control.SelectFeature(       [this.annotationLayer, this.controlPointsLayer], {          hover: true,         renderIntent: "temporary"         }     );     this.map.addControl(select);     select.activate();      But obviously, I can only get it to select the outer feature. What is the typical way of dealing with this? **How can I make it possible for the user to select the inner feature as well?** The simplest solution would probably be to make the outer feature only selectable within the area that's left when subtracting the inner feature – but I have no idea how to go about doing that.


title: Legend of a raster map with categorical data

text: I would like to plot a raster containing 4 different values (1) with a categorical text legend describing the categories such as 2 but with colour boxes. I've tried using legend such as :               legend( 1,-20,legend = c("land","ocean/lake", "rivers","water bodies"))      but I don't know how to associate one value to the displayed color. Is there a way to retrieve the colour displayed with 'plot' and to use it in the legend? Initial raster![enter image description here](http://i.stack.imgur.com/uGuph.png) Raster with legend ![enter image description here](http://i.stack.imgur.com/9wxw4.png)


title: an efficient way to find the pgrouting edges contained within an isochrone

text: So I am working with some pgrouting road network data in postgis2.0/postgresql 9.1. I am interested in routing around certain types of disasters like tornados or such. If a disaster event happens at some spatial coordinate `x1,y1`. Then I would generate an isochrone from that point to every other point that is reachable at some cost `c`. So basically I am creating a travel buffer around the event. Then I would like to loop through all of the road edges that fall within this isochrone and increase the costs of travel along that road segment. I am trying to figure out how to reconcile the pgrouting part with the postgis part--I have less familiarity with pgrouting. Would I create a polygon out of the isochrone and then do some sort of `contains` or `within` query to get the included edges, and then run an sql update to increase their costs. An additional consideration is the efficiency of this operation since this type of spatial operation might present a problem for a large number of network edges.


title: Unknown error concurring when developing QGIS 2.0.1 plugin in Windows 7

text: I am following the tutorial found here for developing a simple QGIS plugin within QGIS 2.0.1 installed on Windows 7. However, as I have began to try developing a link between the GUI and each the mouse click of the map canvas, the main plugin code does not seem to recognize the fields in the GUI. To illustrate with the example in the tutorial. Within in my _plugin.py_ code I have the function:               class vector_selectbypoint:          def __init__(self, iface):              # Save reference to the QGIS interface              self.iface = iface              # refernce to map canvas              self.canvas = self.iface.mapCanvas()              # out click tool will emit a QgsPoint on every click              self.clickTool = QgsMapToolEmitPoint(self.canvas)              # create our GUI dialog              self.dlg = vector_selectbypointDialog()                  def handleMouseDown(self, point, button):              self.dlg.clearTextBrowser()              self.dlg.setTextBrowser( str(point.x()) + " , " +str(point.y()) )      that simply writes the lat and long into a TextBrowser with the GUI by calling the `clearTextBrowser` and `setTextBrowser` functions in a seperate _pluginDialog.py_ script where the following code is ran:               class vector_selectbypointDialog(QtGui.QDialog, Ui_vector_selectbypoint):         def __init__(self):             QtGui.QDialog.__init__(self)             self.ui = Ui_vector_selectbypoint()             self.setupUi(self)              def setTextBrowser(self, output):             self.ui.txtFeedback.setText(output)              def clearTextBrowser(self):             self.ui.txtFeedback.clear()      which interacts with the _ui_plugin_ui.py_ where `txtFeedback` is the TextBrowser in the GUI. Below is a snippet of the _ui_plugin_ui.py_ file:               class Ui_vector_selectbypoint(object):         def setupUi(self, vector_selectbypoint):             vector_selectbypoint.setObjectName(_fromUtf8("vector_selectbypoint"))             vector_selectbypoint.resize(400, 300)             self.buttonBox = QtGui.QDialogButtonBox(vector_selectbypoint)             self.buttonBox.setGeometry(QtCore.QRect(30, 240, 341, 32))             self.buttonBox.setOrientation(QtCore.Qt.Horizontal)             self.buttonBox.setStandardButtons(QtGui.QDialogButtonBox.Cancel|QtGui.QDialogButtonBox.Ok)             self.buttonBox.setObjectName(_fromUtf8("buttonBox"))             self.chkActivate = QtGui.QCheckBox(vector_selectbypoint)             self.chkActivate.setGeometry(QtCore.QRect(20, 250, 81, 20))             self.chkActivate.setObjectName(_fromUtf8("chkActivate"))             self.txtFeedback = QtGui.QTextBrowser(vector_selectbypoint)             self.txtFeedback.setGeometry(QtCore.QRect(80, 20, 221, 151))             self.txtFeedback.setObjectName(_fromUtf8("txtFeedback"))      The issue is that when I run the plugin in QGIS and click on the map canvas the following error is thrown:               AttributeError: 'ui_plugin_ui' object has no attribute 'txtFeedback'      I have verified that there are no syntax issues but the issue persists. If anyone has any insight I would eternally grateful!


title: Enterprise License Server Management

text: I was hoping to get some information from others who have deployed an enterprise GIS License Manager for ArcGIS 10.x software. Honestly, I would like to see if our current setup is adequate or if we should venture into a new setup architecture. Currently, we have multiple departments who each purchase their own licenses. To ensure departments only have access to their purchased licenses, ArcGIS License Managers have been installed for each department and their licenses authorized on each LM respectively. Some of these LM's are on their own dedicated virtual server, and some are installed on client machines in the department (not ideal and we plan to change this immediately). Our idea to both keep the licenses separated, and also not use any of our farm resources is to place physical boxes (ones that have been replaced via PC upgrades, etc.) in each department, install LM on the box and authorize the licenses. Each department would have it's own box, and we can administer these via remove access. Here are questions I have:   1. Is it at all possible to run multiple instances of LM on one machine, and just change the outgoing port so each instance has it's own dedicated I/O access? (i.e. Dept A uses 27000, Dept B uses 27001, etc.)   2. I have heard (and read somewhere), that individual licenses on LM's can be managed down to the user level via changing the LM configuration file. If this is so, have any of you implemented this setup and how successful has it proved to be? Again, we are just wanting to ensure we are implementing a practical architecture. We've used the above setup for years so I believe it's time we revisit this and see if better practices exist. Here are some of the posts/articles we have looked over for ideas so far:   * Managing ArcGIS floating licenses   * hardcoding port number is ArcGIS License Manager Thanks in advance for your input!


title: Finding a Conceptual Site Modeling Software

text: I've been asked by my company to generate a conceptual site model for a project and I have no clue where to start. The purpose is to show groundwater flow through a few layers of rock either in 2D or 3D. I don't need a program like MODFLOW, just something to visualize it. My search always ends at a website that describes how to make a CSM or I read articles with great pictures/figures with no explanation on who/what created it. What I want is to download (doesn't have to be free) software that will help my model look professional. Does anyone have a suggestion for me?


Where can I find in depth documentation on modifying arcgis stylesheets? I have found some of the basics here : Changing tool dialog appearance I want to change the height on one element, the only element in my script. I have a filter value list of 82 counties and I would like to see more than 9 when the tool opens. Usually ESRI documentation is thorough unless it is a topic they expect everyone to know. Is this the case with xsl stylesheets? I think I could figure it out if I just knew what the element was called. does anyone have a suggestion? ![enter image description here](http://i.stack.imgur.com/3Y96V.jpg)


title: Why does GeoServer's GetFeature round my numbers, but only for outputFormat=CSV?

text: When I make a GetFeature request with outputFormat=CSV to my GeoServer v2.4.5 I'm finding that the floating point numbers in my returned data have been rounded. When I make the same request with other output formats the floats are returned correctly. Here's the request returning XML:                http://geothermal.smu.edu:9000/geoserver/gtda-contributions/ows?service=WFS&version=1.0.0&request=GetFeature&typeName=gtda-contributions:UNDRadiogenicHeat&maxFeatures=1      This excerpt from the returned XML shows the data as it was entered into the backing database. There are more floating point numbers in the returned data, but I exclude them for brevity:                <?xml version="1.0" encoding="UTF-8"?>      <wfs:FeatureCollection xmlns="http://www.opengis.net/wfs" xmlns:wfs="http://www.opengis.net/wfs" xmlns:gtda-contributions="http://scr.siemens.com/gtda-contributions" xmlns:gml="http://www.opengis.net/gml" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:schemaLocation="http://scr.siemens.com/gtda-contributions http://geothermal.smu.edu:9000/geoserver/gtda-contributions/wfs?service=WFS&version=1.0.0&request=DescribeFeatureType&typeName=gtda-contributions%3AUNDRadiogenicHeat http://www.opengis.net/wfs http://geothermal.smu.edu:9000/geoserver/schemas/wfs/1.0.0/WFS-basic.xsd">         <gml:boundedBy>             <gml:null>unknown</gml:null>         </gml:boundedBy>         <gml:featureMember>             <gtda-contributions:UNDRadiogenicHeat fid="UNDRadiogenicHeat.2540ea6b-f94a-3088-9b32-762c77a8e57f">                 <gtda-contributions:latitude>47.90752513</gtda-contributions:latitude>                 <gtda-contributions:longitude>-92.79726665</gtda-contributions:longitude>      ...      The same request with outputFormat=application/json:                http://geothermal.smu.edu:9000/geoserver/gtda-contributions/ows?service=WFS&version=1.0.0&request=GetFeature&typeName=gtda-contributions:UNDRadiogenicHeat&maxFeatures=1&outputFormat=application/json      returns similar (correct) data:                {     "type" : "FeatureCollection",     "totalFeatures" : 343,     "features" : [{             "type" : "Feature",             "id" : "UNDRadiogenicHeat.2540ea6b-f94a-3088-9b32-762c77a8e57f",             "geometry" : {                 "type" : "Point",                 "coordinates" : [-92.79726665, 47.90752513]             },             "geometry_name" : "shape",             "properties" : {                 "latitude" : 47.90752513,                 "longitude" : -92.79726665,     ...      However, if the request is for CSV:                http://geothermal.smu.edu:9000/geoserver/gtda-contributions/ows?service=WFS&version=1.0.0&request=GetFeature&typeName=gtda-contributions:UNDRadiogenicHeat&maxFeatures=1&outputFormat=csv      the same portion of data has been rounded off (48 instead of 47.90752513, and -93 instead of -92.79726665):                FID,latitude,longitude,...      UNDRadiogenicHeat.2540ea6b-f94a-3088-9b32-762c77a8e57f,48,-93,...      I've looked for configurations that might cause this, but I see none. I know there are filters that could do this, but I get the impression that these need to be explicitly requested, and the above requests do not include filters. Is there something I'm overlooking? Is there some way that GeoServer is treating the data so glaringly differently just because I'm requesting CSV? Thank you for any insights.


when I try to preview a layer the exception is throw. Internal error occurred javax.imageio.IIOException: Improper call to JPEG library in state 202 Improper call to JPEG library in state 202 logg: org.geoserver.platform.ServiceException: Internal error occurred at org.geoserver.wms.GetFeatureInfo.run(GetFeatureInfo.java:49) at org.geoserver.wms.DefaultWebMapService.getFeatureInfo(DefaultWebMapService.java:346) at sun.reflect.GeneratedMethodAccessor348.invoke(Unknown Source) at sun.reflect.DelegatingMethodAccessorImpl.invoke(DelegatingMethodAccessorImpl.java:43) at java.lang.reflect.Method.invoke(Method.java:606) at org.springframework.aop.support.AopUtils.invokeJoinpointUsingReflection(AopUtils.java:319) at org.springframework.aop.framework.ReflectiveMethodInvocation.invokeJoinpoint(ReflectiveMethodInvocation.java:183) at org.springframework.aop.framework.ReflectiveMethodInvocation.proceed(ReflectiveMethodInvocation.java:150) at org.geoserver.kml.WebMapServiceKmlInterceptor.invoke(WebMapServiceKmlInterceptor.java:34) at org.springframework.aop.framework.ReflectiveMethodInvocation.proceed(ReflectiveMethodInvocation.java:172) at org.geoserver.ows.util.RequestObjectLogger.invoke(RequestObjectLogger.java:54) at org.springframework.aop.framework.ReflectiveMethodInvocation.proceed(ReflectiveMethodInvocation.java:172) at org.springframework.aop.framework.JdkDynamicAopProxy.invoke(JdkDynamicAopProxy.java:202) at com.sun.proxy.$Proxy24.getFeatureInfo(Unknown Source) at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method) at sun.reflect.NativeMethodAccessorImpl.invoke(NativeMethodAccessorImpl.java:57) at sun.reflect.DelegatingMethodAccessorImpl.invoke(DelegatingMethodAccessorImpl.java:43) at java.lang.reflect.Method.invoke(Method.java:606) at org.geoserver.ows.Dispatcher.execute(Dispatcher.java:774) at org.geoserver.ows.Dispatcher.handleRequestInternal(Dispatcher.java:272) at org.springframework.web.servlet.mvc.AbstractController.handleRequest(AbstractController.java:153) at org.springframework.web.servlet.mvc.SimpleControllerHandlerAdapter.handle(SimpleControllerHandlerAdapter.java:48) at org.springframework.web.servlet.DispatcherServlet.doDispatch(DispatcherServlet.java:923) at org.springframework.web.servlet.DispatcherServlet.doService(DispatcherServlet.java:852) at org.springframework.web.servlet.FrameworkServlet.processRequest(FrameworkServlet.java:882) at org.springframework.web.servlet.FrameworkServlet.doGet(FrameworkServlet.java:778) at javax.servlet.http.HttpServlet.service(HttpServlet.java:707) at javax.servlet.http.HttpServlet.service(HttpServlet.java:820) at org.mortbay.jetty.servlet.ServletHolder.handle(ServletHolder.java:487) at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1093) at org.geoserver.filters.ThreadLocalsCleanupFilter.doFilter(ThreadLocalsCleanupFilter.java:27) at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084) at org.geoserver.filters.SpringDelegatingFilter$Chain.doFilter(SpringDelegatingFilter.java:74) at org.geoserver.wms.animate.AnimatorFilter.doFilter(AnimatorFilter.java:70) at org.geoserver.filters.SpringDelegatingFilter$Chain.doFilter(SpringDelegatingFilter.java:70) at org.geoserver.filters.SpringDelegatingFilter.doFilter(SpringDelegatingFilter.java:45) at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084) at org.geoserver.platform.AdvancedDispatchFilter.doFilter(AdvancedDispatchFilter.java:49) at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084) at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:311) at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:68) at org.springframework.security.web.access.intercept.FilterSecurityInterceptor.invoke(FilterSecurityInterceptor.java:116) at org.springframework.security.web.access.intercept.FilterSecurityInterceptor.doFilter(FilterSecurityInterceptor.java:83) at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:72) at org.geoserver.security.filter.GeoServerCompositeFilter.doFilter(GeoServerCompositeFilter.java:91) at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323) at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:68) at org.springframework.security.web.access.ExceptionTranslationFilter.doFilter(ExceptionTranslationFilter.java:113) at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:72) at org.geoserver.security.filter.GeoServerCompositeFilter.doFilter(GeoServerCompositeFilter.java:91) at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323) at org.geoserver.security.filter.GeoServerAnonymousAuthenticationFilter.doFilter(GeoServerAnonymousAuthenticationFilter.java:53) at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323) at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:68) at org.springframework.security.web.authentication.www.BasicAuthenticationFilter.doFilter(BasicAuthenticationFilter.java:150) at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:72) at org.geoserver.security.filter.GeoServerCompositeFilter.doFilter(GeoServerCompositeFilter.java:91) at org.geoserver.security.filter.GeoServerBasicAuthenticationFilter.doFilter(GeoServerBasicAuthenticationFilter.java:82) at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323) at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:68) at org.springframework.security.web.context.SecurityContextPersistenceFilter.doFilter(SecurityContextPersistenceFilter.java:87) at org.geoserver.security.filter.GeoServerSecurityContextPersistenceFilter$1.doFilter(GeoServerSecurityContextPersistenceFilter.java:52) at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:72) at org.geoserver.security.filter.GeoServerCompositeFilter.doFilter(GeoServerCompositeFilter.java:91) at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323) at org.springframework.security.web.FilterChainProxy.doFilter(FilterChainProxy.java:173) at org.geoserver.security.GeoServerSecurityFilterChainProxy.doFilter(GeoServerSecurityFilterChainProxy.java:134) at org.springframework.web.filter.DelegatingFilterProxy.invokeDelegate(DelegatingFilterProxy.java:346) at org.springframework.web.filter.DelegatingFilterProxy.doFilter(DelegatingFilterProxy.java:259) at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084) at org.geoserver.filters.LoggingFilter.doFilter(LoggingFilter.java:75) at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084) at org.geoserver.filters.GZIPFilter.doFilter(GZIPFilter.java:42) at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084) at org.geoserver.filters.SessionDebugFilter.doFilter(SessionDebugFilter.java:47) at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084) at org.geoserver.filters.FlushSafeFilter.doFilter(FlushSafeFilter.java:43) at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084) at org.vfny.geoserver.filters.SetCharacterEncodingFilter.doFilter(SetCharacterEncodingFilter.java:109) at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1084) at org.mortbay.jetty.servlet.ServletHandler.handle(ServletHandler.java:360) at org.mortbay.jetty.security.SecurityHandler.handle(SecurityHandler.java:216) at org.mortbay.jetty.servlet.SessionHandler.handle(SessionHandler.java:181) at org.mortbay.jetty.handler.ContextHandler.handle(ContextHandler.java:726) at org.mortbay.jetty.webapp.WebAppContext.handle(WebAppContext.java:405) at org.mortbay.jetty.handler.HandlerWrapper.handle(HandlerWrapper.java:152) at org.mortbay.jetty.Server.handle(Server.java:324) at org.mortbay.jetty.HttpConnection.handleRequest(HttpConnection.java:505) at org.mortbay.jetty.HttpConnection$RequestHandler.headerComplete(HttpConnection.java:828) at org.mortbay.jetty.HttpParser.parseNext(HttpParser.java:514) at org.mortbay.jetty.HttpParser.parseAvailable(HttpParser.java:211) at org.mortbay.jetty.HttpConnection.handle(HttpConnection.java:380) at org.mortbay.jetty.bio.SocketConnector$Connection.run(SocketConnector.java:228) at org.mortbay.thread.BoundedThreadPool$PoolThread.run(BoundedThreadPool.java:450) Caused by: java.lang.RuntimeException: javax.imageio.IIOException: Improper call to JPEG library in state 202 at com.sun.media.jai.imageioimpl.ImageReadOpImage.computeTile(ImageReadOpImage.java:706) at com.sun.media.jai.util.SunTileScheduler.scheduleTile(SunTileScheduler.java:914) at javax.media.jai.OpImage.getTile(OpImage.java:1129) at javax.media.jai.RenderedOp.getTile(RenderedOp.java:2257) at org.geotools.coverage.grid.GridCoverage2D.evaluate(GridCoverage2D.java:543) at org.geotools.coverage.grid.GridCoverage2D.evaluate(GridCoverage2D.java:468) at org.geoserver.wms.featureinfo.RasterLayerIdentifier.identify(RasterLayerIdentifier.java:178) at org.geoserver.wms.GetFeatureInfo.execute(GetFeatureInfo.java:77) at org.geoserver.wms.GetFeatureInfo.run(GetFeatureInfo.java:43) ... 93 more Caused by: javax.imageio.IIOException: Improper call to JPEG library in state 202 at com.sun.imageio.plugins.jpeg.JPEGImageReader.readImageHeader(Native Method) at com.sun.imageio.plugins.jpeg.JPEGImageReader.readNativeHeader(JPEGImageReader.java:604) at com.sun.imageio.plugins.jpeg.JPEGImageReader.readHeader(JPEGImageReader.java:598) at com.sun.imageio.plugins.jpeg.JPEGImageReader.readInternal(JPEGImageReader.java:1054) at com.sun.imageio.plugins.jpeg.JPEGImageReader.read(JPEGImageReader.java:1034) at com.sun.media.jai.imageioimpl.ImageReadOpImage.computeTile(ImageReadOpImage.java:697) ... 101 more 2014-04-15 17:41:03,734 INFO [geoserver.wms] - Request: getMap Exceptions = SE_XML BgColor = java.awt.Color[r=255,g=255,b=255] Transparent = false SldBody = null ValidateSchema = false Sld = null SldVersion = null TilesOrigin = null FeatureVersion = null SRS = EPSG:3763 Styles = [StyleImpl[ name=raster], StyleImpl[ name=freguesias], StyleImpl[ name=municipios], StyleImpl[ name=distritos]] Layers = [org.geoserver.wms.MapLayerInfo@934b551d, org.geoserver.wms.MapLayerInfo@11d261a9, org.geoserver.wms.MapLayerInfo@5592ff09, org.geoserver.wms.MapLayerInfo@1bd31998] MaxFeatures = null Crs = PROJCS["ETRS89 / Portugal TM06", GEOGCS["ETRS89", DATUM["European Terrestrial Reference System 1989", SPHEROID["GRS 1980", 6378137.0, 298.257222101, AUTHORITY["EPSG","7019"]], TOWGS84[0.0, 0.0, 0.0, 0.0, 0.0, 0.0, 0.0], AUTHORITY["EPSG","6258"]], PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]], UNIT["degree", 0.017453292519943295], AXIS["Geodetic longitude", EAST], AXIS["Geodetic latitude", NORTH], AUTHORITY["EPSG","4258"]], PROJECTION["Transverse_Mercator", AUTHORITY["EPSG","9807"]], PARAMETER["central_meridian", -8.133108333333334], PARAMETER["latitude_of_origin", 39.66825833333334], PARAMETER["scale_factor", 1.0], PARAMETER["false_easting", 0.0], PARAMETER["false_northing", 0.0], UNIT["m", 1.0], AXIS["Easting", EAST], AXIS["Northing", NORTH], AUTHORITY["EPSG","3763"]] Bbox = SRSEnvelope[-35.432430190393 : -2.0255445246072, 28.009986277326 : 44.174608373674] RemoteOwsType = null RemoteOwsURL = null Env = {} FormatOptions = {} Angle = 0.0 CQLFilter = null Elevation = [] FeatureId = null StartIndex = null ViewParams = null Palette = null Tiled = false Buffer = 0 Filters = null Width = 682 Filter = null Time = [] Format = image/png Height = 330 RequestCharset = UTF-8 Get = true BaseUrl =host:8080/geoserver/ RawKvp = {BBOX=-35.432430190393,28.009986277326,-2.0255445246072,44.174608373674, VERSION=1.1.1, FORMAT=image/png, SERVICE=WMS, HEIGHT=330, REQUEST=GetMap, LAYERS=teste_layer, STYLES=, WIDTH=682, SRS=EPSG:3763} Request = GetMap Version = 1.1.1


title: Spatial Interpolation Methods like PRISM

text: I am looking for some kind of software (free software) that allows me to make an interpolation like PRISM software. The interpolation used in PRISM is a > _local regression and accounts for spatially varying elevation > relationships; also accounts for effectiveness of terrain as barriers, > terrain-induced climate transitions, cold air drainage and inversions, and > coastal effects_ What I am aiming to interpolate are climate time-series datasets. source


title: Which pgrouting workshop or tutorial is most current?

text: Please recommend pgrouting tutorials or workshops that are most current. I ask because I just read a post that claims a tutorial is deprecated. (see osm2po vs osm2pgrouting: workshop post) I am working through the pgrouting-workshop (http://workshop.pgrouting.org/) and find it mostly helpful. The OS and environment here: ubuntu 12.04, postgres 9.1, postgis 1.5. While it may be unreasonable to ask that all deprecated tutorials be so marked, it would help if the date of authorship is included prominently in the header of all training materials, then newbies would know to expect some staleness.


title: Programmatically setting ui for QGIS via setEditForm()

text: I've been following this thread in attempts to configure my ui preferences programmatically. I've had no issues with loading styles (which subsequently sets the custom ui if I saved it this way), but I am unable to setEditForm solely. The following things do not work (result is that the properties stay on Autogenerate default option) both in the python console and in my plugin. I've tried with both QGIS 2.0 and 2.4.               layer.setEditForm('C:/my/full/pathname/with spaces/form.ui')     layer.setEditForm('C:/my/full/pathname/withoutspaces/form.ui')     layer.setEditForm("C:/same/as/above")     layer.setEditForm("C:\same\as\above")     layer.setEditForm('C:\same\as\above')


title: QGIS Reshape Tool

text: I just downloaded Quantum GIS and have been playing around with the editing features. For some reason, the Reshape Tool is not working for me. The manual said that at the end, you right-click and it will complete the edit. But I keep getting an Exception error that reads "Unknown exception." Any ideas?


title: Cut polygon shapefile by line shapefile

text: I have two shapefiles, a polygon shapefile (lake) and a line shapefile (boudary) I want to create a new shapefile "polygon to the north of the boundary" ![I want to create a new polygon "area north of the border line"](http://i.stack.imgur.com/zHjqU.png) What tool can I use in QGIS or other free GIS for this task?


title: ArcGIS 10.1 arcpy Update Cursor is returning a list instead of a row

text: I keep getting the following error: > Traceback (most recent call last): File > "R:\GIS_ScriptsandTools\PullHotParametersLabel.py", line 16, in TitleLine = > ""+row.getValue(IDfield)+" ("+ row.getValue(Datefield)+") " > > AttributeError: 'list' object has no attribute 'getValue' When I am trying to run a typical for row in cursor loop when working with an UpdateCursor. For the life of me, I can't figure out why it thinks the row is a list and not a row object. Here is my script:               import arcpy     SourceTable = arcpy.GetParameterAsText(0)     IDfield = arcpy.GetParameterAsText(1)     Datefield =arcpy.GetParameterAsText(2)     SelectedParams = arcpy.GetParameterAsText(3)     OutputField = arcpy.GetParameterAsText(4)     SelectedFields = SelectedParams.split(';')     AllFields = arcpy.ListFields(SourceTable)          FieldDict = {}     for field in AllFields:         FieldDict[field.name]=field.aliasName          with arcpy.da.UpdateCursor(SourceTable,"*") as cursor:         for row in cursor:                  TitleLine = "<und><bol>"+row.getValue(IDfield)+" </bol>("+ row.getValue(Datefield)+")</und>"             LabelList = [TitleLine]             for Field in SelectedFields:                 if row.getvalue(Field) <>0:                     Alias = FieldDict.get(Field)                     ParamValue = row.getValue(Field)                     ParamLine = '{0} = {1:,}'.format(Alias,ParamValue)                     LabelList.append(ParamLine)             if len(LabelList)<2:                 LabelList.append("All < 2L")             Label = '\n'.join(LabelList)             row[OutputField] = Label             cursor.updateRow(row)         del row, cursor      Does anybody have any ideas where I'm screwing up on this?


title: Google Earth Canada NTS 1:50000 Overlay

text: Not that long ago there was a Google Earth kml overlay file (ie. wms-ctis.kml) titled _Canada Geographic Names and Maps_ which included the Overlay Canadian topographic maps. This would allow one to query any geographic name in Canada and/or view either the 1:50K or 1:250K NTS maps in Google Earth. This was incredibly useful and it is very sad to see it disappear and no longer work. Does anyone know of an alternative source or similar file that could provide this functionality?


title: PostGIS Geometry folder is missing

text: I have installed PostgreSQL 9.2 and PostGIS 2.0.3 from Stackbuilder. My template_postgis_20 does not have a geometry table only spatial_ref_sys table. Is that a problem? and if not how can i create table using this postgis template with Geometry. I have both tried Windows 7 and 8. Thank you


title: Qgis: Can prt or python generate high res, scale 1:250 images of entire map?

text: I need to create a very high resolution, scale 1:250, image from my multi- layer map (30 seconds x 40 seconds). The PRT function in Qgis falls short as it reduces the resolution when using PRT. Also the canvas in Qgis only displays a limited portion of my map at the scale I want. My shape files have individual data polygons about 0.4 meter in size. These polygons get consolidated into larger, solid polygons by the PRT function if I change the scale of the map so it fits the canvas (about 1:10000). I need to generate very high resolution images, capable of enlargement by a photo viewer that retain the original scale (1:250) of the map. Size of the image file is not a problem nor is long processing time to create it. Hopefully there is a method to generate an image file directly using python that will solve this issue, but I don't know the command sequence. I can create an example map using fake data if this will help.


title: setStyle() function for GeoJSON features - Leaflet

text: Ok, since I already asked a very long question about this, but since it didn't get any new replies for a while, and not to get confused in details, I will keep this one simple the best way I can. If I'm not mistaken, a `setStyle` function for a named, particular feature would be as follows:               var bounds = [[54.559322, -5.767822], [56.1210604, -3.021240]];     var rect = L.rectangle(bounds, {color: "#ff7800", weight: 1}).addTo(map);          rect.setStyle({color: "#4B1BDE"});      ...which would change the color from orange to blue. I'm also aware of the `resetStyle()` function that will revert the style to the original. This is how I style my GeoJSON:               var everything = L.geoJson(myfile, {         onEachFeature: function(feature){             array_of_layers.addLayer(feature);         },         style: function(feature){                 switch(feature.properties.name){                 case "belgium": return belgium_style; break;                 case "bosnia": return bosnia_style; break;                 case "denmark": return denmark_style; break;                 case "great_britain": return britain_style; break;                 case "greece": return greece_style; break;                 case "italy": return italy_style; break;                 case "serbia": return serbia_style; break;                 case "spain": return spain_style; break;                 }         }          });      **What I want to do is to make just one country blue and the others gray, later in the code.** It's a two-step thing, to paint all countries to gray, and then make one blue. The first thing is, I need such a loop that would iterate over each feature and `setStyle()` for all countries to gray. Does it work if I just `everything.setStyle({color: "#4B1BDE"})` or something? The second thing is, (that's giving me sleepless nights) how do I select **just one feature** out of a group of GeoJSON polygons to work with? Just the country that I need to paint to blue. If it was a matter of mouse hovering, I could place an event listener as are done in Leaflet tutorials. But regardless of user interaction, I want to set and reset the style by calling it with its name, as I did with the rectangle above. Thanks in advance.


title: How to get Last ObjectID of feature class using ArcPy?

text: I have a feature class and I need to know the ObjectID of its last record. What is the best way to find out last ObjectID? Opening a search cursor and reading all the records is one way, but is there any other way?


title: create geometry in arcpad

text: i want to create a point geometry in ArcPad 10 VBScript, to make custom feature drawing. so i want create point geometry then                record.shape = point      or               records.AddNew(point)      i want to create point ,line and polygon geometries. i searched a lot but i haven't found how to create geometry. i don't know if there's another simple manner to do that. thanks in advance.


title: Two machines running the same arcpy based software - only one is leaking memory

text: Ok, so I'm running a rather large Python application based on arcpy on two separate machines - one of them is leaking memory and one is not. Both are running the same operation in the software. The leaking machines is:   * A nearly new Lenovo T420s   * Intel Core i5-2540M processor   * 8 GB DDR 3 RAM   * Windows 7 x64 Professional (non-domain joined)   * 128 GB Toshiba SSD   * running ArcGIS 10.0 SP4 with pyodbc installed into the Python directory(only third party module installed) The nonleaking machine is:   * A 1.5 year old Dell Optiplex 980   * Intel Core i7 870   * 8 GB DDR 3 RAM   * Windows 7 x64 Professional (domain joined)   * 500 GB Western Digital 7200RPM HDD   * running a brand new (today due to another issue with ArcGIS) intallation of ArcGIS 10.0 SP4 with pyodbc installed (also the only 3rd party module I also have this same software running on other machines without leaking, but I haven't run this same operation recently or on SP4, so I'm leaving them out. The software I'm running is here: https://bitbucket.org/nickrsan/pisces - The segment of it that's running is the mapping module (PISCES.mapping). Basically, it runs off the command line and is taking is input data from a Personal Geodatabase, generating layers based upon returned keys from stored queries, and adding those layers into template mxds and outputting them as PDFs/MXDs. All of the leaking portion uses arcpy.mapping and I suspect that's the culprit. The nonleaking machine runs the software utilizing approximate 200 MB of RAM the whole time. The leaking machine starts around 150 MB, and slowly winds its way up to around 1.25 GB where arcpy calls (load layer, output PDF, etc) start failing. Has anyone had any similar experiences? I know arcpy leaks memory sometimes and have had similar (but less severe) situations before. Any suggestions on how to track down a leak like this? Any best practices on the types of things I can do/avoid in arcpy to not have memory leaks? I've been using the gc module, but I think I'm not reading its output correctly. I'm relatively certain that it's not my code, but who knows. Thanks in advance for your help.


title: Rasterband Casting Probem

text: We have a program that was previously using ArcGIS 9.3 and we recently updated it to ArcGIS 10.0. The program was working before with a line of code that cast an object as a RasterBand. Now we are getting a java error that says IRasterBandProxy cannot be cast to RasterBand. Any ideas?


title: How to make a grid with a certain spacing between lines along X and Y axis?

text: We are doing the research about applying remote sensed data on hydro logical model . After delineating, the river basin from DEM data using Argis9.3, We have to map the remote sensed data derived from products like PERSIANN. Our task now is to create the grid with the same resolution with the available remote sensed data( for example 0.25degree x 0.25 degree). and have a certain corner that matches with the obtained data. Does anyone know how to deal with our problem? Your help is highly appreciated.


title: What are the main fields of applications of 3D GIS?

text: I'm doing a scholary research on 3D GIS but i have some problems in finding where they are used. Does anyone know of any real life examples where 3D GIS has been adopted?


title: How can I calculate the orientation of each side of a polygon in ArcGIS 10?

text: I want to examine the orientation of each line in a polygon so that I can calculate their solar exposure. Each polygon represents a building, and has an associated height. For the moment, I just want to consider the orientation, and will later consider shading issues. One approach I thought was to split the polygon into lines, and to calculate the orientation of each line, but the difficulty is that I then have to identify the outer face of that line. Although most of the polygons are simple four-side figures with straight lines, there are a small number where this is not the case (an issue which I just want to consider, but do not have to solve yet). I'm familiar with python and planned on doing this all from a script.


title: How can I change the line width of the outline for polygons in the Mapfile?

text: I display Polygons in the UMN Mapserver. With the `OUTLINE`-keyword I set a color for the outline of the polygons. How can I change the line width of this outline?


title: Has anyone successfully implemented a web service that sits on top of ArcSDE (SQL Server 2008)?

text: Here's the situation. We have planned point locations that lie within predefined polygonal project areas. Our IT department currently has these polygonal project areas in SQL Server 2008 as geography objects. When a user enters a point location in a interface, a query gets issued to the project areas table to return the project area that the new planned point location lies within. Works like a champ. However, these project area boundaries change of course in our GIS and then IT doesn't get the changes and we get data busts. IT has talked of writing a web service that would sit on top of our ArcSDE database (again, SQL Server 2008) so they could just make calls to the fresh data we always have in SDE. This of course would require us to store those datasets in SQL Server geom/geog datatypes. My question is this: Has anyone successfully implemented anything like this? If so, could you share your thoughts. If you have tried and failed, I'd like to hear about that as well. Thanks.


title: What is the FME Equivalent to the ArcGIS Union tool?

text: I'm attempting to union 2 large featureclasses in ArcCatalog, QGIS, and FME. What is the transformer that would be equivalent to the ArcGIS Union tool?


title: Grib2 Libraries For .Net

text: Are there any good libraries for accessing GRIB2 data in .Net? Specifically, I want to be able to quickly pull subsets of weather data from a grib2 file based on time and location.


title: Spatial Database Design (PostGIS) - Point data that doesn't always have spatial data

text: I am creating a spatial database using PostGIS. I am storing the attributes and location of some species records that can be represented as points. The records are recorded from either: a) static survey position (where previously I have just represented the species point as the static position) or b) a transect walked around the site and a position is marked when a record is recorded. The database stores the static positions in a a separate table, so for a) above I don't actually have to store a point against each of these records as I can do a join on the species table and the static location table. If I did store the location, if the static location changed, I would have to also make sure I moved the location of all of the species records associated with that point. The transect results are stored with spatial data as they could be recorded anywhere around the walked route and there is no other way to know where they were recorded. I would all of this data in the same table as it seems silly to have to maintain two tables for records about the same species (would have to make changes to the table structures in two tables). I can see 3 ways to deal with this: 1) Store the species data in two separate tables, one spatial and one not. 2) I am guessing you cant have NULL spatial values in a geometry column? 3) Give all the records spatial data, and move the points from the species table when the static table locations change? Can you have a relationship / procedure in the database that automatically does this for you? Which of these approaches is the most sensible? Or have I missed something. Any thoughts would be much appreciated.


title: How to create simple select script in ArcPy?

text: I am new to Python and want to create a script inside the ArcMap Python window. My current workspace is a series of rivers portrayed in the state of Ohio. I have a list of fish and the rivers (from my series of rivers) they are found in. Ideally, I would like the user to be able to enter the name of a fish and have the script highlight or select the rivers in the map. For example, if I were a fisherman who was not from this area and wanted to catch a certain fish (such as a smallmouth bass) I could type the fish into the script as a raw input and have the map highlight the rivers where this fish could be caught. I have an excel sheet with column titles as fish (i.e. largemouth bass, smallmouth bass, etc) and each row is a river name that they can be found in How might I go about doing this? Excel format:               Largemouth Bass (Col 1. Title)           Smallmouth Bass (Col 2. Title)     Ashtabula River                          Ashtabula River     Chagrin River                            Chagrin River     Cuyahoga River                           Conneaut Creek     Grand River                              Grand River                Huron River                              Huron River     Rocky River                              Rocky River     Sandusky River                           Vermillion River     Vermillion River


title: Select by Attributes Where Clause how to get current date?

text: I am trying to build a where clause for a select by attributes operation in ArcGIS. I have a feature class that has a date field that has dates in it. I wish to construct a where clause that selects records from 65 years ago or more. I want this to be generic, where the clause computes the 65 years from whatever today's date is without actually entering today's date. In other words, I do NOT want to enter something for a where clause like this: "date_field" >= date '1945-09-22' I want to have a where clause like this: "date_field" >= curdate() - 65 years The problem is that I can't figure out how to write this where clause. There is nothing in the ArcGIS documentation on composing where clauses with dates that includes obtaining the current date programmatically as opposed to entering a specific date. Help? Any input is appreciated.


title: QGIS Distance Matrix on large tables

text: I have a table with around 1 million addresses.I am trying to run a 'Nearest Address' on all of them and I understand QGIS has this capability. Yet I am wondering about the most efficient way to run this so that I don't have to have 1 million records searching on 1 million other records? I thought perhaps I could even sort the table somehow on Lat/Long and then only run each record on say the next 1000 records? Any thoughts from anyone who has run this before would be much appreciated.


title: Trying to use SQL API from Oracle PL/SQL to CartoDB

text: I am trying to make an http call from within an Oracle server to CartoDB. To do so I have developed a custom PL/SQL package that will generate the SQL and make the calls. A sample API call is below. The return status looks like a permissions problem ({"error":["permission denied for relation address"]} Any ideas?               http://mycompany.cartodb.com/api/v2/sql?q=INSERT     INTO address(address_id,street_number,city,state,zip_code)  SELECT     'BBE20140321095101832AM5656' address_id, '3469 Shelby Ganges Road'     street_number, 'Shelby' city, 'Ohio' state, '44875' zip_code WHERE NOT EXISTS (     select 1 from ADDRESS WHERE ADDRESS_ID =     'BBE20140321095101832AM5656')&api_key=488e214c52f78b55f2fe9c05daef94074c7e3cd2


title: How can I load VBA customizations stored in a "document-level custom template" into a Map Document created from a "Normal.mxt" template in ArcGIS 10?

text: I have all kinds of custom tools that are called from UI Controls stored in a document-level map document (my custom template). In earlier versions of ArcGIS these were stored in the global-level "Normal.mxt". Since this it is no longer possible to save UI controls within the Normal template in version 10 I use the above mentioned map document that has all of my VBA customization to create most of my maps within so that they may "inherit" the customization. This is a very crude solution since it created a problem with maintaining and upgrading my custom tools until i migrate to .NET. Is there a way to load document-level customization into a map document created from a Normal.mxt template?


title: Is anyone successfully using arcpy.ExportMetadata_conversion?

text: I don't know if it is a bug or my lack of python skills, but I keep getting 'ERROR 000816 The tool is not valid' in my python script when I call ExportMetadata_conversion Initially, I was attempting this on an sde database, but then I tried to take a step back just match the example in the documentation and use a file geodb:               import arcpy     from arcpy import env     env.workspace = "C:/data"     #set local variables     dir = arcpy.GetInstallInfo("desktop")["InstallDir"]     translator = dir + "Metadata/Translator/ESRI_ISO2ISO19139.xml"     arcpy.ExportMetadata_conversion ("data.gdb/roads", translator,          "roads_19139.xml")      But I received the same error. Is anyone successfully calling this method? If so, does your code match the documentation?


title: Is there a tool in arctoolbox for IDistanceOp.CostPathasPolyline?

text: Almost all the methods in IDistanceOp have a corresponding tool in the Distance Toolset Spatial Analyst Toolbox. However, I can't find anything in the toolbox for IDistanceOp.CostPathasPolyline. Does anyone know if such a tool exists?


title: Unzip a file using python caller in FME 2011

text: The idea is to download a zipped file from a website and unzip it using FME 2011. I have attached a screenshot of FME workbench which is not working. Any help is greatly appreciated. Thanks. Sam ![FME Workbench](http://i.stack.imgur.com/AL86f.jpg)


title: File paths in a PostGIS Table in QGIS

text: I have a PostGIS table which I have stored the location of some viewpoint photographs. I would like to store the location of the photo file in a field, but when I add a file path in the attribute table dialog in QGIS, QGIS/PostGIS removes '\' and spaces from the path. Can anyone help me work out why this is? I have done this before without any issues. If I enter the data through pgAdmin it comes across to QGIS correctly, and I can link to the photographs.


title: How to write a SQL WHERE Clause using Multiple Choice values in FME?

text: I published a parameter which is of Choice Multiple type. I wish to use values from this parameter for WHERE Clause when Reading Inputed Data in the begin of Workbench. So there is no problem if I have type of published parameter just Choice, then I use following: WHERE Clause Tile = '$(Name)' but I do not know how to write SQL WHERE Clause statement if I have Multiple Choice and want to read all data for Name of Tiles selected in Multiple Choice, because SQL form above does not work ... Example: In multiple choice i select the names: Anna, John, Tim. Then, I want that data for tiles that have those names are read.


title: Adding text as ICON in mapinfo in mapinfo buttons

text: It's a question with yes or no answer and if yes How? Is it possible to add text label as icon in mapinfo buttons.   * Instead of : > >     Create buttonpad "sitestatusbutton" As >       ToggleButton >         ID 1011 >       Icon MI_CURSOR_FINGER_LEFT <--------------------here >           Calling site_statusAction >           HelpMsg siteStatus(1) >     Show >     * Something like: > >     Create buttonpad "sitestatusbutton" As >       ToggleButton >         ID 1011 >       Icon "MYICON"<-------------here >           Calling site_statusAction >           HelpMsg siteStatus(1) >     Show</li> >   > >


title: GeoServer 2.3.0 Proxy Setting

text: I've upgraded to GeoServer 2.3.0 running on Apache Tomcat 7.0.39 on my host WebFaction. My websites and GeoServer are on same host but different ports. How do I go about setting a proxy so that my site domain can grab data from GeoServer. I remember using a proxy.cgi with openlayers, but I'm now using Leaflet. Is it possible to have my sites with domain capecodgis.com retrieve data from tomcat.capecodgis.com/geoserver. I apologize in advance if I'm not clear. This is very confusing. I've been reading different docs and forums to no avail


title: Route hatch labeling as stations (0+00) with vba

text: I have a route feature which I have set up to label every 500' with the distance along the line. I'd like to use typical stationing format (000+00). I've been using the following VBA script:               Function FindLabel ( esri__measure )     if esri_measure=0 then         FindLabel = "0+00"     elseif   esri_measure<1000 then          varleft = left(esri__measure, 1)          varright = right(esri__measure, 2)          FindLabel = varleft & "+" & varright     else          varleft = left(esri__measure, 2)          varright = right(esri__measure, 2)          FindLabel = varleft & "+" & varright     end if     End Function      each "else" works properly on it's own, but when all the script is entered, the first "else" to be satisfied is applied to all labels. For instance, with the script as copied above, every distance label reads exactly "0+00". To see what happens I try changing the script to:               Function FindLabel ( esri__measure )     if   esri_measure<1000 then          varleft = left(esri__measure, 1)          varright = right(esri__measure, 2)          FindLabel = varleft & "+" & varright     else          varleft = left(esri__measure, 2)          varright = right(esri__measure, 2)          FindLabel = varleft & "+" & varright     end if     End Function      With this script every label is formatted "X+00". When the value is greater than or equal to 1000, I need it to read "X0+00". The problem seems to be that the program tests only the first label (distance 0) against the "if" statement, then applies the outcome to every labeled location. How do I get each label to run through the "if" statement independently rather than only testing the first point?


title: SAGA kernel density error with Sextante Modeler

text: I can succesfully perform **SAGA kernel density** analysis with **Sextante (1.0.7)** in **QGIS 1.8.0.** But i want to built an automated tool for this job with **Sextante Modeler**. But i got foolish errors. This analysis tools works succesfully directly under Sextante but never works with Modeler. First screen shows kernel density analysis screen that called directly Sextante toolbox...![GUI of kernel density analysis](http://i.stack.imgur.com/c1JIN.png) Here is the another screen capture. Second screen show Modeler inputs and the error. i entered every input and output but everytime i got same error... ![Sextante Modeler screen](http://i.stack.imgur.com/XJAMw.png) . . **Here is the Sextante's configuration file for SAGA Kernel Density Analysis tool. Maybe you find the solution...**               Kernel Density Estimation      grid_gridding ParameterVector|POINTS|Points|0|False ParameterTableField|POPULATION|Weight|POINTS ParameterNumber|RADIUS|Radius|None|None|10 ParameterSelection|KERNEL|Kernel|[0] quartic kernel;1 gaussian kernel ParameterSelection|TARGET|Target Grid|[0] user defined Extent USER_XMIN USER_XMAX USER_YMIN USER_YMAX ParameterNumber|USER_SIZE|Cellsize|None|None|1.0 OutputRaster|USER_GRID|Grid **NOTES: ###Maybe you will say that QGIS Heatmap tool can handle this. You are correct but i have to build an automated analysis model so i have to work with analysis tools adapted with Sextante Modeler or Sextante Script but i have no idea with scripts...** **###I got error with Sextante 1.0.9 but i work succesfully 1.0.7...I need your help, thanks in advance...**


title: How to assigning multiple values (for single field) to a polygon layer

text: I am using arcmap 10 to generate a dataset that represents apartment building stock within particular suburbs. I have point data for the apartment buildings, with fields for the number of dwellings and car spaces at each apartment building. And then I have polygon data that maps the extent of suburbs. I'm using a spatial join to encode the apartment building attributes into the suburb polygons - but it seems this can only go so far. I can only get ONE value to transfer across for each field. The Merge Rule gives me some control of what that value is, but the problem is i need all of them. That is, I'd like to end up with multiple values recorded for single fields. I need this to do selections based on attributes. Say I want to select suburbs that have at least 1 apartment building with more than 100 dwellings, the number of dwellings for each and every apartment building needs to be encoded into the field (separated by semicolons say) so that they can all be checked against the selection criteria. Any suggestions how to do this? Is there a nesting function? Thanks in advance!


title: How to save python session in QGIS?

text: After running a lot of scripts in QGIS' python console, I have more than 10 dictionaries and lists, and I need these for further analysis. Is there a way to save python session, with all these variables? In general, my question is how can I export my variables from python console to something, so that they can be loaded in again at a future time? Thanks


title: Convert polygons to a number of randomly located points based on polygon size

text: ArcMap 10.1 on a PC. Working with a large shapefile with ~50,000 poygons with a wide range of size from 1,000+m. to <1m. I would like to convert these polygons into points so I can do some analysis in Maxent. Rather than have just one point for every polygon, I would like a large polygon to have more points than a small one. For example maybe a 100m polygon has 10 points in it, a 50m polygon has 5, and everything smaller has just 1. Ideally these points would just be randomly placed within the boundaries of the polygon it came from. Any ideas? Thanks,


title: Point to Polyline

text: I have shot points that I am trying to convert from Points into a polyline for linear distance for each Tract in ArcView 10.1. What is the best way to do that?   * Aaron, shot points are Seismic shot points (receivers and source points).


title: QGIS not rendering WMS layers

text: I'm a bit of newb' at GIS so be gentle with me. In Qgis I've set up a wms for the British Geological Survey: http://maps.bgs.ac.uk/ArcGIS/services/BGS_Detailed_Geology/MapServer/WMSServer? Now the layers are rendering at 1:25000 but at 1:50000 and above the screen goes blank. According to the BGS website they should be visible between 1:75000 and 1:10000. Does anyone have any suggestions as to what I'm doing wrong? Thanks in anticipation, Neil


title: Shows Netherlands region on map on every zoom level not other region in Open Street Maps

text: I want to show/display only Netherland region on Openstreet map, dont want to show other region on maps. Same thing apply for zoom in or zoom out also. I have tried with applying bounds but does not work,when user decrease the zoom level it shows the whole world map. Can you please help me or guide me in that. Waiting for your response.


title: How do I use ogr2ogr with python GDAL?

text: I've just installed GDAL, but I don't know how to import/use `ogr2ogr`. I've tried:               from osgeo import gdal     from osgeo import ogr      and searched both for `ogr2ogr`, but it wasn't there.


title: How do I use ogr2ogr to convert a GML to shapefile in Python

text: I am trying to convert a GML to an ESRI shapefile using ogr2ogr utility in a python script. I have successfully installed the GDAL/OGR package via osgeo but I am now struggling to find/understand any details on the syntax for using ogr2ogr in python. After importing ogr, all I have found is `ogr2ogr -f "ESRI Shapefile" output.shp input.gml` How do I use this in python? Is it as simple as assigning the input and output files? Everything I try I just get syntax errors. Any pointers in the right direction to get me started would be great. Many thanks


title: How can I get my Qgis area calculations in metres rather than degrees? My shapefiles are projected using WGS 1984 UTM fuseau 35 Sud

text: I am using Qgis 1.7.4 to calculate the area of polygons. My shapefiles have been projected using the coordinate system WGS 1984 UTM fuseau 35 Sud, but my area calculations are still displayed in degrees. how do I get them in metres?


title: Counting how many times a point goes through a loop and write the count to a text

text: I have a several feature classes that go through a loop to calculate points. I want to add up how many go through it each time so I can archive that number. I'm using ArcGIS 10.1. Lets say I have 100 shapefiles with anywhere between 1 - 50 points per shapefile. I want to run a count on the total number of points in those 100 shapefiles. Right now it's printing out the total number of points in each shapefile loop. Is there a way to add the 'count' up to create a total?               SMON = arcpy.ListFeatureClasses("SectMon_PLSS_*")     for fc in SMON:             filename =  "'" + arcpy.env.workspace + os.sep + fc + "'"             arcpy.env.overwriteOutput = True             arcpy.AddField_management(fc, 'shpname','text')             arcpy.AddField_management(fc, 'GIS_ID_PAM','text')             arcpy.AddField_management(fc, 'LATDMS','text')             arcpy.AddField_management(fc, 'LONGDMS','text')             arcpy.AddField_management(fc, 'DESC_','text')             arcpy.AddField_management(fc, 'StatePlane','text')             arcpy.AddField_management(fc, 'State','text')             for fc in SMON:                 result = arcpy.GetCount_management(fc)                 count = int(result.getOutput(0))                 x.append(count)                 c = sum(x for x in x)                 print "The cumulative sum for this iteration is %s" %c      This is what the printouts are giving me: > > > >>>      P:\Scripts\SystemData\TestData\TO PAM\140306_BCR_PNTS102034-102035&103003-103005 (1).shp >>>      This is located in Colorado North >>>      This is located in Colorado North >>>      This is located in Colorado North >>>      This is located in Colorado North >>>      This is located in Colorado North >>>      Fields Calculated >>>      The cumulative sum for this iteration is 5 >>>      #############COLORADO NORTH FINISHED################ >>>      P:\Scripts\SystemData\TestData\5pm\140226_BCR_PNTS100463-100465.shp >>>      This is located in North Dakota North >>>      This is located in North Dakota North >>>      This is located in North Dakota North >>>      Fields Calculated >>>      The cumulative sum for this iteration is 3 >>>      P:\Scripts\SystemData\TestData\5pm\140311_BCR_PNT100466.shp >>>      This is located in North Dakota North >>>      Fields Calculated >>>      The cumulative sum for this iteration is 1 >>>      P:\Scripts\SystemData\TestData\5pm\140312_BCR_PNTS100467-100469.shp >>>      This is located in North Dakota North >>>      This is located in North Dakota North >>>      This is located in North Dakota North >>>      Fields Calculated >>>      The cumulative sum for this iteration is 3 >>>      #############NORTH DAKOTA FINISHED################ >>>   I want to be able write to a log file each time it goes through a loop. I'm trying this, but it's giving me an error.                           txtFile = open(count,"a")                 txtFile.write("This is the total amount " + str(c))      It's telling me "coercing to Unicode: need string or buffer, int found"


title: Multivalue processing in the context of the modelbuilder

text: I am new to python so please bear with me. I have the following script that is supposed to get a multivalue parameter of several shapefiles as an input and then run the unionanalysis on all these shapefules to create a single feature class or shapefile output. The script seems to run fine however only processes up to step 2 of 16 and does not create any output file. I do have the output file specified however and all the parameters are in the right order with the input accepting multivalue variables and both having feature class as their data type. I need to precise that this is within the ModelBuilder framework... Here is the script:                   import arcpy              myList = arcpy.GetParameterAsText(0)              unionResult = arcpy.GetParameterAsText(1)         arcpy.Union_analysis(myList, unionResult, "ALL", "", True)      To further clarify: Parameter (0) is an input, feature class, multivalue (yes) Parameter (1) is an output, feature class, multivalue (No)


title: Problem Docking Add-In Toolbar in ArcGIS 10 SP4/5

text: I have developed an Add-In Toolbar for ArcGIS 10. It docks just fine on my development computer (which is running 10.0, no service packs according to ESRI's patch finder). It has been deployed to users who are running 10.0, SP4 and SP5, and no matter what they (or I) do, I cannot get the toolbar to stay docked where they want it to dock. I do not set the docking behavior in the code so I can't see any way I might be inadvertently affecting it programmatically. This issue seems identical to the one here: http://support.esri.com/en/bugs/nimbus/TklNMDYwNDI2, which was supposedly fixed in SP3. I wonder if it broke again in SP4??? Anyone else notice this? My main tester seems to remember that when she first tested it (on SP3) it docked just fine, but at some point after that it stopped docking - looking at the timeline of her testing, there is a good chance the problem started manifesting after she upgraded to SP4. I have tried every trick I can find online (all to no avail):   * "Save all customizations to the document" is unchecked   * "Create new toolbars and menus in the document" is unchecked   * I have deleted normal.mxt and had ArcGIS recreate it   * I have run AdvancedArcMapSettings.exe and changed monitor dimensions to match the actual dimensions of the monitor  Any thoughts or other suggestions? Thanks!!


title: Error executing CQL on layer other then demo layers

text: I have some sort of wierd problem as I was playing with CQL. CQL works well with given layer but when i try with layer loaded from mapinfo or some views then it does not work. As in the picture when I click the feature it shows the column called **NAME** but when I query using it (or other column) generates error like:               Property 'NAME' could not be found in DISTRICT_BND      Why is the query not recognizing my column? * * * ![enter image description here](http://i.stack.imgur.com/Yfm99.png)


title: How to open .shx files?

text: I want to see the contents of an .shx file to see how the index information is stored in these files so that I can get some help in writing my own shape files. Can any one guide me on how can I open and use the .shx files?


title: How to use polyline Z in shapefiles?

text: I am researching about shapefiles. I am not able to understand the usage of polyline Z parameters. I have searched the internet but I didn't come across any example where this parameter is being used. Please help me to understand this parameter better?


title: Updating polygon geometry: deleting inner rings with python

text: I've searched, but can't find clear explanations how can I simply delete inner rings from existing polygons with python or extract outter rings and overwrite the geometry. Seems GDAL/OGR should be the simplest, but still I need some help... Simple vision of mine would be:               dataSrc=driver.Open(shp)     layer=dataSrc.GetLayer()     for i in xrange(layer.GetFeatureCount()):         feature = layer.GetFeature(i)         geometry = feature.GetGeometryRef()         if geometry.GetArea()<200:             layer.DeleteFeature(i)         if geometry.GetGeometryCount()>1:             geometry=geometry.GetGeometryRef(0) #overwrite geom, getting the first ring outer             feature.SetGeometry(geometry)             layer.SetFeature(feature)      geometry.Destroy() This doesn't loop through, it's just how I would like to see. As I know the first ring or the first part of the geometry should be the outer ring, so I assumed I can just ask for it like that.Suppose that's not a problem. So, probably the question is just about Updating Geometry. I open dataset withouth specific- Reading or Writing, so I assume it should let me write in it not just read. But if I do not want to have a new dataset, can I just update write like that or i should anywhays better create a new dataset? Or maybe there is even more simple way to solve this with python?


title: Geoserver build error

text: i am getting a build error with mvn clean install on geoserver 2.2.x branch saying               Running org.geoserver.platform.exception.ExceptionLocalizeTest     Tests run: 6, Failures: 0, Errors: 0, Skipped: 0, Time elapsed: 0.02 sec     Running org.geoserver.platform.exception.GeoServerExceptionsTest     Tests run: 5, Failures: 0, Errors: 0, Skipped: 0, Time elapsed: 0 sec     Running org.geoserver.platform.GeoServerExtensionsTest     Oct 03, 2012 11:30:39 AM org.geoserver.platform.GeoServerExtensions checkContext     SEVERE: Extension lookup occured, but ApplicationContext is unset.     Oct 03, 2012 11:30:39 AM org.geoserver.platform.GeoServerExtensions extensions     SEVERE: bean lookup error     java.lang.IllegalArgumentException: Argument "value" should not be null.         at org.geotools.util.SoftValueHashMap.ensureNotNull(SoftValueHashMap.java:140)         at org.geotools.util.SoftValueHashMap.put(SoftValueHashMap.java:281)         at org.geoserver.platform.GeoServerExtensions.extensions(GeoServerExtensions.java:113)         at org.geoserver.platform.GeoServerExtensions.extensions(GeoServerExtensions.java:150)         at org.geoserver.platform.GeoServerExtensions.extensions(GeoServerExtensions.java:239)         at org.geoserver.platform.GeoServerExtensionsTest.testExtensionFilterByName(GeoServerExtensionsTest.java:140)         at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method)         at sun.reflect.NativeMethodAccessorImpl.invoke(Unknown Source)         at sun.reflect.DelegatingMethodAccessorImpl.invoke(Unknown Source)         at java.lang.reflect.Method.invoke(Unknown Source)         at junit.framework.TestCase.runTest(TestCase.java:164)         at junit.framework.TestCase.runBare(TestCase.java:130)         at junit.framework.TestResult$1.protect(TestResult.java:106)         at junit.framework.TestResult.runProtected(TestResult.java:124)         at junit.framework.TestResult.run(TestResult.java:109)         at junit.framework.TestCase.run(TestCase.java:120)         at junit.framework.TestSuite.runTest(TestSuite.java:230)         at junit.framework.TestSuite.run(TestSuite.java:225)         at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method)         at sun.reflect.NativeMethodAccessorImpl.invoke(Unknown Source)         at sun.reflect.DelegatingMethodAccessorImpl.invoke(Unknown Source)         at java.lang.reflect.Method.invoke(Unknown Source)         at org.apache.maven.surefire.junit.JUnitTestSet.execute(JUnitTestSet.java:96)         at org.apache.maven.surefire.junit.JUnit3Provider.executeTestSet(JUnit3Provider.java:117)         at org.apache.maven.surefire.junit.JUnit3Provider.invoke(JUnit3Provider.java:94)         at sun.reflect.NativeMethodAccessorImpl.invoke0(Native Method)         at sun.reflect.NativeMethodAccessorImpl.invoke(Unknown Source)         at sun.reflect.DelegatingMethodAccessorImpl.invoke(Unknown Source)         at java.lang.reflect.Method.invoke(Unknown Source)         at org.apache.maven.surefire.util.ReflectionUtils.invokeMethodWithArray(ReflectionUtils.java:164)         at org.apache.maven.surefire.booter.ProviderFactory$ProviderProxy.invoke(ProviderFactory.java:110)         at org.apache.maven.surefire.booter.SurefireStarter.invokeProvider(SurefireStarter.java:175)         at org.apache.maven.surefire.booter.SurefireStarter.runSuitesInProcessWhenForked(SurefireStarter.java:107)         at org.apache.maven.surefire.booter.ForkedBooter.main(ForkedBooter.java:68)      i tried building with both maven 2 and maven 3 on windows 7 with no chance. complete output file is here http://www.4shared.com/office/EvC5AB-w/out.html? edit: compilation fails with this message               Results :          Failed tests:   testDefaultStoredQuery(org.geoserver.wfs.v2_0.GetFeatureTest): expected:<[1]> but was:<[0]>          Tests run: 422, Failures: 1, Errors: 0, Skipped: 0          [INFO] ------------------------------------------------------------------------     [INFO] Reactor Summary:     [INFO]      [INFO] GeoServer ......................................... SUCCESS [0.831s]     [INFO] Core Platform Module .............................. SUCCESS [3.367s]     [INFO] Open Web Service Module ........................... SUCCESS [7.914s]     [INFO] Main Module ....................................... SUCCESS [4:11.637s]     [INFO] Web Feature Service Module ........................ FAILURE [5:37.378s]     [INFO] Web Coverage Service Module ....................... SKIPPED     [INFO] Web Map Service Module ............................ SKIPPED     [INFO] GeoServer Web Modules ............................. SKIPPED     [INFO] Core UI Module .................................... SKIPPED     [INFO] Security UI Module ................................ SKIPPED     [INFO] GeoServer Security Modules ........................ SKIPPED     [INFO] GeoServer JDBC Security Module .................... SKIPPED     [INFO] GeoServer LDAP Security Module .................... SKIPPED     [INFO] Web Coverage Service 1.0 Module ................... SKIPPED     [INFO] Web Coverage Service 1.1 Module ................... SKIPPED     [INFO] GeoWebCache (GWC) Module .......................... SKIPPED     [INFO] REST Support Module ............................... SKIPPED     [INFO] REST Configuration Service Module ................. SKIPPED     [INFO] WMS UI Module ..................................... SKIPPED     [INFO] GWC UI Module ..................................... SKIPPED     [INFO] WFS UI Module ..................................... SKIPPED     [INFO] Demoes Module ..................................... SKIPPED     [INFO] WCS UI Module ..................................... SKIPPED     [INFO] GeoServer Web Application ......................... SKIPPED     [INFO] Community Space ................................... SKIPPED     [INFO] GeoServer Extensions .............................. SKIPPED     [INFO] ------------------------------------------------------------------------     [INFO] BUILD FAILURE     [INFO] ------------------------------------------------------------------------     [INFO] Total time: 10:01.578s     [INFO] Finished at: Wed Oct 03 11:40:36 EEST 2012     [INFO] Final Memory: 25M/67M     [INFO] ------------------------------------------------------------------------     [ERROR] Failed to execute goal org.apache.maven.plugins:maven-surefire-plugin:2.10:test (default-test) on project wfs: There are test failures.     [ERROR]      [ERROR] Please refer to C:\Users\bb\Documents\Projects\geoserver\src\wfs\target\surefire-reports for the individual test results.     [ERROR] -> [Help 1]     [ERROR]      [ERROR] To see the full stack trace of the errors, re-run Maven with the -e switch.     [ERROR] Re-run Maven using the -X switch to enable full debug logging.     [ERROR]      [ERROR] For more information about the errors and possible solutions, please read the following articles:     [ERROR] [Help 1] http://cwiki.apache.org/confluence/display/MAVEN/MojoFailureException     [ERROR]      [ERROR] After correcting the problems, you can resume the build with the command     [ERROR]   mvn <goals> -rf :wfs