title: IDatasetContainer2 CreateDataset method fails after upgrading to ArcObjects 10

text: I have a Network Analyst routine that worked for 9.3 and does not work with 10. I was closely following the example at: http://help.arcgis.com/en/sdk/10.0/arcobjects_net/conceptualhelp/index.html#/How_to_create_a_network_dataset/0001000000w7000000/ It bombs when it gets to the line that actually creates the Network Dataset:               INetworkDataset networkDataset = (INetworkDataset)datasetContainer2.CreateDataset(deDataset);      Did something change that would cause these objects (used with Network Analyst) to not work the same way? This is the exact error message:               System.Runtime.InteropServices.COMException was unhandled Message=Exception from HRESULT:      0x8004228E Source=ESRI.ArcGIS.Geodatabase ErrorCode=-2147212658 StackTrace: at      ESRI.ArcGIS.Geodatabase.IDatasetContainer2.CreateDataset(IDEDataset DataElement) at      NetAnalystTestApp.MainForm.PutRoadLayerOnMap() in C:\@ALPACAS\AlpacasCT\experimental\network      routing\NetAnalystTestApp\NetAnalystTestApp\MainForm.cs:line 352 at      NetAnalystTestApp.MainForm.MainForm_Load(Object sender, EventArgs e) in      C:\@ALPACAS\AlpacasCT\experimental\network      routing\NetAnalystTestApp\NetAnalystTestApp\MainForm.cs:line 83 at      System.Windows.Forms.Form.OnLoad(EventArgs e) at      System.Windows.Forms.Form.SetVisibleCore(Boolean value) at      System.Windows.Forms.Control.set_Visible(Boolean value) at      System.Windows.Forms.Application.ThreadContext.RunMessageLoopInner(Int32 reason,      ApplicationContext context) at      System.Windows.Forms.Application.ThreadContext.RunMessageLoop(Int32 reason,      ApplicationContext context) at System.Windows.Forms.Application.Run(Form mainForm) at      NetAnalystTestApp.Program.Main() in C:\@ALPACAS\AlpacasCT\experimental\network      routing\NetAnalystTestApp\NetAnalystTestApp\Program.cs:line 39 at      System.AppDomain._nExecuteAssembly(RuntimeAssembly assembly, String[] args) at      System.AppDomain.ExecuteAssembly(String assemblyFile, Evidence assemblySecurity, String[]      args) at Microsoft.VisualStudio.HostingProcess.HostProc.RunUsersAssembly() at      System.Threading.ThreadHelper.ThreadStart_Context(Object state) at      System.Threading.ExecutionContext.Run(ExecutionContext executionContext, ContextCallback      callback, Object state, Boolean ignoreSyncCtx) at      System.Threading.ExecutionContext.Run(ExecutionContext executionContext, ContextCallback      callback, Object state) at System.Threading.ThreadHelper.ThreadStart() InnerException:      Thanks, Renee


title: Why would Openlayers always render over other divs?

text: I'm having trouble getting an Openlayers div to allow other divs to appear above it. Having tried the solution proposed here as well as the more classy-looking one here, I haven't had any luck yet. Openlayers always orders itself over other divs. I've also set the z-index for my other div (which I want on top) as high as it goes, still no effect. Is there somewhere it could be getting overridden? Or some setting that needs an extra trick to get rid of? I have asked this question on StackOverflow as well, where it's arguably more relevant. I'm only posting this here because someone here may have specific Openlayers experience on this matter.


title: Update geometry(polygon) of shape file based on a point layer

text: I have two shape files.   1. Polygon shape file   2. Point file contains some points. Now i want to delete all the vertex of polygon shape file which are available in point file. I am using ArcGIS(ArcEditor) platform. Any ideas.


title: Available iOS GIS applications/framworks for mobile data accquisition

text: I'm starting to develop an application for iPad which should enable the users to go to the field and acquire georeferenced information. The idea is to enable acquisition of points and polygons, and the filling of a set of attributes for the geometry. Also important is the ability to connect to externals geo-referenced databases, such as a PostGIS one. Ideally, I would want to be somewhat similar to the ArcGis for iOS application, which has a good user interface. Can you guys point me to other applications similar to what I am developing? I already know the following ones:   * Geo-jot   * pcMapper   * iGIS   * GISRoam   * field_assets   * Geopaparazzi Possibly you can point me to something with similar functionality to the ArcGis for iOS? OR Interesting programming frameworks to achieve the desired functionality on iOS? Thank you.


title: Error when pointing Geoserver to PostGIS

text: I have a PostGIS database that contains all my GIS data along with metadata. I have it published as a Store in GeoServer. They are both installed on the same windows server 2008r2. When i try to add layers from the store resource and publish it, i encounter an error at the bottom of the "Data" tab of the resource configuration page. Any ideas what might be my problem? ![enter image description here](http://i.stack.imgur.com/090IC.png) When i click on the yellow warning icon i get: ![enter image description here](http://i.stack.imgur.com/0sdBe.png) Here is the error:               org.apache.wicket.WicketRuntimeException: Exception in rendering component: [Component id = message] at org.apache.wicket.Component.renderComponent(Component.java:2563) at org.apache.wicket.markup.html.WebComponent.onRender(WebComponent.java:60) at org.apache.wicket.Component.render(Component.java:2361) at org.apache.wicket.MarkupContainer.renderNext(MarkupContainer.java:1298) at org.apache.wicket.MarkupContainer.renderAll(MarkupContainer.java:1428) at org.apache.wicket.Page.onRender(Page.java:1502) at org.apache.wicket.Component.render(Component.java:2361) at org.apache.wicket.Page.renderPage(Page.java:913) at org.apache.wicket.request.target.component.PageRequestTarget.respond(PageRequestTarget.java:63) at org.apache.wicket.request.AbstractRequestCycleProcessor.respond(AbstractRequestCycleProcessor.java:104) at org.apache.wicket.RequestCycle.processEventsAndRespond(RequestCycle.java:1245) at org.apache.wicket.RequestCycle.step(RequestCycle.java:1316) at org.apache.wicket.RequestCycle.steps(RequestCycle.java:1418) at org.apache.wicket.RequestCycle.request(RequestCycle.java:532) at org.apache.wicket.protocol.http.WicketFilter.doGet(WicketFilter.java:356) at org.apache.wicket.protocol.http.WicketServlet.doGet(WicketServlet.java:124) at javax.servlet.http.HttpServlet.service(HttpServlet.java:707) at javax.servlet.http.HttpServlet.service(HttpServlet.java:820) at org.springframework.web.servlet.mvc.ServletWrappingController.handleRequestInternal(ServletWrappingController.java:158) 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:875) at org.springframework.web.servlet.DispatcherServlet.doService(DispatcherServlet.java:809) at org.springframework.web.servlet.FrameworkServlet.processRequest(FrameworkServlet.java:571) at org.springframework.web.servlet.FrameworkServlet.doGet(FrameworkServlet.java:501) 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:502) at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1124) at org.vfny.geoserver.filters.SetCharacterEncodingFilter.doFilter(SetCharacterEncodingFilter.java:108) at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1115) at org.acegisecurity.util.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:265) at org.acegisecurity.intercept.web.FilterSecurityInterceptor.invoke(FilterSecurityInterceptor.java:107) at org.acegisecurity.intercept.web.FilterSecurityInterceptor.doFilter(FilterSecurityInterceptor.java:72) at org.acegisecurity.util.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:275) at org.acegisecurity.ui.ExceptionTranslationFilter.doFilter(ExceptionTranslationFilter.java:124) at org.acegisecurity.util.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:275) at org.acegisecurity.providers.anonymous.AnonymousProcessingFilter.doFilter(AnonymousProcessingFilter.java:125) at org.acegisecurity.util.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:275) at org.acegisecurity.ui.rememberme.RememberMeProcessingFilter.doFilter(RememberMeProcessingFilter.java:142) at org.acegisecurity.util.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:275) at org.acegisecurity.wrapper.SecurityContextHolderAwareRequestFilter.doFilter(SecurityContextHolderAwareRequestFilter.java:81) at org.acegisecurity.util.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:275) at org.acegisecurity.ui.AbstractProcessingFilter.doFilter(AbstractProcessingFilter.java:271) at org.acegisecurity.util.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:275) at org.acegisecurity.ui.logout.LogoutFilter.doFilter(LogoutFilter.java:110) at org.acegisecurity.util.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:275) at org.acegisecurity.context.HttpSessionContextIntegrationFilter.doFilter(HttpSessionContextIntegrationFilter.java:249) at org.acegisecurity.util.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:275) at org.acegisecurity.util.FilterChainProxy.doFilter(FilterChainProxy.java:149) at org.acegisecurity.util.FilterToBeanProxy.doFilter(FilterToBeanProxy.java:98) at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1115) at org.geoserver.filters.LoggingFilter.doFilter(LoggingFilter.java:73) at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1115) at org.geoserver.filters.ReverseProxyFilter.doFilter(ReverseProxyFilter.java:183) at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1115) at org.geoserver.filters.GZIPFilter.doFilter(GZIPFilter.java:41) at org.mortbay.jetty.servlet.ServletHandler$CachedChain.doFilter(ServletHandler.java:1115) at org.mortbay.jetty.servlet.ServletHandler.handle(ServletHandler.java:361) 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:766) at org.mortbay.jetty.webapp.WebAppContext.handle(WebAppContext.java:417) at org.mortbay.jetty.handler.ContextHandlerCollection.handle(ContextHandlerCollection.java:230) at org.mortbay.jetty.handler.HandlerCollection.handle(HandlerCollection.java:114) 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:534) at org.mortbay.jetty.HttpConnection$RequestHandler.headerComplete(HttpConnection.java:864) at org.mortbay.jetty.HttpParser.parseNext(HttpParser.java:533) at org.mortbay.jetty.HttpParser.parseAvailable(HttpParser.java:207) at org.mortbay.jetty.HttpConnection.handle(HttpConnection.java:403) at org.mortbay.io.nio.SelectChannelEndPoint.run(SelectChannelEndPoint.java:409) at org.mortbay.thread.QueuedThreadPool$PoolThread.run(QueuedThreadPool.java:522) Caused by: java.util.MissingResourceException: Unable to find resource: featureTypeReloadWarning for component: resource:tabs:panel:theList:1:content [class=org.geoserver.web.data.resource.FeatureResourceConfigurationPanel] at org.apache.wicket.Localizer.getString(Localizer.java:266) at org.apache.wicket.model.StringResourceModel.getString(StringResourceModel.java:491) at org.apache.wicket.model.StringResourceModel.load(StringResourceModel.java:633) at org.apache.wicket.model.LoadableDetachableModel.getObject(LoadableDetachableModel.java:114) at org.apache.wicket.model.StringResourceModel$AssignmentWrapper.getObject(StringResourceModel.java:230) at org.apache.wicket.Component.getModelObject(Component.java:1565) at org.apache.wicket.Component.getModelObjectAsString(Component.java:1587) at org.apache.wicket.markup.html.basic.Label.onComponentTagBody(Label.java:111) at org.apache.wicket.Component.renderComponent(Component.java:2524) ... 73 more


title: "Thresholding" for Stream Raster

text: Using ASTER DEM, I have created Flow Accumulation. However, to create stream order, I was asked to input flow accumultation and stream raster. But I cannot find how to get the stream raster. Some forums say that I should threshold flow accumulation. How do I do this? Please guide me in the right direction. I am using ArcGIS 10.


title: GIS Buzzwords Topic

text: What are the main buzzwords today in the GIS market? Let's place the meaning of these ones too, so beginners can enjoy the topic as well. Example: **SDI** "A spatial data infrastructure (SDI) is a framework of spatial data, metadata, users and tools that are interactively connected in order to use spatial data in an efficient and flexible way. Another definition is the technology, policies, standards, human resources, and related activities necessary to acquire, process, distribute, use, maintain, and preserve spatial data." wikipedia definition One answer per post :)


title: Python objects in my ArcMap addIn aren't instantiated

text: I've got a weird problem. My addin throws problems intermittently where it doesn't recognize buttons or tools when I try to toggle the visibility in my code. It seems as though the objects aren't instantiated as trying to toggle them works after I try to interact with them. What is the problem and how do I begin to debug it? My code is far too long to post but if there's something I can provide to help can you let me know? Here's a scenario: User starts Arc and clicks a button. The button tries to enable another button with:               anotherButton.enabled = False      Arc throws an error:               global variable "anotherButton" not defined.      User clicks on "anotherButton" and then tries the first button again. This time the command executes and "anotherButton" is disabled.


title: How to use SAGA ShapesToGrid via Processing in a Python plugin?

text: I run these lines in my code:               vectorlayer_vector =  ftools_utils.getVectorLayerByName(inputLayer)     Elevation = self.lineAttrib.currentText()     cellSize = int(self.linePix.value())     outPath = self.inputLayerCombo3.text()     output = ftools_utils.getRasterLayerByName(outPath)           Processing.runAlgorithm("saga:shapestogrid", vectorlayer_vector, Elevation, 0, 0, 4, None, cellSize, output)      When I run, I get an tiff image but with nothing. The .aux file is not created and when I try to open the tiff image that I created appears this message: c:...teste3.tif is not a supported raster data source So, I guess that the output is wrong. What I have to do to get a correct output? What I missing? * * * I update my code but I have no results yet. Please see if you can help me. My code are:               inputLayer = self.inputLayerCombo.currentText()     Elevation = str(self.lineAttrib.currentText())     cellSize = int(self.linePix.value())     outPath = self.inputLayerCombo3.text()          Processing.initialize()     Processing.runAlgorithm("saga:shapestogrid", inputLayer, Elevation, 0, 0, 4, "10000.0,14000.0,50000.0,54000.0" , cellSize,outPath)      Is a simple code but I don't have any results. I don't know what else to do. It doesn't create a raster.


title: ArcGIS Desktop can't return concurrent license

text: We have one concurrent use ArcView license, managed by a license server box in the office. The server reports that one of our desktops has the license 'borrowed', but when looking at the Desktop Admin tool on that box, it doesn't show anything about having a license checked out. It seems the license is somewhere lost in the ether. Is there a way, from the desktop side, to 'force return' or destroy the local license file? There appears to be no way for me to get the license returned from the desktop side. I've filed a report with Esri, but not holding out hope that my issue will be fixed ASAP. Anyone here have experience with this problem? Would a reinstall of the desktop software help?


title: Openlayers WMS getFeatureInfo get features with geometry attribute

text: I am displaying map using WMSlayer and implemented WMSGetFeatureInfo control on wms layer. Tap on a location and get the features using WMS getFeatureInfo service. This service is returning me all the features are located at given Point successfully with feature attributes. My Requirement Along with the Feature attributes, i need to have the feature geometry information too, which not geting as part of feature information from service. So is there any way to capture the geometry information from OpenLayers GetFeatureInfo control? Or Is that a limitation in WMS GetFeatureInfo service?


title: Setting background RGB value to invisible for multiple rasters in ArcGIS Desktop?

text: I'm using ArcGIS Desktop 10.0. I need to set my background rgb value of 255,255,255 to transparent. In the layer-Properties Windows > Symbology I can do it for one layer, but can I batch set this, or use arcpy? I have to do it for around 300 layers.


title: How can I create Multiple joins over 5 tables in ArcGIS?

text: I have a relational database which uses the following relational schema to create a query: Parcel.IsBooked is joined to Bookinglocation.ID Bookinglocation.IsPartOf is joined to Bookingpage.id Bookingpage.id is joined to Namenumber.isPartOf Namenumber.called is joined to Person.id This is the relationship path from the shapefile "parcel" to the table "Person" which contains the name of the parcel owner. My question is....How is it possible with ArcGIS to select parcels from the shapefile and see the selection in the "Person" table? I seem to only be able to join a shp to a table...but not a table to a table. I need to actually remove names out of the "person" table which fulfill a certain criteria...so exporting tables joined to shapefiles won´t do the job. Is this possible with ArcGIS?


title: How would I create an input mask for ArcGIS Desktop?

text: Someone asked me if we can create a simple input mask for geodatabase feature classes so that editing of attributes becomes easier. I am however not sure how this works with ArcGIS 10.x Is this possible with VBA or do I need 3rd Party tools to create masks? The last time I did something similar I could only create tools and toolbars...but no window controls...How could this be done now?


title: Bounding the Google Maps Search area to my Fusion Layer area?

text: I have a searchable Google Fusion Layer of Glasgow Postcodes that highlights the searched areas on a map. The search works fine if I enter a full postcode (for e.g. G1 1BL) that matches the data in the Fusion Table. However, if I search for just the first part of the postcode to display all the postcodes in the sector ('G1'), instead of highlighting all the areas on my fusion layer with G1, it relocates my map to China where there is a highway named 'G1'. This isn't the case for all sectors, for example if I search G84 then it highlights the corresponding areas on my map. This happens with other searches too where unless its a valid postcode I am relocated to the other side of the world to another place with this name as opposed to the fusion layer I am querying. Is there a way to fix this or will I have to limit the search to only full postcodes? Thanks in advance!


title: WFS-T demo with local Geoserver not responding

text:
i have installed the geoserver in my local system and i tried to run the example code from openlayers the code as below               <!DOCTYPE html>                   <meta http-equiv="Content-Type" content="text/html; charset=utf-8">         <meta name="viewport" content="width=device-width, initial-scale=1.0, maximum-scale=1.0, user-scalable=0">         <meta name="apple-mobile-web-app-capable" content="yes">         <title>OpenLayers WFS Protocol with Filter</title>         <link rel="stylesheet" href="http://localhost:8080/openlayers/theme/default/style.css" type="text/css">

                     When i click on the save the request to geoserver is in pending for more then 10 min ![enter image description here](http://i.stack.imgur.com/3CO5m.png)


title: Arcpy scripts fail on big datasets - anyone has experienced similar behavior?

text: Here a question that is not extremely specific but maybe someone has had similar experiences... I have a few Arcpy scripts that work perfectly when used with small datasets but that silently fail when processing big datasets. E.g. I am using Network Analyst tools to determine, all over the Swiss road network, the accessibility of various services (hospital, mall, bakery, school...). I basically create an OD matrix containing the distance between each populated place and the closest service. Apparently nothing rocket- science. When my script is runnung against a small subset, it performs well and spits out a result. When run against the full road network it (usually) doesn't explicitely fail but all the output tables are just empty. I experience a similar behavior when trying to perform a kind of custom filter on a point dataset: when using a few hundreds of points, no problem but when processing the full 4 Million points, the script runs without any warning but outputs an empty table. I am aware that the question is quite vague but it's on purpose! As I seem to experience similar problems in unrelated situations, I would assume that there is a more global issue to address before trying to work on those specific scripts... So, has anyone experienced similar behavior? FYI: ArcGIS 10.1, Python 2.7, Windows Server 2008, 24GB RAM, Intel Xeon 2.67Ghz (2 processors)


title: How to record a base layer change event in OpenLayers?

text: The `OpenLayers.Control` class possesses an `eventListener` property which could be used in the layer switcher control. I am not sure how to make use of it to detect a base layer change event so that I can store the variable name of the chosen layer to be used subsequently when the user navigate back to the map page. Can anyone give me some hints?


title: Find the nearest point on a line

text: I have a line that exists between two latlon positions, and a point on a certain latlon position. What I want to know is what the nearest point on that line is (in respect to the other point). So a latlon position on that line. I know how to do the basic math, but I just can't get my head around the latlong calculation. Can anyone help me or give me some pointers on how to do it? Thanks


title: How to create a SHP file in QGIS using the spreadsheet?

text: I have a spreadsheet which have columns (X, Y, Att1, Att2 ...). I want to create a point shp file, all column in the spreadsheet to be the column in shp file. How can I do this in QGIS? Thanks.


title: On using R data with ArcPy to make histogram in ArcMap

text: I am writing a tool in python which passes parameters to a R script, which in turn creates a DBF file. I read this DBF and want to make a histogram out of the data inside ArcMap as one of the results of the tool. The code snippet is as follows:               import arcpy          simloglkhd = "H:\\CIV\\KULLDORFF\\DATA\\PLOTDATA.dbf"     inputTemplate = "C:\\Program Files (x86)\\ArcGIS\\Desktop10.1\\GraphTemplates\\Default.SimplifiedChinese.tee"     output_graph_name = "VerticalBarGraph1"     verticalBarGraph_bmp = "H:\\CIV\\KULLDORFF\\DATA\\plotdata1.bmp"     # Create the graph     graph = arcpy.Graph()          # Add a vertical bar series to the graph     graph.addSeriesHistogram(simloglkhd, "simloglkhd",10)          # Specify the title of the left axis     graph.graphAxis[0] = "Frequency"          # Specify the title of the bottom axis     graph.graphAxis[2] = "log(lambda)"          # Specify the title of the Graph     graph.graphPropsGeneral.title = "Monte Carlo Distribution of Lambda"          # Output a graph, which is created in-memory     arcpy.MakeGraph_management(inputTemplate, graph, output_graph_name)          # Save the graph as an image     arcpy.SaveGraph_management(output_graph_name, verticalBarGraph_bmp, "MAINTAIN_ASPECT_RATIO", "600", "375")      I am having the following difficulties:   1. I am not able to find a template for this histogram. In the GraphTemplates folder there is PopulationPyramid, but it is still not working. Is there somewhere I can find templates. The create template operation gives a text file instead of a .tee file.   2. I get an error in the line where I addSeriesHistogram saying "Field '10' is not in the dbf"(something like this) while it is only the countBin. Please can someone tell me where I can find good examples of the usage of this function. My dbf file is simple ID(1-999) & simloglkhd(floating point numbers). I really do not know how use this MakeGraph and addSeries_ . I have already read the ArcGIS help. Can you guide me to some more examples and documentation?


title: Create polygon from page size

text: Is there an easy way to create a polygon that matches the data frame in my page layout? I'm trying to set up data driven pages and want to have it key off of polygons that frame the boundaries for each page. I would like to see them as polygons so that I can re-position as necessary and quickly visualize what will appear on each page.


title: How to add LiDAR Intensity attribute to a point feature class?

text: I need to convert LASD (ESRI LAS Dataset) to ASCII text file. I've export the data to point feature class (LAS to Multipoint), converted multipart to single part and calculate X Y and Z.My first question is - How do I add the LiDAR Intensity attribute to the point feature class ? My Second question Is There a better way to convert LASD to ASCII ? i tried the tool las2txt at LasTools but I still looking for an answer for my first question


title: QGIS GPS position after entering attributes

text: I would like to be able to fill a form and, before, during or after, be able to cath a GPS point. I already know the GPS extension but I have to catch the point and enter attributes after... Idea how I can do that? Martin F


title: UK postcodes to census wards

text: I'm planning a project using counts of data aggregated by UK postcodes. I'd like to relate the counts to local small area census population data. Does anyone know is there lookup data available to relate UK postcodes to census wards ?


title: Confidence interval from raster data

text: I have raster data that represents a probability distribution, i.e. each cell has a probability value (in my case the probability that an animal may be found in the cell), and all the cells add up to 100% (I know for certain the animal is within the extents of my raster). I want to be able to generate vector data for confidence values. For example, the 95% line/polygon denotes the boundary wherein I am 95% confident that I will find the animal. Similarly, if I have a kernel density estimate, How do I generate the XX% line/polygon that borders the densest part of the raster containing XX% of the total population? I am willing to use ArcGIS, or open-source software. If there isn't a tool to perform this for me, what is an algorithm I can implement?


title: Exporting PostgreSQL table to FileGDB: wkbUnknown layer geometry type

text: I'm attempting to export a PostGIS (1.5.3) database table as a File geodatabase using the new FileGDB OGR / GDAL 1.9.0 driver:               $ ogr2ogr -overwrite -sql "SELECT * from db_pts LIMIT 1" \         -f "FileGDB" /tmp/test.gdb \         PG:"host='HOST' dbname='DB'"      I get this error:               ERROR 1: Error (-1): FGDB layers cannot be created with a wkbUnknown layer geometry type.     ERROR 1: Terminating translation prematurely after failed translation from sql statement.      The geometry column in the table I'm exported is set up properly and is in the `geometry_columns` table. I think the problem has to do with the way the `-sql` option is returning the layer geometry type. For example, accessing the same database table with `ogrinfo`:               $ ogrinfo -ro -sql "SELECT * from db_allobs LIMIT 1" PG:"host=HOST dbname='DB'"     Layer name: sql_statement     Geometry: Unknown (any)     Feature Count: 1     ...     OGRFeature(sql_statement):0     ...       POINT (-151.472835 69.635245)      It's as though the query doesn't look in `geometry_columns` when building the layer, but later, when exporting the rows it knows that the geometry being exported is a `POINT` (in this case). When I use a similar command to produce a shapefile, all is well, including the layer geometry, which is set to point. Anyone have any ideas?


title: Accessing OpenLayers on "Manage and Install Plugins"

text: I am a new QGIS user and I am trying to add the OpenLayers plugins. I went into 'Manage and Install Plugins', and made sure I had the 'All' tab selected. However, whenever I try to search for the pluggin, it never shows up as on option. In the dialog box I wrote out 'OpenLayers'. Other plugins show up, but not that one.


title: Why is my spatial join count 0 in ArcGIS Desktop?

text: In arcMap, I have a list of locations with lat/long. From a previous post, I learned about the Spatial Join tool to determine the amount of joins that I have. I tried to use the tool with my states layer and it worked perfectly. However, when I tried to do it on the county layer, it is always giving me with 0 joins. What is the cause of this? The county layer that I am using is from http://www.gadm.org/. The one that I used for the states came from arcgis online. There was one for the countries, but it would give me an error when trying to do a spatial join on the data.


title: Is there a way to retrieve/query country data from a shapefile programmatically using map coordinates in GeoTools?

text: I am using Geotools to parse and draw a shapefile. ![Map Drawing](http://i.stack.imgur.com/xemPu.png) Now; when I draw that map I receive coordinates with this bit of code for my `Flight Line` button.               void selectFeatures(MapMouseEvent ev) {                       /*              * Construct a 5x5 pixel rectangle centred on the mouse click position              */                       if(currentClick%2 == 0){             startPoint = ev.getMapPosition();             currentClick++;             }else{              endPoint = ev.getMapPosition();             currentClick++;             }             if(!startPoint.equals(null) && !endPoint.equals(null)){             System.out.println("The two positions are: Start - " + startPoint + "\nEnd - " + endPoint);             System.out.println("Starting point is in ");             }         }      Is there anyway to parse the map/shapefile data for the country that my `startPoint` and `endPoint` are in? ![Show Info Panel](http://i.stack.imgur.com/Udv0M.png) I can see by the `Feature Info` button that was built into the `GeoTools` demo there is a way to parse information from a point; but I am unsure of how to do this programmatically as the `API` is not very well documented.


title: Clipping very large polygon datasets in ArcGIS?

text: I'm attempting to clip a large polygon dataset (~5gb, thousands of features) by a much smaller polygon dataset (~40 features). Is there a best practice or most efficient route for performing this task? The standard geoprocessor clip runs indefinitely on a dataset of this size. Would some form of spatial selection and export be more efficient? **EDIT** : Some great answers below. I selected what I view to be the most thorough response, but each answer provides unique insight into the issue. Thanks!


title: Pan/Scale/Navigation Buttons Disappearing

text: I am working on an OpenLayers 2.10 map using Internet Explorer 8 (used by work), and am encountering some issues with the pan/navigation buttons disappearing after sepcific events:   1. If I change the layer to another background, the pan/navigation panel buttons lose their picture   2. If I resize the browser window, the pan/navigation panel buttons lose their picture When one of these events occurs, the buttons show the broken image link icon. If I refresh the page, the buttons come back. After perusing some answers, the only similar solutions I came up with were at OpenLayers - redrawing map after container resize using the map.updatesize command. These commands did not seem to work, though I am not entirely sure where to place them. Also, I just started using OpenLayers a couple months ago, and a solution may exist using terminology I am not aware of. I would like a solution that keeps the pan/navigation/zoom panel buttons constantly displayed, even if layers are changed or if the page is resized. It is also acceptable for the buttons to refresh after an event that will cause the button to no longer load. Thanks in advance for any guidance.


title: Select parcels from SDE layer using table list and exporting to shapefiles

text: Using a table of county fips codes and state name, I need to select from a large parcel layer stored in SDE and export to individual shapefiles based on the fips code and state name from the table. Can this be done with a Python script or ModelBuilder?


title: Is there demand for a service to publish web maps with lots of shapes?

text: Is anyone in the GIS community trying to publish choropleth web maps with a large (10s to 100s of thousands) number of polygons? and if so what's your current approach and what are the shortcomings of that approach? My company is thinking about reviving an old project, but would like to get a sense of the demand before we invest too much time in something that may not be needed. Features of the service would include:   * Support for millions of polygons.   * Ultra fast updates (recolor the maps in < 1 second, our goal is 60ms).   * Realtime query support (mouseover or click a polygon to display it's info).   * Bring your own data (upload your own shapes, classifications, etc)   * Built in Census Data (all census layers and data preloaded)   * Web access (no coding required)   * API access (REST API with clients in Python and Javascript to start)   * Overlay (display your published maps as overlays in Google Maps, Bing Maps, Open Layers, etc) The core of this system is already working, but it would be a pretty major push wrap it in something useable. We're unsure at this point if there is enough demand to justify the effort. My apologies in advance if this isn't the appropriate forum for this sort of question.


title: Table merging: multiple corresponding objects

text: I'm working with a shapefile which has all the municipalities of France as areas, identificated with their unique code (`INSEE` code). I want to join a table, where every historical monument is included, it has the corresponding `INSEE` column that I may use to merge. But sometimes there can be more than one monument in a municipality (city) identified as its unique `INSEE` code. And when I merge tables together, I would like it to add automatically objects so that every monument is added. If the city of Paris counts 50 monuments for example, I would like to see a pile of 50 "Paris" on my map, each one having a monument. Because when I try to join tables on MaPublisher or Global Mapper, it is not adding objects, it is just merging tables with a unique `INSEE` code. That way, there is only one monument included, because there is no room for others. More understandably, and with simplifications : 1) I have a `CSV` or `TXT` or `XLS`               INSEE;MONUMENT NAME;     10223;Eiffel Tower;     25685;Castle;     25685;City hall;      2) and a Shapefile with objects as areas               INSEE;CITY NAME     10223;PARIS     25685;LYONS      3) I'm getting a shapefile result :               INSEE;CITY NAME;MONUMENT NAME     10223;PARIS;Eiffel Tower     25685;LYONS;Castle      4 ) But what I want is :               INSEE;CITY NAME;MONUMENT NAME     10223;PARIS;Eiffel Tower     25685;LYONS;Castle     25685;LYONS;City Hall      Can you help me and tell me how may I do ? And you are willing to do that for me, I give you the files for you to download : `PROJECTION : LAMBERT 93 (FRENCH)` -it has extra miscellaneous attributes columns- LINK OF THE FOLDER Thank you :)


title: Problems with french coordinate systems (i.e. Lambert Zone I, EPSG 27571)

text: I'm working on a coordinate library implementation and I seem to be stuck with the French Lambert coordinate systems. I'm using this page as a reference (can anybody confirm the results there are accurate?) The individual steps in the conversion seem to fit, but my overall result differs up to 50m from the one computed on that web site. Now step by step: I'm converting as an example the point 2.1°E, 49.8°N from WGS84 to EPSG27571 (NTF Paris / Lambert Zone I)   1. WGS84 -> WGS84 Geocentric: 4122109.79, 151150.90, 4848460.24 (meters)   2. WGS84 Geocentric -> NTF (PM Paris) Datum conversion using Bursa-Wolfe: 4122273.66, 151210.75, 4848135.39 (meters)   3. Geocentric -> Geographic, including meridian shift to Paris: -0.004127387, 0.86917528, -49.58 (radians) I'm quite sure the result is right so far, because these are the same values than said web page returns when using EPSG 4807 as target CS (and converting from grad(!) to radians).   1. Geographic -> Projected coordinates using lambert conformal 1sp, with latitude of origin = 0.86393797974 rad, longitude of origin = 0, False easting 60000, False northing 1200000, Ellipsoid Clarke 1880 IGN: 582974.00, 1233402.04, while expected is 582976.09, 1233397.93. Offset is about 3-5m here.  I have verified the lambert conformal projection using an example for EPSG 24200 JAD96 Jamaica, therefore I _also_ think that the projection itself is right. The error increases the farther away from Paris I get. Not including the false easting or false northing values, the error (quotient between right and wrong values) is in the order of 0.997-0.9998. This looks like some scaling problem, and the WKT description for EPSG 27571 mentions some scale, but I have no clue where to apply that. I do have a scale in the datum, though. Now: Where is the error?


title: I would like to know how to clean topological errors

text: I am a new QGIS 2.0 user. I was so impressed the new default tool, Topology Checker, in this new version. But, hundreds errors are bothering me for cleaning up features. Is it possible to clean errors automatically, not manually? Thank you in advance for your concern. * * * Thank you for your quick response, Mr. Solin and Mr. radouxju. If you mind, let me answer additional 2 questions. 1. Is it possible to clean topological error in QGIS not GRASS. 2. According to "v.clean tool", is it possible to clean overlapped areas. In my understanding, there is no specific tool for overlap error  Thank you so much.


title: How to Map Latitude/Longitude to FAA Chart?

text: How can I map latitude/longitude to an FAA digital chart? Attached is a sample GEOTIFF. More specifically, how can I get a reference points from the digital data? Once I can get reference points, I have no problem translating Pixel-Lat/Long. In my testing, I have been manually identifying the reference points to the the translation. This is a rather tedious process. It looks like I just need a single lat/long-pixel-x/y mapping to derive them all from this data. The GEOTIFF spec suggests the Model Tie Point is the key but the example in the spec does not match this data. -822739.829757951 -64249.12941973 does not give units not does it give offset from what? In this chart the false origin is located at the Pixel coordinates (6521 3865) from the upper left corner. Thanks               ExifTool Version Number         : 9.53     File Name                       : CF-19 WAC 43 South.tif     Directory                       : .     File Size                       : 47 MB     File Modification Date/Time     : 2013:03:15 10:45:54-04:00     File Access Date/Time           : 2014:02:28 13:45:59-05:00     File Inode Change Date/Time     : 2014:02:28 13:06:22-05:00     File Permissions                : rwxr-xr-x     File Type                       : TIFF     MIME Type                       : image/tiff     Exif Byte Order                 : Little-endian (Intel, II)     Image Width                     : 11888     Image Height                    : 4135     Bits Per Sample                 : 8     Compression                     : Uncompressed     Photometric Interpretation      : RGB Palette     Strip Offsets                   : (Binary data 286 bytes, use -b option to extract)     Samples Per Pixel               : 1     Rows Per Strip                  : 128     Strip Byte Counts               : (Binary data 262 bytes, use -b option to extract)     Planar Configuration            : Chunky     Color Map                       : (Binary data 1536 bytes, use -b option to extract)     Sample Format                   : Unsigned     Pixel Scale                     : 127.09668 127.09668 0     Model Tie Point                 : 0 0 0 -822739.829757951 -64249.12941973 0     GDAL Metadata                   : <GDALMetadata>.  <Item name="PyramidResamplingType" domain="ESRI">NEAREST</Item>.  <Item name="RepresentationType" sample="0">THEMATIC</Item>.  <Item name="STATISTICS_MINIMUM" sample="0">0</Item>.  <Item name="STATISTICS_MAXIMUM" sample="0">255</Item>.  <Item name="STATISTICS_MEAN" sample="0">108.47642098115</Item>.  <Item name="STATISTICS_STDDEV" sample="0">83.473159026412</Item>.</GDALMetadata>.     GDAL No Data                    : 256     Geo Tiff Version                : 1.1.0     GT Model Type                   : Projected     GT Raster Type                  : Pixel Is Area     GT Citation                     : PCS Name = CF-19 WAC     Geographic Type                 : NAD83     Geog Citation                   : GCS_North_American_1983     Geog Angular Units              : Angular Degree     Geog Semi Major Axis            : 6378137     Geog Inv Flattening             : 298.257222101     Projected CS Type               : User Defined     Projection                      : User Defined     Proj Coord Trans                : Lambert Conf Conic 2SP     Proj Linear Units               : Linear Meter     Proj Std Parallel 1             : 46.6666666666667     Proj Std Parallel 2             : 41.3333333333333     Proj False Origin Long          : -69     Proj False Origin Lat           : 44.41667     Proj False Origin Easting       : 0     Proj False Origin Northing      : 0     Image Size                      : 11888x4135


title: Geocoding house numbers

text: I have a point feature class with residential information, but all that it has in terms of geographic information is the residence's house number. There is no other information (street, zip, etc). JUST house numbers. The points fall more or less where they are located, but they are somewhat off. I'm working off of a ArcMap basemap and I sometimes I can just move the points to areas that are obvious to me, for example if there's one point in a general area close to a house, and that's the only house, then its obvious to me that that point is that house. However, i need a faster way to do this, and I need help with points that aren't so obvious, like for example when i have to assign a point to a residence and there are several residences to choose from. I thought about geocoding, but I'm not sure how this would work when all I have to work with is a house number. So I thought I could download a shapefile (maybe tiger data) with street names as well as zip codes, and do some kind of spatial join to where the points (house numbers) would be assigned the street name and zip code nearest the point. and then create an address locator to match the points to the actual locations based on those attributes (street name, house number, and zip code). Would this be a good way to tackle this? Does anyone have any other suggestions that might be better? Could you maybe explain to me how I might do this? I'm pretty new to GIS and not very familiar with geocoding.   Also, can someone explain to me more or less how I should perform this spatial join? Thanks!!


title: Calculate lat lon bounds for individual tile generated from gdal2tiles

text: I have multiple tile sources generated using gdal2tiles I'd like to present on the same map. So when presenting a tile I need to determine which source to serve it up from, checking against the bounds. Anyone know how to calculate the lat lon bounds for a single tile based purely on zoom, x, and y (from the file structure generated by gdal2tiles)? BTW: I'm using Google Maps API v3 in case there's a need to call some functionality from the API to help with calculations. The reason I'm restricted to purely zoom, x, and y is because the tiles aren't just getting called by an overlay on the map but also some custom print functionality that allows printing outside of the map. Thanks in advance for any help or suggestions anyone has. Cheers, Gavin


title: How to load Spatialite extension via SQLite3 ODBC Driver?

text: I am having trouble getting the Spatialite extension to load as part of the SQLite3 ODBC driver...[aka the joys of making software somewhat spatial when it tries so hard not to be]. **What works [loading extension into SQLite Driver]:** I have software that is connected using the SQLite3 ODBC driver (http://www.ch-werner.de/sqliteodbc/). I have previously loaded extensions into this driver without any problems. A previous working example was a trigger that updated columns [using triggers] using some trig functions that are not part of base SQLite functionality - I built the *.so extension and pointed SQLite to it as part of the ODBC connection and everything worked fine - when I connect/export data from the software I got no issues. **What works [loading Spatialite extension into SQLite]:** If I use the SQLite shell I can load the Spatialite (or any other extension) using: select load_extension('libspatialite-2.dll');. That works without fail. **What doesn't work:** However I cannot seem to get the Spatialite extension to load properly as part of the ODBC connection - and I am using a trigger that builds the geometry column of a Spatialite database based on the software providing only x and y values. I get an error that says "Extension C:\windows\system32\libspatialite-2.dll did not load: the specified module could not be found". I also tried the spatialite-4.dll and it doesn't work either. The module and its requirements are located there and the environmental path has (always) been there. If I do not load the extension then I will get an error when the trigger calls a function (makegeometry, makepoint or whatever you want) that is not part of the SQLite driver - it doesn't understand the function. My hack workaround at this point are batch files that run what the trigger should do after the data has transferred from the software to the Spatialite database, but I'd prefer it to be a one-stop shop. Any thoughts would be great. Thank you. I haven't seen this much on the web, other than a Google Groups question (https://groups.google.com/forum/#!topic/spatialite-users/3ImfXm7cdpU) that didn't resolve the problem [I haven't tried to build the CppDB yet - I want to see if there is a more ~~lazy~~ straightforward solution].


title: How to perform Intersection using ModelBuilder iteration?

text: ![Pic of Model](http://i.stack.imgur.com/UHZrG.jpg)I would like to take the intersection of two areas (one randomly generated), and do this n times. I can do it and get one result, but I can't save all the intersection areas to a file at the end. The problem is that when I run my model more than once it says that my buffer output file already exists, which I can imagine would happen again for the intersection. What I'd like to happen is to run my model where I produce new random points everytime, then buffer them and intersect the dissolved buffer with another area and then get out the area. Then I want to do this several times. If this isn't worded very well then just let me know and I can explain what I mean in another way. Thanks a lot


title: SQL: Select aggregated Bounding Box of multiple Features

text: I would like to calculate the aggregated bounding box of multiple features using OGC Simple Feature SQL or SQL API of ArcSDE. The following statement will only query the bounding box of each feature:               select ST_Boundary (shape) from polygon_table where country like '%usa%'      What I needed is one bounding box aggregated over all features, resulting to the condition in the where clause.


title: Why am I unable to download Landsat 8 geotiff data from the USGS site?

text: I seem to have an issue with not being able to download geotiff Landsat 8 data from te USGS site. After i find the image i would like to download and proceed to download the Level 1 geotiff option is disabled and I am only allowed to download the LandsatLook jpeg images. I downloaded data from this site before but perhaps I am overlooking something. I tried both http://glovis.usgs.gov/ and http://earthexplorer.usgs.gov/ Screen capture showing the disabled "dimmed" option: ![enter image description here](http://i.stack.imgur.com/dtLjo.jpg)


title: ArcSDE: connecting via LDAP Authentication

text: I was wondering if anyone has a sample connection string for SQLPlus using LDAP Authentication. I am also trying to specify in the connection string the database I wish to use. Thanks, CDB


title: Using pyqgis with existing ArcGIS files

text: I'm pretty new to GIS. I have some files that seem to have been created using ArcGIS (I got them from here http://www.fao.org/ag/againfo/resources/en/glw/GLW_dens.html). I would like to do some analysis on these maps, preferably using python. I'm wondering if it's possible to use pyqgis to do this kind of thing. I'm really just looking for some info to get me started. e.g. - Is it possible to import these files into a pyqgis application? - Are there any tutorials that might help me get started? Thanks. Edit: The file types I see are .lyr, .aux, .avl, .dat, .nit, .adf Here's a specific example of what I'd like to do. Say I have two maps (map A and map B). I think these are raster images, and let's say that the "colour" in a given cell in map A measures quantity qA. And likewise for map B. Something simple I'd like to do is to produce a third map (map C) whose cells have quantities given by qC = qA*qB.


title: Control Map Pan per checkbox

text: i want to add the functionality of switching the option to pan the map and put it in a checkbox. But im not quite sure how i can do that. The reason I want to do this is: I have a layer which displays features and has a select and hover control. So theres no way to pan the map, because the clickEvent triggers direct action from the Handlers (click,hover). I want to add a checkbox where i can activate/deactivate the possibility to pan the map, so i can pan the map when there are features present, which normally would trigger a reaction to the mouseover the feature. The way ive tried to do it is this one but it didnt work. I attached the function to the checkbox but it never worked :/ ` function allowPan() {               if(panAct == false){             highlightCtrl.activate();             selectedControl.activate();             alert("controls activated!!")             panAct = true;         }          if(panAct == true){          highlightCtrl.deactivate();         selectedControl.deactivate();         alert("controls deactivated!!");         panAct = false;         }      } `


title: sa.fill creates NoData raster in standalone script ArcMap 10.0

text: I am trying to put together a geoprocessing workflow in Python, and I've hit a stumbling block with sa.fill. When run as a standalone script, the tool produces a raster composed entirely of NoData values. If I run the same code in the Python window in ArcMap, a correctly filled raster is produced. My code snippet is:               outfill = sa.Fill(raster,"")     outfill.save(fill)     print "Sink fill complete."      I've tried copying my entire script into the Python window and adding it to a toolbox, but I get the following error in both cases:               Runtime error <type 'exceptions.EOFError'>: EOF when reading a line      The standalone script does not produce this error. It actually runs through all the steps without incident, but the output is useless as everything after the Fill step is based on NoData values. Does anyone have any thoughts on this?


title: Creating new shapefile based on its attributes with python

text: I have a shapefile which has an attribute year. I used to generally copy and paste these features in the dataframe and then using the definition query just show polygons which were established before that year. I want to know how can I do it using python. Like if I can give the input feature class, add all the years I want the different shapefiles for and the workspace where it can create shapefile for the mentioned years with the right naming convention. Thank You.


title: How to calculate polygon centroids in R (for non-contiguous shapes)

text: I've spent a little while figuring out the answer to this question. It's not immediately obvious from a Google search, so thought it may useful to post the answer on here. There is also an additional question about non-contiguous polygons. Instant easy answer: **use the command:**               centroids <- getSpPPolygonsLabptSlots(polys)      (This was found in the class description of the SpatialPolygonsDataFrame R data class for the overarching spatial package in R, sp) This seems to do exactly the same thing as               cents <- SpatialPointsDataFrame(coords=cents, data=sids@data, proj4string=CRS("+proj=longlat +ellps=clrk66"))      in the following code, which should be replicable on any R installation (try it!)               #Rcentroids     install.packages("GISTools")     library(GISTools)     sids <- readShapePoly(system.file("shapes/sids.shp", package="maptools")[1],                            proj4string=CRS("+proj=longlat +ellps=clrk66"))     class(sids)     plot(sids)     writeSpatialShape(sids, "sids")     cents <- coordinates(sids)     cents <- SpatialPointsDataFrame(coords=cents, data=sids@data,                        proj4string=CRS("+proj=longlat +ellps=clrk66"))     points(cents, col = "Blue")     writeSpatialShape(cents, "cents")          centroids <- getSpPPolygonsLabptSlots(sids)     points(centroids, pch = 3, col = "Red")      Where cents (blue) and centroids (red) are identical centroids (this should plot should appear after you've run the code): ![centroids calculated by R](http://i.stack.imgur.com/qpCm7.png) So far so good. But when you calculate polygon centroids in QGIS (menu: Vector | Geometry | Polygon Centroids ), there are slightly different results for non-contiguous polygons: ![QGIS generated polygons](http://i.stack.imgur.com/zY5qq.png) So this question is 3-things:   1. A quick and easy answer   2. A warning for people using R to calculate centroids for non-contiguous polygons   3. A question about how it should be done in R to properly account for multi-part (non-contiguous) polygons


title: building opentripplaner gtfs.zip

text: I am building a journey planner app using opentripplanner. i have followed their two minutes tutorial in which we just run webapp which uses pre build graph.obj file, it works fine. But in next five minutes tutorial they shows how can i build .obj from a zip gtfs, while i try to do so i got these errors, i try had to figure them out ! ![enter image description here](http://i.stack.imgur.com/FGDrF.png) ![enter image description here](http://i.stack.imgur.com/qkG0T.png) ![enter image description here](http://i.stack.imgur.com/GkjdW.png)


title: Why doesn't Polygonizer copy the line ids to my polygonized shapefile?

text: I am trying to convert a line shapefile into a polygon shapefile using Polygonizer plugin, and when my new polygons are created all the cells in my id field are "Null", although I have clicked "copy table field from line layer". I have tried both Polygonizer 2 and Polygonizer 2.1. Any idea why the ids from my line layer are not copied to my polygonized shapefile? Thank you!


title: Problem looking at tables within ESRI geodatbase

text: I installed the most recent OSGeo4W with QGIS 1.8.0 and GDAL 1.9.2 on my Windows 7 PC. I can load an ESRI geodatabase (.gdb file) and I can get the window "Select vector layers to add". When I click on any of the layers, however, I get the message "qgis.exe has stopped working" Any idea where the problem might be?


title: How to get the nonspatial tables in the TOC into a combobox in ArcMap?

text: I'm new in using Arcobject . i'd like to set dbf of ArcMap TOC to the Combobox Items. How can i set dbf files to items of combobox using C# programming and Arcobject ?


title: How to add a DEM Band to an existing 3-band TIFF file using QGIS?

text: I'm quite new in the GIS business and I do not figure out how to do the following thing: I have two geoTIFF Files, one is a 3+1 bands file (RGB+alpha), the second one is a DEM geoTIFF file containing one band with elevation data. For a certain application, what I need is a 1+3+1 bands geoTIFF, containing the elevation Data on the first band, and the optical data plus the alpha info in the 2nd to 4th bands. How can I join or add bands from different geoTIFF's using Quantum GIS, for instance? I have found a similar question here, but there does not seem to exist a definitive answer there. So, I would be very happy to get help from the experts on this platform, Thanks in advance.


title: Geoserver - problem with MySQL support (impossible to create MySQL store at all)

text: I need to access spatial data in MySQL database with Geoserver. Official user manual states: just download package, extract it into WEB-INF/lib directory and all will be OK. Instead, just nothing happens - MySQL option doesn't appear when creating new data source, I can't find any relevant messages in logs and Google doesn't provide a hint. What do I do:   1. Download MySQL extension, following recommended link. I have Geoserver version 2.4.3 (uses jre7, downloaded recently). Downloaded archive declares to be of version 2.4.3 (and manual specifically states that versions should match).   2. Unzip archive and copy it contents to `[Geoserver installation directory]\webapps\geoserver\WEB-INF\lib` (also tried other paths mentioned in the manual and extension readme, though they don't make any sense - result the same).   3. Restart Geoserver service (also tried rebooting system), than go to _Data -> Stores -> Add new store_ in Geoserver web interface (logged in as admin). Manual states that MySQL option should be there, but it isn't. As mentioned, I use Geoserver 2.4.3 (latest at current time), OS is Windows 7 x32. Geoserver is installed as a service and starts automatically. Logging profile of server is set to VERBOSE_LOGGING, no traces of "mysql" are found. System logs don't contain anything but messages from Service Control Manager about service being started/stopped. Trying to mess with other JARs in the same directory causes Geoserver to fail. There was installed Geoserver 2.1-beta3 some time ago, but current installation is made into different directory. I'm interested in any advice how to get Geoserver-MySQL combination working, or how to diagnose why it is problematic/impossible.


title: gtfs-editor compile time error

text: I'm trying to build the gtfs_editor project by following the instructions on the given page. But, while executing the command **play run** the log at `localhost:9000` says `Compilation error The file /app/controllers/Check.java could not be compiled. Error raised is : The type Check is already defined` Has anyone gone through this problem before? If yes then your help on this issue is much appreciated. Thanks in advance. p.s.- I'm working on `play framework -1.2.5.3`,


title: Direct Georeferencing with X,Y,Z and Omega, Phi, Kappa

text: I have images taken from plane, it's longitude, latitude, elevation, omega, phi, kappa. I suppose that earth is flat (I don't need any DEM on that stage) and I want to calculate how image is nailed to the ground. I want to find center and 4 corners of nailed image to the ground. Could you advice some simple formula for calculation? Or free tool that has such calculations just to check my own implementation results. I tried to use GRASS but without any luck, it is too complicated :( Thanks in advance! **UPDATE:** Images were taken from around 3.2k meters without physical GCPs. The precision that we got is around 20 meters with camera angles (Omega, Phi) around 2 grad. Do you think it is possible to get better precision? We are using C# and here is our implementation:                   /// <summary>         ///          /// </summary>         /// <param name="omega">camera angle, in grads</param>         /// <param name="phi">camera angle, in grads</param>         /// <param name="kappa">camera angle, in grads</param>         /// <param name="longitude">camera position</param>         /// <param name="latitude">camera position</param>         /// <param name="altitude">camera position</param>         /// <param name="altOverGround">average flight height, it is defined on flight planning stage</param>         /// <param name="imageLength">it is defined on flight planning stage, calculated from altOverGround and focal lenght, it is in meters</param>         /// <param name="imageWidth">it is defined on flight planning stage, calculated from altOverGround and focal lenght, it is in meters</param>         /// <returns></returns>         public static DbGeography CalculateProjection(double omega, double phi, double kappa, double longitude, double latitude, double altitude,             double altOverGround, double imageLength, double imageWidth)         {             // coordinates of the camera             DbGeography cameraCoordinates = CreatePoint(longitude, latitude);                  // converting angles from grad to radian             omega = ConvertFromGradToRadian(omega);             phi = ConvertFromGradToRadian(phi);             kappa = ConvertFromGradToRadian(kappa);                  // get image center (projected to groud) as offset to camera coordinates             // result of imageCenter is in image coordinate system             Point imageCenter = GetImageCenter(omega, phi, kappa, altitude);                  // calculating real image length for this height             var realImageLength = imageCenter.Z * imageLength / altOverGround;             var realImageWidth = imageCenter.Z * imageWidth / altOverGround;                  // corners offset in image coordinates system             var xOfsfet = realImageLength / 2;             var yOffset = realImageWidth / 2;                  // corners list             List<Point> imageCorners = new List<Point>() {                  new Point(xOfsfet, yOffset, imageCenter.Z),                 new Point(-xOfsfet, yOffset, imageCenter.Z),                 new Point(-xOfsfet, -yOffset, imageCenter.Z),                 new Point(xOfsfet, -yOffset, imageCenter.Z)             };                  List<DbGeography> projectionCorners = new List<DbGeography>();             foreach (var point in imageCorners)             {                 // get image corner projected to the ground                 Point projectedPoint = GetProjectionCoordinates(omega, phi, kappa, altitude, point);                      // add projected point to list of projected corners                 projectionCorners.Add(CreatePointWithOffset(cameraCoordinates, projectedPoint.Y, projectedPoint.X));             }                  // add first point to complete the polygon             projectionCorners.Add(projectionCorners[0]);                  return CreatePolygonFromListOfPoints(projectionCorners);                     }              private static Point GetProjectionCoordinates(double omega, double phi, double kappa, double altitude, Point imageCorner)         {             // converting point coordinates from image's coordinate system to normal coordinate system             Point imagePoint = ConvertImageCoordToNormal(omega, phi, kappa, imageCorner);                  // get coordinates of the projected point             return GetProjectionCoordinates(altitude, imagePoint);         }              private static Point GetProjectionCoordinates(double altitude, Point imagePoint)         {             var x = altitude * imagePoint.X / imagePoint.Z;             var y = altitude * imagePoint.Y / imagePoint.Z;             return new Point(x, y, altitude);         }              private static Point ConvertImageCoordToNormal(double omega, double phi, double kappa, Point imageCorner)         {             double x1, x2, y1, y2, z1, z2, x, y, z;             kappa = Math.PI / 2 - kappa;             x1 = imageCorner.X * Math.Cos(kappa) + imageCorner.Y * Math.Sin(kappa);             y1 = -imageCorner.X * Math.Sin(kappa) + imageCorner.Y * Math.Cos(kappa);             z1 = imageCorner.Z;                  x2 = x1 * Math.Cos(phi) - z1 * Math.Sin(phi);             y2 = y1;             z2 = x1 * Math.Sin(phi) + z1 * Math.Cos(phi);                  x = x2;             y = y2 * Math.Cos(omega) + z2 * Math.Sin(omega);             z = -y2 * Math.Sin(omega) + z2 * Math.Cos(omega);             return new Point(x, y, z);         }              /// <summary>         /// Calculates altitude of the image center         /// </summary>         /// <param name="omega"></param>         /// <param name="phi"></param>         /// <param name="kappa"></param>         /// <param name="altitude"></param>         /// <returns>Altitude of image center</returns>         private static Point GetImageCenter(double omega, double phi, double kappa, double altitude)         {             // calculating coordinates of image center             var xCenter = altitude * Math.Tan(phi);             var yCenter = altitude * Math.Tan(omega);             var imageAltitude = altitude / (Math.Cos(omega) * Math.Cos(phi));                  // recalculating center with kappa value             var alpha = Math.Atan(yCenter / xCenter) + kappa;             var aa = Math.Sqrt(xCenter * xCenter + yCenter * yCenter);                  double imageCenterX = Math.Cos(alpha) * aa;             double imageCenterY = Math.Sin(alpha) * aa;             return new Point(imageCenterX, imageCenterY, imageAltitude);         }              private static double ConvertFromGradToRadian(double omega)         {             omega = Math.PI * omega / 200;             return omega;         }      and Point class:               class Point     {         public double X { get; set; }         public double Y { get; set; }         public double Z { get; set; }         public Point(double x, double y, double z)         {             X = x;             Y = y;             Z = z;         }     }


title: Print FeatureLayer labels in Esri JavaScript API 3.5 and earlier

text: I have a need to implement labels on features in a FeatureLayer in version 3.5 if Esri's Javascript API. The labels are from a field returned by a REST feature service. I can't move to 3.7 for various reasons at this time. I have tried using a TextSymbol but my map features just turn to the color of the TextSymbol and no text appears. I may be approaching this in he wrong manner, though. Below is the code I'm attempting to use for labeling with the featureLayer object being my instance of the FeatureLayer I'm adding to the map. Is there a different or proper way to accomplish this task?               featureLayer.on("graphic-add", function (evt) {      var labelColor = new Color([255, 0, 0, 0.25]);     var myLabel = new TextSymbol(evt.graphic.attributes["My Field Name"]);          myLabel.setColor(labelColor);     myLabel.font.setSize("14pt");          evt.graphic.setSymbol(myLabel);          //console.log(evt);               });      Thanks for any help that can be provided!


title: How can I make use of the $scale function in QGIS data defined properties?

text: In QGIS 2.0 I can use with success the $scale function in the labeling expressions. When I tried to use it for symbolizing in the data defined properties it has a value of 0. Is there a way to make it perform like it does in labeling (getting the scale)? As I saw in the expression string builder help that               $scale function     Returns the current scale of the map canvas.      Note: This function is only available in some contexts and will be 0 otherwise.      is there a source where I can read more about this e.g. in what contexts it is available? Thanks, Tudor


title: Are there sources for obtaining pre-made KML data from Google?

text: I have been asked to investigate methods for obtaining KML files. All I know so far is that it is a format based on XML that Google uses to format their maps (or the metadata). However, are there places or methods that I may avail myself of to download KML data without going through Google Earth? For example, is it possible to get the KML for the vector data of a state in the US (roads, boundaries, etc.) ? To restate it simply: I want Google's vector data so I can plop it into ArcMap. Are there repositories for such things? I doubt it, since Google never struck me as being all handy dandy open source randy, but one can hope.


title: Torque cumulative creates grid

text: How can I get Torque's cumulative setting to draw the data points in their correct locations versus artificial grids? http://cdb.io/OAgzRD My dataset of global locations looks fine under non-cumulative setting, but cumulative wants to add the locations in a grid that looks very artificial. C


title: Toggle PopUp Functionality

text: Firstly: I've look through every example on Stack GIS.Stackexchange which might deal with my problem. Of course there was good input, but unfortunately I'm not able to fix it. Admittedly my Java Script skills are not the best, so "easy" thing are also very hard to overcome. Until now I've stuck at the following level: http://jsfiddle.net/84vvR/ I'm using nliu's TOC http://gmaps-utility- gis.googlecode.com/svn/tags/agsjs/latest/examples/toc.html and the editor and popup functionalities overed by esri's API. As you can see, theoretically all work. The mapapp contains an editor an the option for popup. The problem is: If the user wants to add new features by clicking on the map, the popup pops out and gives information about the underlying area, which has the defined popup. The same problem, if the user want to use the measurement tool. On the first click, the popup pops out and so it's not possible to add a starting point for measurement. I've tried to implement a workaround with dojo.connect and a clickHandler or onClick like (button dojoType=”dijit.form.Button” onClick=”....(),...();”>Disable Popup Function I assume that it's not too hard for an advanced developer to solve this problem to bind the Popup Functionality to a checkbox or Button. Would be awesome to get some help. Thanks a lot!


title: Why is the Edit Vertices button (Draw Toolbar) enabled for some annotations but disabled for others?

text: I keep running into this:   * Export Labels in Layer A & B to annotations (In Map, In current extent only)   * Enable Focus Data frame   * Click on annotation exported from Layer A    * **Click the Edit Vertices button - No problem**   * Click on annotation exported from Layer B   * **Edit Vertices button is DISABLED** What's the problem? When I export the same labels to geodatabase annotations, I can easily edit the vertices on any annotation. Figure 1: Click on Layer A annotation - Edit Vertices button enabled ![enter image description here](http://i.stack.imgur.com/Gcibc.png) Figure 2: Click on Layer B annotation - Edit Vertices button DISABLED ![enter image description here](http://i.stack.imgur.com/qlQb3.png)


title: Rotate polygon preserving topology in QGIS?

text: Is it possible in QGIS? What I need is: 1) Rotate a polygon inside another one (bigger), obviously the polygons are not overlapped. 2) Rotate a given polygon in a topologically correct layer (like a cadastral parcels layer) What I would like to avoid is make holes, merge polygon pieces, etc. (the "old way"). Thanks in advance.


title: How to find the locations of the highest values in a raster?

text: Using ArcGIS 10, I have a raster where I would like to find the pixel with the maximum value in the raster and return its location (center of pixel) in decimal degrees. I would like to iterate through this process returning the location of the second highest value of the raster, then third, and so on so that in the end I have a list of N locations that have the highest values in the raster in order. I imagine that this might be most easily done using a Python script but I'm open to other ideas if there is a better way. Thanks!


title: Set Numeric Format for Class Breaks Labels?

text: Why doesn't this code change the numeric format of the labels generated when I build this renderer? When I execute the set_Break() this nnumeric format is correctly set in the renderer but I still get numbers that look like "110.400000" in the generated label.               INumericFormat numericFormat = new NumericFormatClass();     numericFormat.AlignmentOption = esriNumericAlignmentEnum.esriAlignLeft;     numericFormat.AlignmentWidth = 0;     numericFormat.RoundingOption = esriRoundingOptionEnum.esriRoundNumberOfDecimals;     numericFormat.RoundingValue = 2;     numericFormat.ShowPlusSign = true;     numericFormat.UseSeparator = true;     numericFormat.ZeroPad = false;     IClassBreaksUIProperties classBreaksUIProperties =     classBreaksRenderer as IClassBreaksUIProperties;     classBreaksUIProperties.NumberFormat = numericFormat as INumberFormat;     classBreaksRenderer = classBreaksUIProperties as IClassBreaksRenderer;


title: Packaging failed in ArcGIS Server 10.2.1

text: I use ArcGIS Server 10.2.1. on Windows 8 64 bit. ![Packaging failed](http://i.stack.imgur.com/P5bve.png) I'm getting the simple error: "Packaging failed" for each new service I try to publish. Any thoughts what might be causing this? Thanks


title: How to find network distances between two sets of origins and destinations?

text: I'm learning how to use Network Analyst, and I'm not entirely sure what to do (or search for). My objective is to measure the nearest distance from the red circles to the blue circles along a road network: ![enter image description here](http://i.stack.imgur.com/9nNjB.png) I have made a network in ArcCatalog using the network wizard and I'm just interested in the simplest network distance. My network looks like this: ![enter image description here](http://i.stack.imgur.com/h70Bt.png) Is this calculation possible to do? My red points are generated as centerpoints of a grid as I am just interested in the average distance for each grid. I realize that these red points are not on the network, so I'm not sure if that is a problem for the calculation. What should be my next steps? (after reading more of the help documentation!)


title: What is the best way to do reverse geocoding with PostGIS?

text: Is there a way to do reverse geocoding with PostGIS? There's a lot of information on the Web but I found it very confusing... I have imported osm data into my own Postgres with PostGIS database. Now I want to make geocoding and reverse geocoding with my database. I don't know if there are some Stored Procedures, or a project somewhere that takes care of my issue. The next step that I want to do is make a webservice that can handle this kind of requests, but that is another question. Thanks in advance.


title: Is possible to Unit Testing using ArcGISLocalFeatureLayer ?

text: DEAR PEOPLE FROM THE FUTURE: Here's what we've figured out so far… I am developing an application using ArcGIS Runtime for WPF, I want to implement unit test for the functionality in several classes, but some of the workflow includes the use of ArcGISLocalFeatureLayer Class. Trying to analyze the process I came to the conclusion that use this class is going to be an integration test rather than a unit test, due the interaction with the local ArcGIS server that is going to be created. Someone has been able to mock this kind of objects? (please...tell me that you use Rhino Mocks).


title: Finding parcel polygon features with frontage on multiple streets

text: The Department of Public Works would like to know how many properties in the city have a sidewalk at both the front and back. To get at this, I'm looking for a method to find parcels with frontage on more than one road. I'm working in ArcInfo 9.3.1 (file geodatabase). I'm looking for a solution that doesn't require scripting or downloading other tools. Street centerlines and parcel polygons are available; there is no sidewalk layer. Using some of the ideas from this forum, I have developed and tested a process. Some of the results are puzzling (corner lots may or may not be included in the final results-I think DPW doesn't want corner lots included) and I'm concerned about scaling it up to the full 42000 parcels (I tested 1223 parcels). Are there steps to add to this method to improve the results? Is there another way to accomplish the same thing that could be more automated through model builder? • Run Feature to Line tool on parcel polygons feature class (two lines are created for boundaries shared by polygons and lines retain parcel id) • Build a geodatabase topology with this newly created parcel line feature class using the rule “must not have overlaps” then validate topology (interior lines will be in error) • In ArcMap, add topology and feature class and open error inspector then search for errors > select all errors > right click > select features > delete (to delete interior lines) • Run Frequency on remaining features in parcel line feature class with parcel id as frequency field – select records in table where Frequency > 1 and export to new table • Join Frequency>1 table to parcel polygons based on parcel id, keeping only matching records – export join results


title: How to close and fill Gulfs and inlets in QGIS or Grass?

text: I have a raster layer from which I extracted only zones with a certain color, then polygonized it with QGis. I just generalized the new vector layer in GRASS to smooth things a bit. There are many holes, "gulfs", "inlets" and cuts because the color zone in the raster overlays a black and white cartographic map. I would like to fill this spots, but I can't find a easy way. I already closed valid holes, but with "gulfs" and "inlets" I cannot use the same command. ![expample of cuts and inlets](http://i.stack.imgur.com/HuCvz.png) In that image, only the hole on the right is suposed to exist. Is there a easy way to solve this?


title: Why my annotation markers are printed smaller?

text: I annotated one of my maps, for every annotation I set an svg marker whose dimension is either 6 or 8 mm. Everything looks fine in the main window. The problem arises when I try to print my map using the composer, it somehow scales down my markers to a miniscule dimension, not respecting the dimension in mm. This process is glitchy as you can see in the following screenshot: ![composer wrong dimension screenshot](http://i.stack.imgur.com/DZcXN.png) You can notice this particular svg symbol is not drawn correctly: a small portion of it (top left) seems scaled up; this single issue doesn't appear with the other svg markers in my map. EDIT: I just verified that this problem is specific to svg markers, other kinds of indicators are printed just fine


title: GeoServer layer preview shows only missing images for tiles

text: Just getting started with GeoServer and already running into problems. When I click on Layer Preview and view any of the preloaded data in OpenLayers (or any other format) it is just blank and has the missing image icon where the data should be. One interesting thing is that I could click randomly in the blank map and see attributes of selected polygons appear so I think the data is partially loaded. Example screenshot. I've installed it on multiple machines and I am getting the same exact problem. I have connected my own data through PostGIS but I got the same result. That is why I am assuming it is user error. Anyone have any ideas?


title: Multiple field vbscript query

text: I'm fairly new to arc gis and I'm trying to string along several queries. individually they work in same query they don't. background I want to query all permit types rc, re, re/ap, and db; less than 30 days old, and only in certain counties. I'll be making 3 layers with ranging dates <30, 30-60, 60-90, >90 whats working Layer 1: "PERMITTYPE" ='RC' OR "PERMITTYPE" = 'RE' OR "PERMITTYPE" = 'RE/AP' OR "PERMITTYPE" = 'DB' AND "DATEFILED" >= CURRENT_DATE - 30 ** This returns about 60 results in the below counties Layer 2: "COUNTY" IN ( 'BORDEN (TX)', 'GARZA (TX)', 'LYNN (TX)', 'KENT (TX)', 'SCURRY (TX)' ) Layer 3 (looking for different permit type but works) "PERMITTYPE" = 'OD' AND DATEFILED >= CURRENT_DATE - 30 AND "COUNTY" IN ( 'BORDEN (TX)', 'GARZA (TX)', 'LYNN (TX)', 'KENT (TX)', 'SCURRY (TX)' ) New layer (not working ) "PERMITTYPE" ='RC' OR "PERMITTYPE" = 'RE' OR "PERMITTYPE" = 'RE/AP' OR "PERMITTYPE" = 'DB' AND "DATEFILED" >= CURRENT_DATE - 30 AND "COUNTY" IN ( 'BORDEN (TX)', 'GARZA (TX)', 'LYNN (TX)', 'KENT (TX)', 'SCURRY (TX)' ) any help would be awesome thanks in advance


title: Workflow with raster in Qgis 2.2

text: I'm trying to build a workflow with the Qgis 2.2 processing toolbox. I want to clip a raster layer with a vector layer using the GDAL tool "clip with a layer mask". But it doesn't work. The error message in the historic-info is : "GDAL execution console output ERROR 4: `0' does not exist in the file system". The error message mentionned in the historic-algorithm is : "processing.runalg("modeler:model_clipp_raster","D:\Mes Documents\Documents\Qgis processing\Test_LaFleche\la_fleche.shp","D:\Mes Documents\Documents\Qgis processing\Test_LaFleche\72_SCAN_25\SCAN_25_Topo\Extrait_72_Topo.tif",None)". I've done the treatement in the Qgis "classical" interface without problem. It seems to me that problem stand in the fact that the raster can't be loaded. I hope i'm clear enough. If anyone has an idea about the source of my problem and the way to fix it, i would be really grateful. Best regards, Emily.


title: Problems running model in QGIS

text: I am trying to make a batch model to convert band 2 and band 1 MODIS 8 day composites into NDVI. The input files are in tif format, but I want the output to be saved in Idrisi .rst format. I am using QGIS 2.2 64 bit on a Windows 8.1 machine. I installed via the OSgeo4W, and installed everything, so I am using the Orfeo-toolbox for the raster calculator. I am new to QGIS and not very experienced in python either. The simple model is presented below. Ultimately, I want to batch this, but for the time being I’m just testing to see if I can get a simple model to work. I have tested the data, using the raster calculator directly in QGIS to make an NDVI and this worked without any problems. I just have a lot of images I need do this for.... ![screen capture of mode in QGIS][1] Below I have copied the OTB execution command and the execution console output from QGIS. My questions are: What does it mean that I am missing mandatory parameter “il 0”? Is there a better way to make NDVI images than using the model I have made? Any help or advice on this matter is greatly appreciated! Thanks, Brian ++++++++++++++++++++ OTB execution command C:/OSGEO4~1\bin\otbcli_BandMath -ram 128 -exp (NIR - RED) / (NIR + RED) -out "C:/8day_NDVI/NDVI_1.rst" +++++++++++++++++++++++++ OTB execution console output ERROR: Missing mandatory parameter: il 0 ERROR: Troubles loading parameter, please check your line argument... This is the BandMath application, version 4.0.0 Perform a mathematical operation on monoband images Complete documentation: http://www.orfeo- toolbox.org/Applications/BandMath.html Parameters: -progress Report progress  MISSING -il Input image list (mandatory) -out [pixel] Output Image [pixel=uint8/uint16/int16/uint32/int32/float/double](default value is float) (mandatory) -ram Available RAM (Mb) (optional, off by default, default value is 128) -exp Expression (mandatory) -inxml Load otb application from xml file (optional, off by default)


title: arcpy open up a script tool window once user save edits

text: On our map, after user draw a polyline on map and click "Save Edits", we need user to input the attributes related to the polyline drawn. In case user forget or skip this step (input attribute related), I intent to develop a function to open the "Create New Object" window (shows below) automatically, which means once user click "Save Edits" the "Create New Object" window will be open. The "Create New Object" window is script tool in toolbox. It should be auto triggered in the ArcGIS editor, my question is must I write my own editor, or is there a way to use the existing editor provided by ArcGIS? I am using ArcGIS 10.2. ![enter image description here](http://i.stack.imgur.com/oeebo.jpg)


title: python script in toolbox add a button on the window

text: I am getting stuck when I want to create a new python script in toolbox. The function I want to develop is to let user to input the user ID, then retrieve from the database the other values to allow user to do the updating. In this case, I want add a "Retrieve" button in the script to allow user to get their details. Please see the following picture of what I want to give. Is it possible? ![enter image description here](http://i.stack.imgur.com/d03kR.png)


title: FME - Raster Algebra and Multicriteria Analysis

text: I'm trying to use FME to do some raster algebra for some Multi-Criteria Analysis. At the moment, I have three rasters, all with the same extents and resolution: Slope, Hydro, and Infrastructure Mask. In the future, I will be adding additional layers (Vegetation, for example). The slope raster is a classified raster representing low to high slope. The categories are:   * 5: Less than 5 degrees (Preferrable)    * 2: 5-10 degrees    * 1: 10-15 degrees   * 0: Greater than 15 degrees (undesirable - No Go) The Hydro raster is a classified raster of buffer distances from waterways. The categories are:   * 5: Within 50m of a waterway   * 4: Between 50-100m of a waterway   * 3: Between 100-500m of a waterway   * 2: Between 500-1000m of a waterway   * 1: Greater than 1000m of a waterway The infrastructure mask has a value of 0 wherever there is some kind of infrastructure (Road, Pipeline, Powerline, Well, etc) and a value of 1 everywhere else. I would like to combine the Slope and Hydro (and eventually other layers) using FME. Ideally, I would like to add the values together, except, if the slope has a value of 0, the resulting raster should be zero. For the moment, I am using a RasterExpressionEvaluator Transformer to combine the layers, but I'm not sure of the right syntax to get what I want. ![enter image description here](http://i.stack.imgur.com/z84kg.png) This is the expression I'm using:               if(A[0]==0,0,($(Hydro_Weight)*B[0]+$(Slope_Weight)*A[0]))      A is the Slope Raster and B is the Hydro Raster. Hydro_Weight and Slope_Weight are parameters that the user enters before running the workspace. Essentially, what I'm doing is: if the slope is 0, then the result should be 0, otherwise, add the weighted values of the slope and hydro rasters. This seems like the correct way to go but all the cells in my output have a value of -128. What am I doing wrong? Is this the correct way to use the RasterExpressionEvaluator?


title: Attribute table in map composer is limited to 99 rows. How to change it?

text: I am using Qgis 2.2. I am about to create a cadastral index map, which are showing the parcels in specific gridcells, also there is an attribute table connected controlled by an atlas function. So by switching from gridcell to gridcell the attribute table automatically changes. Every gridcell has more than 99 parcels (Sometimes up to a few 100) BUT: The attribute table in the map composer is limited to 99 rows, so not all parcels can be shown in the attribute table. How can I increase the limitation of 99 rows? Thank you very much for your help.


title: Split line at a point with ArcGIS 10.1 Basic level license?

text: I have a point and polyline shapefile and want to split lines using points -- want to emulate the precise functionality of Split Line at Point tool (Data Management). I have ArcGIS 10.1 Basic level license and don't want to use third party tools -- want to achieve this by using Basic level license geoprocessing tools available within ArcGIS 10.1. Any suggestions would be very much appreciative.


title: Using Python in ArcGIS 10, create a line between the middle of the polylines which are intersected

text: I have got a river network as a feature class in a gdb. I need to calculate the topologic attributes via the connectivity graph of the river network. A graph-theoretic representation of river networks uses rivers as vertices and river junctions as links of a connectivity graph to understand the topology of river networks (takes rivers as nodes and river intersections as links of a connectivity graph). It will give information about centrality of lines and also simplify the network. In order to do this, I want to draw lines from the middle of each river branch to the middle of intersected higher order river. By this way, I can calculate the topological centrality of a river branch by using from-to node. When doing this manually, I draw a polyline from the middle of a river line to the middle of intersected higher order river line. Please help me to do this automatically by using Python in ArcGIS 10. If I do it manually it will take too long because of the drainage density.


title: How to resolve an Invalid Topology [Incomplete void poly] error from a geoprocessing service?

text: I am consuming a geoprocessing service from my Flex application. My geoprocessing service takes the input geometry, and intersects another layer. It then runs a Summery stats on the intersection output and returns the output. When I do a query or identify, and send the resulting geometry to the GP service, it runs fine. But when I Keep it in memory and do other things, which take a lot of time, and then send the geometry to the GP service, it errors out, with this error. The only difference I find in the later geometry, is that the rings are not in the same order as in the geometry I get from the server. What can I do to preserve the order of the rings?


title: Creating a Google Maps GroundOverlay with UTM NAD83 data

text: **Edit4:** The solution I posted from Jenner at VBforums was, obviously, in VB, and I used an online converter to port it to C#. Something was lost in translations, and that's why it was 10 miles off. I'm guessing I am just misunderstanding what the algorithms at the Proj.Net discussion boards are for and that's why they weren't doing what I wanted them to. I will attempt to close this question out in two days when I can mark my answer. * * * let me preface this by saying I am by no means a GIS expert. I need to overlay a png image over a Google map. Looking at the v3 api reference, the best way to do this is to use a GroundOverlay, which takes a lat,lng bounds. I am given one UTM NAD83 coordinate - the NW corner. I've found UTM to lat,long conversion algorithms online, but they don't even appear to be very accurate. The lat,lng bounds object in the Google API takes two lat,lng points - NE and SW. (Edit: I found a decent way to get the NE and SW points based on the NW point and the size of the image and the provided meters/pixel, so it appears my only problem is the accuracy of the conversion, if you would like to see the solution I used to get the NE,SW points, look here) (Edit2: I'm using .NET and just found Proj.Net; currently looking for good examples for using it). **Edit3, the current problem** : I am having issues converting the UTM to lat and long. A solution I found here a while ago gives me a point about 10 miles south of where I would expect it. A simple mathematical solution I found here, is giving me unexpected results.               double[] inverseMercator (double x, double y) {          double lon = (x / 20037508.34) * 180;          double lat = (y / 20037508.34) * 180;               lat = 180/Math.PI * (2 * Math.Atan(Math.Exp(lat * Math.PI / 180)) - Math.PI / 2);          return new double[] {lon, lat};     }     double[] toPoint = inverseMercator (686029.702258, 3581213.621173);      I get the following result: long: 6.1627096689832594 lat: 30.602349476368449 I get similar results using the Proj.Net solution provided by D_Guidi in the same thread. Using an online converter, I was able to get something closer to what I am expecting: lat: 32.35238307052292 long: -91.0230710652583 Can anyone shed any light on what I am doing wrong? I'd prefer to stay away from third party software and just find an algorithm to help me accurately overlay the image with the data I'm given, an example of which I've included below:               2.627469     0.000000     0.000000     -2.627469     686029.702258     3581213.621173     UTM     NAD83     15      First row is meters/pixel Any help you may be able to provide is certainly appreciated.


title: How to add feature class from ESRI personal geodatabase (.mdb) in QGIS

text: How do you add a feature class in a ESRI personal gdb in qgis. I understand this can be done through the option add vector layer and select geodatabase. But what after that ? It doesn't allow me to browse through .


title: Creating nanoscale DEM with GDAL

text: A bit of a strange question perhaps, but let me give you a breif explanation of the background before my actual questions: Atomic force microscopy (AFM) is a method that, in short (and to my limited knowledge) lets researchers scan areas at the micro- and nanoscale. It works by "scanning" an area using a probe of sorts. More that is difficult for me to explain, as I have no real understanding of it. What I do know, and what triggered my curiosity was that the result is in fact a "grid" of "height" values (a matrix of say 512x512 values describing the height of the probe at that point). I then thought: Well, apart from the scale, this is in fact a digital elevation model! And, this means that if I can manage to create a DEM file as understood by GIS tools I could apply GIS analysis to it! As it is, my significat other works at a lab that has a AFM-machine, and is using it in one of her projects. I got some scan files from her, and have managed, using Python (struct and numpy), to parse these binary files and what I have now is a numpy array of size 512x512 filled with int16-values. What I'm planning on next, and what I need some help with is the "mapping to a proper DEM"-part. I have some knowledge about DEMS, but when it comes to actual generation of them I am quite new. What I'm thinking is that I have to georeference my data in some way, and for this I need a custom (planar) coordinate system. I envision my coordinate system would use micro- or nano-meters as units. Then it's just a matter of finding the size of the area scanned with the AFM (this I believe is somewhere in the binary file, assume this known). **update** : I also have several scans at different resolutions, but of the same area. For example i have this information about two scans: larger image:               Scan Size: 51443.5 nm     X Offset: 0 nm     Y Offset: 0 nm      smaller (detail) image:               Scan Size: 5907.44 nm     X Offset: 8776.47 nm     Y Offset: 1486.78 nm      What I', thinking is that my custom coordinate system should have an origin in 0,0 and for the larger image I cam assign pixel 0,0 the coordinate value of (0,0) and pixel 512,512 the coordinate value (51443.5, 51443.5) (Guess you get the picture for the other points needed). Then, the larger image would map pixel (0,0) to (8776.47, 1486.78) and (512,512) to (8776.47 + 5907.44, 1486.78 + 5907.44) **1st Question is then** : How do I create a proj4 def for such a coordinate system? I.e: how do I assign these "real world coordinates" to my custom coordinate system (or, if I follow whubers suggestion and using a local coordinate system and lying about units (i.e. treating my nanometers as kilometers) Then I have to transfer my numpy 2-dimensional Array to a Georeferenced DEM file format. I was thinking of using GDAL (or, rather the Python bindings). **2nd Question is then** : How do I create a georeferenced DEM from "arbitrary" data such as mine? Preferably in Python and using open source libraries. The rest should then be fairly easy going, just a matter of using the right analysis tools. The problem is, this task is driven by my own curiosity, so I'm not quite sure what I actually should do with a nanoscale DEM. This begs the **3rd question** : What to do with a nanoscale DEM? What kind of analysis can be done, what are appropriate tools for DEM analysis and finally: is it feasible to make a map with hillshade and contour lines from this data?:) I welcome all suggestions and pointers, but keep in mind that I am looking for free alternativies, as this is a strictly hobby-based project with no budget or funding (and I do not have access to any lisenced GIS-applications). In addition, I do know that Bruker, the company that sells these AFM machines, ships some software, but using that would not be any fun.


title: PostgreSQL, PostGIS and GeoServer publishing 1.7 million records

text: I have a polygon geodata in postgresql and I am publishing out to geoserver the 1.7 million records. How long is too long when processing? EDIT: Forgot to mention that I successfully published a database spatial view subset of this same data. Side question, What is postgresql really doing during this process??? @Devdatta Tengshe, I have a polygon dataset with about 1.7 million records and I am just trying to publish the data to geoserver. I published a spatial view from a small table and a join to the large polygon dataset. This layer seems fine in openlayers. I just published another spatial view from a different small table joined to the large polygon dataset and it published but will not finish loading in openlayers.