title: Convert GML to DXF using ogr2ogr

text: I cannot make ogr2ogr.exe work. I'm using GDAL 1.11. I keep getting the following error message: > ERROR 1: Did not get at least 2 values or invalid number of set of > coordinates 709993.98 6181946.65 32.36 709990.65 6181946.07 32.36 709992.47 > 6181935.8 32.36 709995.78 6181936.38 32.36 709993.98 6181946.65 32.36 ERROR > 1: Invalid exterior ring I run the following command: > C:\Python34\Lib\site-packages\osgeo\ogr2ogr.exe -skipfailures -nlt > MULTILINESTRING25D -dim 3 -f "DXF" -a_srs "EPSG:25832" -spat 711000 6168300 > 713000 6168600 -geomfield posList bygning.DXF Bygning.gml I'm working on a relatively big file 4.2GB. **UPDATE: Example of data to import**               <gml:surfaceProperty>         <gml:Surface srsName="EPSG:25832" srsDimension="3">             <gml:patches>                 <gml:PolygonPatch>                     <gml:exterior>                         <gml:LinearRing>                             <gml:posList>                                 710820.22 6180637.52 36 710818.14 6180633.83 36 710825.47 6180629.71 36 710827.55 6180633.4 38.25 710828.6 6180637.3 38.25 710825.18 6180638.36 38.25 710824.22 6180635.27 36 710820.22 6180637.52 36                             </gml:posList>                         </gml:LinearRing>                     </gml:exterior>                 </gml:PolygonPatch>             </gml:patches>         </gml:Surface>     </gml:surfaceProperty>      Does the GML3 driver in org2ogr support surfaces? **UPDATE 2 - making it simple** I tried making a gml consisting of only one feature, same namespace etc. as before. I then output the feature to WKT. Output WKT "POLYGON ((711383.39 6182528.23,36.49 711379.31,6182528.35 36.48,711379.06 6182520.08,36.48 711383.15,6182519.97 36.48,711383.37 6182527.56,36.49 711383.39,6182528.23 36.49))" Input GML http://pastebin.com/76AGLWH0 Cleary the commas is placed in the wrong positions. The z-coordinate of the first point is parsed as the x-coordinate of the second feature. If I add -dim 3 or -nlt GEOMETRY, it adds zeroes. BR Thomas


title: How to create a time slider using arcgis javascript that only shows a year at a time?

text: I am looking at this example on how to create a time slider using the arcgis javascript api. Instead of showing the data cumulative from the start date. I want to just show the data for that particular year. I have tried changing the ThumbCount to 1, but all that does is set the start time to the start time. Other than that, I can't seem to find any more information on how to get the behavior that I want. I understand that I can use 2 thumbs on the time slider in order to get desired result. However, I am looking to see if it is possible by using only 1 thumb slider.


title: Installation issue with OpenGeo-Suite

text: Ok so last time had a big issue with running WFS. Seems the issue was neither with my OpenLayer code nor Geoserver. The issue was with the servers not able to communicate with each other. The tomcat 6 server on which geoserver was running and apache server which is the default web server for my ubuntu. After some digging around it seems that I need to have tomcat 6 pre-installed on my machine before installing OpenGeo Suite. So having reinstalled the package, at end, the terminal gives me the following message               NOTICE: Please run /usr/share/opengeo-suite/tomcat-admin-setup.sh to complete this      installation.      So i executed the following command (from the required directory),               sh tomcat-admin-setup.sh      And I am getting this error,               tomcat-admin-setup.sh: 143: [: Ubuntu: unexpected operator     Unknown distro. The OpenGeo Suite is supported on Red Hat Enterprise Linux 5 and 6,      CentOS 5 and 6, and Ubuntu Linux 10.04 or later.      Does this mean the OpenGeo Suite is not compatible with Ubuntu 12.04??


title: ESRI Javascript API - Emulate 'Red Glow' Effect From Flex API

text: I am trying to emulate the red glow effect so common in the Flex API Examples. Really, any type of blinking/flashing style marker will do. I am assuming I will have to use a picture marker symbol with an animated gif? I have probably infected my computer with who knows what downloading trial software to convert from swf to gif and all types of other non nonsensical conversions. I am hoping someone can point me to an example or give me some insight on how to accomplish this.


title: How to repeatedly mosaic raster in Erdas?

text: I have found script to mosaic raster in R. But the output raster size is almost 9 times bigger (92mb) than mosaic it in Erdas (12mb). So I prefer mosaic raster in Erdas. But it wastes so much time to mosaic long term (time series) satellite data one by one. So I am looking for possibility to repeatedly mosaic raster to create multiple outputs automatically by single attempt as like batch in Erdas.


title: Opening PDF generated by Python AddIn using Report (*.rlf) file automatcally?

text: As a test I have:   1. Authored a test.mxd with one layer called Localities in a single data frame called Layers   2. Authored a test.rlf using View | Reports > Create Report of ArcGIS 10.2.2 for Desktop   3. Created a Python AddIn which creates the expected PDF file from the Report (*.rlf). This is the code which all works to create the PDF report from the Localities features within the current extent:               import arcpy      import pythonaddins          class ReportForExtent(object):         def __init__(self):             self.enabled = True             self.checked = False              def onClick(self):             mxd = arcpy.mapping.MapDocument("CURRENT")             df = arcpy.mapping.ListDataFrames(mxd,"Layers")[0]             lyr = arcpy.mapping.ListLayers(mxd,"Localities",df)[0]             arcpy.mapping.ExportReport(lyr,r"C:\temp\test.rlf",r"C:\temp\test.pdf","EXTENT",extent=df.extent)             del mxd      Is there a way to have the Python AddIn not just create the PDF file but also open it in Adobe Acrobat Reader?


title: How to install QGIS 2.2 instead of 2.4 to get a working OpenLayers plugin again?

text: I use ubuntu . I don't know how to return installed qgis 2.2 I try uninstall qgis 2.4 and try find how to install old version qgis 2.2 Please .. tell me how to install old version qgis 2.2 for ubunut 12.4


title: Examples of Sanborn map WMS?

text: I'm wondering if anyone has heard of or come across any examples of scanned and georeferenced Sanborn Fire Insurance maps being served as a WMS. I know there are copyright issues, so presumably only for maps before 1922. I'm curious about the possibility of doing this myself... Thanks!


title: Trouble using ogr2ogr w/ -s_srs to combine inconsistent EPSG files into my postgis DB

text: I'm having trouble combining shapefiles from TIGER and different local governments into a single postgis db to query legislative locations. I've imported all data to multipolygon geom, but found some datasets' polygon shapes returning really large floats values lat/lng. Looking more closely, I noticed inconsistent EPSG coordinate systems, based on searching prj2epsg.org with contents of my .prj files:               #Census Data (EPSG:4269 - GCS_North_American_1983)         GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137,298.257222101]],PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]]          #Texas Data (EPSG:3081 - NAD_1983_Lambert_Conformal_Conic)     PROJCS["NAD_1983_Lambert_Conformal_Conic",GEOGCS["GCS_North_American_1983",DATUM["D_North_American_1983",SPHEROID["GRS_1980",6378137.0,298.257222101]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Lambert_Conformal_Conic"],PARAMETER["False_Easting",1000000.0],PARAMETER["False_Northing",1000000.0],PARAMETER["Central_Meridian",-100.0],PARAMETER["Standard_Parallel_1",34.91666666666666],PARAMETER["Standard_Parallel_2",27.41666666666667],PARAMETER["Latitude_Of_Origin",31.16666666666667],UNIT["Meter",1.0]]      I can't seem to get ogr2ogr to explicitly transform the target coordinate systems with `-a_srs EPSG:3081 -s_srs EPSG:3081 -t_srs EPSG:4269`:               ogr2ogr -f PostgreSQL PG:"host=..." file.shp -nlt MULTIPOLYGON -s_srs EPSG:3081 -a_srs EPSG:3081 -s_srs EPSG:3081 -t_srs EPSG:4269 -lco GEOMETRY_NAME=the_geom -nln table_name.      I've tried a number of different variations on the above, but at best I get `ERROR 1: AddGeometryColumn failed for layer table_name, layer creation has failed`. **What am I missing here?** I'm still a bit of a novice on work with shp files. Any help or suggestions would be really appreciated! Example source files:   * From Texas: ftp://ftpgis1.tlc.state.tx.us/DistrictViewer/House/PLANH309.zip   * From Tiger: ftp://ftp2.census.gov/geo/tiger/TIGER2013/CD/tl_2013_us_cd113.zip


title: Connecting Grass with R

text: I need help to connect R with Grass. In RStudio when I send the line: loc <\- initGRASS("C:/Program Files/R/R-3.1.0/bin/i386", home=tempdir(), override=TRUE) It says: Error in system(paste("g.dirseps.exe -g", shQuote(Sys.getenv("WINGISRC"))), : 'g.dirseps.exe' not found Guys, Can you help me?


title: How to translate .bil files in a batch mode with QGIS?

text: I have some .bil files that I would like to convert into ASCII Gridded XYZ format files. It seems to work well with the Raster>Translate tool of QGIS when I process one file at a time. But I would like to do it by batch. So I checked the batch box but got this error message: Processing of the following files ended with error: > > C:/Users/Fall Line Tower/Desktop/Bil To Be Processed/Only > Bill/PRISM_tmax_stable_4kmD1_19810101_bil.bil ERROR 4: `C:/Users/Fall Line > Tower/Desktop/Bil To Be Processed/Onli > Bill/PRISM_tmax_stable_4kmD1_19810101_bil.bil' not recognised as a supported > file format. GDALOpen failed - 4 I have the same error message for each .bil file within my source directory. Any idea where the problem comes from? Does the batch mode take only one type of format? Not the .bil? I am using QGIS 2.2 on Windows 7. Thank you


title: Managing large number of catchments/watersheds

text: Does anyone know of an easy way to manage a large number (10,000's) of catchments/watersheds? I have defined a stream network and for each segment, I have a corresponding catchment (in a unique shapefile). The name of the catchment is the x and y co-ordinate of the catchment outlet which is also assigned to the attribute table of relevant stream segment. The end user will only want to view the stream network and maybe a few of the catchments (certainly not all!) at a time. I need a way to make it easy for the end user to view the catchments of interest (I guess they could search for them). My only idea (which I am a little unsure of how to implement) is to build a tool (in ArcGIS) that loads the catchment when they click on a stream segment (with the tool selected) If anyone has a better idea... For my idea, I think I would have a button similar to the identify tool the end user would use this tool to select a stream segment, the tool would then look in the attribute table to find the file location of the catchment and add this as a shape file. Is such a thing possible to do? I have coded up some python scripts before but have never attempted to integrate such functionality into ArcGIS. Any thoughts or direction would be appreciated. We only have the basic version of ArcGIS, although if this sort of functionality is available in something like ArcHydro I'd be interested to know!


title: What are the advantages/disadvantages of populating a PostGIS database with Raster information?

text: ...as opposed to keeping the rasters as locally stored GEOTIFFs? I have a large project that will require making available approximetely 50,000, 10mb GEOTIFFs as WMS to an OpenLayers based web map. Only 1 WMS will be displayed at any given time so performance is a non-issue. Considering MapServer and GeoServer for this application. I would like to know what are the advantages and disadvantages of storing this type of data. Will updating a PostGIS database when the raster data changes be more difficult, rather than replacing a single (or series of) raster files. WMS load times drastically improved with PostGIS?


title: Looping multiple feature classes in python for ArcGIS 10

text: I currently have a large list of feature datasets and dbf tables in a geodatabase. I would like to join each shapefile with its matching DBF file using the "Make Query Table" tool and then "Copy Features" to produce a new shapefile with specified fields. I am trying to automate this process so I don't have to manually select field names for each set of data (output table would use the same fields). By using the Iterate Datasets in ModelBuilder, I was hoping to do this. However, the Make Query Table will not connect to the geodatabase and still requires me to plug in a specific table. Does anyone know I can use looping to run this process? Are there any similar examples using python? Example of Geodatabase: AAA1_part1.shp AAA1_part1.dbf AAA2_part2.shp AAA2_part2.shp ..... ..... Any information would be greatly appreciated! Thanks


title: Calculating Field based on name of shapefile

text: I have multiple shapefiles to which I have added a field called "gname1". For each shapefile, I want "gname1" to be set to the name of that particular shapefile, for every feature within the shapefile. I realize there are multiple ways to do this, but I want to know how to do this within the "Field Calculator" tool. So, in the "gname1=" box under Field Calculator, what is the expression I use to denote that I want the field to be set to the shapefile name? I am very new to python, and not entirely familiar with the procedure for using python scripts to calculate fields through this tool in Arc10. I have seen other questions similar to this one, but I was confused on how to do this _within_ the field calculator tool. Thanks for any help!


title: Get the different values assigned to a zip code

text: I have a map to which I added a zip code layer and another layer which contains values assigned to every square kilometer within each zip code. I'm trying to find out how many different values there are in each zip code and what they are. So for example, if a zip code has 10 values assigned to it in total, 5 of them being 1.5, 3 of them being 2.1 and the remaining two being 6.4, this zip code basically has **3 different values** , namely **1.5, 2.1 and 6.4**. That is what I want to know. Is there a way to do this without manually checking all the values in every zip code? * * * This is a screenshot of my data. The value that I'm working with is under the column **BI**. I need a table which gives me a list of the long/lat in every zip code and the corresponding values. I can't do anything on excel yet since I don't know which lat/long goes with which zip code. ![data spreadsheet](http://i.stack.imgur.com/wTcJs.png)


title: ArcGIS Javascript API - Refresh

text: The ArcGIS Javascript API has a nice refresh interface for ArcGISDynamicMapServiceLayer:                   var dynamicMapServiceLayer = new esri.layers.ArcGISDynamicMapServiceLayer("http://sampleserver1.arcgisonline.com/ArcGIS/rest/services/Demographics/ESRI_Census_USA/MapServer");         map.addLayer(dynamicMapServiceLayer);         dynamicMapServiceLayer.setDisableClientCaching(true);              if (dynamicMapServiceLayer.disableClientCaching) {             aspect.after(dynamicMapServiceLayer, "onUpdateEnd", function () {                 $.delay(15000, "UpdateDynamicLayer", function () {                     dynamicMapServiceLayer.refresh();                 });             });         }      But this is not generalized across the other layer types. Anyone have a solution for refreshing layers of any type?


title: Is there a simple way to create a new custom vector marker symbol for ArcMap?

text: I would like to create a set of custom marker symbols for environmental sampling locations. I need there to be "vector" marker symbols not picture symbols. I've tried picture marker symbols and there are numerous drawing and exporting issues associated with this approach especially when using transparent colour. I believe I used a font creator in the past and drew a custom true type font but there seem to have been some issues with accessing the custom font from ArcMap as well as some font size problems. Ideally, I would like to draw a vector graphic in Illustrator/ CorelDraw then somehow use this graphic as a marker symbol in ArcMap. Can this be done? I would really appreciate if someone could **describe the process of creating a custom "vector" marker symbol** ; tools used to create the symbol as well as accessing the symbol in ArcMap. Vector graphic approach is preferred but a true type font approach is also fine.   Thanks!


title: leica zeno 5 gps handheld

text: I recently purchased the Zeno 5 gps handheld device form leica geosystems. I would like to ask if anyone knows what other software I can use except the one they offer that integrates with ArcGIS desktop. They suggested the 'Zeno office' and the 'Zeno field' or this software http://www.microsurvey.com/products/fieldgenius/price2007-canada.htm and the 'Zeno connect'. It wasn't clear on the website that you have to buy this extra software so that the zeno 5 gps device will function properly. From my experience I know that when you buy a gps handheld device you dont need to buy extra software to display your position, the satellites, your accuracy, your lat and lon in a user interface. Is there a less expensive alternative software I can purchase and install in this device? Thanks in advance!


title: Print preview showing distortion in layers from GeoServer via MapFish and OpenLayers?

text: I have an application which shows the print preview before printing the map and the problem is that the layer is distorted in the first image but when I zoom into the map the distortion goes away. I am using MapFish as the printing module and I assume that OpenLayers deals with the print preview part. I am new to this so I am not sure if this a code based problem or a setting based problem (either in GeoServer or MapFish) ![enter image description here](http://i.stack.imgur.com/0AIRS.png)


title: Split polygon with holes to get several hole-less polygons

text: Unlike this thread I actually need to split the polygons into several hole- less polygons. I'm not sure on how to approach to this problem. I have polygons with holes of many diferent shapes, like buildings, esplanades, rivers, etc. And I need to split them in a coherent way (if the polygon is a river it should be splitted perpendiculary to his axis) and with less splits as possible. I need to do it programatically. I'm using ArcObjects 9.3, but if there is no other way I could try to use a specific algorithm to solve this problem.


title: Calculate MBR envelope for different geometry shapes in different rows

text: I have some different kinds of geometry (Polygons and line strings) stored in my mysql database. On different rows. I know I can retrieve the envelope of one row with: `SELECT Envelope(zoneShape) FROM geometryTable WHERE refid=4` If there are three different shapes associated with reference id 4, i'm getting three results, obviously. How can I (preferably with one sql query) get the envelope/MBR of all the geometry with e.g. refid=4?


title: How to overwrite a shapefile or feature class using arcobjects?

text: I can copy a shapefile but I want to be able to overwrite the copied shapefile to prevent multiple copies being created each time I run the program. I would also like for the saving and renaming of the copied shapefile to happen behind the scenes i.e. the user doesnt see a gxDialog box. Any suggestions on how to overwrite the copied file, either as a shapefile or feature class would be greatly appreciated.


title: Use output of Add Field as input for Calculate Field in ModelBuilder

text: I have created a simple model in ArcGIS 10 which, after some initial processing of a feature class, outputs `outputfc` and adds a field to it. Regardless of the feature class input (which the user specifies), the field `newField` that is added is always the same. I then need to do some calculations in this newly added field, but because I know this field will only have the name that I specified, I select `outputfc` as the input table for Calculate Field, and an input variable field with the name of the field `newField` as the value. This works fine to achieve the outcome I want, but if I allowed the user to specify the name of the added field, how would I connect this to the next Calculate Field process?


title: Possible ways to display maps offline, inside mobile hybrid app

text: I want to develop a mobile hybrid application to run on both Android & iPhone devices in which i want to display my own map data. I have the SHP files with the data i need. I have Geoserver and PostGIS for loading and serving the map data online. So far everything is fine. My problem is that i still havent figure out the correct walkthrough for displaying the map data OFFLINE. Although i have read many answers here in this website, none of them seems to solve (as per my understanding my issue). The goal is to display points of interrest on map. offline. Some Points of interrest doesnt have address.Just Lon/Lat. I want to cover these requirements:   1. How to offer the map data & the content offline to users device   2. How to offer to user routing options to navigate him to the destination from current location   3. How to add content data from a WordPress CMS, inside the app? (via JSON?GeoJSON?GeoRSS?)  Have been searching and looking many Tile services and many approaches to achieve the requirements, but still havent figure out the correct combination. For the development of app I've been using PhoneGap XDK, Eclipse ADT. And for displaying the maps on mobile devices I'm using OpenLayers mobile. Looking forward your suggestions or walkthroughs. Thank you in advance, with respect


title: QGis not upgrading from 1.9.0. Ubuntu, Raring, using nghtly PPA

text: I'm currently running 1.9.0 on ubuntu raring, through the nightly ppas:               deb     http://qgis.org/debian-nightly raring main     deb-src http://qgis.org/debian-nightly raring main      But it doesn't update with software updater, even though I know there are more recent versions it should be pulling down. Can anyone think off-hand why that might be or what I could do to fix it? I noticed that on the new qgis homepage, the nightly PPAs are not listed for raring... have they been discontinued for some reason, or temporarily broken? Thanks Simon


title: How to set min and max value for stretchType: esriRasterStretch_PercentMinimumMaximum

text: How can I set the min and max values for the stretch type "Percent Clip" in arcobjects? It's possible to set the Stretch type:               rasterStretch.StretchType = esriRasterStretchTypesEnum.esriRasterStretch_PercentMinimumMaximum;      But I have no idea how I can set the min and max values (like in the Properties Dialog of ArcMap). ![Layer Properties with Stretch type Percent Clip and min/max percent values](http://i.stack.imgur.com/T9GLr.png)


title: Is there any "proper" way of setting the viewport for QGis Time Manager?

text: I'm working on video visualisation for some geo data. Everything works fine for me except the viewport that the time manager plugin is using for image frames. It makes image frames size of actual map window size which is something like 1248x636 in my case now. When I resize the QGIS window the time manager plugin starts exporting image frames of different size. Is there any "proper" way of setting the viewport so I don't have to play around with application window size? Or should I perhaps bulk edit (crop?) all the frame images before merging them with mencoder? I'm using Quantum GIS 1.8.0 with Time Manager 0.7


title: retrieve some info from GRIB2 file with JAVA

text: I've just started parsing some info inside GRIB2 file with NOAA weather forecast, and I am realizing this task is not as easy as I could think 'a priori'. I need to process some info inside JAVA routines, and until now I am experiencing with library from Unidata.ucar.edu. The problem is that accessing any specific data inside the file is a little pain in the ass, because of I had never seen a weather interchange file format until now. Therefore, while I read up documentation about GRIB2 file format, if any of you know some JAVA working examples that extract concrete information it would be very appreciated (e.g. retrieve temperature at some level in specific lat/lon)?? I found this library's API just here, but no examples are provided so in the meanwhile it's time to study. I'll update the post with any advance. Best regards. **UPDATE:** I create a Grib2Input object with a RandomAccessFile as param like this:               Grib2Input grib2Input;     String pathfile = "c:/gfs.t06z.pgrb2f03";          //Create RandomAccessFile     RandomAccessFile gribfile = null;     try {       gribfile = new ucar.unidata.io.RandomAccessFile(pathfile2, "rw");     } catch (IOException e) {       e.printStackTrace();     }          //Create grib2Input file     grib2Input = new Grib2Input(gribfile);           try {       grib2Input.scan(false,  false);            } catch (IOException e) {       e.printStackTrace();     }      And here I can get info about GRIB2 file as _Product Status, Grid Name, Grid Units, and so foth_ by getting Grib2GridDefinitionSection of my grib2Input object and then Grib2GDSVariables with getGdsVars() getter on it . Now I suppose **(maybe i'm wrong)** each record is compound ofall sections (1, 2, 3... of GRIB2 spec) First of all: **Is this righg?** For example, record 1 has its GridDefinitionSection, its ProductDefinitionSection, its DataRepresentationSection, and so foth. So, then I try to get one particular record in this way:               List<Grib2Record> records = grib2Input.getRecords();     Grib2Record myRecord = records.get(15); //Get record number 15.       **How can then retrieve information just about this record as Parameter Name, Parameter Units and of course, its data??** The hole file has information about Temperature, Wind, Clouds, and much more, and my problem comes in what I'm seeing that information for each record is the same (no matters wich number of record I get). Of course, I'm doing something wrong or I don't know how the file is format yet. Thankx in advance! **UPDATE 2:** Here are my thoughts: 1) With my _RandomAccessFile_ , I can create a _Grid2Data_ :               Grib2Data gd = new Grib2Data(gribfile);      2) Then I can get data in this way:               float[] data = null;     for (int i = 0; i < records.size(); i++)     {         //First get record information            System.out.println("record " + i +  "-> Param. Categ.: " + records2.get(i).getPDS().getPdsVars().getParameterCategory());         System.out.println("record " + i +  "-> Param. Number: " + records2.get(i).getPDS().getPdsVars().getParameterNumber());         try          {             data = gd.getData(                     records.get(i).getGdsOffset(),                     records.get(i).getPdsOffset(),                      records.get(i).getId().getRefTime() //This param must be new, because it                                                          //is not documented in API                    );         }          catch (Exception e)          {             e.printStackTrace();         }              //Second print some data         for (int j = 0; j < 50000; j+=5000)         {             System.out.println("data[" + j + "] -> " + data[j]);         }              System.out.println("--");     }      **Key question:** If I open GRIB2 file with some software like _tkdegrib_ then I can see information about a **number of records** (pressure, temperature, etc) that **not match with the number of records I can see with the routine above**. For example above, with a GRIB2 from NOAA about weather forecast, _tkdegrib_ shows information about 319 records and I can print information of 315 records ¿? So I can't match my printed data with the extracted with _tkdegrib_ to test this information is the right one. Any suggestion? Am I doing something wrong? Thankx!


title: CartoDB tile proxy

text: I am currently evaluating CartoDB(and really happy with it), and am trying to setup a proof of concept of a secure map. I would like to implement a tile proxy, but really don’t know where to start. I understand we can use this: https://github.com/javisantana/cartodb-tiles- proxy And possibly node.js? I am trying to create a map with an overlay that is accessed through some kind of web portal that users can filter polygons based on numerical sliders/web inputs. It’s important that users can’t see the data in the tables. Any tips or pointers would be greatly appreciated, are there any guides out there?


title: Select area which is at least x meters away from other geometries

text: In my table I have the build-up area from a region as a big number of polygons. > urban_area: id | name | geometry(polygon) I was wondering if and how I would be able to find the area within a bounding box that is at least x-meters away from any geometry in the `urban_area` table. In schematics: ![enter image description here](http://i.stack.imgur.com/0wZRJ.png) With the orange being the geometry in `urban_area`, the grey the offset and the dotted area the desired area.


title: How to zoom to selection on ArcGISDynamicMapServiceLayer?

text: I am using the latest ArcGIS Javascript API. All the examples and code I've found so far only deal with geometries of a `FeatureLayer`. How do I zoom to the visible features of a `ArcGISDynamicMapServiceLayer` or `ArcGISTiledMapServiceLayer` layer? By visible I mean that some of the features of the layers are filtered out thanks to the `layerDefinitions` property. I wish I could keep using FeatureLayers but I can't because I am running in to a wall due the sheer amount of geometry data downloaded, so I have to let the server do the rendering work and use a `ArcGISDynamicMapServiceLayer` to get images instead. Thank you.


title: ArcGIS 10.1 - OD Cost Matrix SOAP Request (Error received)

text: I've been using Network Analyst in ArcGIS desktop (10.1), and using succesfully some SOAP API operations (findAddressCandidates, GetNALAyerNames, GetSolverParameters, GetSolverParameters2) But now want to get OD Cost Matrix (Network Analyst) using SOAP API. I've been following the steps: http://resources.arcgis.com/en/help/soap/10.1/#/Overview/01vp0000001p000000/ Basically, I call getSolverParameters2 to get parameters and content of NAClasses for the OD Cost Matrix Layer. Debugging the code, I can see data returned by getSolverParameters2. When I call solve operation (defined in WSDL) I get and HTTP 400 Error (Bad request) and the following message: "The request sent by the client was syntactically incorrect ()" NAServerODCostMatrixParams params = (NAServerODCostMatrixParams) port.getSolverParameters2("OD Cost Matrix", "", true); NAServerSolverResults results = port.solve(params); (Also, I've been trying fill NAServerODCostMatrix 'by hand', unsuccesfuly ) I trace the SOAP request message and I can see the content. The same results I get when I follow this steps using SOAP UI. Do you have any idea about this problem and/or to solve it? or what can I try? Regards, Hernan


title: Find and list all polygons that overlap with another polygon

text: I am working with two polygon layers which overlap but are unrelated (one is congressional district boundaries, the other is MPO area boundaries). I would like to generate a list of all the polygons in one layer which overlap with each polygon in the second layer. I.e. I want to generate a list of all the Congressional districts which are at least partially within each MPO area. Which vector tools or plugins can I use to get there? I'm using QGIS 1.8 on OSX 10.6.


title: ArcPy - Dividing every raster cell value by the maximum cell value

text: I want to divide every cell of a raster by the raster's maximum cell value using ArcPy. Here is what I have done so far.               densityRaster = LineDensity(shapefileName, None, "", "", "")     maxVal = arcpy.GetRasterProperties_management(densityRaster ,"MAXIMUM")     normalRaster = densityRaster/maxVal #maxValue i get is 5.66      This gives me an error described below. Anybody has clues to where I get it wrong? > Runtime error : ERROR 000732: Input Raster: Dataset 5.666 does not exist or > is not supported If I enter `normalRaster = densityRaster/5.666`, I get the correct result.


title: Openlayers popup - no scroll bars on ios or android

text:
I have an openlayers app with a FramedCloud popup. The popup has a lot of lines of data, so scrolling is needed to get to the bottom. All works well on my pc with Firefox. However, out in the field with an iPhone or Android (running Chrome), I get no side scroll bar, so cannot get down to the end of the form. Here's some popup code that gives the problem:               function onFeatureSelect(event) {              var feature = event.feature;              var content = "" + feature.attributes.art.value +               " (Id: " + feature.attributes.trädnummer.value +               ")
Risk: " + feature.attributes.risk.value +              "
" + feature.attributes.anmärkningar.value +               "
Åtgärd: " +               feature.attributes.rekommendationer.value;         content += "

A string";         content += "

A string";         content += "

A string";         content += "

A string";         content += "

A string";         content += "

A string";         content += "

A string";         content += "

A string";         content += "

A string";         content +=  '
<input type="submit" value="Spara"">';              popup = new OpenLayers.Popup.FramedCloud(             "chicken",             feature.geometry.getBounds().getCenterLonLat(),             null,             content,             null,              true,              onPopupClose);              feature.popup = popup;              map.addPopup(popup);     }      Does anyone have a solution to this? Thanks Mini


title: Arc10.1: intermittent Error using arcpy.Delete() on temporary rasters

text: I've been pulling my hair out with an intermittent error when running from either the MS-DOS prompt, PyScripter, or WingIDE. Haven't bothered to check w/in Desktop apps (but a colleague says it's the same). I run using w/in a 10.1 sp1 Win 7 Desktop application. The code I've mocked up to demonstrate is as follows:               import arcpy     arcpy.env.overwriteOutput=True     arcpy.CheckOutExtension('Spatial')     arcpy.env.workspace = r'd:\wk_geospatialfabric\databin06\lcov_comp'     arcpy.env.scratchWorkspace = arcpy.env.workspace          for i in range(10):         try:             print 'Iteration: ' + str(i)             r = arcpy.sa.Lookup('lulc_mini', 'COVTYPE')             c = arcpy.sa.CellStatistics([r, 1], 'MINIMUM')             arcpy.Delete_management(r)             arcpy.Delete_management(c)         except:             print 'Iteration: ' + str(i) + ': FAILURE'      It fails on the first `Delete()` statement, usually every third iteration. If I do this outside the `try/except`, I see:               Traceback (most recent call last):       File "<string>", line 254, in run_nodebug       File "L:\software\MoWS_NHDPlus_Python\error_cellstatsDelete.py", line 11, in <module>         arcpy.Delete_management(r)       File "C:\ArcGIS\Desktop10.1\arcpy\arcpy\management.py", line 3658, in Delete         raise e     arcgisscripting.ExecuteError: Failed to execute. Parameters are not valid.     ERROR 000732: Input Data Element: Dataset d:\wk_geospatialfabric\databin06\lcov_comp\Lookup_lulc_1 does not exist or is not supported      I've messed with inserting an `arcpy.Exists()` in there prior to `Delete()`. On the problem iterations, it takes an oddly long time to return `True`, and then fails anyway. I'm wondering if there is some latency/asynchronicity that I'm getting hung up on. The IDEs are using EPD's Python 2.7.2 (default, Jun 12 2011, 15:08:59) [MSC v.1500 32 bit (Intel)] on win32. I've been trying to use something like sys.executable to determine what get used w/in Desktop, but can't quite figure that out. Any ideas? **UPDATE** Tried adding `arcpy.RefreshCatalog(arcpy.env.workspace)` both before and after the `Delete` statements, but no impact.


title: Transferring ArcGIS Server directories to another location (D drive) on AMI EC 2

text: Updated question: I am using ArcGIS Server 10.1 + ArcGIS desktop 10.1 on AMI EC2 instance & I am facing space problem so trying to move existing server directories to different drive. Currently our dept is under maintenance so I am copying “arcgisserver” folder (Currently it is on C:\arcgisserver) to my “D drive”. So my question is that which things should I take care (like data ,rest service etc) ? I am familiar with ArcGIS Server 10 and it was easy to perform above settings in 10.0 as compare to 10.1 Currently I am going through ESRI web help and researching on it Any help will be great !!! Thanks in advance : )


title: How to determine average distance from points to line?

text: I have a polyline road dataset and a study area (polygon) covering the the road dataset. I want to determine the average distance for every given point within the study area to the nearest road. i don't have a dataset with the points. Also I don't necessary want to create a new dataset. All I am interested is the average distance. Is there a tool that can do that job (preferable ArcGIS)?


title: PostGIS Union Between Two Tables, Preserving Attributes

text: I'm trying to union US State level upper- and lower-level government boundaries in such a way that generates one PostGIS table with both the upper and lower level names for each shape. In my example below, `U` is a dataset and has a column with rows `UA` and `UB`. `L` is another dataset and has `LA`, `LB`, and `LC`. When I union the shapes together I should get the same shapes as `L` but the new, unioned dataset has both data columns. ![enter image description here](http://i.stack.imgur.com/z8Kog.png) **Edit** : My example image above is pretty crummy. The answers so far have been related to querying the data assuming the geometry is available. This is important but is skipping the more confusing part of my question. A better example: ![image](http://i.imgur.com/ejJHGip.png) The yellow shape is from `U`, the orange shape is from `L`. Notice that they overlap but some pieces of the polygons aren't contained completely inside another shape (like in my example where `LB` and `LC` are completely within `UB`). The question I was trying to ask is about correctly chopping the two layers so the resulting data set has no overlaps. See for example this mailing list post that gets close but doesn't quite work.


title: why sql file generated by postgis was different in notepad and pgAdminIII?

text: When i imported a shapefile which contains some polygons into PostGIS by command line on Win 7, a problem puzzled me: The sql file generated by shp2pgsql opened with notePad has the geom (Mutipolygon) value ,however when put this sql file in pgAdminIII, the value was empty..why??some one can help me??


title: Global Mapper - .SP1 conversion to Shapefile

text: Does anyone know how to convert an .SP1 file to an shapefile from Global Mapper to ArcView 10.1? The .SP1 is a GPX file.


title: Select polylines that are having mutual points from another table with condition

text: I am struggling with this MapInfo SQL (or T-SQL) query too long so I need little help :) I have two tables:   1. PIPES with polylines   2. CHECK_NODES with points that represents starting and ending points from PIPES polylines What I want is to select PIPES that are having mutual points from CHECK_NODES where CNT value is 2 and field values of MATERIAL and DIAM are the same for those 2 pipes wich are joined via that node. ![enter image description here](http://i.stack.imgur.com/UbXcF.jpg) *Note: If CNT=2 that means that only 2 pipes are connected to that node. Its not a problem for me to do first part of the query but I am stuck with a part where I need to further select MATERIAL and DIAM ... So far i could do only               Fields {*}     Tables {CHECK_NODES, PIPES}     Where {CHECK_NODES.obj Intersects PIPES.obj and CHECK_NODES.CNT = 2}     Group {}     Order {}     Into {pipes_for_join}     Browse      p.s Data resides on MS SQL Server so if someone could provide T-SQL its also ok for me. Structure for table: CHECK_NODES               1  NODE_WKT      2  CNT           3  MI_PRINX (NODE ID)      Structure for table: PIPES               1  MATERIAL     2  DIAM         3  SOURCE       4  FROM_NODE     5  TO_NODE       6  MI_PRINX (PIPE ID)      Thanks.


title: loss of attribute data when copying between layers

text: ![enter image description here](http://i.stack.imgur.com/xjDDR.jpg)   I am trying to copy points from one layer to another in QGIS 2.2 The problem I get is that all the attribute data from the copied points is missing when copied onto the new layer, the columns show up as blank. Both layers have the same columns. Is this a problem with 2.2 or is there an easy fix for this?


title: How to set up WebMap

text: I'd like to set up a webmap that clients can drop into their existing website. I'm thinking Geoserver, OpenLayers & GeoExt. It would be very simple slippy map kind of thing. Do I need Geoserver at all then? Just want to be able to use WMS/pan/zoom/turn layers on and off. If I ditch GeoServer, then do I need Apache? And where would be the best place to host this? Could most hosting companies handle it if there was no Geoserver since OL & GeoExt are just JavaScript? I don't want to host it myself. Anyone who could point me in the right direction for a checklist or roadmap to offer up simple webmaps is greatly appreciated. I'm looking for info on setting up the web side of it, since I'm fairly good with GIS already.


title: PyQgis: How can I "append" features to a Vector Layer instead of "stacking" them?

text: As I read data from a text file and add points to a QgsVectorLayer (as well as related attributes), each new record gets 'stacked' on top of the previous one in the Attribute Table (and within the Vector Layer as well) and I end up with the first record read having the highest index and the last record having index 0 and this is not what I want. Short of a solution, I am currently inserting an 'AutoID' field that is incrementing as I read my text file but I find this rather clumsy :               geom = QgsGeometry()     feat = QgsFeature()     layer = QgsVectorLayer("Point?field=AutoID:integer&                            field=Line:string&                            field=Point:string&                            field=Easting:double&                            field=Northing:double&                            field=Elevation:double", file_name, "memory")     layer.startEditing()     i = 0     with open(data_file, 'rU') as filein:         ##read the text file, get data and attributes in corresponding variables         feat.setGeometry(geom.fromPoint(QgsPoint(x, y)))         feat.setAttributes([i, line, point, x, y, z])         layer.addFeature(feat,)         i += 1      This is the end result: the last record read (AutoID 5326) gets index 0 in the table, so compared to my data file, the Attribute Table is just 'reversed' which is no big deal, but annoying. I need the index of my Vector Layer to be the same than that of my data file.   ![Attribute table example](http://i.stack.imgur.com/2MFGK.png)


title: Using ArcPy and Python to export Table View as PDF format?

text: I am using ArcPy to select from a layer and display the result on map as a table view. Now I want to export the table view as a PDF format. Anybody know about it? If it is beyond ArcPy to export as PDF directly, with intermediate format also welcome.


title: How to save separate output files when iterating a model?

text: I need to run a model in model builder 1000 times in order to obtain the average distance of random points in a wildlife species' home range to anthropogenic features. I need to find the average of the random points from EACH of the iterations. However, when running the model 1000 times, the output file is continuously being replaced, so you only end up with one output table, the very last output table. Does anyone know how I can use model tools in order to force the model to save separate files for EACH iteration so that I can find the average distance of the random points to the features for EACH iteration? Thank you very much, any help would be hugely appreciated! I cannot seem to find ANYTHING online about this problem!


title: How to update location via a data streame in openlayers or Leafletjs

text: I'm developing a application where I need to real-time update a Location of a Map using a data stream, I'm planing to use web-sockets, Are there any guides I can follow? Or can you guys direct me in the correct path? Thanks in Advance


title: Slope regions from elevation points using QGIS 2.0

text: I am using QGIS to create a layer of slope region polygons (e.g. less than 30% slope, 30-50%, more than 50% categories). My data comes from point elevations and I interpolated using SAGA b-spline approximation. I then used the Raster Terrain plugin to calculate slope. However, after using the raster calculator and polygonizing the categories of interest, I still have a problem in terms of how neat the borders are. I read about region growing, but am not sure if I have enough information to do so/how to do it. This above image resembles what I'm aiming for. The bottom image is what I currently have. ![enter image description here](http://i.stack.imgur.com/hC0DP.jpg)


title: What is the difference between FID and OID

text: What is the main difference between OID and FID in ESRI software?


title: What software do you use as complementary to your GIS desktop/server?

text: What software do you use as complementary to your GIS desktop/server in your work as GIS programmers/analysts and what programs are an absolute must you can't do without?


title: Maturity of Mapserver mapscript C# bindings?

text: I've been thinking about possibly porting an internal custom GIS application built on top of Python / PyQt / Mapscript / GDAL-OGR / PostGis to C#. I was wondering about the maturity of the C# mapscript bindings (and to a lesser extent, the GDAL/OGR bindings). Can they be considered as mature as the current (very nice) Python bindings since I assume both are SWIG generated? (The most time intensive part will probably be redoing the current PyQt GUI parts) Thanks in advance for any insight / experiences concerning this topic!


title: How do you determine which dbclient connection property to use for an Oracle SQLWorkspace connection?

text: I have an Oracle11g database that I connect to via a SQL Workspace; however the actual dbclient string value that **works** differs depending on machine (and I'm guessing oracle client configuration). On some machines, its `Oracle11g`(like the example on resources.arcgis.com):               // Create the connection properties.     IPropertySet connectionProps = new PropertySetClass();     connectionProps.SetProperty("dbclient", "Oracle11g");     connectionProps.SetProperty("serverinstance", "Riverside");     connectionProps.SetProperty("authentication_mode", "DBMS");     connectionProps.SetProperty("user", "BlackDog");     connectionProps.SetProperty("password", "123456");      Sometimes it's `Oracle`.               // Create the connection properties.     IPropertySet connectionProps = new PropertySetClass();     connectionProps.SetProperty("dbclient", "Oracle");     connectionProps.SetProperty("serverinstance", "Riverside");     connectionProps.SetProperty("authentication_mode", "DBMS");     connectionProps.SetProperty("user", "BlackDog");     connectionProps.SetProperty("password", "123456");      How is that dbclient value used by ArcGIS? Is there a way to detect which is needed at runtime?


title: Is there any library for converting projection transformation in php?

text: I would like to transform data from lat/long ETRS89 to y/x UTM ETRS89 through php. Any idea?


title: How to use ArcGIS Data Interoperability tool in ArcPy?

text: I created an ArcGIS Interoperability Tool, which I like to use in my ArcPY script. In order to find the correct code I created an ArcGIS Toolbox Model including my Interoperability Tool, and exported as an ArcPy script. Everything looks fine, I checked google and stuff. I get an DataInterop license, I import the Toolbox with my Interoperability Tool (which uses an Alias). Then I define the input and ouput. When I try to start the Interopability Tool I get the following message : AttributeError: 'module' object has no attribute I didn't find anything usefull on the internet which can be found with the above error message and Data Interoperability. Does somebody know how to use a self create Data Interoperability Tool with ArcPy?


title: Choropleth map turns one color when filtered, sql error

text: My data is generating a choropleth map perfectly but when I try to filter (both in the filter tab and by writing my own SQL) it turns all one color and tells me I have an SQL error. This happens even when the SQL is generated by CardoDB's filter widget. I have a big file - about 56,000 rows. Is this a file size issue or something else?


title: Openlayers active select control on layer prevents drag and zoom

text: I've build this web-map http://gimoya.bplaced.net/Openlayers/Trails/map.html with bike-trails (retrieved from a KML on my webserver), where several trails from a region are covered by a selectable region's polygon. Selecting a region's polygon will trigger a zoom to the polygons extent. Now, when zoomed inside the selected polygon the active select control prevents dragging and zooming. Where in my script would I have to deactivate/activate the select control to enable dragging and zooming ontop of a selected polygon? Thanks for any pointers!


title: Is there any way to find the shortest distance from a point to a line?

text: I´m working with QGIS 1.8.0. I have a map divided into many features (thousands of small rectangles) and I have a point at the center of each of those features. I also have a map of the transmission lines of the country. Is there any method to measure (not one by one) the shortest distance from the points to the transmission lines? Does the maps have to be projected into any specific CRS? Regards ![enter image description here](http://i.stack.imgur.com/xEqS4.png)


title: Vertex point is returning None

text: I have this multilinestring that has 48 NPoints. When I execute `ST_PointN(line,5)` it returns none. And to continue on, it returns None for everything after index:2 up to index:48. This is how my line looks:               MULTILINESTRING((-3.1642070167833 55.9269296196706,-3.16421351659546 55.9268662214904),(-3.16421351659546 55.9268662214904,-3.16421636372824 55.9268384509897),(-3.16421636372824 55.9268384509897,-3.16422182573761 55.9267851753802),(-3.16422182573761 55.9267851753802,-3.16422870102352 55.926718114886),(-3.16422870102352 55.926718114886,-3.16423309121073 55.926675293667),(-3.16423309121073 55.926675293667,-3.16423565148822 55.9266503211093),(-3.16423565148822 55.9266503211093,-3.16424103159897 55.9265978443265),(-3.16424103159897 55.9265978443265,-3.16424680776317 55.9265415044985),(-3.16424680776317 55.9265415044985,-3.16425267254583 55.9264843002995),(-3.16425267254583 55.9264843002995,-3.16425541048045 55.9264575949012),(-3.16425541048045 55.9264575949012,-3.16426111146586 55.9264019883556),(-3.16426111146586 55.9264019883556,-3.1642667032531 55.9263474469124),(-3.1642667032531 55.9263474469124,-3.16426957768543 55.9263194101362),(-3.16426957768543 55.9263194101362,-3.16427488261739 55.9262676666359),(-3.16427488261739 55.9262676666359,-3.16428009893088 55.9262167875066),(-3.16428009893088 55.9262167875066,-3.164282741107 55.9261910161221),(-3.1642875546472 55.9261440655823,-3.164282741107 55.9261910161221),(-3.1642875546472 55.9261440655823,-3.16429466890915 55.9260746741522),(-3.16429466890915 55.9260746741522,-3.16430092974527 55.9260136069079),(-3.16430092974527 55.9260136069079,-3.16430822838418 55.9259424170929),(-3.16430822838418 55.9259424170929,-3.16431547242401 55.925871759829),(-3.16431547242401 55.925871759829,-3.16431448732505 55.9258328901507),(-3.16431770120536 55.9257327846001,-3.16431547242401 55.925871759829),(-3.164339 55.925777,-3.16431770120536 55.9257327846001))      any ideas anyone what might be going wrong? Thanks in advance!


title: How to download new version of QGIS 2.2?

text: Can anybody tell me how to Download Qgis 2.2 new release version? In QGIS website they mentioned in download window QGIS 2.0 . if u checked released data it was 24 feb 2014. So where i wll get QGIS 2.2 newly released version? Thanks and Regards


title: How to add a table into the TOC in a Python script?

text: I recently experienced some troubles in adding a table from a specific folder to the TOC using a Python script. Any idea on how can I solve? Thanks, Annalisa


title: Halo use background color

text: As you see on the attached image, I labeled contour lines on the map (in meter). Which function will allow the "halo" to take the background color of the map? On the image, the halo is greenish which diminish the look of the map. I use ArcMAP 10. Thanks in advance ![enter image description here](http://i.stack.imgur.com/sI84a.jpg)


title: Postgresql on Ubuntu Server and QGIS

text: Complete newbie... I taught myself the basics of GIS on the job with ArcGIS 9.3. I have been plating with QGIS and have gotten the basics down on both. I am not well versed in the terminology related to all this, hence why I think I am confused. As of now I simply import shapefiles and raster data manipulating them for each project. In reality I am only doing this to create maps for reports, etc and there is not a lot of analysis that needs to be done. Maybe once in awhile. FYI... I am a Landscape Architect, most may understand the type of maps I create. If I understand correctly I can use **Postgresql** on my Ubuntu Server and then connect to it from **QGIS** on a client. In doing this I could set the **data** up so that once I connect to it the **symbology** is saved amongst other things? I hope that is a correct statement. Is someone able to set me straight here? Maybe provide some tutorials, maybe some on proper installation and connections with Ubuntu Server 12.04 LTS and Window 7 and Ubuntu 12.04 Desktop clients. Thanks in advance.


title: How can I suspend default shortcut menu and substitute with a custom menu in the TOC?

text: I looking for solution **(using ArcObject SDK 10.1 (VB.NET) ArcMap Addin template)** to substitute the default shortcut "right-click" menu in the TOC (When a TOC Item is right-clicked) with a custom menu when a custom tool is active. (I would like default functionality when a custom tool is inactive) I already do that in the map - I have a custom tool and when the tool I show a custom menu on right-click. The custom behavior only works within the map frame. How can I do the same in the TOC? ![enter image description here](http://i.stack.imgur.com/c0N4f.jpg)


title: How to avoid duplicate GIDs when copying features from one PostGIS layer to another?

text: When I enter new elements (features) in another Postgres layer, I can do it in two ways:   1. Drawing new element(s) (with 'Add feature') which I do rarely or   2. Copying (or cutting) some elements from another Postgre layer (Source layer), and paste it into target layer which I do frequently In first example, saving of edits works normal because this layer gets **gid** from postgre database sequence *nextval('layer_name_gid_seq'::regclass)* In second example, I got an error during saving edits, because while copying element from source layer to target layer qgis copied **gid** of element from source layer. When attempting to save edits this error is returned: > Could not commit changes to layer „Cjevovodi“   >  Errors: ERROR: 1 feature(s) not added.   >  Provider errors:   >  PostGIS error while adding features: ERROR: duplicate key value violates > unique constraint "cjevovodi_okill_pkey"   >  DETAIL: Key (gid)=(5) already exists.   > I tried to copy *nextval('layer_name_gid_seq'::regclass)* in field gid, but this sequence can`t be pasted in field gid as field is defined as numeric. Does anyone know simple way to copy elements from source layer (with existing gid) assign new gid? Thanks!


title: How to assign correct flow direction to stream network polylines?

text: I would like to create a map showing the results of watershed analysis as well as the resulting stream network with flow direction. I have a stream network which has been created from the flow accumulation raster. I also have a stream link raster which was created from the resulting stream network raster and the flow direction raster. Once I convert the stream network or the stream links to vector polylines the resulting feature class does not retain the correct directions of the various segments. What is the approach? I have ArcGIS 10.2.1 Standard, Spatial Analyst, 3D Analyst _Red lines on the below example show the network I would like to show the map along with arrows indicating direction of flow:_   ![enter image description here](http://i.stack.imgur.com/Wpitu.jpg)


title: draw extend of print composer map as rectangle in canvas

text: After defining the extend of my printed map in the printcomposer I can show the map center in the canvas by hitting "view extend in map canvas". So far so good. But is it also possible to show the exact **extend** with all borders? Or even better: map composer extend is shown as rectangle in the canvas. At the moment I get that result by copying the coordinates of the map-composer extend to a txt-file (;-sparated), import it to the canvas and connect the points as polygon. But if I change the extend of the print composers map, I have to do all steps again. So is there a other way? ![enter image description here](http://i.stack.imgur.com/UcsCA.png)


title: Using pgr_createTopology on only a subgroup of the table

text: Lets say I have a table with the standard id, source, target, and cost needed for pgr_createTopology. Now lets also say I have a 5th column called floor. Is it possible to run pgr_createTopology on just floor = 3 entries. Basically I want to keep the different floors separated. Atm, when i do the standard normal call, it mixes and matches all the floors together. Edit 1: My solution to the problem was to export the the specific floor to a temp table. Run pgr_createTopology on the temp table, and import that date from the temp table to my table. However, this is really ugly. Is there a way to do it without a temp table


title: Lidar and infrared data to detect petroleum releases

text: Would this combo detect differences in soil temperature in proximity to the release due to the higher temperature of pipeline fluids or would the surface temperatures be affected by the volatilization of the released fluids? I know this technique is being used currently and I am just curious about the process by which it would be able to determine a release point.


title: Does Google Maps follow W3C guidelines for accessibility

text: Hello   Does google maps conforms with the w3c accessibility recommendations(A,AA,AAA)? Cheers


title: OpenLayers: show polygon of images in mobile

text: I'm using Phonegap 2.6 and OpenLayers 2.12 for my mobile application. **Feature:** in a panel (jQuery Mobile), I'm clicking in a checkbox and georeferenced images are displayed in the map that is already loaded. These images fits with that we call "blocks" in the map. **With Firefox:** usually, I'm working with Firefox at the beggining of developing stage. So the feature is working, the OpenLayers image is well created and displayed on the map. ![Working with Firefox](http://i.stack.imgur.com/Qg4pU.png) **With the device:** that doesn't work. Instead it draws a polygon around each image. ![Not working on device](http://i.stack.imgur.com/DqX0L.png) The method that provokes that is addLayers (here is the documentation: http://dev.openlayers.org/apidocs/files/OpenLayers/Map- js.html#OpenLayers.Map.addLayers ....) Do you have any idea of what happens? The only difference between the working and not working case, is that with the device my code is encapsulate with phonegap, but I have no clue why that should be a problem. Thanks


title: How to create an accurate 5 km buffer in ArcGIS 10.1?

text: I need 5 km buffer of a basin boundary of Brazil using ArcGIS 10.1. The basin boundary is in GCS. I have tried it giving buffer distance 5 km and 0.0454545 degree decimal (assuming 1degree=110km). I have got two different results for two distances. Which one is more accurate?


title: opening dbf without losing ObjectID and geometry?

text: I found a lot of related questions, but not completely an answer to my struggle. I can open the dbf-file in ACL and save it as dbf. I added some attributes and wanted to open it again in ArcGIS, but apparently I loose 2 columns when opening my dbf in ACL (and it does too when I open it in Excel): the objectid and the geometry. How can I prevent this from happening? Or are there other options I have to join the object ID's with a new attribute that I only have in a table format (no GIS-coordinates). Cheers, Mieke


title: How can I convert WKT to SHP?

text: Using ArcGIS, ArcInfo and fme, I would like to transform data which contains wkt geometries (points, lines and polygons) into SHP. Has anyone a workflow for doing such a transformation?


title: QGIS Atlas generation - crash dumped

text: In QGIS 1.9.0 there's a new funcionality named Atlas generation. Everything works fine until I start generating more 5000 images. Now, after some time (not always the same) the QGIS crashed with the following error: Crash Dumped minidump writtn to C: \users\smart\appdata\local\qgis-20...ac97365etc.dmp Anyone knows whats happening?


title: Raster smoothing

text: I am an italian student. I'm using ArcMap 10. For my thesis I need to simulate a debris avalanche, but in my elevation raster the avalanche deposit is included so I'd need to take it off to simulate on a pre-avalanche-like surface. Unfortunately the only data I have is the average (and supposed uniform) thickness of the deposit of 24m. I removed this thickness from my raster using raster calculator ("my raster" - "raster of tickness created from polygon to raster"). What I obtained is a raster without that thickness, but like a valley with sub-vertical walls! I need to smooth these edges to make it look more realistic (now it looks like a print with the shape of the deposit). Which tool can use? What I was able to find is the Filter tool which smooth these edges too mildly or the Smooth Line tool which could work on the contour lines. Is there a tool to smooth directly the raster? ![Image of raster dem with 24 m thickness removed in each point of the deposit](http://i.stack.imgur.com/hMxWe.jpg) Thank you, Michele


title: What map source/API/tools are recommended to create Web GIS/Mapping application utilising Microsoft technologies?

text: I'm researching about technologies (around Microsoft tech) to create mapping application and would like to have advice/suggestions. Brief requirement to build mapping application with Microsoft Platform:   1. Adding Pins on the map.   2. Clicking on the pins to get more information as popup.   3. There may be links connecting to various pins (i.e. shortest connection, and not actual physical representation of path on the map). In some near future there may be requirement to represent actual physical path between two pins.   4. Clicking on the link to get more information as popup.   5. Base maps can be hosted externally or internally. Which **base map source can a good candidate** for any apparent reasons? (i.e. Google, Bing or OpenStreetMap). **Suggestions for API/Controls** for creating web application using **ASP.NET** or **Silverlight** app or windows based **WPF** application. How **GeoServer or MapServer can be useful** in for such applications? **Suggestions for other notable tools/API/controls** for designing such applications


title: How to determine irrigable acres from a single source/point?

text: I am trying to graphically depict the exact amount of irrigable acres from a single point, in this case a headgate from a drainage ditch. Basically, I have the original water right decrees (i.e. confusing legal documents) from the late 1800's and am interested in accurately illustrating the legally allowable area to which they can irrigate. Example: From the stream a drainage ditch runs west for 1 mile (on contour) before it opens into a small pasture, irrigating 50 acres below its outlet. Is there a way to show this acreage below this specific outlet point? Perhaps a sort of reverse Snap Pour Point operation? Some function that fills in a specified area below a certain point, following topography lines? :![Crude conceptual example](http://i.stack.imgur.com/P36xp.jpg)


title: How do I export attributes with Python scripting in ArcGIS 10?

text: How do I export an attribute table to dbf or txt format by use of Python in ArcGIS 10?


title: GetParameterAsText() error in python script

text: i have a huge problem with an error . when i use getparameterastext and use it for searchcursor , i encounter to an error .for example the below code               import arcpy          featureClass = arcpy.GetParameterAsText(0)          populationField = arcpy.GetParameterAsText(1)          rows = arcpy.SearchCursor(featureClass)      the error is : > Runtime error : ERROR 999999: Error executing function. the error is different in pythonwin . the error is : > Traceback (most recent call last): File "", line 1, in File "C:\Program > Files\ArcGIS\Desktop10.0\arcpy\arcpy__init__.py", line 804, in SearchCursor > return gp.searchCursor(*args) File "C:\Program > Files\ArcGIS\Desktop10.0\arcpy\arcpy\geoprocessing_base.py", line 357, in > searchCursor self._gp.SearchCursor(*gp_fixargs(args))) RuntimeError: ERROR > 999999: Error executing function. please help me . thank you


title: zoomToExtent (getDataExtend) combined CQL filter, does not work

text: I'm using openlayers 2.12 and geoserver 2.1.3. I created a search functionality for a vector layer that contains points/lines/polygons. I used CQL that dynamically gets data from a form. When I use `map.zoomToExtent(pins.getDataExtent());` does not work at all, instead it zooms 100% on an irrelevant place of the map. Never zooms to contains just the newly created features. Map also contains other vector layers with no search functionality and a Google maps layer. All projections are `EPSG 900913`. Check the following code. Form's "name" is a textfield, "era" is a dropdown and "borough" is a dropdown containing string that will turn to `Bounds` , so map will zoom to a borough. Search for : era OR name OR era+name with `zoomToExtent (getDataExtend)` on = zooms 100% at irrelevant place. I want this to zoom at new features. Search for : era OR name OR era+name with `zoomToExtent (getDataExtend)` off = does not zoom, as it should be Search for : era+borough OR name+borough OR era+name+borough with `zoomToExtent (getDataExtend)` on = zooms 100% at irrelevant place. I want this to zoom at new features. Search for : era+borough OR name+borough OR era+name+borough with `zoomToExtent (getDataExtend)` off = zoom,at the borough as it should be. The code...               //first, clear the layer     pins.removeFeatures(pins.features);     pins.refresh({force:true});          //split the chosen borough, to create bounds, so the map can zoom there     var binfo=document.getElementById("boroughselect").value;     if(binfo!=" "){     var split_binfo = binfo.split(",");     var bidchoosen=split_binfo[0];     var zoomTo=new OpenLayers.Bounds(split_binfo[1], split_binfo[2], split_binfo[3], split_binfo[4]).transform(new OpenLayers.Projection("EPSG:4326"), new OpenLayers.Projection("EPSG:900913"));     map.zoomToExtent(zoomTo);}          //create a filter using CQL     var filter_list = []          if (document.getElementById("eraselect").value != " ") {         var filter_b = new OpenLayers.Filter.Comparison({             type: OpenLayers.Filter.Comparison.EQUAL_TO,             property: 'p_era',              value: document.getElementById("eraselect").value         });         filter_list.push(filter_b);}          if (document.getElementById("nm").value != "") {         var filter_c = new OpenLayers.Filter.Comparison({             type: OpenLayers.Filter.Comparison.LIKE,             property: 'p_name',              value: document.getElementById("nm").value         });         filter_list.push(filter_c);}          //create the final filter,combining the above,  to search     var filter_final = new OpenLayers.Filter.Logical({         type: OpenLayers.Filter.Logical.AND,         filters: filter_list});          //pass the filter to the layer     var prot =  new OpenLayers.Protocol.WFS({         url:  "http://localhost/geoserver/wfs",         featureType: "pins",         featureNS: "http://www.cultmap.gr",         defaultFilter: filter_final});          var _CallBack = function(resp) {         pins.addFeatures(resp.features)          var cb = pins.features.length;     if (cb == 0){         document.getElementById("formMsgG").innerHTML = No results, please try again.";                 }};          var response = prot.read({callback: _CallBack});          pins.refresh({force:true});     map.zoomToExtent(pins.getDataExtent());      Thank you for your time slevin


title: Data Transfer of a point on line using the anchor point

text: Can I a line feature, the anchor point touches a point feature transferred the field value, even when the point touches another line(s)? How does this work? ![enter image description here](http://i.stack.imgur.com/Tmzmf.png)


title: How do I describe a special relationship between connected edges?

text: Consider this simple situation where three edges connect at a node: ![edge relationships](http://i.stack.imgur.com/6gGIQ.jpg) I’d like to write a succinct and clear description of the relationship between A and B in such a way that differentiates it from the relationship between A and C. Something like “when traversing the node in the clockwise direction, A is _adjacent?_ to B, but A is not _adjacent?_ to C.” But it’s not really adjacency. Said in a different way: imagine you’re standing on the node and you’re facing towards A. You start spinning yourself clockwise. The next edge you’ll come to is B, not C. Is there a way to describe this relationship between A and B in a more succinct, formal, or correct way than I’ve written above? It must be directional (one relationship of this type exists in the clockwise direction from A, and another exists in the counter-clockwise direction). And it must scale up to cases where more than three edges are connected at the node. Maybe it has something to do with routing? (I am thinking about this in the context of road networks.) Two approaches I’ve tried already but haven't gotten far with:   1. _9IM-like topology references_ : I’ve looked at the DE-9IM, and even though I'm not a mathematician, I think I can still tell from the diagrams and terms that it doesn’t cover this type of relationship. Neither do I find it yet in the topology descriptions at the ESRI help or Oracle help. (Maybe there's something there but I'm just not finding it yet!)   2. _Faces_ : I’ve played around with the fact that the face on the “north” side of A might also bounded also be by B, but not C. However, as you can see in the diagram here, that’s not always true. Imagine my diagram is an extract from a road network where A and C are arterial roads and B is a short dead-end road. I suspect there may not be a single term for what I’m trying to say; at a minimum I’d like to be able to describe such a relationship in a simpler way than I’ve done above. This is a platform-independent question. Right now, I’m just looking for the right words. Later I will try to implement the concept in python (pyqgis or arcpy) on a shapefile, so any answers with that end point in mind will be particularly interesting, but not necessary.


title: gdal_merge to merge hundreds of png files

text: I need to merge about 760 png files into one and am looking at using gdal for this. This is related to the following question Getting a constant background color for transparency               P:\2012\183_TownPlanning_Symbology\Working\Raster_Layer_Creation\tests\BA\25k>gd     al_merge -o ba_25k_gdal -of png -init "255 255 255" *.png      should work but gives > Format driver png does not support creation and piecewise writing. Please > select a format that does, such as GTiff (the default) or HFA (Erdas Imag > ine). When I switch to HFA or GTiff format and try               P:\2012\183_TownPlanning_Symbology\Working\Raster_Layer_Creation\tests\BA\25k>gd     al_merge -o ba_25k_gdal -init "255 255 255" *.png      I get > ERROR 4: `*.png' does not exist in the file system, and is not recognised as > a supported dataset name. > > Traceback (most recent call last): File "C:\OSGeo4W\bin\gdal_merge.py", line > 509, in sys.exit(main()) File "C:\OSGeo4W\bin\gdal_merge.py", line 392, in > main ulx = file_infos[0].ulx IndexError: list index out of range but the format seems correct according to http://www.gdal.org/gdal_merge.html How can I tell it to include all png files in the directory? I can pipe *.png > png_list.txt but can I get gdal to read a list from a text file? So does this mean merging png's in gdal is not possible? Are there other ways of achieving this? I have tried Global Mapper but that estimates a day to complete it and I am just running it. Erdas crashed. I am trying ArcMap 10.1 import into mosaic dataset and export as raster -still running but estimates a 35gb png file. This is from 161 tiles (25k resolution) over 65,000 sq km. But the source png's are only 73Mb.


title: Programatically write values to ArcPad Form control box

text: Using this blog post I've managed to store user names as global variables in ArcPad. The next step is to use those variables in form 'controls' aka form fields. Per the author's comments under his original quote I've attempted to go this route: >   * If you are writing the value to a field that is linked to a control on > the form, set the control value: > pForm.Pages("PAGE1").Controls("EMPLOYEE").Value = > Application.UserProperties("LoggedInEmployee") > Unfortunately I'm not getting anything written into the field when I open the form (this is an onLoad event handling script). I'm not getting an error message at this point, and the value in the control box is changing from to being totally empty. I have managed to 'echo' or 'alert' back the user-name with Application.MessageBox Application.UserProperties("LoggedInEmployee") associated with an onUnload event, so I know that it's stored properly. Any Ideas?


title: Dark Object Subtraction QGIS

text: I'm trying to do Dark Object Subtraction (DOS) for a Spot2 and 4. I have calculated the TOA reflectance and now want to get the at surface reflectance. Using QGIS I would like to know how I can do DOS in QGIS, can I build a model or use an equation to do this. I have tried doing 6S however SPOT2 and SPOT4 HRV1/2 is not listed in the Pre-defined satellite bands.


title: How to split lines on vertex?

text: I've searched a lot but didn't find an answer: When editing lines with QGIS, I often need to split a linestring exactly on one of his existing internal vertices, but I can't find a way; I've tried the tool "Split Features" thinking that it could snap on vertices, but I was wrong, and I've tried several ways of selecting the vertex with no luck... any solution?


title: Openlayers stylemap default style doesn't render until mouseover in IE 7/8

text: I'm using Drupal 7 Openlayers, and adding my own custom behavior. I have no trouble getting my code included, and it works well on other browsers. But in IE 7/8, the layer features render with a null style (solid white with a black outline) until I mouseover them. I have no visible errors in the console, and otherwise the map works fine. On mouseover I (correctly) get the highlight style for that feature, and when I release the hover I (correctly) get the default style for that feature. It's just when the page first loads that it has this problem. Here's the my code for applying the styles, from my custom Behavior's .js file.               Drupal.behaviors.openlayers_ranged_colors = {           attach: function(context) {       var data = $(context).data('openlayers');       if (data && data.map.behaviors.openlayers_ranged_colors) {         var styles = data.map.behaviors.openlayers_ranged_colors.stylemap;         var field = data.map.behaviors.openlayers_ranged_colors.field;         // Collect layers         var vector_layers = [];         for (var key in data.openlayers.layers) {           var layer = data.openlayers.layers[key];           // Only if it's a Vector layer           if (layer.isBaseLayer != true) {             // Apply stylemap from the Behavior             var styleMap = layer.styleMap;             styleMap.addUniqueValueRules("default", field[layer.drupalID], styles['default']);             styleMap.addUniqueValueRules("select", field[layer.drupalID], styles.highlight);             styleMap.addUniqueValueRules("temporary", field[layer.drupalID], styles.highlight);                  layer.redraw();             vector_layers.push(layer);           }         }


title: Openlayers WMS url construction

text: I am working on an extension to openlayers to pre-cache tiles in a local websql datastore to support mapping on mobile device when connectivity is dead. I would like to pre-seed my cache with a set of tiles. Currently this has to be done by visiting the area in question and using the cacheWrite control. This would be done by fetching the tiles in question from my wms mapserver and storing them locally in the websql database. These would then be available for the cacheRead control to use on tileerror. So, my question is: given a bounding box, is there an easy way for me to construct these calls using openlayers and hence populate the cache with the required tiles for a required zoom level? A sample single tile call is:-               "http://mymapserver/pub/mapcache/mytheme?LAYERS=my_layer&FORMAT=image%2Fpng&SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&STYLES=&SRS=EPSG%3A27700&BBOX=533371.66,169851.66,534011.66,170491.66&WIDTH=256&HEIGHT=256"      I can do this by panning the map about automatically, but the user experience is terrible ;-)


title: OpenLayers - changing layer type (wms to OSM) and WKT transformation problems

text: I've created the very simple tool which utilizes `OpenLayers.Format.WKT` with `DrawFeature - polygon` and `ModifyFeature` together, so I can create and modify, load and save polygons in/from WKT format. I've had it working on `OpenLayers.Layer.WMS`, but now I'd like to work on OSM map as I need the city names to be shown. The problem is that after changing the layer type, the polygon is loaded properly - it's displayed somewhere else than it is supposed to be. Here is my loading code (simplified)               var wmsLayer = new OpenLayers.Layer.WMS( "OpenLayers WMS","http://vmap0.tiles.osgeo.org/wms/vmap0?", {layers: 'basic'});          // when using this, it doesnt work properly: var wmsLayer = new OpenLayers.Layer.OSM( "OpenLayers OSM");          var features = wkt.read(data);     features.geometry.transform(map.displayProjection, map.getProjectionObject())     vectors.removeAllFeatures();     vectors.addFeatures(features);      I believe this is due to transformation problem, but I don't know how to fix it...


title: Using stat_density2d with coord_map in ggplot2

text: In my example I am plotting a shapefile for all England & Wales, however I'm interested in the greater London area. If I allow ggplot to simply default and display the entire shapefile then the heatmapping generated by stat_density2d isn't very useful for me, I need to change the scale to show the Greater London area. To address this I need to set the xlim and ylim to effectively zoom in. I can do this in one of two ways:   1. Use the xlim & ylim standalone functions in ggplot   2. Use xlim & ylim within coord_map The former causes the heatmapping to scale correctly (ie size of the tiles reduce to give a meaningful indication of density in a smaller area) but the shapefile data gets truncated to the plot leading to some very odd looking polygons. The latter causes the heatmapping to remain at the scale of the whole country (ie the tiles are massive and so not useful for showing activity density when the map is zoomed in) but the polygons display correctly. Is there a way to get ggplot to change the xlim/ylim values of a plot without truncating the shapefile/polygon data so that it also rescales the tiles generated by geom_density2d? I know I could fudge this by only plotting a subset of the shapefile polygons but this seems...messy. It would be way cooler if I could choose anywhere in the country to zoom/centre on and have everything scale correctly without having to chop and change the shapefile data. This is a follow up to my previous question: How do I create a heatmap in ggplot? **Edit 1:** I'm currently investigating a few undocumented arguments that stat_density2d passes to the kde2d 2d density function from the MASS package. There are a few of these but I get the impression a combination of n and lims may help. Basically if I pass a higher value to n (this dictates the number of rows and columns there are in the matrix produced) then the 'resolution' of the 2d density estimation increases at the expense of processing time (this is applied to the _entire_ bounding area of the shapefile rather than just the window I'm limiting the plot to using coord_map). I then use lims to try and limit the matrix to just the plotting area but that doesn't seem to have much effect.


title: What alternatives do I have for Directions API other than Google Maps and MapQuest?

text: I need an API, other than the two below that, given the starting point and possibly the ending point, will re-order a list of waypoints by the shortest, fastest or optimized route:   * **Google Directions API**   * **MapQuest Directions API**


title: Date field formatting as number in ESRI popup window

text: We are using the Park and Recreation Locator from ArcGIS for State Government and trying to put a date field into the popup info window. The field is pulling the correct data, but it's formatting it as a number (JavaScript milliseconds) instead of a date. There is a variable in the **config.js** file called `FormatDateAs` which has been set as a date according to ESRI's site `"MMM dd, yyyy"`, but has no discernible effect. From experience with some of ESRI's other templates, I tried adding in fields to the `InfoPopupFieldsCollection` such as `isDate` or variations on `Format` but none worked. ESRI's site has no additional documentation for formatting a date in this template. We need the field to remain a date type in the database for future changes we have planned, so simply changing the fields to strings containing the pre- formatted dates won't work. Does anyone know how we can get this field to format correctly in the code?


title: Pgrouting install debian 7.5 whezy

text: I need to install pgRouting on a Debian 7.5 64-bit server. I have PostgreSQL 9.1 and PostGIS 1.5 running. But somehow I can not install pgRouting. 1 - To use add-apt-repository ran the command:               apt-get install python-software-properties      2 - I entered the site georepublic `https://launchpad.net/ ~ / + archive / ubuntu / pgrouting-unstable` ppa to get the source to add to the list but only has to Ubuntu. Can anyone help me???


title: Populate geoext2 gridpanel from openlayer vector features

text: I need to populate a extjs datagrid with a openlayers features. The official example http://geoext.github.io/geoext2/examples/grid/feature-grid.html doesn't help me, beacuse the data is loaded from a JSON file, but I need the features parsed from a WFS layer. Colud anyone help me with a source code example (jsfiddle)?


title: Read large raster into lat/lon value pair to csv file using python and/or gdal?

text: I have rasters for all of the states in US and they are in 10m resolution. So they can be very large, containing a large number of pixels if you think of TX or AK. I want to create a csv file with each line representing the lat/lon value of the pixel, with comma separating lat and lon. So far the methods that I have tried for MN state:   1. gdal2xyz----running time is reasonable, but there is extra column for z value in the csv file. I need to later delete this column through code.   2. iterate each pixel in the raster, get it's lat and lon, and write the value into csv---really takes a long time... To sum up, I just want a csv file with the lat/lon value pair from raster. I feel like it is not hard to achieve. It is just the matter of the running time. I hope there could be some work around. Thank you.


title: QGIS DEM Output cell size does not match DEM input cell size

text: As my first post I offer a humble apology if this is such a basic question as to make any of you cringe. I feel like I have exhausted my search skills and read through the QGIS documentation frequently and thoroughly. My problem is with any of the hydrology analysis in GRASS through QGIS. I have a GRASS mapset that uses a 30m resolution DEM for my area. The extent of the mapset is set to include only the DEM (ie. no whitespace or no-data areas) and yet when I run r.fill.dir, r.basin or r.watershed I get back an image that has a resolution of nearly 73xxmetres. I am most used to ArcGIS where in the hydrology tools the user can define the working environment and set the output resolution to match the input resolution. Is there anyway to set that in GRASS or am I missing a basic step somewhere? I feel like I set up GRASS correctly as other non analytical map work is working fine.


title: Add layers from WMSCapabilitiesStore to MapPanel and select polygons

text: I created a store and added layers to a grid Panel, when you click add layer to map it works fine, the problem is that when I call the function controle(layer) I can't select polygons in the layer that is active. When I add layers manually using OpenLayers.Layer.WMS.... and I send the variable that has the layer to the function controle it works fine but from a store it doesn't. What am I doing wrong? Can anyone please help? The objective that I am trying to achieve is to add layers from several databases to a grid Panel, add the layers to a tree Panel and a Map Panel and click on the states or countries and with shift key pressed be able to select another layer and select another polygon from that layer and obtain a WKT text from all the polygons selected. Thanks :)                       // create a vector layer for selection of polygons             select = new OpenLayers.Layer.Vector("Selection", {                     styleMap: new OpenLayers.Style(OpenLayers.Feature.Vector.style["select"])                 });                  // create a new WMS capabilities store                 store = new GeoExt.data.WMSCapabilitiesStore({                     url: "http://localhost:9090/geoserver/wms?service=wms&version=1.1.1&request=GetCapabilities",                     autoLoad: true                 });                  // create a grid to display records from the store                 var grid = new Ext.grid.GridPanel({                     title: "WMS Capabilities",                     store: store,                     cm: new Ext.grid.ColumnModel([{                         header: "Name",                         dataIndex: "name",                         sortable: true                     }, {                         id: "title",                         header: "Title",                         dataIndex: "title",                         sortable: true                     }]),                     sm: new Ext.grid.RowSelectionModel({                         singleSelect: true                     }),                     autoExpandColumn: "title",                     renderTo: "capgrid",                     height: 300,                     width: 350,                     floating: true,                     x: 10,                     y: 0,                     bbar: ["->",                     {                         text: "Add Layer",                         handler: function() {                             var record = grid.getSelectionModel().getSelected();                             if (record) {                                 var copy = record.copy();                                 // Ext 3.X does not allow circular references in objects passed                                  // to record.set                                  copy.data["layer"] = record.getLayer();                                 copy.getLayer().mergeNewParams({                                     format: "image/png",                                     transparent: "true"                                 });                                 panel.layers.add(copy);                                 panel.map.zoomToExtent(                                 OpenLayers.Bounds.fromArray(copy.get("llbbox")));                             }                         }                     }]                 });                  // create a map panel                 panel = new GeoExt.MapPanel({                     renderTo: "mappanel",                     width: 350,                     height: 300,                     floating: true,                     map: map,                     x: 570,                     y: 0                 });                  // create a treePanel                 tree = new Ext.tree.TreePanel({                     renderTo: "tree",                     root: new GeoExt.tree.LayerContainer({                         text: 'Map Layers',                         layerStore: panel.layers,                         leaf: false,                         expanded: true                     }),                     enableDD: true, //                     width: 170,                     height: 300,                     floating: true,                     x: 380,                     y: 0                 });                  // function to getFeature              function controle(layer) {                 panel.map.addLayers([select]);                 control = new OpenLayers.Control.GetFeature({                     protocol: OpenLayers.Protocol.WFS.fromWMSLayer(layer),                     multipleKey: "shiftKey"                 });                 control.events.register("featureselected", this, function(e) {                     select.addFeatures([e.feature]);                 });                 control.events.register("featureunselected", this, function(e) {                     select.removeFeatures([e.feature]);                 });                 panel.map.addControl(control);                 control.activate();             }


title: ArcGIS - World Imagery Map Service - Force Higher Resolution Imagery

text: I am using the ArcGIS world imagery map service in an ArcGlobe project. Based on the scale and perspective of my view, imagery updates between different resolutions. Sometimes two resolutions are shown at once resulting in choppy imagery. ![enter image description here](http://i.stack.imgur.com/TO2NL.jpg) I can see that this Map Service is made up of many different resolutions that display at different zoom levels. Is there a map service that will only show the more detailed imagery below? I believe this is the I-Cubed 15m data. Is there a way to force arcGlobe to display the higher resolution imagery?


title: How useful is the ESRI Production Mapping extension for cartography?

text: I would like to find out how useful is the $9000? ESRI Production Mapping extension. According to ESRI is helps "Standardize and centralize detailed cartographic production with tools for creating and maintaining derived data, symbology, page elements, and maps." Can someone provide opinions on contour creation and labelling, layout features, styles (apparently stored within a geodatabase?) Area there any enhancements to the legend creation process?


title: How to export a mxd to a image with transparent background

text: I cannot get ArcMap 10.0 or 10.1 to "Export Map" to a png or gif with a transparent background -- no matter what I set the background and transparent colors to. Thus, I created the following script, which does the trick, but I'd rather be able to get the default tool in ArcMap to do it. Any better ideas? * * *               # Purpose: This script exports an MXD view to a png file with transparent background,     # which is problematic to do manually due to a Windows-related bug in the ArcMap UI.     # Warning: It will make any white features (e.g. the default background color) transparent.     #     # INSTRUCTIONS:     # When in ArcMap, viewing a saved mxd...     # 1. From the Geoprocessing menu, select 'Python'.     # 2. In the new Python window, right-click somewhere after >>>, and select 'Load...'     # 3. Select this script.     # 4. When this script appears, press Enter twice to run it.     # 5. A new [mxd-name].png with transparent background will appear in the same directory as the mxd.     # 6. When finished/satisfied, close the Python window.     # http://forums.arcgis.com/threads/65463-Zoom-to-tile-and-Export-to-JPEG-with-World-File     import arcpy     png_resolution = 300    # Can adjust DPI     # Set mxd and png name     mxd = arcpy.mapping.MapDocument("CURRENT")     pngname = mxd.filePath.replace(".mxd",".png")     df = arcpy.mapping.ListDataFrames(mxd)[0]     # Calculate ideal image size from the mxd view and desired DPI     png_width = int((df.extent.XMax - df.extent.XMin) * png_res * 12 / df.scale)     png_height = int((df.extent.YMax - df.extent.YMin) * png_res * 12 / df.scale)     #Export PNG     print "Exporting: {0}\n{1} x {2} pixels, {3}-dpi".format(pngname,png_width,png_height,png_res)     arcpy.mapping.ExportToPNG(mxd, pngname, df, df_export_width = png_width, df_export_height = png_height, resolution = png_res, transparent_color = "255, 255, 255")     #Clean up     del df, pngname, png_res, png_width, png_height, mxd


title: How to export data frame as PNG in ModelBuilder?

text: I've created a model that when finished creates two symbolized raster layers and adds them to the dataframe. I would now like each layer to be exported as a unique png file. Is there a tool that I can add to my model to complete this task? See the model below. I would actually like to create two pngs, one of phaseFinished.lry and another of activityFinished.lyr. ![enter image description here](http://i.stack.imgur.com/8OO2i.jpg)


title: How to programmatically remove a donut hole from Polygon in arcobjects

text: ![Polygon with a hole in it](http://i.stack.imgur.com/bMhEN.jpg) I have a polygon `IPolygon4` with a hole in it, as you can see. How can I get rid of the donut hole, and have the whole thing be one simple polygon with 5 vertices. Here are the constructed WKT points, I projected to GCS_WGS_1984 so it would be more readable. > MULTIPOLYGON (((-93.948870475981934 46.388792874272333, -94.07064721379696 > 46.3298500944495, -94.292075670722767 46.511771884217374, > -94.285856486919386 46.334338634913607, -94.043790864983947 > 46.2678965446348, -93.96530675201366 46.327504695188082, -93.948870475981934 > 46.388792874272333)), ((-93.948870475981934 46.388792874272333, > -93.914198248635017 46.517621156043539, -94.292075670722767 > 46.511771884217374, -93.948870475981934 46.388792874272333)))


title: Why do edits made in attribute table not show in associated layer? QGIS 1.8.0

text: I have compiled a heat map layer to display over a background raster map. I discovered that the heat map layer was incorrectly positioned because I had made an error in the coordinate entry fields of the imported data table. I therefore selected the attribute table of the heatmap layer, toggled editing on and manually amended the coordinates to the correct format. I then clicked save edits. However, although the attribute table now shows the correct coordinates, the heatmap layer is still in it's original incorrect position. Is there a simple solution to this? I am hoping I do not have to amend the original excel spreadsheet and rejoin the table.


title: pulling an excel sheet with =rand(between) in a Python while loop, and exporting results as .dbf

text: To preface my question: I am very new to stack overflow, and relatively new to Python. I am working on setting up a sensitivity analysis. I am working with 40 parameters that range from 0.1 - 1. My analysis requires simultaneously varying these parameters by +-0.1 roughly ~500 times. These values will then be fed into an ArcGIS tool. So, I need 500 sets of random values, for 40 parameters. These parameter values will then be compared to the output of the tool, to see which parameters the model is the most sensitive to. I've set up an Excel sheet that will randomly calculate these values each time it's opened, but the issue is that they need to be in .dbf format to be read by the ArcGIS tool. I have set up a while loop (for 10 iterations to start, but will need to be ~500) and tried two different methods, in hopes that I could automate the process of calling the .xls to generate random numbers, and then exporting it to .dbf. The first, arcpy.CopyRows_management correctly exported to .dbf. The issue was that the output was the exact same for each iteration, and instead of having values of 0.1, 0.2, 0.3 etc, it contained values of 0.22, 0.37, 0.68 etc. It wasn't to the tenths, even though that was specified in the formulas in the .xls. I also tried arcpy.TabletoTable_conversion but that was throwing the ERROR 999999: Error executing function. I am open to all kinds of suggestions. Perhaps there is an easier way to randomly sample and export results to .dbf in Python. This does not need to be done using arcpy, but that is all I've really worked with. I really appreciate any help that is provided! Thanks for your time.               i = 0     while i < 10:         # Set run specific variables         lulc = "D:\\SARuns\\lulc_nosoils_rand.xls\\lulc_nosoils$"         folder = "D:\\SARuns"              print "Reading lulc"              newlulc = "D:\\SARuns\\lulc_nosoils_rand"+str(i)+".dbf"              print "Reading newlulc"              # Copy rows is copying it to a dbf, but the values inside          # are the same for each run. And, none are correct.             arcpy.CopyRows_management(lulc, newlulc)              # Table to table should work. But isn't.             # arcpy.TableToTable_conversion(lulc, folder, newlulc)              print "Converting table"             i+= 1