title: Automating Geoserver & OpenLayers setup

text: I'm using OpenLayers, Geoserver, and PostGIS. For now, everything is fine and happy. I have 2 layers showing on my map through WMS from geoserver. I am taking following steps: Step 1: inserting data into postgis - done manually. step 2: if it is a feature in already present layer then reloading the page add that feature on the map, but for adding the table for a layer into geoserver, I am publishing it manually. step 3: including that layer code into JS manually I want everything to be done automatically. Since my project is to be handed over to a client and they should be given the possibilitz to add features, layers etc. I think automation is necessary. So I want: Step 1: inserting data into postgis - can be achieved by providing a form to user and using a php I can insert them into postgis Step 2: publishing of layer(selecting a style also) in geoserver - has to be done programmatically Step 3: inserting this layer code in JS - to be done problematically. I request your good guidance on this issue. Thank you.


title: How to apply custom OSM map style in OpenLayers?

text: I have a small web mapping application based on PostGIS/GeoServer/OpenLayers. I'm overlaying few WMS layers served by GeoServer and using OpenStreetMap as a background. In OpenLayers OSM is simply called as:               var osm = new OpenLayers.Layer.OSM("OSM");       map.addLayer(osm);      How could I display OSM basemap in shades of grey to accentuate my WMS layers? Is there a way to use existing styles from CloudMade's Style Editor for that purpose?


title: How can I compare two point shapefiles that contain lat/longs and compute distance?

text: Greetings everyone, I have 2 sets of point shapefiles in ArcGIS (group A and B), each set has around 1000 entries that contain lat/long. What I want to do is compare Group A against Group B and figure out for every point in Group A what point in Group B is the closest and what that distance is. What would be the easiest way to do this in arcgis?


title: What tool was used to produce these Cross Sections?

text: To start with the actual question: **What software was used to createthese cross sections? (PDF, aprox. 10MB)** Background: I have been really struggling to find a good software package that can generate professional cross sections. Most software focuses on 3D visualization but analysts want to see actual cross sections. I've already asked a similar question here but this time my requirements are more complex. Where before I only cared about visualization of straight-line boreholes now I need to visualize mining drillhole survey data within section slices of variable thickness. These drillholes are mostly inclined and curved, often slightly spiraled in attempts to intersect subsurface mineralized zones. Basic Parameters:   * Import drillhole and visualize drillhole surveys & data; collars, surveys, assays, lithology, etc.   * Must support deviated and inclined drill holes   * Create subsurface geology models from drillhole data and display on sections   * Import and visualize 3D DXF and other 3D (Leapfrog) subsurface models in sectional view   * Define section slices (variable thickness) in plan view    * able to view sections from various directions (Looking North, South, any azimuth)   * Edit completed section output (adjust labels, colors, layout)   * 3D view and navigation is **NOT** necessary It seems that the above linked sections satisfy all or most of these parameters. I would really appreciate If someone can suggest software that has these capabilities. I've already tried and dismissed: Gemcom GEMS, Geosoft Target, CrossView for ArcGIS Here is an example from the linked document. _(Detour Gold Corporation):_ ![enter image description here](http://i.stack.imgur.com/ze4DQ.jpg)


title: Does Leaflet.js uses WebGL for drawing on HTML5 canvas?

text: If someone could explain how does the Leaflet.js library draws object (and map layers) on HTML5 canvas? Does it uses WebGL and its GPU accelerated usage of physics and image processing in any way?


title: How to use geoserver to display data on hosted website

text: I am using geoserver to display GIS data on my hosted website when geoserver is running I can see my data on my website in both firefox and Internet explorer but only on my laptop where I have geoserver running, when I try to view the website using Firefox or IE on another PC I just get a blank box with error on page warning on the bottom of the web browser. can anyone help with this thanks


title: How to do a search in a line to a specific point

text: I'm using QGIS 2.0.1 Dufour. I'm working with a shape file made by a single line that represents a road with 30km. In that line I want to attach pdf files containing information about some points in that road. I would like to be able to do a search like: 'km = 2,300' and find the files with the information about that point. Is that possible?


title: Drillhole Plan Trace in Arc

text: I'm struggling to find the math to do what I want. I want to plot a top down view of a drillhole in ArcMAP 10.0. I have the beginning UTM coordinate and then for every place in the hole that was measured, a dip (in negative degrees, horizontal being 0), a direction (in degrees, magnetic north being 0, although I can make this true north), and a length or depth. I would just like to figure out the utm coordinates of each measured station so I can make a polyline. This isn't math I know and although I'm finding lots of information about 3d vector math, I haven't been able to figure out how to translate it to my puproses. So even if someone can tell me how to look for the formulas (what would it be called?) it would be much appreciated. Once I know the formulas, I can use it in Python. Thanks!


title: Summarizing data on ArcGIS 10.1

text: I am relatively new to GIS and am facing some difficulties trying to summarize/pool data together. I've attached the file to make it easier. FID_pc11_1 (Refers to postcode identifier) and count refers to the number of each postcode. Basically, there are 2107 objectid_1 (which is the different buffer zones around schools), code represents the different lower super output areas in UK. I would like to summarize my data, such that it would give me data in the format of; in zone 1 (eg.objectid_1= 1),what is the summated count (of postcodes) for each of the code/lower super output layer (eg:E101002711).![enter image description here](http://i.stack.imgur.com/oyCwO.png)


title: Reprojecting a base map in ArcGIS

text: I am new to ArcGIS. I have the EDN license (Desktop and Server); have setup ArcGIS for Server in our offline environment; have published several base maps (with caches that ESRI provided me on a DVD) as services; and can use those services with an ArcGisTiledMapServiceLayer from a Silverlight client we are building. The map displays, I can add our graphics layers to it, and everything works great. But, I do not want to use the default projection that these base map caches have (the 102100 WKID as seen in this arcgisonline service http://services.arcgisonline.com/arcgis/rest/services/World_Street_Map/MapServer). I'd like to use the projection from here: http://services.arcgisonline.com/arcgis/rest/services/ESRI_StreetMap_World_2D/MapServer which is the WKID 4326. Given just the base map caches that I have, is there any way for me to reproject the map? Can the client do it on-the-fly? Or can I configure ArcGis for Server to do that reprojection dynamically? Do I need to rebuild the whole cache some how to the new projection? Or find a new cache entirely?


title: How to get Earth radius in pixels?

text: I need Earth radius in pixels for use in Mercator projection. Well let's say that our map's height is **1478 pixels** and distance between **0 pixel and 1478 pixel is 2.52 degree.** Radius is distance from the equator to pole so the formula for count Earth radius in pixels is **R=90/2.52*1478** which is **52785.71429**. Well then pixel number for **48.54 degree** is **51287.99461** and for **51.06** is pixel number **54885.98194** counted by Mercator so the **diference is 3597.987328** which is not correct because it should be those **1478 pixels**. I think that Mercator formulas are good because when i make inverses of pixel i get the degrees as in input so mistake can be only in Earth radius. Am I missed something, are my calculations right? Thanks in advance.


title: Millions of polygons: displaying layers in Leaflet based on location being viewed?

text: I'm trying to create a map of all parcels in a state. In my state, this is several million parcels. In their current format, they are encoded as topojson and divided up by county. I'd like to show only a certain county's parcels at a time, preferably when the user is zoomed in on that county. Is there a way to achieve this? I'm open to alternative map libraries / ideas.


title: Export arcgis data to xml all at once

text: Is there a way in ArcGIS10 to export all feat classes, tables, datasets all at once as XML? I know there is a batch command but still appears you have to select items. I am in process of writing a script but do not yet know what method/function in ArcGIS to use that would get all at once the FeatureClassList, TableList, and DatasetList instead of getting one by one and then manipulating. Any Help is appreciated! Thanks.


title: How to export a PostGIS raster to PNG?

text: I have a raster which is split into 30 tiles and I want to dump them all into PNGs. I've read the documentation http://postgis.net/docs/RT_ST_AsPNG.html but still don't know how to dump it well. If i use ST_Histogram, i get numbers, think the raster is well inserted,               SELECT band, (stats).*     FROM (SELECT rid, band, ST_Histogram(rast, band) As stats           FROM "corine-g100_06-v15" CROSS JOIN generate_series(1,3) As band           WHERE rid=2) As foo;


title: ITopologicalOperator.Intersect, how does this work?

text: I'm trying to intersect all features of a layer with another's but I'm getting the web ADF application stuck while executing the "ConstructUnion" function Here's the code I've written so far               // create a geometrybag for layer1     ESRI.ArcGIS.Geometry.IGeometry geometryBag = new ESRI.ArcGIS.Geometry.GeometryBagClass();     geometryBag.SpatialReference = m_firstLayerFirstfeatureIGeometry.SpatialReference;     ESRI.ArcGIS.Geometry.IGeometryCollection geometryCollection = (ESRI.ArcGIS.Geometry.IGeometryCollection)geometryBag;     // Cycle all features of first layer and add them to the collection     foreach (DataRow m_dr in m_firstLayerfeatureGraphics.Rows)     {         // Takes this feature's geometry         ESRI.ArcGIS.ADF.Web.Geometry.Geometry adfFeature =             m_firstLayerfeatureGraphics.GeometryFromRow(m_dr) as ESRI.ArcGIS.ADF.Web.Geometry.Geometry;         // convert to a igeometry object         ESRI.ArcGIS.Geometry.IGeometry m_tempfeatureIGeometry =         ESRI.ArcGIS.ADF.Web.DataSources.ArcGISServer.Local.Converter.ToIGeometry(adfFeature, ctx);         System.Diagnostics.Debug.WriteLine("m_firstLayerfeatureGraphics +1 ");         object missing = Type.Missing;         geometryCollection.AddGeometry(m_tempfeatureIGeometry, ref missing, ref         missing);     }     ESRI.ArcGIS.Geometry.ITopologicalOperator mtopo1 = new ESRI.ArcGIS.Geometry.PolygonClass();     // FOLLOWING LINE BLOCKS APPLICATION     mtopo1.ConstructUnion((ESRI.ArcGIS.Geometry.IEnumGeometry)geometryBag);      Any idea why?


title: What does 'building topology' mean for pgRouting?

text: From a couple of answers at gis.stackexchange.com and at the pgRouting mailing list I have found out that pgRouting 'builds topology'. What does that mean? If for example we import an .osm file into PostGIS and then apply all the pgRouting scripts, what does happen? I thought that pgRouting just creates a table of road links (edges) from the OSM nodes, but after following the FOSS4G workshop I see the following tables in my database:               List of relations      Schema |       Name       | Type  |  Owner        --------+------------------+-------+----------      public | classes          | table | postgres      public | geometry_columns | table | postgres      public | nodes            | table | postgres      public | spatial_ref_sys  | table | postgres      public | types            | table | postgres      public | vertices_tmp     | table | postgres      public | ways             | table | postgres      I know what geometry_columns and spatial_ref_sys serve for, but what about the other tables and how do they relate to each other? What is topology in pgRouting's parlance?


title: How to label raster cells in QGIS?

text: I try to show a small resolution raster file in the QGIS Print Manager tool. I would like to show a small label with each cell value (position - _circa_ middle of each cell). How I can do it? Is there any plugin which can help me to accomplish this task? My idea is to use gdal2xyz, load such file using CSV plugin and then save it as a vector file. Unfortunately it is much work, because I have many raster files. Another idea is to use gdal2xyz.py, then ogr2ogr to save it as vector file. Is there any way to batch this process? Is there any open-source program where I can do it much more easily?


title: How do I zoom to the overall extent of two (or more) selected features from two layers?

text: I am creating an application where the user picks a value from a combo box, the application then selects that feature and then selects the corresponding features from a second layer, (both selections are in a different colour). What I want to do is as soon as the selection is made, the map will zoom to the envelope of the selected features. I know how to zoom to 2 individual layers (using pEnvlope.union...), but i need to do it on features. Has anyone got any idea of the best way to do this? (I'm using ArcGIS 9.3.1 and VBA). Thanks


title: Open Layers Plugin and QGis 2.4

text: I have installed Qgis 2.4 with the stand alone installer and the OSgeo4w too on three different computers. Two have win7 64 and one vista 64. It seems that the Open Layers plugin does not work anymore. What I did was the following 1)Open a gpx file with gps tools 1.1)Once I converted the gpx to shp and delcare that I wanted the data in EPSG 4326 2)Tried to use Open layers pluin either with the gpx file or the shp (in a new project) This is what I get 1. Python Warning warning:C:/Users/Mapoteca/.qgis2/python/plugins\openlayers_plugin\openlayers_layer.py:146: DeprecationWarning: QgsMapCanvas.mapRenderer() is deprecated if rendererContext.painter().device().logicalDpiX() != int(self.iface.mapCanvas().mapRenderer().outputDpi()): traceback: File "C:/Users/Mapoteca/.qgis2/python/plugins\openlayers_plugin\openlayers_layer.py", line 109, in draw self.render(rendererContext)   1. Python Error (same message on another computer) An error has occured while executing Python code: Traceback (most recent call last): File "C:/Users/Mapo/.qgis2/python/plugins\openlayers_plugin\openlayers_layer.py", line 109, in draw self.render(rendererContext) File "C:/Users/Mapo/.qgis2/python/plugins\openlayers_plugin\openlayers_layer.py", line 158, in render for res in self.resolutions(): TypeError: 'NoneType' object is not iterable Python version: 2.7.5 (default, May 15 2013, 22:44:16) [MSC v.1500 64 bit (AMD64)] QGIS version: 2.4.0-Chugiak Chugiak, 8fdd08a Morover I uninstalled the Open layers Plugin, since I wanted to reinstallit, but it is no longer on the list. My question is if the plugin is not compatible with Qgis new version or this behavior is just happening on my machines. Ps. I have a 2.3 version and it shows the same symptoms. Thanks before hand Gerardo


title: ESRI JSAPI: How to use GMap3 as 1st layer when using gmaps-utility-gis.googlecode.com?

text: This question relates to ESRI JSAPI and Google Maps V3 integration... I've put together this fiddle http://jsfiddle.net/br4pk/1/ which demonstrates a mismatch between the LOD of the map and gmaps.GoogleMapsLayer instance. The map operates great until you introduce layers which specify an `lod[0]` with a `level:0`. I believe this is because the gmaps layer `lod[0]` looks like this:               {         "level": 1,          "resolution": 78271.5169639999,          "scale": 295828763.795777     }      and the map lod[0] like this:               {         "level": 0,          "resolution": 78271.5169639999,          "scale": 295828763.795777     }      Should the map lod contain a level:0 when the basemap does not specify a level:0 lod?   What is the solution to getting this to work with other layers which do specify a level:0? I did find some success in modifying the lods of the gmap by adding a level 0 detail, and perhaps that is the answer (to modify http://gmaps-utility- gis.googlecode.com/svn/tags/gmapslayer/1.0/src/gmapslayer.js):                   { "level": 0, "resolution": 156543.033928, "scale": 591657527.591555 },         { "level": 1, "resolution": 78271.5169639999, "scale": 295828763.795777 },         { "level": 2, "resolution": 39135.7584820001, "scale": 147914381.897889 },         { "level": 3, "resolution": 19567.8792409999, "scale": 73957190.948944 },         { "level": 4, "resolution": 9783.93962049996, "scale": 36978595.474472 },         { "level": 5, "resolution": 4891.96981024998, "scale": 18489297.737236 },         { "level": 6, "resolution": 2445.98490512499, "scale": 9244648.868618 },         { "level": 7, "resolution": 1222.99245256249, "scale": 4622324.434309 },         { "level": 8, "resolution": 611.49622628138, "scale": 2311162.217155 },         { "level": 9, "resolution": 305.748113140558, "scale": 1155581.108577 },         { "level": 10, "resolution": 152.874056570411, "scale": 577790.554289 },         { "level": 11, "resolution": 76.4370282850732, "scale": 288895.277144 },         { "level": 12, "resolution": 38.2185141425366, "scale": 144447.638572 },         { "level": 13, "resolution": 19.1092570712683, "scale": 72223.819286 },         { "level": 14, "resolution": 9.55462853563415, "scale": 36111.909643 },         { "level": 15, "resolution": 4.77731426794937, "scale": 18055.954822 },         { "level": 16, "resolution": 2.38865713397468, "scale": 9027.977411 },         { "level": 17, "resolution": 1.19432856685505, "scale": 4513.988705 },         { "level": 18, "resolution": 0.597164283559817, "scale": 2256.994353 },         { "level": 19, "resolution": 0.298582141647617, "scale": 1128.497176 }


title: How to eliminate gaps in the ArcGIS union operation?

text: I have a layer that has nested polygons in them. I want to aggregate the polygons. I do this using the Union function. I use the "No Gaps" option when I perform the union. However, when I overlay my unioned layer on top of the original, I can see small gaps in the polygon. Has anyone encountered this before? We also thought that doing a topology would work better than the union. Except when we overlay the polygon count is different. What causes these gaps? And how do I remove the gaps? Thanks!


title: Mapserver tiles - trouble with projection

text: I'm trying to render tiles with **shp2img** tool from Mapserver suite and my tiles aren't matching. Tiles are numbered like in openstreetmaps.org. This is an attempt to draw tile (level 1, x=1, y=0) which has extent [xmin=0, ymin=0, xmax=180, ymax=90]:   ![](http://img69.imageshack.us/img69/4112/z0x1y0.png)   This picture is rendered with extent [xmin=0, ymin=0, xmax=180, ymax=90]. See the problem? Ymin is clearly not zero here.   And this one with extent [xmin=0, ymin=0, xmax=90, ymax=90]:   ![](http://img834.imageshack.us/img834/5489/z0x1y0v2.png)   Now ymin looks zero as it should. It seems like Mapserver tries to keep aspect ratio and is extending rendered area to do that. But this ruins my tile service.   **How do I setup projection in my MAP file to get matching tiles?** Or what else to do to get tiles right. **Update:** this is how tile z=1, x=1, y=0 looks in openstreetmaps.org, I guess I need reproject my map in such projection?   ![](http://a.tile.openstreetmap.org/1/1/0.png)   Some details:   ogrinfo output:               Layer SRS WKT:     GEOGCS["Unknown",         DATUM["MIF 0",             SPHEROID["WGS_84_MAPINFO_Datum_0",6378137.01,298.257223563]],         PRIMEM["Greenwich",0],         UNIT["Degree",0.017453292519943295]]      Projection section from MAP file:               PROJECTION         "init=epsg:4326"     END


title: How can I generate a high resolution rendering of the globe?

text: I am building an ANSI E poster size map collection and would like to add a nice rendering of the globe to the maps indicating the approximate position of the region on the globe. The problem is that I need the globe rendering to be high quality, high resolution image. I think minimum 300 or 400DPI @ about 12" x 12" inches so that I can work with it in CorelDRAW where I will be creating the final product. I have Google Earth (Free), NASA World Wind and ArcGLobe but neither software is capable of producing a satisfactory result. All look OK but the resolution of the image is not a press quality. So far Google Earth seems best (although terrain relief would be better then the spotty imagery), followed by NASA World Wind and legging behind is ArcGlobe. (Given the price tag, one would expect ArcGlobe to deliver but I am fully aware of the ArcScene exporter problems so I am not very surprised.) If someone could recommend a virtual globe (free) program that I can use to do this that would be great. Other options are welcomed as well; the globe does not have to be dynamic but it has to focus on Canada, centered on Ontario. Here are the outputs from the above mentioned programs: Google Earth ![Google Earth](http://i.stack.imgur.com/VI4lg.jpg) NASA World Wind ![NASA World Wind](http://i.stack.imgur.com/VDZNi.jpg) ESRI ArcGlobe 10 ![ESRI ArcGlobe 10](http://i.stack.imgur.com/9LBEo.jpg)


title: Clicking on a polygon-opening a page/link

text: I have a many polygons on my map serving through WMS from geoserver. I want to implement the task such that when the user clicks on a particular polygon , he should be directed to a html page about that polygon. Like if I click on my CSE Dept I should be directed to my cse dept home page. Need help in achieving this... I have followed http://lyceum.massgis.state.ma.us/wiki/doku.php?id=generate_imagemap as but when I try "http://localhost:8080/geoserver/wms?service=WMS&version=1.1.0&request=GetMap&layers=IITKanpur:halls&styles=iitk_halllinks&bbox=80.226,26.504,80.235,26.512&width=512&height=455&srs=EPSG:4326&format=text/html" I am not getting anything. My iitk_halllinks sld is : `<Rule> <Name>href</Name> <PointSymbolizer><Graphic><Size>16</Size></Graphic></PointSymbolizer> <TextSymbolizer><Label><PropertyName>halls_link</PropertyName></Label></TextSymbolizer> </Rule> `


title: Database View in Oracle 9i with SDE 9.3.1

text: If I have a sql statement creating a view against ArcSDE 9.3.1 on Oracle 9i, is the best practice to create the view in ArcCatalog just creating the name (kind of shell of the view) then going in to PLSQL to modify the view? What are the best practices on this? Reading this.. http://webhelp.esri.com/arcgiSDEsktop/9.3/index.cfm?TopicName=Using_database_views ..I am still left without great advice. Thanks. Here is the view               CREATE OR REPLACE VIEW SERIES_DETAIL AS     SELECT LineLoop.LineName, Cont_StationSeries.EventID AS RouteID, Eng_StationSeries.EventID AS SeriesID, Eng_StationSeries.SeriesOrder, Eng_StationSeries.BeginStation, Eng_StationSeries.EndStation, Begin_Cont_ControlPoint.StationValue AS BeginMeasure, End_Cont_ControlPoint.StationValue AS EndMeasure          FROM (((((LINELOOP LineLoop INNER JOIN STATIONSERIES Cont_StationSeries ON LineLoop.EventID=Cont_StationSeries.LineLoopEventID) INNER JOIN STATIONSERIES Eng_StationSeries ON LineLoop.EventID=Eng_StationSeries.LineLoopEventID) INNER JOIN CONTROLPOINT Begin_Eng_ControlPoint ON (Eng_StationSeries.BeginStation=Begin_Eng_ControlPoint.StationValue) AND (Eng_StationSeries.EventID=Begin_Eng_ControlPoint.StationSeriesEventID)) INNER JOIN CONTROLPOINT End_Eng_ControlPoint ON (Eng_StationSeries.EndStation=End_Eng_ControlPoint.StationValue) AND (Eng_StationSeries.EventID=End_Eng_ControlPoint.StationSeriesEventID)) INNER JOIN CONTROLPOINT End_Cont_ControlPoint ON End_Eng_ControlPoint.GroupEventID=End_Cont_ControlPoint.GroupEventID) INNER JOIN CONTROLPOINT Begin_Cont_ControlPoint ON Begin_Eng_ControlPoint.GroupEventID=Begin_Cont_ControlPoint.GroupEventID          WHERE (((Cont_StationSeries.SubTypeCD)=1) AND ((Eng_StationSeries.SubTypeCD)=2) AND ((Begin_Eng_ControlPoint.SubTypeCD)=2) AND ((End_Eng_ControlPoint.SubTypeCD)=2) AND ((Begin_Cont_ControlPoint.SubTypeCD)=1) AND ((End_Cont_ControlPoint.SubTypeCD)=1))          ORDER BY LineLoop.LineName, Cont_StationSeries.EventID, Eng_StationSeries.SeriesOrder


title: how to do warping of shape file in geotools

text: can any body give me simple warping example in geotools. I have non georeference shape file with polygons and i need to georeference this shape file by using warp technique. please help me.


title: Using OpenLayers (with OpenStreetMaps) and PostGis for rendering data acquired from GPS tracking device?

text: I'm researching and planing to acquire data from GPS device and store them in postGIS database, then use geodjango to process/work them out and finally chart them on a map with OpenLayers. Will these technologies be enough for my plan, since I don't want my maps, I want to use OpenStreetMap and plot my points on an OpenLayers' vector layer. I also want to implement real time tracking on a map, I'll query the PostGIS database in realtime using AJAX and then recreate objects on a map. I just can't find any good examples to do this with OpenLayers. Can I use Jquery for this and combine it with OpenLayers since I'm good at using Jquery's AJAX methods. P.S. Please bear in mind that I'm still researcing these technologies, but OpenLayers/geodjango/postGIS seems like a great combo for my needs. I found literature for almost everything that I'm interested in these 3, only that I lack the examples of real time charting and refreshing objects on OpenLayers' vector layer. TIA


title: How to make a route event layer?

text: I'm using ArcGis 10.1 but I'm a beginner in this program. I want to make a route event layer based on a route I have already defined. My question is in the third line when the program ask me for an input event table. I read a couple of explanations about this tool but any of them explains what table is this. Do I have to create an excel table? If yes, how does it have to look like?


title: How does Nearest Neighbor Image Resampling work in ArcGIS?

text: I am trying to understand how nearest neighbor for the resampling of Image datasets works in ArcGIS. The output raster cell value is the value of the closest cell value in the input raster: ![enter image description here](http://i.stack.imgur.com/CbXZj.png) In this case the center of each output cell is the middle cell of each 3x3 input cells. what happens if they are all at the same distance? if the output has half the dimension of the input, the center of the output will have the same distance to the 4 closest neigh input cells ?               InRas1=6x6     OutRas=3x3      He gets the majority of cell values? **No** Or am I missing something here?


title: MODIS MOD09A1 processing

text: modis state flag for MOD09A1 data (Surface Reflectance 8-Day L3 Global 500m) is different for each 8 day composite. I have 10 years data how to extract the cloud information in batch? and is their any process to interpolate and fill the gaps after removing cloudy pixels? I need a urgent help


title: Best OpenSource Database to serve images on web?

text: I am using a mysql database where I store images as blobs. From here, I serve them on the web. Does anyone know if there is another opensource option that works better? I would need something to serve them faster maybe? If anyone had some experiences on this issue, please drop an opinion. Thank you,


title: How to split circles in 8 quadrants

text: I am using ArcGis 9.3 and I need to split a large number of circles in quadrants. The circles corresponde to a 200m buffer of a point layer and, for each point, I need to split the buffer in 8 quadrants. Is there any ArcGis tool or extension that does this for all points/buffers at the same time?


title: Is Clip (Data management) or Extract By Mask (Spatial Analyst) more efficient?

text: I'm trying to clip down an Orthophoto to a county boundary. I've been using the Extract By Mask Tool, but this process has been running for about 12 hours! It's gone through 2 cycles already, so I'm assuming this is each band. Can someone correct me if I'm wrong on that. Would the Clip in Data Management work better (faster)? How would the accuracy be? I've always tried to use Extract By Mask, but this is taking too long.


title: Openlayers dynamically label polygon segments

text: My requirement is to show dynamic label on polygon feature's segments with some offeset. The label text should have the rank of each segment. If it is a square which starts at top left corner it will look as                    1           -------       |      |      4|      |2       |      |       --------          3      Can anyone give me somoe pointers in achieving this using openlayers? Thanks, Vijay


title: Zoom to Shapefile attribute table row feature

text: I am working with the python code below in ArcGIS to zoom into a shapefile's attribute table row features without selected until the end of table one by one. I am trying to use this code but this one requires that a row is selected.               import arcpy          mxd = arcpy.mapping.MapDocument('CURRENT')          df = arcpy.mapping.ListDataFrames(mxd, "Layers") [0]          df.zoomToSelectedFeatures()          arcpy.RefreshActiveView()      Any Help? Thanks Thank you for your reply. lets say, I have mapindex shapefile which has a 10 record in the attribute table. I want to zoom into first row automatically, then save the ArcMap .mxd or export as pdf. then move to next row, and do the same process(zoom in, save and export) until the end of record(row 10) in mapindex's attribute table. Again Thank you


title: Why is there a small white space between tiles in Chrome?

text: I'm using openlayers to display a google base-map layer and it works fine on Firefox/IE/Safari. I'm having problems with Chrome as the tiles have a small white space between them shown in this image. ![enter image description here](http://i.stack.imgur.com/TlrpM.png) I found this discussion about invasive CSS which suggested the problem was due to > Your theme is "helpfully" adding margins to img tags. you can override this > with               .openlayers-map img{ margin:0px }      I'm unclear where I put this line. I read through the openlayers CSS file but I don't have a good idea what is going on, on where I should be editing the `img` tag.


title: Keeping package dependencies for QGIS master up-to-date

text: Not for the first time, QGIS has broken due to unsatisfied package dependencies. This is not surprising, as I'm using the 1.9 Master version which is released nightly, and the package dependencies seem to update rapidly. Most recently, for example, QGIS could not reinstall due to the following package dependencies.               Depends: libc6 (>= 2.14) but 2.15-0ubuntu10.3 is to be installed           Depends: libgcc1 (>= 1:4.1.1) but 1:4.6.3-1ubuntu5 is to be installed           Depends: libqt4-network (>= 4:4.5.3) but 4:4.8.1-0ubuntu4.3 is to be installed           Depends: libqt4-sql (>= 4:4.5.3) but 4:4.8.1-0ubuntu4.3 is to be installed           Depends: libqt4-svg (>= 4:4.5.3) but 4:4.8.1-0ubuntu4.3 is to be installed           Depends: libqt4-xml (>= 4:4.5.3) but 4:4.8.1-0ubuntu4.3 is to be installed           Depends: libqtcore4 (>= 4:4.8.0) but 4:4.8.1-0ubuntu4.3 is to be installed           Depends: libqtgui4 (>= 4:4.8.0) but 4:4.8.1-0ubuntu4.3 is to be installed           Depends: libqtwebkit4 (>= 2.2~2011week36) but 2.2.1-1ubuntu4 is to be installed           Depends: libspatialite3 (>= 3.0.0~beta20110817) but 3.1.0~rc2-2~precise1 is to be installed           Depends: libsqlite3-0 (>= 3.5.9) but 3.7.9-2ubuntu1.1 is to be installed           Depends: libstdc++6 (>= 4.6) but 4.6.3-1ubuntu5 is to be installed           Depends: qgis-providers (= 1.9.0+git20121208+7841975~precise1) but 1.9.0+git20121208+7841975~precise1 is to be installed           Depends: qgis-common (= 1.9.0+git20121208+7841975~precise1) but 1.9.0+git20121210+df2c1c5~precise1 is to be installed      The strange thing for me is that more up-to-date versions of the dependent packages already seem to be on my machine, based on the above warning. I know the simple answer to this question is **"use 1.8.0 if you cannot sort your dependencies out"** , and this is precisely what I've done for now - (see very useful answers here if you need to fix broken packages). However, I really like using some of QGIS's new features such as the new labelling engine and would also like to contribute to the development process by submitting up-to-date bug reports. So if anyone has faced similar issues and knows of a reliable way to keep abreast with all QGIS 1.9 dependencies, please let me know.


title: Why can't I edit polygons once they have been saved to MySQL?

text: I have a simple table with an id field and a polygon field in mySql. Using QGIS I can add new polygons and give them an id. Before saving them I can also edit them. But once a polygon has been saved it cannot be edited any more: every time I move it, it simply moves back to the original position when I save the edits. Added points disappear when saving. The same happens with a point field in a different table. What could be the problem?


title: Should the geometry column always be the last column of a postgis table?

text: When we need to add an attribute to a table (e.g. a table of 350,000 polygons), my instinct is to always recreate the table so that the geometry column is last. I think this is a best practice for optimal performance, but I can't find any source to back this up. All imported shapefiles have the geometry column last, but I don't know if this is for performance reasons or not. I believe column order is important for some database platforms, but I'm not sure it is for PostgreSQL. Does anyone have an opinion?


title: Looking for TSP with Target different from Source

text: I know the TSP (traveling sales person) in osm2po is beta but it works quite well. Is there a way to designate a final target location that is different from the initial source location. In other words, I would like to designate:   1. Start location   2. Several intermediate locations   3. End location (different from start location)   4. Determine the best way to begin at the Start location, travel to all the intermediate locations while hitting the End location last. Currently, I see the TSP only does a round trip from an initial start location.


title: Adding images to placemarks multiple times

text: I have a route that I travel that is saved by coordinates to a .csv file for use in GE. I have images (pictures) that I would like to add to placemarks I have along the route. Is there a script or macro I can run that would allow me to add the url of the images to the .csv file so that I could have the images appear inside the placemarks balloon without having to add the images individually to the balloon?


title: Add Google Basemap to ArcGIS

text: I am looking for a way to add Google Earth imagery to ArcMap 10. Over a year ago, a GIS.SE member mentioned Arc2Earth was the only way to accomplish this (04/2011 thread). Another somewhat convoluted approach is described here. Has there been any progress in bringing Google Earth imagery to ArcGIS?


title: I want create a crystal report in gis with writing code

text: I am working in ArcGIS 9.3.1 and I try to write code for creating a crystal report and work with it in VBA, but I recieved an error when the program arrived at this line of code:               Dim pProj As CRAXDDRT.Application     Dim pReport As CRAXDDRT.Report      The error is "User defined type not defined". Why does the program give this error? Please help me with this problem.


title: How do I calculate curvature using GIS co-ordinates

text: I have GIS data for which I would like to calculate where curvature is present. Anyone has done this before, please help.


title: Connecting to WMS - problem

text: I have a problem connecting to my WMS service. In ArcMap i can connect but what i see is only blank page (cannot see the layers). In Qgis it also connects with service but when i choose add layers i got this message: ![alt text](http://i.stack.imgur.com/LCssR.jpg) Some of my WMS properties: Geographic Coordinate System: GCS_WGS_1984 Datum: D_WGS_1984 Prime Meridian: Greenwich Angular Unit: Degree And i have additional question: Can i set the version of my WMS when creating in Arcgis Server ? Because what i can see is only possibility to set a WMS capability (default version is set up to 1.3.0). * * * I think now i know what is my problem. Look this screen from ArcInfo. Here we got (below) my WMS service and above we have got our governmental geoportal. As i can see that my WMS have no style how it should be display and i think that is why i get white page. ![alt text](http://i.stack.imgur.com/LCssR.jpg) Do you know how i can set up a styles of display for my layers? * * * I do not know if it is what you expect (response from fiddler), because i never used this software but this is what i get: HTTP/1.1 200 OK Server: Apache-Coyote/1.1 Content-Type: application/vnd.ogc.wms_xml Content-Length: 4807 Date: Fri, 20 Aug 2010 09:43:47 GMT WMS karpaty OGC:WMS application/vnd.ogc.wms_xml image/bmp image/jpeg image/tiff image/png image/gif image/svg+xml application/vnd.ogc.wms_xml text/xml text/html text/plain application/vnd.ogc.se_xml application/vnd.ogc.se_inimage application/vnd.ogc.se_blank Layers EPSG:4326 EPSG:32634 0 KEO_boundaries KEO_boundaries EPSG:4326 EPSG:32634 default KEO_boundaries image/png 1 KEO_cities KEO_cities EPSG:4326 EPSG:32634 default KEO_cities image/png * * * Ad Simon Or you mean this? ![alt text](http://i.stack.imgur.com/LCssR.jpg) Here is some problem...


title: How to georeference a LANDSAT scene in QGIS 1.8?

text: I tried to georeference a LANDSAT scene in QGIS 1.8 but I failed. My first question is Can I use geografical coordinates as x,y coordinates of control points? or I need projected coordinates (from Campo Inchauspe/argentina 4 for example)? I used the geographical coordinates (lat; long in dd mm ss) and I obtained a black small (less than 1 meter) image.


title: How to Append to rather than overwrite single table using ArcPy?

text: The code thus far works to create multiple output files, but I'd like to set it up so that the data is all exported onto a single table across multiple rows. The issue is I believe with the part in bold below:               #script to calculate home range overlap (hectares) between two individuals          # import modules     import arcpy     import time     import os     #import tabletoexcel     arcpy.ImportToolbox("xxx\Additional Conversion - Generic Tools.tbx")          from arcpy import env          env.overwriteOutput = True               # local variables          #hardcode test     ou1 = r"xxx\bbmmAnsPoly90.shp"     arcpy.env.workspace = r"xxx\ninety"     outputfolder = r"xxx\ninety\test"     #logfile = r"xxx\ansoverlap.txt"     #ou1 = arcpy.GetParameterAsText(0)     #arcpy.env.workspace = arcpy.GetParameterAsText(1)     #outputfolder = arcpy.GetParameterAsText(2)     #int_output_Statistics = arcpy.GetParameterAsText(3)     logfile = outputfolder + "\\" + ou1 + "_" + "overlaps.txt"     print logfile          try:         # process: intersect         #intersect_analysis (in_features, out_feature_class, {join_attributes}, {cluster_tolerance}, {output_type})         ou2s = arcpy.ListFeatureClasses()              for ou2 in ou2s:               fcname1 = arcpy.Describe(ou1).basename             fcname2 = arcpy.Describe(ou2).basename             **int_output = os.path.join(outputfolder, str(fcname1 + fcname2) + ".shp")             #int_output = os.path.join(outputfolder, str(fcname1) + "_overlaps.shp")             int_output_Statistics = r"xxx\Default.gdb" + "\\" + str(fcname1 + fcname2)             #int_output_Statistics = r"xxx\Default.gdb" + "\\" + str(fcname1) + "_overlaps"**             excel_output = outputfolder + "\\" + str(fcname1 + fcname2)               #excel_output = outputfolder + "\\" + str(fcname1) + "_overlaps"                  logfile = outputfolder + "\\" + fcname1 + "_" + "overlaps.txt"             if arcpy.Exists(logfile):                 f = open(logfile, "a")             else:                 f = open(logfile, "w")                      f.write(time.strftime('%c'))             f.write("\n")             f.write("------------------------------")             f.write("\n")                  arcpy.analysis.Intersect([ou1, ou2], int_output)             f.write("\n" + "Calculating home range overlap between " + fcname1 + " & " + fcname2)             print "Calculating home range overlap between " + fcname1 + " & " + fcname2                  # process: add field             arcpy.management.AddField(int_output, "Overlap", "DOUBLE", "", "", "", "", "NULLABLE", "NON_REQUIRED", "")             f.write("\n" + "Adding \"Overlap\" field to " + int_output)             print "Adding \"Overlap\" field to " + int_output                  # process: calculate field             arcpy.management.CalculateField(int_output, "Overlap", "!shape.area@HECTARES!", "PYTHON_9.3", "")             f.write("\n" + "Calculating geometry of overlap between " + fcname1 + " " + fcname2)             print "Calculating geometry of overlap between " + fcname1 + " " + fcname2                  # process: delete identical             arcpy.management.DeleteIdentical(int_output, "Overlap", "", "0")             f.write("\n" + "Deleting identical values from \"Overlap\" field")             print "Deleting identical values from \"Overlap\" field"                  # process: summary statistics             f.write("\n" + "Determine overlap value in hectares")             print "Determine overlap value in hectares"                  arcpy.analysis.Statistics(int_output, int_output_Statistics, "Overlap SUM", "")             arcpy.management.AddField(int_output_Statistics, 'orangutan1', 'TEXT')             arcpy.management.AddField(int_output_Statistics, 'isopleth1', 'DOUBLE')             arcpy.management.AddField(int_output_Statistics, 'orangutan2', 'TEXT')             arcpy.management.AddField(int_output_Statistics, 'isopleth2', 'DOUBLE')                  search1 = arcpy.SearchCursor(ou1)             for scrap in search1:                 ou1_name = scrap.getValue("ou")                 ou1_iso = scrap.getValue("iso")                  search2 = arcpy.SearchCursor(ou2)             for scrap in search2:                 ou2_name = scrap.getValue("ou")                 ou2_iso = scrap.getValue("iso")                 #row_values = [(scrap.getValue("ou"), scrap.getValue("iso"))]                  cursor = arcpy.UpdateCursor(int_output_Statistics)                  for row in cursor:                 row.setValue("orangutan1", ou1_name)                 row.setValue("orangutan2", ou2_name)                 row.setValue("isopleth1", ou1_iso)                 row.setValue("isopleth2", ou2_iso)                 cursor.updateRow(row)                       f.write("\n" + "Writing to:" + " " + int_output_Statistics)                  #export_to_xls(dataset, output):             arcpy.TableToExcel_conversion2(int_output_Statistics, excel_output, "XLS")             print "Writing to excel file"             f.write("\n" + "Writing to excel file:" + " " + excel_output)             f.write("\n" + "------------------------------")             f.write ("\n")                             f.close                    except Exception, e:         # If an error occurred, print line number and error message         import traceback, sys         tb = sys.exc_info()[2]         print "Line %i" % tb.tb_lineno         print e.message      So, is there a method (perhaps similar to f = open(logfile, "a" ?), whereby I could append to rather than overwrite a table?


title: convert GPS to XY coordinates

text: if i have obtain GPS coordinates as Latitude/Longitude, then i need to convert them to XY coordinates, i see this link actually i need to get the Trilateration , in its equation i need the XY coordinates from GPS coordinates. i need obvious and clear equation for conversion. also how the conversion will be if i assume there is reference node !! my application is : i need to make Trilateration for sensor nodes and in the middle of the there is pc which get positions of nodes(depends on gps) then pc can calculate its position depend on Trilateration. So Trilateration equation use XY rather than lat/lon ?? can we use the pc as reference node to conversion.


title: ArcGIS opposite to pivot table

text: fellow GIS-ers :)!   I've been trying to find a tool to do the opposite to what the "pivot table" tool does but I dont thin there is any tools for that. I have a raster data set derived from Landsat imagery with a only a single value per pixel. In Arcmap the table for this file is organize in a way that it has two fields (apart from the ID field). 1. count, and 2. value, or in other words all of the one-to-one relations were joined into one-to-many relations. i need to find a way to go back to the one-to-one relations "form". Which by looking at how a pivot table tool works is kind of the opposite. please any help would be much appreciated.                            **how it is**                                **how i need it**          ID    Count     Value                         ID      Value          1      3          75                          1        75          2      4          65                          2        75                                                        3        75                                                        4        65                                                        5        65                                                        6        65                                                        7        65      Fabs


title: What do ArcGIS regression analysis warnings 000853, 000916 & 000981 mean?

text: I'm trying to do regression analysis using a shapefile which contains dependent variable (NDVI) and explanatory variable (prec) in ArcGIS. First, I run GWR analysis and then run spatial autocorrelation Moran I on residual of GWR result. Finally, I wanted to see summary statistics on result window in ArcGIS but it just gives me 3 warnings related with analysis. Do you have any ideas about why these warnings appear? The shapefile has already been projected into Geographical lat/long which makes me wonder whether the result of regression is correct or not. ![result of GWR shown in result window](http://i.stack.imgur.com/RCoG0.jpg)


title: QGIS rounds up the coordinates to six decimals when saving layer changes

text: I have the problem that when I create a new point (in a point layer) or a new line (in a line layer of course) and I save the changes to the layer (so NOT the project save) then QGIS rounds up the coordinates of the new element to six decimals. I can see that in the "identify results" window by clicking on "(Derived"). So as a result, I don't get the required position of the new element. The layers are DB tables that I received and I just open them in QGIS. I don't know the source of the layers. The CRS of the project is the same as the CRS of the layer. The thing is, if I create a new shapefile-layer myself (point for example), the new element remains in the wanted location even after saving layer changes. So it has something to do with the layer I received. But I don't know what that is. I am trying for 3 days to find a solution and can't seem to solve the problem. Does anyone have a hint for me? Bevore saving the layer changes: ![bevore save](http://i.stack.imgur.com/gNuhN.jpg) After saving the layer changes: ![after save](http://i.stack.imgur.com/NXGiT.jpg)


title: How to automatically create a map for each column of attribute table?

text: I have an atribute table with yearly data for 10 regions, is there a way how to automaticaly create a set of maps (one for each year) with fixed scale? I tried to search the web but perhaps I didnt ask the correct question Thanks a lot for any help Jakub Picture - how I would like the result to look like ![result](http://i.stack.imgur.com/zJlbI.jpg)


title: How to view the whole KML

text: When I load a KML in Google maps, it doesn't show all polygons, instead there are Google "next-previous" pages, each containing parts of the KML. How can I view all the information in one page? Thanks, Tamir


title: Extracting points bordering shape extent

text: I'm using ArcGis 9.31.   I have a set of points as elevation and in order to get a soft clip area for the TIN, I would like to extract only all the points bording the extent of this shape, or to get a shape deliminating the extent of the set of points. I don't have any other data than this point shape. Is there a way to select them? Many thanks!


title: Is there a workaround for GRASS i.maxlik 'class' being a keyword in Python?

text: I am running i.maxlik in python - attempting to automate a processing chain. I use the command:               grass.run_command('i.maxlik', group= group, subgroup=subgroup, class='classification', sigfile = sigfile, verbose = True)      Infuriatingly, however **'class'** is a python keyword and it doesn't run due to **'invalid syntax'.** Would you know a way to work around this? I have tried emailing the GRASS developer list but no response. Regards Becky


title: Using Google Maps Offline in QGIS

text: I am working on a QGIS project that involves me using the Google Hybrid Map Plug-in. I am due to spend time travelling over the next few days and won't have an internet connection. I was hoping someone could let me know if it's possible to continue using the plug-in without internet? For example by saving the relatively small I need for offline use. Thanks Sarah


title: How to tell if GeoWebCache REST API is working?

text: In my web application, i wanted a service to cache updated data using geowebche. so i worked on seeding tiles automatically with GeoWebCache as in this post: custom time interval for geowebcache then i replaced `wms` with `gwc/service/wms` in layer requesting in OpenLayers. i didn't do anything more than what's in the post. here it said i have to install the extinction but i didn't. http://docs.geoserver.org/2.0.0/user/extensions/rest/installation.html#rest- installation Now, how can i tell if it works? do the tiles get refreshed by themselves without me refreshing the page or what? and I use WMS to display my points on the map, is there any way to cache the layer where it doesn't need to reload points whenever i navigate on the map?


title: Mapping the European Debt Crisis

text: I'm a journalist based in Paris, and for a project I'm working on, I would like to display how public debt rose dramatically across Europe during the crisis. My plan is to show each country's ratio of debt/GDP between the years 2000 and 2010 inclusive. Any country with a ratio under 50% will be shown in green, anything from 50-60% in lighter green or yellow, etc. I'm just not sure of the most effective way to portray this, using Cartodb. I'm really new to the software and to geographic visuals in general. Thanks for your help!


title: Which edition of ArcGIS for Server and which capacity level?

text: I need to define the edition of ArcGIS for Server product to invest for a small scale oil & gas company. What we need is to be able to:   1. Work on a geodatabase with 2-3 concurrent desktop users while a allowing access to a few more desktop users,   2. Serve maps from geodatabase (and satellite images and documents from the file server) to internal non-GIS users through web browsers and enabling specific tools like querying and basic geoprocessing on the client side.   3. Allow some users feed the geodatabase using their tablets or smartphones, I read the ArcGIS for Server functionality matrix and got confused. The first part of the document tells me that I need ArcGIS for Server Enterprise Advanced (full product) but the case studies made me think that a workgroup- standard license seems to be enough. Similarly I first got the impression that I need the Portal for ArcGIS product to serve any maps (do not want to use cloud) but then I thought ArcGIS for Server alone can be enough for services like WMS, WFS etc.. My company can invest on these products so I do not consider Open Source platforms but just want to make sure that we are not over investing. Any advise on this?


title: Why do the values calculated by pyproj and the ones by a website differ greatly?

text: By this website http://www.movable-type.co.uk/scripts/latlong.html, I would like to calculate distances between latitude 37.919 and 37.465, longitude 15.358 and 14.601, the results are 50.48 km and 66.81 km, respectively. e.g. Point1: 37.919, 14.601 and Point2: 37.465, 14.601, the website gives 50.48 km, Point1: 37.465, 15.358 and Point2: 37.465, 14.601, I got the distance among them, 66.81 km. However, I use the pyproj.Proj to calculate the distances, which differs greatly from the results above. Below is my python code:               import pyproj.Proj     p=Proj(proj='eqc',lat_0=37.692, lon_0=14.980)     p(14.601,37.465) #return (-42190.08701065082, -25269.524410073795)     p(15.358,37.919) #return (42078.76751985713, 25269.52441007167)      The distance between 37.919 and 37.465 is 25269.52441007167-(-25269.524410073795), roughly 50539 m = 50.539 km, which matches approximately the value by the website: 50.44 km But, the distance between 15.358 and 14.601 is roughly 84269 m, 84.269 km, which differs greatly from the one by the website: 66.81 km. I don't know why?


title: Applying different labeling options at different scales to different points (that aren't easily separated from the other points)

text: I am creating maps (in ArcMap 10.1) for a web-based application. There are more than 500 places to be labeled at some scale (the map scales range from 1:32 million to 1:1.5 million). As you can imagine, some of the places are very near one another and some are far away from other places. I would like those places that are farther away from other places (where there is no labeling conflict) to be labeled when the map is at a smaller scale than those for which zooming in to a larger scale is needed to avoid a cloud of overlapping names. (I don't mind them overlapping at the largest scales...) I am using the Maplex labeling but don't see how to easily accomplish this. At some point they may still overlap, so I cannot use the "never overlap" option... Thank you for any suggestions or guidance you can offer!. Sandy


title: Printing map to scale

text: How to print this map (page) to 1:500 without georeferencing or without using any gis software. It is better to say that the marked line of 3 km should be 6 cm on paper (A4). See image: ![enter image description here](http://i.stack.imgur.com/Ob5dR.png)


title: Summarizing different (local) shapefiles/sources into one centralized database using ArcGIS for Desktop?

text: Using ArcMap/ArcCatalog 10.2, is there a way, to summarize a huge number of widely spread shapefiles/sources (different servers, different hard drives etc.) into one centralized database, in order to be able to transfer this database as well as a corresponding mxd-file e.g. to a notebook or another computer without internet connection?


title: How to batch create filenames in ModelBuilder?

text: I am running a model in ModelBuilder. The input are raster files in *.tiff format and the output are *.asc format. I am using iterate raster and there will be 365 input and output files using %name%.asc e.g input file : LST01011990.tiff --> output file LST01011990.tiff.asc How I can output the filename to LST01011990.asc (without repeating the *.tiff)?


title: Obtain all elevations points from a raster DEM

text: I need to obtain all elevation grid points from a raster DEM. I see the steps required as:   1. Read in raster DEM (Arc Info Binary Grid format).   2. Create subset area using Points2One (see triangle area here: http://imgur.com/aeOPBBD).   3. Write subset-grid elevations to a file (I can't do this yet). I only need elevations (ie. Z only, not XY). **Any ideas for step 3?** I'm happy with getting all grid elevations or a histogram of elevations (with frequencies for (say) 100 bins). I will later manipulate the elevation data in EXCEL or R or Fortran or Python or ??? Note: I'm very new to using GISs and am using QGIS 1.8.0 on Windows (since I found that the easiest program to use). Thus, I'd prefer a solution using QGIS, but I'm willing to use another free GIS program if required. Also, I'll likely have to automate/program this extraction procedure later (ie. steps 2 and 3), but would prefer a GUI-based solution first for simplicity. Thanks, Tom.


title: qgis, autocad and google layers

text: Ive downloaded qgis and loaded the google layers plugin. I have imported a dxf cad drawing. Question is how do i place the shapefile in the correct place on the map? Also how do i tell google layers to return a particular are of the map? Thanks in advance


title: Which registry key holds the license type in ArcGIS 10.x?

text: I'm unable to change my license type (ArcView, ArcEditor or ArcInfo) in ArcGIS 10.0 using the Desktop Administrator due to permissions restrictions. I'm trying to use the back-door method of changing the registry key in Regedit. In older versions of ArcGIS this was found under `> HKEY_LOCAL_MACHINE > Software > ESRI > License` but I'm not seeing this in ArcGIS 10.0 The post Opening ArcGIS on a specific license level mentions a similar problem on Citrix, however I'm unable to set a system variable due to the same permissions issues. Which registry key holds the license type in ArcGIS 10.0, and is it still possible to change the license type using this method?


title: Migration from Oracle to PostgreSQL

text: I emphasise that I am a newbie in databases. Is it possible to migrate from Oracle database to Postgres? If yes do you know any tested solutions?


title: Why should a user choose to use open source software in their remote sensing research?

text: If you judge by the amount of the questions regarding FOSS software that are being asked in gis.SE many users seems to prefer FOSS over proprietary software. I've read some articles - more precisely some personal blogs - supporting this choice. Furthermore corporate giants like Esri seem to acknowledge the open source development movement. So, I am asking your position in this matter. Why do you use open source tools? What are the advantages or disadvantages of your choice, if any?


title: Why does OpenLayers load OSM-like tiles, but does not show them?

text: I created a tile server with OSM-data in the German style. The data is quite similar to the official tile server for openstreetmap.de, provided by Fraunhofer IOSB.               http://tile.iosb.fraunhofer.de/tiles/osmde/${z}/${x}/${y}.png      If I include the standard OSM-server               http://a.tile.openstreetmap.org/${z}/${x}/${y}.png      with the following constructor               new OpenLayers.Layer.OSM("OSM",myOSMURL,{...});      I get a valid map. The OpenCycleMap tile server works as well. If I want to use my own tile server or the one from Fraunhofer, the browser successfully gets an respond to the GET-request from the tile server and loads the tiles needed, but doesn't show them in the map. Following the steps here, I generated a map with local data successfully. Nevertheless, if I replace the link of the local data with a link to Fraunhofer or my tile server, the browser still loads the tiles, but doesn't show them.               // This is the layer that uses the locally stored tiles     //var newLayer = new OpenLayers.Layer.OSM("Local Tiles", "tiles/${z}/${x}/${y}.png", {numZoomLevels: 19, alpha: true, isBaseLayer: false});     var newLayer = new OpenLayers.Layer.OSM("Local Tiles", "http://tile.iosb.fraunhofer.de/tiles/osmde/${z}/${x}/${y}.png", {numZoomLevels: 19, alpha: true, isBaseLayer: false});     map.addLayer(newLayer);     // This is the end of the layer      ![enter image description here](http://i.stack.imgur.com/73jOX.jpg) Image in full resolution: http://i.stack.imgur.com/73jOX.jpg How do I get the tiles from my tile server or Fraunhofer to my map?


title: How to copy annotation feature and paste within same layer?

text: I am having trouble executing a very basic command which is to copy and paste an annotation within the same layer in ArcGIS 10.1. First, I created an annotation class, then created one label (a balloon call-out with a number) whose attributes I want to copy for all of my locations. When I went to Editor > Edit Annotation Tool > selected my annotation label > then under Attributes (in my Editor toolbar) > right-click my label > Copy Attributes > then right- click on my Layer > Paste Attributes, nothing happens. Can someone explain why I cannot paste attributes in this way? Thanks for your help.


title: only one popup displayed on the openlayers

text: The issue is that i only one popup on the map is dispalyed. I can see all the points but popup can see only on the last point of the list. Please check my code http://jsfiddle.net/3HSjY/24/ I'm stack with rhis problem from Sunday and can not find any info that could solve my problem on the internet :(


title: arcpy searchCursor numeric variable in where clause

text: I know this question has been answered but I am not seeing an adaptation that I can understand. I am trying to use a numeric variable in a where clause and just having a hard time figuring it out. Your help is much appreciated               fc1 = "L:\ENGINEERING\collind\ArcGIS\Default.gdb\sedaSampleHistory"          fieldname = "pocid"          delimfield = arcpy.AddFieldDelimiters(fc1, fieldname)          pocidCnt = 3          cursor = arcpy.da.SearchCursor(fc1, ["pocid", "cl2total"], delimfield + " = pocidCnt ")      I have tried all kinds of ways to get pocidCnt to work in the where clause and I am not having any luck. If I write it as               cursor = arcpy.da.SearchCursor(fc1, ["pocid", "cl2total"], delimfield + " = 3 ")      it works fine. What am I doing wrong?


title: Leaflet Draw - add title to marker or shape?

text: I am trying to create a participatory mapping site where users can draw pts, lines and polys using Leaflet draw. Is it possible to create a pop window where users can add a title or description of the feature they draw -or- after they draw, click on the feature and add a label or description? My mock-up is here: http://geointerest.frih.org/smap/ Thanks, GV


title: Trouble installing PostGIS PlugIn for OpenJump

text: I’m trying to install PostGISPlugin version 1.52 in OpenJUMP 1.7.1 release 4004. I moved the JAR file into my installation’s lib/ext directory, but I’m not sure how to access the plugin-specific UI components. I’m not even sure it’s been properly installed. Can anyone point me in the right direction?


title: How to show average walking distance to closest point in census tract in ArcMap?

text: Given a map layer of point locations, and a second map layer of census tract polygons for the same area, how would you create a map layer that shows the average walking distance to the closest point for each census tract? Euclidean distance?


title: How to insert a map into another map?

text: If you look at this site it has a main map window but it also has a small window map at the top right corner.   1. What is the technical term for the smaller window?   2. How is that window created?  I am currently using leaflet js and I wanted to insert that into my map here.


title: Defining areas for inclusion in polygon conversion using ArcGIS for Desktop?

text: Using ArcGIS 10.2 for Desktop, I am trying to create a polygon shape(s) from a raster image. Is there a way to indicate areas I do not want included in the shape. Similar to a buffer I suppose, but that will not do what I want it to do. I have the image below, top is an inset. I want the "main" areas to be within a polygon, but if a point is too far out it will be left out, this can be arbitrary depending on the location and position of other points. Is there any way to control this? I have looked at using Minimum Bounding Geometry, Convex hulls, etc. The main problem is I need to do this for over 200 files and aside from hand drawing a polygon around each shape, is there any better way to do it? ![Russia](http://i.stack.imgur.com/cW7L8.jpg) ![enter image description here](http://i.stack.imgur.com/MIEOI.jpg)


title: cadtool and spline plugins not working

text: Am trying to edit a line shapefile using a spline tool in QGIS 2.4.0 .Its drawing but the ok button is inactive. what could be the problem? Generally am trying to create a curved line.


title: Scatterplot, correlation of dependence and r-square between two rasters

text: I have two rasters of the same region and of same resolution and I want to plot in a scatter plot their values and compute the best - fit line with r-square. Any idea how can I do this on Arcmap or QGIS or extract the values of the rasters in a file so I can import it into Excel? Unfortunately I don't use R neither python.


title: Buffering algorithm

text: Does anyone know the algorithm of buffering? I had difficulty in drawing buffers. So I tried many values to acquire 1km buffer. I found value 0.008966 produces 1km buffer. Then I calculated inverse number. 111.5324 was the result. Circumference of the earth is 40077km. 1 rad is 57.29578. 40077*57.29578/360=6378.4527 6378.4527 was divided by the 111.5324. Answer was 57.18923 . I am sure there are errors in the calculation, but the final result was close to 1 rad. Is there any possibility that a squared value was used? Thank you.


title: Drawing outline of points in ArcGIS for Desktop

text: Simple question, can Arc take a shapefile of points, calculate which ones are at the "edges" and draw a polygon or polyline around it. I know there are masks and bounding geometry and all of that, but it just needs to draw a line around and can intersect the exterior points.


title: Join by attribute when unique identifier has multiple rows using ArcMap 10.2?

text: I would like to join the attributes of one polygon layer to another polygon layer by their common ID. The problem is that in layer A, each ID has multiple rows, because each represent shapes from the layer B cut in multiple pieces (intersected with other information) and I need that information to be replicated each time within the boundaries of the layer B. Is there a way to join by attribute, so that the target layer B repeats rows of the adjoining Layer A, when the ID in Layer A occurs multiple times?


title: Data structure: mongoengine to GeoJson

text: I recently posted my first question on StackOverflow, but I think the question was quite extensive and full of potential errors. So I'd begin again, step-by- step, with a much shorter and simplier question : **The purpose is to structure my models.py in my django website to store datas in a Mongo database asGeoJson objects.** Do you think this code is correct, I'm not very sure about the syntax ? Thanks a lot ! * * * ## **models.py :**               # -*- coding: utf-8 -*-          from mongoengine import *     connect(‘mongodb_jsons’)          import datetime          class GeoJson(Document):      # from a geojson object      # save it in Mongodb with a geojson structure              Date_created = DateTimeField(default=datetime.datetime.now)         Latitude = FloatField()         Longitude = FloatField()         Content1 = CharField()         Content2 = CharField()              meta = {'db_alias': 'mongodb_jsons', # save in DB ‘mongodb_jsons’                 'indexes': [                 # the geojson structure                     {'type' : 'Feature', {                        'geometry':{                             'type' : 'Point',                             'coordinates' : ['Longitude','Latitude']},                              },                         'properties':{                            'content1' : 'Content1',                            'content2' : 'Content2',                                'date_creation' : 'Date_created'                             }                         }                     }]                 }      * * *


title: How do i get country boundary and names to show up in postgis backed mapnik rendered map?

text: I followed this to setup postgis server with planet-latest.osm. I used generate_xml.py to generate stylesheets. Using this style sheet with mapnik produces a world map with no country boundaries or country name on it. What is the best way to render tiles that looks similar to official tiles.openstreetmap.org tiles. For example, for same zoom, x, y, this is what i get: ![enter image description here](http://i.stack.imgur.com/tl3jm.png) And this is what tile.openstreetmap.org has: ![enter image description here](http://i.stack.imgur.com/NLTS0.png)


title: What is the correct projection to use with tm_wgs84_sans_antarctica.shp?

text: I am using tm_wgs84_sans_antarctica.shp with mapnik to generate tiles. I am using the tiles with Google Maps Javascript API V3. But the projection is way off. Here is roughly what my code looks like:               import mapnik     m = mapnik.Map(           tile.size()[0],           tile.size()[0],           )     s = mapnik.Style()     r = mapnik.Rule()          polygon_symbolizer = mapnik.PolygonSymbolizer(mapnik.Color('#f2eff9'))     r.symbols.append(polygon_symbolizer)     s.rules.append(r)          m.append_style('My Style',s)          ds = mapnik.Shapefile(file='tm_wgs84_sans_antarctica.shp')     layer = mapnik.Layer('world')     layer.datasource = ds     layer.styles.append('My Style')     m.layers.append(layer)          bbox = tile.bounds()     bbox = mapnik.Box2d(bbox[0], bbox[1], bbox[2], bbox[3])     m.zoom_to_box(bbox)     im = mapnik.Image( *tile.size() )     mapnik.render(m, im)      As you can see I am not setting any srs either in map or in layer. I tried different combination of projection strings without really understanding what they mean, but either got no map or wrongly projected map. Please let me know if any other information is required. Thanks.


title: Print Library for Open Layers?

text: I realize this is a bit of an open question but are there any known open javascript libraries which facilitate printing an OpenLayers map? I'm looking to produce the content on the client and possibly upload it to convert to PDF or simply set the innerHTML on some dom element. I started a proof-of-concept but was looking for something which honors various tile sizes and vector markup (IE 8) and all the other complexities that are sure to arise.


title: How to convert Geostatistical analyst (GA) file to Raster FAST?

text: I have been trying to convert a number of Geostatistical analyst files that i were created using Geostatistical wizard in ArcGIS 10. I want to save these files as Raster and use the File Am i missing something, is there a way to do it faster ?


title: OpenLayers: How to get *real* centroid of a Line Geometry/Feature?

text: I have declared a line geometry as:               var lineString = new OpenLayers.Geometry.LineString([source,target]);      then I am trying to get the centroid of this lineString as `lineString.getCentroid()`. But what's returning is a Start Point of the lineString. So, is there any way to get the _real_ centroid of a LineString? P.S: I have also looked through this discussion. But couldn't figure out the solution.


title: How to change DEM's gradient colour?

text: I'm a newbie about GIS. I'm from Italy and I've got a *.tif file of my region. I've loaded it has raster but I have a question: The dem has a gradient colour magenta. How change this colour to typical "blu- green-yellow-brow" of Greografic representation? Thanks


title: Border style in MapServer

text: I have WMS polygonal layer and want to draw polygon's border with the following style: ![enter image description here](http://i.stack.imgur.com/J2A7x.png) My symbol's description:               SYMBOL        NAME "double_line"        TYPE vector        POINTS          0 0 1 0 -99 -99 0 1 1 1       END      END      CLASS section in mapfile:                   CLASS         STYLE             OUTLINECOLOR 255 0 0             SYMBOL 'double_line'             SIZE 3             WIDTH 1         END      But in this case you can see self intersection in polygon's edge: ![enter image description here](http://i.stack.imgur.com/5Etem.png) What is the righ way to get correct border style?


title: Conversion of GRIB and NetCDF to my database

text: . ![enter image description here](http://i.stack.imgur.com/RfSQo.png) I have downloaded "High Resolution Initial Conditions" climate forecast data for one day, it was in extension .tar.gz so I extracted it in my local directory and I get the files like in the attached image. I think, that the files without extension are GRIB data (because first word in them is "GRIB"). So i want to get data from the big files (GRIB and NetCDF formats containing climate data like temerature & pressure in grid) to my database, but they are binary. Can you recommend me some easy way for getting data from these files? I can't get any information about handling their datasets on their website. Converting these files to .csv would be nice, but I can't find any working program to convert the GRIB files... please help me.


title: IDW Interpolation Raster clipping

text: I used the CLIPPER tool under the 'Raster' menu to clip my interpolation raster with a polygon. the resulting output is not displaying anything even though the process completes without any errors. Also I tried looking other posts in this forum and on google too, but What i got are the posts where some error occurs. i even did the compute histogram, which gives a histogram for my input raster but after clipping the output raster has no histogram. I don't know why this is happening. Probably something is there I missed out.


title: Using TIN - interpolation plugin - to calculate volume

text: I am having trouble calculating the volume of my TIN. I have attempted to create a TIN using the interpolation plugin in QGIS. Here are the steps that I followed.   1. I added a Zvariable row in the attribute table of the vector layer with the value of 500 (I am assuming this corresponds to 500 meters).   2. I used the interpolation plugin in QGIS to create a TIN layer. I created it following steps advised through the qgis website which recommended however that I use a cell size of 5000 for both both X and Y cells but this came up with an error message for me so I left it as it was set by default (please see the screenshot below for details of exactly what I entered).  ![enter image description here](http://i.stack.imgur.com/Rrn7O.png)   3. After, I defined the coordinate system as WGS84, and the attached layer appeared (please see in the below attachment, the layer titled testTIN which is shown in black shading). ![enter image description here](http://i.stack.imgur.com/M8Iej.png) Do you believe that my TIN layer is correct and if so, how do I now go about calculating the volume of the area? I have also attached for you a screenshot of the attribute table of the layer used in the interpolation plugin to create the TIN. ![enter image description here](http://i.stack.imgur.com/6jWKD.png) Thank you in advance for your time. My apologies if my questions are basic but I have not used ArcGIS or QGIS in years. If you prefer writing in German please do so, I am able to read it.


title: OpenLayers Control.WMSGetFeatureInfo does not send request to geoserver

text: I'm refactoring the code used to getfeatureinfo from geoserver running on port 8080 at http://www.edpp.cz/zul_mapa-povodnoveho-planu-mesta/. The old version uses OpenLayers.loadURL and successfully requests the response from the server. The code is like:               map.events.register('click', map, function (e) {                         var params = {                     REQUEST: "GetFeatureInfo",                     EXCEPTIONS: "application/vnd.ogc.se_xml",                     BBOX: map.getExtent().toBBOX(),                     SERVICE: "WMS",                     VERSION: "1.1.1",                     X: e.xy.x,                     Y: e.xy.y,                     INFO_FORMAT: 'text/html',                     QUERY_LAYERS: layers_on,                     FEATURE_COUNT: 1,                     Layers: layers_on,                     Styles: '',                     WIDTH: map.size.w,                     HEIGHT: map.size.h,                     format: format,                     srs: map.layers[0].params.SRS};                      if(offline)                     OpenLayers.loadURL("http://localhost:8080/geoserver/wms", params, this, setHTML, setHTML);                 else                     OpenLayers.loadURL("http://www.edpp.cz/mapinfo/geoserver/wms", params, this, setHTML, setHTML);                  OpenLayers.Event.stop(e);         });      With the new code I try something like this:               popup = new OpenLayers.Control.WMSGetFeatureInfo({                 url: offline ? 'http://localhost:8080/geoserver/wms' : 'http://www.edpp.cz/mapinfo/geoserver/wms',                 layers: unbaseLayers,                 layerUrls: offline ? ['http://localhost:8080/geoserver/wms'] : ['http://www.edpp.cz/mapinfo/geoserver/wms'],                 title: 'Identify features by clicking',                 queryVisible: true,                 maxFeatures: 1,                 eventListeners: {                     getfeatureinfo: function(event) {                         if (console) {                             console.log(event);                         }                         if (event.text.length > 1088) { // nezobrazí prázdné bubliny                             map.addPopup(new OpenLayers.Popup.FramedCloud(                                 "popup",                                  map.getLonLatFromPixel(event.xy),                                 null,                                 event.text,                                 null,                                 true                             ));                         }                     },                     nogetfeatureinfo: function(event) {                         if (console) {                             console.log(event);                         }                     }                 }             });      If I look into Firebug network log, I see that the clicking the map doesn't send any request to the geoserver (it runs fine on localhost though). Is this just a proxy related problem? I've already read about need to set up a proxy when running geoserver on the different port than a website. But it's really puzzling me there is no request being sent. Might it be a typo? I can't see any. What if I am not able to set up a proxy, do I have any other option? PS: All I get as a result is object of nogetfeatureinfo type.


title: Is there anywhere that I can obtain global UV data in either raster or vector format?

text: I'm looking for some UV data (currently with no spatial or temporal resolutions in mind, but ideally split into UV-A and UV-B) for the whole earth. Can anyone recommend a source for data like this or point me in the right direction? Any help much appreciated.


title: Android app: GPS Logger and POI recorder, with offline map

text: I am looking for a good app that can do one or many of my below listed requirements. They are listed with the most fundamental functions first. Preferably the app can do all of it, but it should at least pass the first requirements.   * Have track logger and capability of making own POI   * Use offline maps to display position/track/POI using OSM or something else   * For the twist: I want to use the app for outdoor surveying of habitat/wildlife. Where a POI is recorded I want to be able to select habitat/species from a list of predefined values, so I dont need to type in the value myself. Later I want to export that information into a shapefile and build up a database with observations.  I have not been looking into it too much, but I think ArcMap collector could do this? But I am looking for a free app to easy distribute it to volunteers etc. Any ideas for this? Cheers


title: What software do I need?

text: I run large site focused on history. Our current project expects the use of maps (for simplicity, we use current maps, not historical) with quite many points (custom SVG markers), lines, polygons on it with date (and if possible, time) parameter. I need a tool to create such items on existing map layer, to be able to place this map with features in form of slippy map on my website (without editing, just viewing), and to link clicks on the objects to some function that will load additional data (by some value field taken from the geo database). So far, I was experimenting with QGIS desktop, where I was able to achive - use background map layer (via openmap plugin, google, openstreepmap data) - put on svg markers, lines and polygons, save locally in ?shape file? - add value fields and date to marker table so lets say, I am able to CREATE the features on the map. But, I dont like the look of background map layer, as google is not showing the map 100% correctly in areas where I need it, on the other hand openstreetmap is great but it is too detailed for me (main roads, cities with names, water/rivers, country borders would be sufficient). I'd love to have height shading added to the map, to create look simillar to Google Physical map. While reading a tons of information on the internet, I came to following requirements :   * I need a geo app server with database - doable with geoserver (java) and postgis   * I need to grab data from geofabrik (I plan to use it for the whole world with most of our focus will be Europe, Asia, Africa and North America) and load it into postgis (using tool like osm2pgsql)   * I need to render the tiles from this data with tool like mapnik   * and I need to publish my layer from QGIS to Geoserver, then use some javascript library (leaflet) to show it on my page Am I right ? And how do I pick only those information from the openstreetmap that I need on my map ? I expect it to look similar to this http://www.geofabrik.de/maps/rendering.html Geofabrik TOPO STYLE Footnote : I do not need any special GIS measurement tools, geocoding, path finding, no no. Just map and items to be shown on the map, in various layers (able to turn them on and off in browser), browsable by date, with possible thousands of objects all around the globe. I cannot predict the usage, hardware is no limit for me, so I prefer to create the full stack on our server and run it from there, not to pay for some ready service (that might or might not fully fit into our needs). Thank you


title: What index to use in this situation

text: I am working with 3 tables with 3 x 500 000 rows. There is one column that gives each object a unique ID (UIDN), there also is a column with the geom data (the_geom). There also are some other columns but the 2 above are the most important. When I run my code, it takes a very long time, beceause these 3 tables are so big (I tested the same code with smaller columns and there it runs much faster, so I know for sure that the problem is the size of these columns). PostGIS automatically made a GIST index for the 'the_geom' columns. Beceause it was real slow, I added a btree index for the 'uidn' column. I hoped this would speed things up, but it doesn't. What do I need to do? Do I need another type of index? EDIT: Here is an example of a real query found in my function                           SELECT wbn.uidn, wgo.uidn                 FROM wbn, wgo                 WHERE ST_Intersects (wgo.the_geom, wbn.the_geom) = 't'                   AND ST_Length (ST_Intersection (wgo.the_geom, wbn.the_geom)) > 2                   AND wbn.uidn = f_uidn                   AND wgo.type = f_type      wbn and wgo are two of the big tables I was talking about


title: Why is arcpy.StackProfile_3d -> arcpy.SaveGraph_management creating empty graph?

text: In a python script I'm trying to save a graph to a .png file with graph output from arcpy.StackProfile_3d but the .png file shows an empty graph. The table created from arcpy.StackProfile_3d shows values. Does anybody know why the graph is still empty?               import arcpy     # Set Local Variables     profileLine = "C:/temp/data.gdb/profile"     profileTargets = "C:/temp/data.gdb/raster"     profileTable = "C:/temp/value.dbf"     graphName = "Sample Graph"     outGraph = "C:/temp/graph.png"          # Execute StackProfile     arcpy.StackProfile_3d(profileLine, profileTargets, profileTable, graphName)          # # Execute SaveGraph     arcpy.SaveGraph_management(graphName, outGraph, "MAINTAIN_ASPECT_RATIO")      ![enter image description here](http://i.stack.imgur.com/MvwS3.png)


title: Correcting a table's projection in MapInfo

text: Occasionally I have a MapInfo TAB file that doesn't match up with others. Presumably this is because MapInfo is interpreting it as a different projection to the one it was saved as. I am also assuming that it is possible to have two TAB files of different projections in the same workspace and have them match up. As far as I can see the only way to change a tables projection is to save a copy which doesn't seem to help me. I there a way to switch projections on the fly? Or at least to find out which projection the TAB file thinks it is and which projection MapInfo is interpreting it as? Am I just being thick?


title: How to split a polygon into multiple equal-width polygons in ArcMap 10?

text: I am adding parking spaces into a parking lot, and I need them to be exactly the same shape and not overlap each other or have spaces between. In ArcMap 9.3, there was an editing task that allowed you to split a long rectangle into multiple, equal-width polygons. I am having trouble finding this tool in ArcMap 10. Does anyone know how to do this?


title: Is it feasible to reclassify a raster which has been resampled/smoothed for display purposes?

text: **Background:** I have a Google Earth KMZ file which shows ground movements after an earthquake, arranged into discrete classes: ![enter image description here](http://i.stack.imgur.com/dk9hz.png) From the KMZ file I'm able to download the original PNG images, which are multi-band (RGBA). I used the mosaic tool to create a single, seamless raster over my area of interest. I need to reclassify this image to match the legend above. **The problem:** The problem is that while there are 13 discrete classes in the legend, there is a far greater range of values in the mosaicked raster. I believe this has occurred due to resampling and smoothing of the images (which occurred before the data provider gave them to me) which has resulted in interpolation between the known values. Here is a screenshot of the cells at close magnification: ![enter image description here](http://i.stack.imgur.com/O4e5D.png) For example, the yellow in the legend corresponds to RGB values of 255,255,0 but the yellows in the mosaicked raster vary considerably (eg 239,239,15 or 246,246,8, etc). **My question:** Is it feasible to "resurrect" the mosaicked image, and convert it into something approximating the values in the legend? I realise that this approach isn't ideal, and that it would be far better to start with the original source dataset - but this isn't possible in this case.


title: How to add custom Feature attributes to ESRI Shapefile with Python?

text: I am seeking a way to take an existing ESRI Shapefile that has a Feature set of 200 countries. Each country Feature has an attribute of "NAME." My objective is to create a Python script that adds an arbitrary (for now) additional attribute, say, "POPULATION". Of course I have the OSGeo and GeoDjango modules installed. I'm as far as:                from osgeo import ogr              infile = ogr.Open('sample.shp', 1) #'sample.shp' is a pre-existing ESRI shapefile described above         inlyr = ogr.GetLayerByIndex(0)      Am I missing an OGR function that will allow me to insert Feature attribute fields into an existing Shapefile?


title: Similar color on points and lines

text: I have two feature classes/shapefiles. One containing points and one containing lines. They have both a field called Gate_number. I can use Symbology to give every point with the same gate_number similar color, and i can do the same with the shp containing lines. My problem is that I want to have the same color on the lines and ponits which have similar gate_number. Every point with gate_number=2155 should be green and ever line with gate_number=2155 should be the same color green. Can theis be done in ArcMAp og FME?


title: ArcGIS modelbuilder: Iterate through Excel sheets?

text: I am trying to use modelbuilder (ArcGIS 10.1) to add XY data to ArcGIS, then export the events as Feature Classes. The data comes in one xlsx document with a large number of sheets, each with an X and Y column. So in the end I would have a set of point feature classes. So far I haven't been successful in automating this procedure with modelbuilder. My idea was to use the "Iterate Tables" tool and choose the xlsx document as input workspace; then continue with "Make XY Event Layer" and finally "FC to FC". In order to not overwrite previous results, I have added "%i%" to the output name. To be more precise, the model does run through the first sheet; the second run fails because the model tries to overwrite the first result (despite the %i% appendix). Using "%n%" instead, I get an infinite number of results all using the first sheet. I suspect that either it is not possible to iterate through the sheets of an Excel table, or else maybe something is wrong with my naming system. Any help on this / workarounds much appreciated.


title: How to Enable DriveTime (Minutes) in Network Analyst Service Area Creation?

text: I'm using ArcGIS 10.1 Network Analyst Extension to try and do a drive time analysis for fire station service areas. Using this tutorial I am going through step by step to make this happen, however, I am unable to create the analysis using drive time as an impedance and can only choose distance as an option: ![enter image description here](http://i.stack.imgur.com/SPba6.png) My road network dataset has a field that shows the time it would take to traverse the segment at 40mph. How do I enable DriveTime as an option in my analysis settings?


title: Is it possible in QGIS print composer to include a numeric scale in an html box?

text: For work I always have to create a table with my name, project name, date etc. into my map layouts. I also need to include the scale like "1:150000" in that table. So far I am solving that by including an html table with the html box function of the print composer. Is there a way to automatically include the right scale? So far I only managed to automatically include the date by using a little javascript. Maybe something similar is also possible for the scale?


title: How to clean/edit GPX (tracks) files recorded from GPS device?

text: During a long trip, I have recorded tracks of my journey with the GPS of my smartphone. I have ended up with about 30 GPX datafiles. Now that I am back, I am trying to draw a map of this trip. I imported each file in QGIS as vector layer, and visualized the tracks with _OpenLayers Plugin_. Now I need to edit some of these files:   * basically trimming the end and beginning of many files   * merging some tracks   * adding a title or a description for each track What is the way to do it with QGIS ? I would like to keep as much data as possible from the source data (timestamp, elevation)


title: ArcGIS Online map labels do not appear in ArcGIS app for Android?

text: I have a map on ArcGIS online. The map shows labels just fine as you can see on my web map (http://bit.ly/1q7hc0h). The problem is that when I view the map on may Android device using the ArcGIS app the labels do not appear. I have the most recent version of the Android app and I do not have scale dependent labels set in ArcGIS Online. Any ideas on why the labels are not showing in the mobile version?


title: Subquery in SearchCursor with ArcPy?

text: I'm searching for the value from the field "demA_demA" which has the lowest elevation by having the value 8 from the field "outTable_2". ![enter image description here](http://i.stack.imgur.com/AdMmz.jpg) How can I do this in the following script.               import arcpy     import os     from arcpy import env     from arcpy.sa import *     arcpy.CheckOutExtension("Spatial")     env.overwriteOutput = True          env.workspace = r"D:\Julia_T\projekt\aufgabe2\testdatenC"          inASCII = r"D:\Julia_T\projekt\aufgabe2\testdatenC\2007040_cloud_free.asc"     outRaster = inASCII + "_out.img"     rasterType = "INTEGER"     # Process: ASCII in Raster     arcpy.ASCIIToRaster_conversion(inASCII, outRaster, rasterType)     # Process: Projektion definieren        arcpy.DefineProjection_management(outRaster,"GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHER OID['WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]]")     # Process: Sample     outTable = r"D:\Julia_T\projekt\aufgabe2\ergebnisse3\outTable.dbf"     Sample(outRaster, outRaster, outTable, "NEAREST")     #works     # Join the table     # Table view     arcpy.MakeTableView_management(r"D:\Julia_T\projekt\aufgabe2\ergebnisse\outTable.dbf",     "out_tview") #works          joinField = "Y"     joinTable = r"D:\Julia_T\projekt\aufgabe2\ergebnisse3\demA.dbf"     arcpy.AddJoin_management("out_tview", joinField, joinTable, joinField, "KEEP_ALL")     # Persist the layer on disk     arcpy.CopyRows_management("out_tview",     r"D:\Julia_T\projekt\aufgabe2\ergebnisse3\joinTestC.dbf")     # Execute CreateFileGDB     arcpy.CreateFileGDB_management(r"D:\Julia_T\projekt\aufgabe2\result", "test1.gdb")     #works     outLocation = r"D:\Julia_T\projekt\aufgabe2\result\test.gdb"     # Execute TableToTable     arcpy.TableToTable_conversion(r"D:\Julia_T\projekt\aufgabe2\ergebnisse3\joinTestC.dbf",     outLocation, "jointest1") #works     # search for the lowest elevation with the value 8     with arcpy.da.SearchCursor (r"D:\Julia_T\projekt\aufgabe2\result\test1.gdb\jointest1",     ["outTable_2", "demA_demA"], """""demA_demA" = min("demA_demA")AND "outTable_2" = 8""") as     cursor:         for row in cursor:             print row [0], row [1]      Do I have to add the geodatabase because of an subquery? How should I write the subquery? I'm pretty new and working with ArcMap 10.2.1


title: no route detected in network

text: why this warning box appear when i click at solve button in network analyst toolbar ? - No route from location "Graphic Pick 46" to location "Graphic Pick 61". Error: No solution found.


title: Get CartoDB query by distance to work

text:
Has anyone tried this tutorial: Query by distance ? This is a copy of the final code which they use for the map:                           <link rel="stylesheet" href="http://libs.cartocdn.com/cartodb.js/v3/themes/css/cartodb.css" />

         500          100          50          10
                I can see all of the points, but nothing else works, when I make the following changes:   * I replace var layerUrl with my viz.json link.   * I replace the name of their table, with the name of my table in the sql inside var subLayerOptions   * I replace the name of their table, with the name of my table in the css inside var subLayerOptions   * I replace the name of their table, with the name of my table in the cartocss inside sublayers[0].set Am I missing something? Does anyone have an example of this that they have got to work, with their data? Could I see the code? Thank you!!


title: How to measure sinuosity in ArcGIS Desktop?

text: I am trying to quantify the sinuosity of a GPS track (shows a loop pattern) and cannot seem to find any way to do this. I need to be able to quantify sinuosity for the removal of some unwanted data. Is there a tool to measure sinuosity in arcmap, i am very unfamiliar with the codes in ArcMap 10?


title: Mapping between two different agencies roadway shapefiles of Las Vegas in ArcGIS

text: I have two roadway shapefiles of Vegas. But they do not overlap. One of the file has a certain shift from its exact location. Can anyone help me how can I move the location so that they overlap exactly and I can map those two different shapefiles in ArcMap?


I'm new when it comes to GIS and PostgreSQL. I have searched everywhere for the cause of this problem but I can't find the solution. I was following the documentation here First i created a table:               CREATE TABLE tableA (id integer, polygon geometry);      I couldn't use `st_polygon` because that was causing a "`type not exist`" error though GIS is fully updated. (through Ubuntu repos) Then I inserted a row. I couldn't do it like how it was described in the documentation because that was giving me a weird errors:               INSERT INTO tableA VALUES (1,         st_polygon ('polygon ((40 40, 40 60, 60 60, 60 40, 40 40))', 0)         );              ERROR: lwline_deserialize: attempt to deserialize a line which is really a Invalid type         SQL state: XX000         Context: SQL function "st_polygon" statement 1               INSERT INTO tableA VALUES (1,         st_polygon ('polygon ((40 40, 40 60, 60 60, 60 40, 40 40))', 0)         );                ERROR: function polygon(unknown, integer) does not exist     SQL state: 42883     Hint: No function matches the given name and argument types. You might need to add explicit type casts.     Character: 44      In the end I was able to insert it like this:               INSERT INTO tableA VALUES (     1,     ST_PolyFromText ('polygon ((-30.8947 138.6522,     -34.8951 138.6522,     -34.8952 138.6512,     -34.8947 138.6512,     -30.8947 138.6522))', 1)     );      Now when I perform               SELECT * FROM tableA WHERE ST_OVERLAPS(ST_BUFFER(GeomFromText('POLYGON((-34.8947 138.6522,     -34.8951 138.6522,     -34.8952 138.6512,     -34.8947 138.6512,     -34.8947 138.6522 ))'),0), polygon);      It outputs the error:               ERROR: Operation on two geometries with different SRIDs     SQL state: XX000      Anybody got a clue what it could be now and how to fix it?