title: How to change default zoom in the Google embedded map?

text: I am having difficulties to change the default zoom level when I open the iframe below in my site. Can anyone please help me to set zoom level when following iframe is opened?               <iframe src="http://mapsengine.google.com/map/u/0/embed?mid=zYk_BmXXTdxw.k83snDbX5SJc" width="745" height="480"></iframe>      Thank you so much everyone in advance. PS: This is the newest version of Google Maps.


title: How to read vector datasets using GDAL library

text: Up to now, I have used GDAL in order to read raster information and geotiffs and show them with the help of OpenGL. Now I have a vector dataset, and by that, I mean I have the files with types: **DBF File , PRJ File , SBN File , SBX File , SHP File , LOCK File , XML Document , SHX File** I want to read the vector information **(e.g vertexes of the geometries and geometrical and topological relationships between feature)** and then show them using OpenGL.   * What part of GDAL library should I study?   * Can you introduce me a tutorial link that has fully covered what I want to read from the vector dataset?


title: please Introduce me a free 3rd party library appropriate for showing Geospatial information

text: Up to know I have used GDAL to extract raster data from geotiffs and then opengl's used to show that raster data.   But when it comes to extract vector data with the help of ogr library, what library would you suggest to show these vector data in visual-studio picture controls opengl that is more appropriate than opengl? * * * I've found out from your comments that you have thought I'm searching for a way to read vector data. this is the subject of my another question. Here, I'm just asking suppose that we have extracted all the geometries. Which environment is the best to show these geometries?


title: Is there a work-around for a completely_within spatial join in ArcGIS? This is bugged (confirmation by Esri Inc)

text: I tried to perform a 1-to-many spatial join between polygons (target feature) and lines (join features) to join the lines that are completely covered by a polygon to that polygon. (the polygon is actually a result of buffering the lines with a dissolve option) I used the option 'completely contains', but this returned a background processing error. After sending parts of my data to Esri Inc, I got the confirmation that there is indeed a bug they need to fix (#NIM086228). Anyway, my project needs to be completed. In search of a work-around, I already tried to turn the lines into polygones (small buffer around it), but it returns the same error. When trying the contains_clementini (or actually all contains and within options) suddenly give me an error even before running it. The error is 'Relationship invalid for selected layers'. Here is what I need to do: > I have (pipe)lines that I want to 'merge' when they have common > characteristics. So I create a buffer and let it dissolve when the value of > those characteristics is equal. A polygon layer is created, but there is not > link back to the underlying lines that are the basis of the polygon. In > order to have these ID's, I wanted to perform the spatial join with match > option 'completely contains Has someone any experience with the last type of invalid error messages? Or another work-around for me? Thank you already! Mieke


title: Scalar relationships ArcGIS

text: I have point data which has a number of fields/columns, each of which is at a different scale. I'm looking to find a way of constructing a very simple index which would combine this data. Say I have 3 levels of scale (so 3 fields) - small, medium, large. Each of the points either has or hasn't got a text string at each of those levels. I'm looking to assign a number (or not) depending on whether the text string is present or not, and aggregate that data. Basically I want to see if there are patterns in the presence/absence of text strings at these different scales (possibly re-combining the value given to each of them in more than one way). I'll then figure out a way of showing this map format (Arc). I don't think I've been very clear but please ask if you need clarification. Here's a sample table with something along the lines of what I was thinking of doing (in grey). I'm not sure this will give me the outcome I need - ideally each column should have a different weighting reflecting relative size. But should give a rough idea of the sort of thing I was thinking of doing. ![enter image description here](http://i.stack.imgur.com/y01Xt.png)


title: How to Interpolate through ModelBuilder?

text: I would like to create an iterated interpolation tool through ModelBuilder. I have an excel file containing different columns and each column represents a data that will be interpolated using Spline Tool. How can I do this through ModelBuilder? I am using ArcGIS 10. ![enter image description here](http://i.stack.imgur.com/d0t0m.jpg) Image contains point data with lat/long and sample# are the columns to be interpolated.


title: GPS suitability

text: I was wondering what type of GPS/GNSS receiver should I use for recording and mapping trees in a particular farm (inventorying trees). Will handheld GPS like Garmin do the job? Or should I go for the expensive RTK system? Please help. *Remark: The area is a mountainous region.


title: How to package up a Python script for an end user to run

text: I have a Python script that I have made to some basic data management tasks. I want to package up a Python script in a nice user friendly way for an end user to run, they do not know how to run Python scipts and I need them to input some parameters in an easy way i.e. say a filename - whats the best way to do this. Cheers


title: How to cache online tiles on a development machine?

text: I am developing a slippy map web application using openlayers. Right now my tiles come straight away from OSM and google, etc. For the moment I always work on the same area for development purposes and I would like to have those tiles cached on my dev machine so I can work on the web app offline yet still see some map and to reduce network consumption. What I would like is to set up some sort of proxy which would just cache tiles I watch and if they are not on the local machine download them or show a dummy tile if I am offline. I thought tilecache was the tool for this, but I get quite confused trying to set it up. It seems more appropriate for mapserver generating tiles from map data. I know firefox does some caching, but I would like some more control and be able to clear that cache to make sure I see the last update of my work. What is the most straightforward way to achieve this?


title: Elevation statistics for a given Area

text: I am attempting to determine the best way to calculate the Minimum and Maximum elevation for a defined area. The area is defined by an `IPolygon` The elevation is defined as an `ISurface` ## **Possible Solutions:**   1. The only way I can think of is to check every cell of the surface. Is there a method that performs this?   2. In the comments, @whuber, suggested a Zonal Maximum. Expanding on this idea, I've begun looking into the `IZonalOp`. _Are there other options?_


title: Feature linked annotation subtype removal in migrating to ArcGIS 10

text: I have a question about emergent properties of deleting an annotation subtype. If we delete an annotation subtype how will this effect the annotation feature class, will there be any changes to the visible annotation the active data frame? I haven’t seen any changes with the removal of the subtype, but it seems too good to be true. Background info: In our enterprise implementation, we are in the process of upgrading from a 9.2 to 10 data format, to help achieve this we came up with a Python script, turns out there are quite a few databases that fail at the arcpy.UpgradeGDB_management() method. The exception is:"The feature class' default subtype code cannot be retrieved or is invalid” In order to resolve this issue we removed the Annotation’s subtype and the script will run just fine. Thanks in advance.


title: Why are ArcGIS 10.2.1 for Server services stuck in perpetual Stopping... state?

text: Running into a problem repeatedly with our ArcGIS Server installation where all services are in a perpetual state of "Stopping..." whenever we do an upgrade (just upgraded from 10.2 to 10.2.1). The same thing occurred when we changed the ArcGIS Server account from a local account to a domain account. The services cannot be stopped and started manually. Restarting the ArcGIS Server service does not seem to help either. Rebooting the server that ArcGIS Server is running on does not help with this issue either.


title: How can I change the size of the infowindow in CartoDB?

text: I'm working with info windows, but I'd like them to be a bit bigger than the defaults in CartoDB. To display pictures and additional information. How can I make it so they are bigger? Is there any way in the GUI to do this? Is there a more advanced way that has some template I could easily apply?


title: Multiple geometry in a cursor SQL

text: I have a cursor which I create from a table which has one geometry field only. But in my cursor I would be using two geometry fields as I need to calculate two different areas.Below is the sql statement which I am using . After that I need to insert those two geometries in the table which also has two geometry columns. But when I try to do the same it gives me an error **"inconsistent datatypes: expected NUMBER got MDSYS.SDO_GEOMETRY""** How could I fix this particular geometry issue.                 Declare          Cursor c1 is           select 0 as BG_ID,store_number ,geometry , geometry bg_geometry,0 as total_area          from   TMCS_All_Stores_TA  where CLIENT_ID= 1 and ORG_ID = 1 and  rownum <= 1;         Type C1_TAB_TYPE is table of c1%ROWTYPE;                c1_list c1_TAB_TYPE;               Cursor c3 is           select BG_ID,0 as BLK_ID ,store_number,0 as total_area_blk           from   TEMP_PROP_AREA_BGLEVEL ;          Type C3_TAB_TYPE is table of c3%ROWTYPE;                c3_list c3_TAB_TYPE;                      Begin                     For r1 in c1                       Loop                 Select  BG_ID,store_number,geometry,bg_geometry,total_area                  Bulk Collect into c1_list from(                      SELECT a.BG_ID,r1.store_number,(SDO_GEOM.SDO_INTERSECTION(a.geometry,r1.geometry, 0.005)) GEOMETRY,r1.bg_geometry bg_geometry,                      ( (SDO_GEOM.SDO_AREA(SDO_GEOM.SDO_INTERSECTION(a.geometry,r1.geometry, 0.005), 0.005, 'unit=sq_mile') /SDO_GEOM.SDO_AREA((a.geometry), 0.005, 'unit=sq_mile'))*100) Total_area                      FROM   TMCS_BG_STI_DATA a                       Where sdo_relate(a.geometry, r1.geometry, 'mask=anyinteract') = 'TRUE');                                             For j in 1..c1_list.count loop                                        DBMS_OUTPUT.PUT_LINE('block group geometry--> '||c1_list(j).store_Number);                                        insert into   TEMP_PROP_AREA_BGLEVEL                                        values c1_list(j);                                         End loop;                             For j in 1..c1_list.count loop                         DBMS_OUTPUT.PUT_LINE('area--> '||c1_list(j).total_area);                       if c1_list(j).total_area !=100 then                                             select BG_ID,BLK_ID ,store_number,total_area_blk                                 BULK COLLECT INTO c3_list from                              (                              SELECT c1_list(j).BG_ID BG_ID,a.BLK_id,c1_list(j).store_number store_number,                              ( (SDO_GEOM.SDO_AREA(SDO_GEOM.SDO_INTERSECTION(a.geometry,c1_list(j).geometry, 0.005), 0.005, 'unit=sq_mile') /SDO_GEOM.SDO_AREA((r1.bg_geometry), 0.005, 'unit=sq_mile'))*100) total_area_blk                              FROM   TMCS_BLK_STI_Q213_DATA a                               Where sdo_relate(a.geometry, c1_list(j).geometry, 'mask=anyinteract') = 'TRUE'                              And BG_ID = c1_list(j).BG_ID);                                               For i in 1..c3_list.count loop                                       --   DBMS_OUTPUT.PUT_LINE('geo_id--> '||c3_list(i).Blk_id);                                       --    DBMS_OUTPUT.PUT_LINE('area--> '||c3_list(i).total_area_blk);                                          insert into   temp_Prop_area_BG                                          values c3_list(i);                                          End Loop;                            End if;                   End Loop;                        Insert into   temp_prop_cal                   (bg_id,  store_number, total_area)                   SELECT bg_id,  store_number,  SUM(total_area)                   FROM        (                   SELECT bg_id,  store_number,  total_area as total_area          FROM   TEMP_PROP_AREA_BGLEVEL a          WHERE store_number='336655'          UNION ALL          SELECT bg_id,   store_number,  total_area_blk as total_area          FROM   temp_Prop_area_BG b          WHERE store_number='336655')        GROUP BY bg_id, store_number;              End Loop;           End;


title: How to populate ESRI GeoPortal programmatically?

text: I'm trying to find either some documentation or an example of how to manually populate the GeoPortal metadata datastore. I'm using a MySQL backend. Given a ISO XML file I can manually harvest it but this is tedious and doesn't fit with our model. I can easily write a script to populate entries in the SQL database but it appears that there is more to it than just adding entries to SQL, there is a part that goes into the Lucene engine as well. If anyone has any experience or can provide a simple example, how-to or preferably some documentation on populating the datastore it would be greatly appreciated. The goal is to have GeoPortal as an external catalog of locally registered datasets, the trouble is registering the datasets with the GeoPortal catalog in a preemptive manner rather than a reactive (having an administrator manually "harvest" or ingest resources).


title: Trying to write a python script to crawl a directory and list all datasources

text: I'm editing a scrip from ESRI to try and create a list of all datasources within MXD's. I've got the script to crawl the file tree and find all the info, but I would like it to output the data for each layer/datasource on a single line, so it is easy to sort in excel. If anyone could take a look at the code I have so far and let me know how to accomplish that, I would appreciate it.               import arcpy, datetime, os          try:          arcpy.gp.overwriteOutput = True          #Read input parameters from GP dialog     folderPath = arcpy.GetParameterAsText(0)     output = arcpy.GetParameterAsText(1)          #Create an output file     outFile = open(output, "w")          #Loop through each MXD file     count = 0     for root, dirs, files in os.walk(folderPath):         for filename in os.listdir(folderPath):           fullpath = os.path.join(folderPath, filename)           if os.path.isfile(fullpath):             if filename.lower().endswith(".mxd"):                      #Reference MXD                 mxd = arcpy.mapping.MapDocument(fullpath)                 count = 1                      #Write MXD data to file                 MapDoc = os.path.basename(mxd.filePath)                 MapDocPath = mxd.filePath                       #Reference each data frame and report data                 DFList = arcpy.mapping.ListDataFrames(mxd)                 for df in DFList:                     #Format output values                     if df.description == "": descValue = "None"                     else: descValue = df.description                               #Reference each layer in a data frame                     lyrList = arcpy.mapping.ListLayers(mxd, "", df)                     for lyr in lyrList:                         lyrName = lyr.name                          if lyr.supports("dataSource"):                             lyrDatasource = lyr.dataSource                                                      else: lyrDatasource = "N/A"                      outFile.write("Map Document, MXD Path, DataFrame List, Description, Layer name, Layer Datasource")                 str = ",";                 seq = (MapDoc, MapDocPath, DFList, descValue, lyrName, lyrDatasource);                 outFile.write(str.join(seq))                      del mxd                                outFile.close()          #Open resulting text file     os.startfile(output)          #Delete variables that reference data on disk     del folderPath, output, outFile, fullpath          except Exception, e:     import traceback     map(arcpy.AddError, traceback.format_exc().split("\n"))     arcpy.AddError(str(e))      The error I recieve when I try to run this is as follows:               \MultiMXDReport_crawl.py", line 60, in <module>     outFile.write(str.join(seq))     TypeError: sequence item 2: expected string or Unicode, list found     <type 'exceptions.TypeError'>: 'str' object is not callable     Failed to execute (MultiMXDCrawler).


title: Select by Location tool

text: I am using the Select by Location tool to select certain polygons on my shape- file layer. I start by having a given point as my Selecting Feature and that its relationship is completely_ contains. After I run the tool the first time, I get the polygon feature that my point is in (lets call it Polygon A). When I run another select by feature tool that its relationship is boundary_touches my Polygon A, I get a layer that also contains Polygon A and the features that share boundaries with it. Is there a way I can prevent this from happening? Meaning that every time I run the select by location tool it doesn't select polygons from the previous time. Or is there a tool or a way that will allow me to join multiple layers together and remove duplicate features? I have tried the Union tool as well as the Intersect tool. Any suggestions? Thanks


title: Using American Fact Finder to download data for all block groups in a state at at time?

text: I am having issues using the advanced search on the American Fact Finder to download census block data for a whole state at a time. The best option that it gives me, when trying any route, is to select geographies of all census block groups within a county within a state. This leaves me with the option then of selecting every county, but in some states (like Texas) that have more than 250 counties, I am out of luck with this method. Does anyone have a solution for this? To duplicate my error:   * Visit http://factfinder2.census.gov/faces/nav/jsf/pages/index.xhtml   * Then select Advanced Search,   * Then select on the left "Geographies",   * Then select "all geographic types" radial button,    * Then in the dropdown, select Block Group or part,   * Then select state in the next drop down, and you will start to see the issue. Also, if instead of using the list search and you use the tab that says "name" I run into the same issue. Does anyone have any experience with American Fact Finder downloads? I would appreciate any suggestions, Thanks! Also, I am not asking the same question here as I am looking to see if this is possible just by accessing AFF not their ftp. Gathering data for a large number of census block groups


title: topology problem

text: Can anyone tell me how to fix an error after running the topology rule in ArcGIS? I tried to subtract and create a new layer, but a failed warning appears every time.


title: Symbology of contour lines in ArcMap

text: I have contour lines with height values from 0 to 1030 (sorted by height 110, 115, 120, 125). Values with in-between number (like 5, 15, 115, 135...) should be grouped and colored with same color, and other values should be grouped this way 10, 60, 110, 160, 210, 260... (yellow), 20, 70, 120, 170, 220... (red), 30, 80, 130, 180, 230... (cyan) etc. I hope you get it despite my bad English :) So, I don't know how to change symbology in ArcMap to color them that way. Can somebody help me? Thank you.


title: ArcGIS Server 10 and Desktop 9.3 interoperability

text: I want to install ArcGIS Server 10 but our users only have ArcGIS Desktop 9.3.1. So does ESRI support:   * a) Access/use AGS services from 10 in 9.3.1?   * b) Manage (create services etc) AGS 10 from ArcCatalog 9.3.1?  I know you can use the web manager but our users are used to manage AGS from desktop.


title: How can I use fortify() to create a filtered R data frame from a shape file?

text: I am in the process of building a `ggplot` choropleth map of population in administrative areas in Wales. I have downloaded the Boundary-Line data from the Ordnance Survey and extracted what seems to be the right shapefile (community_ward_region.shp). Using R, I have got as far as reading in the shapefile.               require(maptools)     shape <- readShapePoly(wards)     str(shape)      Which gives me this promising output:               Formal class 'SpatialPolygonsDataFrame' [package "sp"] with 5 slots       ..@ data       :'data.frame': 1690 obs. of  4 variables:       .. ..$ NAME      : Factor w/ 1507 levels "Abbey Cwmhir",..: 969 90 111 200 441 477 1455 249 255 305 ...       .. ..$ DESCRIPTIO: Factor w/ 4 levels "COMMUNITY","COMMUNITY WARD",..: 2 2 1 1 2 2 2 2 1 1 ...       .. ..$ COMMUNITY : Factor w/ 858 levels "Abbey Cwmhir",..: 67 67 81 128 152 152 152 152 157 190 ...       .. ..$ FILE_NAME : Factor w/ 23 levels "ABERTAWE_-_SWANSEA",..: 1 1 1 1 1 1 1 1 1 1 ...       .. ..- attr(*, "data_types")= chr [1:4] "C" "C" "C" "C"       ..@ polygons   :List of 1690       .. ..$ :Formal class 'Polygons' [package "sp"] with 5 slots       .. .. .. ..@ Polygons :List of 1       .. .. .. .. ..$ :Formal class 'Polygon' [package "sp"] with 5 slots       .. .. .. .. .. .. ..@ labpt  : num [1:2] 259009 188524       .. .. .. .. .. .. ..@ area   : num 1923892       .. .. .. .. .. .. ..@ hole   : logi FALSE       .. .. .. .. .. .. ..@ ringDir: int 1       .. .. .. .. .. .. ..@ coords : num [1:1629, 1:2] 259413 259420 259427 259427 259432 ...       .. .. .. ..@ plotOrder: int 1       .. .. .. ..@ labpt    : num [1:2] 259009 188524       .. .. .. ..@ ID       : chr "0"       .. .. .. ..@ area     : num 1923892      Now if I do this:               bar <- fortify(shape, region = "NAME")      I get a nice data frame called `bar` that looks pretty much as I expected:               > str(bar)     'data.frame':   4744053 obs. of  7 variables:      $ long : num  302962 302970 302974 303013 303015 ...      $ lat  : num  280066 280076 280078 280097 280105 ...      $ order: int  1 2 3 4 5 6 7 8 9 10 ...      $ hole : logi  FALSE FALSE FALSE FALSE FALSE FALSE ...      $ piece: Factor w/ 29 levels "1","2","3","4",..: 1 1 1 1 1 1 1 1 1 1 ...      $ group: Factor w/ 1762 levels "Abbey Cwmhir.1",..: 1 1 1 1 1 1 1 1 1 1 ...      $ id   : chr  "Abbey Cwmhir" "Abbey Cwmhir" "Abbey Cwmhir" "Abbey Cwmhir" ..      However, this is a large data frame and `ggplot` runs out of puff when trying to display it. In reality I only want to look at one area at a time. It looks as if the `FILE_NAME` factor in the shape object is what I want as it mostly corresponds to counties and the major conurbations.               > unique(shape@data$FILE_NAME)      [1] ABERTAWE_-_SWANSEA      [2] BLAENAU_GWENT_-_BLAENAU_GWENT      [3] BRO_MORGANNWG_-_THE_VALE_OF_GLAMORGAN      [4] CAERDYDD_-_CARDIFF      [5] CAERFFILI_-_CAERPHILLY      [6] CASNEWYDD_-_NEWPORT      [7] CASTELL-NEDD_PORT_TALBOT_-_NEATH_PORT_TALBOT      [8] CONWY_-_CONWY      [9] GWYNEDD_-_GWYNEDD     [10] MERTHYR_TUDFUL_-_MERTHYR_TYDFIL     [11] PEN-Y-BONT_AR_OGWR_-_BRIDGEND     [12] POWYS_-_POWYS     [13] RHONDDA_CYNON_TAF_-_RHONDDA_CYNON_TAFF     [14] SIR BENFRO - PEMBROKESHIRE     [15] SIR_BENFRO_-_PEMBROKESHIRE     [16] SIR_CEREDIGION_-_CEREDIGION     [17] SIR_DDINBYCH_-_DENBIGHSHIRE     [18] SIR_FYNWY_-_MONMOUTHSHIRE     [19] SIR_GAERFYRDDIN_-_CARMARTHENSHIRE     [20] SIR_YNYS_MON_-_ISLE_OF_ANGLESEY     [21] SIR_Y_FFLINT_-_FLINTSHIRE     [22] TOR-FAEN_-_TORFAEN     [23] WRECSAM_-_WREXHAM     23 Levels: ABERTAWE_-_SWANSEA ... WRECSAM_-_WREXHAM      Q. How can I select only a subset of the data from the `shape` object I extract from the shapefile? For example, only the `POWYS_-_POWYS` parts? If I can somehow include the `FILE_NAME` in the data frame that is created with `fortify` then I could easily subset the `bar` data frame but I don't know how to do that. Or is there a way to use `fortify` to extract only parts of the object?


title: How to pass an output stored in a single variable into another script using arcpy

text: I am running set of processes each with different outputs in a python script. Is there a way to set a variable as global so that just the output of that single step can be accessed from other scripts?


title: Python Error while installing Cartogram

text: I am quite new to GIS and qGIS especially. I would like to install cfarmers cartogram plugin. From what I understand, this has to be done manually - so I 've followed the instructions from RobinLoveLace (see link below). Where has the cartogram plugin for QGIS gone? When I try to click the icon for the Cartogram Creator, I get an error message (see below). Any tips on how to get it to work? > An error has occured while executing Python code: > > Traceback (most recent call last): File > "C:/PROGRA~2/QGISCH~1/apps/qgis/./python/plugins\cartogram\plugin.py", line > 136, in run d = doCartogram.Dialog(self.iface) File > "C:/PROGRA~2/QGISCH~1/apps/qgis/./python/plugins\cartogram\doCartogram.py", > line 61, in **init** for name, layer in layermap.iteritems(): File > "C:/PROGRA~2/QGISCH~1/apps/qgis/./python/plugins\cartogram\doCartogram.py", > line 61, in **init** for name, layer in layermap.iteritems(): File > "C:\PROGRA~2\QGISCH~1\apps\Python27\lib\bdb.py", line 49, in trace_dispatch > return self.dispatch_line(frame) File > "C:\PROGRA~2\QGISCH~1\apps\Python27\lib\bdb.py", line 68, in dispatch_line > if self.quitting: raise BdbQuit BdbQuit > > Python version: 2.7.4 (default, Apr 6 2013, 19:54:46) [MSC v.1500 32 bit > (Intel)] > > QGIS version: 2.4.0-Chugiak Chugiak, 8fdd08a > > Python path: ['C:/PROGRA~2/QGISCH~1/apps/qgis/./python/plugins\processing', > 'C:/PROGRA~2/QGISCH~1/apps/qgis/./python', u'C:/Users/HEC/.qgis2/python', > u'C:/Users/HEC/.qgis2/python/plugins', > 'C:/PROGRA~2/QGISCH~1/apps/qgis/./python/plugins', > 'C:\PROGRA~2\QGISCH~1\bin\python27.zip', > 'C:\PROGRA~2\QGISCH~1\apps\Python27\DLLs', > 'C:\PROGRA~2\QGISCH~1\apps\Python27\lib', > 'C:\PROGRA~2\QGISCH~1\apps\Python27\lib\plat-win', > 'C:\PROGRA~2\QGISCH~1\apps\Python27\lib\lib-tk', 'C:\PROGRA~2\QGISCH~1\bin', > 'C:\PROGRA~2\QGISCH~1\apps\Python27', > 'C:\PROGRA~2\QGISCH~1\apps\Python27\lib\site-packages', > 'C:\PROGRA~2\QGISCH~1\apps\Python27\lib\site- > packages\GDAL-1.11.0-py2.7-win32.egg', > 'C:\PROGRA~2\QGISCH~1\apps\Python27\lib\site-packages\PIL', > 'C:\PROGRA~2\QGISCH~1\apps\Python27\lib\site- > packages\jinja2-2.7.2-py2.7.egg', > 'C:\PROGRA~2\QGISCH~1\apps\Python27\lib\site- > packages\markupsafe-0.23-py2.7-win32.egg', > 'C:\PROGRA~2\QGISCH~1\apps\Python27\lib\site- > packages\python_dateutil-2.2-py2.7.egg', > 'C:\PROGRA~2\QGISCH~1\apps\Python27\lib\site- > packages\pytz-2014.2-py2.7.egg', > 'C:\PROGRA~2\QGISCH~1\apps\Python27\lib\site-packages\win32', > 'C:\PROGRA~2\QGISCH~1\apps\Python27\lib\site-packages\win32\lib', > 'C:\PROGRA~2\QGISCH~1\apps\Python27\lib\site-packages\Pythonwin', > 'C:\PROGRA~2\QGISCH~1\apps\Python27\lib\site- > packages\Shapely-1.2.18-py2.7-win32.egg', > 'C:\PROGRA~2\QGISCH~1\apps\Python27\lib\site-packages\six-1.6.1-> > py2.7.egg', 'C:\PROGRA~2\QGISCH~1\apps\Python27\lib\site- > packages\wx-2.8-msw-unicode', > 'C:\Users\HEC\.qgis2\python\plugins\mmqgis/forms', > 'C:\PROGRA~2\QGISCH~1\apps\qgis\python\plugins\fTools\tools'] > The plugin for Cartogram Creator for qGIS is downloaded from here: https://github.com/cfarmer/cartogram-plugin/tree/master


title: Cartogram plugin for leaflet or OSM

text: Does anyone know if there is a plugin to make cartograms for leaflet ? I have sets of geotagged Twitter datas and I'd like to distord the map according to the density of tweets of every area of the world. As cartograms are mostly changing the tiles I'm not sure if Leaflet is the most interesting plugin to do so (as it's mostly done to display datas on a static map layer like an OSM one), but as I begun to use it a lot for other purposes it would be easier to implement in my code... Another approach would be to find a plugin to transform directly the OSM layer without using Leaflet, but I don't know if that exists too. Thx in advance !


title: Cannot view published map content on ArcGIS Online outside of firewall

text: I successfully published and added a map onto ArcGIS Online using ArcGIS for Server, which is located on a secure server. I can view the map just fine within the server environment (with localhost:6080) and from my work desktop (with [servername]:6080). However, when I try to view the same map from home, it fails to find http://[servername]/arcgis/rest/services/[mapname]/MapServer. Is there any way to bypass this firewall issue without having to obtain a web server? I do not have server admin rights, and obtaining permission for installing a web server would take longer than the timeframe that I'm working under. (Would a web server even resolve this issue?)


title: GIS Analysis - best way to combine and display demographic data?

text: I am currently undertaking a demographic research project for my employer. We have the data on income, age, sex, language, transport etc and I have mapped most of it and created some thematics of said variables. Whilst the maps are perhaps adequate from a raw data representation perspective, I am looking for a technique or thematic style which displays statistical and demographic data in a more intuitive and interpretive fashion. Here is an example of what i would like to improve: ![http://i894.photobucket.com/albums/ac142/mellingham/Income_zps95372b84.png](http://i.stack.imgur.com/aXzp1.png) **Is there a specific technique or method that people are using which displays thematic data which is more easily interpreted by a non GIS user?** Similar to this **question** but perhaps with more technical descriptions on how to achieve some of the better cartographic outputs. **EDIT: Thanks Ian who suggested using rasters to generate a heat map type approach:** ![Heat Map of Non-English speaking families](http://i.stack.imgur.com/mWlnq.png)


title: How to validate geometry in openlayers?

text: How can I validate the geometry in Openlayers? I need to check " **MULTILINESTRING()** " (which is empty) is valid/not valid. In postGIS, for Geometry validation we are using " **isValid()** " function. Like this, is there any options in Openlayers?


title: Efficient algorithms for recovering watershed boundaries?

text: This is somewhat related to a previous question I asked: What Algorithm is used by ArcGIS Watershed tool? I was wondering if there is a more efficient way to calculate catchment areas if you are interested only in the polygon defining the catchment area (rather than all cells contributing to the catchment area). Doing a recursive search seems very expensive and would return a lot of information we don't want (namely each cell within the area). Assuming the same starting point - a raster grid and a pour point/outlet (I leave it open for now whether the raster is a DEM, Flow Direction or Flow Accumulation). My first thought was to simply follow the 'steepest path' upstream of the outlet point, which would quickly take us to the 'edge'. But I could not think of a way to follow the 'edge' around. Are there any other possibilities?


title: Linking up GPS points taken from rail network

text: Bit of a newbie here, trying to work out whether any PostGIS functions can help with my task. Apologies in advance if my use of GIS terminology is not up to scratch yet. I have a large database of GPS points, ie lat/long coordinate pairs (large being "millions" of points). Currently they have no relational information, so they are not linked together. However, they were collected over the UK rail network, typically every few metres. So all points map to rail tracks. I need to process these points and link them together, so they map onto real railway lines. I want to identify which GPS points are buffers (ie ends of track) and points (in the rail sense - ie where track splits) as well as those that are "just sections of track". This needs to be a process that runs automatically with no user intervention. There are many areas where multiple tracks run in parallel, or areas such as large stations with lots of points & sidings, where there are large numbers of GPS points in a small area. Being rail, there is a limit to how much the line will curve, so that must be a key feature to aid analysis. As the GPS points are close together, it should be possible to see which form lines that could be track, and dismiss those that turn too sharply (which must be points from parallel tracks). Initially I had been planning to write some java code that would pick GPS points from the DB that fall within a limited geographical area of each other, and look for GPS points within that region that continued in a straight-ish line. The add/update "previous node id" and "next node id" fields for each point. Then repeat, using the last linked node as the centre of the next geographical region to analyse. When the end of a track was reached, then I'd pick the next "untraversed" point in the DB and start again. Then repeat this cycle until every point in the DB had been traversed. Initially the data is in a mySQL database table. I'd added an extra field, in which I'd converted the lat/long pairs to a geo "Point" field. Then used some rudimentary "MBRContains" queries to return points within a certain area (eg 100m2). I was then going to write some java code that explored/traversed the points, with repeated calls to get points within small areas. So my question(s) are: Does PostGIS offers any functions that could assist with the task of identifying related points ? Or are there any better approaches out there to achieve this ? Maybe there are tools/routines to do just this kind of thing, but I don't know the terminology enough to know about it yet :) Thanks for any advice.


title: Problem with taverna

text: I want to use taverna to call OGC WFS,and the file returned from WFS is GML, but what I want is SVG,so I have to convert GML to SVG by using taverna.The Method I used ExecuteProcess_gml2svg,which comes form a WSDL,.The process and the GML datasource is as shown ![enter image description here](http://i.stack.imgur.com/pJL3z.png) The problem is when I use the datasource url in the picture (that is the 'random_points' gml datasource is http://rsg.pml.ac.uk/wps/testdata/random_points.gml ),it works well,but when I download the datasource ,and use Appache server in my pc to publish the same gml datasource(the url is like http://chenhongyan.com/random_points.gml ),the process does not work,and the errors are: Processor 'ExecuteProcess_gml2svg_2' - Port 'ProcessOutputs': An error occurred invoking the WSDL service object of type 'exceptions.IOError' has no len() and I just don't know what can I do to solve it.Can you help me?


title: Record Set (Recordset) for Join in ModelBuilder

text: I am trying to use a recordset to develop a Table View to use to join to a feature layer in model builder. I can export the table view directly with the user input included. However, when I attempt to use this table view in a join i get no features in the resulting dataset. ![My Model](http://i.stack.imgur.com/KE7kk.png) FYI. The ultimate goal is to allow users to input multiple attributes from a feature layer (Using domain in the field of the recordset) and assign a unique value (i.e cost) for each of those attributes. I would then use the join to filter the table on those results and include the user added value in the resulting feature layer. So I am also open to another approach. Thanks!


title: Is there a way to force a featureclass to prevent the creation of non-simple shapes when using native geometry storage?

text: After spending a good chunk of the morning trying to figure out why a feature was displaying a circle when I was setting it as a rectangle polygon via SQL, I found the SDE_GEOMETRY table and some documentation about how ArcGIS stores curves in a side table. So even though I was setting the shape field correctly, ArcGIS was using the original geometry that was stored in the SDE_GEOMETRY side table (for those features that had been created with curves). Is there a recommended way to prevent a this side table from being used and have ArcMap disable any tools that would create the specialized geometry. My current workaround is to delete the row in an after insert/update trigger (seems to work but I haven't really tested it much).


title: Creating a new time column in PostGIS / CartoDB

text: I received a CSV file that has year (just two digits, '70'), month (as 'JAN'), day and time (as '0310') as separate columns. Bringing in to CartoDB / PostGIS they all appear as strings. I'd like to convert them to a single column that is a real 'date' type so I can display with Torque. What's the best way to go about doing that?


title: Projection problem? -- OpenLayers, Drupal, and MapBox

text: I've installed OpenLayers in Drupal 7 and created a view that pulls geofield data for a "Geography" vocabulary (states, tribes, and U.S. jurisdictions) and shows points on a map that I'm storing on MapBox's server. The map tiles I created in TileMill. The map displays fine, the points appear in the proper locations, but I can't get the map to center properly -- it appears somewhere off the coast of western Africa. I'm guessing that this is a "projection" problem, right?   * TileMill maps are always projected to "Web Mercator" (EPSG:3857)   * So, for the layer that I've added, I've used layer type "XYZ" and set the "Projection" to "EPSG:3857"   * In the "Layers & Styles" tab for the map, I've selected "EPSG:3857" for the "Map Projection"   * I've also set the "Display Projection" to "EPSG:3857"   * Because I'd like the map to be centered over/near California, I've made the "Centerpoint" in the "Center & Bounds" tab equal to "-122, 33"   * I've also tried "Zoom to Layer" in the "Behaviors" tab -- but that doesn't seem to work either. I'm guessing there's something conceptual that I'm missing here -- that there's something about projections that I don't understand, right? Or something else? I'd appreciate any help anyone can offer... I'm so confused!


title: FME remove polygon holes of under size ### sqft and chop the ends off

text: Quetion #1   I have a wokspace that takes the latest version of my data and creates buffers and then merges those to one polygon.   When 2 buffer polygons are near each other and close to create a hole. I would like to find a reasonable area size and remove only the holes. If there are small polygons for one reason or another I need to keep them.   ![remove small hole](http://i.stack.imgur.com/tuZtM.png)   Another question not exactly related mentions the donuthole extractor. Do I need to extract these and then delete them?   **EDIT:**   I used the solution with a tester and WAS successfull. The part I was missing was I had to put all of this "behind" my first disolver.   It wasn't an actual donut until the two parts were together. I suppose I could get fancy and put some other locgic in front to make it that way.   The other problem I was having was it is using sq meters not hectares or acres. :)   On to the second problem which I have not had any luck with yet.   ![image3](http://i.stack.imgur.com/gyRWm.png) ![workspace](http://i.stack.imgur.com/BlGYw.png) Question #2   I like the rounded version not the butt of bufferer. However would like to limit the ends (after the merge) to being only a specified distance from the end of the actual line.   ![bufferchop](http://i.stack.imgur.com/D3pEg.png) ![image6](http://i.stack.imgur.com/SbNHX.png)   OK I figured out I need to use the parameter created when I set the AreaCalculator. But I still don't see proper results. I would like to chop off the end flat at maybe 20' offset from the actual end of the line. Can I get a direction or pointer on accomplishing this?


title: Modelbuilder Parameters - cannot select input files when model dialogue box is open

text: I created a model in modelbuilder and set the inputs to be parameters. When I open the model and the tool dialogue box shows up, the input files that I wish to select are not shown. These inputs are used for the "Select layer by location". Why is that and how do I fix it? I have all the input data set up in a geodatabase, and I can see it when I navigate to the .gdb in ArcCatalog.


title: ArcGIS: Rapid digitizing of Point Data

text: I have several big maps with a lot of points that I have to digitize. All these Points have different Values (depth in Meter, double) and will be stored in 1 Shapefile. Is there a way to speed up my digitization process? e.g. Create Point with my Mouse in right Hand and directly Type the Value into a NumPad in my left Hand? Right now I always have to select the Item in the Attribute Table, where i can edit the Value. (Too much clicks) Thanks!


title: Collector cable optimizer

text: I'm hoping to find some pointers for a project that i'm supposed to be working on as part of a team. We're going to create an extension within ArcMap that uses python to process information within a shapefile that the user drags into a toolbox. The program is supposed to determine as final output the best routes for laying collector cables. The model will need to account for, among other things, road allowances, properties signed onto the project, the maximum number of Wind Turbines per electrical circuit, environmental features, cost for running x number of collector cable circuits overhead, cost for running x number of collector cable circuits underground etc It's a pain because i've got to learn a new language from scratch. Anyway, i understand that such models exist already and i was wondering if someone could shoot me off in the right direction by pointing one or more out and possibly any other useful resources as well. Thanks


title: QGis points distance from shore in metres

text: I have a grid of depth points in the Irish Sea that i'm trying to add other data fields to, in this case distance from shore. I've tried both of these suggestions to no avail, but am having some success with this one. My problem is that the resulting raster ![resulting raster](http://i.stack.imgur.com/tpNgk.png) has distances as decimal degrees rather than metres. I've changed my project properties > canvas units to metres, re-done the rasterise, re-done the proximity, but it's still decimal degrees. Any ideas how I can get the results in metres?


title: Null Values After "toreal" Command

text: I am very new so help is greatly appreciated! I have joined a simple CSV file to my project and am attempting to change a string value to a real value. In editing mode, I bring up the field calculator and create a new field from the first field. I try to change the string to real by doing the command "toreal(nameofmyfield)". When I create the new field though, I get NULL for every value, even though I see the normal numbers in the column I am trying to change. Any idea what's going on?


title: Calculate attributes to table based on other attributes

text: I am going through thousands and thousands of road lines and giving them a value based on QC with imagery. So for instance for the add_data field I give it one of these four numbers 0 = not examined 1 = include 2 = consider excluding 3 = exclude 4 = DC But for each of these I also need to fill out comments and type of road data. So I was wondering if there is a way to calculate the comment and other fields based on the add_value field? So maybe it might look like this for the add data: "1" = "include" comment = looks good RD_Typ = Major Road.   "2" = "include" comment = looks good RD_Typ = Dirt Road.   "3" = "include" comment = Joined data to other data RD_Typ = Major Road.   "4" = "include" comment = Joined data to other data = Dirt Road. ..... and so on you get the point. So the question is can I calculate the other fields based on the number I calculate in the add data field. I know I could do this with a join and then calc field but that just does not seem to be as much fun. Thanks for your help.


title: How do I limit update/insert/delete capabilities for a WMS layer in a QGIS project?

text: In QGIS 2.0 I have a layer based on a query to a PostgreSQL/PostGIS database and would like to serve this as a WFS layer. Users of the layer should be able to select subsets of the data based on properties such as dates or unique IDs, but should not be able to update, insert, or delete records. The Creating a WMS/WFS from a QGIS project section of the documentation simply says, (emphasis mine) > **In the WFS capabilities area you can** select the layers that you want to > provide as WFS, and **specify if they will allow the update, insert and > delete operations**. How, exactly, do I specify this? What is the correct way to deny, and to allow, these actions? * * * Here's what I've done so far. I added the layer I want to publish to a new QGIS project. Then in the Project Properties window, under OWS server, I checked the box in the Published column for the layer. Now I don't know what to put in the free-form text input boxes in the Update, Insert, and Delete columns. I've tried putting values of "0", "NO", and "FALSE" in here, as well as more facetious entries such as "Sure" and "Why not?". No matter what I enter though, when I click Apply or OK, QGIS gives no error or feedback of any sort. When I return to the project properties later, the text boxes are all blank as they were initially. Oh, and in case it makes a difference, I'm running this under Windows (currently 7, but it will be deployed on Server 2008 or 2012 eventually).


title: How to create Esri Flex Geometry from WKT

text: I am working with Esri flex 2.3 API, I am also using Abode BlazeDS for server side binding. I am sending a geometry as WKT from server into my flex client and I am looking for a way to convert the string into a flex geometry so that i can make a Graphic object and add to Feature Layer! or do I have to write a parser myself to accomplish this?


title: Symbology or Styles

text: Anyone know where to find me some more symbology or styles that I can use for my ArcMap 9.3.1 ? I have lots of mxd files on my computer and how do I apply all of my symbology or styles turn on ? Every time I open my mxd files I have to turn on each in the Style Reference dialog box because they are not enable check.


title: Looking for an annotated vectorial map of the world

text: I'm looking for an annotated vectorial (PS/SVG/...) map of the world (+countries). Wikipedia contains many vectorial maps (e.g: http://commons.wikimedia.org/wiki/Category:Blank_SVG_maps_of_the_world )but the paths are defined like that:               <path fill="#ccc" d="M98.171,199.174c-0.505-0.076-0.938-0.322-0...      ... and I'd like to find an **identifier** so I can colorize a given country (using XSLT or another method....). Something like that:               <path id="JAPAN" fill="#ccc" d="M98.171,199.174c-0.505-0.076-0.938-0.322-0...      thanks


title: Where can I find the various Esri tools to design a geodatabase?

text: There is a confusing multitude of tools written by Esri and others (some supported, some unsupported) for designing and diagramming a geodatabase. Here is an attempt at creating a definitive list of where to find them. Related to What are convenient uml tools to create geodatabase models?


title: Is it possible to create color tables for floating point raster maps in GDAL?

text: I've got a script that dumps a floating point raster that I'd like to automatically color. Currently I load this raster in Arc or QGIS and manually create a color table every time I re-run my model. I'd like my model to automatically create a nice color table, but the only functionality I can find is to create a color table for integer only rasters. My approach so far is to output 2 rasters, one with floating point values rendered in grayscale and a second integer raster that has the colors on it. I can flip between the two to see values, but it's far from optimal. Is there a technique to procedurally create floating point color tables in GDAL that automatically load with QGIS or Arc?


title: ArcGIS 10: AddIn Extension loading toolbox issue

text: I have an extension class in an ArcMap addin project of which loads a custom toolbox into ArcToolbox when the user opens ArcMap or opens a previously created map. My toolbox gets added to ArcToolbox fine except when the user creates a new map from File > New. My toolbox just disappears (no messages). I'm running AG10 SP2 but it doesn't work in AG10 base either. Anyone else experience this? Solutions?               protected override void OnStartUp()     {        AddToolboxes();        #region Wire events        ArcMap.Events.NewDocument += delegate() { ArcMap_NewDocument(); };        #endregion     }     void ArcMap_NewDocument()     {        AddToolboxes();     }     void AddToolboxes()         {                  // Load toolbox into ArcToolbox             IArcToolboxExtension arcTbxExt;             try             {                 arcTbxExt = (IArcToolboxExtension)ArcMap.Application.FindExtensionByName("ESRI ArcToolbox");                 if (arcTbxExt == null) { throw new Exception("ArcToolbox extension could not be found."); }             }             catch (Exception ex)             {                 System.Windows.Forms.MessageBox.Show("Error finding ArcToolbox extension." + Environment.NewLine + " - " + ex.Message);                 return;             }                  try             {                 IGPToolbox2 toolbox = Maersk.Functions.Toolbox.OpenToolbox("C:\\Toolboxes\\MyToolbox.tbx");                 arcTbxExt.ArcToolbox.AddToolbox((IGPToolbox)toolbox);                 arcTbxExt.ArcToolbox.Refresh();                  }             catch (Exception ex)             {                 System.Windows.Forms.MessageBox.Show("Error launching  toolbox." + Environment.NewLine + " - " + ex.Message);                 return;             }              }


title: Moving a Staging ArcGIS Server to the DMZ?

text: I am trying to move my EDN Arcgis Server with all the mapservices and caches to the DMZ .. The move is from our intranet to the DMZ to provide access to mapservices to outside world using the internet port 80 or 8080.. Everything is fine when we try accessing applications and mapservices from within our network to this box in the DMZ!! When it comes to external access it cannot access any ?? **I have tried a similar setup using windows xp box in the DMZ! Both internal and external machines can access and works well.. This should account for Network configuration correctness, hopefully..** _I am wondering if it has to do with ArcGIS Server on Windows 7 ?? Am I missing something wrt settings with ArcGIS Server or IIS ??_


title: Icon missing when adding custom combobox command to Toolbar

text: I have a custom command:               class ATL_NO_VTABLE CSliceSelector :         public CComObjectRootEx<CComSingleThreadModel>,         public CComCoClass<CSliceSelector, &CLSID_SliceSelector>,         public ISliceSelector,         public ICommand,         public IComboBox     {        ...     }      I also have a custom Toolbar where I try to add my custom command:               class ATL_NO_VTABLE CSeismicToolbar :         public CComObjectRootEx<CComSingleThreadModel>,         public CComCoClass<CSeismicToolbar, &CLSID_SeismicToolbar>,         public ISeismicToolbar,         public IToolBarDef     {        ...        STDMETHOD(GetItemInfo)(long pos, IItemDef * itemDef) {             if (itemDef == NULL)                 return E_POINTER;                  if (pos == 0) { // CSliceSelector                  itemDef->put_Group (VARIANT_FALSE);                 HRESULT hr = itemDef->put_ID(CComBSTR(L"{96A2B223-BF0D-410C-B22B-892FED3D9A15}"));                 if (FAILED (hr)) {                     return hr;                 }             }            }     }      However when I start ArcScene; I see a red circle icon and the text "Missing": ![enter image description here](http://i.stack.imgur.com/N7W7K.png), instead of my command :-(. However if I add it manually from Customize toolbars; it works fine. I have tried both returning a null pointer and a valid bitmap:               CSliceSelector()     {          m_hBitmap = ::LoadBitmap(_AtlBaseModule.GetModuleInstance(),MAKEINTRESOURCE(IDB_SLICE));     }     STDMETHODIMP CSliceSelector::get_Bitmap(OLE_HANDLE * Bitmap)     {          if (Bitmap == NULL)             return E_POINTER;        *Bitmap = (OLE_HANDLE) m_hBitmap;          return S_OK;     }      I have other commands that are not combobox; and they work fine.


title: Load partial shapefile into Postgis / GeoDjango project

text: I have a shapefile with Canadian postal codes, but am only looking to load a small subset of the data. I can load the entire data file and use SQL or Django queries to prune the data, but the load process takes about 2 hours on the slower machines I'm using. As the data I'm actually after is about 10% of the dataset, this isn't a very efficient process. I'm following the instructions in the Geodjango tutorial, specifically the following code:               from django.contrib.gis.utils import LayerMapping     from geoapp.models import TestGeo     mapping = {'name' : 'str', # The 'name' model field maps to the 'str' layer field.                    'poly' : 'POLYGON', # For geometry fields use OGC name.                    } # The mapping is a dictionary     lm = LayerMapping(TestGeo, 'test_poly.shp', mapping)     lm.save(verbose=True) # Save the layermap, imports the data.      Is there a way to only only import data with a particular name, as in the example above? I'm limited to the Linux / OS X command line, so wouldn't be able to utilize any GUI tools.


title: QGIS proximity raster with extent beyond original coverage?

text: I need a raster coverage of distances from points outside the geographic region of interest, and am having trouble doing this with QGIS. Questions like this one (Calculating distance to points in QGIS) seem close to what I want, but I'm getting hung up on the first step. The raster that I create with the QGIS "Rasterize" tool has an extent that does not overlap with the region of interest (because all of the points I want the distance from are outside of that region). I don't see any option for changing the extent of the generated raster. I am then unable to generate a Proximity raster with the extent I need. How can I get the correct extent? In addition, when I run "Proximity", I consistently get an error that my raster cells are not square and so distance calculations may not be accurate. I don't know how to make square cells because the tool asks for how many cells I want in the x and y directions, not the x and y dimensions of each cell. Is there a way of fixing this? Ultimately, I'm a bit uncomfortable with the necessity of converting point data to raster in order to calculate proximity, and wonder whether there is an open source tool that can do what I want without the "Rasterize" step...


title: Set symbols for a layer in ArcMap with two classes (greater than zero and less than zero)

text: I need to symbolize features in a feature class with two colors thus defining what features have positive (>0) and negative (<0) values. The only two ways I found this to be possible are: To use a static field with the values calculated beforehand OR Use definition queries using two identical layers in the TOC. Since it is more time consuming to update labels for two layers, I have approached this problem in another way by using Graduated colors symbol setting with two classes (Manual classification). I have added two extra features in the feature classes with really big and really small values (to cover possible values that can come later on while the feature class will be edited) and then edited the intervals in the Range columns in the Symbology tab of the Layer Properties. Afterwards, those two features are deleted. Now features' values will fall within these two ranges and will be symbolized properly. ![enter image description here](http://i.stack.imgur.com/qkBMb.png) Considering the need to update multiple documents with these settings and having different datasets as layers, are there any other approaches (without using Python to automate the hack described above)? All other GIS.SE posts I've found so far refer to the methods outlined above.


title: Arcmap point symbol/marker size, How it works

text: I am trying to figure out how the size of markers in arcgis works. I understand that a square marker set to 1 pt is 1 pt for specific scale. but when I import emf or font, How does the change size works then? if like the symbol is rectangular or a picture of a bus? if we take the bus, I set 20 pt size does it take 20 point of the smallest part of the picture or what? please help! :) and also, how can I set a specific size to be correct for illustrator export? the size gets different, is there any script for arcmap or illustrator to set correct size of symbols?


title: How do I format a text box in MapInfo

text: I want to import a table (created and formatted in MS word?) to appear as text in the MapInfo layout window. Obviously this isn't individual data records which appear in the MapInfo browser window. So would I need to format the MS Word table, save the file as a JPG/png format then import into MapInfo as a raster file? Is there another way in MapInfo of doing this - just type the table contents as text into the layout window then draw a rectangle around the text and make this box transparent but with a border?


title: How to load interlis1 files directly into QGIS 2.0.1 Dufour 64Bit (W7Pro)

text: with earlier version of QGIS (1.8) it was possible to load data in interlis 1 format directly. If I try now it says something like unknown format. Is this because I'm using the 64Bit version of the installer? A possible workaround could be using the examples at http://giswiki.hsr.ch/HowTo_OGR2OGR "interlis to shp". It would be very nice to be able to open Interlis Files directly again. Thank you for any hints. Regards, Adrian


title: Why does one part of GML -> KML conversion not work on Google Maps?

text: Brand new to all of this, not even sure this is the right forum to ask in... I have a client who wants to visualise cancer incidence and mortality rates by local/unitary authority across Great Britain. I've used Ordnance Survey data to get the boundaries of the authorities in GML format and saved all of them as one large GML file. I used QGIS to convert the data to KML format. I have the data in Fusion Tables now and it all looks great except for one blatant gap in the coverage. When I look at the KML data for the missing local authority, it hasn't imported at all (ie the cell is empty). I've tried re-uploading the geometry (and re-downloading and re-converting) and Fusion Tables is definitely rejecting it as KML. I've read somewhere that Fusion Tables can mess around with your geometry somewhat, so I uploaded the KML file and pointed Google Maps directly at it and that didn't work either (I don't know how good a test that is). I don't know that I originally created a valid GML file (I don't know the standard) but they seemed to import ok to QGIS and the final result looks great apart from this one place. The area _is_ visible in QGIS when I load the GML file, so I don't know if QGIS is producing invalid KML, or mayble Google Maps/Fusion Tables have limitations? **Edit:** I've downloaded the OS Open data pointed to by nhopton, loaded the appropriate layer into QGIS ( _Layer_ -> _Add vector layer_ ), saved it as KML ( _Layer_ -> _Save as_ ) and uploaded it to Fusion Tables. (Note that I didn't explicitly do anything with the CRS). Generally it's worked well (and importantly the particular area that had failed with my first method worked with this, so I have the missing geometry!). However, there's a number of missing areas again. And in this case, there's much less that _I've_ done to the data myself, so it's not so easy to assume I've just messed up the data. Any ideas why my final results look like this? I'd really appreciate some guidance from anyone who knows this area a bit better:) TIA


title: How to find the elevation differences between adjacent pixels in ArcGIS?

text: Is there a way to determine what the difference in elevation is between two adjacent pixels? I am using an SRTM DEM to map faults. I am using the vertical error (>16m) as a criterion for the difference in elevation to accurately map these faults. There must be a more efficient way to do this than by using the i-tool. Cheers


title: SQL expression for finding strings including letters in ArcGIS 9.3.1

text: I'm working on ArcGIS 9.3.1. I have a shapefile with a string field containing point names like these:               10     12     12/13     t/9     6577y     1234     12t9     124y67     u9      and so on (and yes, I know it looks like a colossal mess. It's not my design). Using `WHERE MyColumn LIKE '%Search String%'` I can find (using select by attribute) the ones containing a specific letter or a backslash (/). MY QUESTION IS THIS: How can I select all values containing a letter (in this example: t/9, 6577y, 12t9, 124y67, u9), not caring what letter it is or it's position? (is there a way to specify a position?)


title: QGIS converting to ASCII question

text: I am attempting to do some analysis on a raster file and I need to convert it to ASCII. In QGIS, I used [Raster -> Conversion -> Translate] with the asc file format, and I suspect it is simply copying the file without converting anything. Are there any other suggestions on how to do the conversion?


title: points on google layer using openlayers

text:
HI, I want to load some points, whose lat-long values are saved in a .txt file over google maps. The following is the code which I am using for the same. The google map layer is succesfully loaded if I keep lines 1-10 & 13-17. For loading the data points I am including 11-12 lines but I am not getting any output at all. I request you to please find the error and suggest how to do the same. "data.txt" contains the lat-long values with tab delimited present in my server.               1.     2.      3.<title>OpenLayers Example</title>      4.      5.      6.      7.

      8.      16.      17.       When I am replacing above 11-12 lines with               11.var pois = new OpenLayers.Layer.Text( "My Points",                      { location:"./data.txt",                        projection: map.displayProjection                      });      12.map.addLayer(pois);      Even then I am not able to get the output. I require some needful help from any one of you


Does anyone know if there are any sample questions available for the ESRI Web Application Developer Certification exam? Thanks, Seth _**_*** _UPDATE_ ** _*_** ESRI has released 3 sample questions for this test: http://downloads2.esri.com/campus/downloads/certification/samplequestions/EWDA10_Sample_Questions.pdf


title: WMS layer disabled when i grouped

text: i am using geoserver 2.4.0 for my GIS based web application.My Wms layer disabled when i grouped more than 3 layer in my application.Previously it worked fine.This problem occurs in real time only,in geoserver this problem does not exist.i don't know! any updation happened in geoserver? How can i fix it?kindly provide me solution for this problem?


title: Why is the geonames not loading?

text: I want to import a geonames file. I downloaded the geonames .txt file from http://download.geonames.org/export/dump/ I used: `v.in.geonames input=PH.txt output=ph` but it displayed this error: `ERROR: File 'C:UsersRhenDesktopPHPH.txt' not found` I don't know why it cannot be imported..


title: create temporary point geometry feature from coordinates in OGR with Python

text: I have two shapefiles containing two types of geometry features: one contains polygons and the other linestrings. Because I'd like to use the Within method to check whether the end(or start) of a linestring falls within a polygon located in the other shapefile, I need to access the first and last points (as geometry features). At the moment I iterate through the linestring feature and can identify the points using geometry.GetPoint(x) which returns the x and y coordinates of that point. But when I'd like to use within to check whether a given point falls in the polygon object I can't because as Paolo explained in this question I need the geometry reference (using GetGeometryRef())of a point which can be retrieved like:               ...     pointshp = driver.Open('test.shp')     pointlyr = pointshp.getLayer(0)     point = pointlyr.GetNextFeature()     point_geom = point.GetGeometryRef()     point_geom.Within(polygon_geom)      but I'm actually opening a shp that has linestring features. So after skimming through Utah State University GIS slides, I came up with this:               def within(pointcoord):     driver = ogr.GetDriverByName('test.shp')     ds = driver.CreateDataSource('test.shp')     layer = ds.CreateLayer('test',geom_type=ogr.wkbPoint)     fieldDefn = ogr.FieldDefn('id',ogr.OFTInteger)     layer.CreateFiled(fieldDefn)     featureDefn = layer.GetLayerDefn()     feature = ogr.feature(featureDefn)          point = ogr.Geometry(ogr.wkbPoint)     point.AddPoint(pointcoord[0],pointcoord[1])     feature.SetGeometry(point)     feature.SetField('id',1)     layer.CreateFeature(feature)     pointlyr = layer.GetLayer(0)     point = pointlyr.GetNextFeature()     point_geom = point.GetGeometryRef()     for x in xrange(len(polygons))             point_geom.Within(Network.polygons[x].GetGeometryRef())      Which looks really inefficient to me. So just for the sake of using a method(within) I need to create not only a layer, but a shapefile, when all I really need is a simple point feature from which I can get the geometry reference to compare it with the one of the polygons. The question is, how to create a point feature (temporary one, to be destroyed) from x and y coordinates? Or is there a simpler way of doing that than the example above? On another hand, is there any other way of getting the geometry reference of the start or end point (as a geometry feature) of a linestring feature?


title: How to select a raster file that is under a polygon in a vector file?

text: I am using: v.db.select map=PHL_adm1 where=NAME_1='Bulacan' My vector map is over the raster map. I want to select the part of the raster map under the polygon(NAME_1='Bulacan') of a vector. Is it possible to do?


title: calculate kernel density for polygons in ArcGIS

text: I have a population data set at the census block level. My goal is to use the data set to identify some population centers in the region based on population density. The simplest way to do it is to calculate the density by census blocks and use an arbitrary density cut-off value to determine those centers. One shortcoming of this method, however, is that it relies on census blocks, which have irregular size and shape. The other approach I am thinking of is to covert census blocks into points based on their centroids, and calculate kernel density based on these points. But I am not sure if this makes sense. In other words, does this approach help me overcome the problems that are associated with the first approach? Any advice is appreciated. Also, if it makes sense to do that, how the output raster size should be determined?


title: Water Quality Data Acquisition Sources

text: I was wondering what sources are available, or if anyone has water quality data for the Great Lakes and specifically Lake Ontario. We are interested in mapping Phosphorus, Organic Matter, Nitrogen and Dissolved Oxygen concentrations in Lake Ontario as they relate to algal blooms and hypoxia. Thank You


title: PostGIS: convert point data to raster

text: I have a table of points and I'm trying to visualize their distributions. Each point has a POINT geometry (WGS84) and some numeric attributes: a, b, c.               | id | geom       | a  | b | c |     |-----------------|----|---|---|     | 1  | POINT(...  | 12 | 9 | 2 |     | etc ...                      |      The first thing I tried was to project to an equal area projection (EPSG: 3175), use `ST_SnapToGrid` and then group on the resulting point geometry.               SELECT        avg(a), avg(b), avg(c),        ST_Transform(ST_SnapToGrid(ST_Transform(geom, 3175), 50), 4326) AS bin     FROM points     GROUP BY bin      This groups the points into 50x50 meter bins with their aggregate values which I can import into qgis to play with. I optimized this by adding another column to the points table called `bin_id` which was a foreign key to a `bins` table:               | id | geom      |     |----|-----------|     | 1  | POINT(... |      Then I could write queries like this:               SELECT avg_a, avg_b, avg_c, geom     FROM (       SELECT          avg(a) AS avg_a,         avg(b) AS avg_b,         avg(c) AS avg_c,         bin_id       FROM points       GROUP BY bin_id     ) AS tmp     JOIN bins ON tmp.cell_id = bins.id      This all works fine but i want to use Rasters instead. All the examples I've seen show how to import existing rasters or shapefiles into PostGIS. But I need to know how to generate one from data already in PostGIS. Ideally I want a raster with 3 bands: `avg_a`, `avg_b` & `avg_c` Can someone point me in the right direction? EDIT 1: I'm not opposed to doing it outside of the database.   EDIT 2: Seems like I can use the GDAL API to do it.


title: How to set custom break points based on formula?

text: I work on projects that each require the creation of about 40 maps. In all of these maps, I use Graduated Color symbology to represent concentrations of a different population in area Census tracts. The break points I use are not equal interval--basically the first break point is the mean, the second is the mean+10 (the "disproportionate share threshold"), and the third and fourth break points are equal intervals between the disproportionate share threshold and the maximum value in the dataset. The maximum value is the last break point. I have a workaround that I'm using. Well, I've got the "around" part--the "work" part is not quite there: I use SPSS to collapse population per tract into either the minimum value, the mean, the disproportionate share threshold, then the two break points between that threshold and the maximum value, and the maximum value itself. Once I've done that, I join the table to a Census tract shapefile, and open up symbology. When I set the classification to "Natural Jenks", the break points are ALMOST exactly what I want, except that the minimum value is in a class by itself, and one of the upper break points are omitted. So, I have to adjust it by hand anyway. Is there some way to have the break points "skip" minimum value and snap to all of the values above that? Ideally, each of the five breaks should include a range between:   1. the lowest value and the mean,   2. the mean and the disproportionate share threshold (DSP),   3. the disproportionate share threshold and ((maxvalue-DSP)/3+DSP,   4. that last one and (2((maxvalue-DSP)/3)+DSP,   5. that last one and the maximum value. Thanks in advance for your help--I hope this was clear. I'm pretty new to ArcGIS, and just getting to the point where I "know what I don't know".


title: how to create a map with oneway flows using FlowMapper for QGIS

text: my plan is to create a map that shows the top 50 (by amount of passengers) international flights departing the United States (e.g. New York -> London, Honolulu -> Tokyo...). I want to use the plugin FlowMapper for QGIS to display those flights and have the width of the flows correspond with the amount of passengers. I have read the tutorial and had a look at the examples given (all in the download file of FlowMapper). The problem is, in those examples FlowMapper is only used for back and forth flows. Thus, the data matrix in the examples looks like this, each row and column have the same: ![enter image description here](http://i.stack.imgur.com/NxBTQ.jpg) But I want to use one way flows, so row and column have different geographical locations, like this: ![enter image description here](http://i.stack.imgur.com/khsth.jpg) So my question is, how does the flow matrix- and coordinates-file have to look like in order to be correctly displayed by the FlowMapper plugin? This is my complete speadsheet that I wanna use as basis for my map: http://www.sendspace.com/file/ckaejd


title: ArcPy Cursor iteration update

text: I'm trying to get the cursor to iterate through rows and fields looking for blank or null values and change them to `999999` but it's not changing anything.               fc = "C:\Users\\bbrock\Documents\ArcGIS\Ports.shp"          # Create a search cursor      #     rows = arcpy.SearchCursor(fc)           # Create a list of string fields     fields = arcpy.ListFields(fc, "", "String")          for row in rows:         for field in fields:             if field.type != "Geometry":                 if row.getValue(field.name) == '':                     row.setValue(field.name, '999999')                                   print "%s: Value = %s" % (field.name, row.getValue(field.name))                      if row.isNull(field.name):                     row.setValue(field.name, '999999')                                   print "%s: Value = %s" % (field.name, row.getValue(field.name))


title: QGIS/MATLAB add geo-coordinates flag

text: I've created some polygons in QGIS with which I'd like to create some thematic maps. Unfortunately, if I want to import the shapefile in MATLAB with: shaperead('nameoffile, 'UseGeoCoords') i get an error which states: Expected geo-coordinates flag (T/F) to follow parameter name 'UseGeoCoords'. Below you'll find a screenshot from MATLAB: ![MATLAB Screenshot](http://imgur.com/uNQCTKu.png) As far as I see there should be an option in QGIS, with which I can say that programs can use GeoCoords and not only X,Y coordinates. Has anyone an answer for this specific problem? Thanks and best regards!


title: Tools for planning a building

text: I'm planning to build a cabin, and I'd like to combine all the GIS data available from my home country with my own measurements for placing and the cabin layout. That is, I need a tool which can:   * save data in layers,   * import third-party GIS data for elevation, infrastructure, etc., and   * import a PDF, SVG or other vector graphics format to get the cabin layout. Is there any software useful for this available for Linux (at a pinch Windows XP), preferably for free?


title: Publish MXD to server in ArcObjects C#

text: I have an MXD that I would like to publish to an ArcGIS server using C# ArcObjects (preferably as an MSD). I have been unable to find anything regarding this. Essentially I would like to implement the following ArcPy code in ArcObjects:               arcpy.mapping.ConvertToMSD()     arcpy.mapping.PublishMSDToServer()      Any help would be appreciated! Thanks!


title: Table analog of MakeQueryLayer

text: As of 10.1, `arcpy` provides the rather handy `MakeQueryLayer` function. This allows us to use `arcpy` to create query layers programmatically instead of having to go through MXDs (as was the case in 10.0). This is extremely useful when working with a database that does not contain the SDE structure (formerly created via post-install and now created by Enable Enterprise Geodatabase, I believe) or where the tables you want to interact with have not been registered with SDE. However, now I want to create a _table_ instead of a feature layer/class, as in the query will not contain a shape column. How can I go about this? ## What I've tried: I considered `MakeQueryTable`, but this takes an existing layer or table as an input, rather than a query. I also just tried to use `MakeQueryLayer` and got some... confusing results. `conn` below is a `str` that contains the file path of my database connection file.               >>> r = arcpy.MakeQueryLayer_management(conn, 'TEST', 'SELECT ID AS OBJECTID, MY_TABLE.* FROM MY_SCHEMA.MY_TABLE', oid_fields=['OBJECTID'])     >>> r.status     4     >>> r.outputCount     1     >>> r[0]     u'TEST'     >>> type(r[0])     <type 'unicode'>      The status code `4` corresponds to Succeeded, but the result object contains a Unicode string as its only output. So I'm not sure if it successfully created a layer or table and I just can't find it or what. Normally, `MakeQueryLayer` outputs an actual `Layer` object:               >>> r2 = arcpy.MakeQueryLayer_management(conn, 'TEST2', 'SELECT ID AS OBJECTID, MY_TABLE_WITH_SHAPE.* FROM MY_SCHEMA.MY_TABLE_WITH_SHAPE', oid_fields=['OBJECTID'])     >>> r2.outputCount     1     >>> r2[0]     <map layer u'TEST2'>     >>> type(r2[0])     <class 'arcpy._mapping.Layer'>      I tried to do some work with `arcpy.mapping.ListLayers` and `arcpy.mapping.ListTableViews` to see if I could find the output layer, but those require an input document (and I'm not using a document). My back end database happens to be Oracle, but I imagine that shouldn't matter.


title: Map-based online questionnair/survey

text: thank you for your advice and I made changes to my question; I am planning to create map-based online questionnair/survey for my master dissertation on bicyle safety. What I want to do is, asking my potential respondants to mark their answers by adding pins/mark on Online map I create. For example, I ask them a question like, " where did you have a near-miss incident while you are on a bicycle?", then respondant can answer that question by locating pins/mark on where they had nearmiss. My initial investigaton suggested that I can achieve this with Google map API, ArcGIS online or Crowdmap. I have seen some crowd-sourced maps with a lot of pins located by many people. But I am not sure if I can;   1. download all these marked locations in csv/text file with xy coordinate   2. make my survey map accessable only for my "group" Can anybody suggest me what can I use (Google API, crowdmap, ArcGIS online or else) to achieve above 1 and 2 without any or with little programming? Thank you.


title: Positional error of line features

text: Does anyone know how to model the error in the spatial location of a line? I think Survey Analyst extension in ArcGIS had a possibility to construct ellipse errors but it is not available any more. Is there any extension for this or someone can describe any other way how this could be done?


title: Zoom to Selected Feature for point layers in QGIS 2.0

text: I find the QuickFinder plugin for QGIS 2.0 very useful to quickly search, find, and select data. However, I'm under the impression that checking the "scale" option doesn't seem to do anything. I thought that perhaps it would change the scale of the map and zoom it to a certain level, but my results are the same whether or not this option is checked. It seems as though Zoom to Selected Feature doesn't change the scale of the map either. I'm wondering if this might be why the scale feature for the Quickfinder doesn't do so. I'm trying to confirm if these are the general effects of these tools or if it is simply my implementation. After finding this very old ticket: https://trac.osgeo.org/qgis/ticket/154 I still feel confused about the current status of this topic. Does anyone have any ideas of what I may be missing? Thanks in advance.


title: What graph traversal algorithm should I use?

text: I would like to write an algorithm which can traverse a graph and hopefully later I can implement it to use for an indoor navigation system. The graph would come from floor plans of a building and the graph nodes represent the building objects such as doors, corridors, stairs, rooms, etc. The navigation system will accept different users with different capabilities, these capabilities are stored in user profile. For each individual user (or a group of users of the same type, e.g. staff member) the system returns a path specific to that user/s considering his personal functionality. For example if the user of the system is using a wheelchair the system should exclude all the stair nodes from the path (they will be marked non-traversable). Similarly if in order to enter into an indoor space, the user needs to have a key or be an staff member but he is NOT, the system should mark the node (e.g. the door which requires the key) non-traversable and should search for those paths which don't require these specifics requirements. This was the first part and here I use a function to check these conditions for a given user. The inputs to this function are a node and a user. If the conditions that are required by the node are NOT fulfilled (means user profile doesn't match the condition required by the node) the function returns false= non-traversable. If the condition holds true, the node stays isTraversable and will be stored somewhere and will be considered for the final path calculation. The second part is that I have a list that contains all the nodes that are traversable and now the shortest path will be calculate based on 3 different costs, Energy, time and money. The final cost can be a combination of these costs and it sets the weights for the edges. My problem is with the first part. I want to know how can I make an algorithm that produces the subgraph (the list with the traversable nodes in it). And later how can I use this algorithm to modify Dijkstra to fit this system. Can anybody help with writing a pseudocode for the algorithm so that I can understand how would it work in clear steps or any hints on how can I make these steps into pseudocode? Please note that at this stage the performance is not very important for me. Thanks in advance!


title: Python Def Query Help for multiple layers

text: I have Three layers in a .mxd, all three layers have a common value I can use, but a different Label in each Layer. As an example, the Value I want to use is "100" In Layer A, the Value is in [Name]   In Layer B, the Value is in [Series]   In Layer C, the Value is in [SeriesName] The Def Query for Layer A is "Name = '100'"   The Def Query for Layer B is "Series = '100'"   The Def Query for Layer C is "SeriesName LIKE '100-%'" What I am looking for is a python script that will allow me to enter in the number, and apply that number to all three Def Queries. Appreciate any an all help!


title: the difference between pattern recognition and image classification in remote sensing?

text: Does anyone know the difference between pattern recognition and image classification?


title: Setting zoom level

text: Is there any way I can set the zoom level between lets say zoom level 5 and 6? I can' select 5.5.                map = new Map( "map" , {                     basemap: "gray",             center: [-85.416, 49.781],                     zoom : 5,                     logo: false             });


title: Access atributte processing algorithm

text: I have a doubt. In my plugin I use algorithms from processing sextante. For example:               Processing.runAlgorithm("grass:v.to.rast.attribute", None, inputLayer, 0, Indexes, extent, cellSize, -1.0, 0.0001, outPath)      inputlayer = directory of inputlayer extent = defined extent of the new raster file cellsize = option for user to choose the cellsize outPath = output raster directory and other parameters that doesn't matter here but my question is in the attribute field (that is defined as Indexes in the example). This parameter is the attribute field to rasterize. So, usually, I used a combobox to define this parameter but in this case I create, internally, a new column with values and I want to use this column. But I don't know how this parameter are defined, as a list, as a stringm...? Does anyone know that? Or how I can access the attribute created by me?


title: Select distinct values from a single column of an attribute table (or layer)

text: Is there a way to select distinct values from a column in ArcMap? I have the data in both GDB and SHP formats. I have searched for ways to select using SQL, QueryLayers, ModelBuilder and individual toolboxes and it appears as all selection options always SELECT * FROM tableName WHERE ... In SQL I would write SELECT DISTINCT columnName FROM tableName. Thanks for your help.


title: ArcMap 10.1 SP Spatial Join Problem (lines within polygon)

text: I have a line layer (faults) and a polygon layer (fishnet grid). The fishnet has a risk of earthquake attributed to each polygon. I want to take that risk value and give it to the fault lines which pass through that polygon. I am attempting to use the spatial join tool to do this. I have the target layer as the fault lines and join layer as the risk polygons. I have selected the 'Within' operator and one-to-one join operation and set the risk values to be added as a field to the new layer with the merge rule set to sum (I want to add the risk of all the polygons that the line passes through). The tool works for lines which are completely within a polygon and returns the risk value correctly but any lines which pass through 2 or more polygons are returned with a null value instead of the summed risk. If I uncheck the keep all target features then it leaves them out, implying that the tool doesn't think the lines are 'within' any polygons. Is this an error in Arcmap? and how can I fix this? Or is that not how the 'within' operator works for lines within polygons? My assumption is that the 'Completely_within' operator should return what I am getting from the 'within' operator (also this is how it's described in the tool help). Green yellow and red lines have had the risk values correctly attributed but as you can see its only the lines which fall completely within a polygon. The Purple liens are the ones which aren't being computed with the summed risk but are returned with null. ![The output](http://i.imgur.com/vU6OlYZ.png) A screen cap of the tool before running (the warning is just the because the file name exists since i took it afterwards...) ![tool](http://i.imgur.com/IjrZhOM.png)


title: Can shaded relief be "blended" with DEM colour ramp in ArcGIS to produce a more visually dramatic result?

text: Not being able to use true transparencies with ArcGIS colour ramps where only solid colours are allowed has me always wondering about how to improve the look of the surface in my maps without resorting to time-consuming post processing techniques. (For example the basic graduated black colour ramp typically used with hillshade does not have levels of transparent black but rather shades of grey and solid white where 100% should be. A as a result, it is not possible to overlay a hillshade without impacting the colours of the features under the hillshade; overlaid colours will appear washed out) In this discussion @matt wilkie provides a link to GDAL's gdaldem for producing a hillshade without Spatial Analyst using. He also includes a link to Terrain Data Tutorial which shows how to how to use gdaldem to **build a "composite hillshade with hypsometric tinting".** **Can ArcGIS (Spatial Analyst, 3D Analyst) be used to produce a similar "blended" product?**


title: Why is ArcMap deleting data from a field when I add a new one?

text: My copy of ArcMap 10 is acting pretty strangely. I have a shapefile and all I've done is added a new field as a FLOAT and for some reason all data in another one of my fields is being deleted! Does anyone know why this is happening? Thanks, Alan! ![enter image description here](http://i.stack.imgur.com/Mxm8y.png) ![enter image description here](http://i.stack.imgur.com/7Mjvk.png)


title: Merge overlapping GeoJSON GeometryCollections

text: I need to merge three overlapping GeoJSON GeometryCollections of Polygon. Ultimately, the end result needs to be a geom type for PostGIS (e.g. MULTIPOLYGON) so I can easily check if particular points fall within it (e.g. with ST_Intersects). The GeoJSON files come from OpenTripPlanner, and they represent isochrones. The reason that there are multiple, overlapping files is that we use the same starting point but more than one starting time. So actually they are _very_ similar but not exactly. We tried mapshaper for merging but it does not really "union" the polygons, instead it puts them all in one geometry so effectively that just puts them all on top of each other. We were able to convert each file individually into multipolygon using ST_CollectionHomogenize however still missing what function or tool we can use to merge these into a single multipolygon. Note that preservation of outer rings is a requirement, which seemed to rule out any standard PostGIS function we tried.


title: QGIS 1.8.0-Lisboa error message when converting DEM to Contours

text: I am trying to convert a small DEM (1 megabyte) into contours. When I go to - Extraction - Contour- and run it I get the following message: "ERROR 1: CPLRealloc(): Out of memory allocating 416 bytes." It doesn't seem to matter what the contour interval is, although I will get a slightly different error message (for example 40 bytes). And, I have successfully converted a much larger DEM in the past. Now I get the same error message when I attempt to convert that one to contours as well. Computer doesn't have the most power, but it worked in the past. Windows 7 Professional, 2 gigabytes of RAM, I have tried different contour intervals, checked the projection, made sure no other applications were running, and can't think of what else might be going wrong.


title: shp2kml license expired, looking for free solution

text: So I work for a small African non-profit. I do all the GIS and data work. Most of the layers I create, I convert to kml layers so my boss can then use them quickly and easily in Google Earth. But now in the last months my shp2kml program says its license is expired and to download the newer version. I've downloaded the most up to date version I can find which isn't that new Feb 2011 and I'm still getting the same error. Anyone know another freeware program that does the same thing or know a way around this or to get a more up to date version of shp2kml. My operating system is windows 7 and i'm using QGIS 1.8.0. The non-profit can't afford Google Earth Pro edition. I need free or practically free solutions. Thanks for your help.


title: ESRI's WmsLayer failure

text: I've been using the WmsLayer for my WMS Service needs, rather than my old way of doing things by sussing out a GetMap URL and posting that image to a GraphicsLayer. It's been working well but I seem to have found a service that breaks it. Here is the getCapabilities URL I'm using:               http://ags.pdc.org/arcgis/services/global/pdc_active_hazards/MapServer/WMSServer?SERVICE=WMS&REQUEST=GetCapabilities      I'm able to to get the Name of the layers available (they're using simple numbers to ID them). I provide the WmsLayer the getCapabilitiesURL, the layer ID number & a Proxy Url based upon the proxy handler from the documentation here at ESRI. Unfortunately that doesn't seem to be enough. When debugging, I crash to an unhandled exception as soon as I add it to my map control.               Category: ManagedRuntimeError     Message: System.InvalidOperationException: Sequence contains no elements      Using Firebug to look at my hosted application, it reports this error:               Unhandled Error in Silverlight Application Data at the root level is invalid. Line 1, position 1.   at System.Xml.XmlTextReaderImpl.Throw(Exception e)        at System.Xml.XmlTextReaderImpl.Throw(String res, String arg)        at System.Xml.XmlTextReaderImpl.ParseRootLevelWhitespace()        at System.Xml.XmlTextReaderImpl.ParseDocumentContent()        at System.Xml.XmlTextReaderImpl.Read()        at System.Xml.Linq.XDocument.Load(XmlReader reader, LoadOptions options)        at System.Xml.Linq.XDocument.Parse(String text, LoadOptions options)        at ESRI.ArcGIS.Client.Toolkit.DataSources.WmsLayer.client_DownloadStringCompleted(Object sender, DownloadStringCompletedEventArgs e)        at System.Net.WebClient.OnDownloadStringCompleted(DownloadStringCompletedEventArgs e)        at System.Net.WebClient.DownloadStringOperationCompleted(Object arg)      The XML returned by the GetCapabilities call seems sound to me - at least it appears like all of the other WMS services I've been communicating with - so I'm not sure what parsing errors it is running into. Any ideas? I can provide further info if there is something I haven't thought to check.


title: Open multiple hyperlinks in one point directly from QGIS?

text: I have an enormous shapefile and each point has about twenty different attributes. About five of these twenty are hyperlinks- is there anyway to open these directly from selecting a point in QGIS? For example, in ArcGIS Explorer, you can click on a point, look at the point's attributes, click on the hyperlink, and it will open adobe/autocad/etc. I installed hot links plugin but couldn't figure out how to use it, let alone any information about it online. Any advice would be much appreciated! Thanks.


title: Using python to Raster clip the extent to be the same as the display?

text: I've been trying to figure this out for almost a week now. I am creating a tool which first needs the raster to be clipped out to the current display. I have tried all sorts of things. Below is the current code I've been trying to get working. You can also see the previous code commented out.               import arcpy     view = arcpy.mapping.MapDocument.activeView     #view = arcpy.env.extent          # Local variables:     inputRaster = "D:\bt_gis\Contours\InstaContours\LA_OC_DEMs.img"     outputRaster = "C:\\Users\\JZimmer\\Documents\\ArcGIS\\Default.gdb\\LA_OC_DEMs_Clip2.img"          # Process: Clip     arcpy.Clip_management(inputRaster, view, outputRaster, "", "-3.402823e+038", "NONE")     #arcpy.MakeRasterLayer_management (inputRaster, outputRaster, "", view, "")      Any input is appreciated.


title: Styles in WMTS layers

text: Served through GeoServer 2.4.1. This behavior has been observed trying to view the layer both through OpenLayers and in QGIS. I can display multiple styles on the same layer using WMS, but only get the default style in WMTS. If I specify a style that does not exist, the layer won't display at all, but if I give it a valid style it returns the tiles for the default. Any ideas? My suspicion is that either styles are a feature not yet implemented in Geoserver for WMTS, or that there's some setting that I'm missing on the server side since I get the same behavior on two clients.


title: Is it possible to use the ArcObjects Runtime SDK on a Windows Store Application (windows 8, 8.1)

text: I'm looking to develop a simple windows 8 tablet application with arcobjects. Is it possible to use the ArcObjects Runtime SDK? If not, is there an ArcObjects SDK for Windows store Applications?