title: computeDistance() does not work in Google Maps API

text:
I am trying to learn Google Maps API. When I try the below code, execution stops at `var km = computeDistance(position.coords, ourCoords);`. I get to see the `"before compute distance"` alert but not `"after compute distance"` alert.               <!doctype html>                           <meta charset="utf-8">             <title>Where am I?</title>
                 Your location will go here

                 Distance

                    How would I fix this?


title: Can I just swap out a .prj file to reproject a shape file?

text: I'm an experienced software developer -- but a GIS noob. I have a shape file that I am trying to turn into a SVG map with Kartograph. When I try to run Kartograph I get this error:               ERROR 6: No translation for Lambert_Conformal_Conic to PROJ.4 format is known.      Because I don't want to delve into Kartograph, it seems like I need to reproject the shape file. The wikipedia page for shape files says that the .prj file describes the projection.               .prj — projection format; the coordinate system and projection information, a plain text file describing the projection using well-known text format      Can I just swap out the prj file to reproject the file? If so, where can I find a .prj file that will work to reproject the map so that Kartograph can turn it into an SVG?


title: Can input to ModelBuilder be selected features in ArcMap?

text: I want to find out if the manually selected features in ArcMap can be the input to a ModelBuilder tool. If not, is it possible to copy or export selected features in ArcMap using ModelBuilder? Once in ModelBuilder, the data will be manipulated and exported to a dwg file but I first need to figure out how to get the selected features in ArcMap as the input to ModelBuilder.


title: Filter in Qgis Web Client url

text: I don't know if this is possible but I read examples in qgis web client documentation to filter layer data like FILTER=Country:...... Can i do something like:               mylocalserver/QGIS-Web-Client/site/qgiswebclient.html?map=/usr/lib/cgi-bin/argentina/argentina.qgs&FILTER=localidades:LOCALIDAD='AGUARAY'      localidades is the layer name in my qgis desktop project. When I do a "filter" en the qgis desktop application with the "Query" panel, I have the information just like I want. Sorry for my English and thanks in advance!


title: Drawing polygons with Google Maps using local shape data

text: I am trying to draw a census tract on a Google map. It looks like Google Maps lets you draw polygons using longitude and latitude coordinates. However, the data that I have is not longitudes/latitudes. The specific dataset I have looks like this. The "spatial reference" for this data is 3857. But after looking up this spatial reference on this page, I still am not quite sure how to go about drawing the shape to the map using longitudes/latitudes. I assume that this sort of thing can be done with GIS software or command line tools, but I am more interested in code libraries or APIs that could make this possible. Also, if anyone knows of a better source of census tract shape data which is more compatible with Google Maps, that would be great too.


title: NoData specification via IEEE NaN in GeoTIFF for ArcGIS

text: I am a programmer exporting 32-bit floating point GeoTIFF rasters from my own application that my users need to read into ArcGIS. Unfortunately I am not an Arc user myself and have _zero_ expertise with it, so please target comments accordingly. Arc is reading and locating my rasters properly, but I am having problems making it understand which grid cells are empty. Since all positive and negative floating point numbers are legal data values in my rasters, I am trying to use IEEE NaNs as my NoData values, which I understood to be an accepted practice. Furthermore, as my rasters are often dimensionally large but sparsely populated, I'm trying to get maximal compression out of the TIFF PackBits algorithm (which really likes to see sequences of identical bytes) by making an optimal choice with regard to my particular NaN value. Since any 32-bit float whose exponent bits are all 1 is considered a NaN, I'm using 0xffffffff (32 1-bits) as my NaN value since it's both a valid IEEE NaN and is also 4 identical bytes. Arc seems to treat my 0xffffffff NaN values as low but legal values, i.e., it renders them as black similarly to the actual low legal data values in my file. How can I get Arc to treat this particular constant (which is of course _not_ a value that can be expressed as a floating point number since it's an actual IEEE NaN) as my NoData value and not render it at all? Not being an Arc user it's quite possible that I'm just not setting the right display/import options within Arc to get what I want. If someone could explain how to do this that may well solve my problem.


title: How to load multiple images with the same CRS?

text: When I add several images to QGIS 1.8 it asks the coordinate reference system for each one. There are several seconds in between the dialog and it is not easy to just continue working as the dialog keeps popping up with focus. Is there any way to check "the same" CRS for all selected (inserted) rasters?   With 10 it is not that big of a deal but when it is 100 it seems to occupy a lot of time. I normally use the gui so writing a batch program is not what I am looking for.   I select NOT all of the images in a given folder.


title: Using onRectangle for user input in ArcMap/ESRI addin

text: In an effort to try to add user input to my project, I'm developing an esri python addin to generate rectangles to my map. I'm on step one now, just exploring the addin functionality. I tried making an addin 'tool' precisely per the ESRI sample code, but the button isn't working. It appears on my ArcMap desktop, but when I click it, nothing happens. I just really want to know what's happening with the onRectangle, and how I can get it to accept user input. Here's the sample code that I essentially stole from ESRI:               import arcpy     import pythonaddins          class fishnetsExample(object):         def __init__(self):             self.enabled = True             self.cursor = 3             self.shape = 'Rectangle'              def onRectangle(self, rectangle_geometry):             """Occurs when the rectangle is drawn and the mouse button is released.             The rectangle is a extent object."""             extent = rectangle_geometry             # Create a fishnet with 10 rows and 10 columns.             if arcpy.Exists(r'in_memory\fishnet'):                 arcpy.Delete_management(r'in_memory\fishnet')             fishnet = arcpy.CreateFishnet_management(r'in_memory\fishnet',                                 '%f %f' %(extent.XMin, extent.YMin),                                 '%f %f' %(extent.XMin, extent.YMax),                                 0, 0, 10, 10,                                 '%f %f' %(extent.XMax, extent.YMax),'NO_LABELS',                                 '%f %f %f %f' %(extent.XMin, extent.YMin, extent.XMax, extent.YMax), 'POLYGON')             arcpy.RefreshActiveView()             return fishnet


title: Arcpy ERROR 000732 using tools from Toolbox (Copy, FeatureToFeature etc.)

text: I have got a lot of feature classes in my sde postgreSQL database. When I check the existence of sde database using arcpy.Exists(inGDB) I get the True result. In the case when I check the existence of one feature class using the same function I get False. Any attempt to copy feature class from database is ended by ERROR 000732 (Input Data: Dataset does not exist or is not supported) This situation is weird because when I copy code with function into empty script file and perform checking with the same parameters I get True result for feature class. I do not understand why everything is fine in standalone script but in the case more complex scripts it does not work. Any idea? EDIT: I think I have found the reason of my errors but it is very strange behavior of Python. At the beginning of complex script I import a few my other scripts as modules. Inside of each script I perform import of arcpy. During different attempts to solve problem I have noticed that when I commented import in my other modules in the main script everything is ok, when I have commented out (import arcpy) again I get errors (ERROR 000732) in the main script. I do not know how to solve this problem, because without imports in the other modules they do not work (I get error “global name 'arcpy' is not defined). SOLUTION: The error was caused by lack of path in environment variables. Because my scripts should work at ArcGIS Server I use different version of GDAL. One of my scripts that imported as module in main script caused my error. Initially in my code I use only path to GDAL catalog. I solved that by this code:               if 'ArcGIS10.2' in str(os.path):         path = r'C:\Program Files (x86)\GDAL         path1 = r'C:\Program Files (x86)\ArcGIS\Desktop10.2\arcpy'         path2 = r'C:\Program Files (x86)\ArcGIS\Desktop10.2\ArcToolbox\Scripts'         path3 = r'C:\Program Files (x86)\ArcGIS\Desktop10.2\bin'     if 'ArcGISx6410.1' in str(os.path):         path = r'C:\Program Files\GDAL'         path1 = r'C:\Program Files\ArcGIS\Server\arcpy'         path2 = r'C:\Program Files\ArcGIS\Server\ArcToolbox\Scripts'         path3 = r'C:\Program Files\ArcGIS\Server\bin'     os.environ["PATH"] += os.pathsep + path     os.environ["PATH"] += os.pathsep + path1     os.environ["PATH"] += os.pathsep + path2     os.environ["PATH"] += os.pathsep + path3


title: organize the map theme and style configuration file

text: I am working on an android application which render some features directly, but these different features should be rendered differently at different zoom level or only be visible from a certain zoom level. So I think I need a config file to keep this information. During the render I will read the file, and get the matched style for a given feature, and the render it. Since all of these are done at the client side, so I think the file size should be as less a possible(consume less memory), and the speed of querying for a matched style should be as quick as possible. I have no experience about organize this kind of file, so I wonder if someone can give me some suggestion? BTW,all the data no matter points or lines have an property named "uid" to identify them from others.


title: How to extend OSM2PO with ferry durations?

text: I would like to have the ferry duration [1] taken into account, when transforming OSM files to SQL. But I'm not sure if I need to write my own waytagresolver, or it is possible to simply add it to the osm_meta tag by adding a parameter or by config? I'm hoping for the latter, as this of cause is the simples by far :) [1] http://wiki.openstreetmap.org/wiki/Tag:route%3Dferry


title: Can QGIS Sextante or ArcGIS ModelBuilder batch extract vertices from KML polygons and calculate coordinates

text: I want to know how to extract the vertices of any KML polygon and also calculate the co-ordinates of such vertices in Degree Decimal format. I am looking at a batch process. Can we use the QGIS Sextante or ArcGIS ModelBuilder to do such a thing?


title: Error when converting .shp to .dxf

text: I have created a polygon shape file in Erdas Imagine 9.2 from a raster image. I want to export it to .dxf or .dgn for further processing. When I export the .shp to .dxf in Global Mapper, large polygons sharing the border of an image are missing. I tried in Erdas Imagine but the result was in Line String. Would appreciate any help. Thank you


title: What DEMs does QGIS support? + DEM data sample

text: I am working on QGIS rasters. I want to know about DEM file formats supported by QGIS. I know that QGIS can read the GDAL supported formats. but what is GDAL supported format ?   Does QGIS support the DEMs supported by ArcGIS ?   I want to have some samples of raster and DEM files for QGIS.   Thank you for your help.


title: Point Layer WMS GetFeatureInfo failing when zoom in

text: I have a problem that i think is related with the interaction between OpenLayers and my Google map Base Layer (inconsistencies with coordinates systems, or something like that). Maybe also this could be related with the request that im sending to GeoServer through WMS. But i'm confused. My setup: Like i said, i have an olMap with a base layer that comes from google. Here is the map init:               function initMap() {              // World Geodetic System 1984 projection (lon/lat)         var WGS84 = new OpenLayers.Projection("EPSG:4326");              // WGS84 Google Mercator projection (meters)         var WGS84_google_mercator = new OpenLayers.Projection("EPSG:900913");              var options = {            projection: WGS84_google_mercator,            displayProjection: WGS84,         };                   window.map = new OpenLayers.Map('map_canvas', options);              var google_default =  new OpenLayers.Layer.Google("Google Default", {wrapDateLine: false, numZoomLevels: 30, sphericalMercator: true})         var google_satellite = new OpenLayers.Layer.Google("Google Satellite", {type: google.maps.MapTypeId.SATELLITE, wrapDateLine: false})              map.addLayers([google_default]);                   window.map.zoomTo(2);         window.map.setCenter(0);                  // build up all controls         window.map.addControl(new OpenLayers.Control.Zoom());         window.map.addControl(new OpenLayers.Control.MousePosition())         window.map.addControl(new OpenLayers.Control.Navigation());         window.map.addControl(new OpenLayers.Control.LayerSwitcher());          }      When the zoom is small the server send me the response with the feature info, but when zoom in the server no longer responds. I tried to fix that passing buffer param region to make the intersection in the server side. But i have no good results when zooming even with an excessive value for that:               vendorParams: {        buffer: 100, // one hundred is an excessive value     },       Also i'm printing out the event.xy that is different when zooming but that make sense, isn't it? An Alert is launched in the beforegetfeatureinfo method of the WMSGetFeatureInfo Control listener:               beforegetfeatureinfo function(event) {                       // Code here to set the content of queryableMapLayers                       // The event object will contain xy of mouse click                       alert(event.xy);                     }      Without zoom i get (1136, 273) values and (600, 304) with large zoom approximately in the same map location. Like i said, i'm confused about the issue. I spend much time googling it but i didn't found any guideline. Any thoughts?


title: GeoDatabase topology and rules with a *.rul file

text: I have been working for some time on a solution to a problem much like this question.   I have several layers that I want to resolve together. I have them all in topology (geodatabase) and have been trying to decide which rules and in what order would be best practice or even would accomplish.   For instance the same data in the overlapping sections needs to not overlap or have gaps.   That data should not move where a county polygon is near and should be moved to it.   As should the school districts.   All data layers are polygonal.   Ending with all features being moved and made exact (within tolerance) to the plss layer which is first repaired for gaps and slivers. The topology tool allows you to load a *.rul file, which I assume is an ordered XML. I haven't found one anywhere to examine it.   To my question does anyone hold a *.rul file that would begin to give me insight into the rules and the order (the tolerances I will most likely have to tweak for my data)? If user823 would ever come back by I might convince him to post his complex rule set.


title: Postgis Raster merge geometry

text: I have a raster table and a geometry table in postgis. My goal is to create a single raster merging elements from both:-               WITH geom_raster AS (              SELECT wkb_geometry, ST_AsRaster(wkb_geometry, 200,200,ARRAY['8BUI', '8BUI', '8BUI'], ARRAY[118,154,118], ARRAY[0,0,0]) as rast FROM geometry_table     ),     raster AS (              SELECT ST_Union(rast) from raster_table,geom_raster where wkb_geomtry && rast     )        SELECT ............      What is the last step to merge the two into a single raster so I see the geometry on top of the raster? I tried ST_Union but get an error about alignment which I suspect means the grids or SRIDS or something else are not in unison?


title: Search Functionality using Cartodb API

text:
I posted a similar question about this a month or so ago and I was still in early stages of researching the problem thus not to many results came from it .After some time and effort I have come close to having my desired results. The goal is to have a search box that queries the table on my cartodb account and then results get sent back to the API. The code below works as described however that is because I have the value of the form hard coded. I did this to show that the function is correct and that the value does get passed on the the query.My problem thought is once i delete the hard coded value and try to type the text in, the query results on my Console are blank. So some where the value is getting dropped and I am not quite sure where. (I posted the whole code so that you can copy and paste to see the results...If desired)               <!DOCTYPE html>                 <title>Bounds example | CartoDB.js</title>     <meta name="viewport" content="initial-scale=1.0, user-scalable=no" />     <meta http-equiv="content-type" content="text/html; charset=UTF-8"/>     <link rel="shortcut icon" href="http://cartodb.com/assets/favicon.ico" />                  <link rel="stylesheet"            href="http://libs.cartocdn.com/cartodb.js/v3/themes/css/cartodb.css" />

     <form action="#" method="post"  id="form">         <input type="Text" name="userInput" id="test" value="ARNOLD RD">          <input type="submit" value="Submit" onclick="processFormData();" />     </form>


title: FME : How can I use multiple Workspacerunners in one Workspace?

text: I created two Workspaces, the first one does a clipping process and the second has to run it. It contains a "Creator" and 4 Workspace-runners. I have 4 different Inputs and outputs paths, how can I launch my main workspace and only specify in each workspacerunner the INPUTS and the OUTPUTS? Thanks !![enter image description here](http://i.stack.imgur.com/i9Nok.jpg) * * * I'm sorry I wasn't very explicit. What you said describes what I wanted to do. I have a Dataset that I use each time "communal areas" and 4 others "pipelines" that I change each time. What I want is that I don't have to change my INPUTS/OUTPUTS directories each time. This is why I used Workspace- runners but it didn't work. I think I have a problem in my main workbench because I had this Error Message : * * *               2013-09-13 17:04:31|   0.4|  0.0|INFORM|WorkspaceRunner_2: Running FME Workspace W:\LAB\...decoupage_canalisations.fmw with command line:     2013-09-13 17:04:31|   0.4|  0.0|INFORM|WorkspaceRunner_2: {C:/Program Files/FME_2013/fme} W:/LAB/.../Decoupage_cana_dep/decoupage_canalisations.fmw --SourceDataset_SHAPE {W:\LAB\...\INPUTS\DIFF_MPC\*.shp} --SourceDataset_SHAPE_3 {W:\LAB\..._REGION\*.shp} --DestDataset_SHAPE {W:\LAB\...OUTPUTS\DIFF_MPC} --LOG_FILE {W:\LAB\...\LOGS\DIFF_MPC} 2>@1     2013-09-13 17:04:35|   0.4|  0.0|WARN  |WorkspaceRunner_2: Failed to run workspace W:/LAB.../decoupage_canalisations.fmw -- Trace was:                           FME 2013 SP2 (20130605 - Build 13499 - WIN32)                       FME Desktop Professional Edition (floating)                                    Permanent License.                             Machine host name is: L50026151                     Copyright (c) 1994 - 2013 by Safe Software Inc.                                    Safe Software Inc.               Unable to open log file `W:\LAB\...\LOGS\DIFF_MPC'     Program Terminating      * * * Translation FAILED. I Put another screenshot of the main workbench.! ![enter image description here](http://i.stack.imgur.com/i9Nok.jpg) ![enter image description here](http://i.stack.imgur.com/q6Uze.jpg)


title: Copy features from feature class to an existing feature class with arcpy

text: I see CopyFeatures_management, but this creates a new feature class. http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//001700000035000000 Is there a way to copy features from a feature class into an existing feature class?


title: How can I use ogrinfo to reach information anout a .shp?

text: How can I use ogrinfo?   where should I write it?   In the visual studio's command prompt or somewhere else?   Can you make some examples? * * * when I use `ogrinfo` in the windows or visual studio's command prompt, I get the result: ![enter image description here](http://i.stack.imgur.com/kL70C.jpg)   This is the same result that I get, when I use `gdal_translate`.


title: Has ArcMap's background progress bar been fixed in 10.2?

text: I don't know if any other geoprocessing tool authors are pondering whether to invest the time and effort in a workaround for this "feature", but   * in 10.0 background processing tasks had a progress bar indication both at the bottom of the main window and as a little 'clock' in the results window   * in 10.1 the progress indication is gone, for background processes Has anyone installed 10.2 yet? Have we got it back?


title: Memory considerations for cloud based GeoStack

text: Some time ago I asked about hardware considerations of cloud based GeoStack implementation. As a follow up to this question There are some good pointers in related question already, but they mainly refer to PostGIS. So I wanted to ask you for some advice regarding memory allocation for a simple Web mapping Ubuntu box, based on PostGIS, GeoServer and OpenLayers/GeoEXT. By simple in this case I mean serving maximum 10-20 WMS layers, caching, getFeatureInfo, and getting legends. Underlying dataset will consists of ~3000 polygons of administrative areas. It will be hosted on cloud VPS - either AWS or rackspace. I have three questions regarding this setup:   1. Does it make sense to keep PostgreSQL and Geoserver on two separate instances? Or is it better to focus on one, more powerful server?   2. Cloud servers are scalable. So assuming 100% is the total amount of memory available - what would be your suggestions of setting memory parameters in PostgreSQL? And how about java memory allocation for GeoServer?   3. Any other hints for tweaking memory? Thanks for help. PS I asked this question on GeoServer Users mailing lists but haven't recieved any suggestions so far.


title: Path to proxy.cgi for GeoServer is not accessible in browser

text: There have been several questions very similar to this asked already (including ones not involving GeoServer), but their solutions have not solved my problem. I am trying to use getFeatureInfo on my GeoServer WMS layers, which are hosted on a different domain than the website on which they will be used, and have since learned that I need a proxy to accomplish this. I have followed the instructions on how to set this up in Tomcat 6 (except for renaming .jar file as it's unnecessary in v6), added www.greengeography.org and the IP of my local network to the proxy.cgi file, then changed the top line to `#!/usr/bin python` as that's where it's at in Amazon's EC2 Linux (and also tried a slash instead of space as I've seen both ways). I am confused on URL mapping, servlet mapping, and where to create a cgi directory, so I tried it in just about every location within the /tomcat6/webapps, including the main geoserver directory. The only place I am able to access the cgi file in a browser is at the root. Why can I not make this happen in the geoserver directory, which is where I think it needs to be? I did not see instructions about changing the URL mapping, so I'm not positive it's the problem. As suggested in several other questions posted on this site, I attempted to install the GeoServer Proxy Extension instead of configuring this stuff manually, but every time I tried to access it from the GS admin page, I got java errors (no matter which version of the extension I tried), so I gave up and went back to manual. If it helps, here is some of my javascript, and I have tried just about every URL variation I could think of:               OpenLayers.ProxyHost="http://50.112.123.161:8080/geoserver/cgi-bin/proxy.cgi?url=";      And the GeoExt part:               var featureInfo = new OpenLayers.Control.WMSGetFeatureInfo({                  url: 'http://50.112.123.161:8080/geoserver/gwc/service/wms',                  queryVisible: true,                      exceptions: "application/vnd.ogc.se_xml",                         layerUrls: ["http://50.112.123.161:8080/geoserver/OR_wild_wkspc/wms","http://50.112.123.161:8080/geoserver/gwc/service/wms"],                     eventListeners: {                         getfeatureinfo: function(e) {                              new GeoExt.Popup({                              title: "Feature Info",                              width: 400,                              height: 300,                              autoScroll: true,                              queryVisible: true,                              map: map,                              feature_count: 50,                              lonlat: mapPanel.map.getLonLatFromPixel(e.xy),                              html: e.text                             }).show();                         }                     }             });               mapPanel.map.addControl(featureInfo);          featureInfo.activate();      I have also been getting these errors in my GeoServer logs:               2012-02-13 23:52:22,922 ERROR [config.XMLConfiguration] - Configuration file cannot be read or does not exist!     2012-02-13 23:52:22,923 ERROR [layer.TileLayerDispatcher] - Configuration /usr/share/tomcat6/webapps/geoserver/../resources contained no layers.     2012-02-13 23:52:22,924 ERROR [layer.TileLayerDispatcher] - Error reading service information from /usr/share/tomcat6/webapps/geoserver/../resources: Error parsing file /usr/share/tomcat6/webapps/geoserver/../resources/geowebcache.xml      I also get ones when a file cannot be found. This one happens to be an html file since it was just a test to see if I could figure out the paths, but the rest of them are for proxy.cgi:               2012-02-13 19:50:22,785 WARN [servlet.PageNotFound] - No mapping found for HTTP request with URI [/geoserver/OR_wild_wkspc/test.html] in DispatcherServlet with name 'dispatcher'      Last but not least, I also tried more paths to point to the .cgi file, as suggested here. But no matter what, I get this error in Chrome Developer: > XMLHttpRequest cannot load > http://50.112.123.161:8080/proxy.cgi?url=http://50.112.123.161:8080http%3A%2F%2F50.112.123.161%3A8080%2Fgeoserver%2FOR_wild_wkspc%2Fwms%3FSERVICE%3DWMS%26VERSION%3D1.1.1%26REQUEST%3DGetFeatureInfo%26LAYERS%3DOR_wild_wkspc%253AExistingWildWebMerc%2COR_wild_wkspc%253AProposedWildWebMerc%26QUERY_LAYERS%3DOR_wild_wkspc%253AExistingWildWebMerc%2COR_wild_wkspc%253AProposedWildWebMerc%26STYLES%3D%2C%26BBOX%3D-14061041.537506%252C5096118.80202%252C-12776899.462494%252C5842144.19798%26FEATURE_COUNT%3D10%26HEIGHT%3D610%26WIDTH%3D1050%26FORMAT%3Dimage%252Fpng%26INFO_FORMAT%3Dtext%252Fhtml%26SRS%3DEPSG%253A900913%26X%3D409%26Y%3D296. > Origin http://www.greengeography.org is not allowed by Access-Control-Allow- > Origin. OpenLayers.js:516Uncaught TypeError: Cannot read property 'x' of > null If you feel I should've been able to solve this given that there have been so many similar questions asked already, that's fine, but please don't downvote it, because I did do my research. Thanks!


title: cartodb runs out of space *after* shrinking the data in the tables?

text: I try to work in CartoDB with data I import from CSVs (or xlsx, see below). The data is 24.1 MB in ASCII CSV, 10 MB in xlsx. I spent quite some time shrinking my data's size, and I was happy to upload less than 25 MB (the limit of my free academic account). I am puzzled by what happened next. I could delete two longish string columns (latitude and longitude, as cartodb georeferenced it automatically). I could convert a column of years to numbers from strings, which should only shrink the size of the table. I could convert a column of dates string to date, which should shrink it further. But when I would do one final operation I need, the conversion of an intensity measure of my temporal heat map from string to number, which can only shrink the size further, THEN cartodb complains that I am over the limit of my account? I spent hours deleting my tables and importing them again, twice in xlsx where I thought numbers will be numbers from the get-go, but that import even fails completely. I would rather not pay for an upgrade until I understand what is going on and why it is necessary. I have also mixed feelings about going over the limit because of the created_at and updated_at columns, which cartodb creates automatically and does not let us delete… Any thoughts on how to work around this? Any experience with this behavior?


title: Exporting large high resolution map to PDF in ArcMap errors?

text: I am attempting to export a large E size map at 300 DPI in ArcMap 10.2.2 but each time I export there are errors on the map. Once the PDF is generated there are areas where symbols are stretched out in a line across the width of the map. I have tried it on ArcMap 10.1 and 10.2.2 and I get the same errors. I know I can avoid these errors if I export at a lower DPI, but I want a really high resolution print of this particular map to be prominently displayed. All other parameters for the export to PDF are the default values. Here is an example of the errors: ![example of the errors.](http://i.stack.imgur.com/pLsqI.png)


title: Personal Geodatabase: Access speed up due the use of RelationshipClasses?

text: I've a special question about the effects of RelationshipClasses. For some reasons I've created a Personal Geodatabase and imported two features from Shape-Files and created an own table. For example we have the following structure feature_I = ID_K, Cnt, Closed feature_II = ID_S, Name, Number Table = ID_X, ID_Y, ID_Z,... The table contains 101046 rows, the ID_K is equal with ID_X, ID_S is equal with ID_Y and ID_Z. The value of ID_K is at least available twice in the table (ID_X). Now I've written some algorithms which use cursors to get through the all entries in feature_I. Then the algorithms will look at the table and filter all rows for which the statement ID_K = ID_X is true. For each of this rows two queries are created. This first query delivers the object of the feature_II (ID_Y = ID_S) and the second query will deliver a 2nd object of the feature_II (ID_Z = ID_S). I've measured the time for the whole operation. The result varies between 80 and 90 min. Now I've created three RelationshipClasses between the Table and the two features. All are one-to-many relations (ID_K with ID_X, ID_S with ID_Y and ID_Y with ID_Z). If I rerun the whole algorithm the execution time drops to 19 min. Now the questions is why? Whats the reason for the increase of performance by accessing the tables in the geodatabase when RelationshipClasses exist? Thanks for the clarification.


title: How to properly add a Python Module to a Plugin?

text: I've been trying to use a module that is installed in my local version of python, but whenever I open the plugin itself, it does not manage to locate said module. The module in this case is Noise 1.1.1. The error I get is this, is there somewhere I need to declare the module so it allows its use in qgis?                   Couldn't load plugin PointGenerator due an error when calling its classFactory() method          Traceback (most recent call last):       File "C:/PROGRA~2/Quantum GIS Lisboa/apps/qgis/./python\qgis\utils.py", line 164, in startPlugin         plugins[packageName] = package.classFactory(iface)       File "C:/Users/Chris/.qgis//python/plugins\PointGenerator\__init__.py", line 52, in classFactory         from pointgenerator import PointGenerator       File "C:/PROGRA~2/Quantum GIS Lisboa/apps/qgis/./python\qgis\utils.py", line 309, in _import         mod = _builtin_import(name, globals, locals, fromlist, level)       File "C:/Users/Chris/.qgis//python/plugins\PointGenerator\pointgenerator.py", line 26, in          from noise import snoise2       File "C:/PROGRA~2/Quantum GIS Lisboa/apps/qgis/./python\qgis\utils.py", line 309, in _import         mod = _builtin_import(name, globals, locals, fromlist, level)     ImportError: No module named noise


title: Remove Shapefile lock

text: I have an ArcObject Dissolve procedure that creates a shp in a directory. This dissolved layer is a temp layer and is not added into ArcMap. I would like to delete it when the user clicks a button again, however it will not delete because their is a lock file (Dissolved.shp.6H8Z5D1.3448.1520.sr.lock) in the directory that happens after the dissolve gp runs:               pSumTable = pGeoProc.Dissolve(pInputTable, True, "ROUTE", "Dissolve.Shape, Minimum.ROUTE", pdatasetname2)      I looked at Using Scheme Locks, however I'm not using this layer to perform any edits. I also found this post on GIS SE, however it too does not seem to apply for my situation. How do I remove this sr lock file (my entire code below), while keeping the ArcMap session open? Thanks                Try                      'delete all temp files                 Dim s As String                 For Each s In System.IO.Directory.GetFiles("J:\20111030\GIS\Scripts\TempData2")                     System.IO.File.Delete(s)                 Next s                      'dissolve route layer                 Dim pMap As IMap                      pMap = mxdoc.FocusMap                      Dim pInputFeatLayer As IFeatureLayer                      pInputFeatLayer = featurelayer                      If pInputFeatLayer Is Nothing Then                          MessageBox.Show("The testinline layer was not found")                     Exit Sub                      End If                      Dim pInputTable As ITable                 pInputTable = pInputFeatLayer                      If pInputTable Is Nothing Then                          MsgBox("Table QI failed")                     Exit Sub                      End If                      If pInputTable.FindField("ROUTE") = -1 Then                          MsgBox("There must be a field named MUSYM layer ")                     Exit Sub                      End If                      Dim pInputFeatClass As IFeatureClass                 pInputFeatClass = featurelayer.FeatureClass                      Dim pfeatureclassname As IFeatureClassName                 pfeatureclassname = New FeatureClassName                      With pfeatureclassname                          .FeatureType = esriFeatureType.esriFTSimple                     .ShapeFieldName = "Shape"                     .ShapeType = pInputFeatClass.ShapeType                      End With                      Dim pNewWSName As IWorkspaceName                 pNewWSName = New WorkspaceName                      pNewWSName.WorkspaceFactoryProgID = "esriCore.ShapefileWorkspaceFactory.1"                 pNewWSName.PathName = "J:\20111030\GIS\Scripts\TempData2\"                      Dim pdatasetname2 As IDatasetName                      pdatasetname2 = pfeatureclassname                 pdatasetname2.Name = "Dissolved"                      pdatasetname2.WorkspaceName = pNewWSName                      Dim pGeoProc As IBasicGeoprocessor                 Dim pSumTable As ITable                 pGeoProc = New BasicGeoprocessor                 pSumTable = pGeoProc.Dissolve(pInputTable, True, "ROUTE", "Dissolve.Shape, Minimum.ROUTE", pdatasetname2)                 Dim pOutputFeatClass As IFeatureClass                 pOutputFeatClass = pSumTable                      If pOutputFeatClass Is Nothing Then                     MsgBox("FeaturesClass QI Failed")                     Exit Sub                 End If                      Dim pOutputFeatLayer As IFeatureLayer                 pOutputFeatLayer = New FeatureLayer                 pOutputFeatLayer.FeatureClass = pOutputFeatClass                 pOutputFeatLayer.Name = pOutputFeatClass.AliasName                      ' CLEAN UP                 featurelayer = Nothing                 pInputFeatLayer = Nothing                 pInputTable = Nothing                 pInputFeatClass = Nothing                 pfeatureclassname = Nothing                 pNewWSName = Nothing                 pdatasetname2 = Nothing                 pGeoProc = Nothing                 pSumTable = Nothing                 pOutputFeatClass = Nothing                 pOutputFeatLayer = Nothing                      pMap = Nothing                 mxdoc = Nothing                  Catch ex As Exception                 MessageBox.Show("Caught an unspecified error in the calling code: " & vbCrLf & ex.ToString)             End Try


title: No intersection found between polygons I know intersect!

text: I'm using two geographies: one is made up of small areas (about 2200), the other large areas (about 70). Both cover all of Australia, and the small areas are not necessarily nested inside large areas. For the large areas, I was given a kml file, and a mapinfo file (both come from the Australian Bureau of Statistics). I converted both into .shp format, for use in R (using rgeos). For each of the small areas, I would like the proportion that fits in each large area. This can be done with gArea(gIntersection(small area, large area))/gArea(small area) Which I execute inside a foreach loop that returns a 2200 by 70 matrix of weights. My problem is that for some of the small regions, this method does not find any intersection (as in it returns zero for all larger areas). I have run it on shape files constructed from both the kml file and the mapinfo files, and areas with no intersection _are different_! What would be the most likely cause of this sort of problem? * * * Edit- Thanks guys for your help. I have put together a little zip file containing my two shape files and the R script I'm using. I'm pretty fresh to this (3 days) so reading up on projections. Thanks again. (The file is about 60 meg) https://dl.dropbox.com/u/63100926/Example.zip * * * Edit 2: Thanks to your advice, I now know that my small regions are in GRS80, and my big regions are in WGS84. So thanks! I am having a bit of difficulty finding out how to change them. Any pointers on where to look? Cheers.


title: Simulate flood to assess Asset inundation footprints

text: I am trying to evaluate a number of flood gauges to assess at what levels of river height different flood risks would need to be issued. The gauges are currently set to warn of flood a couple of hundred km downstream but those levels don't affect my AOI. For example a "Major" flood downstream will only result in a "Minor" flood in my area etc. I am trying to work out a workflow for this work which needs to be delivered in the next 3 days. We have Arcgis 10.1 (Advanced with 3D Analyst/Spatial Analyst), Global Mapper and QGIS/GRASS. Datasources are Lidar point clouds (4ppm), 1m DEM, 125mm imagery, gauge locations (3D vector) and streams. What's the best pathway? So far, ArcGlobe lets me easily fly through the terrain and visualize all but the lidar datasets. Is there someway I can add the buildings as 3D surfaces into this? I am evaluating tools such as LP360, Lidar Analyst and E3D3 to do this (will add results as I work through them) ![enter image description here](http://i.stack.imgur.com/uTxjU.jpg) What I am planning to do is create inundation 3D polygons based on the Gauge height plus 1m, plus 2m etc (up to 10m) and then intersect them with the asset vertical profiles to get an idea of what's affected at different flood levels and how this relates to a Minor-Major local flood. Global Mapper does the inundation as below but I can't figure out how to specify what height levels I want or extract a 3D polygon of the inundation zone. ![enter image description here](http://i.stack.imgur.com/UfbOF.jpg) My intention is to use this post as a thread for people to easily compare the real world functionality of the different options so your comments/suggestions/critiques are welcome.


title: PgRouting Driving Distance crashing postgis server

text: I haven't used PgRouting much so I'm a little in the dark here. However, I am attempting a driving distance calculation like the one described by Underdark here. My SQL is identical, I have a routable network etc etc. I have also successfully run this in the past on another database. However, when I run the simple piece of sql below, it completely crashes the Postgresql server with no error other than connection to server has been lost:               SELECT * FROM driving_distance('           SELECT gid AS id,               start_id::int4 AS source,               end_id::int4 AS target,               shape_leng::float8 AS cost           FROM network',           494902005,           50000,           false,           false)      I first spotted the crash when trying to run this from Python. I restarted the computer and tried to run the stripped down version above from PGAdmin. It still crashes. BTW the '494902005' value is my node id and I want to limit the search to 50km. Because I am not getting any feed back, I don't know where to start diagnosing the issue. I am using PostGIS 2/PostGreSQL 9.1 on 32bit Windows 7. Any ideas gratefully received!


title: How to use QGIS Time Manager sample data?

text: I've downloaded the Time Manager plugin for QGIS for Windows, version 1.7.1. I'm not sure whether I should really be asking a question on here, but since there doesn't seem to be an issue tracker on the new git, I want to post pictures, and I know the developer frequents this list, I figured I may as well. My problem is I can't seem to make it work. Using the test data "tweets.shp", I can get an animation, however it is backwards, i.e. the screen is full at the start of the timeline, and tweets are removed at the time they are posted, leaving nothing at the end of the timeline. Graphics attached: ![Start of timeline](http://i.stack.imgur.com/OA6VS.png) ![End of timeline](http://i.stack.imgur.com/iVujY.png) The settings I am using are: **start** : t **end** : none **index** : tweets20120109101002200 **timeformat** : %Y-%m-%d %H:%M:%S **offset** : 0 I installed it manually (due to firewall issues I think) into " _\Program Files\Quantum GIS Wroclaw\apps\qgis\python\plugins\timemanager_ " - wonder if it's something to do with that. Alternatively, is there a way I could debug the queries it is making to QGIS for each time point?


title: Exception from HRESULT Error - Update Attributes

text: I have a form that allows users to update attribute information for selected features...but I'm recieving the error: "Exception from HRESULT: 0x80041051" when it tries to update the data. I'm using VB.net, Visual Studio Express 2008 and ArcGIS 10. Any suggestions on where I went wrong with updating the data?               System.Runtime.InteropServices.COMException was caught       ErrorCode=-2147217327       Message="Exception from HRESULT: 0x80041051"       Source="ESRI.ArcGIS.Geodatabase"       StackTrace:            at ESRI.ArcGIS.Geodatabase.IFeatureCursor.UpdateFeature(IFeature Object)            at ImprovementProjectsLines.ImprovementLines.cmdUpdate_Click(Object sender,     EventArgs e)       InnerException:      CODE:                    Private Sub cmdUpdate_Click(ByVal sender As System.Object, ByVal e As System.EventArgs) Handles cmdUpdate.Click         Try          Dim pMxDoc As IMxDocument             pMxDoc = My.ArcMap.Document             Dim pMap As IMap             pMap = pMxDoc.FocusMap             'Dim pFLayer As IFeatureSelection                       'Get a reference to the editor.             Dim editor As IEditor             editor = My.ArcMap.Application.FindExtensionByName("esriEditor.Editor")                  'Get Crossings FeatureLayer             Dim pFeatLyr As IFeatureLayer             pFeatLyr = GetLayerByTOC("TransportationImprovementProjects_Lines")                  'Get Dataset from FeatureLayer             Dim pEditDataset As IDataset             pEditDataset = pFeatLyr                  'If not editing, start editing dataset's Workspace             If editor.EditState <> esriEditState.esriStateEditing Then                 editor.StartEditing(pEditDataset.Workspace)             End If                  'Get SelectionSet             Dim pSelSet As ISelectionSet             pSelSet = GetSelection(pFeatLyr)                  'Get FeatureCursor from SelectionSet             Dim pFCursor As IFeatureCursor = Nothing             pSelSet.Search(Nothing, False, pFCursor)                  'Initialize first (only) feature             Dim pFeat As IFeature             pFeat = pFCursor.NextFeature                  'Update attributes with user input             'Fill in Attributes of selected feature in User Form             Do Until pFeat Is Nothing                 pFeat.Value(pFeat.Fields.FindField("YEAR")) = txtPrjYear.Text                 pFeat.Value(pFeat.Fields.FindField("PROJECTYPE")) = txtPrjType.Text                 pFeat.Value(pFeat.Fields.FindField("LOCATION")) = txtLocation.Text                 pFeat.Value(pFeat.Fields.FindField("AGENCY")) = txtAgency.Text                 pFeat.Value(pFeat.Fields.FindField("DESCRIPTION")) = txtDescrip.Text                 pFeat.Value(pFeat.Fields.FindField("COST")) = txtCost.Text                 pFeat.Value(pFeat.Fields.FindField("DATEEDITED")) = txtDate.Text                 pFeat.Value(pFeat.Fields.FindField("USEREDITED")) = txtUser.Text                 pFCursor.UpdateFeature(pFeat)                      pFeat = pFCursor.NextFeature             Loop                       'Stop editing and save edits             editor.StopEditing(True)                  'Close form             Me.Close()              Catch ex As Exception             Windows.Forms.MessageBox.Show(ex.ToString, "error")                   End Try     End Sub          Public Function GetSelection(ByVal pFeatLyr As IFeatureLayer) As ISelectionSet              'Gets selection of feature layer and returns selection set              'Initialize the required variables         pMxDoc = My.ArcMap.Application.Document         Dim pFeatSel As IFeatureSelection         pFeatSel = pFeatLyr              GetSelection = pFeatSel.SelectionSet          End Function          Public Function GetLayerByTOC(ByVal lyrName As String) As IFeatureLayer              'This function finds a feature layer based on its TOC Name              'Initalize global variables         pMxDoc = My.ArcMap.Application.Document         Dim pMap As IMap         pMap = pMxDoc.FocusMap              Dim lyrCntr As Integer         Dim pFeatLyr As IFeatureLayer              For lyrCntr = 0 To pMap.LayerCount - 1                  'Ensure that the layer is valid             If pMap.Layer(lyrCntr).Valid = True Then                      'Ensure that the layer is a feature layer                 If TypeOf pMap.Layer(lyrCntr) Is IFeatureLayer Then                     pFeatLyr = pMap.Layer(lyrCntr)                          If UCase(pMap.Layer(lyrCntr).Name) = UCase(lyrName) Then                         GetLayerByTOC = pFeatLyr                         Exit Function                     End If                 End If             End If         Next lyrCntr              Return Nothing         'Layer not found. Show a message         MsgBox("Layer " & lyrName & " not found !", vbCritical, "Error")          End Function      End Class


title: Zipcode overlay on OSM

text: I have to display an overlay with different colors depending on the zipcode in Europe. Can I get this information from OSM shapefile ? If yes, how exaclty ? Is there a tutorial somewhere ? Thanks


title: How to access tiles in Tilecache using QGIS?

text: I've Mapserver set up with a mapfile that points to a tileindex. In QGIS I can call the mapfile and view the rasters dynamically. I've got access to tilecache and I've seeded the same mapfile/rasters that I view in QGIS. My question is, what do I need to do to get QGIS to access the tilecache as opposed the dynamic rasters? I'm guessing that I need to do either something to my mapfile to get it to point to the seeded rasters as opposed to the dynamic rasters or repoint the WMS in QGIS


title: PyQGIS + PostGis: FATAL: sorry, too many clients already

text: I searched for questions regarding this problem and hadn't find one. I'm developing a custom application with a map canvas to show some geographic information from a PostGis database. After some testing I began to receive the following error: > FATAL: sorry, too many clients already This error means that the maximum number of database connections was reached. The layers I get with instances of QgsVectorLayer, which does not have any method to close the connection, unlike QgsRasterLayer which has the method closeDataProvider(). Should I assume that after getting the data from the database, the QgsVectorLayer instance automatically closes the connection? Or is there a proper method to do it at the end of an application? Thanks in advance.


title: How to compute surface of raster or any point using ARCGIS?

text: I have a raster image whose values are a digital elevation model of terrain. I need to know the surface type (topographic characterization) of particular points along transects of this DEM. How to find it? (In my application I have two points: a transmitter end and receiver end. I find the highest point between the transmitter end and receiver end: see the figure below. I would like to know whether the highest point is smooth or is a knife edge; that is, to within the resolution of the image it appears to be non-differentiable.) ![smooth or knife edge surface calculation](http://i.stack.imgur.com/4uFlj.jpg)


title: Loading Shapefiles into PostGIS using shapefile import/export manager

text: I am a complete beginner with PostGIS. I've started learning it using the tutorials provided by Boundless. In section 4 of this tutorial where I want to import shapefile in PostGIS, I do not have the PostGIS shapefile Import/Export Manager tool in the PgAdmin III. I've followed the tutorial step by step. They have never pointed if I need to install such tool or how to activate it, etc. Any idea what the problem could be? and how I can solve this? I've installed Postgresql 9.3, PostGIS 2.1 , pgAdmin3 1.12.3 on Windows Server 2012 R2 Standard. by the way, the relevant page for this issue in the tutorial is at: http://workshops.boundlessgeo.com/postgis-intro/loading_data.html Thanks!


title: How to translate 1980 US Census Place Code to post 1990 Place FIPS?

text: US Census PLaces are coded differently pre-1990 and post-1990, after the introduction of FIPS Code. For instance, for Albertville, AL: Pre-1990 : Place Code - PLACEA : 0030 Post-1990: Place FIPS - PLACEFP : 00988 Is there an easy way to find a correspondance between these two codes?


title: How to prevent symbol displacement when loading graphics on ArcGISDynamicMapServiceLayer?

text: I have a map with a ArcGISDynamicMapServiceLayer. In response to the user click on a button, some features are loaded on a GraphicsLayer. The GraphicsLayer renderer is set by Renderer.fromJSON(descService.drawingInfo.renderer), using the renderer of the mapService. When the loaded graphics are displayed on the GraphicsLayer, they are not located exactly where they should. They are displaced. In order to place them where they should be, a zoom is necessary. Do you know how to resolve this problem?


title: How to calculate the incidence angle of a polyline?

text: How can I calculate the angle of each polyline in ArcGIS?


title: How to remove coincident points with z values greater than min for Spline

text: I am trying to run a Spline interpolation on some terrestrial LiDAR data, but it keeps failing (I'm getting error 010235 - returned error from spline process and error 010067 - error executing grid expression). I think it's due to the fact that there may be coincident points (ie. points with the same xy but different z values). Can anyone suggest a way to delete the duplicate points with z values greater than the minimum? I can't do it manually because I think there may be a lot of duplicates (I tried the collect events tool, and there were many coincident points).


title: How can I find an OSGB six-digit code's bounding box?

text: How can I convert Ordnance Survey six-digit codes to the latlng bounding box of the related 1km grid square? For example, NH0325: according to Wikipedia this means a 1km square whose south-west corner is 3 km east and 25 km north from the south-west corner of square NH. Could anyone point me at some code for converting this code to its 1km grid square, ideally in Python? I already have Paul Agapow's OSGB library - I think this returns the south- west corner of the grid square, in OSGB36 coordinates. I could potentially use this to find the north-east corner, but I don't know how. I'm displaying the results on an WGS84 map, so it would be better to have the results in WGS84, if possible. This is all far beyond my knowledge of map projections. Any suggestions very much appreciated.


title: How can I use QGIS to convert points with elevations to DXF?

text: I'm trying to convert points (with their elevations) into a dxf file for use in AutoCAD LT. I've tried "save as" to dxf for the layer, which appears to almost work...the dxf has x and y location information, but the elevations associated with the points are not present. I also tried the "dxfauth" openware tool, but that did the same as the above. Thanks!


title: adding WMS service to application ArcGIS JS api

text: I've got map in projection: 4326 (WGS-84) when I try to add WMS layers which is in different projection I always get `connection timeout` Do you know if I could change WMS projection to set it up compatible with my projection? When I was using JAVA sdk I could connect WMS in different projections and arcgis server changed it projection ad hoc.


title: GeoServer WFS Row Level Security?

text: I'm trying to implement row level security with GeoServer. I know out of the box it supports service, layer and REST security and not row level. However I'm sure there must be a work around. Has anyone successfully implemented it?


title: Get all columns from postgis table when doing a getFeatureInfo request from mapserver's WMS layer

text: ## some background:   * im using openlayers WMSGetFeatureInfo control to get feature info from a WMS layer.   * the data is in postgis (its a table of points).   * im using mapserver, and you need to define a template file for getfeatureinfo request. * * * if i put in the template file: > any string it works. when clicking a feature i get back this string (for any feature of course). if i put in the template file: > [date] it returns the **current** date, (not a column in the table) if i put in the template file: > [id] it returns a number that is not the id column of the feature, though it is different for each feature. * * * how can i request all columns for that feature, as if i would do a "select *" quary ? thanks allot!


title: Problems converting CSV to Shapefile using pyshp

text: I have a python script which converts csv data to shapefiles. It converts but won't get angle from csv. Can anyone help me? The CSV data looks like this:               fid_nr;x;y;angle     1;55.474040000;13.020460000;90     2;55.474040000;13.020460000;90      And the Python script is:               import shapefile as shp     import csv     import os     import sys          fileToOpen = open('test.csv', "rb")          #Set up blank lists for data     fid_nr,x,y,angle=[],[],[],[]          #read data from csv file and store in lists     with open(fileToOpen, 'rb') as csvfile:         print "File found and will be now converted"         r = csv.reader(csvfile, delimiter=";")         for i,row in enumerate(r):             print row             if i > 0: #skip header                 fid_nr.append(row[2])                 x.append(float(row[0]))                            y.append(float(row[1]))                          angle.append(row[3])                   #Set up shapefile writer and create empty fields     w = shp.Writer(shp.POINT)     w.autoBalance = 1 #ensures gemoetry and attributes match     w.field("fid_nr","N")     w.field("x","F",10,8)     w.field("y","F",10,8)     w.field("angle","N")              #loop through the data and write the shapefile     for j,k in enumerate(x):         w.point(k,y[j]) #write the geometry         w.record(k,y[j],fid_nr[j],x[j],y[j],angle[j])     #Save shapefile     w.save("Resultat/" + out_file3)     print "Done!"


title: QGIS 2.2 documentation

text: I thought I would try to write a plugin for QGIS (some Python code that I hoped I could easily port) I started following http://www.qgisworkshop.org/html/workshop/plugins_tutorial.html but that turned out to be hopelessly out of date. I turned to http://www.qgis.org/en/docs/pyqgis_developer_cookbook/intro.html instead and started going through the intro for the Python console and stumbled at the second query.               layer = qgis.utils.iface.activeLayer()      worked fine but then               layer.getLayerID()      threw up               Traceback (most recent call last):     File "<input>", line 1, in <module>     AttributeError: 'QgsVectorLayer' object has no attribute 'getLayerID'      and sure enough code completion in the console would seem to suggest that this doesn't exist. So my question is "where do I find up to date, usable instructions for developing plugins for QGIS 2.2?" It would seem that unless you are already "in the know" there is no way in. Or am I just being paranoid?


title: Finding fields in Feature

text: I have a program that has to loop throug a very large amount of features. I have created a wrapper for the feature. This made me think: How expensive is the Fields.FindField method? Would it be faster - and if so; how much - to store the field index in a static field in the wrapper class. An example:               internal class XFeatureWrapper     {         private static int XColumnIndex =-99;         private string _xValue;               public FeatureWrapper(IFeature xFeature)         {             if (XColumnIndex == -99)             {                 XColumnIndex = xFeature.Fields.FindField("XXX");             }                  _xValue = xFeature.get_Value(XColumnIndex).ToString();              }              public string XValue         {             get { return _xValue; }         }     }      Would this be faster than finding the field index in the feature every time? Morten


title: Geotriggers multiply when using app

text: I managed to get the ESRI Geotrigger example running on my Android device. Created a small number of geotriggers (tag, location, boundary and rules) using the Geotrigger Editor. Works great. One thing is weird though: every time a trigger is fired a new geotrigger is created for the same location. Somewhere in the code a part is relevant for this problem, but it seems to contradict:                   // Create the trigger if we haven't done so already.         if (mShouldCreateTrigger) {             // Set create trigger flag here so that we don't create multiple             // triggers if we get a few initial updates in rapid succession.             mShouldCreateTrigger = false;      OK, Done. The next part is causing the problem I guess.                       // The TriggerBuilder helps build JSON parameters for use with the             // 'trigger/create' API route.             JSONObject params = new TriggerBuilder()                     .setTags(TAGS[0]) // make sure to add at least one of the tags we have on the device to this trigger                     .setGeo(loc, 100)                     .setDirection(TriggerBuilder.DIRECTION_LEAVE)                     .setNotificationText("You left the trigger!")                     .build();                  // Send the request to the Geotrigger API.             GeotriggerApiClient.runRequest(this, "trigger/create", params,                     new GeotriggerApiListener() {                 @Override                 public void onSuccess(JSONObject data) {                     Toast.makeText(GeotriggerActivity.this, "Trigger created!",                             Toast.LENGTH_SHORT).show();                     Log.d(TAG, "Trigger Created");                 }      Should I delete this part ?


title: QGIS Arch Linux compilation error

text: Just pulled QGIS from git and am trying to compile it, here's the relevant error: ` [ 89%] Building CXX object python/CMakeFiles/python_module_qgis_core.dir/core/sipcorepart3.cpp.o In file included from /home/noah/qgis/qgis-git/qgis/Quantum- GIS/python/core/qgsmaplayerregistry.sip:10:0: /home/noah/qgis/qgis- git/qgis/Quantum-GIS/python/../src/core/qgsmaplayerregistry.h: In member function ‘virtual void sipQgsMapLayerRegistry::connectNotify(const char*)’: /home/noah/qgis/qgis-git/qgis/Quantum- GIS/python/../src/core/qgsmaplayerregistry.h:252:10: error: ‘virtual void QgsMapLayerRegistry::connectNotify(const char*)’ is private void connectNotify( const char * signal ); ^ /home/noah/qgis/qgis-git/qgis/Quantum- GIS/build-master/python/core/sipcorepart3.cpp:21714:46: error: within this context return; ^ make[2]: *** [python/CMakeFiles/python_module_qgis_core.dir/core/sipcorepart3.cpp.o] Error 1 make[1]: *** [python/CMakeFiles/python_module_qgis_core.dir/all] Error 2 make: _*_ [all] Error 2 ` I'm currently downgrading python as suggested on the AUR page-- I don't understand the problem, I hope this is the right thing to do.


title: Installation question: initdb is missing after install of postgis2 via RPMS on RHEL 6.4

text: I'm attempting to install postgis2 on RHEL 6.4 by using the instructions postgis.net/install, but after the postgis2 install, neither initdb nor /etc/init.d/postgreql92 is present. The same thing happens with the 9.3 beta. Is a source installation my best option? **DETAILS**               [root@pgsrv ~]# cat /etc/issue     Red Hat Enterprise Linux Server release 6.4 (Santiago)          [root@pgsrv ~]# rpm -qa | grep post     postfix-2.6.6-2.2.el6_1.x86_64          [root@pgsrv ~]# rpm -qa | grep pgdg          [root@pgsrv ~]# rpm -ihv pgdg-redhat92-9.2-7.noarch.rpm     Preparing...                ########################################### [100%]        1:pgdg-redhat92          ########################################### [100%]          [root@pgsrv ~]# cd /etc/yum/pluginconf.d/     [root@pgsrv pluginconf.d]# ls     changelog.conf               product-id.conf          security.conf              versionlock.conf     downloadonly.conf            rhnplugin.conf           subscription-manager.conf  versionlock.list     epel-release-6-8.noarch.rpm  rhnplugin.conf-original  verify.conf     [root@pgsrv pluginconf.d]# cat rhnplugin.conf     [main]     enabled = 1     gpgcheck = 1          # Added "exclude=postgresql*" for PostGIS      exclude=postgresql*          # You can specify options per channel, e.g.:     #     #[rhel-i386-server-5]     #enabled = 1     #     #[some-unsigned-custom-channel]     #gpgcheck = 0               [root@pgsrv pluginconf.d]# yum install postgis2*     Loaded plugins: changelog, downloadonly, product-id, rhnplugin, security, verify, versionlock     This system is receiving updates from RHN Classic or RHN Satellite.     Setting up Install Process     Resolving Dependencies     --> Running transaction check     ---> Package postgis2_92.x86_64 0:2.0.3-2.rhel6 will be installed     --> Processing Dependency: postgresql92 for package: postgis2_92-2.0.3-2.rhel6.x86_64     --> Processing Dependency: libpq.so.5()(64bit) for package: postgis2_92-2.0.3-2.rhel6.x86_64     --> Processing Dependency: libgdal.so.1()(64bit) for package: postgis2_92-2.0.3-2.rhel6.x86_64     ---> Package postgis2_92-debuginfo.x86_64 0:2.0.3-2.rhel6 will be installed     ---> Package postgis2_92-devel.x86_64 0:2.0.3-2.rhel6 will be installed     ---> Package postgis2_92-docs.x86_64 0:2.0.3-2.rhel6 will be installed     ---> Package postgis2_92-utils.x86_64 0:2.0.3-2.rhel6 will be installed     --> Processing Dependency: perl-DBD-Pg for package: postgis2_92-utils-2.0.3-2.rhel6.x86_64     --> Running transaction check     ---> Package gdal.x86_64 0:1.7.3-15.el6 will be installed     ---> Package perl-DBD-Pg.x86_64 0:2.15.1-4.el6_3 will be installed     ---> Package postgresql92.x86_64 0:9.2.4-1PGDG.rhel6 will be installed     ---> Package postgresql92-libs.x86_64 0:9.2.4-1PGDG.rhel6 will be installed     --> Finished Dependency Resolution          Dependencies Resolved          ========================================================================================================================      Package                          Arch              Version                       Repository                       Size     ========================================================================================================================     Installing:      postgis2_92                      x86_64            2.0.3-2.rhel6                 pgdg92                          2.5 M      postgis2_92-debuginfo            x86_64            2.0.3-2.rhel6                 pgdg92                          1.5 M      postgis2_92-devel                x86_64            2.0.3-2.rhel6                 pgdg92                          172 k      postgis2_92-docs                 x86_64            2.0.3-2.rhel6                 pgdg92                          3.7 M      postgis2_92-utils                x86_64            2.0.3-2.rhel6                 pgdg92                           31 k     Installing for dependencies:      gdal                             x86_64            1.7.3-15.el6                  epel                            3.0 M      perl-DBD-Pg                      x86_64            2.15.1-4.el6_3                rhel-x86_64-server-6            196 k      postgresql92                     x86_64            9.2.4-1PGDG.rhel6             pgdg92                          970 k      postgresql92-libs                x86_64            9.2.4-1PGDG.rhel6             pgdg92                          185 k          Transaction Summary     ========================================================================================================================     Install       9 Package(s)          Total download size: 12 M     Installed size: 53 M     Is this ok [y/N]: y     Downloading Packages:     (1/9): gdal-1.7.3-15.el6.x86_64.rpm                                                              | 3.0 MB     00:00          (2/9): perl-DBD-Pg-2.15.1-4.el6_3.x86_64.rpm                                                     | 196 kB     00:00          (3/9): postgis2_92-2.0.3-2.rhel6.x86_64.rpm                                                      | 2.5 MB     00:00          (4/9): postgis2_92-debuginfo-2.0.3-2.rhel6.x86_64.rpm                                            | 1.5 MB     00:00          (5/9): postgis2_92-devel-2.0.3-2.rhel6.x86_64.rpm                                                | 172 kB     00:00          (6/9): postgis2_92-docs-2.0.3-2.rhel6.x86_64.rpm                                                 | 3.7 MB     00:00          (7/9): postgis2_92-utils-2.0.3-2.rhel6.x86_64.rpm                                                |  31 kB     00:00          (8/9): postgresql92-9.2.4-1PGDG.rhel6.x86_64.rpm                                                 | 970 kB     00:00          (9/9): postgresql92-libs-9.2.4-1PGDG.rhel6.x86_64.rpm                                            | 185 kB     00:00          ------------------------------------------------------------------------------------------------------------------------     Total                                                                                9.2 MB/s |  12 MB     00:01          Running rpm_check_debug     Running Transaction Test     Transaction Test Succeeded     Running Transaction     Warning: RPMDB altered outside of yum.       Installing : postgresql92-libs-9.2.4-1PGDG.rhel6.x86_64                                                           1/9        Installing : postgresql92-9.2.4-1PGDG.rhel6.x86_64                                                                2/9        Installing : gdal-1.7.3-15.el6.x86_64                                                                             3/9        Installing : postgis2_92-2.0.3-2.rhel6.x86_64                                                                     4/9        Installing : perl-DBD-Pg-2.15.1-4.el6_3.x86_64                                                                    5/9        Installing : postgis2_92-utils-2.0.3-2.rhel6.x86_64                                                               6/9        Installing : postgis2_92-devel-2.0.3-2.rhel6.x86_64                                                               7/9        Installing : postgis2_92-docs-2.0.3-2.rhel6.x86_64                                                                8/9        Installing : postgis2_92-debuginfo-2.0.3-2.rhel6.x86_64                                                           9/9        Verifying  : postgis2_92-debuginfo-2.0.3-2.rhel6.x86_64                                                           1/9        Verifying  : postgis2_92-docs-2.0.3-2.rhel6.x86_64                                                                2/9        Verifying  : postgis2_92-devel-2.0.3-2.rhel6.x86_64                                                               3/9        Verifying  : postgresql92-libs-9.2.4-1PGDG.rhel6.x86_64                                                           4/9        Verifying  : postgis2_92-2.0.3-2.rhel6.x86_64                                                                     5/9        Verifying  : postgresql92-9.2.4-1PGDG.rhel6.x86_64                                                                6/9        Verifying  : postgis2_92-utils-2.0.3-2.rhel6.x86_64                                                               7/9        Verifying  : gdal-1.7.3-15.el6.x86_64                                                                             8/9        Verifying  : perl-DBD-Pg-2.15.1-4.el6_3.x86_64                                                                    9/9           Installed:       postgis2_92.x86_64 0:2.0.3-2.rhel6                       postgis2_92-debuginfo.x86_64 0:2.0.3-2.rhel6                       postgis2_92-devel.x86_64 0:2.0.3-2.rhel6                 postgis2_92-docs.x86_64 0:2.0.3-2.rhel6                            postgis2_92-utils.x86_64 0:2.0.3-2.rhel6                          Dependency Installed:       gdal.x86_64 0:1.7.3-15.el6                               perl-DBD-Pg.x86_64 0:2.15.1-4.el6_3                                 postgresql92.x86_64 0:9.2.4-1PGDG.rhel6                  postgresql92-libs.x86_64 0:9.2.4-1PGDG.rhel6                           Complete!          [root@pgsrv pluginconf.d]# cd /usr/     [root@pgsrv usr]# cd pgsql-9.2     [root@pgsrv pgsql-9.2]# cd bin     [root@pgsrv bin]# ls     clusterdb   createuser  dropuser       pg_dump         pg_restore     psql          shp2pgsql     createdb    dropdb      pg_basebackup  pg_dumpall      pgsql2shp      raster2pgsql  vacuumdb     createlang  droplang    pg_config      pg_receivexlog  pg_test_fsync  reindexdb          [root@pgsrv bin]# cd /etc/init.d     [root@pgsrv init.d]# ls     abrt-ccpp         cgconfig   iptables      network                 postgresql-9.3.rpmsave  rpcsvcgssd      trace-cmd     abrtd             cgred      irqbalance    nfs                     psacct                  rsyslog         tuned     abrt-oops         crond      kdump         nfslock                 qpidd                   sandbox         udev-post     acpid             cups       killall       nscd                    quota_nld               saslauthd       wdaemon     arptables_jf      dbora      ktune         nslcd                   rdisc                   single          winbind     arpwatch          dovecot    lvm2-lvmetad  ntpd                    restorecond             spamassassin    xinetd     atd               firstboot  lvm2-monitor  ntpdate                 rhnsd                   spice-vdagentd     auditd            functions  mcelogd       oddjobd                 rhsmcertd               sshd     autofs            haldaemon  mdmonitor     portreserve             rngd                    svnserve     avahi-daemon      halt       messagebus    postfix                 rpcbind                 sysstat     blk-availability  ip6tables  netconsole    postgresql-9.3-data-pg  rpcgssd                 systemtap     certmonger        ipsec      netfs         postgresql-9.3-orig     rpcidmapd               tcsd          [root@pgsrv init.d]# yum provides */initdb     Loaded plugins: changelog, downloadonly, product-id, rhnplugin, security, verify, versionlock     This system is receiving updates from RHN Classic or RHN Satellite.     elrepo-s6-x86_64/filelists                                                                       | 132 kB     00:00          epel/filelists_db                                                                                | 7.6 MB     00:07          pgdg92/filelists_db                                                                              | 157 kB     00:00          rhel-x86_64-server-6/filelists                                                                   |  25 MB     00:03          rhel-x86_64-server-optional-6/filelists                                                          | 9.3 MB     00:01          rhel-x86_64-server-supplementary-6/filelists                                                     | 1.4 MB     00:00          rhn-tools-rhel-x86_64-server-6/filelists                                                         |  57 kB     00:00          postgresql92-server-9.2.4-1PGDG.rhel6.x86_64 : The programs needed to create and run a PostgreSQL server     Repo        : pgdg92     Matched from:     Filename    : /usr/pgsql-9.2/bin/initdb          postgresql92-debuginfo-9.2.3-2PGDG.rhel6.x86_64 : Debug information for package postgresql92     Repo        : pgdg92     Matched from:     Filename    : /usr/src/debug/postgresql-9.2.3/src/bin/initdb          postgresql92-server-9.2.3-2PGDG.rhel6.x86_64 : The programs needed to create and run a PostgreSQL server     Repo        : pgdg92     Matched from:     Filename    : /usr/pgsql-9.2/bin/initdb          postgresql92-debuginfo-9.2.4-1PGDG.rhel6.x86_64 : Debug information for package postgresql92     Repo        : pgdg92     Matched from:     Filename    : /usr/src/debug/postgresql-9.2.4/src/bin/initdb          [root@pgsrv init.d]# rpm -qa | grep post     postgis2_92-docs-2.0.3-2.rhel6.x86_64     postgresql92-9.2.4-1PGDG.rhel6.x86_64     postgis2_92-utils-2.0.3-2.rhel6.x86_64     postgis2_92-devel-2.0.3-2.rhel6.x86_64     postfix-2.6.6-2.2.el6_1.x86_64     postgis2_92-2.0.3-2.rhel6.x86_64     postgresql92-libs-9.2.4-1PGDG.rhel6.x86_64     postgis2_92-debuginfo-2.0.3-2.rhel6.x86_64          [root@pgsrv init.d]# rpm -ql postgis2_92| grep initdb     [root@pgsrv init.d]# rpm -ql postgis2_92 | wc         239     239   14920     [root@pgsrv init.d]#      [root@pgsrv init.d]# rpm -ql postgresql92 | wc         327     327   15799     [root@pgsrv init.d]# rpm -ql postgresql92 | grep initdb     /usr/pgsql-9.2/share/locale/cs/LC_MESSAGES/initdb-9.2.mo     /usr/pgsql-9.2/share/locale/de/LC_MESSAGES/initdb-9.2.mo     /usr/pgsql-9.2/share/locale/es/LC_MESSAGES/initdb-9.2.mo     /usr/pgsql-9.2/share/locale/fr/LC_MESSAGES/initdb-9.2.mo     /usr/pgsql-9.2/share/locale/it/LC_MESSAGES/initdb-9.2.mo     /usr/pgsql-9.2/share/locale/ja/LC_MESSAGES/initdb-9.2.mo     /usr/pgsql-9.2/share/locale/ko/LC_MESSAGES/initdb-9.2.mo     /usr/pgsql-9.2/share/locale/pl/LC_MESSAGES/initdb-9.2.mo     /usr/pgsql-9.2/share/locale/pt_BR/LC_MESSAGES/initdb-9.2.mo     /usr/pgsql-9.2/share/locale/ro/LC_MESSAGES/initdb-9.2.mo     /usr/pgsql-9.2/share/locale/ru/LC_MESSAGES/initdb-9.2.mo     /usr/pgsql-9.2/share/locale/sv/LC_MESSAGES/initdb-9.2.mo     /usr/pgsql-9.2/share/locale/tr/LC_MESSAGES/initdb-9.2.mo     /usr/pgsql-9.2/share/locale/zh_CN/LC_MESSAGES/initdb-9.2.mo     /usr/pgsql-9.2/share/locale/zh_TW/LC_MESSAGES/initdb-9.2.mo


title: Error message starting gis.m in GRASS

text: I finally broke down and bought Neteler and Mitasova's book, Open Source GIS: A GRASS GIS Approach, because I need to learn command-line GRASS. On trying to work through the book's examples, I'm finding that there may be a problem with my GRASS installation, or more likely with the installation of the tk framework. Here's the code and the error message that I'm getting: > gregory@Gregorys-MacBook> gis.m   >  ~ gregory@Gregorys-MacBook> /usr/local/bin/wish: line 2: > /usr/local/bin/../../../Library/Frameworks/Tk.framework/Versions/8.5/Resources/Wish.app/Contents/MacOS/Wish: > No such file or directory I installed GRASS and its dependencies using the binary installers from KyngChaos, so I'm pretty confident that the TK framework got installed correctly, but the error message seems to indicate otherwise. Also, I don't know how to interpret the section of the path in the error message that reads "/../../../". Can anyone suggest a fix for this?


title: Is ArcMap's Fixed Zoom In/Out algorithm documented anywhere?

text: I only occasionally use the **Fixed Zoom In/Out** buttons on the **Tools** toolbar of **ArcMap** and had always assumed that they zoomed in/out by a set scale factor. However, using **ArcGIS for Desktop 10.2** , I just set my scale to be **1:10,000,000** and then hit the **Fixed Zoom In** button **10 times** while recording the subsequent scale each time and was surprised to see:               10,000,000     8,000,000     6,000,000     5,000,000     4,000,000     3,000,000     2,500,000     2,000,000     1,500,000     1,250,000     1,000,000      Dividing each scale denominator by its successor shows a sequence of: 1.25, 1.33, 1.2, 1.25, 1.33, 1.2, 1.25, 1.33, 1.2 Curious about what happened in between the scales above I experimented a little to discover that the first issuance of the Fixed Zoom In button simply takes the scale to the first of those listed above and then starts following the same pattern. Try starting with 1:8,000,001 to see ...               8,000,001     8,000,000     6,000,000     5,000,000     4,000,000     3,000,000      This seems like quite a neat solution, but not one that I had realised was in place - does anyone know if it has been documented anywhere? I've looked in the Online Help without success. In answer to @DevdataTengshe, the scales currently in my available list (see below) seem to bear no relation to those used by Fixed Zoom In/Out. ![enter image description here](http://i.stack.imgur.com/ue1v1.png) However, when I add 1:7,500,000 into my available list it does seem to insert that as an additional scale step ...               10,000,000     8,000,000     7,500,000     6,000,000     5,000,000     4,000,000      This makes the algorithm being used even more mysterious to me.


title: Mapping Geocoded addresses onto Census shapefiles in R

text: I have a list of addresses that I have geocoded to get the corresponding longitude and latitude.               > AddressDetail <- geocode(match_address)     > coordinates(AddressDetail) <- c("lon", "lat")      Now I believe I should assign a projection so that R knows how the lon/lat coords are set up               > proj4string(AddressDetail) <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs")      Now I read in the shapefile that uses the NZMG projection               > Shape_MB = readOGR("/Work in progress/MB06_LV2/", "MB06_LV2")        Now this is where it starts to go wrong. I try transforming the projection of the lon/lat addresses into NZMG projection.               > AddressDetail <- spTransform(AddressDetail_v2, CRS("+proj=tmerc +lat_0=0 +lon_0=173 +k=0.9996 +x_0=1600000 +y_0=10000000 +ellps=GRS80 +units=m +no_defs")                                                           non finite transformation detected:     lon lat       Error in spTransform(xSP, CRSobj, ...) : failure in points      In addition: Warning message:     In spTransform(xSP, CRSobj, ...) : 9 projected point(s) not finite      So when I continue on to use the over function to find the corresponding census meshblock index of the polygon that the address is within I get the following error.               > temp <- over(AddressDetail_v2,Shape_MB,MB06)     Error: identicalCRS(x, y) is not TRUE      I'm a bit of a novice with geospatial stuff in R so would appreciate any advice!


title: What do these curly graphics represent on this map?

text: What do the two curly graphics to the left and right of this map represent? They are the only two on the map and seem to be attached to the border of an area. ![enter image description here](http://i.stack.imgur.com/dVlwD.jpg) I'm sorry I haven't posted the full map, I'm not sure I'm allowed.


title: How to batch clip to single-part shapefile using ArcGIS 9.3 and Python 2.5?

text: How do I clip a shapefile to a single-part shapefile using ArcGIS 9.3.1 and Python 2.5? Example: I have a set of lines (lines.shp) I am trying to clip to each record in a shapefile called Grids.shp _In grids.shp there is a field called [NAME] with a record for each grid; grid-1, grid-2, ..._ The grids represent Map extents, and the lines will be roads/streams/other base layers (could not post base layers publicly so have provided Lines.shp as an example). I have been trying to use the **SearchCursor** method but have some trouble defining objects and whatnot as I am new to Python. I got this same process to work in V10 but the V10 syntax is way different/simpler, especially for cursors and uses Python 2.6 EDIT (working script! thank you so much everyone): Clip_wCursor_v2a.py ESRI Forum with python script attached.               # ---------------------------------------------------------------------------     # Clip_wCursor_v2.py     # Created for ArcGIS 9.3.1     #   (clip generated by ArcGIS/ModelBuilder)     # Attempted cursor mods by a newb, Donovan     # ---------------------------------------------------------------------------          # Import system modules     import sys     import string     import os     import arcgisscripting          # Create the Geoprocessor object     gp = arcgisscripting.create(9.3)          # Load required toolboxes...     gp.AddToolbox("C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Analysis Tools.tbx")          # Designate a workspace     gp.workspace = r"C:\gis\data"          # Assign variables...     clipto = "grids.shp" # grids to clip to     clipit = "lines.shp" # features to be clipped          # Cursor time...     rows = gp.SearchCursor(grids)     row = rows.Next()     feat = row.Shape # used as the clip_features in the Clip_analysis     while row:         n = str(row.NAME) # assign a variable for the processing message based on a field         print "clipping to: " +n # tells you what grid was clipped         # Clip_analysis <in_features> <clip_features> <out_feature_class> {cluster_tolerance}         # output format, grid_1_c.shp, grid_2_c.shp, etc...         gp.Clip_analysis(clipit, feat, str(row.NAME)+"_c"+".shp", "") #update str if necessary         row=rows.next()          # reset the array...     del rows     del gp          print "...fin..."     print "no error checking happened, examine your workspace..."      appreciate all the input, I could only mark one as the correct answer, and i considered the fact that my original error was because of the lack of a `while` and improper use of `row` in the clip function, although the makefeaturelayer method is definitely valid!


title: Poor Temporal Raster Performance over REST

text: When trying to consume temporal rasters (float values) served from ArcGIS Server over a REST endpoint, the results served up are "wonky". One or two time steps will deliver data as expected but in general, it crashes and returns missing data or inverted values, etc. If I consume the data in ArcMap directly from SDE, it behaves correctly. If I consume the data in ArcMap via REST, it is a crap-shoot on what I'll get back, sometimes missing data. The data is symbolized using "Classified" with the data broken into 7 categories with zero values excluded. I am currently in the process of converting to integer-based rasters and also reclassifying to see if the issue is corrected but the end-goal is to put this into a web-viewer that I can also query against to get depth values at the x/y click. The server environment is as follows: **Two Separate Servers**...   * Windows 2008 R2 - ArcGIS Server 10.1 SP1   * Windows 2008 R2 - ArcSDE 10.1 SP1 with Microsoft SQL Server 2012 R2 PS: If I consume the raster using stretched symbology, it returns the data as expected although I can't control the color returned by the server after posting to AGS. I've tried using FGDB's and registered folders but achieved the exact same issues with that setup. Thoughts?


title: Error assessment of calculation done using unprojected vs. projected data

text: This question builds off of the one with the subject line "Calculating Flow Direction and Delineating Basins from Projected vs. Unprojected Data.": Calculating Flow Direction and Delineating Basins from Projected vs. Unprojected DEM Data This is an entirely separate question, however, as the aforementioned question has established that there are problems with using algorithms (e.g., ArcGIS Flow Direction) that assume Euclidean distance on data in a spherical/unprpojected geographic coordinate system. We know that map projections are sort of like taking an orange peel and attempting to flatten it out on a desk - you will have some error inherently introduced by the map projection. But, it seems that the benefits of projecting offset any error introduced, particularly when you're running calculations that assume a Cartesian/projected planar surface. In this case, the algorithm I'm interested in is the ArcGIS Flow Direction algorithm which does assume that your data are projected (and this is the assumption taken by most applications based on my research) as it uses a Euclidean approach for calculating distance. **My question is** : how could one quantify the error that might be introduced with calculating flow direction in a given study area using unprojected DEM data (DEM data in a geographic coordinate system) vs. projected data (DEM data in an appropriate projection such as a UTM or something conformal)? Granted, you can derive a flow direction raster using unprojected and then the same DEM data projected. But what then? Since our goal is to model the earth's surface as accurately as we can (and we're not addressing any errors that may be introduced in the process of creating the original DEM etc. - those are a constant as far as I'm concerned)....do we just assume the flow direction data derived from the projected DEM is better, and then compare the individual cell values of the two rasters to identify which cells have different directional values (in the context of the normal D-8 model)? I guess to do this then you'd have to take the flow direction raster derived from unprojected data, and then apply the same projection used with the projected flow direction raster. What would make the most sense, and what should the unprojected DEM be compared to as a benchmark of accuracy? Getting into the nitty details of the mathematical equations might, to those that understand it, give you proof at the ground level and be enough for some, but that as well as something that could convey the error to someone that does not have an in-depth understanding of the math but may just know enough geography/GIS to be dangerous would be great (ideally both levels would be good which would resonate with the hardcore geography geeks and the average GIS dabbler). For the higher level folks, saying that the proof is in the math possibly leaves it somewhat open for argument - I'm looking for something more tangible (e.g., akin to attaching a dollar figure to some sort of inefficiency in government). Any thoughts or ideas on how one might quantify this would be greatly appreciated. Tom


title: Cannot get areas for polygon layers in ArcMap 10

text: For an ArcGIS project opened in ArcMap 10, we have added polygon layers representing different geographical features. After all the polygons were created, we opened the attribute tables for those layers, and to our shock/dismay, the cells in the "shape_area" column all showed zero! (The "shape_length" field was correctly filled.) We also tried the Measure tool in the toolbar, but it could only measure distances, and the Area button was greyed out. We checked to make sure that all layers were in the correct geographical coordinate system (WGS 1984). What did we do wrong? Is there any way to obtain the areas of the polygons we created? Thanks! P.S. As Brad suggested, I tried to right click on the "SHAPE_Area" column (both on the column label and in the column), but the "Calculate geometry" (and "Field calculator") options were also greyed out.


title: Options for displaying PostGIS rasters in OpenLayers

text: This question is similar to an earlier question about displaying vector data in a web-map. I want to have a basic web interface which uses OpenLayers and can display raster data that is currently in PostGIS, using a basemap such as google. I have installed PostGIS2.0 with the raster library and see from the documentation that there are several options to output data, such as JPEG, GeoTIFF and PNG. Is their a recommend format to use? To display vector data I am currently using javascript to send coordinates from the user's browser to my server, and then performing a query using PHP/SQL on PostGIS. The results are then returned as GeoJSON which can be overlaid on the basemap. I wondered if a similar approach could be applied using rasters but I don't know what raster format is best, or if there are any other constraints I should consider. I have seen a lot of references to GeoServer but I'm not sure if it is necessary for this project (I'm not entirely sure when GeoServer should be used or not).


title: Reading palsar k&c 50 data in GDAL or other open source application

text: I posted this inquiry in the grass-user list. JAXA is providing free access to PALSAR 50m Orthorectified K & C Mosaic. According to the website, data are provided in: > Data files (Raw format) and image files (PNG format) However, I can't get GDAL to access the data (Raw format). I've tested other sample PALSAR data and GDAL can easily access them using the PALSAR driver I'm looking for open source apps that can read and export this data to the more common raster formats like geotiff, etc. The header file is described in this pdf: http://www.eorc.jaxa.jp/ALOS/en/kc_mosaic/map/PALSAR_MSC_dsc.pdf **UPDATE:** I found a non-foss application that can directly read the kc_mosaic format. It is called MapReady developed by the Alaska Satellite Facility. It can convert to various image format including geotiff. Another tool is the Next ESA SAR Toolbox (NEST). However testing with NEST4C-1.1, the data format is not recognized.


title: Server side clusterer for Google Maps

text: After doing a lot of research on different types of clustering coordinates (server side) I am still having problems with choosing the best approach for my project. My requirements:   1. Ability to work with more than 1.000,000 coordinates.   2. Be able to filter coordinates by point of interest.   3. Support map zooming and dragging.    4. Fast   5. I can't use any third party services. Here is what I found:   1. Region quadtree seems to the most suitable algorithm.   2. Geo hashing coordinates + Solr for quick retrieval/filtering of points (might only work with small set of data since the clustering will have to happen on the fly) I would like to know how to deal with map zooming & dragging while maintaining fast response from the server. How can clusters be cached if the maps is dragged, zoomed? Some clusters can be pre-clustered for large areas (continents) but what if there are 10,000 points within once city? My software stack is postgresql, python, django.


title: Correct network analysis settings needed for shortest route problem

text: Ok 1st sorry my English is not very good. What i want to ask is why the route analysis takes long path to reach the point when i click at solve button even though I can see a shorter path with my naked eye? I'm trying to digitize again and make a new network dataset, but it's still same like before. Am I missing any steps during creating the network?


title: How to create SDE GDB without post installation?

text: Is there any way to create a SDE sqlserver2008 geodatabase without :   1. SDE post-installation execution on the server   2. having sde service running on the server and using :   1. execution of some DDL scripts to create database on the SQLServer   2. execution of "sdelayer -o register" command to register some tables (containing Geometry type fields) as feature class.


title: How to create adjacency matrix of topology graph of arc map for shortest path?

text: I have created a topology graph of a polygon feature and I want to find the individual FIDS of the lines of the polygon. I want to construct an adjacency matrix to create the shortest path. How can I do this with ArcObjects programming with VBA?


title: call Sql Server stored procedure form arcpy

text: Is there a way to execute a SQL Server stored procedure with arcpy or Python and retrieve the results returned by the stored procedure? There is a very similar question posted on ArcGIS forum with no answer so far.


title: How to export Raster to JPEG in Python?

text: I built a model to import table data, create xy event, kernel density and clip density raster. I would like to export the clipped raster as a JPEG but I am not sure how to write the ExportToJPEG code. Any help would be greatly appreciated. Below is the model exported as a Python script.                   # Import arcpy module     ... import arcpy     ...           ... # Check out any necessary licenses     ... arcpy.CheckOutExtension("spatial")          ... # Local variables:     ... Input_Table = "C:\\Users\\cbriglia\\Desktop\\GIS\\Worksheets\\accidents_jan_2013.xls\\Layer1$"     ... state_border_shp = "C:\\Users\\cbriglia\\Desktop\\GIS\\Shapefiles\\state_border.shp"     ... Layer1_pytest4 = "Layer1_pytest4"     ... Ouput_Feature_Class = "C:\\Users\\cbriglia\\Desktop\\GIS\\Shapefiles\\layer1_pytest4.shp"     ... Hot_Spot_Output = "C:\\Users\\cbriglia\\Desktop\\GIS\\Shapefiles\\pytest5"     ... Clip_Output = "C:\\Users\\cbriglia\\Desktop\\GIS\\Shapefiles\\pytest6"     ...           ... # Process: Make XY Event Layer     ... arcpy.MakeXYEventLayer_management(Input_Table, "longitude", "latitude", Layer1_pytest4, "GEOGCS['GCS_North_American_1983_CSRS',DATUM['D_North_American_1983_CSRS',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]];-400 -400 1000000000;-100000 10000;-100000 10000;8.98315284119521E-09;0.001;0.001;IsHighPrecision", "")     ...           ... # Process: Copy Features     ... arcpy.CopyFeatures_management(Layer1_pytest4, Ouput_Feature_Class, "", "0", "0", "0")     ...           ... # Process: Kernel Density     ... arcpy.gp.KernelDensity_sa(Ouput_Feature_Class, "NONE", Hot_Spot_Output, "4.10886000000002E-03", "3.42405000000002E-02", "SQUARE_MAP_UNITS")     ...           ... # Process: Clip     ... arcpy.Clip_management(Hot_Spot_Output, "-73.7275781441499 40.9867050023656 -71.7867901391991 42.0506901260974", Clip_Output, state_border_shp, "-3.402823e+038", "NONE")


title: Simplify by number of points in QGIS

text: I have a fairly simply question about something I just don't know how to do. I am aware you are able to simplify polygons in QGIS based on distance tolerance but my questions is, is there a way to simplify based on a set maximum number of desired points? For example, if the number 1000 was entered, the polygon would be trimmed to 1000 points regardless of any distance thresholds. Does anyone know of any ways to do this?


title: Why does export to AI rasterize all my vector layers?

text: I have no raster layers in my map document and no transparencies. Yet the map is rasterized when exported to AI. (Map frame and coordinate grid are vector graphics but everything in the map itself is a single raster image) I've exported maps to AI before and I did not have this problem. This time the map is quite large (ANSI E @ 34"x44") Perhaps that is the problem? I know that all layers below a raster layer get rasterized in the output, any transparent layers and layers below a transparent layers. What are the other factors that will cause the AI exported to rasterize all the layers?


title: PRJ file for Stereo70 S-42 projection

text: I receive lots of Romanian maps with the projection STEREO70 with Datum "S-42 for ROMANIA". I have only found the "Stereo 1970.prj" by ESRI, and when I open the maps with my GIS software the coordinates do not match. As you can see, the PRJ by ESRI does not work with "S-42", but with "Dealul Piscului 1970": `PROJCS["Stereo_70",GEOGCS["GCS_Dealul_Piscului_1970",DATUM["D_Dealul_Piscului_1970",SPHEROID["Krasovsky_1940",6378245.0,298.3]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Double_Stereographic"],PARAMETER["False_Easting",500000.0],PARAMETER["False_Northing",500000.0],PARAMETER["Central_Meridian",25.0],PARAMETER["Scale_Factor",0.99975],PARAMETER["Latitude_Of_Origin",46.0],UNIT["Meter",1.0],AUTHORITY["EPSG",31700]]` I don't know if anybody is working with Stereo70 + S42 projection and can share their PRJ. Thank you very much.


title: Python script - how to avoid exclusive schema lock error?

text: I have a script for updating some features on my database every night (just for copying and replacing some features). This features are "read-only". My problem is that I can't avoid that these features are opened by users, and my script can show the following error:               ExecuteError: ERROR 000464: Cannot get exclusive schema lock.                    Either being edited or in use by another application.      Can I force the phyton script, via some command, copying the files, even been opened by some user? Can I fall down all the connections on my database before running the script? Thanks!


title: Topographic correction (for land use classification) with SAGA or GRASS

text: Is there a module in GRASS or SAGA (or a alternative open program) to do a topographic correction of a landsat image to enhance the land use classification? Thanks! Martin


title: Convert 32 bit image to 8 bit image in QGIS

text: I want to convert my images from 32 bits to 8 bits. I'm using otb to calculate texture and the raster output is on 8 bits, I'm going to use this information as a layer stack for my multispectral image at 32 bits! Please tell me how to convert from 8 bits to 32 bits and if it possible in QGIS! Thanks!