title: How to normalize choropleth maps of census information?

text: I'm trying to improve on the census-based maps on this site, which shows census values at 5 scales from Collection Districts through to State. Currently the maps show whole numbers (eg total population per polygon) which is not optimal. How should they be normalized? **Option 1: % of Total** I was advised to normalize by % of the total. This results in legend ranges such as "0.0004% to 0.0006% of the Total Population" which (IMO) make the legend hard to decipher - it's not an intuitive fraction, and is hard to put into perspective. **Option 2: Population Density** Another option could be to normalize based on area. Since the datasets are stored in meters, this also results in values like "0.000000019 to 0.000000035 persons per square meter". I could flip this to show "area / person" - but this wouldn't make sense for datasets such as housing. And would it be offensive when applied to datasets such as religion (1km2 per X person)? Is there a standard methodology to convert either of these values into something more meaningful? What do other people do on their census maps? Thanks, Steve


title: Shp's will NOT line up  (coordsys) - shps made using GPSeismic sotftware (client made shp's)

text: OK, I asked this previously but solved it on my own last time. Same project, but differant shp files. I recieved from a client shp files that were made in GPSeismic 2011. Upon recieving the shp files from client, he sent the 5 files (dbf, shp, shx, txt and sp1) as you can see in the screen shot. The POREC did not have the prj (coordinates) so I put in the correct coordinates that are being used. I've tried the coords that we are currently using, "undefined" coords, even did the same with the surface shp file prj, switched that to 'undefined'. Nothing seems to work. The postplot keeps falling into the same place, just to the East even with differant coords. What is going on with this, I make the coordinate changes in ArcCatalog as I have ArcView 9.3. I have two projects as this is the 2d one, the first one I had no problem with the coords as the client sent the shp files made with the GPSeismic software program. My Proj Manager could not get this to line up on her computer either. Waht could be the problem AND SOLUTION? Could the SP1 file be the problem? ![enter image description here](http://i.stack.imgur.com/RtA6b.png) * * * ## sp1 file below (partial) ![enter image description here](http://i.stack.imgur.com/hLHeJ.png)


title: How to add super/sub script on a QGIS text annotation?

text: How do I add super/sub script on a QGIS text annotation? For example, 44.95 mi^2 ![example image with superscript 2](http://i.stack.imgur.com/q1Fr5.png)


title: Natural Earth: How to get the position of all provinces that belong to a state

text: This question is about Natural earth 1:10 states and provinces.shp file. I'm using ogr2ogr and want to know the name of the id i need to pass . In purpose of pulling out the right states regions. I've use couple of properties name that didn't work (`-where "adm0_a3 IN ('MLI')" \` , `-where "iso_a2 = 'ML'" \` ). Can anyone help me i'm running out of places to find the answer.


title: Why don't libgrass.so.1.8.0 and libqgissqlantconnection.so.1.8.0 load?

text: I have compiled Qgis 1.8 at Ubuntu 12.10, but some phyton plugins does'nt work, like sextante. Trying to solve this problem, some files don't have been load. They are: libgrass.so.1.8.0 and libqgissqlantconnection.so.1.8.0. Any ideia about this?


title: Conversion from 3 dimension geometry to 2 dimension

text: I have a Point table in which the geometry column contains both 3 dimension and 2 dimension values.I want to convert this 3 dimensional geometry to 2 dimension.Is there any function in PostGIS? Please help me to solve this problem. Regards Sreesha.T.A


title: POI Datasets for North America

text: I use foursquare because it contains a lot of POIs. Which database do you use when you need a point of interest in North America?


title: Call ArcGIS Server Buffer command without using SOM

text: All, I wanted to see what the best way was of calling the geometry service to carry out a Buffer operation on ArcGIS Server 10, without using the SOM. Ideally I want to consume the REST API, using a wrapper so that it was easy to call from c#. Is the GISClient the way forward? From what I can determine, it looks like you only get stateless (and thus low footprint) access to the Mapserver and Geocode services. I know that I could construct a call using a URL, but I wondered if there was C# API that wrapped the REST API for you. Cheers, Matt


title: features supported by WCS in geoserver?

text: i have implemented a geotiff image layer using WCS service in geoserver 2.2.4. I have searched on net for operations supported by WCS...it tells me about 3 operations: GetCapabilities, DescribeCoverage and GetCoverage I have used these operations using http request and all are working perfectly. My questions are:   * is there any other way to use these operations other than http request (if yes, please explain with some example)   * is wcs support any other feature/facility/operation other than the above mentioned operations. Is it possible to apply stretch operations(standard deviation/histogram etc) and other operations on geotiff image in geoserver using WCS?


title: how to add GML files in geoserver 2.2.2

text: How can I add GML files to a GeoServer 2.2.2 instance? How can edit them "on the fly" and save the results back to the server using OpenLayers? I would appreciate examples if possible. sir i m adding a vector layer. in which we want to use gml. we want to add feature like point ,line etc. to gml. when we do editing on line that features will be saved on this layer. i want to insert that feature using openlayers on line like in wfs layer we do using transanction. similar type i want to do using gml .... please provide some solution..


title: .NET And ArcGIS 10x - Where have all the wrappers gone?

text: With prior versions (to ArcGIS 10x) of ArcGIS, the .NET wrappers used to be placed in a dotnet directory. Where are they located with 10x? With prior versions, I used the absence or presence of this directory to tell whether or not the .net framework was installed in the proper order with installation of ArcGIS. How can I tell if the .net framework was installed in the correct order with respect to ArcGIS???


title: ArcGIS 10.1 Network Analysis- How do I get appropriate network dataset?

text: I have acquired London road dataset from the Open Street Map however after creating and building the network dataset in ArcCatalogue the solve tool could not return the routes because of errors no valid stops. That after I had enter the addresses as stops, No result was returned. Olawumi


title: How to create a webmap showing point density?

text: I am trying to develop hotspot map showing point density in our gis based web application. After searching for an option I have decided to use v.kernel command. I have a table containing some point data ie crime data. I need to display the hotspot of crime for a particular area. In v.kernel vector file containing points and polygon are used for displaying hotspots/heatmaps. In qgis there is an option for selecting some rows and exporting it to shape file. Is there any option for running this command from a python file or is there any better option for displaying hotspot/heatmaps? I have tried openlayers heatmap, python's heatmap and django-heat app but nothing seems to work for me. I am using openlayers tms(Openlayers-2.12) for displaying maps generated through mapnik and tilecache. Thanks in advance.


title: Finding outliers with wrong Z values in LIDAR data

text: I have more than 2000 XYZ (only ground points) .txt files and I have to do Quality Control (QC). Let's say I have file, which I converted to multipoint (or even singlepart feature). I have mean Z value for a file, but if most of the points are 20-40 metres, and two, three points per file are 70-80 metres (wrong Z values). How to catch these points and remove? On another hand, it could be that in another file 70-80 metres is common Z value and some points wrong Z value 120 metres. In another file 120 metres is common and correct. How to find and remove points with wrong Z values in LIDAR data? Standard deviation? Raster cells?


title: Plotting two shapefiles using R?

text: I'm not familiar with using R for GIS. I was wondering if you can help me figure out how to plot different shapefiles at the same time. I'm using the "map" function to generate a map in order to get the boarders of the country I'm working on. Then I would like to add on top of it 20 shape files representing each a state in a way that they look adjacent to each other and fit the country's borders.


title: How to draw graphics in a child window?

text:
I am adding some graphics to the map in a parent window, and I would like to show the same graphics in a child window. I am using the ArcGIS Javascript API (version 3.1). I have written code as below: Complete code in parent window:               <%@ Page Language="C#" AutoEventWireup="true" CodeBehind="Test1.aspx.cs" Inherits="ePermit.Test.Test1" %>          <!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">                        <meta http-equiv="Content-Type" content="text/html; charset=utf-8" />         <meta http-equiv="X-UA-Compatible" content="IE=7" />                  <meta name="viewport" content="initial-scale=1, maximum-scale=1,user-scalable=no" />         <title>GIS Search </title>         <link rel="stylesheet" type="text/css" href="http://serverapi.arcgisonline.com/jsapi/arcgis/3.1/js/dojo/dijit/themes/claro/claro.css">         <%--<link href="Styles/claro.css" rel="stylesheet" type="text/css" />--%>                                                                     <form id="form1" runat="server">

                         Content for the second tab                          Layer List :

                         Content for the third tab

                 <%--

--%>

                                                             <%-- <button onclick="map.panUpperLeft();"></button>--%>                                                             	                  <%--<button onclick="map.panUp();"></button>--%>                                                              	                 <%--<button onclick="map.panUpperRight();"> </button>--%>

                                                              <%--<button onclick="map.panLeft();"></button>--%>

                             	                                                  <%--<button onclick="map.panRight();"></button>--%>

                                                              <%--<button onclick="map.panLowerLeft();"> </button>--%>                                                              	                 <%--<button onclick="map.panDown();"></button>--%>                                                              	                 <%--<button onclick="map.panLowerRight();"></button>--%>

                 This is the results section
               <asp:HiddenField ID="hfxmin" runat="server" />          <asp:HiddenField ID="hfxmax" runat="server" />          <asp:HiddenField ID="hfymin" runat="server" />          <asp:HiddenField ID="hfymax" runat="server" />              </form>                     Complete code in Child window:               <%@ Page Language="C#" AutoEventWireup="true" CodeBehind="WebForm3.aspx.cs" Inherits="ePermit.WebForm3" %>          <!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">                   <meta http-equiv="Content-Type" content="text/html; charset=utf-8" />         <meta http-equiv="X-UA-Compatible" content="IE=7" />                  <meta name="viewport" content="initial-scale=1, maximum-scale=1,user-scalable=no" />         <title>GIS Search </title>         <link rel="stylesheet" type="text/css" href="http://serverapi.arcgisonline.com/jsapi/arcgis/3.1/js/dojo/dijit/themes/claro/claro.css">         <%--<link href="Styles/claro.css" rel="stylesheet" type="text/css" />--%>                                                                    <form id="form1" runat="server">

            <asp:HiddenField ID="hfxmin" runat="server" />          <asp:HiddenField ID="hfxmax" runat="server" />          <asp:HiddenField ID="hfymin" runat="server" />          <asp:HiddenField ID="hfymax" runat="server" />          </form>                       When user click on print button child widow should open with the same graphics as in parent window. Thanks, Uday


title: How to center a map in Print Composer to the same center as in map window?

text: I can't figure out how can I **center the map in Print Composer to layer/object/map center** in the same way as in map window in **QGIS 1.8**. Let's say I have layer which I want have in the center of the map. In map window it is pretty easy - I just click on " **Zoom to layer extent** ". Unfortunately when I create map frame in Map composer, the layer is centered on width (good) but not centered on height (I guess it is aligned to bottom). If I use button _Set to map canvas extent_ in _Print Composer / Item Properties / Extents_ it is again centered on width but this time aligned to top. I think in both cases it would be much better to have it centered. **Is there any way to have map centered in Print Composer to the same center as map in map window?**


title: Batch Export feature attribute to ascii

text: I have around 160 feature classes that I need to convert into ascii files. Spatial statistics tool ---> Utilities ---> Export feature attribute to ascii only allows converting one file at a time. Is there any way to export all my features to ascii in one go? Any script written? Thanks


title: How to create Sequential Animation of Polylines in ArcGIS10?

text: I have a data set consisting of sets of polylines which reflect movement of individual organisms in my sample within their habitat each day. I am trying to create an animation which cycles through the days, displaying only the movement for one day at a time. When I try to do this in ArcGIS10, it adds the lines from each day onto a map containing the lines from all previous days. How can I get the previous days' lines to disappear when the time slider moves to the next day?


title: How to run two versions of QGIS simultaneously

text: I really enjoy seeing the new features and icons of QGIS 1.9 nightly versions. It's looking fantastic, and some features are already very useful for my work. However, it is inherently less stable than QGIS 1.8 as it's a development version. To deal with this, I generally keep QGIS 1.9 installed and updated on my system, and use it until something I need does not work (this happens every fortnight or so - most commonly with plugins such as ftools and openlayers). At this point, I uninstall v1.9, remove the ubuntugis repository, and proceed to install v1.8. This solution works fine, but I was just wondering if there's a better way? How would I go about running 2 versions of QGIS simultaneously (with updates), one for stability, the other for features, fun and bug-checking? And is it worth doing? One solution I tried a while back was to clone the github version into a new folder on my machine and run that binary when I wanted the latest version...


title: Create Polygon Overlay Handling Holes

text: How can I create a polygon overlay that correctly handles holes? This answer provides a method of creating a polygon overlay that does not handle holes. Can it be modified to correctly handle holes? Software Versions:   * PostgreSQL 9.2.1   * PostGIS 2.0.1 r9979   * Windows Server 2008 R2


title: How to create a point cloud from DEM?

text: Please forgive my ignorance if I am asking a totally irrelevant question, I am still a GIS newbie. I looked in Google and found a lot of links as of how to convert point clouds to DEMs. I am trying to achieve exactly an opposite: to create a point cloud from a DEM of the specific location and then to create a 3D model of this location from this point cloud. I didn't find one single link about it. So is all my thinking wrong? Is there other way to achieve what I am trying to achieve?


title: How to create shapefile from OS MasterMap data using ArcGIS Desktop 10?

text: How do I create a shapefile from **OS MasterMap** data in Arcgis 10? I have tiles from the Ordnance Survey Site and would like to create a topomap from the existing tiles. These tiles contains the data I require for the creation of the shapefiles. I will be grateful if someone could help me out with this. Thanks


title: Is it possible to have a dynamic table in a data driven page environment?

text: I am creating a map series that has many polygons per sheet. There is annotation for the polygons however some polygons are very small with no annotation for them (annotation would completely hide a polygon). Thus, I have an attribute table where lables are stored for those specific polygons. How can I have that table displayed on the map, with only the attribute values per sheet, rather than all 2000+ values? I am looking for a simple, dynamic solution.


title: Export polygon to points (vertices)

text: This is the first time that I have worked with geospatial data and QGIS, so please be kind. I was given a shape file that contains house districts for the state of FL. This link provides instructions (for ArcMap users) to import the house districts into Tableau, which is my goal. I was able to generate a file that contains the district, and coordinates, but I need the order column, and am not sure how to do this using QGIS (assumedly using an add-in).               "shapeid","x","y"     "0","-87.634896","30.866261"     "0","-87.634871","30.866344"     "0","-87.634851","30.866536"      The format that I am looking for is ...               "shapeid","x","y","order"     "0","-87.634896","30.866261","1"     "0","-87.634871","30.866344","2"     "0","-87.634851","30.866536","3"


title: How to change selection tolerance for WMS/WFS in QGIS?

text: Does anyone know if there is any equal function to the selection tolerance (ArcMap) in QGIS. I have a problem activating points in my WMS/WFS. The "selection tolerance" function solves the problem in ArcMap (If I raise the pixel numbers there it works), but I would like to find a similar function in QGIS. Does any one know?


title: How to integrate .NET application in QGIS using Python?

text: In our company we have developed a .NET application that's now integrated with 3 different GIS platforms. Now we would like to make it work with QGIS but we need some help on getting started. I've been thinking of using callbacks and I've looked at this **posting** on NathanW's blog but I can't make it work in QGIS. Any ideas?


title: Anything similar to TMS or WMTS for WFS?

text: I've been working with OpenLayers for a while, but problem got a question with WFS. Not sure if it's supported by WFS, or OpenLayers, I want to achieve something like this, mostly for cache purpose: As user pans/zooms the map, OpenLayers keeps sending out WFS requests with BBOX. I want to make coordinates in such requests keep only 3 decimal digits, in a sense, aligning to 0.001 degrees. Do I have to tinkle with OpenLayers, or is it built-in? **EDIT** I'm thinking about tiling for WFS, like TMS and WMTS do. Just we don't render image here, instead we "render" XML files.


title: Integer rasters referring to the same scale

text: I'm using ArcMap 10. I want to display the area of simulated avalanches. I want this area to be divided in different colors depending on the thickness of the deposit. I managed to do that. But I want all the simulated avalanches to refer to the same scale and same legend of thickness. This means that for every simulation a certain range of values of thickness (e.g. from 2m to 4m) falls within the same color. The problem is I need to transform my raster into polygons (for future calculations of area ratios) and to do this it must be integer. So I guess I have two choices: "raster calculator" or "reclassify". With "raster calculator" I see I obtain my integer raster with the number of classes of thickness I want but the areas are changed a bit (I don't know why). With "Reclassify" I obtain an integer raster where my areas are essentially the same as before, but it can have a number of classes not higher than the value of thickness it reaches. In other words using the "reclassify" tool I would have each simulated avalanche with its own scale of thickness and not a common thickness for all my simulations. Is there a way to accomplish what I'm looking for? I hope I explained it clear enough (It's a bit hard to explain this even in my own language). Thank you, Michele


title: Need to give same color for Points and its label

text: ![like this I have to give](http://i.stack.imgur.com/tZF4L.jpg)I have 10000 points, using Graduate Color I have given different color for each points in 10000. I want to give a same color for each label according to each points colors using ArcGIS. is it possible, How to do it. the above image is model for my work


title: Query Postgis Table with PHP

text: Would you be so kind as to help me with the following in OpenLayers. I need to be able to use the coordinates from a user click in a postgis query. To catch the coordinates of a click I built a click handler to use the function `map.getLonLatFromViewPort(ext.xy)`, which delivers me the coordinates of my click in the following format: lon = 568417.33622114, lat = 681985.1902844. These are sent to a php-script using a GET request. I then want to query my postgis db using `ST_GeometryFromText`, which requires that I put my coordinates in the following syntax `((568417.33622114 681985.1902844))`. Any help on how to get my coordinates in the correct format would be appreciated!


title: How do I change geotiff "no data" color to white/transparent rather than black?

text: I am preparing data for inclusion into a geonetwork node, where the maps are displayed via geoserver. The problem is that cells that have"no data" are displaying as black in the geotiff and I am would like the color to be set to white or transparent. I am using ArcView and do not have access to either Spatial Analyst or 3D analyst. I would prefer not to have to assign a new sld file and am presently using "Use colormap" when exporting from Arc. Any help with this would be greatly appreciated. Thanks


title: Can Geoserver return the raster value of a lat/lon point

text: I have a GeoTIFF published by Geoserver 2.2 and being used by Leaflet via WMS. How can I get Geoserver to return the raster value at a lat/long point? I'm guessing it has to do with using WCS, but no examples seem to exist for doing this! This will be similar to querying PostGIS with `ST_Value(raster, ST_SetSRID(ST_Point(lon,lat),4326))`. I dont really want to import the same raster into PostGIS just to find the raster value at a point. Or is this the recommended solution?


title: Opening GeoServer Admin returns HTTP ERROR: 503?

text: I try to open Geoserver Admin and returns the following error. > HTTP ERROR: 503 > > Problem accessing /geoserver/web/. Reason: > > Service Unavailable What is causing this?


How can I convert a set of ShapeFiles to an OSM file?


title: Arcgis JS api WMS layer on top of basemap

text: I'm trying to place WMSlayer on top of basemap "streets", but have ran into problems. Code:               // mapHandler.mapextent = new esri.geometry.Extent(350000,6370000,750000,6630000,new esri.SpatialReference({"wkid":3301}));          mapHandler.map = new esri.Map('map', {     // extent: mapHandler.mapextent,     center: [25.88270117, 58.70311058],     basemap: 'streets',     spatialReference: new esri.SpatialReference(4326),     zoom: mapHandler.defaultZoom     });     mapHandler.wms1layer1info = new esri.layers.WMSLayerInfo({title:"Maa-ameti fotokaart"});     mapHandler.wms1resourceinfo = {extent: new esri.geometry.Extent(300000,6300000,800000,6700000,{wkid:3301}),layerInfos:[mapHandler.wms1layer1info],version:["1.1.1"]};     mapHandler.wms1layer1 = new esri.layers.WMSLayer("http://kaart.maaamet.ee/wms/fotokaart?",{resourceInfo:mapHandler.wms1resourceinfo,visibleLayers:["MA-FOTOKAART"]});     mapHandler.wms1layer1.setImageFormat("jpg");     // mapHandler.wms1layer1.setVisibleLayers(["0", "1", "2"]);     mapHandler.wms1layer1.setImageTransparency(false);     mapHandler.map.addLayers([mapHandler.wms1layer1]);      With the above code only the basemap is displaying. I'm getting the following warning in console:               Resource interpreted as Image but transferred with MIME type application/vnd.ogc.se_xml: "http://kaart.maaamet.ee/wms/fotokaart?SERVICE=WMS&REQUEST=GetMap&FORMAT=ima…=2445863.801486572,7992358.782554072,3272606.6994188153,8240626.2504242575".       If I comment out:               basemap: 'streets',      Will result in WMSlayer displaying correctly and also prints the following the console:               Map: Geometry (wkid: 4326) cannot be converted to spatial reference of the map (wkid: 3301)       Any help would be appreciated.


title: Calculate Field using Esri Python Addin

text: I am trying to select a polygon then it must automatically change (calculate) the status to "C". The tool seems to work in the background but nothing actually happens (the status is not updated). I have used the code from this post but nothing actually gets selected. If I have a feature in the layer selected already it does deselect it though. Below is the code that i am using. Is there something that I have gotten wrong? Thanks for any help.               import arcpy     import pythonaddins          def getSearchDistanceInches(scale, selection_tolerance=3):         """Returns the map distance in inches that corresponds to the input selection         tolerance in pixels (default 3) at the specified map scale."""         return scale * selection_tolerance / 96.0 # DPI used by ArcMap when reporting scale          class tlStatToC(object):         """Implementation for PythonAddin_addin.tlMapCoods (Tool)"""         def __init__(self):             self.enabled = True             self.shape = "NONE" # Can set to "Line", "Circle" or "Rectangle" for interactive shape drawing and to activate the onLine/Polygon/Circle event sinks.         def onMouseDownMap(self, x, y, button, shift):             mxd = arcpy.mapping.MapDocument("CURRENT")             pointGeom = arcpy.PointGeometry(arcpy.Point(x, y), mxd.activeDataFrame.spatialReference)             searchdistance = getSearchDistanceInches(mxd.activeDataFrame.scale)             lyr = arcpy.mapping.ListLayers(mxd)[0] # assumes you want to select features from 1st layer in TOC             arcpy.SelectLayerByLocation_management(lyr, "INTERSECT", pointGeom)             arcpy.RefreshActiveView()             arcpy.CalculateField_management(lyr,"STATUS","C")


title: How to get a message when the map is moved in QGIS?

text: I want to get a message when every time I move the map. I had:               def info(self):         QMessageBox.information( self, self.tr( "Info" ), self.tr("info") )      I suppose that used to be something like that...               iface.mapCanvas().???.connect(self.info)


title: Retaining original thematic colors in a USA source shapefile

text: > **Possible Duplicate:**   >  Adding Layerfiles (.lyr) to Qgis New to QGIS and need a lot of help! I am using QGIS with MAC - I have downloaded a USA thematic map but when I open it in my system it changes the colours. I need the original colours! There are over 30 categories and they are very difficult to match, so manually changing them is not an option. How can I retain the original colours? The downloaded zip file includes a .shp, .dbf, .sbn, .prj, .sbx, .lyr, .xml as well as a .kml file. Thank you any help you can provide


title: Getting MapInfo instance from C# code

text: Please sorry for my bad english :( Here is the situation: I have a .Net DLL written on C#, witch contains windows form. MapInfo calls function from this DLL witch opens form. Is there any way to get instance of MapInfo, that called this function? (Not just first opened instance, because i can have multiple opened)


title: QGis 1.7.4 on Ubuntu 12.04

text: Will there be available stable version of Qgis 1.7.4 in Ubuntu 12.04? Version 1.9 is working fine at the moment, but as it is unstable yet, dont want to risk with my work. Thanks.


title: How do you change the prime meridian of a map?

text: I am working with the Natural Earth world map as a SHP file, but when I try to change the prime meridian (`PRIMEM`) of the coordinate system > > GEOGCS["GCS_WGS_1984",DATUM["D_WGS_1984",SPHEROID["WGS_1984",6378137.0,298.257223563]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]] I get this: ![enter image description here](http://i.stack.imgur.com/vGGVi.png) Projections seem to give the same results. What is happening here?


title: QGIS: WMS misalignment

text: I am experiencing a small screen space offset when using a WMS layer (origin is GeoServer 2.1.4): ![Alignment visualisation](http://i.stack.imgur.com/lj6ZQ.png) The image shows a red square which is the top left pixel of a tile (GeoTIFF). An exact copy of the same tile is among the tiles in the background provided by my GeoServer. The four images represents the same view (the red square is the same pixel in all four images) just panned a bit inside my QGIS window. The yellow line outlines the corresponding pixel in my WMS layer - added to illustrate the misalignment. I would expect the red square to always align to the yellow line. **Can anyone provide me with some knowledge to why I'm having this misalignment depending on the extent shown on my screen?** _I have posted a similar questionhere proposing that the misalignment could be from a difference in the aspect ratio between the width/height and the bbox in the WMS request. But I'm not sure this is the reason and thus the reason for posting two different questions on the same topic._


title: Summary statistics on attribute table

text: I've got a polygon feature class with attributes fields `y2010`, `y2011`, `y2012`, `y2013`, `y2014`, `country` The 'year' fields contain counts of an event occurrence, which can be 0 or > 0. The 'country' field' is an identifier for country A or B, so the table looks like this:               y2010 | y2011 | y2012 | y2013 | y2014 | country             0  |  10    |   1   |   0   |   5   |   A        1  |   0    |   3   |   0   |   5   |   A        0  |   1    |   1   |  45   |  20   |   B      I'd like to reformat this table to show _per year_ the _number of features_ that fit certain criteria. The criteria are:   1. More than zero event occurrences for countries A+B;   2. More than 1 event occurrences for countries A+B;   3. More than zero event occurrences for country A. so that the final table using the example above would looks like:               criterium | 2010 | 2011 | 2012 | 2013 | 2014 |        1           1      2      3      1     3        2           0      1      1      1     3        3           1      1      2      0     2         and the table content of columns 2:6 would be the number of features of the original table, that meet the criterium _for that specific year_. I've been toying around with the Statistics toolset, but so far without success. Any tips are welcome. Edit: Opted for Q&D solution by exporting the attribute table to R and wrote a script to calculate the cell values.


title: Creating a new raster based on pixel value of original raster using raster calculator but new layer is not evident

text: I am trying to create a new raster based on a current raster (w001001 - NVIS) which has one paletted band, with palette values ranging from 1 to 29, with each palette/value representing a vegetation type (e.g. 1 = rain forests). I want to create a new raster (actually, preferably a vector, but I cross that bridge later) for the palette of value 3 (this one:) ![enter image description here](http://i.stack.imgur.com/XIcWR.jpg) I have followed the answers to the questions posted in this forum, and have used the Raster Calculator to do this: ![enter image description here](http://i.stack.imgur.com/4JXm0.jpg) It all seems to process fine, and the new layer appears in my layer list. However, I can't see the layer at all and when I use the ![enter image description here](http://i.stack.imgur.com/NCMam.jpg) nothing comes up. When I open properties for the layer, metadata tells me this: ![enter image description here](http://i.stack.imgur.com/hnu8W.jpg) Is there something obvious that I am missing? This is my first post and I am a fairly new user of QGIS so apologies if anything is unclear or nonsensical. Thanks in advance.


title: Which is the best way to store roads for the shortest path finding application?

text: I would like my application which is simply a shortest path finder in a city to meet the following requirements:   * the city roads network should be stored persistently on disk and loaded into memory on demand (to execute the shortest path finding algorithm);   * I would like to use something like OpenLayers to view the results that have been yielded by the shortest path algorithm;   * the 'OpenLayers view' should contain not only the results, but the map (uninteractive) of the city itself. As for now I am thinking of the following architecture:   * PostGIS enabled database to store roads ('roads' table);   * PostGIS enabled database to store the actual map data to render it in OpenLayers;   * a table for the results of shortest path algorithm in the PostGIS enabled database (I would like to view a couple of shortest paths simultaneously);   * a web app with OpenLayers which shows the map from the 'map' table, shows the roads layer from the 'road' table and the shortest path from the 'shortest path' table. Does this sound correct or should I change anything?


title: QGIS Relation Manager

text: I have a PostgreSQL database with a one-to-many relationship. Here is a reduced example:               CREATE TABLE buildings.building     (       objectid serial NOT NULL,       building_no character varying(40) NOT NULL,       geometry geometry(Polygon,32638),       CONSTRAINT building_pkey PRIMARY KEY (objectid)     )          CREATE TABLE buildings.flat     (       objectid serial NOT NULL,       fk_property_group integer,       fk_building integer NOT NULL,       CONSTRAINT flat_pkey PRIMARY KEY (objectid),       CONSTRAINT flat_fk_building_fkey FOREIGN KEY (fk_building)       REFERENCES buildings.building (objectid) MATCH SIMPLE      )      I defined a relationship in the QGIS (2.2 and 2.4) Project Properties and it looks good in the forms. But when I enter a new building I will have the value `nextval('real_estate.property_group_objectid_seq'::regclass)` for the objectid of the building and the flats that I want to add won't have the correct foreign key. Its actually NULL. Is there any workaround/solution for this issue?


title: How to create a 9.3 file geodatabase with ArcObjects SDK for ArcGIS 10?

text: In ArcGIS 10 I could use the Create File GDB geoprocessing tool to create a 9.3 file geodatabase or the Create Personal GDB geoprocessing tool to create a 9.3 personal geodatabase. But I cant't find a similar method in the ArcObjects SDK for ArcGIS 10. The _IWorkspaceFactory2.Create_ method seems only to create geodatabases compatible with ArcGIS 10. Possible workarounds are:   * using the geoprocessing interface   * centrally storing one empty 9.3 geodatabase and using that as template geodatabase But is there really no possibility to create a 9.3 Geodatabase with the ArcObjects SDK for ArcGIS 10?


title: Resume Tilemill MBTiles export after crash

text: Is there a way to resume a Tilemill export (MBTiles) after it has crashed? Today I just got the following crash:                FATAL ERROR: Evacuation Allocation failed - process out of memory      This was just a little sample export, but I will need to make a pretty huge export (size of Europe with Zoom Level from 0 - 18) and if that would happen inbetween this export I'd like to know that I don't loose everything, but just a little part and that Tilemill will continue afterwards... Is there a button for resuming that I am missing? Or what would be the most reliable solution for this process? Thanks


title: How to change shape file's extent in QGIS

text: I have created a network from 3 shapefiles using GRASS. While haveing problem to perform a network analysis with it, I've been told that the extent of the network shouldn't exceed that of the original shape file. Since the extent of each original shape is smaller than that of the network, I am looking for a way to modify the extent of the original shape files in Qgis. Can someone help me with it? Many thanks, Dor


title: spatialpolygons heatmap in R - NAs in grey?

text: I have a spatial polygon data-frame and I am trying to plot a heatmap via spplot. It's a map of Europe and the values for the heatmap are in a column of the data-frame. Plotting it works fine but I have some NAs in the column used for the colors and these countries are just not colored, so they are white. I'd really like them to be grey, however.. I have no idea how I can do this in spplot.. can anyone help?


title: add attributes to newly created shapefile from existing shapefile using arcpy

text: I have created a new shapefile and want to add shape attribute from another existing shapefile. It will be a recursive process with update of shape attributes of new shapefile (ofc) more than 2 times. Any help is appreciative. Following is what I am trying to make work:               Import arcpy     ofc = "C:/outFC.shp"     ifc = "C:/inFC.shp"     fields = arcpy.ListFields(ifc,"*","ALL")     for f in fields:         shape = f.Shape         query = '"Shape" = \'' + str(shape) + '\''         cur = arcpy.UpdateCursor(ofc,query)          del f, fields      Diagram: I want to get FID (0) from first table and add it to second table of new shapefile and only need FID, Shape, and id fields in second table. ![enter image description here](http://i.stack.imgur.com/jSjTv.png) ![enter image description here](http://i.stack.imgur.com/qPsVt.png)


title: Determining Pixel-by-pixel Catchments in ArcHydro

text: I used ArcHYDRO to prepare HydroDEM (5m) without sinks, Flow direction (FDR) raster, flow accumulation raster (FAC). As I understand FAC, value of pixel shows number of pixels that contribute to this pixel. Now, let say I derive drainage line from FAC (in raster form), I would like to link somehow all pixels that contribute to particular pixel in my drainage line, in way that I have separate contributing areas for each pixel in that drainage with clear relation to this pixel. How do I do that? The main goal is to derive height difference between pixel in the drainge and any pixel that contribute to it and only to it.


title: Conditional field calculator expression

text: I need to calculate the PIPEID field to equal the G_KEY_COMP field if PIPEID is null or PIPEID and G_KEY_COMP do not match. I think the first condition of the PIPEID being null is covered by the second condition of the fields not matching? I'm a novice on this topic and can't get the syntax I have to work. Here's what I've piecemealed together: Pre-Logic Script Code:               def update(PIPEID,G_KEY_COMP):      if PIPEID != G_KEY_COMP       return G_KEY_COMP     else:      return PIPEID      PIPEID = update(!PIPEID!, !G_KEY_COMP!) I'm posting a screenshot of the field calculator window since I'm not able to get the indentation right in the codeblock. ![enter image description here](http://i.stack.imgur.com/4rY91.jpg)


title: Mysterious Invalid Syntax in Field Calculation

text: I'm trying to calculate a field in Model Builder using the Calculate Field tool. I'm using the following code block and expression in the calculation, but am getting an "invalid syntax (, line 1)" error.               # Code Block:     def get_block(block_id):         block_only = str(block_id)[1:]         if block_only[-1].isalpha():             return block_only.zfill(3)         else:             return block_only.zfill(2)          # Expression:     get_block(%block_id%)      **%block_id%** is an in-line substitution that points to a string variable I create as the output of a custom Python script tool. Any ideas what the error is?


title: ESRI Javascript API: setting feature layer selection symbol

text:
This seems like it should be a relatively simple thing to do, but for whatever reason I can't figure it out. I'm trying to add functionality to select a group of parcels as laid out in this example, and display queried information for those parcels in a popup info window. I've got it displaying the queried data in an infowindow correctly, and it highlights the correct parcel for the info being displayed, but I'm having issues highlighting every selected parcel. In the above example it sets the selection symbol with setSelectionSymbol(), but for whatever reason this isn't working for me, and I'm not exactly sure what I'm doing wrong. I've included the parts of my code that are relevant to my problem, including parts that deal with toggling the parcel layer on and off in a legend. Thank you for any help! If i can help clarify anything at all let me know!               function init() {              //initialize popup window          var popup = new esri.dijit.Popup({         fillSymbol: new esri.symbol.SimpleFillSymbol(esri.symbol.SimpleFillSymbol.STYLE_SOLID, new esri.symbol.SimpleLineSymbol(esri.symbol.SimpleLineSymbol.STYLE_SOLID, new dojo.Color([255,0,0]), 2), new dojo.Color([255,255,0,0.25]))         }, dojo.create("div"));              //add geometry service from server         esri.config.defaults.geometryService = new esri.tasks.GeometryService("http://2501s-jjohnson:6080/arcgis/rest/services/Utilities/Geometry/GeometryServer");              //initalize map         app.map = new esri.Map("map", {           basemap : "satellite",           center : [-111.924531, 40.611871],           zoom : 13,           infoWindow: popup         });              //query parcel and show info in popup window          function executeIdentifyTask(evt) {           var query = new esri.tasks.Query();           query.geometry = pointToExtent(app.map,evt.mapPoint,10);           var deferred = featureParcel.selectFeatures(query,esri.layers.FeatureLayer.SELECTION_NEW);           app.map.infoWindow.setFeatures([deferred]);           app.map.infoWindow.show(evt.mapPoint);           console.log(evt.mapPoint);              }              //trigger feature layer parcel query if parcels are turned on           function clickConnect(connect){             if(connect){                 //perform the identify task on click                  clickHandler = dojo.connect(app.map, "onClick", executeIdentifyTask);             }else{                 //disconnect the click handler                     dojo.disconnect(clickHandler);               clickHandler = null;               }               }              //add parcel layer          var parcels = new esri.layers.ArcGISDynamicMapServiceLayer("http://2501s-jjohnson:6080/arcgis/rest/services/Parcels_UT/MapServer", {           id : 'parcels',           visible : false         });         parcelLayer.push({           layer : parcels,           title : 'Parcels'         });              //add parcel feature layer query info for popup window         var content = "Address: ${ADDR}" + "
Owner Name: ${OWNER}" + "
Parcel ID: ${APN}" + "
City: ${CITY}" + "
Acres: ${TOTAL_ACRES}"+ "
County Assessor Site"            var popUpTemplate = new esri.InfoTemplate("Parcel", content);                 //add parcel feature layer         featureParcel = new esri.layers.FeatureLayer("http://2501s-jjohnson:6080/arcgis/rest/services/Parcels_UT/MapServer/0",{           mode: esri.layers.FeatureLayer.MODE_SELECTION,           outFields: ["*"],           infoTemplate:popUpTemplate         });              //set selection Symbol         var selectionSymbol = new esri.symbol.SimpleFillSymbol().setColor(new dojo.Color([255,255,0,0.5]));         selectionSymbol.setOutline(new esri.symbol.SimpleLineSymbol("solid", new dojo.Color([255,0,0]), 2));         featureParcel.setSelectionSymbol(selectionSymbol);              //create legend         dojo.connect(app.map, 'onLayersAddResult', function(results) {           var legend = new esri.dijit.Legend({             map : app.map,             layerInfos : legendLayers           }, "retail");           legend.startup();         });         //add all layers to map         app.map.addLayers([parcels]);         //add layer arrays into  legend object         legendLayers["parcel"] = parcelLayer;              dojo.connect(app.map, 'onLayersAddResult', function(results) {             //for each layer array in legend object             for (var k in legendLayers){           //add legend check boxes for each layer           //to toggle layer visibility           dojo.forEach(legendLayers[k], function(layer) {             var layerName = layer.title;             var checkBox = new dijit.form.CheckBox({               name : "checkBox" + layer.layer.id,               value : layer.layer.id,               checked : layer.layer.visible,               onChange : function(evt) {                 var clayer = app.map.getLayer(this.value);                 clayer.setVisibility(!clayer.visible);                 this.checked = clayer.visible;                     //toggle the feature layer parcel query                //parcelClick indicates whether the parcel layer is on                 if (this.value ==="parcels" && this.checked === true){                     clickConnect(true);                     parcelClick = true;                 } else if (this.value ==="parcels"){                     clickConnect(false);                     parcelClick = false;                 }               }             });                  //add the check box and label to the table of contnents             dojo.place(checkBox.domNode, dojo.byId(k), "after");             var checkLabel = dojo.create('label', {               'for' : checkBox.name,               innerHTML : layerName             }, checkBox.domNode, "after");             dojo.place("
", checkLabel, "after");           });           }         });              dojo.connect(app.map, "onLoad", initSelectToolbar);            }                  function initSelectToolbar(map) {       selectionToolbar = new esri.toolbars.Draw(map);       var selectQuery = new esri.tasks.Query();            dojo.connect(selectionToolbar, "onDrawEnd", function(geometry) {         selectionToolbar.deactivate();         selectQuery.geometry = geometry;         var deferred = featureParcel.selectFeatures(selectQuery, esri.layers.FeatureLayer.SELECTION_NEW);         app.map.infoWindow.setFeatures([deferred]);         app.map.infoWindow.show();            });     }            function pointToExtent(map, point, toleranceInPixel) {        var pixelWidth = map.extent.getWidth() / map.width;        var toleraceInMapCoords = toleranceInPixel * pixelWidth;        return new esri.geometry.Extent( point.x - toleraceInMapCoords,                     point.y - toleraceInMapCoords,                     point.x + toleraceInMapCoords,                     point.y + toleraceInMapCoords,                     map.spatialReference );                                  }                dojo.ready(init);


title: Changing kml marker in esri javascript api

text: I can load my kml file but I want to be able to define a marker for my points. I can't find anywhere in the api which allows me to define the marker. How can you define your own marker using a kml file with the esri javascript api?


title: Implementing custom IRenderer for ArcGis 10.0/WPF

text: **EDIT:** Updates in the bottom As a followup to my earlier question, I'm trying to implement custom symbol rendering for my WPF/ArcGIS application. I'm accessing layers that have been published on the ArcGIS server, with pre- existing symbology already in place. I want to keep this symbol, and add text to it, since FeatureLayers don't automatically display their data labels. I'm running into several problems. I've defined my renderer thus:               public class EntityLabelRenderer : IRenderer     {         public Symbol GetSymbol(Graphic graphic)         {             var objectId = graphic.Attributes["OBJECTID"]; // always set.             var graphicFromServer = //Use a QueryTask to fetch the object;                  return graphicFromServer.Symbol; //add text         }     }      When I load my layer, GetSymbol is called for each object, but the `Graphic` object I get in the call is very barebones - the Symbol property is null, and the only values I can get are the ID. So I use that ID, as I show above, to fetch the full Graphics object with all its properties using a QueryTask. But even the one I get now has its Symbol set to null. Additionally, I don't know what the default data label for this object is, as set on the server. I have all the fields, but have to decide myself which to display. So my question is: 1) Can I get, for a Graphic object fetched from an ArcGIS server, the default data label as defined in the MXD? 2) Get I get the default symbol for an object, as defined in the MXD on the server? Can I get that, then create a composite symbol comprised of both that symbol and an additional TextSymbol? 3) Is there another, better way? The only one I could think of was to create a custom GraphicsLayer and add TextSymbols to it manually for each object of mine, but this seems insanely complicated. ----- **Update** ----- I've gone over this and tried some things, and here's where I am in my attempts to get this working: If I wait until my FeatureLayer's `Initialized` event is fired, I see that the layer comes populated with an IRenderer for the server-configured symbols. This is great, since I now have the symbol I want to display, as well as the layer's definition of the default text (to display as a label). So now I can create my own IRenderer which wraps this server-generated IRenderer (usually a UniqueValueRenderer or UniqueMultipleValueRenderer) and adds the label. I can see two approaches: 1) Hacking into the existing renderer, changing its ControlTemplate to add my label, or 2) Somehow returning a composite Symbol that contains both the existing Symbol and my own TextSymbol. Option #2 is better, but it seems that unlike the Flex and Java SDKs, the WPF SDK doesn't have a CompositeSymbol object. I'd rather not go for option #1 before I'm sure that a CompositeSymbol can't be implemented somehow. Does anyone know how I can go about this?


title: How to display intersecting roads and create multiple labels for long roads in QGIS?

text: I am new to GIS and I am using Qgis 1.7.3 on Windows. Can any one please tell me how to display a map from line shapefile with these two characteristics:   1. Roads are displayed with borders at the sides and intersect at crossroads   2. Name label for long roads are duplicated at specified intervals Here is the link to an image that illustrates my idea: ![https://picasaweb.google.com/lh/photo/MlvnVJmY1MMXqI5AR-B66fiYNl9YxKwsJ_DSpLaCxjs?feat=directlink](http://i.stack.imgur.com/4P5zB.jpg) (The map was taken from _underdark_ 's blog) This is the best I can achieve by drawing lines with two overlapping simple line, one slightly bigger than the other: ![https://picasaweb.google.com/lh/photo/80tHFuFnw0SrXG2Us8jLgPiYNl9YxKwsJ_DSpLaCxjs?feat=directlink](http://i.stack.imgur.com/csMlN.jpg) Thanks in advance!


title: is there any direct access to Mapbox data? (similar to CartoDB's sql api)

text: I've been asked to add a form to a map that will write data in the background (e.g., ajax to a magic file that does datawrites and some emailing, notifications etc) and a new marker will be added to the map being displayed. **I have done this before** with CartoDB SQL api. Sent all the form data to a magic ajax script that wrote the new info to the database, then the front-end would refresh the map to show the new marker. So that's fine, and I had intended upon copying my code over wholesale, but then I found out this is in MapBox ... Is there any way to write to the MB database? or will i need to re-work the entire map storing the data elsewhere ?


title: Renderd outputs "failed to initialize projection 4326"

text: I'm setting up a tile server using the guide provided here. I have everything set up. `proj4` is installed in version 4.8.0. The espg file is existant, I ran an `xmllint` over my xml style file without errors. `mod_tile` is also running without problems. I also have a fully populated PostGIS database running. I used `osm2pgsql`. The problem is starting `renderd`. It gives me the following error message:               renderd[25257]: An error occurred while loading the map layer 'default':      failed to initialize projection with: '+init=epsg:4326'       encountered during parsing of layer 'necountries' in Layer at      line 58 of '/src/mapnik-style/osm.xml'      (I inserted line breaks for legibility) It's the only error message. Everything else works fine. I checked for the espg file at `/usr/local/share/proj/epsg` which is present and contains (of course) the WGS84 definition. I googled around but couldn't find something. There were hints to the proj4 library. But since it's installed in version 4.8.0, i can't see the problem. Thanks for the help!


title: Is there an option to remove duplicate polygon labels in QGIS?

text: I would like to know if there is an option in QGIS labeling engine which would remove all the duplicate labels. My polygon layer contains a few features with same name. So when labeled it becomes a cluster which does not look good.


title: Grass algorithms doesn't appear in the sextante toolbar in QGIS

text: I was having problems with some sextante algorithms in my QGIS 1.8, something related with wrong "os" name, I found a post in which reccomend to install the 1.9 master version. I clean uninstalled the QGIS and then I installed the 1.9 version. When I activate the sextante toolbar, the Grass, TauDEM and even the QGIS algorithms doesn't appear in it. I have tried activating them in the configuration of sextante, reinstalling the plugin but the tools still missing. Any ideas of what is happening?


title: Qgis. I got a error message about sqlite3_step while I was doing tutorial. Help me please

text: I just started using Qgis. Eventually, I got successful installation of QGIS but I'd been through hell of installation problems. Though, I got this a message from doing tutorial. I was converting landuse.shp to spatialite (See below) ![Error occured](http://i.stack.imgur.com/gUKgj.png) Does anyone have any idea what the problems is ? Suggesting solution would be thankful. Here is where I got tutorial. http://www.qgis.org/en/docs/training_manual/foreword/preparing_data.html


title: Visualize given temperatures on image/map/plane

text: I am looking for a way to visualize some data of the form :               latitude longitude temperature      **as image(s)**. Have many of these 3-tuple lines uniformly distributed on a region. Ultimately, I am looking for a dynamic way, but even a statically creation of an image would be fine. I need the image in order to use it as an overlay for google maps (or something similar) and looking for an end result like the one here.


title: Rasterize vector layers with different extents into identical raster grids with QGIS

text: I am working on a quite simple spatial potential analysis using QGIS. I have some ArcGIS knowledge, but I am new to QGIS. My task is to identify the functional inner city of a region based on different criteria like density or distance to public transport. Therefore I created polygon layers representing different criteria. All the polygons have values between 0 and 1. The next step is to rasterize the vector layers and use the raster calculator to make an additions for each raster cell. The problem I am facing is that the layer extents are different. Therefore all the raster outputs also differ in extent and cell size. As a result I cannot use the raster calculator to make the cell calculations. So what I need to do, is to find a way to rasterize the polygone layers into raster sets with the same extend and the same cell size. ArcGIs offers the options to set a specific “Analysis extent” and a specific “Analysis cell size” for further raster calculations, which is basically what I am looking for in QGIS. I hope you guys can help me!


title: EPSG: 900913 x Display EPSG:4326

text: I have a question when the projected 4326 and 900913. The open layers uses by default 4326 using the projection coordinates in decimal fomato already uses google maps with coordinates 900913 meters. I saw several examples of creating the map in the projection and display property 900913 projection in 4326. What difference does it make on the map? If I want to record a point in the database will be in the map projection 900913 in this case?


title: ST_Distance is not delivering right results

text: I want to calculate the distance between lat/lon points. I tried the following queries:               Select      ST_DISTANCE(ST_GeomFromText('POINT(-73.97923 40.75571)', 4326), ST_GeomFromText('POINT(-73.974509 40.757108)', 4326))       which gives 0.00492364143698933 as a result. and               Select      ST_DISTANCE_SPHERE(ST_GeomFromText('POINT(-73.97923 40.75571)', 4326), ST_GeomFromText('POINT(-73.974509 40.757108)', 4326))      which gives 426.951968421872 as result. If I look at the distance between this point in Google Maps http://goo.gl/maps/GvhFM I get 1,1 miles (~1770m) as a result. Why this mismatch? How can I get the real (Google Maps) distance?


title: Showing a point shapefile with raster layers in geoserver

text: I'm new to geoserver, and I am fiddling around with the geoserver suite for Windows. I have two aerial images georeferenced and imported as GeoTIFF which worked out fine. I am also trying to import a point shapefile with two points in it, but it wont show. The CRS of all three layers are the same, and if I just preview the shapefile, it looks good. You can see how my setup looks here: ![enter image description here](http://i.stack.imgur.com/XTagd.png)


title: QGIS OS grid overlay

text: I'm totally new to GIS. However, I've proudly managed to create a map, which is correctly geo-referenced with a Great Britain Ordnance Survey (OS) background (streetview level) and I've added some point data. I would like to know how to create or import a layer with the OS grid lines. I'm sure this has been done a million times, but I can't seem to find the answer on the web - I suspect I'm not using the correct search terms. Any help appreciated.


title: ArcGIS map service not resizing for print-friendly page

text: ( I asked this **same question on ArcGIS forums** , with no response...I hope someone here can help!) I am using ArcView, ArcServer 10.0 and the JS API v3.1 (so no print function) ... I have an app where the user chooses parameters from a few dropdown menus, and a print-friendly map page opens in a new tab. From this page, the user can click a 'print' button, which uses window.print() to open the browser's print dialog. I have a print.css page, which is working well with my 'normal' divs, text and images. I have this tag in my HTML:               <link rel="stylesheet" href="print.css" type="text/css" media="print" />      My problem is that I cannot seem to reset the size of the Dojo map object. The map in my printed page (or PDF, if I choose CutePDF) is the same size as whatever the size was in the print-friendly page. I would like my print.css to hard define the printed version / PDF as an exact height and width, regardless of the window size in the print-friendly page. I cannot define the size of page to be printed in the print-friendly page, b/c for one, my client would like to pan and zoom the map after the print-friendly page is generated. I have no problem setting the size of the other elements in my print.css rules (text box, hard coded legend with images)...but any rules I try to set for the map itself in my print.css page are ignored. I tried using 'exportMapImage', thinking that maybe I could make an image of the map and its contents, but this loses the basemap which I need, and also I do not want the user to open yet another page or tab. (I am having the same problem with my scalebar, but I this isn't so important for this project). My goal is to create a PDF or printed page that is the same size no matter what the screen size or current window size., and I have it working except for the map itself. Any ideas? I hope I explained my problem clearly. Thanks , J Also, I am hoping to find a client-side solution; I don't have much experience with server-side programming. If this can't happen, please let me know any server-side ideas, thanks!)


title: Convert .tiff data to bitmap just when rendering and at runtime by visualC++ code

text: I need to convert a .tiff image to a bitmap one without saving the bitmap. I know that I should use GDAL_translate but because I'm new, this function's interface is a little complicated for me. Could you please tell me how should I call the function in this special case?   And if you have any better suggestion than using GDAL_translate, please allow me be informed? ## Note that I don't want to save the bitmap. It's just for displaying purposes in mfc and my superior has told me not to save the bitmap just convert data to bitmap at run-time. well @user1240592,@Mapperz,@Alex Marcov,@Jamo,@Tonny,@johanvdw and @soheil I've seen the question you've answered here. my problem is I don't need a bitmap file to be created on the hard disk. suppose that I have this MFC application ![enter image description here](http://i.stack.imgur.com/z479M.png) already I have written the code to print header information of a geotiff in the rich edit box now I need something to convert a 15931x16997 geotiff to bitmap at runtime and show it in a picture control ( I need a high performance way) and be able to do various image processing operations on this bitmap data at run-time (using opencv) and immediately show the result on the picture control and also be able to do these tasks when the image is viewed:   1-zoom in   2-zoom out 3-pan 4-fixed zoom in 5-fixed zoom out 6-go back to previous extent   7-go next extent   and then if needed I can convert the bitmap data created by image processing operations to a geotiff image file and convert it in the hard disk.   I don't have problems with image processing parts,opening and saving geotiff images, extracting raster data and write raster data.   my problem is converting data to bitmap at runtime and use a high performance way (like openGL) to display the image with high performance and do various navigation tasks on it.   **Do you anyway or any library that have implemented these tasks Notice just need bitmap data on the buffer not saved on the hard disk?**


title: Create .avl file Table

text: i would like to create a table of the .avl file i use as my symbology layer but cant come up with a way to do it in arcmap 10. Any thoughts?


title: Creating Periodic table with geoserver

text: I'm trying to create a periodic table using geoserver and openlayers. I'm trying to do this in geoserver as I want the periodic table to change colors. But I'm new to geoserver and not sure which configurations(like layers) to use and how to create shapefiles for this.Any help would be really appreciated. Thanks for your time and help


title: Intersecting Layers Returns Empty Collection Despite Visible Intersection

text: When I perform an intersection of a polygon I draw in OpenLayers and a postgis database layer, it seems that I am getting incorrect results. Intersection works correctly on some layers. For instance, if I intersect a triangle with a layer of polygons that represent crop fields, I get the following: ![enter image description here](http://i.stack.imgur.com/GLxna.png) The query my app generates to produce the above result is:               SELECT ST_AsText(ST_Intersection(%(geometries_0)s::geometry, %(geometry)s::geometry))      where `geometries_0` is my triangle:               POLYGON((-104.84928345939991 40.518951354186285,-104.82319093011056 40.51953858115158,-104.83700967095314 40.50707521626648,-104.84928345939991 40.518951354186285))      and `geometry` is my layer of crop fields, as well-known text:               MULTIPOLYGON(((-104.841309611298 40.5075331998226,-104.84173356681 40.5069932245841,-104.842041204329 40.50640946683,-104.842224948796 40.5057962996657,-104.842280275816 40.5051688207073,-104.842205823049 40.5045424803865,-104.842003423773 40.5039327015263,-104.841678061729 40.5033544995574,-104.841237748411 40.502822112724,-104.840693325791 40.5023486513933,-104.840058199365 40.5019457751149,-104.839348008051 40.5016234053897,-104.838580239118 40.5013894812384,-104.837773797582 40.5012497635973,-104.836948540713 40.501207693373,-104.836124789073 40.5012643066572,-104.83532282616 40.5014182091969,-104.834562398965 40.5016656107496,-104.833862231727 40.5020004184754,-104.833239564888 40.5024143870601,-104.832709730574 40.5028973218633,-104.83228577506 40.5034373300773,-104.831978137541 40.5040211136997,-104.831794393074 40.5046342970926,-104.831739066055 40.5052617810515,-104.831813518821 40.5058881146554,-104.832015918097 40.5064978757385,-104.832341280141 40.5070760506105,-104.83278159346 40.5076084036796,-104.833326016079 40.5080818278834,-104.833961142505 40.5084846673069,-104.834671333819 40.5088070040565,-104.835439102753 40.5090409023397,-104.836245544289 40.5091806037522,-104.837070801158 40.5092226689759,-104.837894552799 40.5091660624086,-104.83869651571 40.509012177646,-104.839456942906 40.5087648031902,-104.840157110143 40.5084300292289,-104.840779776982 40.5080160977709,-104.841309611298 40.5075331998226)))      However, if I perform the same query with a different layer ("soils"), I get an empty result: ![enter image description here](http://i.stack.imgur.com/91qJS.png) The query is the same:               SELECT ST_AsText(ST_Intersection(%(geometries_0)s::geometry, %(geometry)s::geometry))      with a polygon `geometries_0` that should overlap:               POLYGON((-104.84627938530097 40.54511058649626,-104.83460641167578 40.545175808723876,-104.84070039055733 40.537283458057615,-104.84627938530097 40.54511058649626))      and a `geometry` layer representing soils, similar to the crop fields in the above query:               MULTIPOLYGON(((-104.939716 40.258166,-104.939775 40.258174,-104.939963 40.258159,-104.940159 40.258065,-104.940039 40.257671,-104.939917 40.25749,-104.939928 40.257419,-104.94003 40.257404,-104.940265 40.257641,-104.940632 40.257902,-104.940826 40.258061,-104.941051 40.258188,-104.941123 40.258235,-104.941205 40.258283,-104.941246 40.258275,-104.941287 40.258212,-104.941186 40.258094,-104.941186 40.258007,-104.941167 40.257921,-104.941105 40.257858,-104.941044 40.257786,-104.941045 40.257716,-104.941127 40.257676,-104.94122 40.257653,-104.94141 40.257731,-104.941559 40.257671,-104.941255 40.257181,-104.940857 40.256794,-104.940644 40.256478,-104.940319 40.255997,-104.940003 40.255728,-104.939676 40.255561,-104.939419 40.255544,-104.938895 40.255529,-104.938287 40.255512,-104.938046 40.255528,-104.937733 40.255549,-104.937322 40.255533,-104.937012 40.255577,-104.936947 40.255593,-104.936623 40.255774,-104.936581 40.255924,-104.93658 40.256042,-104.936661 40.256223,-104.93671 40.256436,-104.936842 40.256618,-104.937262 40.256753,-104.937662 40.256818,-104.937897 40.257,-104.938181 40.25745,-104.938374 40.257742,-104.938465 40.257931,-104.938782 40.258051,-104.939121 40.258092,-104.939439 40.258133,-104.939716 40.258166)))      I use the postgis `ST_AsText` function to convert database layers to well- known text, and I checked to make sure that all layers have the EPSG:4326 projection (using the `Find_SRID` function). Why would one layer (crop fields) intersect correctly, and another (soils) not? I've tried the same query using geographies instead of geometries, with the same result.


title: Python: using CalculateField_management for duplicating attributes with a double quotation mark

text: I am writing a python script which duplicates a field. The field contains names of pipelines but some of these names has double quotation (ex. 8" NGL PIPELINE). I have created a new field and trying to copy the attributes using: arcpy.CalculateField_management(PL_XING_SUMMARY,"PLNAME", "'!FIRST_PLNA!'", "PYTHON_9.3"); But the script double quotes whole attributes (ex. "8" NGL PIPELINE"). Is there a way to remove the double quotation marks at the start and end of each attribute? P.S. I am using ArcGIS version 10.1 Service Pack 5


title: How to set Zero values to NULL using QGIS pycalc?

text: I'm using QGIS 2.4 pyclac and windows 7 OS. I have columns that should show NULL value but instead of that they show Zero (0) value, I need those fields to be NULL. This far I have seen a script that will solve my problem: ![enter image description here](http://i.stack.imgur.com/MP5sB.png) TIPO_USO is a string of characters, but VOL_COMER is double, so I need to pass it to NULL as the first field. I'll use the next script, using PyCalc: ![enter image description here](http://i.stack.imgur.com/iDGuK.png) Then I get this message:               Field code block does not return 'value' variable! Please declare this variable in your code!      Now I don't know what to do because I don't work with python but I'm sure this will solve my question, I hope :D


title: How to force polygons to the bottom, without compromizing drawing speed?

text: I have a large polygon dataset and no time to rebuild it to regoranize the Object IDs. Many polygons are drawing above others, blocking their symbol (for example on one side of the polygon). I tried using the 'advanced' symbol level drawing however this resulted in the polygons drawing at a glacial pace. For my clients & users, this is not cool. Is there any way to force the polygons I don't want up front to the bottom, without compromizing drawing speed?


title: How to implement VisibleBoundsUpdatedEventHandler

text: I have an ArcObjects add-in button that adds attribute value info as a text graphic into a map document layout. Now, instead of using the add-in button click event, I would like the code to run on VisibleBoundsUpdatedEventHandler or when the active view has been updated. Can this be implemented through add- ins? The link below uses this for a panning tool. I'm not sure how and where to begin to implement this. Any suggestions or code samples would be most appreciated! http://forums.esri.com/Thread.asp?c=93&f=993&t=145041 Thank you


title: Point-in-polygon using Python and OGR

text: I want to use OGR in Python to write a simple point-in-polygon test:               def point_in_polygon(point, polygon):         """         point : [longitude, latitude]              polygon : [(lon1, lat1), (lon2, lat2), ..., (lonn, latn)]              """         # Create spatialReference         spatialReference = osgeo.osr.SpatialReference()         spatialReference.SetWellKnownGeogCS("WGS84")         # Create ring         ring = osgeo.ogr.Geometry(osgeo.ogr.wkbLinearRing)         # Add points         for lon, lat in polygon:             ring.AddPoint(lon, lat)         # Create polygon         poly = osgeo.ogr.Geometry(osgeo.ogr.wkbPolygon)         poly.AssignSpatialReference(spatialReference)         poly.AddGeometry(ring)         # Create point         pt = osgeo.ogr.Geometry(osgeo.ogr.wkbPoint)         pt.AssignSpatialReference(spatialReference)         pt.SetPoint(0, point[0], point[1])         # test         return pt.Within(poly)      However, it doesn't seem to work properly:               In [22]: polygon = [(-10., -10.), (-10., 10.), (10., 10.), (10., -10.), (-10., -10.)]          In [23]: point_in_polygon([0., 0.], polygon)     Out[23]: True          In [24]: point_in_polygon([359., 0.], polygon)     Out[24]: False      Any ideas what I'm missing? Thanks for your insight :)


title: Alias Name Issue ArcGIS Explorer

text: Question regarding arcgis explorer I seem to be running into an issue related to obtaining and using alias names in a datagridview for the columns. I cannot seem to use or even aquire the Alias names but clicking on a point on the map does open the popup with the names intact. I even have tried just looping through the original table and it still does not seem to display the alias names event from there (When using columns.getaliasnames) the names that appear are the original column name. Has anyone else encountered this issue or something similar to this? Thanks


title: Is it possible to use Python to perform network analysis on only part of a network?

text: I have a large network dataset and want to perform an analysis on it which would output the start/end point, the distance between them and whether they're first, second or third (etc.) nearest neighbours. There are approximately 30,000 nodes in the dataset, so I thought it might be possible to process 1,000 at a time (still comparing to all 30,000 nodes, but only doing the analysis on part of the dataset). I have an idea of how to do it (using a for loop, creating new files for each part of the analysis and putting it together) but I'm not certain as to how to create a for loop to use only part of the data and to go through and eventually finish analysis on all data. Does anyone know if this is possible? **Edit:** I've been having thoughts about the partitioning by rowID and was wondering then if I were to be making a service area analysis layer if then the layer would be impossible to create due to the partitioning or otherwise incorrect because of it?


title: How to tag schools in OpenStreetMap?

text: Hello and good evening dear GEO-friends, How can we tag schools? One problem is that they are based on large areas, and there are entrances and barriers. How can we handle footways or all the parts of a footway or service road from the main road network so that it is possible to route to the entrance of the school. Any ideas here?! Is it right that generally speaking school buildings can be tagged with building=school?


title: Win7 Senior users_Which pgRouting and postGIS versions are stable with postgreSQL8.4?

text: I'm working on a postgreSQL ,postGIS and pgRouting project on Windows 7   I need to work on a stable version (no time for risks), and I can see that many people are asking the same question many places..   I think because there is not a clear instruction about the compatible and stable versions anywhere, so we (Win-pgRouting beginners) would appreciate a complete answer from senior users..      **Here are the tools I would like to use:**   PostgreSQL 8.4   postGIS ??   pgRouting ??   QGIS ??   OpenLayers ??   GeoServer ??   osm2po ??      **Here are some facts I found on during my search about stable versions compatible with pgRouting**   1) pgRouting on postgreSQL 9 is not a stable version (it is in beta form). postgis.net/windows_downloads      2) Expert users confirming that version 8 is used by most of them who are working on the collection of tools that I'm interested in rndxpress.wordpress.com/2012/07/11/…         3)Problems with the new version collections gis.stackexchange.com/questions/26816/…      4)This is a pgRouting4Windows mail list response, it is written couple of days ago it is stating explicitly by the Win pgRouting binaries' developers that it does not work properly with PostgreSQL 9 and they are at this moment forming a group of developers who are willing to build a stable version pledgebank.com/pgRoutingWindows **THE QUESTION:**   From your experience using this collection of tools, what versions do you suggest that I should use? I installed postgreSQL8.4, but I’m not sure which compatible versions I should install for postGIS and pgRouting and if you have any suggestions for QGIS ,OpenLayers and osm2po, because I'm new to the area so I'm not sure if compatibility even matters for some of those tools,         I appreciate your support,   sara


title: QGIS automatically inserts segments when editing polygon nodes

text: I have created a shapefile with multiple polygons. When I attempt to move one of the veritices with the node tool new sements automatically appear that cross the polygon and even cross the opposite poygon border to a vertex on another polygon. I can sometimes clean this up and toggle editing off. Then I "Check plygon validity" and there are no errors, but when I toggle editing again the offending segments reappear. I started again with a new QGIS project and a new shapefile layer (not snapping to the "bad" layer) and got the same problem at another location. Same problem occurs if "avoid intersecting polygons" is on or off. Snapping distance is set to 10 m, distance is >100m to the vertices that the new lines automatically go to. One, two, or three new lines may radiate from the same "offending" vertex. Working with this new layer and project, if I close QGIS and reopen, the shapefile looks good, but the probelm reappears when I start to edit the node. I was able to delete the offending node, but it will reappear when I edit vertices weith the node tool. Help! (V 1.8, Windows). I don't recall this happening with v 1.6.


title: import libtiff error

text: I am unable to use libtiff from Python on windows 7. Here is what I get when I try from libtiff import TIFF I do not understand the statement " Make sure that libtiff is installed and its location is listed in PATH|LD_LIBRARY_PATH|.. PythonWin 2.7.6 (default, Nov 10 2013, 19:24:18) [MSC v.1500 32 bit (Intel)] on win32. Portions Copyright 1994-2008 Mark Hammond - see 'Help/About PythonWin' for further copyright information. > > > from libtiff import TIFF Traceback (most recent call last): File "", > line 1, in File "D:\Python27\lib\site-packages\libtiff__init__.py", line 4, > in from .libtiff import libtiff, TIFF File "D:\Python27\lib\site- > packages\libtiff\libtiff.py", line 35, in raise ImportError('Failed to find > TIFF library. Make sure that libtiff is installed and its location is listed > in PATH|LD_LIBRARY_PATH|..') ImportError: Failed to find TIFF library. Make > sure that libtiff is installed and its location is listed in > PATH|LD_LIBRARY_PATH|..


title: Geocoding in 64 bits - background geoprocessing

text: I have 10.1 SP1 installed plus the background geoprocessing patch for 64 bit geoprocessing. Im trying to geocode a large file, so wanted to take advantage of 64 bits processing. My Surprise when i try to change properties of locators to make use of runtime memory limits to more the 2024 it gives me an error saying out of range limits. How can i process set up locators as described in this blog post I also have enabled the background geoprocessing in ArcCatalog > Geoprocessing options and created a model to run the job, but when looking at taskmanager at run time its always a process of arccatalog (32bits) running and nothing related to 64 bits. Thank you in advance all the support.


title: What are the units of the calculated area from a raster in decimal degree

text: I'm trying to calculate area for a raster file which contains 2 layers (0 for unsuitable habitat and 1 for suitable). It's been calculated successfully but I don't know what is the unit of the calculated values. Calculated Area for 43873037 counts:365608.662141 and for 3241800 counts: 27015.001513 Spatial reference: WGS84_world_Mercator, Linear unit: Meter, Angular unit: Degree, Cell size:0.0083333338, GRID format


title: Using PNG/JPEG Compression in GeoServer

text: Is there any use in setting the PNG and JPEG compression in GeoServer (Services - WMS - PNG Options/JPEG Options)? The default is 25, but we heard rumours that having compression resulted in a time penalty for decompression on the way out, so we set it to 0. We deliver PNG8 and JPEG tiles, and these are relatively small but if we can make them smaller without significant quality or time loss then that would be great. We just aren't sure whether sticking compression on would reduce the image sizes but result in a (significant) time penalty or not and would like to know for certain. FYI We can't really just test it as currently we do not have a Dev GeoServer environment! Thanks!


title: GeoServer REST API: Creating Database Data Store does not persist password

text: I'm passing the XML created by the method below to a GeoServer REST API request that successfully creates a postGIS data store. However, password is not persisted. I have to update the password manually in the data store. I need the password to be saved for subsequent requests using the data store to be successful. Any suggestions on what I can do for the password to be saved?                   private string getDbXml(string dsName)         {             string dbHost = "localhost";             string dbPort = "5432";             string dbName = "testgis";             string dbUser = "mydbuser";             string dbPassword = "mydbpassword";             string dbType = "postgis";             string exposePKsParameter = "<entry key=\"Expose primary keys\">true</entry>";                  string xml = "<dataStore>" +                             "<name>" + dsName + "</name>" +                              "<enabled>true</enabled>" +                              "<connectionParameters>" +                                  "<host>" + dbHost + "</host>" +                                  "<port>" + dbPort + "</port>" +                                  "<database>" + dbName + "</database>" +                                  "<user>" + dbUser + "</user>" +                                  "<password>" + dbPassword + "</password>" +                                 "<dbtype>" + dbType + "</dbtype>" +                                  "<namespace>" + dsName + "_NS</namespace>" +                                  exposePKsParameter +                               "</connectionParameters>" +                            "</dataStore>";                  return xml;              }      I'm using GeoServer 2.1.3 and PostGIS 1.5.


title: How to convert .geotiff --> .tiff?

text: How to convert .geotiff files to .tiff files (ordinary plain images)? Basically, I have cropped the image to a rectangle of interest and knowing its left-lower location and having chosen the pixel size, I could combine several data sources as ordinary .tiffs (and could finally forget the tiresome ArcGis world..) Any working method goes. We will employ usual image processing tools, and the location information is useless at that phase. Thank you for any help! =) UPDATE: This is what I got when I tried to read GeoTiff in GIMP 2.8.4: ![enter image description here](http://i.stack.imgur.com/3LLIu.png)


title: How can I hide a tool from command window under category

text: I have created a tool, a button and a dockable window under one project. To open/close the dockable window, button is used. For some specific work, by clicking some specific button from dockable window, I have activated the tool. So, basically the tool is used from the dockable window. So, the tool shouldn't visible to the user in the command window under the category. I couldn't find any document to hide the tool. Is there any process that can help me to hide the tool?


title: R : convert a spatial polygon objet to spatial polygon data frame

text: My goal is to modify an existing shapefile by merging certain polygons. After importing the shapefile and using the UnionSpatialPolygons command, I get the polygon outline that I want. However, this is now a SpatialPolygons object and not a SpatialPolygonsDataFrame, so I'm unable to export it to a shapefile using writeOGR. How can I get around this problem? Thanks in advance, John


title: Finding Ring Smaller than a certain tolerance

text: In Simplifying a shapefile for world timezones http://efele.net/maps/tz/world/ all generalize/simplify algorithms fail on 5 countries/states, shown in this shape file http://enclick.com/tz_world_handpicked.rar The failure of the simplifying scripts is due to small rings below the size of the tolerance used for simplifying. The entire multi-polygon is deleted for each of these countries How can I find and delete rings below a certain size with QGIS?


title: What kind of coordinates are this ones?

text: I have to map a point using this coordinates: 38 R MF 18852 15799 I can't understand the coordinate system. For sure '38' is the UTM zone number and 'R' is the designation (indicating southern Iraq)... but the rest?


title: Pipe Network: Computing the inflows for systems with multiple reservoir (input source)

text: Let's say, for the sake of simplicity, we have the following Y Shape water pipe network: ![](https://docs.google.com/drawings/pub?id=1lnPJkaP80nkzdHSB9YfHapth_tlhInnxF- UhPFhQN3s&w=960&h=720) `S1` and `S2` are the reservoirs ( input source), `Sout` is the water demand where water is drawn out for commercial use or whatever use. We know the height of the reservoirs at `S1` and `S2` relative to the ground and we also know the height of the `Sout`. All `S1`, `S2` and `Sout` are in`m^3/s`. For the pipe, we know for each pipe their corresponding diameter `d`, length `L`, roughness coefficient `C`. My question is, given the above information, is there anyway we can uniquely compute `S1` and `S2` for the Y shape network above? From what I know, we have one equation but two unknowns ( the conservation of source): `Sout=S1+S2` But this is not sufficient for us to uniquely determine `S1` and `S2` unless we have other equations here. I think the height of the nodes are also playing a role in the solution, but I just don't know how it can be done. Any ideas? On a side note, if our water pipe network is a loop, we can use the hazen william formula and hardy cross method to compute the head loss. This means that whatever that is required to be known about the pipe is already known, as long as they are used in hazen william formula and hardy cross method.


title: Can we merge shapefiles with different geometries (point and polygon) to get a single layer which shows both geometries?

text: I want to merge two shape files in which one is of point geometry type and other is polygon type in order to get a single layer containing both points and polygon in a single map.


title: Speed up pgRouting point_to_id using new <-> operator in PostGIS 2.0.0

text: The point_to_id function, which assigns a node id to a point based on existing nodes and tolerance, is using a query as follows:               SELECT     id, the_geom     FROM vertices_tmp WHERE         ST_DWithin(_p,the_geom, tolerance)      ORDER BY ST_Distance(_p,the_geom) LIMIT 1;      This example is very similar to the one used to explain the new KNN index based queries. I modified the query like this:               SELECT     id, the_geom     FROM vertices_tmp WHERE         ST_DWithin(_p,the_geom, tolerance)      ORDER BY _p <-> the_geom LIMIT 1;      However, as the documentation for '<->' says, the index kicks in only if one of the two operands is constant, as in ST_MakePoint(10,10). Is there any trick to get the index to be used also in a function?


title: Geodesic buffers for polygons

text: How can I create geodesic buffers around polygons? I'm on 9.3.1 ArcMap, but I need a programatic solution. IBUfferConstruction only does geodesic buffers for points. What are my options?


title: Records in related table do not appear when using the 'Related Tables' button in attribute table

text: ArcMap 10 - Geodatabase with related tables & relationship classes. I've written a script that parses data between two databases. When I go into the attribute table of the point feature class and select a record, I have the option to click the "Related Tables" button. When I use this button, the resulting table is either empty, or has only 1 record. The original data has many related records, and they have been parsed correctly into separate feature classes. The records are ALL in the related table, but they're not able to link for some reason. The relationship is through the Date/Time stamp, and I've checked to make sure that the values are the same. I tried to troubleshoot by adding a join, and everything worked as expected. Obviously with Joins you can only have a 1-1 relationship as opposed to a 1-Many. ![enter image description here](http://i.stack.imgur.com/Je3Qx.png)


title: what format distanceTo function returns in OpenLayers?

text: i am using distanceTo function in OpenLayers Geometry. http://dev.openlayers.org/docs/files/OpenLayers/Geometry/Point- js.html#OpenLayers.Geometry.Point.distanceTo here am confused in which format it will return the value. It is in Kilometers or meters?


title: How to set interactivity on multiple columns in CartoDB?

text: I'm trying to create a sublayer with interactivity on it for a tooltip. The tooltip is going to need data from a column named plant_name. So I am setting the interactivity property as I create the sublayer, like this:                       sublayer = layer.createSubLayer(               sql: sql,               cartocss: css               interactivity: ["cartodb_id", "plant_name"]             )      The problem is I get an error saying "layer1: Postgis Plugin: ERROR: column "plant_name" does not exist." even though the column is in the table. What am I doing wrong?


title: How to dynamically disperse point labels in ArcMap, using image files and leader lines?

text: My GIS group needs to find, or somehow develop, a labeling tool for points in ArcMap. Ideally, it would work like this:   1. Each point would be labeled using an image file (jpg or png) so that a location, such as a McDonald's restaurant, would get labeled with thumbnail-sized image of the golden arches, etc.   2. A leader line would extend from the image to the point location on the map   3. The labels would position themselves as close to the points as possible, without overlapping each other.   4. When you zoom in or out on the map, the labels would reposition themselves, moving closer or farther from the points, in order to keep from overlapping. Basically, we're looking for a method to label points with graphics, instead of text, that we can then throw onto a tiled base map online. If we can use any open-source or non-propriety methods, Javascript libraries, HTML5 tools, etc, it would be best for the budget. But, if need be, we are open to any alternative software options We also have access to ArcGIS for Server, so we could publish a map service if there is a way to accomplish this in ArcMap for Desktop. We have been looking for some kind of way to do this in ArcMap for more than a year, without any luck. ESRI has a few articles on how to use PNG files as point symbols, but we're looking for a way to use PNG/JPG files as point labels. Maybe it's too much.


title: Videos (as headers) in CartoDB Infowindows

text: I need to know how to embed a video (YouTube or Vimeo) into an infowindow in CartoDB in the same way that an image appears as an image header. Potentially more complicated, when there is no video link I'd like an image to appear instead. Can anyone assist me (as a non-programmer) with the code for this and WHERE to insert it in the HTML? I have viewed the "How do you embed iFrame content into a CartoDB infoWindow?" question, but I cannot make that work for either YouTube or Vimeo videos. Also what I am asking for is for it to be a part of the "image header" in CartoDB which is not what that question answers. Thank you!


title: filtering GeoRSSLayer contents in arcgis javascript API?

text: I would like to make use of a public geoRSS feed with content from 10years (200+ items), however since it is too much information, i would like to show only a subset of this content in my map... I have been unable to find any reference of how to filter/query the content of an esri georsslayer with the JS API any pointers on how this might be done, would be much appreciated!


title: How to handle simultaneous ArcGIS flex API GeometryService calls?

text: In the code below, when responses come back from the server for the b and c calls, I have the impression that the resposes are randomly given to b and c. The geometry service doesn't know which object is the owner of each response. Do I need to declare a GemetryService in B and in C?               Class A{         public static const SERVICE:GeometryService = new GeometryService(url);     }          Class B{         function callService(){             ...             SERVICE.buffer(geometries, responseHandler);            }               function responseHandler(event:GeometryServiceEvent){             ...         }     }          Class C{         function callService(){             ...             SERVICE.buffer(geometries, responseHandler);            }               function responseHandler(event:GeometryServiceEvent){             ...         }     }          Class D{         function method(){             var c:C = new C();             c.callService();                  var b:B = new B();             b.callService();         }      }


title: ArcGIS - Extract by Attibutes in Model Builder

text: I have a model which at one point uses the Extract by Attributes tool to produce a new raster with just the areas containing the maximum values of the first. I would think that the SQL expression would go something like: "VALUE" = (SELECT MAX("VALUE") FROM in_raster) No matter what combination I try of expression it always comes up invalid. I always have to run the model up to the point of this tool, then GET UNIQUE VALUES by hand in the sql builder window to input something like "VALUE" = some_number Then it works fine. The rest of the model runs and I get the correct output. I did try with Zonal stats but the output at the end of the model is not correct when I do. Is there a way of automating this perhaps with some python or another set of tools? I have a number of initial values I need to iterate through for this model and I want the whole thing to be able to run without me having to do that step manually. Thanks in advance


title: matching open layers plugin Google earth to my data

text: I'm novice to the GIS and am trying to figure out how to keep my datapoints matching to the same projection on Google Earth satellite that I downloaded through the Open Layers Plugin. My data is in a gpx format, downloaded from my garmin. When I first open the google earth image the location of my tracks make sense, until I zoom in our out and it changes completely. I understand this is a projection problem, but cant figure out how to fix - ive tried different projections. My current projection is WGS 84 If anyone can help me with this it would be much appreciated!


title: Loop - convert layer to raster coming up blank

text: I am looping through files, selecting out building footprints one by one and converting the selected building footprint to a raster. The script worked fine for about 2 hours this morning then stopped working. Suddenly all the rasters started coming up blank and I can't figure out why. The script says its running successfully. Thoughts?               import arcpy, os, shutil     from arcpy import env     from arcpy.sa import*          #set workspaces     arcpy.env.workspace = "C:\\ArcGIS\\SurfCityData\\Data" #arcpy.GetParameterAsText(0)     arcpy.env.overwriteOutput = True     outputWorkspace = "C:\\ArcGIS\\SurfCityData\\Data\\Scratch"             # Check out the ArcGIS 3D Analyst extension license     arcpy.CheckOutExtension("3D")     arcpy.CheckOutExtension("Spatial")          #Variables     Year = "2005" #arcpy.GetParameterAsText()     Elevation = "C:\\ArcGIS\\SurfCityData\\Data\\LiDAR_DEM_2005\\mosiac_2005"          ObsPts = "C:\\ArcGIS\\SurfCityData\\Data\\test_pts2005_prj.shp"      footprint =  "C:\\ArcGIS\\SurfCityData\\Data\\Footprints\\BuildingFootprints_SC.shp"     Ocean = "C:\\ArcGIS\\SurfCityData\\Data\\AtlanticOcean.shp" #arcpy.GetParameterAsText()          #Set analysis extent to elevation raster     arcpy.env.extent = Elevation     #Create Featue Layers     PointsFL = outName(ObsPts,"_Layer")     footprintFL = outName(footprint,"_FL")     arcpy.MakeFeatureLayer_management(ObsPts, PointsFL)     arcpy.MakeFeatureLayer_management(footprint,footprintFL)          #Select observation points one by one     RangeCount = int(arcpy.GetCount_management(PointsFL).getOutput(0))     for points in range (0,RangeCount):         FID = "FID=%s" % (points)             #Change elevation of building footprint to zero for a selected obs. point             arcpy.SelectLayerByAttribute_management(PointsFL,"NEW_SELECTION",FID)        print "\n", arcpy.GetMessages()        arcpy.SelectLayerByLocation_management(footprintFL,"INTERSECT",PointsFL,"","NEW_SELECTION")        print "\n", arcpy.GetMessages()        RastFootprint = outputWorkspace+"\\fp_"+str(FID)[4:].split(".")[0]        arcpy.PolygonToRaster_conversion(footprintFL,"FID",RastFootprint,"MAXIMUM_AREA")        print "\n", arcpy.GetMessages()      I haved tried coping the selected footprint feature to a new feature class and then converting that to a raster as well. The footprint copied to a feautre but still would not produce a raster even though the script says is running fine. I've tried this on two different machines, a personal laptop and a work machine with the same results. I should also add that I'm testing the script on a set of 4 points so nothing that should be bogging down my system.


title: Extract Raster from Raster using Polygon shapefile in R

text: I am new to R and using the raster package. I have a problem extracting polygons from an existing raster file. If I use               extract(raster, poly_shape)      function on the raster it always creates a list with the data. What I really want is to extract another raster file that I can load with ArcGIS again. After reading a bit more I think the crop function is what I really need. But when I try to use this function               crop(raster, poly_shape)      I get this Error:               Error in .local(x, y, ...) : extents do not overlap     In addition: Warning message:     In intersect(extent(x), extent(y)) : Objects do not overlap      The files raster and poly_shape are the same for both functions. Can you tell me what could be wrong here? Is it even right that the crop function creates another raster and not a list? Thanks for help!