title: Geoserver+postgis+leaflet

text: I have been using the getfeatureinfo operation on layers published from shapefile and it works however the data which is being served from within postgis throws an internal server error. I am not sure why it is doing it. Has anyone else experienced similar issue. Can it be compatibility issue between versions of Postgis and geoserver. The layers display correctly in both cases, the issue is just with getfeatureinfo. Version postgres 9.3 postgis 2.1 Geoserver 2.5 Latest version of leaflet. Any help will be appreciated. EDITED: Request: (Works fine when data is published from a shapefile but fails when layer is from postgis. SRID is 4326 and has been checked in the query and it is correct) var URL="http://localhost:8080/geoserver/Victory/wms?INFO_FORMAT=application/json&REQUEST=GetFeatureInfo&EXCEPTIONS=application/vnd.ogc.se_xml&SERVICE=WMS&VERSION=1.1.1&WIDTH="+ WIDTH+"&HEIGHT="+HEIGHT+"&X="+X+"&Y="+Y+"&BBOX="+BBOX+"&LAYERS=Victory:topolinefinal&QUERY_LAYERS=Victory:topolinefinal"; var content="";               $.getJSON(URL,function(data){         alert(data);     });      FireBug: Internal error occurred null Details: org.geoserver.platform.ServiceException: Internal error occurred at org.geoserver.wms.GetFeatureInfo.run(GetFeatureInfo.java:49) at org.geoserver.wms.DefaultWebMapService.getFeatureInfo(DefaultWebMapService.java:346) at sun.reflect.GeneratedMethodAccessor411.invoke(Unknown Source) at sun.reflect.DelegatingMethodAccessorImpl.invoke(Unknown Source) at java.lang.reflect.Method.invoke(Unknown Source) at org.springframework.aop.support.AopUtils.invokeJoinpointUsingReflection(AopUtils.java:319) at org.springframework.aop.framework.ReflectiveMethodInvocation.invokeJoinpoint(ReflectiveMethodInvocation.java:183) at org.springframework.aop.framework.ReflectiveMethodInvocation.proceed(ReflectiveMethodInvocation.java:150) at org.geoserver.kml.WebMapServiceKmlInterceptor.invoke(WebMapServiceKmlInterceptor.java:34) at org.springframework.aop.framework.ReflectiveMethodInvocation.proceed(ReflectiveMethodInvocation.java:172) at org.geoserver.ows.util.RequestObjectLogger.invoke(RequestObjectLogger.java:54) at org.springframework.aop.framework.ReflectiveMethodInvocation.proceed(ReflectiveMethodInvocation.java:172) at org.springframework.aop.framework.JdkDynamicAopProxy.invoke(JdkDynamicAopProxy.java:202) at com.sun.proxy.$Proxy64.getFeatureInfo(Unknown Source) at sun.reflect.GeneratedMethodAccessor402.invoke(Unknown Source) at sun.reflect.DelegatingMethodAccessorImpl.invoke(Unknown Source) at java.lang.reflect.Method.invoke(Unknown Source) at org.geoserver.ows.Dispatcher.execute(Dispatcher.java:774) at org.geoserver.ows.Dispatcher.handleRequestInternal(Dispatcher.java:272) at org.springframework.web.servlet.mvc.AbstractController.handleRequest(AbstractController.java:153) at org.springframework.web.servlet.mvc.SimpleControllerHandlerAdapter.handle(SimpleControllerHandlerAdapter.java:48) at org.springframework.web.servlet.DispatcherServlet.doDispatch(DispatcherServlet.java:923) at org.springframework.web.servlet.DispatcherServlet.doService(DispatcherServlet.java:852) at org.springframework.web.servlet.FrameworkServlet.processRequest(FrameworkServlet.java:882) at org.springframework.web.servlet.FrameworkServlet.doGet(FrameworkServlet.java:778) at javax.servlet.http.HttpServlet.service(HttpServlet.java:620) at javax.servlet.http.HttpServlet.service(HttpServlet.java:727) at org.apache.catalina.core.ApplicationFilterChain.internalDoFilter(ApplicationFilterChain.java:303) at org.apache.catalina.core.ApplicationFilterChain.doFilter(ApplicationFilterChain.java:208) at org.apache.tomcat.websocket.server.WsFilter.doFilter(WsFilter.java:52) at org.apache.catalina.core.ApplicationFilterChain.internalDoFilter(ApplicationFilterChain.java:241) at org.apache.catalina.core.ApplicationFilterChain.doFilter(ApplicationFilterChain.java:208) at org.geoserver.filters.ThreadLocalsCleanupFilter.doFilter(ThreadLocalsCleanupFilter.java:27) at org.apache.catalina.core.ApplicationFilterChain.internalDoFilter(ApplicationFilterChain.java:241) at org.apache.catalina.core.ApplicationFilterChain.doFilter(ApplicationFilterChain.java:208) at org.geoserver.filters.SpringDelegatingFilter$Chain.doFilter(SpringDelegatingFilter.java:74) at org.geoserver.wms.animate.AnimatorFilter.doFilter(AnimatorFilter.java:70) at org.geoserver.filters.SpringDelegatingFilter$Chain.doFilter(SpringDelegatingFilter.java:70) at org.geoserver.filters.SpringDelegatingFilter.doFilter(SpringDelegatingFilter.java:45) at org.apache.catalina.core.ApplicationFilterChain.internalDoFilter(ApplicationFilterChain.java:241) at org.apache.catalina.core.ApplicationFilterChain.doFilter(ApplicationFilterChain.java:208) at org.geoserver.platform.AdvancedDispatchFilter.doFilter(AdvancedDispatchFilter.java:49) at org.apache.catalina.core.ApplicationFilterChain.internalDoFilter(ApplicationFilterChain.java:241) at org.apache.catalina.core.ApplicationFilterChain.doFilter(ApplicationFilterChain.java:208) at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:311) at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:68) at org.springframework.security.web.access.intercept.FilterSecurityInterceptor.invoke(FilterSecurityInterceptor.java:116) at org.springframework.security.web.access.intercept.FilterSecurityInterceptor.doFilter(FilterSecurityInterceptor.java:83) at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:72) at org.geoserver.security.filter.GeoServerCompositeFilter.doFilter(GeoServerCompositeFilter.java:91) at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323) at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:68) at org.springframework.security.web.access.ExceptionTranslationFilter.doFilter(ExceptionTranslationFilter.java:113) at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:72) at org.geoserver.security.filter.GeoServerCompositeFilter.doFilter(GeoServerCompositeFilter.java:91) at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323) at org.geoserver.security.filter.GeoServerAnonymousAuthenticationFilter.doFilter(GeoServerAnonymousAuthenticationFilter.java:53) at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323) at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:68) at org.springframework.security.web.authentication.www.BasicAuthenticationFilter.doFilter(BasicAuthenticationFilter.java:150) at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:72) at org.geoserver.security.filter.GeoServerCompositeFilter.doFilter(GeoServerCompositeFilter.java:91) at org.geoserver.security.filter.GeoServerBasicAuthenticationFilter.doFilter(GeoServerBasicAuthenticationFilter.java:82) at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323) at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:68) at org.springframework.security.web.context.SecurityContextPersistenceFilter.doFilter(SecurityContextPersistenceFilter.java:87) at org.geoserver.security.filter.GeoServerSecurityContextPersistenceFilter$1.doFilter(GeoServerSecurityContextPersistenceFilter.java:52) at org.geoserver.security.filter.GeoServerCompositeFilter$NestedFilterChain.doFilter(GeoServerCompositeFilter.java:72) at org.geoserver.security.filter.GeoServerCompositeFilter.doFilter(GeoServerCompositeFilter.java:91) at org.springframework.security.web.FilterChainProxy$VirtualFilterChain.doFilter(FilterChainProxy.java:323) at org.springframework.security.web.FilterChainProxy.doFilter(FilterChainProxy.java:173) at org.geoserver.security.GeoServerSecurityFilterChainProxy.doFilter(GeoServerSecurityFilterChainProxy.java:134) at org.springframework.web.filter.DelegatingFilterProxy.invokeDelegate(DelegatingFilterProxy.java:346) at org.springframework.web.filter.DelegatingFilterProxy.doFilter(DelegatingFilterProxy.java:259) at org.apache.catalina.core.ApplicationFilterChain.internalDoFilter(ApplicationFilterChain.java:241) at org.apache.catalina.core.ApplicationFilterChain.doFilter(ApplicationFilterChain.java:208) at org.geoserver.filters.LoggingFilter.doFilter(LoggingFilter.java:75) at org.apache.catalina.core.ApplicationFilterChain.internalDoFilter(ApplicationFilterChain.java:241) at org.apache.catalina.core.ApplicationFilterChain.doFilter(ApplicationFilterChain.java:208) at org.geoserver.filters.GZIPFilter.doFilter(GZIPFilter.java:42) at org.apache.catalina.core.ApplicationFilterChain.internalDoFilter(ApplicationFilterChain.java:241) at org.apache.catalina.core.ApplicationFilterChain.doFilter(ApplicationFilterChain.java:208) at org.geoserver.filters.SessionDebugFilter.doFilter(SessionDebugFilter.java:47) at org.apache.catalina.core.ApplicationFilterChain.internalDoFilter(ApplicationFilterChain.java:241) at org.apache.catalina.core.ApplicationFilterChain.doFilter(ApplicationFilterChain.java:208) at org.geoserver.filters.FlushSafeFilter.doFilter(FlushSafeFilter.java:43) at org.apache.catalina.core.ApplicationFilterChain.internalDoFilter(ApplicationFilterChain.java:241) at org.apache.catalina.core.ApplicationFilterChain.doFilter(ApplicationFilterChain.java:208) at org.vfny.geoserver.filters.SetCharacterEncodingFilter.doFilter(SetCharacterEncodingFilter.java:109) at org.apache.catalina.core.ApplicationFilterChain.internalDoFilter(ApplicationFilterChain.java:241) at org.apache.catalina.core.ApplicationFilterChain.doFilter(ApplicationFilterChain.java:208) at org.apache.catalina.core.StandardWrapperValve.invoke(StandardWrapperValve.java:220) at org.apache.catalina.core.StandardContextValve.invoke(StandardContextValve.java:122) at org.apache.catalina.authenticator.AuthenticatorBase.invoke(AuthenticatorBase.java:501) at org.apache.catalina.core.StandardHostValve.invoke(StandardHostValve.java:171) at org.apache.catalina.valves.ErrorReportValve.invoke(ErrorReportValve.java:102) at org.apache.catalina.valves.AccessLogValve.invoke(AccessLogValve.java:950) at org.apache.catalina.core.StandardEngineValve.invoke(StandardEngineValve.java:116) at org.apache.catalina.connector.CoyoteAdapter.service(CoyoteAdapter.java:408) at org.apache.coyote.http11.AbstractHttp11Processor.process(AbstractHttp11Processor.java:1040) at org.apache.coyote.AbstractProtocol$AbstractConnectionHandler.process(AbstractProtocol.java:607)               at org.apache.tomcat.util.net.JIoEndpoint$SocketProcessor.run(JIoEndpoint.java:314)     at java.util.concurrent.ThreadPoolExecutor.runWorker(Unknown Source)     at java.util.concurrent.ThreadPoolExecutor$Worker.run(Unknown Source)     at org.apache.tomcat.util.threads.TaskThread$WrappingRunnable.run(TaskThread.java:61)     at java.lang.Thread.run(Unknown Source)      Caused by: java.lang.NullPointerException Experimented a bit more with projections: For shapefile based WMS: '..../wms?INFO_FORMAT=application/json&REQUEST=GetFeatureInfo&EXCEPTIONS=application/vnd.ogc.se_xml&SERVICE=WMS&VERSION=1.1.1&WIDTH=1958&HEIGHT=626&X=613&Y=337&BBOX=0.37724465131759644,52.74812277455308,0.3824964165687561,52.74913913313306&LAYERS=Victory:TopolineFinal&QUERY_LAYERS=Victory:TopolineFinal&styles=&srs=EPSG%3A4326' The above doesnt return any features when srs is included 4326 or 3857, so removing srs, it return the features '.../wms?INFO_FORMAT=application/json&REQUEST=GetFeatureInfo&EXCEPTIONS=application/vnd.ogc.se_xml&SERVICE=WMS&VERSION=1.1.1&WIDTH=1958&HEIGHT=626&X=613&Y=337&BBOX=0.37724465131759644,52.74812277455308,0.3824964165687561,52.74913913313306&LAYERS=Victory:TopolineFinal&QUERY_LAYERS=Victory:TopolineFinal' The Layer has native SRS of 4326 and declared of 3857. Reproject native to declared selected as option. Can anyone explain why the url with srs doesnt return any feature when srs is included?


title: What is the best raster format for large NoData areas?

text: Is there a good raster format for dealing with a large number of NoData cells? I do a lot of flood modeling and mapping, and recently I've been using rasters to show predicted flood depth and/or elevation for large areas, often for an entire county. I've just been using the default ESRI binary grid format, but I'm wondering if there might be something more appropriate. These rasters use a huge amount of space (especially at cell sizes of 5-10 ft), and the vast majority of the cells are NoData, so they're just going to waste.


title: Getting an error while trying to export a table

text: I'm trying to export a table to one of my folders but everytime i try to do it, i get this error:               An error occured trying to save the object.       What could cause this error? There are no umlauts in the path of the workspace.


title: OS X Python vs QGIS Python

text: I'm wondering why after I updated to Python 2.7.6 on my OS X 10.9 Mac, QGIS Python console says that it is still running the old 2.7.5 version. This is not affecting my work but I'm curious about why this happens. I though QGIS used the Python that comes with OS X but it seems it doesn't. Do any of you know something about it? I'm using KyngChaos QGIS 2.0.1 installers. Thanks!


title: Increasing OpenLayers Max Zoom Level

text: I develop HTML pages by using OpenLayers, which pointing to a local GeoServer WMS service. On displaying a layer, it stops at its available max zoom level, but it is still too small to manipulate (a polygon of 10 meter square). Is there a way to increase the zoom level? Where should I do the workaround -- on the OpenLayers or the layer setting in GeoServer? This is my OpenLayers code:                   var xwms = new OpenLayers.Layer.WMS( "Central Java",         "http://localhost:8080/geoserver/myneighborhood/wms/",          {             layers: 'basemap',             styles: 'tlx_map_land',             bgcolor: '0xBED0DE'         } );         var map = new OpenLayers.Map({            div: "map",            layers: [xwms],            controls: [                new OpenLayers.Control.Navigation({                    dragPanOptions: {                        enableKinetic: true                    }                }),                new OpenLayers.Control.PanZoom(),                new OpenLayers.Control.Attribution()            ],            center: [110.248191,-7.477177],            zoom: 13         });      Anyway, thanks for your attention.                        ============================ solved ===============================      Finally, I decided to utilise Klokantech's tool to generate custom bounding box(es), set it into both layer publishing in GeoServer and OpenLayers map options. To everyone who participate, thanks for your time and attention.


title: How can I create my own minifed openlayers js file?

text: I want to create my own minified openlayers js file with only the relevant classes my app used (including some custom class). Is there a way to know by using the full openlayers.js file which are the ones I need to include in my custom build (instead doing trial and error and start deleting files and see which breaks my app) ?


title: Extra row and column in a raster after GDAL grid import

text: I'm exporting an ESRI grid to TIFF format and after that I'm reading this exported raster file into SAGA GIS. However when Saga reads this file in there appears to be an additional row and column in the same file and thus changed raster extent. Also cellsize has got some pseudo decimals (2 m -> 2.0000000000000022 m). ArcMap (10.1) reports 24000 and 21000 while Saga (2.1.0) reports 24001 and 21001 columns and rows for the same file. Coordinate system is EUREF FIN TM35. Update: Saga uses GDAL import grid -module as a default to reading rasters in. If one selects this option from module tree (not by load -> file) there are two options to select:                 select from multiple bands [ ]        transformation [ ]       Under transformation there are different interpolation options, which seem to cause the problem. If none is selected, Saga reads my raster in with correct dimensions, but changes cellsize to 0.5. Any ideas on this one?


title: qgis polygon to coordinates

text: I have the ordnance survey postcode polygons, i can get the centre point on the postcode (Using polygon centroids) and all the corner coordinates (using extract nodes) Does anyone know how to export the coordinates of each postcode boundary to csv so it can be drawn on a google map?


title: Metadata in Dimap Generic 1.1

text: Does anyone know anything about the Dimap Generic structure for metadata? I have checked in http://www.spotimage.fr/dimap/spec/dictionary/dictionary.htm but I can't find anything about how to implement it on other data. I am supposed to deliver a DSM in GeoTiff with XML metadata in Dimap structure. What does that really mean? //Jenni


title: Does anyone know what a geoid10 and its numbers mean?

text: I am doing home for my GIS class where we had to download census data from our county and to define what geoid10 and its numbers mean. I cannot find the definition anywhere.


title: How to create a GeoJSON that works with D3?

text: I'm simply trying to convert a .shp file to a geoJSON format using:               ogr2ogr -f geoJSON output.json input.shp      After executing the command nothing seems to be wrong. Here's an excerpt from the output.json                   {         "type": "FeatureCollection",              "features": [         { "type": "Feature",          "properties": { "ID_0": 86, "ISO": "DEU", "NAME_0": "Germany", "ID_1": 1, "NAME_1": "Baden-Württemberg", "NL_NAME_1": null, "VARNAME_1": null, "TYPE_1": "Land", "ENGTYPE_1": "State" },          "geometry": { "type": "MultiPolygon", "coordinates": [ [ [ [ 8.708400940398242, 47.715581894910606 ], [ 8.713716147005524, 47.701734382960055 ],      ...      But when i try to use de JSON file in d3 (http://d3js.org/) to draw SVG polygons the output is just wrong. Since the shp files are displayed correctly in QGIS I think there has to be something wrong with the way i use `ogr2ogr`. The SVG I get is not completely wrong but there seems to be a detail I can't find. It seems like it has been turned upside down and somehow distorted into two separated parts. Here's the javaScript i used to generate the svg:               //dimensions     var w = 2000;     var h = 2000;          var svg = d3.select("#chart").append("svg")         .attr("width", w)         .attr("height", h);              d3.json(         "http://localhost:8888/data/data.json",         function (json) {              var path = d3.geo.path();              svg.append("g")             .attr("class", "black")             .selectAll("path")             .data(json.features)             .enter()             .append("path")             .attr("d", path);      Does anyone have an idea what went wrong here? I also tried to convert the shp-file using Qgis and myGeodata (http://converter.mygeodata.eu/vector). But neither of them work the way they should. I'm very new to this whole cartography stuff. So I would be very happy to get some advice. Thanks very much!


title: getResultImageLayer is not working

text: I am trying to obtain the shapefile output that I am receiving from the gpTool, which is hosted on the ArcServer, and placing it on the Map using,               gp.getResultImageLayer(jobInfo.jobId,"final_shp", imageParams, function(gpLayer) {                 map.addLayer(gpLayer);               });      where `gp` is,               gp = new esri.tasks.Geoprocessor("http://domain/ArcGIS/rest/services/webGIS/Shock_Models/GPServer/Income_Log");      This is how it has been explained in the resource but not sure why it is not working. The firebug gives me, > NetworkError: 500 Internal Server Error - > http://webgis.arizona.edu/ArcGIS/rest/services/webGIS/Shock_Models/GPServer/Income_Log/jobs/j10905b1a50f141e0b62f1ab8ef230f51/results/final_shp?dpi=96&transparent=true&format=png8&imageSR=102100&f=image&bbox=%7B%22xmin%22%3A-14519891.438356264%2C%22ymin%22%3A637618.0139790997%2C%22xmax%22%3A-6692739.741956295%2C%22ymax%22%3A6507981.786279075%2C%22spatialReference%22%3A%7B%22wkid%22%3A102100%7D%7D&bboxSR=102100&size=800%2C600 And when I goto that particular link, I get this error, > Server Error - Object reference not set to an instance of an object. Any idea how I can correct it? **UPDATE** I was told to use fiddler and see the error details that is occurring over the network and this is the output that I got, > SESSION STATE: Done. Response Entity Size: 849 bytes. > > == FLAGS ================== BitFlags: [ClientPipeReused, ServerPipeReused] > 0x18 X-CLIENTPORT: 2010 X-RESPONSEBODYTRANSFERLENGTH: 849 X-EGRESSPORT: 2023 > X-HOSTIP: 128.196.53.161 X-PROCESSINFO: firefox:2248 X-CLIENTIP: 127.0.0.1 > X-SERVERSOCKET: REUSE ServerPipe#2 > > == TIMING INFO ============ ClientConnected: 15:53:51.383 > ClientBeginRequest: 15:53:51.494 GotRequestHeaders: 15:53:51.494 > ClientDoneRequest: 15:53:51.494 Determine Gateway: 0ms DNS Lookup: 0ms > TCP/IP Connect: 0ms HTTPS Handshake: 0ms ServerConnected: 15:52:45.077 > FiddlerBeginRequest: 15:53:51.495 ServerGotRequest: 15:53:51.495 > ServerBeginResponse: 15:53:51.679 GotResponseHeaders: 15:53:51.679 > ServerDoneResponse: 15:53:51.679 ClientBeginResponse: 15:53:51.679 > ClientDoneResponse: 15:53:51.679 > > Overall Elapsed: 00:00:00.1850106 > > The response was buffered before delivery to the client. > > == WININET CACHE INFO ============ This URL is not present in the WinINET > cache. [Code: 2] * Note: Data above shows WinINET's current cache state, not > the state at the time of the request. * Note: Data above shows WinINET's > Medium Integrity (non-Protected Mode) cache only. But I am still confused as to what the error is? This is the application.


title: Ogr2ogr: Unable to open datasource problem

text: I downloaded the Precints shapefile from the City of Boston's data hub, and I'm trying import it into SQLServer (I'm running the Express verson). I've already created the database in the SQLServer (but not the tables), so now I'm using OGR to import the data. I keep getting this error:"FAILURE: Unable to open datasource 'precints.shp' with the following drivers." Here's what I've written in the OSGeo4W command prompt: ogr2ogr –f “MSSQLSpatial” “MSSQL:server=Comp1\SQLExpress;database=Precints- ogr2ogr; trusted_connection=yes” “Precincts.shp” I've tried downloading the data again, but no luck. I tried getting rid of the XML file that has the same name as the .shp file, but it still didn't work. Any idea on how to fix this? Thanks for your help!


title: Question about licence

text: I would like to ask about the licence. This software is presented as "Open Source Geographic Information System". But it is free also for company? Can I use it for creation of commercial maps in my business? It seems to be good alternative to ArcWIE. Thank you Vasek


title: MCD43A3 MODIS products

text: I have MCD43A3 product and I am using that to get the albedo values(Black sky and White Sky). I have some in-situ snow depth data, corresponding to which I want to get these albedo values. For most of the in-situ data I'm getting a value of 32767, which is a fill value for the product. For some data I'm getting, a value like 908,889,976 which is multiplied by scale factor of 0.0010 to get the albedo values like .908,.889,.976. My question is what should I do about the 32767 values which I am getting for most part of my data? Does getting this type of value mean that the image has noise or error? I have little experience with MCD products, any kind of help would be appreciated. Thanks, GBh


title: How to present three spatial variables in one map?

text: I want to assign 3 category of variables in a single map. In other words, I have to produce a cartographic visualization of the damage including a legend listing:   * damaged surface area in the residential areas (in km2),    * number of living units and number of people affected, and    * damage in terms of euroes. I just can present two of them by "graduated colors" and "graduated symbols". I cannot add third variable, in a meaningful way. How I can present all the 3 variable in a same snapshot simultaneously?


title: ArcGIS.com and related tables

text: I haven't found anything on the ArcMAP forums. Is there any way to represent related tables on arcGIS.com web maps? I have a geodatabase that has related tables, but so far I am only able to show the feature classes and not the related tables.


title: How to cluster 2 datasets in Google Maps API V3?

text: I'd like to produce a map that shows thousands of Markers.   I'm using Markersclusterer to show all these markers.   Yesterday, I received another huge data set (completely different) and I'd like to see its markers with the ones that belong to the previous data set (using MarkerClusterer.   So, I can have different icons for those two data sets.   Can you help me do this?   thanks


title: Is there a way to programatically trigger popups using the operationallayer config?

text: I have a widget that queries a group of features, and I'd like to be able to trigger a popup that looks and behaves the same as the one that gets generated when the user clicks on the map. Is there a call to do this in the flexviewer itself?


title: Clipping .mdb in place and schema locks

text: I have an .mdb with a feature class that I want to clip in place - that is, I want to clip it with another polygon, and then store the results back in the original .mdb. I am doing this with arcpy. My current strategy is to clip it to an in-memory feature class, then copy that feature class back over the original.               in_mem = "in_memory/temp_clip"     arcpy.Clip_analysis(src_fc, clip_fc, in_mem)     arcpy.env.overwriteOutput = True     arcpy.CopyFeatures_management(in_mem, src_fc)      However, the CopyFeatures call fails with `ERROR 000601: Cannot delete /path/to/src_fc. May be locked by another application.`. As best as I can tell, there is no lock on `src_fc`. There is no .ldb with it, there's nothing in `GDB_DatabaseLocks`, no program has it open (according to Process Explorer), I've tried restarting my computer, and running compact and repair both from MS Access and `arcpy.Compact_management()`. I can edit the database fine in MS Access, it is only Arc that claims there is a lock on it. Therefore, I suspect that there's some obscure lock arc is holding that it doesn't want to be released, or there's a bug somewhere. So, my question is twofold: First, is there a better way to do what I'm trying to do (replace a feature class with the result of clipping it), and second, if not, how can I resolve this locking issue?


title: Using ArcGIS to interpolate X, Y, Coordinates and ordinary kriging of point data

text: I have been using Geosoft's Target to interpolate X, Y coordinates then use them to create grids using ordinary kriging. I am now working with ArcGIS which I am new to and cant seem to be able to interpolate my coordinates and create the grids. Any ideas?


title: Need help in how to map ZIP codes data by year

text: I'd like to build a graphic showing where city employees live, based on their date of hire. I have a list showing the hire date and the ZIP code for each employee. There are about 4,000 employees. I can't figure out how to do that, despite it being very easy to explain!


title: Is there a way to automatically snap loose points of lines to close by vertices with Quantum GIS?

text: I have a vector layer of lines that were drawn with the snapping option off. I need to have them snapped but, there are hundreds of lines. Is there any solution to do this automatically or semi-automatically? Thanks.


title: Problems with transformations in ArcMap

text: I have files projected in British National Grid.   After I have transformed them into WGS_1984.   I noticed that if I were to transform them back to BNG the original would not match the new.   Is there a way to make sure this does not happen especially when moving from BNG to WGS 1984 multiple times?


title: How can I use Openflights database in QGIS?

text: I would like to have a shape file with the airports position and the routes of the planes I've found Openfilghts but I can't use it in QGIS. Seems that the format of the DB `.dat` is not supported by QGIS. How can I overcome this problem?


title: Finding routes in which more than one car passed (gps data)

text: I'm using ArcGIS 10.1. I have point GPS data of routes of several cars over a few days, and I want to find only the routes where more than one or two different cars passed through (Each GPS points has the car id in it). How can I find a thing like that without going through the data manually ? EDIT: The time interval (and the distance) between points differs - it moves from minutes to hours and from a few meters to tens of kilometers apart. I'm comparing routes right now by viewing the points by the eye (which isn't accurate or useful in any way). I thought about creating routes using network analyst, but as I have a lot of points I kind of hoped to fan out the less popular ways first, and calculating routes for only the popular ones...


title: Adjusting Z Units on a DEM

text: I am learning how to use ArcHydro and attempting to change the Z Units because I get a warning message saying the Spatial Reference Doesn't Match. The DEM was in Meters, and I adjusted the Z Unit to 100. My DEM is in a GCS and when I project it to the same PCS my vector data is in, it puts them in two very different places, even though they are in NC State Plane 3200 FIPS (feet). If someone can give me a recommendation to fix this, I would greatly appreciate it.


title: OpenLayers: Combine multiple OGC filters

text: I have two _OGC_ filters client side:               <ogc:Filter xmlns:ogc="http://www.opengis.net/ogc">         <ogc:PropertyIsEqualTo matchCase="true">             <ogc:PropertyName>code</ogc:PropertyName>             <ogc:Literal>secret</ogc:Literal>         </ogc:PropertyIsEqualTo>     </ogc:Filter>             <ogc:Filter xmlns:ogc="http://www.opengis.net/ogc">         <ogc:Intersects>             <ogc:PropertyName>geom</ogc:PropertyName>             <gml:Box xmlns:gml="http://www.opengis.net/gml">                 <gml:coordinates decimal="." cs="," ts=" ">1298543.0367787,7107400.531374 1299312.4070614,7203464.9016567</gml:coordinates>             </gml:Box>         </ogc:Intersects>     </ogc:Filter>      I have both filters represented as strings and I want to combine them into one filter. I can see in the documentation that OpenLayers can create a _Filter Object_ from a string, but it doesn't seem to work for me. **Am I doing something wrong and/or how should I do it?**


title: How to use for loop to add and calculate new field?

text: I am trying to use for loop to add and calculate a newfield for 2 shapefiles and a different new field for 2 other shapefiles. When I tried to do these 2 loops as one, it did not work so I separated them into 2 loops and it worked. I am using python 2.7 stand alone (arcmap 10.2). more details: So I want to make a new field and add 'CA' to shapefiles containing _CA.shp (*_CA.shp). I also want to make a new field and add 'AZ' to shapefiles containing _AZ.shp (*_CA.shp). When I created a loop (the first code), I got my roads_CA and railroads_CA to work but roads_AZ and railroads_AZ did not create a new field containing 'AZ'. I did not get an error which means the loop ran through but I guess the elif came out false so it did not create anything. here is my code:                   CAshp = arcpy.ListFiles("*_CA.shp")              for CA in CAshp:             if (CAshp == 'roads_CA.shp', CAshp == 'railroads_CA.shp'):             newField = 'StateAbbre'             arcpy.AddField_management(CA, newField, 'TEXT')             arcpy.CalculateField_management(CA, newField, "'CA'", "PYTHON_9.3")         elif (CAshp != 'roads_AZ.shp', CAshp != 'railroads_AZ.shp'):             AZshp = CAshp.extend(arcpy.ListFiles("*_AZ.shp"))             newField = 'StateAbbre'             arcpy.AddField_management(AZshp, newField, 'TEXT')             arcpy.CalculateField_management(AZshp, newField, "'AZ'", "PYTHON_9.3")


title: What is the most efficient way to iterate through folders and geodatabases?

text: I am looking at cataloging some data, and I have code that works, it just seems to take a lot of time to run through all of the items. The way this data is structured is a bit of a hindrance, particularly with the amount of directories I am looking to analyze. The data structure is similar to this:               Working Folder      Property Folder (All with a specific label ogl###)           Geodatabase(all the same name in each different folder)                Feature Dataset                     3 Shapefiles (boundary, parcel, corner)                60 Additional Shapefiles of varying types           Misc. Folders, Tiffs, and geodatabases      Misc. Folders      This may not have been the best structure to set up our data library, but that is above what I am trying to do, and I have no influence over that. My goal is to do some quality control on this data, making sure that each Property Folder, Geodatabase, and Feature Class is structured correctly. To do this I need to check each Property Folder\Geodatabase individually, write to excel the names of the corresponding Feature Classes in each database, and then within each feature class the correct fields. I have written working code to check the existence of the correct feature classes and export to excel the results of that check. To me, however, it seems like it is taking longer than it should. Here is an example of my code to accomplish that (minus the extraneous excel stuff):               import arcpy           arcpy.env.workspace = r'E:\Working'     working_folder_list=arcpy.ListWorkspaces('ogl*','Folder')           ##Generate list of Property                                                                         #Folders      ######This block is cataloging a blank reference folder######          arcpy.env.workspace = r'E:\Working\Reference GDB\OGL.gdb\Admin'     ##looking into the reference      fc_ref_ad=arcpy.ListFeatureClasses('*')                             ##folder Feature Dataset     fc_ref_ad_dic={}     for fc in fc_ref_ad:         fc_ref_ad_dic[fc]=[f.name for f in arcpy.ListFields(fc)]          arcpy.env.workspace = r'E:\Working\Blank DB 4-05-11\SGL.gdb'        ##looking into the reference     fc_ref=arcpy.ListFeatureClasses('*')                                ##geodatabse     fc_ref_dic={}     for fc in fc_ref:         fc_ref_dic[fc]=[f.name for f in arcpy.ListFields(fc)]           ##creating a dictionary of fields                                                                         ##per FC     for fc in fc_ref_ad:                            fc_ref.append(fc)                                               ##Merging both lists and dicts          for fc in fc_ref_ad_dic:         fc_ref_dic[fc]=fc_ref_ad_dic[fc]     fc_ref.sort()     #############################################################               ######This block looops through each Property Folder#########          for folder in working_folder_list:         arcpy.env.workspace = folder + "\\ogl_data0.gdb\\Admin"         ##Listing FC in Feature Dataset         admin_ds=arcpy.ListFeatureClasses('*')          arcpy.env.workspace = folder + "\\sgl_data0.gdb"                    ##Listing FC in Geodatabase     current_gdb=arcpy.ListFeatureClasses('*')     for fc in admin_ds:         current_gdb.append(fc)                                          ##Merge the 2 lists          if len(current_gdb) > len(fc_ref):                                  ##Dealing with uneqaul len(lists)         for j in range(len(current_gdb)-len(fc_ref)):             fc_ref.append('')          ##block of code to write lists to corresponding cells in excel          #############################################################               ######This block goes through each FC and lists fields#######     ######Approximately 4000 directories(66 Property Folders and 60 FCs######          for folder in working_folder_list:          arcpy.env.workspace = folder + "\\sgl_data0.gdb\\Admin"     fc_cur_ad=arcpy.ListFeatureClasses('*')     fc_cur_ad_dic={}     for fc in fc_cur_ad:         fc_cur_ad_dic[fc]=[f.name for f in arcpy.ListFields(fc)]          arcpy.env.workspace = folder + "\\sgl_data0.gdb"     fc_cur=arcpy.ListFeatureClasses('*')     fc_cur_dic={}     for fc in fc_cur:         fc_cur_dic[fc]=[f.name for f in arcpy.ListFields(fc)]          for fc in fc_cur_ad_dic:         fc_cur_dic[fc]=fc_cur_ad_dic[fc]          ##########################################      I know it is a lot to look through, and I may not be the cleanest coder. Truth be told I am still learning this stuff, and most of it I am teaching to my self on the fly. I am just wondering if there is a more efficient way I could be doing this apart from restructuring the data itself. I am aware of the arcpy.da.Walk, but I am unfamiliar with it, and not sure if it would do me any good. Perhaps I am wrong. Either way, I am hoping someone may be able to give me some tips. I want to thank you in advance for looking. Edit: After trying to learn arcpy.da.Walk, I created a sample script that averages the times it takes for each method. I am a bit shocked that after 20 iterations, that the Walk() function took an average of 1 minute longer per run through my data. Here is that script:               import arcpy     import os     import time          worktimes=[]     walktimes=[]          counter=0          while counter <20:     t_start=time.clock()          arcpy.env.workspace = r'E:\Working'     working_folder_list=arcpy.ListWorkspaces('ogl*','Folder')          for folder in working_folder_list:         arcpy.env.workspace = folder + "\\ogl_data0.gdb\\Admin"         admin_ds=arcpy.ListFeatureClasses('*')         for fc in admin_ds:             print fc             fields = arcpy.ListFields(fc)             for f in fields:                 print f.name              arcpy.env.workspace = folder + "\\ogl_data0.gdb"         current_gdb=arcpy.ListFeatureClasses('*')         for fc in current_gdb:             print fc             fields = arcpy.ListFields(fc)             for f in fields:                 print f.name          workspacetime=time.clock()-t_start     worktimes.append(workspacetime)          t_start2=time.clock()     workspace = 'E:\Working'     for dirpath, dirnames, filenames in arcpy.da.Walk(workspace):         for filename in filenames:             fields=arcpy.ListFields(os.path.join(dirpath,filename))             print filename             for f in fields:                 print f.name          walktime = time.clock() - t_start2     walktimes.append(walktime)     counter+=1          worksum = 0     walksum = 0          for time in worktimes:     worksum += time          for time in walktimes:     walksum += time          workavg = worksum/len(worktimes)     walkavg = walksum/len(walktimes)          print 'Average workspace time: '+str(workavg)     print 'Average walk time: '+str(walkavg)      Is there something I'm missing, or is Walk() not as efficient as I assumed. The average times were 104 seconds for changing workspaces, 167 using the Walk() function.


title: From a shapefile, copy specific rows (defined by a separate table) to a new shapefile

text: I have a polygon shapefile of origin zones. I also have a polygon shapefile of destination zones. I also have a table that lists origin-destination pairs. I would like to create a multi-part shapefile where each row represents an origin-destination pair. So when you select a single row, two separate polygons would be highlighted (one polygon representing the origin zone, the other representing the destination zone). I don't know how to go about creating the multi-part shapefile. Does anyone have any suggestions?


title: Geoserver Status Screen

text: Can anyone tell me what all the information on this Geoserver status this means and if there are any issues with it? I cannot find anything in the documentation that explains what all this is actually telling us. We have figured out the reload button, but not sure if it is required for EVERY configuration change or not, but we generally press it after changing anything, including styles, GeoWebCache changes, layer changes, etc. Thanks in advance Mark ![enter image description here](http://i.stack.imgur.com/m5f20.jpg)


title: Drawing a map of a country with specific boundaries/zones using postcodes latitude and longitude

text: I'm quite new in geocoding - it has been fun! My project is to draw Australia with specific boundaries/zones using R program. I was able to successfully upload a shapefile of Australia in R program, which can be seen below. The file can be found in: http://www.abs.gov.au/AUSSTATS/abs@.nsf/DetailsPage/1259.0.30.001July%202011?OpenDocument I want to be able to edit the boundaries/zones in Australia. My question: 1) Can I edit the zones/boundaries in a shapefile? 2) I've got a list of postcodes with latitude and longitudes in Csv format. Is there a way of grouping postcodes to form zone/boundaries (polygon) using R program? If not, is there a freeware software to do this? I haven't had much lucky with QGIS & ArcGIS. Thank you for your help.


title: Why am I missing the the geometry column (geometry_columns)?

text: Can anyone tell me what could I be doing wrong? I have installed OpenGeo and started doing the exercises. I noticed that every time I create a database through the pgAdmin, the geometry column (geometry_columns) is not created. Did I missed a step?


title: Geoserver 2.2.5 and Integrated GeoWebcache caching Defaults not saved over restarts

text: I am experience an annoying issue with Geoserver 2.2.5 and the integrated GeoWebCache running as a service on Windows. The Caching Defaults Configuration is not retained over GeoServer Restarts or Server Reboots. This is particularly annoying in a production Environment as a server reboot or _GeoServer Service Auto Restart on failure_ ends up with GeoServer running with no cache at all. We have our Service set to restart the service on failure, then reboot the computer automatically on the third failure. The documentation points to a configuration file in the Geoserver data directory named _gwc-gs.xml_ that holds the configuration details. No matter what we do by editing this file, we cannot get the configuration to persist. Interestingly, I noticed that a file named _gwc-gs.xml.tmp_ was also in the same location. When I delete this .tmp file and then restart Geoserver, it appears again, always containing the default values So my assumption is that this _.tmp_ file is being created on Geoserver start up and it is being read by the GeoWebCache code, not the actual _.xml_ file. Question: Has anyone else experienced this or know of this behavior? Question: Does anyone know of a work around or know if this behavior has been corrected in a later version of Geoserver? Thanks in advance Mark


title: How can you conditionally halt an FME translation?

text: I have an FME workspace that takes a set of street centerlines and converts it into address points. I have been working my way across the state of Missouri with this translator one county at a time, since each county can take up to an hour to process. I select a county at the start of the translator as a choice parameter. My problem is that I keep having to double check to make sure of what counties I have processed so I do not process a county twice. When I have a file of several million address points, I do not want to be removing duplicates. What I would like to have FME do is take the county I have selected, and check against a shapefile of counties to see if the county has been processed yet. I have a flag attribute, "COMPLETE" that is set to 0 if the county has not been processed and 1 if it has. The county shapefile has nothing to do with any other readers or writers in the translator. Edit: To clarify, the counties are not being processed at all by the translator, only the centerlines. Since there are over 2 million centerlines, I do not want to read the centerlines at all (much less maintain an attribute on the centerlines) if the county has already been processed. Edit: Here is the solution. I moved up the st_county10 reader (the one in the image) to the top of the reader list. This executes before anything else. By using the feature holder, I avoid writing over the st_county10 shapefile if the translation is terminated. One last piece, since I had a python shutdown script, I had to add this line:               import pyfme     if FME_Status == True:      To the front of the shutdown script to stop execution of the shutdown script too. ![If COMPLETE=0, translation proceeds, otherwise it goes to the Terminator](http://i.stack.imgur.com/Ro3XG.gif)


title: How to create a circular sector in Mapinfo?

text: I want to build a tool that creates a cell sector from data like that below. The sector diagram should be like the second picture. The first picture shows a Mapinfo tab file with the required fields Antenna beam, Azimuth, Band, (the radius), Latitude, and Longitude. How could I achieve this? I had a look at a link in mapinfo-L but the solution was not conclusive! ![Different Sites](http://i.stack.imgur.com/qoh6D.png) ![Result Cell Sector](http://i.stack.imgur.com/LDukn.png)


title: GeoNetwork opensource - Geoserver/Google Maps integration

text: I'm using GeoNetwork 2.6.4 and I would like to change the default base map to Geoserver, adding a Google Maps Satellite Layer. Which is the best way? Any suggestions?


title: Basic .NET mapping components?

text: Which .NET mapping components are currently available? Focus is on desktop use and fairly basic mapping of SHP files (thematics etc.), not the more advanced GIS functionality (reprojections, clipping, point-in-polygon etc). Royalty free distribution is needed. I am aware of TatukGIS, AspMap, Map Suite and MapXtreme. I believe ESRI has nothing to offer in this area any longer, now that mapobjects is gone. Other options? Any stable open source offerings available?


title: How to change the appearance of a Feature.Vector that's already displayed?

text: I've created a bunch of Feature.Vectors from Point objects and assigned them a Style. Now I want to dynamically change the appearance of one of the features in real time; e.g., change it from red to blue. (I am NOT hovering over the feature.) I've tried changing its Style, then doing refresh() on the layer, or redraw(), etc, and either nothing happens or the feature disappears. Any suggestions? Thanks.


title: What's wrong with my label formatting?

text: I have some formatting issues with my labels. I am using a fairly straight forward expression (see below). It sort of works fine but **some formatting is being ignored, namely the ALL CAPS ACP tags**. All owner names are lower case with first letter capitalized. I want to display them as all caps and only SOME labels come out with owners capitalized. I am also using italics tags when a claim number value exists and that works fine. Another issue is that although I am **NOT spreading characters** in Maplex, **some labels but not all have Owner's names spread** , but strangely not the PIN numbers and claims numbers all of which are a part of one single string variable. There are other issues with Stacking Labels in Maplex which doesn't seem to stack properly when using Stacking Separator but let's ignore that for this question. Is there a way to fix this or do i need to export to annotations and fix manually?               Function FindLabel ([Owner], [MNDMF_Active_Mining_Claims_CLAIM_NUM], [PIN], [PIN_Surface], [PIN_Mining] )     FindLabel = ""          if not isnull([Owner]) then     if trim([Owner]) <> "" then     FindLabel = FindLabel & trim("<ACP>" & [Owner] & "</ACP>")     end if     end if          if not isnull([MNDMF_Active_Mining_Claims_CLAIM_NUM]) then     if trim([MNDMF_Active_Mining_Claims_CLAIM_NUM]) <> "" then     if len(FindLabel) > 0 then     FindLabel = FindLabel & vbnewline     end if     FindLabel = "<ITA>" & "Claim: " & [MNDMF_Active_Mining_Claims_CLAIM_NUM] & "</ITA>"      end if     end if            if not isnull([PIN]) then     if trim([PIN]) <> "" then     if len(FindLabel) > 0 then     FindLabel = FindLabel & vbnewline     end if     FindLabel = FindLabel & trim([PIN])     end if     end if          if not isnull([PIN_Surface]) then     if trim([PIN_Surface]) <> "" then     if len(FindLabel) > 0 then     FindLabel = FindLabel & vbnewline     end if     FindLabel = FindLabel & trim([PIN_Surface])& " S"     end if     end if          if not isnull([PIN_Mining]) then     if trim([PIN_Mining]) <> "" then     if len(FindLabel) > 0 then     FindLabel = FindLabel & vbnewline     end if     FindLabel = FindLabel & trim([PIN_Mining])& " M"     end if     end if          End Function


title: Log file for ArcGIS geoprocessing?

text: Whenever I run a geoprocess overnight, it seems my computer reboots. Does ArcGIS save some log file information? Where does it save it? For example, when I run a large dataset or image process I need to perform overnight, does it save the logfile that tells you what it performs from start to finish? Thanks. PS Do I make myself clear? If not, perhaps someone in here can help me correct what I am trying to say.


title: Query for Selecting only Selected Fc's from a GDB

text: I need to version n number of fc's which start with sde."ENT-QA". I need to select only those feature classes and run versioning for it. Can anyone help me in writing the query for selecting only fc's which start with sde."ENT-QA". Remeber ENT-QA is mentioned has quotes. Thanks in advance.               import arcpy     arcpy.env.workspace= r"xyz.sde"     fdList = arcpy.ListFeatureClasses("","","")          for fc in fdList:         print fc         arcpy.RegisterAsVersioned_management(fc, "NO_EDITS_TO_BASE")         print "fc registered"


title: QGIS Dev install issues

text: I need to use some of the new functionality in the Dev build and just can not install it on our windows 7, 64 bit machines. It works fine on my personal windows 8 laptop but I can't get it to work at work. I have tried unistalling / reinstalling many times and Lisboa is fine but Dev won't work. I tried the fix in Why is QGIS not launching? but no luck.               Problem signature:       Problem Event Name:   APPCRASH       Application Name: qgis.exe       Application Version:  0.0.0.0       Application Timestamp:    4fdce845       Fault Module Name:    QtCore4.dll       Fault Module Version: 4.7.1.0       Fault Module Timestamp:   4d530715       Exception Code:   c0000005       Exception Offset: 000ee27b       OS Version:   6.1.7601.2.1.0.256.4       Locale ID:    3081       Additional Information 1: 0a9e       Additional Information 2: 0a9e372d3b4ad19135b953a78882e789       Additional Information 3: 0a9e       Additional Information 4: 0a9e372d3b4ad19135b953a78882e789      This is on multiple machines with Win 7 -64 bit. I have been trying to install qgis-dev but I keep getting qgis_core.dll not found issues and I have looked online but can't find anything that works. You get the Dev screen (image) for a millisecond and then no output. The paths seem to load fine.               OSGEO4W home is C:\OSGeo4W          C:\>path C:\OSGeo4W\apps\grass\grass-6.4.3RC3\bin;C:\OSGeo4W\apps\msys\bin;C:\OS     Geo4W\bin;C:\Program Files (x86)\Common Files\Intergraph\Grid Analysis SDK\1.0\P     rogram;C:\Windows\system32;C:\Windows;C:\Windows\System32\Wbem;C:\Windows\System     32\WindowsPowerShell\v1.0\;c:\lastools;C:\apps\FME\;C:\Program Files (x86)\Quick     Time\QTSystem\;C:\OSGeo4W\apps\msys\bin;C:\OSGeo4W\apps\Python27\Scripts;C:\OSGe     o4W\apps\qgis-dev\bin;C:\OSGeo4W\apps\grass\grass-6.4.3RC3\lib          C:\>set QGIS_PREFIX_PATH=C:/OSGeo4W/apps/qgis-dev          C:\>rem Set VSI cache to be used as buffer, see #6448          C:\>set VSI_CACHE=TRUE          C:\>set VSI_CACHE_SIZE=1000000          C:\>start "QGIS" "C:\OSGeo4W"\bin\qgis-dev-bin.exe          C:\>pause     Press any key to continue . . .


title: How draw a Rectangle with INewRectangleFeedback?

text: I´m developing in ArcGIS Engine and my question is how to draw a rectangle with `INewRectangleFeedback` interface. Here is my code but it does not create the rectangle. If someone could help me it would greatly appreciate it.               Public Sub mousedownFBRectanglePoly(ByVal px As Long, ByVal py As Long)              feedbackStartPoint = New ESRI.ArcGIS.Geometry.Point         feedbackMovePoint = New ESRI.ArcGIS.Geometry.Point         feedbackScreenDisplay = Principal.m_mapControl.ActiveView.ScreenDisplay         feedbackFillSymbol.Color.CMYK = RGB(255, 0, 0)              'Initialize a new RectangleFB object         NewRectangleFeedback = New NewRectangleFeedback         NewRectangleFeedback.Display = feedbackScreenDisplay         NewRectangleFeedback.Symbol = feedbackFillSymbol 'QI              'pass the start point from the mouse position, transforming it to map coordinates.         feedbackStartPoint = feedbackScreenDisplay.DisplayTransformation.ToMapPoint(px, py) 'QI         NewRectangleFeedback.Start(feedbackStartPoint)     End Sub          Public Sub mousemoveFBRectanglePolygon(ByVal px As Long, ByVal py As Long)         '    'pass X and Y to feedbackMovePoint to transfer to NewEnvelopeFeedback         feedbackMovePoint = feedbackScreenDisplay.DisplayTransformation.ToMapPoint(px, py) 'QI         NewRectangleFeedback.SetPonint(feedbackMovePoint)     End Sub          Public Sub DoubleclickFBRectanglePolygon(ByVal px As Long, ByVal py As Long)         'feedbackStopPoint = feedbackScreenDisplay.DisplayTransformation.ToMapPoint(px, py)         'when mouse comes up, end the new envelope and pass it to feedbackEnvelope.         feedbackRectangle = NewRectangleFeedback.Stop(feedbackStopPoint) 'QI         'pass the new rectangle element, the geometry defined by our feedback object         feedbackElement = New RectangleElement         feedbackElement.Geometry = feedbackRectangle         'make sure the element is activated in the current view         feedbackElement.Activate(feedbackScreenDisplay)         'now add the newly created element to the ActiveView with default symbology.         Principal.m_mapControl.ActiveView.GraphicsContainer.AddElement(feedbackElement, 0)         'and refresh the view so we can see the changes.         Principal.m_mapControl.ActiveView.Refresh()     End Sub      And also I have the problem that when I call my methods in the mouse events it does not enter the method; it skips even if all the parameters we pass are good.


title: Query/find features in ArcPad, hide the rest in map?

text: Trying to help out another department that does quarterly assessments of field data. Their symbology is all set by the status of the assessment (the only way I could figure out how to hide something was to hide symbology) so that doesn't help much. Their checkouts contain all data for the year but they want to be able to filter/limit what is visible in the map based off the quarter (and potentially other data columns). Been trudging through documentation but can't really find anything. Is this possible?


title: Distance between GPS coordinates

text: Given a database of GPS points, a center GPS location, and a radius R (in meters), how can I find all points in the database that are within a distance of R of the center point? Will the distance iso-lines be an ellipse?


title: How do you merge small gaps into the larger surrounding polygon?

text: I have a layer which has small rectangular gaps, I need to merge the gaps into the larger polygon that surrounds it, has anyone got any idea of how to do this, I know i can fill in the gaps with a new feature then go to Editor> Merge, but i have hundreds of these. If theres a way of doing this via vba or python, it would be great. I'm using ArcGIS 9.3.1 (ArcEditor) Heres what it looks like ![gaps is layer](http://i.stack.imgur.com/Bx1VZ.jpg)


title: Label should appear inside the polygon only

text: I'm working on vectorLayers of OpenLayers and I'm labeling my layer based on some attribute, while the label appear at some certain zoom level. The problem I'm facing is that though label appears at certain zoom in, yet there are some polygons with enough small size that the labels get conjusted there. Well, what I want is that the labels should appear only inside the polygon, if the polygon size is small, the label may not appear. I searched for this context and I found OpenLayers function of "Strategy" and "clusters", but I couldn't find a way for my problem... Any sample code or relevant link...


title: How to reproject geodatabase (*.gdb) with datasets and featureclass using custom transformation?

text: I have a file geodatabase in arcMap which in in EVEREST_1830 datum and Modified UTM projection. Now I need to use OpenStreetMap or googlemaps imagery WMS layer as base layer and overlay the datasets. I know that I can Project each feature class in the Gdb to WGS 84 and do Overlay. But I will take time to project each featureclass. So the short-cut might be to bring the WMS imagery layer to local (.gdb) coordinate system or to bring the .gbd to WGS 84. Is it Possible? In Short:   Is it possible to project a whole geodatabase including all datasets and feature class (*.gdb) in ArcMap to a another projection system using a user defined custom transformation? Or is there a tool, python script etc...


title: is there significance difference between qgis version 1.8.0 & 2 to install plugin?

text: I am a Qgis learner on widows xp and I want to know whether I can install semiautomatic classification plugin onto Qgis 1.8.0 in the same way for Qgis 2. is there significance difference between the two version in accepting the plugin?


title: How to extract all lines that do not intersect with buffer?

text: I am working in ArcGIS Desktop 10.0 and/or 10.1. I have two road network feature classes in a file geodatabase. I have created a 5 meter buffer around the first one, so it is now a polygon feature. The other feature is a line feature. I want to intersect the line feature road network with the buffered one and generate an output feature class that includes only the road lines that fall **outside** of the buffer. When I run Intersect, the output is only the portion of the lines that fall inside the buffer. Is there a process to "clip off" the lines that fall inside the buffer polygons and just output a line feature containing all the lines outside the buffer? I am open to using other software if necessary. Thanks for your help!


title: alternative to mapClickEvent

text: I use gwt-openlayers lib in my GWT app. I have added a drawFeature control and now I want to show a radius of circle while I draw a feature. Now it looks like this:               vectorLayer.addVectorFeatureAddedListener(new VectorFeatureAddedListener() {             @Override             public void onFeatureAdded(FeatureAddedEvent eventObject) {                 final VectorFeature feature = eventObject.getVectorFeature();                 if (feature.getGeometry().getClassName().equals("OpenLayers.Geometry.Polygon") && !hasTemplate) {                     if (currentAreaFeature != null) {                         vectorLayer.removeFeature(currentAreaFeature);                     }                     if (templateAreaFeature != null) {                         removeTemplateArea();                     }                     currentAreaFeature = feature;                     Geometry geometry = currentAreaFeature.getGeometry().clone();                     geometry.transform(new Projection(map.getProjection()), DEFAULT_PROJECTION);                     geometry.calculateBounds();                     Bounds bounds = geometry.getBounds();                     currentCenter = bounds.getCenterLonLat();                     currentRadius = calculateDistance(currentCenter, new LonLat(bounds.getUpperRightX(), currentCenter.lat()));                          style.setLabel(currentRadius.intValue() + "m");                     feature.setStyle(style);                     feature.redrawParent();                 }                 hasTemplate = false;             }         });      In this case, a label with radius represents after the feature added. I want to represent a radius always `if(drawFeature.isActive)`. The first idea to realize this - using the MousePosition control. Thus I can take a second LonLat parameter to calculate a radius using my custom method `calulateRadius()`. I hope that I can take the first Lonlat (= future center of circle) using mapMouseClickEvent. But if `drawFeature.isActive` it is impossible to handle this event. So, I am at a deadlock.   This is my solution that doesn't work:                  map.addMapClickListener(new MapClickListener() {                         @Override                         public void onClick(final MapClickEvent mapClickEvent) {                             MousePositionOutput mpOut = new MousePositionOutput() {                                 @Override                                 public String format(LonLat lonLat, Map map) {                                     String out = calculateDistance(mapClickEvent.getLonLat(),lonLat)+"m";                                          return out;                                 }                             };                                  MousePositionOptions mpOptions = new MousePositionOptions();                             mpOptions.setFormatOutput(mpOut); // rename to setFormatOutput                                  map.addControl(new MousePosition(mpOptions));                         }                     });


title: Where can I get documentation for Netezza Spatial Functions

text: I am looking for some documentation or manual or reference guide for Netezza Spatial Functions. Any idea where can I get the same. Thank you. Dharmesh


title: Google Maps for other planets

text: I have visited, but do not remember where, Google maps for the Moon and Mars (I think). I am wondering if anyone knows what you have to do to generate the necessary GIS data (coordinate system, etc) for a completely new planet complete with aerials? I would like to take other worlds and generate a Google Map complete with topography and draggability and more using the Google Maps API but I cannot find any information on how to do this online. I would like to do this for both planets that are real and for fictional worlds. Think Narnia as a Google Map.


title: ArcGIS Raster Symbology - Break contour at 0

text: In ArcGIS 10, how can I break raster symbology (stretched) at a particular value (i.e. 0). For example, I have a DEM representing both topography and bathymetry and I want a break at the 0 m contour line (shoreline) so the green colors represent positive elevation (land) and the blue colors represent values less than 0 (water depths). Similar to what Global Mapper offers.


title: 'Out of bounds' Error in profile51 function of matlab

text: I am trying to Extract long profile data in Matlab. I have followed the steps of this pdf. Now, in Matlab I passed the function `profile51(‘file_name’, ‘arcmap_directory’, ‘matlab_directory’, ‘n’);` and I am experiencing the error "starting point chosen geographic coordinates, out of bbounds of DEM , aborting". Where is the problem. I am using ARCGIS 10.


title: Is it possible to put scale(ex; 1:10000) at print composer on QGIS 1.8.0

text: I would like to display a cartographic scale (e.g. 1:10,000) in the print composer. I can set the scale from map tab. After setting, I can type by myself as Label. Is it possible to put scale automatically??


title: Is there a way to tell the geodatabase version using Python?

text: I want to find out if either a file or personal geodatabase is 9.3, 10.0, 10.1, etc Is there a way to do this using Python?


title: Do I have to transfer gwc-layers from my GEOSERVER_DATA_DIR directory to a production server

text: I am working with Geoserver and I know that it is possible to transfer the contents to the GEOSERVER_DATA_DIR between servers. I have already moved my GEOWEBCACHE_CACHE_DIR to a different location so that the tiles are not transferred. I am using the embedded geowebcache in Geosserver version 2.2.5. I however have a number of xml files in GEOSERVER_DATA_DIR/gwc-layers:               LayerGroupInfoImpl--570ae188_124761b8d78_-7fab.xml     LayerGroupInfoImpl--570ae188_124761b8d78_-7fac.xml     LayerGroupInfoImpl--570ae188_124761b8d78_-7fad.xml     LayerInfoImpl-18e1ce0d_13ef93f68ce_-7ffd.xml     LayerInfoImpl-1c3a926c_13ef90929b2_-7a75.xml     LayerInfoImpl-2fbcd275_13e632965e7_-7feb.xml     LayerInfoImpl--4a9eec43_132d48aac79_-7fff.xml     LayerInfoImpl--570ae188_124761b8d78_-7fae.xml     LayerInfoImpl--570ae188_124761b8d78_-7fb0.xml     LayerInfoImpl--570ae188_124761b8d78_-7fb4.xml     LayerInfoImpl--570ae188_124761b8d78_-7fb6.xml     LayerInfoImpl--570ae188_124761b8d78_-7fb8.xml      So my question is - are these files generated automatically if I tranfer the contents of GEOSERVER_DATA_DIR or do I have to include them in the transfer?


title: Prism mobile phone metadata: does it include location?

text: I've read that Prism collects metadata, but no details of what that metadata includes. Does anyone know if location is part of the metadata?


title: Can I view the Geoserver zoom level in QGIS

text: I am creating some sld's using QGIS for use in Geoserver and wanted to view the current zoom level so that I can then define it in the zoom level. Is this possible? For example                   <!-- Local Roads 1787 - 3573-->         <FeatureTypeStyle>             <Rule>                 <ogc:Filter>                     <ogc:PropertyIsEqualTo>                         <ogc:PropertyName>State_cont</ogc:PropertyName>                         <ogc:Literal>Y</ogc:Literal>                     </ogc:PropertyIsEqualTo>                 </ogc:Filter>                 <MinScaleDenominator>1787</MinScaleDenominator>                 <MaxScaleDenominator>3573</MaxScaleDenominator>                 <LineSymbolizer>                     <Stroke>                         <CssParameter name="stroke">#CDAA66</CssParameter>                         <CssParameter name="stroke-width">7</CssParameter>                         <CssParameter name="stroke-linecap">round</CssParameter>                     </Stroke>                 </LineSymbolizer>             </Rule>         </FeatureTypeStyle>         <FeatureTypeStyle>             <Rule>                 <ogc:Filter>                     <ogc:PropertyIsEqualTo>                         <ogc:PropertyName>State_cont</ogc:PropertyName>                         <ogc:Literal>Y</ogc:Literal>                     </ogc:PropertyIsEqualTo>                 </ogc:Filter>                 <MinScaleDenominator>1787</MinScaleDenominator>                 <MaxScaleDenominator>3573</MaxScaleDenominator>                 <LineSymbolizer>                     <Stroke>                         <CssParameter name="stroke">#FFFFFF</CssParameter>                         <CssParameter name="stroke-width">6</CssParameter>                         <CssParameter name="stroke-linecap">round</CssParameter>                     </Stroke>                 </LineSymbolizer>                 <TextSymbolizer>                     <Label>                         <ogc:PropertyName>road_name2</ogc:PropertyName>                     </Label>                     <Font>                         <CssParameter name="font-family">Arial</CssParameter>                         <CssParameter name="font-size">10</CssParameter>                         <CssParameter name="font-style">normal</CssParameter>                     </Font>                     <LabelPlacement>                         <LinePlacement/>                     </LabelPlacement>                     <Halo>                         <Radius>                             <ogc:Literal>1</ogc:Literal>                         </Radius>                         <Fill>                             <CssParameter name="fill">                                 <ogc:Literal>#FFFFFF</ogc:Literal>                             </CssParameter>                             <CssParameter name="fill-opacity">                                 <ogc:Literal>1</ogc:Literal>                             </CssParameter>                         </Fill>                     </Halo>                                             <Fill>                         <CssParameter name="fill">#000000</CssParameter>                     </Fill>                     <VendorOption name="followLine">true</VendorOption>                     <VendorOption name="maxAngleDelta">30</VendorOption>                     <VendorOption name="maxDisplacement">800</VendorOption>                     <VendorOption name="repeat">1000</VendorOption>                 </TextSymbolizer>             </Rule>         </FeatureTypeStyle>              <!-- Local Roads 3573 - 28573-->         <FeatureTypeStyle>             <Rule>                 <ogc:Filter>      where               <!-- Local Roads       and                            <MinScaleDenominator>3573</MinScaleDenominator>                 <MaxScaleDenominator>28573</MaxScaleDenominator>      defines what to do at different zoom levels. I need QGIS to display the above scale denominator values. Over 150k ![over 150k zoom](http://i.stack.imgur.com/XNyMU.jpg) 30k-150k ![enter image description here](http://i.stack.imgur.com/YpAMR.jpg) <30k ![enter image description here](http://i.stack.imgur.com/lQiCF.jpg) TOC ![enter image description here](http://i.stack.imgur.com/Lgwv1.jpg) The idea is to do rapid prototyping of different style options in QGIS and then export the SLD's to Geoserver for production. So people can zoom in out and turn on/off the relevant symbology. This would be easier if you can set zoom levels for different symbology from within the properties-->style tab like in arcgis. What I found for QGIS http://qgis.spatialthoughts.com/2012/02/styling-vector-data-in-qgis-using- size.html doesn't do the job as I want the layer turned on/off at different zoom levels.


title: Merge Spatially Coincident Features in the Same Feature Class

text: I need a way to programatically merge grid cells (polygon 15x15m) that are directly on top of each other and add together some of their attributes. I am working in ArcGIS 10.1 with the data stored as a polygon feature class. The data is the result of pedestrian survey across a pre-defined grid. On occasion the folks conducting the pedestrian survey get off course and record something in a column or row adjacent to the one they are working in (or the GPS position wanders). This results in two sets of data being recorded for each cell when the data from multiple GPS units are appended into a single feature class. I need a way to go through each cell, check to see if there are any duplicate features, if there are then merge their attributes (some integer, some text), and then delete the one feature that didn't get the merged data. The 'Delete Identical' and 'Find Identical' tools don't do what I'm looking for - and doing it manually isn't an option as the feature class often has upwards of 10,000 features in it. I know how to iterate through the feature class using python and check to determine if there are any duplicate features in that location. What I don't know how to do is deal with merging the features and updating the attributes. If anyone can offer some direction on how to accomplish this I would greatly appreciate it. EDIT - More Description of Before/After State: In the picture below there are two features selected - identical 15x15 meter square polygons, I need to add together all the integer fields (the ones with TOT_ as a prefix) and append the GPS_UNIT, INITIALS, and REC_DATE string fields. It's difficult to show since the features directly overlap. ![enter image description here](http://i.stack.imgur.com/1TMOA.png) **EDIT 2:** I have found that I can easily identify the duplicated features by calculating their XY centers, concatenating them into a string field as [X] & ',' & [Y] then opening the database in access and using a Find Duplicates query to identify which of the cells have been duplicated. I have not however, been able to figure out a good way to merge the attributes for each of the cells into a single row. Help with how to do that, either in python or in a database environment, would be very welcome!


title: What are my options for Server side clustering using open source GIS tools?

text: I need to make a webapp that will show over 40K Points on a map. These points are all across the globe. Using WFS Service and OpenLayers, it is possible to do clustering on small set of points, but it is not really useful to do this clustering on the client side for my case. How can I do server side clustering, in a WMS or tiled WMS service? There are several such options for Google maps. What are the corresponding tools in the OSGeo field?


title: How can I get a window by clicking on a mapsubject in the QGIS web client?

text: Is it possible in Qgis-web-client to show a new little window (with custom informations and maybe a little picture) by clicking on a mapsubject? How to modify Qgis-web-client if possible?   Maybe its possible to modify the attribute-window, that i get on the right side by using the identification cursor-mode? Thanks


title: how to use a particular output from previous code in my next code in matlab

text: I have written a code in matlab and there are several outputs. I have to use only one of the output(containing latitude, longitude, and one parameter) of these several outputs in my next code. Can anyone help me how to use only that output in my next code......and also it has a large amount of data but i am interested only in some particular range of data. How to select only that range..


title: World territorial units for statistics

text: My current assignment required to do research and produce world territorial units for statistics. In Europe it easy thanks to Eurostat. Even non-commercial GIS data is also available here. Now outside Europe I am completely lost. I know there are some regional organizations exist. What i dont know how they organize their territorial units for the member countries. Any ideas???


title: Why is torque not working with large polygons?

text: It seems that torque isn't working with a "big" polygon. By "big" I mean the size of a french department. Could you confirm that? I would like to show the evolution of jobless rate during the last 30 years. Thanks a lot.


title: How to get the distances (shortest path) between all pairs of points in the river network?

text: I have 2 vector shapefiles (one line, one point), and I have tried at least 4 ways to try to get the distance matrix at one time (and most concern the directions of each river segment)   1. quantum navigator   2. GRASS tool - v.net.allpairs   3. pgrouting   4. roadgraph plugin In the end I think that roadgraph plugin is the one which is the closest to my goal, and it also can set the directions of the river, but I can't modify it to get the distances between all pairs of points(and I want to get this in matrix form). What/How could I do? Is there any other way to get this distance matrix on qgis？


title: Arcgis Stop Working with my Loop Scrip (Could Be Memory)

text: I am having issues with my scrip, because it makes my computer very slow after 15 loops arcgis 10.1 and explorer.exe crashes. But instead of this, i can see my 15 pdf files saved correctly and i need to move the cursor to 16 for begin the loop. I am using ddp, dynamic tables (with script), GDB, everything is loaded in memory and then i delete memory every time the cicle begins. Here is my Scrip, can anyone tell me what i am doing wrong?.               import arcpy, os, sys              relpath = os.path.dirname(sys.argv[0])          #Create final output PDF     pdfPath = relpath + "\FinalOutput.pdf"     if os.path.exists(pdfPath):         os.remove(pdfPath)     finalResult = arcpy.mapping.PDFDocumentCreate(relpath + "\FinalOutput.pdf")          #Reference MXD and layers     mxd = arcpy.mapping.MapDocument(relpath + r"\MAPABASE.mxd")     LiderLayer = arcpy.mapping.ListLayers(mxd, "LS2")[0]     queryLayer = arcpy.mapping.ListLayers(mxd, r"SUBSECC_VIC")[0]     BFLayer = arcpy.mapping.ListLayers(mxd, "BF1")[0]     BFquery = arcpy.mapping.ListLayers(mxd, r"BF1")[0]          #Reference page layout elements     for elm in arcpy.mapping.ListLayoutElements(mxd):       if elm.name == "bar1":  bar1 = elm       if elm.name == "bar2":  bar2 = elm       if elm.name == "bar1txt":  bar1txt = elm       if elm.name == "bar2txt":  bar2txt = elm       if elm.name == "NoGrowth": noGrowth = elm       if elm.name == "horzLine": horzLine = elm       if elm.name == "vertLine": vertLine = elm       if elm.name == "cellTxt":  cellTxt = elm       if elm.name == "headerTxt": headerTxt = elm       if elm.name == "bar1b":  bar1b = elm       if elm.name == "bar2b":  bar2b = elm       if elm.name == "bar1txtb":  bar1txtb = elm       if elm.name == "bar2txtb":  bar2txtb = elm       if elm.name == "NoGrowthb": noGrowthb = elm       if elm.name == "horzLineb": horzLineb = elm       if elm.name == "vertLineb": vertLineb = elm       if elm.name == "cellTxtb":  cellTxtb = elm       if elm.name == "headerTxtb": headerTxtb = elm          #Reference DDP object     ddp = mxd.dataDrivenPages     #Loop through each DDP page     for pageNum in range(44, mxd.dataDrivenPages.pageCount+1):       mxd.dataDrivenPages.currentPageID = pageNum            #Graphic table variable values for LS and BF(b)       tableHeight = 3.0       tableWidth = 4.32       headerHeight = 0.16       rowHeight = 0.10       upperX = 11.18       upperY = 7.60            tablebHeight = 3.0       tablebWidth = 4.32       headerbHeight = 0.16       rowbHeight = 0.10       upperbX = 11.18       upperbY = 6.63            #Build selection set       seccnum = ddp.pageRow.seccioncom       arcpy.AddMessage("Processing: " + seccnum)            arcpy.SelectLayerByAttribute_management(LiderLayer, "NEW_SELECTION",                                               "\"seccioncom\" = " + "'" + seccnum + "'")                 #selection       arcpy.Select_analysis(LiderLayer, "in_memory\select", "\"seccioncom\" = " + "'" + seccnum + "'")       numRecords = int(arcpy.GetCount_management(LiderLayer).getOutput(0))       arcpy.Select_analysis(BFLayer, "in_memory\select2", "\"seccioncom\" = " + "'" + seccnum + "'")       numRecordsb = int(arcpy.GetCount_management("in_memory\select2").getOutput(0))            #Show selected features               #Add note if there are no Lider .................................................................       if numRecords == 0:         noGrowth.elementPositionX = 11.20         noGrowth.elementPositionY = 7.50       else:              #if number of rows exceeds page space, resize row height         if ((tableHeight - headerHeight) / numRecords) < rowHeight:           headerHeight = headerHeight * ((tableHeight - headerHeight) / numRecords) / rowHeight           rowHeight = (tableHeight - headerHeight) / numRecords              #Set and clone vertical line work         vertLine.elementHeight = (headerHeight) + (rowHeight * numRecords)         vertLine.elementPositionX = upperX         vertLine.elementPositionY = upperY              temp_vert = vertLine.clone("_clone")         temp_vert.elementPositionX = upperX + 2.00         temp_vert = vertLine.clone("_clone")         temp_vert.elementPositionX = upperX + 4.32              #Set and clone horizontal line work         horzLine.elementWidth = tableWidth         horzLine.elementPositionX = upperX         horzLine.elementPositionY = upperY              horzLine.clone("_clone")         horzLine.elementPositionY = upperY - (headerHeight - 0.01)              y=upperY - headerHeight         for horz in range(1, numRecords+1):           y = y - rowHeight           temp_horz = horzLine.clone("_clone")           temp_horz.elementPositionY = y              #Set header text elements         headerTxt.fontSize = headerHeight / 0.0200         headerTxt.text = "NAME SUPERVISOR"         headerTxt.elementPositionX = upperX + 0.80         headerTxt.elementPositionY = upperY - (headerHeight / 2)              newFieldTxt = headerTxt.clone("_clone")         newFieldTxt.text = "ADDRESS"         newFieldTxt.elementPositionX = upperX + 2.70              #Set and clone cell text elements         cellTxt.fontSize = rowHeight / 0.0200              y = upperY - headerHeight          rows = arcpy.SearchCursor("in_memory\select")         for row in rows:           x = upperX + 0.02           col1CellTxt = cellTxt.clone("_clone")           col1CellTxt.text = row.getValue("NOMBRE_LS")           col1CellTxt.elementPositionX = x           col1CellTxt.elementPositionY = y           col2CellTxt = cellTxt.clone("_clone")           col2CellTxt.text = row.getValue("DOMICILIO")           col2CellTxt.elementPositionX = x + 2.10           col2CellTxt.elementPositionY = y           y = y - rowHeight     #End of Lider Loop.................................................................................     #Add note if there are no Beneficiarios............................................................       if numRecordsb == 0:         noGrowthb.elementPositionX = 11.20         noGrowthb.elementPositionY = 6.30       else:              #if number of rows exceeds page space, resize row height         #if ((tablebHeight - headerbHeight) / numRecordsb) < rowbHeight:         #  headerbHeight = headerbHeight * ((tablebHeight - headerbHeight) / numRecordsb) / rowbHeight         #  rowbHeight = (tablebHeight - headerbHeight) / numRecordsb              #Set and clone vertical line work         vertLineb.elementHeight = (headerbHeight) + (rowbHeight * numRecordsb)         vertLineb.elementPositionX = upperbX         vertLineb.elementPositionY = upperbY              temp_vertb = vertLineb.clone("_clone")         temp_vertb.elementPositionX = upperbX + 2         temp_vertb = vertLineb.clone("_clone")         temp_vertb.elementPositionX = upperbX + 4.32              #Set and clone horizontal line work         horzLineb.elementWidth = tablebWidth         horzLineb.elementPositionX = upperbX         horzLineb.elementPositionY = upperbY              horzLineb.clone("_clone")         horzLineb.elementPositionY = upperbY - (headerbHeight - 0.01)              yb=upperbY - headerbHeight         for horzb in range(1, numRecordsb+1):           yb = yb - rowbHeight           temp_horzb = horzLineb.clone("_clone")           temp_horzb.elementPositionY = yb              #Set header text elements         headerTxtb.fontSize = headerbHeight / 0.0200         headerTxtb.text = "NAME OF BENEFIT"         headerTxtb.elementPositionX = upperbX + 0.80         headerTxtb.elementPositionY = upperbY - (headerbHeight / 2)              newFieldTxtb = headerTxtb.clone("_clone")         newFieldTxtb.text = "ADDRESS"         newFieldTxtb.elementPositionX = upperbX + 2.70              #Set and clone cell text elements         cellTxtb.fontSize = rowbHeight / 0.0200              yb = upperbY - headerbHeight         rowsb = arcpy.SearchCursor("in_memory\select2")         for row in rowsb:           xb = upperbX + 0.02           col1CellTxtb = cellTxtb.clone("_clone")           col1CellTxtb.text = row.getValue("NOMBRE_BF")           col1CellTxtb.elementPositionX = xb           col1CellTxtb.elementPositionY = yb           col2CellTxtb = cellTxtb.clone("_clone")           col2CellTxtb.text = row.getValue("DOMICILIO")           col2CellTxtb.elementPositionX = xb + 2.10           col2CellTxtb.elementPositionY = yb           yb = yb - rowbHeight            #Export to PDF Append Disabled       #ddp.exportToPDF(relpath + r"\temp.pdf", "CURRENT")       #finalResult.appendPages(relpath + r"\temp.pdf")       #os.remove(relpath + r"\temp.pdf")       #Export to Single PDF Document Each       arcpy.mapping.ExportToPDF(mxd, relpath + "\\" + seccnum + ".pdf")             #Clean-up before next page       for elm in arcpy.mapping.ListLayoutElements(mxd, "GRAPHIC_ELEMENT", "*clone*"):         elm.delete()       for elm in arcpy.mapping.ListLayoutElements(mxd, "TEXT_ELEMENT", "*clone*"):         elm.delete()       noGrowth.elementPositionX = -6       cellTxt.elementPositionX = -6       headerTxt.elementPositionX = -6       horzLine.elementPositionX = -6       vertLine.elementPositionX = -6       noGrowthb.elementPositionX = -6       cellTxtb.elementPositionX = -6       headerTxtb.elementPositionX = -6       horzLineb.elementPositionX = -6       vertLineb.elementPositionX = -6       del row       del rows       del rowsb       del horz       del horzb       arcpy.SelectLayerByAttribute_management(LiderLayer, "CLEAR_SELECTION")       arcpy.Delete_management("in_memory\select")       arcpy.Delete_management("in_memory\select2")          finalResult.updateDocProperties("MAPABASE", "AREAS", "Esri", "map sheets" , "map book", "USE_THUMBS")     finalResult.saveAndClose()


title: GeoServer subdomain, where to place a html file

text: i would like to host a website, which is addressable as a subdomain of my GeoServer. For example http://mygeoserver.com:8080/geoserver/web/my_website.html. My question is: where do i have to place the html file in my folder structure so that i can access and run it with the http request above? How can i make the web page available in the web directory of GeoServer? The GeoServer is installed under Ubuntu. Thanks in advance * * * Edit (1): Thanks for your comment BradHards. Background: The background to my question is that i want to have an html file with an OpenLayers WMS map, which shows my GeoServer data. In this map i need WFS as an overlay to the WMS map. Up to now, i can run the html file with WMS fine, but i cant get the WFS overlay to work. In this forum thread (WFS layer not showing?) someone mentioned a hint according to which the OpenLayers/WMS/WFS html file must be located within the GeoServer web directory in order to get it to work. My steps: First i installed Tomcat6 under Ubuntu 12.04 LTS (sudo apt-get install tomcat6). After that, i installed the GeoServer version 2.2-SNAPSHOT (apt-get install opengeo-geoserver). It is running well. My tries: I tried a bit for myself and i am close to a solution. In the folder "var/lib/tomcat6/webapps/geoserver/data/www" are two html demo files, ol- demo.html as demo for openlayers and wfs-t.html as demo for wfs. I can successfully call these files in the browser (localhost and from remote) by (http://myaddress.com/geoserver/www/ol-demo.html). But when i additionaly place my html file in this folder i can't reach them from the browser in an analogous way. Instead: HTTP 404 the requested resource() is not available. Hope this helps to be able to better answer my question. Thanks in advance


title: Fiona: get each feature extent (bounds)

text: I was looking in Fiona to get for each feature its extent but I didn't find how. I have try to far to do something like below               import fiona          with fiona.open('countries/ne_10m_admin_0_countries.shp', 'r') as source:         for f in source:             geom = f['geometry']             print geom      I was expecting to get a method for this at the f variable level. After some documentation reading, I've seen that f is a pure python record. So at the end, how with Fiona can I get the extent (or bounding box) of each feature geometry? PS: I already know the pure GDAL/OGR python solution so I expect a Fiona solution please


title: How to avoid creating corrupt Shapefiles during Editing?

text: I have one of my GIS technicians digitizing some lines in QGIS in shapefile format. I don't know how he did it (and neither does he), but somehow the shapefile became corrupt. It was creating random lines or some of the lines he created would just disappear. I went in into ArcCatalogue to see how it looked in ArcGIS and this is what I saw: ![enter image description here](http://i.stack.imgur.com/NItlM.jpg) Notice the question mark icon where I should see a shapefile 'line' icon. Obviously ArcCatalogue cannot read this file. Also, a second dbf file seemed to have been created with the '_packed' attached to the end. When I look at the shapefile using Windows explorer, I see that there already is a .dbf for the shapefile 'M3_PRE_SMU_lines_10Apr13_SMC.dbf', so I don't know where this _packed shapefile came from and I can't seem to find anything online that speaks to it. I tried to add this file into ArcMap and I receieved the following error: ![enter image description here](http://i.stack.imgur.com/Ljbhm.jpg) The error is pretty self explanatory...number of shapes does not match the number of records. I just don't know why that is occuring. There doesn't seem to be anything online that explains how this is occurring in QGIS, but I do see a couple repair tools. I actually repaired this myself by just opening QGIS, adding the layer, and then right clicking on the layer and 'saving as' another shapefile. So, I've figured a work around, but I'm hoping to find a solution that will stop this from occurring in the first place. Thanks, Mike


title: QGIS Equivalent of MapInfo's 'Combine Objects using Column'

text: MapInfo has a function that takes a table, and creates a new table populated combines them by a specified value. So, combining on the name field, and specifying that the Area field should be summed...               Name                  Area     Ben's Field           100     Ben's Field           50     Tom's Field           25      Becomes               Name                  Area     Ben's Field           150     Tom's Field           25      How can one achieve this in QGIS?


title: How to combine (join?) two overlapping polygons to determine combined area they cover?

text: How do I combine (join?) two overlapping polygons to determine the combined area they cover in ArcGIS 10? I am sure the Geoprocessing tools can do it, but just how eludes me.


title: Grid registration issue converting from GeoTIFF to ARC ASCII Grid

text: I've got a 3 arc second digital terrain model represented as a 32-bit GeoTIFF that looks like this:               $ gdal_info gom03_v1_0.tif          Driver: GTiff/GeoTIFF     Files: gom03_v1_0.tif     Size is 10200, 7800     Coordinate System is:     GEOGCS["WGS 84",         DATUM["WGS_1984",             SPHEROID["WGS 84",6378137,298.257223563,                 AUTHORITY["EPSG","7030"]],             AUTHORITY["EPSG","6326"]],         PRIMEM["Greenwich",0],         UNIT["degree",0.0174532925199433],         AUTHORITY["EPSG","4326"]]     Origin = (-71.499580383300781,45.999583333330001)     Pixel Size = (0.000833333333333,-0.000833333333332)     Metadata:       TIFFTAG_SOFTWARE=nc2geotiff       AREA_OR_POINT=Area     Image Structure Metadata:       INTERLEAVE=BAND     Corner Coordinates:     Upper Left  ( -71.4995804,  45.9995833) ( 71d29'58.49"W, 45d59'58.50"N)     Lower Left  ( -71.4995804,  39.4995833) ( 71d29'58.49"W, 39d29'58.50"N)     Upper Right ( -62.9995804,  45.9995833) ( 62d59'58.49"W, 45d59'58.50"N)     Lower Right ( -62.9995804,  39.4995833) ( 62d59'58.49"W, 39d29'58.50"N)     Center      ( -67.2495804,  42.7495833) ( 67d14'58.49"W, 42d44'58.50"N)     Band 1 Block=10200x1 Type=Float32, ColorInterp=Gray      Note that the `AREA_OR_POINT` attribute is set to `AREA`. When I convert this to ARC ASCII Grid format using GDAL:               gdal_translate gom03_v1_0.tif -of AAIGrid test.asc      and look at the top of the resulting ASCII file, I see:               $ more test.asc          ncols        10200     nrows        7800     xllcorner    -71.499580383301     yllcorner    39.499583333335     cellsize     0.000833333333      315.814208984375 314.712371826171875 ...      These `xllcorner`,`yllcorner` values are the same as the GeoTIFF `Lower Left` values above. From the official description of the ESRI ASCII Raster Format, however, either `xllcorner`,`yllcorner` or `xllcenter`,`yllcenter` can be specified. Since the GeoTIFF values represented AREA (not POINT), isn't the specification of `xllcorner` and `yllcorner` incorrect? Shouldn't the ASCII file in this case really read as below?               $ more test.asc          ncols        10200     nrows        7800     xllcenter    -71.499580383301     yllcenter    39.499583333335     cellsize     0.000833333333      315.814208984375 314.712371826171875 ...


title: Should point data be equidistant projected when using ArcGIS IDW spatial interpolation?

text: I am wondering if anyone can clarify whether point data should be projected to an equidistant projection when using ArcGIS IDW spatial interpolation? I am working on a dataset from Western North America spanning about 30 degrees of latitude. The data are currently in Lat/Long (NAD83). Does ArcGIS "project on the fly" or somehow adjust for latitude when calculating distances to my sample points during the interpolation procedure or should I be supplying everything in a projection that preserves distances?


title: esri.tasks.PrintTask won't print WMS layers hosted by non-ESRI servers

text: I'm using the ESRI Javascript API. My esri.Map is successfully displaying WMS layers via esri.layers.WMSLayer. Some of the WMS layers are hosted by an ArcGIS Server and some are hosted by non-ESRI servers. But when I print using esri.tasks.PrintTask, only the WMS layers hosted by ArcGIS Server are visible. The WMS layers hosted by non-ESRI servers don't show up. Has anyone succeeded in printing a WMS layer which is hosted by a non-ESRI server? -Rick


title: Syntax in Qgis "Rule properties?

text: I struggle with Qgis and some drawings of water pipes for a water company. According to the program it is possible to place different signatures for different components (stopcocks, T-pieces, etc.) From the various data located in the file "attribute table" I have to define that: the symbol is used if the value in column UNDERKLASS is equal to U22 thank you in advance Uffe


title: Has ESRI made any definitive announcements about 10.1 SP2?

text: My organization is awaiting SP2 before we can move our primary GIS team from 9.3.1 to 10.1. I have been testing 10.1 and have found a particularly bad bug that we've been told will be addressed in 10.1 SP2. However, when I look at it on ESRI's site, the bug is set to be resolved in 10.2. When is that coming out? Will there be a 10.1 SP2 in the meantime?


title: How do I make an icon with some text to popup when I hover over it and for it not to be slow?

text: I would like to make an icon on the screen and have hover popup text. Ideally the text will be black with a white background. Right now my icon is a bitmap created like this               GraphicTrackerClass incidentIcons = new GraphicTrackerClass();     incidentIcons.Initialize(mapControl.Map as IBasicMap);     int graphicID = incidentIcons.Add(thisIncident.point, alertIcon);      Then I set up a hit test and would like a text popup when I hover over the GraphicTracker (should I use a different class for displaying the icon?). `AxMapControl.DrawText` is fast, but I can't figure out how to customize the text and background MapTips look really cool but I can't find a function to plot a maptip with a point and text - I've only found `layer.ShowTips = true;` I tried creating a balloon callout and activating my IElement like this (but it is really slow about 2 seconds for the popup to appear) :               pGraphicsContainer = activeView as IGraphicsContainer;      pGraphicsContainer.AddElement(pElement, 0);     pElement.Activate(activeView.ScreenDisplay);     activeView.PartialRefresh(esriViewDrawPhase.esriViewGraphics, pElement, null);      I also had a bad experience experimenting with dynamicdisplay. _**Below I've left the original post in case anybody wants to read it:_** What are the main steps to use to improve performance for an arcengine hosted in WPF with C#. It's super slow. I made a bunch of layer invisible which helped but its still no where near being snappy / usable. Edit: I load DataMaps10/streets_na/StreetMap North America.mxd and turn make the following layers invisible CanadaBoundaries, CoastalWater, Parks, Airports, Landmarks, UrbanAreas, USBoundaries(World-25m), USBoundaries(25-1.5m), USBoundaries(1.5m-0), CanadaBoundaries, CanadaCoastlineandProvinces(World-25m), CanadaCoastline(25m-0), CanadaProvinces(25-1.5m), CanadaProvinces(1.5m-0), Exits(Numbers)(15k-0) I haven't tried switching to Winforms, I don't exactly know how to do that. Edit 3: SO the problem is DynamicDisplayEnabled = true ruins EVERYTHING! so I turn it off... but now I want text that pops up when you hover over something... I was going to use dyanmic display for that, but dynamic display turns everything to crap ! I'm super frustrated.


title: How to drag a point and then use savestrategy in OpenLayers

text: Im trying to use the wfs-t in openayers using the savestrategy. I have created a savestrategy and it works when I for example uses the DrawFeature Control to add a point. A new point is added to the database. However, when I use the DragFeature control to move a point (in the same vector layer as I have added using DrawFeature) it doesnt work. The response from Geoserver is               <wfs:InsertResult>        <ogc:FeatureId fid="none"/>     </wfs:InsertResult>     <wfs:TransactionResult>        <wfs:Status>          <wfs:SUCCESS/>        </wfs:Status>      </wfs:TransactionResult>      And to me that looks like for some reasons no information sent regarding the fid that should have been moved. What could I havce done wrong? BR Mike