title: Create easy Annotations from point features with python in ArcGIS 10.0

text: Because in ArcGIS 10.0 there is still no simple toolbox function to convert labels to Annotation, I try to go a substitude way with the callable Toolbox function **[TiledLabelsToAnnotation][1]** on a polygon shape, which covers all of my data area (= 1 polygon). Thank you for your help - the script is now running, because of the simplifications for the posting!               import arcpy          # Settings     scale = "30000"       folder = "D:\Temp\LBL.gdb"       arcpy.env.workspace = folder            # Mapdocument     map = r"D:\Temp\MAP-LBL.mxd"            # Local variables:     dataFrame = "Layer"       grid = "Grid"       group = "Group"       suffix = "Anno"       fid = "tile"       link = "STANDARD"  # not linked = "STANDARD" / linked = "FEATURE_LINKED"        unplace = "GENERATE_UNPLACED_ANNOTATION"            # Process: Tiled Labels To Annotation     arcpy.TiledLabelsToAnnotation_cartography(map, dataFrame, grid, folder, group, suffix, scale, "", fid, "", "", link, unplace)


title: Is there a way to know if an error occurred while loading a WMS layer, in OpenLayers 2?

text: First of all, I already asked this question in stackoverflow. I don't want to spam by asking again here. But I think the question fits better in gis.stackexchange. I am sorry for the replication. Feel free to delete one of the two questions. Let's start! I am working with OpenLayers 2. I am using the events "loadstart" and "loadend" from here. When I use OpenLayers.Layers.WMS, i want to have some kind of callback when there is a loading error. For example, when I enter an incorrect SLD filter, I want to have some feedback about the layer not being loaded correctly. Unfortunately, "loadend" doesn't provide this kind of error detection. I already found this question where someone else has the same problem. Is there even an event that signals this kind of error? Is there something like a "layerfail" or "layererror" event? How could I implement this? Thank you.


title: Break line into 100m segments

text: I have a line geometry in PostGIS and I am wondering how I would break the line into 100m segments. Is there a function for this?


title: Orthogonal Height Profile along line (sextante -> v.segment, r.profile?)

text: While trying to find a way to find "valleys/ridges" within a DEM, I'd like to know if there are tools, plugins or "methods" known or available to automatically generate height profiles along a given path. Or maybe a tool that (allows to) calculates the height profile (UNIX preferred, no Windows at all) at a certain position, takin coords and azimut as input? **\--- / edit 19.12. --- /** Even though I'm yet not used to sextante and GRASS at all, I've searched and read a lot. Shouldn't it be possible to solve above discribed problem by using the sextante modeller? If there is a single line within a line-shape plus a raster elevation map... Shouldnt it be possible to...   1. Step: generate points using _v.segment_ with direction (within line) as attributes for the resulting points **P**.     2. Step: generate 2 new points for every **P** within a fixed given distance, which mark the start and end of an orthogonal line (angle from "directon" at step-1), where **P** is the middle resp. point of intersection.    3. Step: calculate profile of this orthogonal using _r.profile_   4. Step: write/append elevation profile to a CSV with ~~lat/lon/height~~ **x/y/elev** , so that it can be further processed or simply beeing 3D illustrated via Speadsheet, gnuplot,.. The lan/lot is delete cause I didnt point out I would like to have it "straightned". Otherwise I could simply wrap the area from the DEM. An please - as mentioned - I'm new to GRASS and even less when it comes to the modeller :-/ It's not too long ago I started using the sextante toolbox to improove the possibilities within QGIS. Below show my two rasters examples. First one is straight ASTER elevation data overlayed by a rugh drawn line from a geological map, respresenting a tectonic fault. Second on is the terrain convergence index (tci) of the DEM, where I get the impression that I can see this line. ![ASTER elevation data](http://i.stack.imgur.com/f0BFe.jpg) ![terrain convergence index](http://i.stack.imgur.com/5XPYI.jpg)


title: tif file creation

text: Does anyone have any idea how to get around this problem? I've been looking for a solution to no avail. I'm trying to convert a shp file to a grid, and the tif is hanging up. Output raster layer: H:/reclassify/test2.tif The above files could not be opened, which probably indicates that they were not correctly produced by the executed algorithm Checking the log information might help you see why those layers were not created as expected This algorithm requires GRASS to be run. A test to check if GRASS is correctly installed and configured in your system has been performed, with the following result: GRASS seems to be correctly installed and configured


title: Default workspace in QGIS 2.0

text: I have been using QGIS for quite some time. It always happens that I change the position of toolbars and panels to make the workflow more efficient. But then when I close QGIS and I have to start a new project, the position of toolbars and panels remain the same. I have to put them back individually to their position. My question is, is it possible to create a default workspace? I mean, if I want to put back at once all the tools and panels in their position, I just have to click one or two buttons. Like in photoshop, they have several workspaces. Thank you for your help! Jake


title: ArcGIS JavaScript API Auto Panning to Cursor Position

text: Does anyone know if it's possible using the ArcGIS JS API to make it so the map auto-pans based on the user's cursor position if they get within a close distance to the edge of the map viewing area. This would allow the user to continue editing and not have to do some other action (such arrow keys, or a pan icon) to get the map to pan to an area they want to continue to edit in. I've been looking for some examples of this and I can't seem to find any. Any help or direction would be appreciated.


title: How to best intersect vector data with raster data?

text: I have Open Street Map [OSM] data imported so that I can access the geometry of buildings. I have also loaded SRTM digital height data from NASA from a GeoTIFF raster. I have projected the raster into EPSG:3857 (the same as OSM) from the original EPSG:4326. For the SRTM tiles that I am interested in, I have warped them into 1 meter resolution with cubic spline interpolation via ST_Resample and produced smaller tiles from that with ST_Tile(as a single INSERT .. SELECT which is very handy). **Now I run into a problem:** The ST_Intersection() of the resampled tile and the building polygon simply produces a square bounding box, which is not what I want. I want to get at the exact x,y points on the raster for the intersection because I'm using this data for generating 3d voxel scenes (1 pixel = one voxel) and I want to produce a 2d map of the rendered area to pixel-point precision as well. **What I am currently doing:** I use ST_Clip() to clip the raster by the geometry of the buildings, which puts NODATA data into the returned clipped raster everywhere except where the building is actually situated. This is like using an template to "burn" the image of the building into the raster, leaving the points inside the building completely intact while everything else is burned away. I dump the resulting clipped raster with ST_PixelAsPoints() or ST_PixelAsCentroids() which gives me x,y positions on the raster where the building is located. **Is this the right approach?** I have heard of something called a 'pointcloud' but PostGIS doesn't seem to have this available. Converting the raster into polygons doesn't work, because the x,y data from the raster is lost.


title: setting up php_mapscript in portable GIS

text: I'm new here, and I'm also new in webgis, I have download portable gis, and I love this new version with lot more tools, I have a question regarding to mapserver that comes with portable GIS v2, how to set up or add php_mapscript extension to php in xampplite, I have try to do like this but failed :   1. Add the necessary files (libmap.dll, gdal17.dll, etc.) to a new folder (foldername = A ). I get those files from new version of MS4W.   2. Add the folder path to the windows path environment.   3. Modified php.ini to add php_mapscript When I try to launch apache server it saya something like "couldn't load msvcr90.dll", I copied the file from windows folder and add it to folder A, new error says "wrong msvcr version". Can anyone help me, so I can use php_mapscript extension. Regards, Zery


title: ArcGIS Server 10 - Licensing by Core

text: I know awhile back they changed the licensing model to license by the core. So if you have one server license, and you want to install it on a four core server?   * Do you need four licenses to be legit?   * If you want to restrict Server from using the other 3 cores, does ESRI provide information on how to do this?


title: ArcGIS server on Amazon EC2 vs buying your own machine. Pro's and Con's

text: What are the pros and cons using this virtual machine compared to buying your own?


title: arcgis - updating joined excel file when adding a new sheet

text: When adding a new sheet to a joined Excel file in ArcGIS, the new sheet will not show till I exit both arcgis and excel. Is there a better way to update the connection...? edit: sorry for the delay and thanks: I have some projects based on an excel workbook joined to a specific feature class (based on an identical field). That workbook has many worksheets, each one joined to a different field in the feature class, and sometimes there is the need to add new data, in a new worksheet (I wouldn't want to create a new workbook for each new table that needs to be joined), and create a join. That's where my problem is: I can't find the way to refresh arcmap, so it recognizes a newly added worksheet in the excel, without restarting arcmap. Any advice? Thank you, Chaim


title: Are there any convenient sources of color palettes to use in ArcGIS software while producing maps?

text: Since it was asked in Beta and I want to know the answer, I post it here.


title: Eliminating/minimizing distortion in 200 NM EEZ buffers

text: I frequently need to calculate buffers of 200 nautical miles from country or baseline polygons in order to map maritime Exclusive Economic Zones (EEZs). However, I've noticed that when generating the buffers in QGIS, I get different results depending on the projection specs of the layer I'm buffering from. For example, in a recent project I'm using an azimuthal equidistant projection, and the buffer turns out differently depending on where the projection is centered. This also means that even if I center the projection on the country, there is often still significant distortion at more distant corners of the same country. Is there any way to generate buffers without any distortion (i.e. equidistant from each node on the polygon from which they're being calculated)? Is there any tool (or different projection) that would generate this result? How do others usually deal with this issue? Thank you very much in advance!


title: How to convert zip codes to latitude and longitude?

text: I have a data set that includes a list of zip codes. Is there a way in R to convert these zip codes to latitude and longitude coordinates?


title: Calculating the minimum distance between points and polygons in QGIS

text: I have two vector layers: one point and one polygon and I want to calculate the minimum distance of each point from the polygons (ie the distance from the closest point of the closest polygon). Is there any available plugin in QGIS for doing this? Thanks!


title: Sum a field in buffer intersection

text: I have looked and do not see this answered in any other discussions. Purpose: Site suitability analysis for brewers spent grain processing facility Method: Locate intersecting buffers with greatest amount of spent grain - need field to sum barrels based on intersections (ex. Mountain Brewing Co. produces 3,000 barrels a year and Sierra Nevada produces 800,000 barrels a year - if the radii from the two breweries intersects the intersection would have a sum of 803,000 a year) Using ArcMap 10.1, I have one shapefile that I have a 15 mile radius around each point (brewery). I am trying to identify all intersecting buffers by summing a specific field. I would eventually like for the intersections with the greatest number to be a 3D representation of barrels per year based on the sum field (I don't think I can do this in ArcMap 10.1 but need to transfer everything over to ArcScene, but I'm hoping not).


title: How can I style OSM data in geoserver?

text: I've been trying to style OSM data imported using osm2pgsql and served through geoserver. The tables generated by osm2pgsql are: planet_osm_line, planet_osm_point, planet_osm_polygon, apart from the others generated through osm2pgsql slim mode. Default line/point styles work well but I'm unsure how to achieve the effect like OSM's example map. I have tried using other style files I managed to find, but they are outdated and don't validate with geoserver. Others like "OSM in a box's" source or http://svn.osgeo.org/osgeo/foss4g/benchmarking/wms/2011/geoserver/styles/ have separate sld's for roads, railways etc more than the layers created by osm2pgsql. How can I style geoserver's output to look similar to OSM ? Any help would be greatly appreciated. Thank You! I'm open to using SLD's which are available, but older ones don't validate correctly.


title: excel column does not appear when doing join after doing modification to the column

text: Using ArcGIS 10.1 I am trying to add " **/1** " my column provided Shp File how it appears in the county supplied excel spreadsheet under "parcel #" 43002/1 18524/1 1666/1 4270/1 Currently, the county provided shapefile only displays 43002 18524 1666 4270 How do I ADD the /1 to the already current number within my shapefile to make it look like 43002/1 etc? There are over 13,000 tracts so doing it by hand it out of the question. When I go into excel and delete / REPLACE, when I do the join, that parcel column does not show up, where before it did before I did the delete of the /1. I went into format, number and used zero decimal placement. Still nothing. 1 - If someone can tell me how to add " /1 " to the parcel numbers 2 - Is there something I can do with excel so when I do the join, the parcel num column shows up and not disappear


title: Issue with using TopoJSON server API to return TopoJSON object with properties

text: I am creating a node application and using TopoJSON server API reference (see the link below) to convert a GeoJSON object to a TopoJSON object. TopoJSON server API reference I used these three lines of code below which returns the TopoJSON topology but without returning the properties from the GeoJSON object.               var collection = {type: "FeatureCollection", features: […]}; // GeoJSON     var topology = topojson.topology({collection: collection}); // convert to TopoJSON     console.log(topology.objects.collection); // inspect TopoJSON      However, I also want to include the properties of each topology to the final TopoJSON object, and according to the API, I need to use this function below copy all properties from the GeoJSON object:               function propertyTransform(properties, key, value) {       properties[key] = value;       return true;     }      But my issue is I don't know how to use this propertyTransform function as an option inside of the topojson.topology() function in order to return the TopoJSON object with properties. Does anyone know how to do it or provide me some example code? Thanks a lot! Vanni


title: Splitting polyline by points using ArcGIS 10.2 for Desktop?

text: I have a layer containing a curved polyline and I would like to split it at multiple point locations, which I have in another layer. In ArcGIS 10.2, I can do this with the tool "Split Line at Point". It is working, but now I have a chaotic attribute table. The order of each feature in the table is not the same as it is geographically. How can I order the polylines geographically?   1. The line is curved at many sections, so order only by latitude or longitude is not possible. (I think that this method is applied by ArcGIS, because the order of the points for splitting is the correct.)   2. I have tried it with searching for the same End_X and Start_X, but it seems that they have not always the same coordinates. ![Splitted Polyline with Object ID of sections](http://i.stack.imgur.com/N3yga.jpg) ![Splitted Polyline with Object ID, splitting points with ID](http://i.stack.imgur.com/a2yzU.jpg)


title: Unable to install ECW support on lubuntu 14.04

text: I try unsuccessfully to install the ECW support (mainly to use ECW file in QGIS 2.2). We attempt to follow and adapt the howto provided by makina-corpus. My first supprise was the repository ubuntugis-unstable is unavailable for trusty and the repository http://qgis.org/debian doesn't have libgdal-ecw-src package. Finally in desperation we add ubuntugis-unstable for the raring distribution in source.list. And apply the following commands               chmod +x ERDAS-ECW_JPEG_2000_SDK-5.1.1.bin     ./ERDAS-ECW_JPEG_2000_SDK-5.1.1.bin # Select option desktop read only     sudo cp -r ERDAS-ECW_JPEG_2000_SDK-5.1.1 /usr/local/     sudo ln -s /usr/local/ERDAS-ECW_JPEG_2000_SDK-5.1.1/Desktop_Read-Only/lib/x86/release/libNCSEcw.so /usr/local/lib/libNCSEcw.so     sudo ldconfig     sudo apt-get install libgdal-ecw-src     sudo gdal-ecw-build /usr/local/ERDAS-ECW_JPEG_2000_SDK-5.1.1/Desktop_Read-Only      The ECW format seems to be installed               gdalinfo --formats | grep -i ecw       ECW (rw+): ERDAS Compressed Wavelets (SDK 5.1)       JP2ECW (rw+v): ERDAS JPEG2000 (SDK 5.1)      But .... Gis crash when a ECW raster file is open with the following message :               QGIS died on signal 4Abandon (core dumped)      If anyone have an idea


title: TIF format to view a vegetation Map with QGIS

text: I'm a new user of QGIS and I'm trying to view a map present in an online GIS dataset. There are several formats for this dataset and one of them is a TIF folder with several files: .tif, .tif.aux, .tif.rrd, tif.vat.dbf and I have no idea how to handle these files to view the map with all its informations about the vegetations of a region. * * * Thank you for your answer. Effectively it was easy to do. In my research I'm tring to view a vegetation map and to have also the possibility to separate each layer. So I went in the website of vegetation mapping project to download what I would need: http://www.vegmad.org/default.html In this way I have tried to find the right file to download, and every time (always using QGIS) I can see only the image of Madagascar but I can't distinguish each layer. Do you know how to do? * * * I have tried to do what you said, but it's not so easy as I expected. I'm showing you what I have obtained downloading MODIS Dataset Geotiff, and which classes I have opening Classes.xls file. And I think that before I separate each layer in different raster files I need to give the exact classification to each pixel value range, but I can't treat those data... ![enter image description here](http://i.stack.imgur.com/335Q5.png) In QGIS first of all I have opened Modis_classification2.tif and what I observe is more pixel values than what that I have in this classification file.xls. This is shown below... ![enter image description here](http://i.stack.imgur.com/4HGEN.png) So what you think that i should do? Thank you again!


title: Calculating distances between consecutive points along rivers (polygon) using ArcGIS Desktop?

text: I have radio-tracking locations of semi-aquatic mammals. Now, I would like to calculate the distance between the consecutive locations using the river as the animals would swim (opposite to as the crow flies). I have two different files: a) Polypoint layer with all the tracking locations b) Polygon layer with the rivers, streams and ponds (and some buffer around it) These approaches seem not to work for me: Network analyst (which seems to be most likely the easiest way): no license to create layer The linear referencing tool to locate features on a route: My animals move once in a while a bit off the beaten water track, thus I added a buffer to the water bodies = polygon layer I'm using ArcGIS Desktop 10.0.


title: Exporting QGIS shp file to dxf-Which dxf format?

text: As I understand it, dxf comes in a range of format types. QGIS (v2.0 'Dufour') can export any layer to a dxf file. However I have a hydraulic modelling package which relies upon importing the dxf file as the R12 format. This import does not seem to work. I was wondering which dxf format does QGIS export to? If it is not R12, this may be the reason why this is not working.


title: Converting KML To Layer using ArcPy losing objects?

text: I am trying to a convert kml file to a layer using python. This is my code:               import arcpy      arcpy.KMLToLayer_conversion('C:\kml\\result.kml', 'C:\kml', 'result')      But the layer has lost a lot of objects as compared to the layer converted by _Conversion tools_ -> _kmltolayer_ I attach 2 screenshots: ![layer, which was converted by Convertion tools -> kmltolayer](http://i.stack.imgur.com/dcPwu.png) ![layer, which was converted by python script](http://i.stack.imgur.com/df3HU.png) What could be the problem?


title: Edit OSM Map Graph

text: Can I use pgrouting (http://pgrouting.org/index.html) to change OSM data? My goal is to change the cost of some ways (increase and decrease) in order to change the routing calculation process. If possible, after edit, is it possible to save as a .osm file to use on NavIt? **Resume:** Instead of create Tags inside of OSM data, I want to alter his graph structure. It's like a filter of the OSM original data. Thanks for reading that ;-)


title: How to create an offset polyline in arcpy?

text: I am trying to create an offset polyline using arcpy, identical to the function in ArcObjects (constructoffset). After trying to access ArcObjects in arcpy, which I found difficult for those not well-versed in COM languages, and after much searching for functionality in the regular geoprocessing toolboxes, I discovered that this functionality seems to only exist in ArcObjects. There may be ways to program this functionality in arcpy using trig, but that seems really hard. Essentially, the product I'd like is an offset polyline that preserves all original geometry and curvature. Simply using a SearchCursor on the input and building a new feature class based on its geometry requires some algorithm I don't know. This is because without an algorithm based on slope, etc., a one-size-fits-all addition or subtraction will not preserve an offset shape for all slopes. Any help?


title: Simplifying Shapefile without risking Spatial Accuracy?

text: I am building a feature for a web application. Basically would use the ABS (Australian Business Statistics) data to render the regions on the Google Maps. I am able to implement the basic working prototype however the ABS dataset (65 MB Shape File) is huge and contemplated to use the tools   1. Map Shaper   2. QGIS  for simplification. The rendering time is poor without simplification as the data set is huge. Map Shaper was not really helpful because the shape file (65 MB) didn't get imported into the database after simplification because of error 333 object not found (something like this). QGIS helped me to simplify the shape file from 65 MB to 4 MB, 2 MB, 1 MB depending upon the tolerence level (0.0001, 0.0020,0.0080 etc). I was able to successfully import this into the database and used it to render on maps. The rendering time has been significantly improved, however the spatial accuracy is going for toss because of the level of simplification. The functionality here is to store regions into the database for a user and then notify the users when a delivery matches. We are using Grails, Postgres (using PostGIS), GeoJSON here to implement it. Could anyone here assist me in simplifying the shape file (Statistical Area Level 2) of size 65 MB to atleast 1-2 MB without risking the spatial accuracy when we zoom in to the lowest level?


title: Is there a good library for WMS "Tile Servers"?

text: I am looking for WMS Tile Servers that I can consume with "OpenLayers". **One:** Are there any good library / resources online that collate WMS tile servers or where is a good place to start looking? (Street maps & others etc) **Two:** Are there any **free** tile servers available?


title: Define a circular neighborhood for a custom raster calculator function?

text: I have a 1 meter resolution raster containing values from 0-3 which represent increasing levels of habitat suitability. I need to identify circular regions that have contiguous pixel values. For example, I define a 50 meter radius, and I am looking for regions with values of 3. I need to know if every cell within a 50 meter radius of each target cell has a value of three. If true, then I need to return a true value (1) to the output target cell, if false, then return a false value (0). This is similar to a focal statistics function using a circular neighborhood. Region group will identify regions of contiguous cells, but I need the regions of contigous cells to be circular and of a minimum diameter. I'm using ArcMap 10.2 and any help is appreciated. Thanks


title: Is there a way to shift data in an attribute table so that all attributes are moved to a different coordinate?

text: I have several shapefiles where the attribute table is composed of timestamped data that is appended to the coordinates. Because of the nature of my data, I need to apply a time lag so that the data from a given time are placed a few minutes before. In other words, in Excel I would just delete the desired number of cells and have everything shifted up to the new coordinates. The reason I do not do it in Excel is that I need to perform 'trial and error' until I find the right time lag to apply and it would be a lot faster if I could do it straight in QGIS where I can visually explore my data.


title: Setting Point Cloud Data Colour from Raster

text: I was reading this article: http://blog.safe.com/2013/10/14-ways-to-take-charge-of-lidar-data/ Number 12 states that you can change the colour of your point-cloud data by using a raster image. Does anyone know the process for doing this within ArcGIS. I am looking for a way to make the visualization of lidar data really pop and I have hi-res aerial photography that I could use.


title: How can I get proj4text from srtext?

text: I have a custom projection used in some shapes on a prj file. I want to insert this projection in the table spatial_ref_sys PostGIS. The data in the prj file will SRTEXT table field, how can I get PROJ4TEXT from SRTEXT??               INSERT into spatial_ref_sys      (srid, auth_name, auth_srid, proj4text, srtext)     values      ( 101010, 'epsg', 101010, '?????????????',      'GEOGCS["Carta Geral",DATUM["Not specified (based on International 1924 ellipsoid) (EPSG ID",SPHEROID["International 1924 (EPSG ID 7022)",6378388.0,297.0]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",200000.0],PARAMETER["False_Northing",5000000.0],PARAMETER["Central_Meridian",-51.0],PARAMETER["Scale_Factor",1.0],PARAMETER["Latitude_Of_Origin",0.0],UNIT["Meter",1.0]]');


title: Using Flex what event is fired when loading the popup defined in my layer tag?

text:
Here is my Config.xml ...snip...                   <layer label="Pop-up (turn on/off)" type="dynamic" visible="true"               alpha="0.8"               url="xxx.xxx.xxx">               <sublayer id="0" popupconfig="popups/Parcels.xml"/>  -->              </layer>      </operationallayers>      So, when I click a point on the map a popup appears. That popup displays information about a parcel that was clicked, provides a "zoom to" button. What event fired to make that occur? I'd like to heavily customize what occurs when I click the map. Bonus points if you can verify which object is being used to render the popup ## Parcels.xml               <?xml version="1.0" ?>     <configuration>         <title></title>                 <description>             <![CDATA[Account# {Strap}
Owner: {Name1}
Address: {Locn} {Locd} {Locs}

**2014 Values**
Just(Market)Value: ${Just}
Land Value: ${Land_Val}
Improvement Value: ${Improvemt}
Assessed Value: ${Assessed}

Homestead(x=yes): {Homestead}
Land Sq Ft: {Lsqft}
Total Building Area: {GrossArea}
Total Living Area: {Living}
Year Built: {Yrbl}

Additional Information
]]>         </description>                        <fields>             <field name="Strap"/>  Used in title, do not repeat in content                  <field name="Yrbl" alias="Year Built">                 <format usethousandsseparator="false"/>             </field>                  <field name="Land_Val" percision = "2">                  <format usethousandsseparator="true"/>             </field>             <field name="Improvemt" percision = "2">                  <format usethousandsseparator="true"/>             </field>             <field name="Just" percision = "2">                  <format usethousandsseparator="true"/>             </field>             <field name="Assessed" percision = "2">                  <format usethousandsseparator="true"/>             </field>             <field name="Lsqft" percision = "2">                  <format usethousandsseparator="true"/>             </field>             <field name="GrossArea" percision = "2">                  <format usethousandsseparator="true"/>             </field>             <field name="Living" percision = "2">                  <format usethousandsseparator="true"/>             </field>                  <field name="SALE_DATE">                 <format dateformat="shortDate"/>             </field>            </fields>                   </configuration>


title: How to Dynamically refresh / reload a KML layer in OpenLayers

text: I have previously asked this question to no avail so I thought putting this up here would be a good idea. **Here is an example on:** How to Refresh / Reload a KML layer in OpenLayers. Dynamic KML Layer. See my answer below. * * * **TLDR:** See my answer below.


title: Callout box labels converted to annotation justifying strangely

text: I am helping a co-worker out and she had some callout box labels that she converted to annotation. Most of them look normal - they were originally right justified and maintained that justification - but when she edited the label information and added an additional line, the additional bottom line looks full justified. Kind of. The word spacing is huge but the character spacing appears to be the same. The other lines still look right justified just as far as the box used to extend (meaning they are both in line on the right side but there is a ton of space to the right of them as well). It's a little hard to describe so I'll include a screenshot. We tried even moving all of the lines into one line and nothing changed in the appearance of the box/layout. All of her settings look correct, however. I inserted a new callout box as text and mimicked everything there but the box justification didn't change at all. Any ideas what might be happening and how to fix it without having to insert new labels? ![Top Label is the Weird Updated One](http://i.stack.imgur.com/YIPYD.jpg)


title: Collaborative map editing using local openstreetmap server

text: I tried to make a collaborative web-map solution. that is allowing multiple users to edit map features and enter data collaboratively. I read that OpenStreetMap give this solution. I've installed a local osm server using switch2osm tutorial. Now I can see a map canvas on a web page. **But** , How can I manage map data. I mean:   * Defining features categories.   * Defining users as editors.   * Assigning privileges to users according to geographic area.  Is there any other open source solution, more easy to set up and meet my requirement?


title: Labeling expression to label one value in a attribute table column

text: I am using QGIS 2.0 on Windows. I have a shapefile for a point layer that has the the following attribute table: ![Attribute Table](http://i.stack.imgur.com/iRjcs.png) I am trying to label using the Status column and only want the value that equals Office Location to show on the map. I tried the expression Status = 'Office Location' and I get either a 1 or a 0 show on the map. I cannot figure out the appropriate expression to do what I want and I haven't been able to find anything specific to this on here or through a Google search. I could just make another column an only provide data for the specific values that I need, but that seem counterintuitive. **How do I label a single value and not show the other values?**


title: Problem with XElement usage in ArcGIS Runtime SDK for WPF Sample C# code

text: I was attempting to use some of the code from an ArcGIS Runtime SDK for WPF Sample that adds a shapefile or raster data to a map. I've not been able to get past this one diagnostic error so that I can test the code . . .               // Read the JSON response as XML     XmlReader reader = System.Runtime.Serialization.Json.JsonReaderWriterFactory.CreateJsonReader(downloadDataEventArgs.Result, new XmlDictionaryReaderQuotas());     // Get the root XML element      XElement root = XElement.Load(reader);     // Query for the "geometryType" element     XElement geometryType = root.XPathSelectElement("//geometryType");      When it gets to the last line, I get the error message "'System.Xml.Linq.XElement' does not contain a definition for 'XPathSelectElement' and no extension method 'XPathSelectElement' accepting a first argument of type 'System.Xml.Linq.XElement' could be found (are you missing a using directive or an assembly reference?)" How do I get the code to recognize the Extension Method XPathSelectElement? By the way, I have included the System.Xml and System.Xml.Linq libraries. Thanks, Renee


title: .rul topology rule

text: How can I create * _.rul_ *file for topology in ArcGIS?


title: Use csv file as raster file and stitch them together

text: I'm having an issue trying to stitch together 100 32bit Tiff files (coming from thermal imagery video and extracted as 32bit Tiff). I used Hugin to mosaic the pictures extracted from the video which give good result but has the bad habits to distort the finale stripes. I can as well extract as CSV, the file is then a matrix of 640 columns and 480 lines of temperature value. I was wondering if there is anyway to stitch those csv together in a way to obtain a no distorted finale file using R by exemple. Is a script could be able to say ok I recognize this combination of value in those two csv and I can stitch them together. We considered some structure from motion software to do this work but they usually don't handle the 32bit Tiff. It is a bit unusual question but the development of UAV technology will lead to the multiplication of data as movie that will need to be transfer in GIS application so it might be interesting to find a solution.


title: Inverse distance weighted interpolation for categorical and/or binary datasets

text: I want to perform an `Inverse distance weighted (IDW) interpolation` to a series of points. However my data set can be `categorical` or `binary` and just like the `IDW` I want to interpolate the new values (categorical or binary) based on the highest frequency of closest points values surrounding the prediction location: > A if the majority of the closest points are A > > C if the majority of the closest points are C > > ... In addition I want to avoid using close points that are behind a barrier, that can be a wall or just an elevation. This is how the problem/dataset looks like: ![enter image description here](http://i.stack.imgur.com/dd0js.png) The grey points is the point I want to interpolate and the points behind the barriers can't be used in the process. The output of the interpolation process would be A or B, in this case probably A, since the closest points are all A. What I am looking for is some kind of work-flow that would allow me to do this in ArcGIS, if possible I would prefer R (if anyone has any experience in doing this kind of processes in R). This is just the example how my problem looks like, I still have to implement it to a bigger extent and to many layers of points. Bear in mind that I have little experience with interpolation and that this might be a easy/dumb question. Any help/guidelines are deeply appreciated.


title: Running Regular Queries Through ArcSDE

text: Working on an .net arcobjects app that's connecting to a Geodatabase (Oracle) through ArcSDE. I make a bunch of edits inside a EditSession but the data I am touching have a bunch of foreign Key constraints that I have to enforce by inserting/updating a bunch of records (non spatial) before I commit/flush/close my edit session. Is there a way to execute those queries through the same session being used by ArcSDE without having to deal with Table feature classes? Something that allow me to run pure SQL through an edit session?


title: Center a selected point on map. Openlayers + Geoext1 + ExtJs3.4

text: I've got an ExtJs combobox next to the map with several points of the map that are provided from a Postgis database through Php/XmlHttprequest. At the same time I've got a map with all points containing their respective attributes. I need to select a point from the list and then this point must be searched on the map an centered on it. Can somebody tell me which is the best way to proceed with this case? Relating combobox valuefield with map feature id maybe? Retrieving coordinates directly from database and centering the map on this coordinates? Thank you


title: How to calculate point from lat/long, angle, and distance?

text: I am trying to calculate the point from a given lat/long, a bearing, and a distance. This is for plotting weather given a aeronautical way point. For example: Given a waypoint with a known lat/long (LAX airport control tower), with a bearing of North by Northwest and a distance of 45 nautical miles, where is the point? I could treat the world as flat, but I think the resulting calculations would be inaccurate the further north the point is since longitudes get more narrow near the poles. I am mainly interested in plotting these points in the US. Is there a projection that would help me do this? This is for a Python application I'm writing, and the libraries I have access to are GDAL and Proj4 for projections. So no fancy tools like ArcGIS are available to me. Ultimately I want to combine 4 or these calculated points to form a polygon that should roughly be about 100-200 miles in width and height. Thanks!


title: How to estimate surveying territory

text: I am helping set up a project to gather survey data from small retail outlets in an emerging market using Open Data Kit. Is there any established methodology to estimate how much territory a surveyor can cover in a given time? I'm finding it difficult to estimate time and manpower requirements. Any help or links would be greatly appreciated.


title: How To Understand Geological Directions?

text: I was reading a "Guide to Some Minerals and Rocks in Indiana," which was published in the 1950s, and in it there were these bewildering directions: "Bluff along Mud Pine Creek 1.5 miles west of Rainsville, SW 1/ 4NE1/ 4sec. 29, T. 23 N., R. 8 W." It is on p.38 of this pdf: https://scholarworks.iu.edu/dspace/bitstream/handle/2022/235/C04.pdf?sequence=1 How may those directions be made understandable?


title: Select points within a relative number range

text: I need to select points that have attributes that are within a relative number range of each other. For instance; given the following numbers: 95,705,139,65,537,266,525 and 188; I need to find all of those numbers that are within a range of 50, 100 and 150 of each other. I then need to find those that are within a distance of each other. I am using QGIS. Any ideas gratfully accepted. Thanks.


title: Extracting Cities, Towns and Villages boundary coordinates

text: How can I extract the polygon coordinates of cities,towns and villages from OpenStreetMap? Are there any other free databases available?


title: GeoDjango Order by distance for Model that is linked to User who has a location on their user profile

text: I need to order my query model by distance, however the point field in question is not on the model class but on the user profile and the i am getting an error when attempting the following.               g = geocoders.Google(domain='maps.google.co.uk')     place, (lat, lng) = g.geocode('SW1')       MyModel.objects.distance(Point(lat, lng), field_name='user__userprofile__location').order_by('distance')      The error i get is "ST_Distance output only available on GeometryFields" Both MyModel and UserProfile have a GeoManager but obviously User does not. Filtering works fine using               MyModel.objects.filter(user__userprofile__location__distance_lte=(Point(lat, lng), D(mi=distance)))      But i cant understand why sorting is not working. Any help or workaround would be really appreciated.


title: In what data format should I export a (lat,long,alt) or (x,y,z) grid to be compatible with the largest number of GIS software?

text: I am new to GIS and developing an application that generates a 3D grid with (latitude, longitude, elevation) or (x, y, z) metric triplets, and I would like to add an exporting feature that ensures compatibility with the most GIS software. In case of several possibilities, I'd love to understand the pros and cons considering the current industry standards and also the "efficiency" of the format, as I'd like to reach the largest number of users. **UPDATE:** I did a bit of homework but I am still a bit lost with the variety of formats. What I have come accross which seem interesting: GRD (Surfer), KML (Google), DXF (AutoCad), WRL (seems cool but not used anymore), SHP (ArcView - closed format?) Anyone could help me choose?


title: Arcmap 10 feature with undefined coordinate system projects incorrectly when projecton is defined

text: I'm fairly new to GIS, and am slowly feeling my way through ArcMap 10. I have run across a problem defining the Geographic Coordinate System and projection for a feature (the feature currently has an Undefined GCS) that I created from a raster file using the Contour tool in Spatial Analyst. (Although the Undefined data lines up with other (defined) features that I have available, as I understand it, it is important to have all features within a map using the same GCS and projection? Is this correct, or can I work across multiple GCS?) I attempted to set the GCS to GDA_1994, and the projection to Lambert Conformal Conic (these are the GCS and projection of the other data sets that I have available). The data frames GCS is set as above (GDA_1994, Lambert). I also know the approximate location in Decimal Degrees (~ 150.9, -25.0 to 153.9, -28.5 => South East Queensland). When I attempt to use the 'Define Projection' tool (I believe this is the one I should be using, as 'Project' is only used when a GCS is already specified?) the scale and position of the data is completely thrown, with the entire feature fitting within 131.00, 0.00 DD (the measure tool gives the entire data set a size of ~ 3m x 3.5m - pretty much the same values as the true size of the data's range in DD - Coincidence?). Clearly I'm doing something wrong, or I'm missing a step that I have to do before using Define Projection, but I have no idea what..? Do I have to somehow set the datas units, or convert them from the angular units(DD) to linear units(m)? Any help would be appreciated.


title: get data from WFS and present in using d3js - correct projection

text: I'm trying to get data for map of Poland. I found WFS service that provides GML files.   It is located under: http://sdi.geoportal.gov.pl/wfs_prg/wfservice.aspx?request=getCapabilities&service=wfs   I'm able to all provinces of Poland using this request: http://sdi.geoportal.gov.pl/wfs_prg/wfservice.aspx?request=getFeature&version=1.1.0&service=WFS&typename=gmgml:WOJEWODZTWA My problem is that data from this server is in `EPSG:4326` projection, when d3js is using U.S.-centric (d3.geo.albersUsa).   I can convert gml to shp (using QGIS) and then to TopoJSON or GeoJSON (using http://www.mapshaper.org/), but I need to know what projection I must use to display my map correctly.


title: setSubsetString with limit?

text: I have a layer with couple thousands objects in QGIS and I would like to create a filter of the layer where result would be only one object.   I would like to do it in console or plugin. I am using `setSubsetString` function now. My current non-effective code looks like this:               self.vectorlayer.setSubsetString( "importance <= 10" )      and then I have for loop and select first loaded record and break. However, does exist any effective way? For instance, use a limit cause of other trick?


title: For a SDE FeatureClasses list all Roles that have been Granted any Privileges on it, and which Privileges each Role have on it

text: I have searched through the `ArcObjects 10 SDK` API and googled, but havnt been able to find what I am looking for, so thats the reason for me posting this question. **Q: How can I for a SDE FeatureClasses list all Roles that have been Granted any Privileges on it, and which Privileges each Role have on it?** The equivalent in normal _Oracle_ would be something like this:               SELECT dtp.PRIVILEGE,            dtp.grantable ADMIN,            dtp.grantee,            dtp.grantor,            dbo.owner,            dbo.object_type,            dbo.object_name,            dbo.status,            '' column_name       FROM ALL_TAB_PRIVS dtp, ALL_OBJECTS dbo      WHERE dtp.TABLE_SCHEMA = dbo.owner AND dtp.table_name = dbo.object_name            AND dbo.object_type IN ( 'TABLE', 'VIEW' )     UNION     SELECT dcp.PRIVILEGE,            dcp.grantable ADMIN,            dcp.grantee,            dcp.grantor,            dbo.owner,            dbo.object_type,            dbo.object_name,            dbo.status,            dcp.column_name       FROM ALL_COL_PRIVS dcp, ALL_OBJECTS dbo      WHERE dcp.TABLE_SCHEMA = dbo.owner AND dcp.table_name = dbo.object_name            AND dbo.object_type IN ( 'TABLE', 'VIEW' )      The query would list out all object of the types in the list provided (in the example I have listed `TABLE` and `VIEW`) which privileges have been granted on them like `SELECT`, `UPDATE` etc... and to which Roles. **Related Q:** SDE FeatureClass Granting and Revoking Privilege with ArcObjects


title: ArcGIS maps center to country

text: Is it possible to set center and zoom/fit bounds to country by name with ArcGIS Javascript API? I mean, with Google Maps you can do something like this: http://stackoverflow.com/a/6086437/1745228 How can we get data from ArcGIS and set extent on country?


title: How to save 'No data value' to a QGIS layer style properties?

text: I have many raster tif files (100+) I found this link to the MultiQml plugin very helpful, from this post: Can I change properties of several layers simultaneously? However the 'No data' values does not seem to be saved with the style? Is there a way to do this? Any suggestions are very welcome.


title: Leaflet's worldCopyJump breaks map, zooming fixes it up (updated)

text: I had my map working with `worldCopyJump` set to `true`, but then it stopped working. I'm not sure exactly what broke it (the original form of this question said there was a connection to a plugin, but that doesn't seem to be the case), but it might be that I updated to Leaflet 0.62. In any case, this is the problem: when the map loads, dragging the map would throw a bunch of errors claiming to have invalid LatLong objects (with a `NaN` instead of the longitude), and the map would freeze. Oddly enough, if _before dragging_ the map I zoomed in or out, and _then_ dragged, everything would work fine. Whatever the map is doing during `onzoom` seems to be somehow fixing whatever problem it had. Also peculiar is the LatLongs that are reported to be invalid: they are nowhere in the code. I have no idea where those numbers might be coming from. I traced the error back to the `worldCopyJump` option. If this is set to `false` then the map works from the get go. But I need that option turned on. Can anybody think of what might be going on? You can find a (non-)functional version of the map online here. In case that does not last, here's the code:                           <link rel="stylesheet" href="http://cdn.leafletjs.com/leaflet-0.6.2/leaflet.css" />


title: How can I delete a shapefile using OGR?

text: In my QGIS Python plugin, I use shapefile layers and I am having a problem rewriting existing shapefiles. Is it possible to delete a shapefile from disk using OGR?


title: Geoserver-WMS-problem when zooming point layer

text: I need help. I would like to have no blank(white) space between points when zooming in WMS geoserver. ![This is OK, whitout blank space between points at the smallest zoom](http://i.stack.imgur.com/xFEMf.png) ![When I zoom in, blank spaces appear, but I don't want that](http://i.stack.imgur.com/KK7UT.png) Any idea? Greetings, Andrew


title: Projection of Shapefiles for Spain (INE)

text: I downloaded the official shapefiles for the municipalities of Spain from INE here. It has a lot of shapefiles, but I found out that the right one for what I need is **esp_muni_0109.shp**. I would like to convert these shapefiles to **WGS84** with **ogr2ogr**. The problem is that these shapefiles have no .PRJ files, so I don't know from what projection system to convert them from. Does anyone know the right projection system for these shapefiles?   If so, can anyone tell me what parameters to use with ogr2ogr? (I'm pretty new in all of this). Thank you.


title: How to install Sextante in QGIS 1.8 without Plugin Installer?

text: How can I install sextante (and where can I download it) if the Plugin Installer doesn't works? As I work for a corporation, I have some problems conecting to the official repository, and I cannot solve them ...


title: Parsing field in field calculator using Python

text: I have field1 which has a values of "A_B" (example). I created two fields out of this, values to the left of the underscore go in fielda and to the right of the underscore go in fieldb. I used the code below in field calculator to separate them out into two fields. fielda =               field1.split("_")[-1]      fieldb =               `field1.split("_")[-1]`      I now have field1 with values of A_B_C_D. I need to separate the A_B_C_D into four separate fields using field calculator. fielda = A, fieldb = B, fieldc = C and so on. I'm not sure how to go about this. I basically need to parse the values by the underscore.


title: Field Calculator - replace string function

text: Within an attribute table I'm trying to replace one string with another using the 'replace' function . As an example, there are a number of features whose field 'NAME' contain the string '(B)' and I'm wanting to replace this with the word 'County'. I'm selecting the 'update existing field' checkbox and using the following expression                replace( '(B)','(B)','County')      The end result is that the field 'NAME' for every feature is replaced by 'County' irrespective of whether or not the field 'NAME' originally contained the string '(B)'. Any help would be appreciated.


title: What GIS applications and services have been affected by the Heartbleed vulnerability?

text: Heartbleed is the name given to a recent, very severe security vulnerability in the OpenSSL library, used in a large number of websites and applications. It is being called one of the worst computer security vulnerabilities in recent history. Even StackExchange was vulnerable as of April 8th (since resolved according to this answer). You can read more about it in these Security.SE questions:   * How exactly does the OpenSSL TLS heartbeat (Heartbleed) exploit work?   * What should a website operator do about the Heartbleed OpenSSL exploit?   * What clients are proven to be vulnerable to Heartbleed? What GIS applications and services have been affected by this vulnerability? Link to any official statements by the application maintainers about the vulnerability and/or how to resolve it. **Note to moderators:** Please mark this question and answers as community wiki, as this is more of a PSA than a Q&A.


title: Question about ERDAS TOC

text: On occasionally when I open a raster in ERDAS 2013, I get white or yellow circle icons on my ERDAS Contents and tried to find information on ERDAS. I want to know what does the icon mean ? Thanks appreicate it. Her is an attachment of what it looks like![Look I mark circle here](http://i.stack.imgur.com/rSJO2.jpg)


title: Animate a Choropleth Map in CartoDB?

text: I used Torque to animate points over time in CartoDB I have a map of polygons and I want to animate the changing colors. Over five years the intensity of the type of incident I'm mapping changes. I can make a slideshow that shows a visualization of each year in an individual choropleth map, but I'd love to be able to animate it in torque instead. Is this even possible?


title: Calculating scale from tiles using their envelope

text: I have a few problems I need to solve, and it comes down the the below example. It may be my understanding or the way I am doing it. Can someone please explain why my logic does not work. I am taking an XYZ parameters from an openlayers 3 client call, then creating a referenced envelope for this tile. I would then expect that if I calculated the scale from this envelope and the image size (plus assuming a dpi of 96) I would get an expected scale.               int x = 242784;     int y = 152342;     int z = 18;     double north = tile2lat(y, z);     double south = tile2lat(y + 1, z);     double west = tile2lon(x, z);     double east = tile2lon(x + 1, z);          ReferencedEnvelope envelope = new ReferencedEnvelope(new Envelope(west, east, north, south), DefaultGeographicCRS.WGS84);          double scale = RendererUtilities.calculateScale(envelope, 256, 256, 96);     System.out.println("scale = " + decimalFormat.format(scale));      So running the above I get a scale of 1988.682461 I was expecting something close to 2256.9988666883. Also different x,y gives vastly different scales, ie x=500, y=500 gives 197. I would have thought for the same zoom level the scales would all be the same. I'm pretty sure it's my understanding, so any help to clarify this would be much appreciated.


title: Binary insert into Postgis geography column results in value being inserted as geometry

text: When writing WKB as a value (and as a binary parameter) to a geography column, the persisted value is not geography. Example Java is also available here: https://github.com/ayuudee/issue- pad/blob/master/src/com/jesusthecat/im/pggeog/BinaryGeogTest.java **Here's what it does:**   1. Creates a table with an ID and Geography(Point, 4326)   2. inserts a row using WKT.   3. inserts a row using WKB (as bytes).    4. Prints out id, point, and ST_SUMMARY(point). **SQL result of step #4 is:**                id |                         pt                         | st_summary      ----+----------------------------------------------------+------------       1 | 0101000020E610000009C6C1A5E3E662406BB75D68AEED40C0 | Point[GS]       2 | 0101000020E610000009C6C1A5E3E662406BB75D68AEED40C0 | Point[S]      **Example Log of the binary insert (where pt is Geography(Point, 4326))** > LOG: execute : insert into px(pt) values($1) > > DETAIL: parameters: $1 = > '\x0101000020e610000009c6c1a5e3e662406bb75d68aeed40c0' **You'll notice that:**   1. The WKB for both #1 and #2 are the same; and    2. The flags in the ST_Summary result for #1 are [GS], where for #2 (the binary) they're [S]. The documentation for ST_Summary would indicate that #2 has spatial information, but is not geodetic (i.e. not geography). I'm writing a Java library that seeks to persist geography as binary, but this would seem to indicate that it isn't possible. Also, is it normal that this should be allowed to happen anyway (i.e. writing a value to a geography column that's not geography)? **Versions:** POSTGIS="2.1.3 r12547" GEOS="3.3.3-CAPI-1.7.4" PROJ="Rel. 4.7.1, 23 September 2009" GDAL="GDAL 1.9.0, released 2011/12/29" LIBXML="2.7.8" LIBJSON="UNKNOWN" TOPOLOGY RASTER POSTGRES PostgreSQL 9.3.5 on x86_64-unknown-linux-gnu, compiled by gcc (Ubuntu/Linaro? 4.6.3-1ubuntu5) 4.6.3, 64-bit Would appreciate any help.


title: How to select point feature rendered by picture

text: I use GeoServer as my gis service and render some points with pictures, SLD code like this:               <PointSymbolizer>         <Graphic>            <ExternalGraphic>              <OnlineResource xlink:type="simple" xlink:href="http://xxx/xxx.jpg" />              <Format>image/jpeg</Format>            </ExternalGraphic>            <Size>32</Size>         </Graphic>     </PointSymbolizer>      and in openlayers all work fine,i can see picture is rendered normally. But problem is when i want to click(select) this picture , it seems openlayers using WFS/WMS generated click bbox coordinates and send them to geoserver, and if this bbox intersects with the Point(picture) coordinates,the point is selected only if i click the very close to point coordinate, but when i click the edge of this picture or not very close to the point(still inside this picture), the click query bbox generated by openlayers do not intersect with picture rendering point and this feature will not be selected finally. Simplly said: there is a picture which 256*256 size, rendered in geoserver by PointSymbolizer,and this point's coordinates is [128,128](screen coordinate,center of picture), when you click close to 128,128 this picture can be selected, but when you click [250,250](still inside picture)this point cant be selected. Any advice ? thanks


title: How do you remove the title in the banner from the esri geoportal extension?

text: I would like to remove the catalog.site.title text from the banner area in the geoportal page. However, I would still like to keep the title of the page in the browser bar/window. Does anyone know in which file the catalog.site.title is located for just the banner text? If I remove the text from the catalog.site.title from the site resources in the gpt.properties file the banner title and browser title are replaced with `???catalog.site.title`. Thank you in advance. ![Title to remove and title to keep](http://i.stack.imgur.com/nYHOQ.jpg)


title: Issue blocking QGIS feature form dialog with Python

text: I'm not sure where I'm going wrong with this python logic for a custom form in QGIS 2.0.1. There are a few things required by the validation function: that index of combo boxes aren't 0 (first index is empty), the line edits aren't empty or NULL, or Autre. It's very puzzling, because blocking validation works except for when it comes to the last 2 or 3 items, the dialog accepts even when the if statements are true. When some of the form is isolated, it works:               def formOpen(dialog, layerid, featureid):         global myDialog         myDialog = dialog         global Naturefield         Naturefield = dialog.findChild(QComboBox, "Nature")         buttonBox.accepted.disconnect(myDialog.accept)         buttonBox.accepted.connect(validate)         buttonBox.rejected.connect(myDialog.reject)          def validate():         global Naturefield         #intending to block dialog accept if 'Autre' was not changed         if Naturefield.currentText() == 'Autre':             msgBox = QMessageBox()             msgBox.setText("Define autre")             msgBox.exec_()         else:             myDialog.accept()      However, it seems as though I run into problems when too many conditions are required. Entire thing:               from PyQt4.QtCore import *     from PyQt4.QtGui import *               Diametrefield = None     Naturefield = None     PressCons1field = None     PressCons2field = None     SensFermfield = None     SensFluidefield = None     Fournsrfield = None     Commentfield = None     OriginInfofield = None     Nat_autrefield = None     Four_autrefield = None     myDialog = None          def formOpen(dialog,layerid,featureid):         global myDialog         myDialog = dialog         global Diametrefield, Naturefield, PressCons1field, PressCons2field, SensFermfield, SensFluidefield, Fournsrfield, Commentfield, OriginInfofield, Nat_autrefield, Four_autrefield         Diametrefield = dialog.findChild(QComboBox,"Diametre")         Naturefield = dialog.findChild(QComboBox,"Nature")         PressCons1field = dialog.findChild(QLineEdit, "PressCons1")         PressCons2field = dialog.findChild(QLineEdit, "PressCons2")         SensFermfield = dialog.findChild(QLineEdit, "SensFerm")         SensFluidefield = dialog.findChild(QLineEdit, "SensFluide")         Fournsrfield = dialog.findChild(QComboBox,"Fournsr")         Commentfield = dialog.findChild(QLineEdit,"Comment")         Nat_autrefield = dialog.findChild(QLineEdit,"Nat_autre")         Four_autrefield = dialog.findChild(QLineEdit,"Four_autre")         OriginInfofield = dialog.findChild(QLineEdit,"OriginInfo")              Naturefield.currentIndexChanged.connect(validate)         Fournsrfield.currentIndexChanged.connect(validate)              buttonBox = dialog.findChild(QDialogButtonBox,"buttonBox")         # Disconnect the signal that QGIS has wired up for the dialog to the button box.         buttonBox.accepted.disconnect(myDialog.accept)              # Wire up our own signals.         buttonBox.accepted.connect(validate)         buttonBox.rejected.connect(myDialog.reject)               def validate():         global Diametrefield, Naturefield, PressCons1field, PressCons2field, SensFermfield, SensFluidefield, Fournsrfield, Commentfield, OriginInfofield, Nat_autrefield, Four_autrefield         # Make sure that the Nature field isn't empty.         if Fournsrfield.currentText() == '..Autre':             QMessageBox.warning(None,'Probleme de saisie',"Definir 'Autre' Fournisseur")             Fournsrfield.setLineEdit(Four_autrefield)         if Naturefield.currentText() == '..Autre':             QMessageBox.warning(None,'Probleme de saisie',"Definir 'Autre' Nature")             Naturefield.setLineEdit(Nat_autrefield)         if Diametrefield.currentIndex() == 0:             Diametrefield.setStyleSheet("background-color: rgba(255, 107, 107, 150);")             QMessageBox.warning(None,'Probleme de saisie',"Veuillez choisir 'Diametre'")         if Naturefield.currentIndex() == 0:             Naturefield.setStyleSheet("background-color: rgba(255, 107, 107, 150);")             QMessageBox.warning(None,'Probleme de saisie',"Veuillez choisir 'Nature'")         if not len(PressCons1field.text()) > 0 or PressCons1field.text() == 'NULL':             PressCons1field.setStyleSheet("background-color: rgba(255, 107, 107, 150);")             QMessageBox.warning(None,'Probleme de saisie',"Veuillez remplir 'Pression de Consigne 1'")         if not len(PressCons2field.text()) > 0 or PressCons2field.text() == 'NULL':             PressCons2field.setStyleSheet("background-color: rgba(255, 107, 107, 150);")             QMessageBox.warning(None,'Probleme de saisie',"Veuillez remplir 'Pression de Consigne 2'")         if not len(SensFermfield.text()) > 0 or SensFermfield.text() == 'NULL':             SensFermfield.setStyleSheet("background-color: rgba(255, 107, 107, 150);")             QMessageBox.warning(None,'Probleme de saisie',"Veuillez remplir 'Sens de Fermeture'")         if not len(SensFluidefield.text()) > 0 or SensFluidefield.text() == 'NULL':             SensFluidefield.setStyleSheet("background-color: rgba(255, 107, 107, 150);")             QMessageBox.warning(None,'Probleme de saisie',"Veuillez remplir 'Sens Fluide'")         if Fournsrfield.currentIndex() == 0 or Fournsrfield.currentText() == '..Autre':             Fournsrfield.setStyleSheet("background-color: rgba(255, 107, 107, 150);")             QMessageBox.warning(None,'Probleme de saisie',"Veuillez choisir 'Fournisseur'")         if not len(Commentfield.text()) > 0 or Commentfield.text() == 'NULL':             Commentfield.setStyleSheet("background-color: rgba(255, 107, 107, 150);")             QMessageBox.warning(None,'Probleme de saisie',"Veuillez remplir 'Commentaire'")         if not len(OriginInfofield.text()) > 0 or OriginInfofield.text() == 'NULL':             OriginInfofield.setStyleSheet("background-color: rgba(255, 107, 107, 150);")             QMessageBox.warning(None,'Probleme de saisie',"Veuillez remplir 'Origin d'Information'")         else:             # Return the form as accepted to QGIS.             myDialog.accept()      I've now attempted to clear up the signal connect by adding another function for the "autre" requirement. This has cleared up MOST of my problems, but the last few validation requirements still slip through the cracks.               def formOpen(dialog, layerid, featureid):         Naturefield.currentIndexChanged.connect(autre)         Fournsrfield.currentIndexChanged.connect(autre)         Typefield.currentIndexChanged.connect(autre)          def autre():         if Diametrefield.currentText() == '..Autre':             QMessageBox.warning(None,'Note de saisie',"Definir 'Autre' Diametre")             Diametrefield.setLineEdit(Diam_autrefield)         else:             Diametrefield.editTextChanged.connect(validate)          def validate():         if Diametrefield.currentIndex() == 0 or Diametrefield.currentText() == '..Autre':             Diametrefield.setStyleSheet("background-color: rgba(255, 107, 107, 150);")             QMessageBox.warning(None, 'Probleme de saisie',"Veuillez choisir 'Diametre'")


title: Determine when I have crossed state line?

text: We are currently using the Google Maps API to get route information and total mileage but we would like to determine when we are crossing into a new state. When we return the route information we can be handed 2000+ Latitude and Longitude points. Initial thought was to do a lookup on each point to determine what state it was in but that would be very intensive. Another option would be to try and find a way to determine if the Lat & Long where inside a polygon or KMZ. Has anyone seen code for this? Otherwise what other options should I be investigating?


title: strange behaviour of piechart function of QGIS

text: I am using QGIS to visual health status of elderly in an elderly home, whether they have medical devices or not is one attribute that I am trying to visualize. I have create the map using `mapwindows`, and use `R` to create a tab- delimited txt file to carry the medical device information (7 columns, 6 columns of `0` or `1`, and one uid that will be used to join the bed layer in QGIS. I used the diagram tab in the properties menu of the bed layer, and pie-chart is being used for visualization: if the elderly have 3 devices, the pie will be presented with 3 parts with 3 different colours. And strange things are visualized. A majority of elderly do not have any devices (i.e. that should be straight zeros for the 6 column in the txt file). I support no pie chart should be visualized, but QGIS still show a solid pie using the last attribute I have inserted in the diagram tab of the properties menu. I have tried to insert a `csvt` file to explicitly state that those medical device columns are `integer`, `real`, but it still fails. QGIS visual an attribute with value 0 as 1, if all of the attributes used for plotting are straight zeros. Can anyone help? Thanks.


Does the a wms layer in openlayers have "maxfeatures" option? I'm asking this because my wms layer contains allot of points ( > 100.000). When zoomed in I don't want the rest of the features (not sure if you call it features when using wms) to be rendered. When I was using WFS, this was my code:                   wfs = new OpenLayers.Layer.Vector("WFS", {             strategies: [new OpenLayers.Strategy.BBOX({                             resFactor: 1,                             ratio:1                         })],             protocol: new OpenLayers.Protocol.WFS({                             maxFeatures:1000,                             url:  "my url",                             featureType: "my_featuretype",                             featureNS: "my_feature ns",                             version: "1.1.0"                         })     });      How can I accomplish the same render functionality with a wms layer? **Edit** Solution: Hide the layer at a certain zoomlevel. This prevents a big amount of features to be loaded. code:                   map.events.on({ "zoomend": function (e) {             if (this.getZoom() > 7) {                 layer.setVisibility(true);             }             else {                 layer.setVisibility(false);             }         }     });


title: Where to look for GIS Internships for Summer 2012

text: I am currently looking internships for this summer mainly in the Washington, D.C. areas (although, any area would be alright). Currently I am getting my masters in GIS and have a B.S. in anthropology. I have fair amount of background in Python and SQL. Does anyone know for a good place to look for internships besides usajobs.gov and indeed.com ?


title: Need help calculating forest cover within a buffer when it overlaps with other buffers

text: I'm new to ArcGIS and having real trouble getting this done… Hopefully someone can help. Thanks in advance! I have a layer showing forest cover in a landscape. I'm trying to create a 5km buffer around every forest patch in the landscape and then calculate the amount of forest within that buffer. I've done this before with non- overlapping buffers, just by dissolving all forest patches into a single multipart feature and then intersecting this with the buffers. However, in the landscape I am working on now many of the buffers overlap with each other, and wherever a forest patch within the buffer of interest is also intersected by another buffer it gets split into several features. So, instead of a single value for forest cover per buffer I'm getting value for every little piece of patch that has been created by the overlaps. Does anyone know how to get around this, so that I have a single value of forest cover within each buffer? Seems like it should be easy, but I've been fiddling with spatial join and merge all day but haven't gotten anywhere. Using V10. Cheers, Jay


title: How to create Network Dataset for Network Assistant using ArcPy?

text: Is it possible to create a network dataset with ArcPy to be used with the network assistant? In the documentation I found the **Build Network** and **Dissolve Network** functions, but it looks like both require an already existing network dataset as input parameter. In this **thread** people say it is not possible with Python but with **.NET**. In this **thread**, somebody says it is possible to create geometric networks with arcpy in **10.1**. Is it possible to use a geometric network with the network assistant? In the documentation it says: Only point and line feature classes are supported as input. The feature classes cannot participate in another geometric network or other advanced geodatabase dataset such as a topology, network dataset, terrain, or fabric.


title: Google Map Projection Skew - Need correct Projection

text: I'm having a really difficult time resolving skewed grids when projecting onto Google maps using Alteryx outputs. Within the module, I'm taking GRD files and aggregating rows and columns to re-size from 30m to 250m grids. Then projecting using EPSG:3857 which I understand to be the correct Web Mercator Aux Sphere for Google/Bing. However, no matter what I research and attempt to implement, I always end up with skew in my squares. Any suggestions would be very helpful. Thanks ![Alteryx Process Flow](http://i.stack.imgur.com/TVmQU.jpg) ![enter image description here](http://i.stack.imgur.com/9iFc6.jpg)


title: How to use buffer() in QGIS Value Relation field

text: When adding a point, I'd like to auto-populate one of its fields (the street address) setting it to a field in the nearest point in another layer (a layer containing all existing addresses). It looks like the Value Relation edit widget could do the trick, if I could figure out how to use the `buffer()` function. Right now, I can get the list of all addresses by selecting the `Layer`, `Key` column and `Value` column. I thought I could filter that list so that it has only nearby points (ideally just one), using something like: `within($geometry_A, buffer($geometry_B,50)) = 1` where `$geometry_A` is the new point and `$geometry_B` is the existing point. `$x`, `$y` and `$geometry` seem to refer to the existing points in the selected `Layer`. Assuming I am on the right track, how do I refer to the geometry of the new point (i.e. `$geometry_A`)? ![screenshot](http://i.stack.imgur.com/rcwGv.png)


title: Flood modelling insights

text: I am trying to create a flood modelling system. The system consists of two parts mainly a river flooding simulation and rainfall simulation. To do this I am running the latest version of GRASS 6.4. Also I am using the SRTM obtained from CGIAR as my DEM. I would like to ask what are the most appropriate commands to run a simulation. Ex. Terraflow Also since we are modelling flooding based on a river, is it advisable to use HEC-RAS? It was created by the same folks who made GRASS. :) Sorry again for the dumb question. I am more of a programmer rather than an engineer. It is my first venture into simulation systems as part of my degree. Thanks for your insights. BTW as of now Here is our plan for the rainfall (pseudocode):               int i = 0, j = 0;     int passes, numberOfmesh;     double rainAmount;     double store1[][];     double store2[][];          numberOfmesh = getTotalMeshCells();     passes = computePasses(numberOfmesh());     rainAmount = getRainAmount();          /* rainAmount corresponds to the input of the user in the textfield.     store1 contains the elevation and the number of times an amount of water     has been passed to it (per pixel)     store2 contains the elevation and the total accumulated water (per pixel)     Both store1 and store2 have 'n' number of rows and 2 columns;     ASSUME THAT THE ARRAY HAS ALREADY BEEN POPULATED */          while(i != passes)     {         while(j != numberOfmesh)         {             direct = getDirection(store[j][0]);             /* direct corresponds to the direction of a mesh cell;             I don't know what id the datatype for direct; */                         store2[j][6] += rainAmount;             gotoDirection(direct);              }        }     computeAccumulated();          //Here is what gotoDirect is supposed to do:          gotoDirect(direct)     {        /* How can we determine the elevation of the mesh cell to where the          'direct' points to?          In the following code, j refers to the elevation of the mesh cell where          the 'direct' point to.        */             store1[j][7] += 1;        /* increases the number of times that an amount of water        was passed to it */     }          //Here is what computeAccumulated is supposed to do:          computeAccumulated()     {         while(j != numberOfmesh)         {             store2[j][8] += store1[j][9] * rainAmount;                     }                    }


title: How to import raster data in GPS with Arc GIS

text: How do I see the map on GPS as a background layer? In QGIS, I could see the map which is converted from shape file & raster data on QGIS with "Garmin Custom Map". I need to do the similar thing through Arc GIS. GIS: Arc info 10.1 GPS: Garmin GPS map 62sc Any comment or advice is very helpful. Best Regards,


title: How to create a single contour defining the boundary between sea and land

text: Hello!   I am still learning so forgive me if my question is very simple... I am trying to create a map with ArcGIS. The area I am working on is a seasite and I found that in some parts the cells appear as white. Do you have any idea how I can fix it? Moreover, I would like to create a single contour along the seacoast just to highlight the boundary between the land and the sea. Does anybody know how I can do it? Thank you in advance! :)


title: How to save R raster plot for ArcGIS raster file

text: I want to save my raster plot that produced with R to ArcGIS raster (.img, .tif and ASCII) file.


title: Using kriging with data in GDB

text: I have series of data in my GDB I am running the program under VBA and what I want to do is using the data by the tool Kriging. This is what I thought: First, Kriging needs data that must be a point shape file and I am getting the data out of my GDB (one point in one table but there is a lot of data)   1. I get the data I need and make it in a shape file   2. I use the Kriging tool to calculate the shape file   3. I use the raster to point tool to change the raster Kriging make to a point layer   4. I get the point data I want! Here is my question: **When I get the data from the GDB (using sql), how to change it to a shapefile in VBA?** I have only found the page to display XY data using VBA: http://edndoc.esri.com/arcobjects/9.0/samples/tables/create_a_layer_based_on_xy_data.htm but it's can't save as shape file! How to save it automatically? Is there any wrong with my steps? If need more information please tell me! @artwork21 An error occurred in the line gp.FeatureClassToShapefile_conversion "'C:\Temp\New File Geodatabase.gdb\point'", Temp_ _2_ says automation error unspecified error how to fix it? Call tools from ArcToolBox in VBA form in this post about my problem I fix it like follow gp.Kriging "C:\temp\123.shp", "F30", "C:\AA\Kriging_123", "Spherical 249.101224", "250" but this time no matter what I try it's don't work....


title: How do I parse complex json?

text: My json response:               {"type":"FeatureCollection", "features":[{"type":"Feature","id":"SPATIAL_VIEW.fid-54038ee6_13f98ee9f52_4049","geometry":{"type":"Polygon","coordinates":[[[7.2970943E7,-4.69048595E8]]]},"geometry_name":"GEOM","properties":{"ID":"123","GIS_CODE":"AP","NAME":"Andhra Pradesh"}}],"crs":{"type":"EPSG","properties":{"code":"42106"}}}      In this I want to get the GIS_CODE as an alert. What I tried is: That json response is in "jsoncontent"                  var content = eval('('+jsoncontent+')');        alert(content.features.properties.GIS_CODE);      I am getting this error Uncaught TypeError: Cannot read property 'GIS_CODE' of undefined How to resolve it? Any help is greatly appreciated .. Thank you


title: Order for inserting data into postgis

text: I am using postgis for my project along with geoserver and openlayers. I give user an html form with drop down menu list of available layers(each layer corresponds to a table in db).The user can select any layer , based on which I am retrieving the fields of that particular from db to insert the values. As soon as the user inputs data and clicks on submit then the data should be inserted into db. Inserting data into db is a trivial thing as far as you know the order in which to insert into the db. But I have tables with different number of fields and in different order. I want to know how can I insert values into particular table when I don't know the order in which I have to store. I shall explain with an example: Table 1 of layer1: id,name,point Table 2 of layer2: id,polygon,name,area When user selects layer1 in dropdown menu , the form will provide you id,name,point text fields in that order to insert data. If the user selects layer2 then form will give you id,polygon,name,area text fields in that order. Suppose say user has selected layer1 then after inserting the values into text fields and submitting , the php file should perform a query "insert into (id,name,point) values('$_post[id]','$_post[name]','st_geomfromtext('POINT($_post[point])',4326));" Same way for inserting into layer2 the query should be "insert into (id,polygon,name,area) values('$_post[id]','st_geomfromtext('POINT($_post[point])',4326),'$_post[name]','$_post[area]');" But how can I know that I have to use only this particular query??? I can know the table name from hidden field. For now my php file is like this: $tablename=$_POST[tablename]; if($tablename=='layer1') { perform query1; } else if($tablename=='layer2') { perform query2; } But this works only if I know the tables already present. Is there anyway to write the php file provided I know only the table name but no idea of table description in db?? (I can get the table description through a query if I want).


title: Visualization of 3D shapefiles with Mapserver and Postgis

text: I'm wondering is there a way to visualize 3d data with mapserver? 3D .shp files are in postgis database, and I already connected mapserver to database, but with simple layer in mapfile.


title: What is Lanczos resampling useful for in a spatial context?

text: GDAL includes a resampling method beyond the normal mix of nearest neighbor, bilinear, cubic and splines: "Lanczos windowed sinc resampling". I understand that its a convolution filter, but unlike images where results tend to be subjective, the resampling used for spatial data has other implications. What's Lanczos and how does using it affect the output?


title: how to overlay GPX data on rasters in QGIS

text: > **Possible Duplicate:**   >  Overlaying GPX data on rasters in QGIS When both layers are displayed in QGIS, they are not even visible in the same display. The coordinates displayed for each are wildly different and are not of the same magnitude. how to solve that problem?


title: Replacing NULL values in geodatabase table?

text: How do you update NULL values within a geodatabase table using the ArcGIS Field Calculator (service pack 4) without changing other values to NULL? I've tried the suggestions on this previous question: How do you test for a NULL string field in the Field Calculator with Python? When I run the script it changes everything to NULL.               def RemoveNULL(x):         if x is None:             return '0'      * * *               RemoveNULL(str(!myxField!))      e.g. **before running script**               1     23     NULL     5      **after running script**               NULL     NULL     NULL     NULL      Can you not use python script on mdb table?


title: What does map scale refer to in map making / field survey metadata?

text: If I have, say, a printed map at the scale of 1:20000 and the legend says that it was mapped at a scale of 1:10000 does this simply mean, that the surveyors/cartographers made drafts & sketches or used aerial photographs at that scale or can there be some other meaning? Also, if the metadata of vector data say that its scale/resolution is e.g. 1:5000, can it have any other meaning besides having been digitized from a source (aerial photograph, printed map) of that scale? **edit** See my comment below for specific context (general answers are still fine).


title: TileStache:TileStache.Core.KnownUnknown: Couldn't get a layer from data source!

text: I want to use tilestache to get geojson files. I hava run tilestache-server.py -c tilestache.cfg ,and it works well. However,I want to get tiles from data stored in postgis. Maybe it will need to run tilestache-seed.py. So I run`               tilestache-seed.py -c ~/TileStache-1.49.8/cfg1 -b 32.000 114.000 31.000 115.000 -l postgis 12  -e geojson      then I got               Traceback (most recent call last):      File "/usr/local/bin/tilestache-seed.py", line 331, in mimetype, content = getTile(layer, coord, extension, options.ignore_cached) File "/usr/local/lib/python2.7/dist-packages/TileStache/ **init**.py", line 64, in getTile status_code, headers, body = layer.getTileResponse(coord, extension, ignore_cached) File "/usr/local/lib/python2.7/dist- packages/TileStache/Core.py", line 414, in getTileResponse tile = self.render(coord, format) File "/usr/local/lib/python2.7/dist- packages/TileStache/Core.py", line 500, in render tile = provider.renderTile(width, height, srs, coord) File "/usr/local/lib/python2.7/dist-packages/TileStache/Vector/ **init**.py", line 583, in renderTile layer, ds = _open_layer(self.driver, self.parameters, self.layer.config.dirpath) File "/usr/local/lib/python2.7/dist- packages/TileStache/Vector/ **init**.py", line 453, in _open_layer raise KnownUnknown('Couldn\'t get a layer from data source %s' % source_name) TileStache.Core.KnownUnknown: Couldn't get a layer from data source PG:dbname='osmgis' user='postgres' host='127.0.0.1' password='xxxxxx'port='5432' here is my cfg1               {      "cache": { "name": "Disk", "path": "~/TileStache-1.49.8/stache", "umask": "0000" }, "layers": { "postgis": { "provider": { "name": "vector", "driver": "PostgreSQL", "parameters": { "dbname": "osmgis", "user": "postgres", "password":"xxxxxx", "host": "127.0.0.1", "port": 5432, "query": "SELECT ST_AsGeojson(geom) from wuhan_osm_roads"}} } I am a beginner ,and I don't know how to solve the problem. by the way, if i delete "12" behind "postgis" in code,it works but output nothing.


title: Where is list of Geographic Transformations in ArcGIS, and their areas of applicability?

text: I have some data in Provisional South American Datum 1956, and I need to project it to Web Mercator with WGS84 as the datum. When I use the project tool in ArcGIS desktop, it requires a Geographic Transformation, and offers me 13 different ones: from `PSAD_1956_to_WGS_1984_1` to `PSAD_1956_to_WGS_1984_13` How do I find out which is the transformation meant for my area? I remember finding a pdf on the Esri Forums several years ago, which showed the various transformations between NAD 83/27 & WGS 84. Is there something similar for all transformations in ArcGIS?


title: Fastest way to create an ESRI Grid and populate cell values?

text: In terms of processing speed what's the fastest way of creating and then adding values to a raster? The values will originally come from ASCII files which can be about 500MB in size - ideally I would like to create an Esri GRID but that is not essential. I have Python (GDAL) and/or ArcOjects (9.3.1) in my arsenal I cannot use the ASCII to Raster tool due to the format of the values in the ASCII file . I'm looking for a solution of beaing able to write values to the raster directly. As an example I may have an XYZ asccii as like:                   @File_Version: 3     @Coordinate_Type_is: 3     @Export_Type_is: 1     @Number_of_Projects 1     @Project_Type_Name: 3 cstone25      @Project_Unit_is: meters     #File_Version____________-> 3     #Project_Name____________-> cstone25     #Project_Type____________-> 3     #Project_Unit____________-> meters     #OpenWorks_Project_______-> EXPLORATION_UTM31     #Master_Project_______->      #Coordinate_type_________-> 3     #Number_of_points_in_hzd_-> 1     zpxc1_ufa_nra_env_sele+20+100_Max_Peak_Amp     #Horizon_internal_id_____->_120     #Horizon_extremes_are____->   7245798912.00000     #Horizon_onset_is__Minimum_____-> 1     #Horizon_type_is_Time______-> 1     #Horizon_color_is________-> 0 255 0      6148.0 17392.0     529940.17623    6274994.18251 394123616.00000      6148.0 17396.0     529940.83286    6275019.17388 326897280.00000      6148.0 17400.0     529941.48949    6275044.16526 567459712.00000      6148.0 17404.0     529942.14613    6275069.15663 429482368.00000      6148.0 17408.0     529942.80276    6275094.14801 133182280.00000      6148.0 17412.0     529943.45939    6275119.13938 32126824.00000      6148.0 17416.0     529944.11603    6275144.13076 17940496.00000      6148.0 17420.0     529944.77266    6275169.12213 320429248.00000      6148.0 17424.0     529945.42929    6275194.11351 453848192.00000      6148.0 17428.0     529946.08592    6275219.10488 336978496.00000      6148.0 17432.0     529946.74256    6275244.09626 290245120.00000      6148.0 17436.0     529947.39919    6275269.08763 323590976.00000      6148.0 17440.0     529948.05582    6275294.07901 249979616.00000      6148.0 17444.0     529948.71246    6275319.07038 221433360.00000      6148.0 17448.0     529949.36909    6275344.06176 244277504.00000      or I might have a header than values:               !     !     FILE NAME :  Rob_Grid                     !     FORMATTED FILE CREATION DATE:  APR 25 2011     !     FORMATTED FILE CREATION TIME:  17:11     !     @Rob_Grid HEADER                         ,      GRID,         5           15,  0.1000000E+31,            ,      7,      1          245,    232,   262410.0    ,   273960.0    ,   6426990.    ,   6439190.             200.0000    ,   0.000000    ,   0.000000         @       0.1000000E+31  0.1000000E+31  0.1000000E+31  0.1000000E+31  0.1000000E+31       0.1000000E+31  0.1000000E+31  0.1000000E+31  0.1000000E+31  0.1000000E+31       0.1000000E+31  0.1000000E+31  0.1000000E+31   1023.936       1021.771            1020.476       1019.986       1019.994       1020.656       1022.158            1024.154       1026.193       1027.803       1028.656       1028.808            1028.548       1028.080       1027.858       1028.684       1030.780          **NOTE: this question is about how to create and populate the raster values the fastest** Thanks Rob


title: How to limit Create Fishnet extent by mask & carry over attributes?

text: I have a 5km polygon grid of an area which is not a standard square/rectangle, it is in fact a rough ovoid made up of around 500 5km2 cells. What I want to do, eventually, is have a 100m polygon grid covering the entire area that is currently represented by the ovoid(ish) 5km grid - about 1.25 million cells. The caveat is that the 100m polygons should be linked to a 'parent' 1km grid square: i.e. you have 12,500 1km cells with a given ID (1KM_ID) and then the 100m cells are given an ID too (100M_ID) which goes from 1 to 100 for each 1km cell. Using Create Fishnet will create 1km cells based on the bounding box of the 5km cell ovoid, which will greatly increase processing time, i just want the fishnet based on the cells that exist. Same applies for then creating 100m cells from the 1km cells. This can be done really nicely in FME but for some reason there are processing issues and duplicates are created around the edge of each 1km cell, messing the whole thing up. Any ideas for AcrMap?


title: Loop over features of a GeoJSON object from CartoDB

text:
I have just little knowledge about HTML and JavaScript. Nevertheless I was playing, copying and pasting around with Leaflet and CartoDB. I want to display the features GeoJson in a table. This is what I tried.               var map = L.map('map').setView([51.0, 10.0], 6);          var base = {     "S/W" : L.TileLayer.provider('OpenStreetMap.BlackAndWhite'),     "TOP":  L.TileLayer.provider('Thunderforest.OpenCycleMap'),     "OpenStreetMap": L.TileLayer.provider('OpenStreetMap.Mapnik')     };               window.backupAlert = window.alert;     window.alert = function () {return true};     var alle = L.cartoDB('http://leak-my-brain.cartodb.com/api/v2/sql?q=SELECT * FROM jobsnew', {     pointToLayer: function(feature, latlng) {     return L.marker(latlng, {     icon: L.icon({     iconSize: [4, 4],     iconAnchor: [0, 0],     popupAnchor: [0, 0],     iconUrl: "http://www.zoowelt.org/wp-content/uploads/2013/02/viv_point.png"     })     })     },      onEachFeature: function(feature, layer) {     layer.bindPopup('' + feature.properties.title + '
' + feature.properties.kw);     }     });          map.addLayer(base.OpenStreetMap);      map.addLayer(alle);               map.on( "zoomend", function( e ) {         console.log( "zoom level is " + map.getZoom())     }); //zur ausgabe des Zoomlevels in der Konsole     var layersControl = L.control.layers(base,alle);      layersControl.addTo(map);       for (var key in alle._layers) {       if (alle._layers.hasOwnProperty(key)) {         backupAlert("hi comment");       }     };     function writeTable() {         var out = '';         for (var row=0; row<1000; row++) {         if ( typeof(alle._layers[row])=='undefined') continue;             out += '
'              + '' + alle._layers[row].feature.properties.title + '	' +              + '' + alle._layers[row].feature.properties.kw + '	'              + '
'          }                document.write(out);      };      This somehow works but there are some issues. The main problem is that the indices in the `alle` are not starting at `0` and they are not increasing by the step size of `1`. ![](http://www.digital-geography.com/wp- content/uploads/2013/05/leaflet_jobs.png)


title: "Add Delimited Text Layer" plugin is enabled, but does not show up in the plugin menu

text: I'm using QGIS 1.7.1 on Windows 7. I want to import a delimited text layer, but that option does not appear in the "plugins" menu. This problem persists, despite the following steps I've taken to remedy it:   1. Disable then Enable "Add Delimited Text Layer" plugin    2. Restart QGIS/Restart computer    3. Uninstall QGIS and GRASS using add/remove programs, then re-install  Is there any way to completely remove all plugins, then re-install and re- enable "Add Delimited Text Layer?" Thank you ![Plugins](http://i.stack.imgur.com/dPRYm.png)


title: Software program (or spreadsheet calculations) to find Borehole Deviation?

text: I was given a bunch of surveyed borehole data (azimuth, depth, dip, geology, etc) and am required to calculate a deviation variable. Besides querying the data and averaging the surveyed borehole results, I was wondering if anyone had any specific calculations or recommended software programs in order to calculate a deviation variable and to visualize it in a 3D view around a planned borehole?


title: What is the difference between mosaicking and merging?

text: What is the difference betweem image Mosaicking and Merging?


title: Applications of GIS in geomorphology

text: I'm taking a geomorphology class right now that is theoretical, learning about geomorphic processes and formations. There is no computer or quantitative component. I'm enjoying it as far as it goes, but I'm left wondering: what does GIS get used for in modern geomorphology. What kind of data gets collected? What do people do with the data?


title: Getting/downloading shapefile from OpenLayers/Web page and saving to local disk

text: I have developed a web application where i am using GeoServer and OpenLayers. I have displayed some Layers over Google map in a ASP.Net Web Page. I want to know is there any Method or Script available in OpenLayers for downloading and saving shapefile (the shapefile which is displayed over Google map in Web page) from web interface. Actually the user wants to download the shapefile from Web interface and open it in his Dektop GIS (ArcView) system.


title: Can QGIS read ESRI TIN file directly?

text: Can QGIS open directly a TIN file generated with ArcGIS?


title: How to return data layer using geoserver rest api?

text: I'm starting to use the Geoserver Rest Api, testing by the browser. I had success with the following commands:               http://localhost:8080/geoserver/rest/workspaces               http://localhost:8080/geoserver/rest/workspaces/webGis/datastores.json          http://localhost:8080/geoserver/rest/workspaces/webGis/datastores/TransitFeed          http://localhost:8080/geoserver/rest/workspaces/webGis/datastores/TransitFeed/featuretypes/Shapes.json      The latter returns me a jsom with infomações the Shapes layer, but not only is the data attributes.               {"attribute":[{"name":"geom","minOccurs":0,"maxOccurs":1,"nillable":true,"binding":"com.vividsolutions.jts.geom.MultiLineString"},{"name":"id","minOccurs":0,"maxOccurs":1,"nillable":true,"binding":"java.lang.Integer"},{"name":"shapeid","minOccurs":0,"maxOccurs":1,"nillable":true,"binding":"java.lang.String"},{"name":"sentido","minOccurs":0,"maxOccurs":1,"nillable":true,"binding":"java.lang.Integer"},{"name":"lenght","minOccurs":0,"maxOccurs":1,"nillable":true,"binding":"java.lang.Double"},{"name":"servico","minOccurs":0,"maxOccurs":1,"nillable":true,"binding":"java.lang.String"},{"name":"tipo","minOccurs":0,"maxOccurs":1,"nillable":true,"binding":"java.lang.String"}]}}}      ![enter image description here](http://i.stack.imgur.com/RJBwj.png) How do I return the data in json?? Can I make a query type id = 1?


title: Why does QGIS generate Shapefiles without DBF?

text: I need to output shape files for archiving, which requires generating .shp, .shx and .dbf files. When I save a layer as a shape file it generates the first two, but not the dbf. All the help I can find tells me how to import a dbf file but not export. The dbf files will need to contain ID, X and Y fields. Any ideas?


title: Data visualization techniques for web maps using ESRI and open source suits

text: I want to represent statistical data over the web map and looking the examples/techniques to represent data on the web map (apart from Pie-chat, bar graph) with new tools like gauge widget etc. Therefore, does anyone know other techniques to show information interactively on the web map? Currently I am going though following links & question : 1.Examples of Data Visualization 2.15 fantastic data visualization 3.Good examples of interactive maps and associated technologies 4.What’s new with JavaScript and geospatial - wrapup from the js.geo event Software Platform: ESRI Suit 10.1(ArcGIS Server, Desktop & ArcGIS Java script API) and Open source (Leaflet, Open layer, QGIS Server) Development platform: HTML,Java script, Dojo and Jquery **EDIT :** Looks like D3-based visualizations is suitable for Open source still I am trying to integrate this with ESRI JSAPI too.. Any more tools similar to this will be great !!