title: How to create a turn restriction in ArcObjects

text: I am having difficulties creating turn restrictions in ArcObjects. I have found a way to Add a a network turn to restrictions through the `INetworkForwardStarEx.AddTurnRestriction()` interface. But how to I create the INetworkTurn that this method takes.


title: Fiona error on installation

text: Receiving the following error when installing Fiona.Does anyone can guide ? error: Unable to find vcvarsall.bat


title: Merge Zip Codes to Create new sales Markets

text: Let me preface this question by saying I am very novice with respect to GIS. With this in mind, I am looking to merge zip codes into sales territories. Instead of plotting data by zip code, I want to have a new file that lets me plot data for the respective territories. Currently, I get around this by plotting the same value for all zip codes in the territory. What is the best/easiest way to do accomplish this task? I am quite familiar with the programming language `R`, but usually can work through examples if this requires programming in other languages. Many thanks in advance


title: How to convert/typecast from Polygon to Multipolygon in PostGIS?

text: Is there any way to typecast Polygons to Multipolygon type? For example, I have created a circle using ST_Buffer() on a point. The return type of that function is of Polygon type. The problem is that on the geometry column there is a restriction limiting it to Multipolygon type. So, I cannot insert a simple polygon. I have used ST_MPolyFromText() function, but as usual it is returning null for that and ST_Dump() is not giving the proper output. The constraint is:               CONSTRAINT enforce_geotype_geom CHECK (geometrytype(geom) = 'MULTIPOLYGON'::text OR geom IS NULL),


title: Search Raster Files from Folder that Have Same Spatial Coordinates Using Python

text: I have approximately 700 raster files in a single folder and I need to find files that have the same spatial extent and copy to a different folder. The raster files measure daily snow cover and therefore have the exact same spatial extent.


title: Can ArcGIS Desktop 9.3.1 open Geodatabase created using ArcGIS for Server 10.1?

text: can it possible to use old version of arcgIS desktop, say ver 9.3.1, if we develop a geodatabase using arcGIS 10.1. The above answers were very useful will make me to understand more on arcgis server, since we are planning Thanks _Suplimentary Question:_ may be question is not understood well..we have old version of arcgis desktop 9.3.1, now we proposing and wanted to develop geodatabase system using new version of arcgis server i.e 10.1, will it compatible.


title: How do I control what features get rendered with an ArcGIS Server service?

text: **Question:** Using a mapping service set up with ArcGIS Server, how can I set up a parameter that limits which features are being rendered? e.g. "/service?organization=megacorp"? **Background:** I am working on a web mapping application which will be used by organizations. They all have their own data, but it is stored in the same tables. When they log in, they will be brought to a map which shows their data points. Due to existing features, our preference is to store this data in a single table:               points_table     id | organization | name            ------------------------------     1  | the state    | poor house      2  | megacorp     | tenements       3  | union carbide| bhopal      It looks like the proper way to secure this service will be to use windows integrated authentication, and although I'm not clear on the details, use that to control which features are being rendered.


title: How to perform Ordinary Kriging in SAGA via QGIS SEXTANTE Plug-in?

text: I'm trying to learn how to apply the kriging analysis in QGIS using the SEXTANTE plug-in and SAGA tools. It appears that I have my SEXTANTE and SAGA folders configured properly according to the user guide. I've never used a kriging function before and I can't find a good resource from SAGA explaining the different input parameters. Does anyone have experience using the ordinary kriging tool in SAGA that can provide advice or a reference? I'm currently getting the following error message in QGIS:               An error has occured while executing Python code:      Traceback (most recent call last):        File "C:/Documents and Settings/PS/.qgis/python/plugins\sextante\gui\ParametersDialog.py", line 137, in accept if self.setParamValues():        File "C:/Documents and Settings/PS/.qgis/python/plugins\sextante\gui\ParametersDialog.py", line 91, in setParamValue if not self.setParamValue(param, self.paramTable.valueItems[param.name]):        File "C:/Documents and Settings/PS/.qgis/python/plugins\sextante\gui\ParametersDialog.py", line 129, in setParamValue return param.setValue(widget.getValue())        File "C:/Documents and Settings/PS/.qgis/python/plugins\sextante\parameters\ParameterExtent.py", line 11, in setValue if value is None: NameError: global name 'value' is not defined


title: Can't import OSM data to PostGIS, permission denied

text: I use Postgres 9.1 on Xubuntu 11.10. I've created a db 'aus' and user 'aus', granted ALL on the db to the user. Now I try to import the data, and it outputs this error:               $ osm2pgsql --create --database aus --username aus --prefix planet australia.osm -S /usr/share/osm2pgsql/default.style     osm2pgsql SVN version 0.70.5          Using projection SRS 900913 (Spherical Mercator)     Setting up table: planet_point     NOTICE:  table "planet_point" does not exist, skipping     NOTICE:  table "planet_point_tmp" does not exist, skipping     SELECT AddGeometryColumn('planet_point', 'way', 900913, 'POINT', 2 );     failed: ERROR:  permission denied for relation geometry_columns     CONTEXT:  SQL statement "DELETE FROM geometry_columns WHERE         f_table_catalog = '' AND f_table_schema = 'public' AND f_table_name = 'planet_point' AND f_geometry_column = 'way'"     PL/pgSQL function "addgeometrycolumn" line 133 at EXECUTE statement     SQL statement "SELECT AddGeometryColumn('','',$1,$2,$3,$4,$5)"     PL/pgSQL function "addgeometrycolumn" line 5 at SQL statement      What's wrong with my database? **edit:** the solution:               $ sudo su -l postgres     $ psql -d aus     # grant all on database aus to aus;     # alter database aus owner to aus;     ^D     ^D     $ osm2pgsql ...     here we go...


title: Finding Two locations with a conditional search on DEM

text: I want to find two locations H difference in elevation between two locations on a DEM. Ex: point A and B are my output after search on DEM. conditions are (1)buffer R distance around point A and (2) difference of H elevation to point A I am very poor in programing and looking for a tool or progra/script to achieve this. Your expert suggestions are really appreciated and useful. Thanks for your support. Eagerly waiting for your inputs. Thanks, Prasad.


title: Why does a layer imported from GPS have lat/lng rather than project coordinates?

text: I am importing tracks from my GPS using the GPS Tool in QGIS. I check the layer CRS and it shows WGS 84 but when I try and view the layer in my project which has a CRS of EPSG:2193 the coordinates are still lat and log (i.e. native WGS84). I have the same data on another computer and that works fine. I have compared all the setting that I can think of (Project Properties and Setting->Options) between the two system and all the ones I looked at are the same so far as I can tell. In particular I have made sure that "settings->options->CRS->'on the fly projection' is checked. I am pretty sure that I have had this issue before but can't remember how I solved it. What setting do I need to get QGIS to automatically convert the wgs84 layer to project projection? Alternatively what setting might I have inadvertently turned on that might have broken this? later: I worked through the various things that people suggested (thanks very for the suggestions BTW) I and still could not figure out what the problem was so I then started from scratch again and now everything works as expected. Sigh... I have been in IT long enough to know that when things just don't work as they should cycling power (or the software equivalent there of) should be tried before shouting "HELP!" not sure of the etiquette here -- should I mark this as answered?


title: How do I use the SearchCursor where_clause to get a max value from a column?

text: This is my first time asking a question on StackExchange. I hope I phrase my question and describe my attempts to understand my problem well enough! While working with SearchCursor, I'm trying to specify the SQL where_clause from a table called "Parcel" within the featureClass, and I want the MAX and MIN value from the column LAND_VALUE. Below is what I try just for MAX.               import arcpy          feature_class = "C:\foo.ext"     where_clause = "SELECT MAX(LAND_VALUE) AS MaxLandVal from Parcel"          rows = arcpy.SearchCursor(feature_class, where_clause, "", "", "LAND_VALUE D")      I'm told an invalid SQL statement was used. However, I'm not sure if the SQL statement is invalid, or if perhaps I'm stating the SQL query within SearchCursor incorrectly. Despite using the reference material for SearchCursor, SQL MAX() and SQL MIN(), and despite referencing a few different examples, I can't get this to work! Please advise! Finally, how would I display MaxLandVal if the query were to be correct? Perhaps something like...?               for blah in rows:         print blah.MaxLandVal      Thank you for your help!


title: OpenLayers LatLon conversion

text: I have an OpenLayers map with a settings like a below. All LatLon points are represented like a `POINT(3684182 5207957)`. In a more usual form (for me) it looks like a `(104.287422, 52.284126)`. How can I "convert" "POINT(…)" to this form?               projection: "EPSG:4326"     actualBounds: "0,0,4938083,5375780"     bounds: "0,0,9248768,9248768"     initCenterCoord: "3791668,5149639"     initCenterZoom: "5"     maxExtent:         bottom: -90         left: -180         right: 180         top: 90     maxResolution: 36128     numZoomLevels: 10     resolution: 141.125     restrictedExtent:         bottom: 1185098         left: -4938083         right: 9876166         top: 11936658     utm:         cityy: 9248768         hemisphere: "N"         inverseflattering: 298.2571644         majorsemiaxis: 6378137         offsetx: -1902571         offsety: -2422059         scalex: 100         scaley: -100         shiftx: -41454992         shifty: 583328750         zone: 48      Update: Some JS examples about a conversion between a projections:               > var projWGS84 = new OpenLayers.Projection("EPSG:4326");     > var proj900913 = new OpenLayers.Projection("EPSG:900913");     > var point = new OpenLayers.LonLat(3684182, 5207957);     c       lat: 5207957       lon: 3684182     > point.transform(proj900913, projWGS84);     c       lat: 5207957       lon: 3684182     > point.transform(projWGS84, proj900913);     c       lat: 5207957       lon: 3684182      Update 2: PostGIS conversion from 900913 to 4326 do it wrong:               SELECT ST_AsText( ST_Transform( ST_GeometryFromText('POINT(3684182 5207957)', 900913) , 4326 ) );     ------------------------------------------      POINT(33.0955700007803 42.3128389015332)      Screeshot: ![enter image description here](http://i.stack.imgur.com/soZMm.png)


title: How to report issues with QGIS documentation?

text: I stumbled across a gotach while installing mod_fastcgi on MacOs 10.8 (mountain Lion) -- there is a bug in 10.8 in the apache build env which breaks the install instructions here: http://hub.qgis.org/projects/quantum- gis/wiki/QGIS_Server_Tutorial The fix is to put a symlink                sudo ln -s /Applications/Xcode.app/Contents/Developer/Toolchains/XcodeDefault.xctoolchain/ /Applications/Xcode.app/Contents/Developer/Toolchains/OSX10.8.xctoolchain      as documented here: http://apple.stackexchange.com/questions/58186/how-to-compile-mod-wsgi-mod- fastcgi-etc-on-mountain-lion-by-fixing-apxserror It would be good to get this noted in the docs... Should I open standard bug report or is there a better way to report 'problems' with a document?


title: Openlayers & WMS - GetFeatureInfo problem

text:
I'm making a small project but I encounter some problems. Let me explain that : Having generated mapfile with mapserver and they are viewable by the openlayers interface, I would like to make me GetFeatureInfo on features. Subsite only problem is that I can not select the correct coordinates of a feature x and y (i and j for the url) because the x and y coordinates correspond to the "frame" openlayers and not of the layer displayed (problem with zooms, good feature pointing, etc.).               http://xxx/cgi-bin/mapserv.exe?map=C:/ms4w/Apache/htdocs/mapfiles/picc.map&SERVICE=WMS&Version=1.3.0&REQUEST=GetFeatureInfo&bbox=217711.311089,129521.320380,230314.065583,139320.681620&CRS=EPSG:31370&query_layers=arc&layers=arc,polygone&info_format=mime&i=913&j=162&WIDTH=1000&height=1000      img : http://i.stack.imgur.com/vNYm4.png **Do you have any idea to solve my problem?** Attached my html and mapfile code:               <!DOCTYPE html>               <meta charset="utf-8">     <meta name="WMS" content="WMS">     <meta name="MA2 Géomatique">     <meta name="Copyright" content="Copyright (c) 2013 MA2 Géomatique, Ulg">     <title>WMS </title>          <link rel="stylesheet" href="xxx/OpenLayers-2.12/OpenLayers-2.12/theme/default/style.css" type="text/css">     <link rel="stylesheet" href="style.css" type="text/css">          <link rel="stylesheet" type="text/css" href="xxx/OpenLayers-2.12/GeoExt/resources/css/popup.css">
WMS

PICC

URL UTILES :

             URL pour accéder au WMS dans une application client SIG :	             xxx/cgi-bin/mapserv.exe?map=C:/ms4w/Apache/htdocs/mapfiles/picc.map

                 URL pour accéder aux métadonnées du MapFile du PICC - routes et bâtiments (XML) :	             xxx/cgi-bin/mapserv.exe?map=C:/ms4w/Apache/htdocs/mapfiles/picc.map&SERVICE=WMS&REQUEST=GetCapabilities

                 URL pour accéder aux métadonnées du MapFile du MNT (XML) :	             xxx/cgi-bin/mapserv.exe?map=C:/ms4w/Apache/htdocs/mapfiles/mnt.map&SERVICE=WMS&REQUEST=GetCapabilities

xxx

          </meta>           Mapfile :               OUTPUTFORMAT     NAME agg     DRIVER AGG/PNG     IMAGEMODE RGB     END     # Legend     LEGEND       IMAGECOLOR 255 255 255     STATUS ON     KEYSIZE 18 12     LABEL       TYPE BITMAP       SIZE MEDIUM       COLOR 0 0 89     END     END          # Web interface definition. Only the template parameter     # is required to display a map. See MapServer documentation     WEB     # Set IMAGEPATH to the path where MapServer should     # write its output.     IMAGEPATH '/tmp/'          # Set IMAGEURL to the url that points to IMAGEPATH     # as defined in your web server configuration     IMAGEURL '/tmp/'          # WMS server settings     METADATA       'ows_title'           'WMS-PICC'       'ows_onlineresource'  'http://xx/cgi-bin/mapserv.exe?map=C:/ms4w/Apache/htdocs/mapfiles/picc.map'       'ows_srs'             'EPSG:31370'       "wms_enable_request" "*"       "wms_include_items" "all"     END          #Scale range at which web interface will operate     # Template and header/footer settings     # Only the template parameter is required to display a map. See MapServer documentation     TEMPLATE 'fooOnlyForWMSGetFeatureInfo'     END          LAYER     NAME 'arc'     TYPE LINE     DUMP true     TEMPLATE fooOnlyForWMSGetFeatureInfo     EXTENT 217711.311089 129521.320380 230314.065583 139320.681620     CONNECTIONTYPE postgis     CONNECTION "dbname='projetsig' host=xx port=5432 user='xx' password='xx' sslmode=disable"     DATA 'geom FROM "arc" USING UNIQUE gid USING srid=31370'     FILTER ( "rue" IS NOT NULL )     METADATA       'ows_title' 'arc'       "wms_enable_request" "*"       "wms_inlude_items" "all"       "ows_featureinfo_textplain_include_items" "*"     END     STATUS OFF     TRANSPARENCY 100     PROJECTION     'proj=lcc'     'lat_1=51.16666723333333'     'lat_2=49.8333339'     'lat_0=90'     'lon_0=4.367486666666666'     'x_0=150000.013'     'y_0=5400088.438'     'ellps=intl'     'towgs84=-106.869,52.2978,-103.724,0.3366,-0.457,1.8422,-1.2747'     'units=m'     'no_defs'     END     CLASS        NAME 'arc'         STYLE          WIDTH 0.91           COLOR 54 54 54        END     END     END          LAYER     NAME 'polygone'     TYPE POLYGON     DUMP true     TEMPLATE fooOnlyForWMSGetFeatureInfo     EXTENT 217711.311089 129521.320380 230314.065583 139320.681620     CONNECTIONTYPE postgis     CONNECTION "dbname='projetsig' host=xx port=5432 user='xxx' password='xxx' sslmode=disable"     DATA 'geom FROM "polygone" USING UNIQUE gid USING srid=31370'     METADATA       'ows_title' 'polygone'       "wms_enable_request" "*"       "wms_include_items" "all"     END     STATUS OFF     TRANSPARENCY 100     PROJECTION     'proj=lcc'     'lat_1=51.16666723333333'     'lat_2=49.8333339'     'lat_0=90'     'lon_0=4.367486666666666'     'x_0=150000.013'     'y_0=5400088.438'     'ellps=intl'     'towgs84=-106.869,52.2978,-103.724,0.3366,-0.457,1.8422,-1.2747'     'units=m'     'no_defs'     END     CLASS        NAME 'polygone'         STYLE          WIDTH 0.91           OUTLINECOLOR 0 0 0          COLOR 109 188 214        END     END     END          END      Thx.


title: Why is "Convert Graphics toe Features" grayed out in a dataframe that has georeferenced layers?

text: I have brought two graphics, a species range and a heatmap, into a dataframe and georeferenced them using control points. I rectified them and now I want to clip one with the other but I can't get them converted to features because that function is grayed out. So what I think is a graphic is not? When I open the clip tool it doesn't show these layers because, I assume, ArcMap thinks they are not features? Thanks, Steve


title: Can't select layer to upload to GPS

text: I'm trying to upload a vector layer to my Garmin GPS unit. I've saved it as GPX (WGS 84); opened the GPS Tools dialogue in QGIS 1.8; selected the "Upload to GPS" tab; and clicked on the "Data layer" dropdown arrow. However, the box just turns blue; nothing appears for me to navigate through. Nor can I enter anything into the box. It just turns blue, but will not take any input. Any idea how I select the layer to upload? ![enter image description here](http://i.stack.imgur.com/2Nt5x.jpg)


title: How to show attribute table in the form

text: Is there any way to show the attribute table of a layer in a QtextBrowser or a QtableView on form? I did this but,(shows attribute table in a separate dialog) I have no idea, how to show it in my form. iface.showAttributeTable(iface.activeLayer) Thank You


title: Python GDAL/OGR: check lines geometries connectivity

text: I'm working with shapefiles representing roads, so basically it's stored as (poly)lines.   I need to perform a road connectivity checking. In other words, I want to ensure that roads meant to communicate are really in contact, because I don't wanna see a little gap between my 2 segments of road.   I don't really see how to proceed. Should I compute the distance from the current line to the closest ones and then look if this distance is under a threshold? Or is there an easier solution?   Thanks,   eo


title: Geoserver SQL view parameterizing

text: I've got a postgresql table that has the following attributes:               name        String     begintime   String     endtime     String     the_geom    Point      Now I'd like to create an sql view in geoserver (2.2 snapshot). In short, i only want to map the features where begintime equals a certain param. Looking at the geoserver documentation and the questions here at the forum, I came up with this:               select name, begintime, endtime, the_geom      from n80_tellingen_merge      where begintime = %utct%      This gets me the following error:               ERROR: syntax error at or near ")" Position: 111      So basically I'm asking two things:   * What am I doing wrong?   * Since I don't have brackets or more than 111 positions, where can I find the actual query sent to postgresql? Can I consult a log somewhere? (Ubuntu server 13.04) My apolagies if I'm missing something obvious.


title: Why does map display like it is in Isometric when everything is in EPSG:4326?

text: I'm beginning with Geoserver and try to display a map of my country (Belgium) is an openlayer map (embedded into a MapPanel object). I've downloaded the municipalities limits from gadm.org and copied the file into a POSTGIS database (using the tools provided in the Geoserver package). I guess it is a projection issue, but I don't get why since everything is in EPSG:4326. Thanks in advance! Here is my code:               panel = new GeoExt.MapPanel({        title: "MapPanel",        renderTo: "map-id",        layers: [             new OpenLayers.Layer.WMS("Administrative borders",                 "/geoserver/wms",                 {layers: "bel_adm4", projection:"EPSG:4326", format: "image/png"},                 {transitionEffect:'resize', isBaseLayer:true}        )],        map:           {             projection:"EPSG:4326",              maxExtent: new OpenLayers.Bounds(2.5412,49.5100,6.4013,51.502)            }      });      Here's a screenshot of the generated map: ![enter image description here](http://i.stack.imgur.com/VkpRt.png) Here's how Belgium should be: ![enter image description here](http://i.stack.imgur.com/dTXOl.gif)


title: Polygon to fill gap

text: Using postgis I want to fill a gap between other polygons. Can you help with some directions what functions to use ? Thank you and sorry if the question is rather noobish. Image to explain from where the to extract the poylygon: http://i.stack.imgur.com/PjcTl.png


title: Upload geotiff file using GeoServer Rest API

text: I can upload shapefile to GeoServer via REST api using cURL. I used this code:               curl -v -u admin:123 -XPUT -H "Content-type: application/zip" --data-binary `@F:\temporary\shp.zip` `http://192.168.79.131/geoserver/rest/workspaces/arc/datastores/testShape/file.shp`      Now I want to upload georeferenced **TIFF files** in same way. I zipped the tiff file and then, used the following command to upload:               curl -v -u admin:123 -XPUT -H "Content-type: application/zip" --data-binary  @F:/testdem.zip http://192.168.79.131/geoserver/rest/workspaces/nurc/coveragestores/testRaster/file.geotiff      But I got errors: > Unsupported format: geotiff* we are done reading and this is set to close, > stop send I tried other formats like: file.imagemosaics | file.tif | file.tiff | file.geotif, but nothing changed. I use GeoServer 2.3.3 on WinXP SP3.


title: How to install pgrouting 2.x on ubuntu 12.04 LTS

text: I try to install pgrouting 2.0 on ubuntu 12.04 LTS with postgresql 9.2 I have ppa in the system by `sudo add-apt-repository ppa:gwibber-daily/ppa` answer: `You are about to add the following PPA to your system: Daily builds of Gwibber trunk More info: https://launchpad.net/~gwibber-daily/+archive/ppa Press [ENTER] to continue or ctrl-c to cancel adding it gpg: zbiór kluczy ,,/tmp/tmpO_cWbg/secring.gpg'' został utworzony gpg: zbiór kluczy ,,/tmp/tmpO_cWbg/pubring.gpg'' został utworzony gpg: zapytanie o klucz 72D340A3 z hkp serwera keyserver.ubuntu.com gpg: /tmp/tmpO_cWbg/trustdb.gpg: baza zaufania utworzona gpg: klucz 72D340A3: zaimportowano klucz publiczny ,,Launchpad PPA for gwibber-daily'' gpg: Ogółem przetworzonych kluczy: 1 gpg: dołączono do zbioru: 1 (RSA: 1) OK` next I've typed `sudo apt-get update` looks like everything OK, but what next? on the topic i haven't found much informations. Thanks for answer.Tom


title: Plotting a map of New Zealand with regional boundaries in R

text: I have been asked to make a map of New Zealand using R _with regional boundaries_. I need something exactly like this example of a map for Switzerland: http://blog.revolutionanalytics.com/2009/10/geographic-maps- in-r.html I tried copying the same commands but after I downloaded the data from gadm and tried plotting it with spplot, I ended up with ten tiny plots with New Zealand on the extreme right hand edge. I can't paste the output into this question because R will not let me save it as a file. But what I did was simply download the "level 1" data from gadm in R format and then `spplot(gadm)`. Please could someone help me to make a map of New Zealand showing regional boundaries, so that I can make the regions different colours? It looks like the main problems at the moment are   1. I would like the North and South Islands to be in the middle.   2. I don't want a super-detailed map as I would like to plot it quickly.   3. I just want a single map, not ten of them. Alternatively, any pointers to tutorials would be very welcome. I don't necessarily want to use the gadm data if there is a better alternative.


title: How to install PostGIS correctly with OpenGeo Suite 4.0 community edition on Ubuntu?

text: I have installed the OpenGeo Suite community edition, version 3 and 4, onto Ubuntu Server 10.04 and 12.04 multiple times in the past, with virtually no problems. With the latest version of OpenGeo Suite (which includes PostGIS 2.1), I am able to install almost all of the components successfully, except PostGIS, on ubuntu 10.04 and 12.04. With earlier editions, PostGIS was installed as part of the full Opengeo installation process, but now, although PostGreSQL 9.3 is successfully installed, there are no PostGIS or opengeo database templates installed, nor are the spatial functions installed. According to apt-get, however, PostGIS 2.1 is installed. I have also attempted to install OpenGeo v4 community edition on a local ubuntu server on a vm. Same issue... geoserver, geowebcache, the opengeo dashboard, PostGres 9.3, and PostGIS 2.1 all installed, but no opengeo or postgis databases, and no functions. I'm wondering if this is an issue with the Opengeo installation package?


title: Set default raster output format in qgis

text: When I save a map as image, either from the main window or from a print composer in qgis (v 2.0.1) on debian, the default format is BW-files. Never heard about them before, never used them.. Is there any way to set another (i.e. .png) as the default image format? I have not been able to find anything in Settings / Options. Is it possible to set it directly in some configuration file? (I am using qgis both on linux and windows on the latter, the file format "sticks")


title: Find ordinal polygon points

text: How might I go about finding the following points of a `MULTIPOLYGON` in PostGIS?   * Northeasternmost   * Northwesternmost   * Southeasternmost   * Southwesternmost Can't seem to find a `ST_` function to accomplish this.


title: Label size in map units, geoserver .sld file

text: Is it possible to define labels in map units in a sld file? Thanks.


title: How to setup a debugger for QGIS Plugin development?

text: As the title of the question implies, I wish to know how to set up a QGIS Python Plugin Debugger so that I can develop plugins more efficiently. Any information regarding this would be much appreciated as I am fairly new to QGIS Python Plugin Development. Thanks


title: Is code completion available in QGIS Python Console?

text: I wanted to know if the autofill functionality is available in the python console for QGIS as the one available in ArcGIS Desktop's Python Console or am I supposed to hardcode it.


title: Is it possible to show changes in database from geowebcache?

text: Is it possible to show changes in database from geowebcache (used in conjunction with geoserver) which is showing cached data instead of reflecting changes?


title: OSGB conversion error

text: So I am using `ArcGIS 10` and trying to work in British National Grid. My shape file is in WGS84 UTM30N. By chance I passed the shape file across to a colleague and we both transformed into BNG. When we compared the data there was a significant difference in our results. I have since cross checked with QGIS, Global mapper and a website or two and got a whole range of results but in general the most common is the result my colleague got. My question I guess is three-fold:   1. What have I done wrong?   2. Why is there such a variation (~150m) in values? (1 or 2m error is probably ok for what I am looking at.)   3. Which value is "right"? Bit of a beginner-ish question I fear but help would be much appreciated!


title: Applying rotation before FeatureToRaster_conversion?

text: I am converting polygon shapefiles to rasters. With arcpy I am using the function `FeatureToRaster_conversion`.               gp.FeatureToRaster_conversion(InFeatures, InField, OutRaster, InCellSize)      This function gives control of the cell size but nothing concerning the rotation. I would like the output raster with a given orientation other than North. Say for example having the picture at 45 degrees. Any ideas how I can achieve this? With or without using FeatureToRaster_conversion.


title: Display a WFS layer with OpenLayers

text:
Good Morning! I have some troubles displaying a WFS layer with OpenLayers and hope someone here can give me a hint. Before moving on, I precise that I have already checked other posts, like here and here but, for any reason, it hasn't solved my problem yet (and by the way, thank you to the contributors to those posts, they made me understand a little bit better already)... The WFS layer appears in the layer switcher, but not at all on the map. So, here is the configuration of my data:   * Workspace: luma_project (this is the default workspace)   * Workspace URI: http://localhost:8080/luma   * URL for WFS: http://localhost:8080/geoserver/wfs   * Layer name: adresspt   * Srs: 4326 and here is my OpenLayers code:               var adresspt_wfs = new OpenLayers.Layer.Vector("Adress points WFS", {       strategies: [new OpenLayers.Strategy.BBOX()],         protocol: new OpenLayers.Protocol.WFS({         url: "http://localhost:8080/geoserver/wfs",         featurePrefix:"luma_project",         featureType: "adresspt",         featureNS: "http://localhost:8080/luma",         srsName: "EPSG:4326",         geometryName: "the_geom",         version: "1.0.0"       })     });     map.addLayer(adresspt_wfs);      Do you see anything missing here? As my data is stored in a (local) PostGIS database, should I refer to it as _workspace:store:layer_? or _workspace:layer_? Also, I precise that I have no trouble to display this layer as WMS in OpenLayers and that I can even display it as WFS in Quantum GIS. The problem thus doesn't come from the data or the Geoserver, but really from my OpenLayers code. Any help will be greatly appreciated! Thanks a lot in advance, Stéphane EDIT: Here is the full code:               var map;          function init() {     var bounds = new OpenLayers.Bounds(         13, 55.67, 13.5, 55.8       );           var options = {         projection: new OpenLayers.Projection("EPSG:4326"),         displayProjection: new OpenLayers.Projection("EPSG:4326"),         maxExtent: bounds,         units: 'degrees'       };            map = new OpenLayers.Map('map');            var transportation = new OpenLayers.Layer.WMS(         "Lund transportation network", "http://localhost:8080/geoserver/wms",         { layers: 'lu_transportation',           format: 'image/png',           srs:'EPSG:4326',           transparent:'true',           zoomOffset: 11,         },        {isBaseLayer: true, opacity: 0.4});            var adresspt_wfs = new OpenLayers.Layer.Vector("Adress points WFS", {         strategies: [new OpenLayers.Strategy.BBOX()],         protocol: new OpenLayers.Protocol.WFS({           url: "http://localhost:8080/geoserver/wfs",           //featurePrefix:"luma_project",           featureType: "luma_project:adresspt",           featureNS: "http://localhost:8080/luma",           srsName: "EPSG:4326",           geometryName: "the_geom",           version: "1.0.0"         })       });              map.addLayer(transportation);         map.addLayer(adresspt_wfs);            map.addControl(new OpenLayers.Control.Navigation());       map.addControl(new OpenLayers.Control.LayerSwitcher());       map.addControl(new OpenLayers.Control.Scale($('scale')));       map.addControl(new OpenLayers.Control.MousePosition({element: $('location')}));            map.zoomToExtent(bounds);          }      And then the html code:
         Coordinates of the mouse:
...

                  EDIT (2) Just to try out, I have added another WFS layer, from a remote server, and still nothing displays... I took the exact code of this page, which is:               var test_wfs = new OpenLayers.Layer.Vector("WFS", {       strategies: [new OpenLayers.Strategy.BBOX()],       protocol: new OpenLayers.Protocol.WFS({         url: "http://demo.opengeo.org/geoserver/wfs",         featureType: "tasmania_roads",         featureNS: "http://www.openplans.org/topp"       })     });      I really don't understand... this layer displays perfectly on this page, but not on mine... with the same code...


title: Why does ArcGIS Server not return data when querying joined fields?

text: I've been trying to help troubleshoot an issue that is very strange, and I think we've boiled it down to an issue with ArcGIS Server. We have a map service that shows up in the services directory as:               Layers:         Some-LayerA (0)         ...         Some-LayerE (4)          Tables:         databaseName.SchemaName.TableName (5)      And the layer we are querying against looks roughly like this:               Fields:         databaseName.SchemaName.LayerName.ObjectID (type: esriFieldTypeOID, alias: OBJECTID)         databaseName.SchemaName.LayerName.SomeField (type: esriFieldTypeString, alias: SomeField, length: 3)         databaseName.SchemaName.TableName.ProjectKey (type: esriFieldTypeString, alias: ProjectKey, length: 18)         databaseName.SchemaName.TableName.Spent To Date ( type: esriFieldTypeDouble , alias: Spent To Date )          ... More joined fields (No duplicate names or aliases, all strings, dates, doubles)      When issuing a query with only a where clause **against the joined table** , I can do something like this and have it bring results up:               Where: databaseName.SchemaName.TableName.ProjectKey = 'ValidProjectKey'          --> Returns 2 results     -----> databaseName.SchemaName.TableName.ProjectKey: ValidProjectKey          -----> databaseName.SchemaName.TableName.ProjectKey: ValidProjectKey      However, if I add some extra parameters, then I get no results at all:               Where: databaseName.SchemaName.TableName.ProjectKey = 'ValidProjectKey'     Out Fields: *          ---> Returns 0 results      Apparently, if specifying the out fields manually, it will give results:               Where: databaseName.SchemaName.TableName.ProjectKey = 'ValidProjectKey'     Out Fields: databaseName.SchemaName.TableName.ProjectKey          ---> Returns 2 results     -----> databaseName.SchemaName.TableName.ProjectKey: ValidProjectKey          -----> databaseName.SchemaName.TableName.ProjectKey: ValidProjectKey      But, if I change that second query to go against one of the **layer's native (non-joined) fields** instead, then it seems to return an appropriate number of results, but with all of the join fields set to null:               Where: databaseName.SchemaName.LayerName.SomeField = 'SomeValue'     Out Fields: *          ---> Returns 514 results     ------> databaseName.SchemaName.LayerName.ObjectID: 553     ------> databaseName.SchemaName.LayerName.SomeField: SomeValidValue     ------> databaseName.SchemaName.TableName.ProjectKey: null     ------> databaseName.SchemaName.TableName.ProjectManager: null          ------> databaseName.SchemaName.LayerName.ObjectID: 1938     ------> databaseName.SchemaName.LayerName.SomeField: SomeOtherValidValue     ------> databaseName.SchemaName.TableName.ProjectKey: null     ------> databaseName.SchemaName.TableName.ProjectManager: null      Any clues why going against the joined table would lead to such different behavior? I don't have much control over the actual query, because it is being issued through geocortex.


title: Understanding How ArcGIS Displays Non Spatial Data Created in Microsoft Access

text: I am working on getting American Community Survey(ACS) data into SDE table format. I have managed to import all the text files downloaded from ACS into excel and then imported to Access for further extracting. Anyway, so all the tables have an ID field. When I try to preview the data in ArcCatalog it would not show the data in proper order shown below. The problem is when I register this tables as a geodatabase tables it is not using an ID field to create OID field but it assigns OID to the order that it was displayed in preview tab. I want to **understand how it creates and assigns the OID** to each row. I know I can right click on the field and sort by ascending order and it would display how I want it but it is a temporary solution. This does not solve the problem when I convert the access table to a SDE table. ![enter image description here](http://i.stack.imgur.com/kFGJt.png) Thank you in advance.


title: What does MODIS data reflectance values (MOD13Q1) mean?

text: Several publications suggest filtering MODIS data (in my case MOD13Q1) prior to use. They recommended to filter out pixels having blue reflectance more than 10%, but what 10% here means? is it 10% of the total possible value for blue reflectance or 10% or the total reflected over the whole spectrum? Any advices?


title: QGIS Digitizing in Cartesian

text: I have a client, that needs some geological maps digitized (polygons) in a Cartesian grid. It would be a basic, square grid with North up, and coordinates in x and y. Some of the maps will need to be digitized in meters and some in feet. I've looked around a lot for a way to set this up in QGIS, but I can't find anything that I can understand. I think if my technical knowledge was better, I could make a custom projection that would work. I'm baffled that this answer isn't easier to find. Although I wouldn't be surprised if I'm overlooking something. I'm not interested in placing this spatially. It can stay in the Cartesian grid. I'm doing the digitizing in QGIS and he is placing the SHP files into MapInfo for plotting. Edit: A question that I've been wondering is if it is fine to digitize in a UTM grid as a SHP file, then just delete the PRJ? Would the SHP file then have the right coordinates in a flat plane? Or is there more to it than that? Thanks,


title: Connect not working for custom QGIS form

text: So I have a custom form for editting feature attributes. However, I can't seem to correctly disconnect and reconnect certain bindings. The init code is below and should work for any simple QT widget with buttons. It doesn't throw an error. The OK button is disconnected from the original connect that QGIS assigns. But it does not reconnect to the validate function. Can anyone point me in the direction of the mistake here?? Thank you in advance if you can!               from PyQt4.QtCore import *     from PyQt4.QtGui import *          # QGIS passes in the dialog, he layer_id and the feature_id     def form_open(myDialog,layer_id,feature_id):         # QGIS passes in the dialog, he layer_id and the feature_id         # I send all this to be handled by its own class. So I can use          #'self' and avoid having to use globals...         return MyClass(myDialog,layer_id,feature_id)          class MyClass(QDialog):         def __init__(self,myDialog,layer_id,feature_id):             self.ui = myDialog             self.buttonBox = self.ui.findChild(QDialogButtonBox,"buttonBox")                  # The problem is here...             self.buttonBox.accepted.disconnect(self.ui.accept)             self.buttonBox.accepted.connect(self.validate)              def validate(self):             msgBox = QMessageBox()             msgBox.setText("I validate here...")             msgBox.exec_()


title: Working with Interrupted Goode Homolosine data

text: I'm trying to use Global Land Cover Characterization data from the USGS: http://edc2.usgs.gov/glcc/fao/data_specifications.php It is in the Interrupted Goode Homolosine projection. As I am working on a continental scale (South America) I need to reproject to WGS84. I've tried using the Goode Homolosine projection (land) in ArcGIS but the resulting raster was offset to the East by around 11km. I shifted this using the header file of the original .bil file. I need to clip it down to the area of my study site and so used a snap[ raster of existing data (worldclim data). The snap did not work and my GLCC raster is still warped and does not line up with underlying rasters. I have seen the related questions on GIS stack exchange but none of those appear to work either. I'm running ArcGIS10 sp1. _**hacky fix_** Okay, here's what I did to get the forest cover Global Land Cover Characterization Data to line up with Worldclim data (for Latin America in WGS84 (geographic) projection). Took the global file from the USGS website, added the header to the same folder and then projected to Goode Homolosine (land) in Arc Catalog. Modified the Goode Homolosine projection to use the Sphere (Arc Info) as it's spheroid. Clipped the bit I needed and projected it to WGS84 (Geographic). I then resampled the GLCC data as it was at 1000m resolution whereas the worldclim stuff is in 30 arc seconds. Before resampling I checked all the environments and added a snap raster (the worldclim ascii I wanted it to match), extents, output coords, cell size, etc). They now line up. Hacky, as I said but it worked. And just for the record, I'm working with MaxEnt (hence the ASCII requirement). There is word on the wire that you need to turn off background geoprocessing in order to work around a bug in 10 with the snap raster option. Mine was off anyway. In addition, you can resample first and then add the snap raster to a "copy raster" command. Might be safer that doing it all in one.


title: How to create an accurate Tissot Indicatrix?

text: A Tissot Indicatrix is useful method for communicating at a glance the kinds of distortion a given projection is prone to (in the figure below, each of the red circles occupies the same area). I've been told that popular methods for generating TI's have their own problems, to the point of being sometimes woefully inaccurate. What is the problem with the popular methods, and what is most correct way of generating a TI that is accessible to your average GIS dude(ette)? ![Mercator and globes with tissot's](http://i.stack.imgur.com/r1TPf.jpg)


title: How to add a field to a .zip shapefile?

text: I am exploring the possibility of using ArcGIS Online (Free Public Account and Explorer) with my new Global Studies students. (So far, the Free Public Account appears to be more user friendly and stable.) I found a nice shapefile on Natural Earth that shows all the countries of the world. I would like to add some fields so that students could add data and manipulate the map. I want to create blank fields for students to enter data (Govt, GDPperC, Generic1, Generic2). I have tried taking the dbf file out of the zip and editing it in Access, but when I put it back in the zip, ArcGIS Online couldn't load it. I tried loading the .zip file as a layer in QGIS, but it won't let me edit anything. :( Is there a way to add the 4 fields (ideally at the beginning, so they don't get lost in the plethora of data fields in the shapefile)? I am not a programmer, so this has to happen via some GUI method.


title: How to georeference an image in ArcMap, using coordinates from a KML file

text: My starting point is KML file containing a PNG image, which displays correctly in Google Earth: ![enter image description here](http://i.stack.imgur.com/3l3OH.jpg) From Google Earth I choose `> Save Place As > KML`, which gives me the unprojected PNG file, plus the KML file telling Google Earth how to display the image. This includes the coordinates in the **LatLonbox** section:               <GroundOverlay>       <drawOrder>1</drawOrder>       <Icon>         <href>image.png</href>       </Icon>       <LatLonBox>         <north>-43.35635491009509900</north>         <south>-43.74776058465478500</south>         <east>172.75212961661239000</east>         <west>172.35586651444777000</west>       </LatLonBox>     </GroundOverlay>      The PNG file is 399px x 544px in size: ![enter image description here](http://i.stack.imgur.com/LFFQ9.png) ... so I input the coordinates manually in the ArcMap Georeferencing toolbar, specifying 0,0 for the top left of the image, and 399,-544 for the lower right of the image, matching these with the north/west and south/east coordinates from the KML file. The data frame is set to WGS84, which is apparently the CRS used by Google Earth. The resulting geocoded image does not correctly match real-world coordinates (the image should be slightly offset as shown by the red arrow - and possibly rotated so the coastline matches): ![enter image description here](http://i.stack.imgur.com/XBNKy.png) The image correctly displays in Google Earth, so the coordinate reference must be correct - presumably my logic of using that bounding box with the image bounding box is flawed. Any ideas what I'm doing wrong here? Thanks


title: What are the set of projections that can be accurately georeferenced using polynomial 2?

text: **Background** I'm starting with a map from the internet (Encyclopedia Britannica) and georeferencing it using QGIS 2.0.1, I get really accurate georeference solution using polynomial 2 transformation: very small residuals, most less than 5 pixels. (This is surprising because the previous time I tried this, I got pretty poor results: How can I improve my Georeferencing results? though now I am using more points (31) and was quite careful with them.) The georeference window: ![The fit](http://i.stack.imgur.com/JwoNY.jpg) (note the very slight residuals, the biggest of which is the small line emanating from Beijing). This produces the following overlay, which you can tell is quite good. ![Georeferenced map](http://i.stack.imgur.com/xWwJw.png) **Goal** I'd like to identify the projection and its parameters used by the cartographer. **Questions** Can the fact that a polynomial 2 fits really well be used to narrow down the list of projections that the map uses? Or do many/most projections georeference well with polynomial 2 when applied over such a large geographic extent and with so many control points? (I can disable most of the control points and as long as I include at least one in the northern regions, I can get the same very good georeference. And of course polynomial 3 works very well also.) (I've also tried to fit most of the pseudocylindrical projections mentioned in The Proj Manual (1996, revised 2003) to see if one can be made to match the control points as well as the polynomial 2 georeference fit, but the latter is much better than any of those fitted projections.) **The GCPs** I have put the GCPs at http://pastebin.com/pWti0W8k. The map x/y are in EPSG:4326 (I used the Google Physical OpenLayer map as the reference but then reprojected those EPSG:3857 to 4326 so I could get lon/lat).


title: resampled thematic map accuracy

text: I've generated a map of ecoregions based on some environmental data using a semi-supervised classification where I combined a statistical clustering with random forests to generate my map. I created the map at a 30m resolution and to remove some of the small scale variability I've re-sampled (with the majority filter) to 2.5km. I was thinking of using McNemars test to investigate for any significant differences between my original (30m) and re-sampled (2.5km) data, but my validation dataset becomes significantly clustered at less than 2.5km (according to Ripley's K). Also, my validation dataset is quite large (~2000 samples) and from what I've noticed so far is that the difference is coming back as significant (p<<0.001) using McNemars test even though % agreement is ~95% with a kappa of .96. I think it's due to the large sample size of my validation set. Would it be suitable to compare the differences in the map with something like a Mann-Whitney or Wilcoxon test on an independent random sample of pixels? i.e. create _x_ amount of sampling points, extract class info from the 30m and 2.5km maps and then compare? The other thing I've thought about doing, was to look at lumped class accuracy with McNemars test i.e. lump the class information from my validation sites on both maps and compare those counts across classes to the validation dataset. An issue I can see here though (apart from the clustering of my validation samples at less than 2.5km) is that there only needs to be a difference of 1 pixel between the validation set and my classification sets for it to counted as misclassified. Any help would be greatly appreciated.


title: Convert points to lines with Python GDAL

text: I have a dictionary with points. The points are derived from a raster (see image below). The points are not sorted and are not in order. The dictionary looks for example like this:               pointDict = {0: (345645.1276541934, 1267223.104499615), 1: (345626.87681620114, 1267223.2540966477), 2: (345645.2772512261, 1267268.581997563), 3: (345617.751397205, 1267223.4036936804), 4: (345654.1034761568, 1267259.306981534), 5: (345636.15183223, 1267231.781127513), 6: (345636.30142926273, 1267268.2828034975), 7: (345626.87681620114, 1267259.306981534), 8: (345617.90099423775, 1267259.306981534), 9: (345608.7755752416, 1267259.6061755994), 10: (345599.7997532782, 1267250.1815625378), 11: (345590.6743342821, 1267250.4807566034)}      ![enter image description here](http://i.stack.imgur.com/IMQ9T.png) I want to create a multiline. The maximum distance from point to point is 14 m. If the points are further apart, they go to a new line. So far I got the following code. It works, but the problem is, that the points are connected in the wrong order, as you can see in the image below.               import ogr, gdal, os     from math import sqrt          pointDict = {0: (345645.1276541934, 1267223.104499615), 1: (345626.87681620114, 1267223.2540966477), 2: (345645.2772512261, 1267268.581997563), 3: (345617.751397205, 1267223.4036936804), 4: (345654.1034761568, 1267259.306981534), 5: (345636.15183223, 1267231.781127513), 6: (345636.30142926273, 1267268.2828034975), 7: (345626.87681620114, 1267259.306981534), 8: (345617.90099423775, 1267259.306981534), 9: (345608.7755752416, 1267259.6061755994), 10: (345599.7997532782, 1267250.1815625378), 11: (345590.6743342821, 1267250.4807566034)}          multiline = ogr.Geometry(ogr.wkbMultiLineString)          i = 0     lineDict = {}     for item in pointDict:         stop = False              x = pointDict[item][0]         y = pointDict[item][1]              if item != 0:             xPrevious = pointDict[item-1][0]             yPrevious = pointDict[item-1][1]             distance = sqrt((y-yPrevious)**2+(x-xPrevious)**2)              for line in multiline:             if line.GetPointCount() > 0:                 j = 0                  for j in range(0, line.GetPointCount()):                     point = line.GetPoint(j)                     xExisting = point[0]                     yExisting  = point[1]                     distance = sqrt((y-yExisting)**2+(x-xExisting)**2)                     j += 1                     if distance < 14:                         line.AddPoint(x,y)                         stop = True              if not stop:             lineDict[i] = ogr.Geometry(ogr.wkbLineString)              lineDict[i].AddPoint(x,y)               multiline.AddGeometry(lineDict[i])             i += 1          for line in multiline:         print line          outSHPfn = 'test1.shp'        shpDriver = ogr.GetDriverByName("ESRI Shapefile")     if os.path.exists(outSHPfn):         shpDriver.DeleteDataSource(outSHPfn)     outDataSource = shpDriver.CreateDataSource(outSHPfn)     outLayer = outDataSource.CreateLayer(outSHPfn, geom_type=ogr.wkbMultiLineString )     featureDefn = outLayer.GetLayerDefn()     outFeature = ogr.Feature(featureDefn)     outFeature.SetGeometry(multiline)     outLayer.CreateFeature(outFeature)      ![enter image description here](http://i.stack.imgur.com/VbJEu.png) Any ideas are appreciated. I want to do this in Python ideally with GDAL/OGR (no ArcPy).


title: Drawing a Satellites Coverage FootPrint on an OpenLayers Vector Layer

text: We have been discussing offering at some point in the future satellite footprints on some of the maps we offer. Initially the ISS but some Satellites may be appropriate as well. We have also been looking at this in the context of disaster management, and providing some tools for disaster planners to visualize where different useful satellites may be located to assist. These two A and B excellent questions provide a great deal of information but are way too complex for my old brain to put into practice. Would someone be able to show me the correct Openlayers JScript Code to draw the shape of the footprint on a vector layer, given the Satellites position and height. Accuracy would not be critical as it is a visualization. Maps would be EPSG:4326 and/or EPSG:900913


title: Choosing correct proj4string in R for a set of points

text: I was given a set of points which I eventually will overlay with a polygon shapefile in R. The issue is that the person who provided me with the points did not give them to in long/lat coordinates but the points were already projected (I will not be able to ask them to give me the points in a better format for a while). I was given the following information about the projection of the coordinates:               Top: 583309.373443 m     Bottom: -1996690.626557 m     Left: 221989.400576 m     Right: 3572989.400576 m          Projected Coordinate System: World_Sinusoidal     Projection: Sinusoidal     False_Easting: 0.00000000     False_Northing: 0.00000000     Central_Meridian: -76.00000000     Linear Unit: Meter      Here is the head of the csv I was given:                  [,1]     [,2]     1 2413489 -1995191     2 2416489 -1995191     3 2419489 -1995191     4 2422489 -1995191     5 2425489 -1995191     6 2428489 -1995191      I need to assign the correct projection now to read these points in correctly. Here is the pro4string and SpatialPoints code I have so far:               sinu_prj <- "+proj=sinu +x_0 +y_0 +units=m +no_defs"     brazil_fish_pts <- SpatialPoints(brazil_fish_grid2,  proj4string=CRS(as.character(sinu_prj)), bbox = NULL)      I know that the proj4string code is not correct but I don't know how to fix it? Can anyone be of help?


title: Manipulating distances in python. Problem with ShapeLength

text: I have a question more than a problem. I'm getting the distance of a line using python with arcpy.da.SearchCursor or FeatureClassToNumPyArray. The result I get is a double. Example: arcpy.da.FeatureClassToNumPyArray(Line, "SHAPE@LENGTH", None, SpR) Result: 1.184971580000024 And when I use the Measure tool get the distance in meters: 124,571.381742. The length from Search I use it to add points using 'positionAlongLine()' based on a distance for example: 10,000 meters. How can I convert 10.000 meters to the format required for my points to be inserted precisely on the line?. Example: Meters(Arcmap) => Shape@Length(SearchCursor) 124,571.381742 => 1.184971580000024 10,000 => ? I'm trying to divide my distance of 10,000 m/100, 000 and I getting 0.1, however when inserting and create my points with my python script when I see on arcmap the result feature, I measure the distance between the points and I have a distance of 10,600m.               distance_m = 10000     distance_double = float(distance_m/100000) #How convert meters to double LENGTH          with arcpy.da.SearchCursor(Line, ("SHAPE@", "SHAPE@LENGTH"), "", spRef) as cursor:                 for row in cursor:                     leng = row[0].length                     pts.append(arcpy.PointGeometry(row[0].firstPoint))                     i = 0                     while i < row[0].length:                         pts.append(row[0].positionAlongLine(i))                         i += float(distance_double)          segPts = arcpy.CreateFeatureclass_management(outPath, outName+"_FC", 'POINT','','','',spRef)          edit = arcpy.da.Editor(outPath)     icursor = arcpy.da.InsertCursor(segPts, ('SHAPE@'))     edit.startEditing(False,False)     for newPt in pts:         icursor.insertRow((newPt,))          edit.stopOperation()     edit.stopEditing(True)      Thanks for helping a beginner in python and GIS


title: arcpy.Clip_analysis() Automatically Affixing .shp.shp to Target Feature Class Name?

text: I created a Python script that was previously working fine, but is now giving me an error when I try to use the arcpy.Clip_analysis() function in arcmaps 10.2.1. I am not adding a filetype extnsion to any of my feature class names when creating them. Previously, this worked fine. Now however i looks like ".shp" is being affixed to my feature class name twice, making the name invalid. Here is my code:                   in_features = os.path.join(arcMapsAuto.pc_gdb, arcMapsAuto.merged_feature_class )         clip_features = os.path.join(arcMapsAuto.shapes_gdb, "Terminal_Polygon_Merge")          out_feature_class = os.path.join(arcMapsAuto.tc_gdb, "All_Terminal_Clip")               # Check to see if feature class already exists         if arcpy.Exists(out_feature_class):             print "All_Terminal_Clip creation: FAIL."             print " --All_Terminal_Clip already exists."             #return               print in_features         print clip_features         print out_feature_class                      # printed out:  {geodatabas}\All_Terminal_Clip                   # Execute Clip         try:             arcpy.Clip_analysis(in_features, clip_features, out_feature_class)             print " All_Terminal_Clip creation: SUCCESSS."              except Exception as err:             arcpy.AddError(err)             print "All_Terminal_Clip creation: FAIL."             print err              # printed error:  ERROR 000210: Cannot create output {geodatabase}\All_Terminal_Clip.shp.shp ERROR 000354: The name contains invalid characters                  return      I am missing something? Is arcpy.Clip_analysis() affixing .shp.shp to my target feature class name, and is this a new issue? I swear I've used the same code before without a problem? Any ideas on how I can fix this? > {geodatabase}\All_Terminal_Clip > > All_Terminal_Clip creation: FAIL. > > ERROR 000210: Cannot create output {geodatabase}\All_Terminal_Clip.shp.shp > ERROR 000354: The name contains invalid characters


title: ArcSDE 10.2: How to enable and create Geodatabase in Oracle without ArcGIS Desktop

text: Since there is no Post Installation Wizard anymore in ArcSDE 10.2 - how can i enable und create a geodatabase in oracle without ArcGIS Desktop? We have only an EDN description for ArcGIS Server/SDE, but no desktop license. Is there any documentation about it? Any help is appreciated!!


title: How to create a range domain in ArcGIS with an interval of 10

text: I need to create a range domain in a file geodatabase. The user will be able to select a value ranging from 0 and to 200 with an interval of 10. For example, 0, 10, 20, ..., 200. How can I make and configure a domain type like that?


title: Is Lock Graphics functionality from ArcGIS 9 ArcScript available for ArcGIS 10.1?

text: I have created grids and have some layout graphics and I dont want to accidentally move these elements. I found a old script **here** which does what I want, and I understand that the Production Mapping extension does this but is more advanced. I cannot recompile from ArcGIS 9 to 10 (not smart enough), but do you know if there is an finished DLL file like the one in the link but for 10.1? I also want to know if it is possible to lock Map frame element and/or save layout settings in ArcMap 10.1


title: How to change rasterlayer color?

text: I have rasterdataset in gdb. I create layer from this rasterdataset.                   private static void CreateLayer()         {             ...             IRasterLayer rlMosaicDataset = new RasterLayer();             rlMosaicDataset.CreateFromDataset(rasterDataset);// name in gdb                  IRasterRGBRenderer stretchRenderer = new RasterRGBRenderer();             IRasterRenderer rasterRenderer = (IRasterRenderer)stretchRenderer;             IRaster raster = rasterDataset.CreateDefaultRaster();             rasterRenderer.Raster = raster;             rasterRenderer.ResamplingType = rstResamplingTypes.RSP_BilinearInterpolation;  // smooth             rasterRenderer.Update();             IRasterStretch stretchType = (IRasterStretch)rasterRenderer;             stretchType.StretchType =                 esriRasterStretchTypesEnum.esriRasterStretch_ESRI;             stretchType.StandardDeviationsParam = 2; // what is?             rlMosaicDataset.Renderer = rasterRenderer;                  rlMosaicDataset.MaximumScale = 0.1;             rlMosaicDataset.MinimumScale = 500000000;             ...         }      This layer is black-and-white imagery. I'd like that its color was green and red. (In ArcMap I can do: Layer Properties -> Symbology -> Stretched -> Color Ramp -> here I can choose strip with green, yellow, orange, red colors)


title: How do I copy or edit a cell of attribute table and paste in text file?

text: I'm new in VBA. I want to copy one cell or optional cell in attribute table and paste in a text file. I'm using ArcGis 10.1.


title: Area Boundaries, London

text: Hope someone can help me with this. I'm having a little trouble finding a data set which has polygons of London districts as listed here: http://en.wikipedia.org/wiki/London_districts I've checked Ordnance Survey and none of their boundary data sets corresponds to the above list. I have Borough polygons & ward polygons but need the above as (based on my knowledge of a couple of boroughs) these equate to sub-borough areas people would commonly refer to in everyday life. Can anyone help? Thanks


title: Formatting OSM Data For Use in ESRI Address Locator

text: I've been given the task of researching the possibility of setting up an ArcGIS Geocoding Service that can work with OSM data. I am aware of Nominatim, but since these service(s) will be running on Windows Server 2008, and the customer has already invested big $$ on the ArcGIS Server license - I don't think its going to be an option. From what I've read, installing Nominatim on Windows is a daunting task in itself. If someone knows of any tutorials of how to accomplish this, please let me know. I am far from an expert in either geocoding or OSM, but I would think the biggest barrier to using OSM data in an ESRI address locator is lack of address ranges. Is there any way at all to get address ranges for a given OSM road segment? Any other tips or suggestions welcome.


title: Convert Shapefiles to GeoJSON?

text: I'm trying to convert Shapefiles from the State of Texas Senate Districts website into a GeoJSON file for use with Google Maps API. I'm using QGIS to facilitate but I can't get the projection/coordinates right. So far I've tried Google Mercator (EPSG:900913) and WGS 84 (EPSG:4326) with no luck. I'm using the OpenLayers to verify results and it keeps placing Texas near the Gulf of Guinea (0 lat, 0 lng).


title: What is the most efficient way to publish different thematic versions of the same large feature class?

text: I have a large feature class containing many thousands of polygons. For several reasons, at least part of the time in our development/deployment process, it lives in a file GDB. It is part of a map service that we publish where there are different layers based on multiple (4-6) thematic versions of this same feature class (rendered, filtered, and labeled based on different attribute values). Basically, I've been adding the feature class to the map, renaming the layer, setting the display attributes, then adding it again and doing the same things. This works just fine, but I have a nagging feeling that it's inefficient, particularly for a map service. The consumer of the map service has limited abilities to change rendering.


title: How to export x/y coordinates from geodata?

text: I am a novice about Geographic Information Systems. Now I want to get real line segment dataset to do an experiment. So I searched it, and found the site that includes the real data: http://www.wsdot.wa.gov/mapsdata/geodatacatalog/ From this, I downloaded railroad data and opened it with Quantum GIS. I want to get the end points of each line segment (or feature) in the data. When I click the one feature in vector layer, I can see the firstX, firstY, lastX and lastY in 'Identify Results' window. I think this information is what I needed. However I cannot find the way to export this information to csv file or text file. I searched, but all of these solutions are not working for me: Get list of coordinates for points in a layer? How to get X and Y coordinate values of the map? How can I export end points of each line segment in such dataset?


title: How to decrease GeoTIFF file size by eliminating or negating the blackfill?

text: I have a 5 band, 16-bit orthorectified GeoTIFF from a satellite image that contains a significant amount of blackfill and has a file size of over 21GB. I would like to decrease the file size by eliminating or negating the blackfill. I have attempted to use the SPARSE_OK option in GDAL as such: gdal_translate -of GTiff -co "BIGTIFF=YES" -co "SPARSE_OK=TRUE" -co "TILED=YES" input.tif output.tif This has had no effect on the file size. Is my understanding of how SPARSE_OK works incorrect or is there a potential problem with the SPARSE_OK option? Is there a better way of reducing the file size? I have tried "COMPRESS=DEFLATE", which works, but then the product is veryyyy slow is loading into the viewing software.


title: constant deviation projecting lat lon to itm(israel tranverse mercator

text: i use two different ways to project lat lon data taken from google earth to ITM (Israel Transverse Mercator) and export the data to arcmap   same ways i try to unproject ITM data and export the data to google earth   both ways i get a deviation of about 90 m   first way: cs2c2   \- ` cs2cs +proj=longlat +ellps=WGS84 +datum=WGS84 +to +proj=tmerc +lat_0=31.73439361111111 +lon_0=35.20451694444445 +k=1.000007 +x_0=219529.584 +y_0=626907.39 +ellps=GRS80 +units=m file.txt>output.txt `   second way: DNRgarmin : i downloaded DNRgarmin open a shp file or a kml file and try to convert it   both ways i get a deviation of about 90m does anyone have an idea wht i am doing wrong?


title: Urban Heat Island USING ENVI

text: I have downloaded MODIS images for a specified area.. Its my first time that I am working on UHI (Urban Heat Island).. I need to know the methods to work on it.. I am using ENVI for this.


title: find route distance & time taken to travel between two points using leaflet

text: I am able to find route between two points using cloud mode & leaflet using this. CloudMade Routing with Leaflet: route is not displayed But I also want to find distance & time taken to travel from first point to last point.


title: ArcGIS for Server 10.3 and Portal for ArcGIS

text: I just had a meeting with ESRI distributor and learned that Portal for ArcGIS will not be an extension but be a part of ArcGIS for Server beginning from 10.3. We currently do not have either Server or Portal products and will invest on both of them. In current version (10.2) if we buy ArcGIS Server and the Portal Extension what we will pay for Portal will be as much as what we pay to Server. Below is the question they can not answer; Is there a possibility to have the price of the Server 10.3 product (which will include Portal as well) remain the same with 10.2? Because if that will be the case, investing for both products now will be a total loss of money. But I do not really think they will double the Server price in 10.3 version since it will then be a problem for the customers that require Server but not Portal. Anyone had similar past experiences with ESRI and give me and idea how to proceed?


title: Selecting a Polygon based on the location of a point

text: Currently I have a polygon layer made up of three shapes. All the shapes are listed in the attribute table. One of these three shapes intersects a point. I would like to select just the shape that has an intersection of that point. I am guessing I use the select tool with a piece of code to select the correct shape? Thanks for reading Daniel


title: Arcpy Table To Table Conversion Error

text: I am running a script that uses an excel spreadsheet to define and output maps, kinda like my own data driven pages set up. The only input into the tool is the path of the spreadsheet. When I run the script for the first time in a session I get an error stating I have used an invalid SQL statement, however this is not the problem I cant fix...when I try and run again I get Error 00732 Telling me that the Table (Spreadsheet) does not exist or is not supported. If I run the same excel sheet through the stand alone TableToTable geoprocessing tool it works fine. Almost like the Python session is locking itself out of the excel spreadsheet or something along those lines? Here is the portion of the script where the error occurs:               XLSxPath = arcpy.GetParameterAsText(0)          #Check Temp PDF Folder exists, otherwise create     TempPdfFolder = r"C:\Temp"     if not os.path.exists(TempPdfFolder):         os.makedirs(TempPdfFolder)     arcpy.AddMessage("First level")          #Read in admin project details     AdminTable = XLSxPath + "\\" + "Admin"     AdminDBF = arcpy.TableToTable_conversion(AdminTable,r"C:\Temp", "KBM_Admin" + ".dbf")     AdminSC = arcpy.SearchCursor(AdminDBF)     for vars in AdminSC:         Client = vars.Client         Project = vars.Project_Na         JobDir = vars.GIS_Direct      Has anyone else had a similar problem? I am running Arc 10.2.1


title: qgis how to assign a field to map tips

text: I am a beginning user of QGIS 1.8. I cannot find the way to assign a field to 'map tips'. There seems to be no way in 'properties' of the layer to do so.


title: Road network Analysis

text: I have a Road layer and Point features of some Bus stops from the same area.I want find the Road distance between the Bus stop.Please help me find it. I have QGIS,ArcGIS,and FME. Thanks, RKM


title: Unable to connect PostgreSQL to ArcMap 10.1

text: I'm trying to connect my local PostgreSQL database to ArcMap 10.1 but for some reason I'm unable to connect to the server. Any ideas of what might be causing the issue. For the instance option in ArcMap I'm using the name of my server which is PostgreSQL 9.3 but i get "the unable to connect to data server" error. I have also tried localhost:5432 but i get the same error. Attached there is a pic showing the error the inputs and PostgreSQL server properties. Thanks for the help ![enter image description here](http://i.stack.imgur.com/81wSh.jpg)


title: R - Efficient way to segment rasters into sub-regions and sample... without using extract()

text: My question is: Is there a more efficient way to sample from a sub-region (defined by a shapefile read into a spatial object) of multiple rasters without having to use `extract()` for each region on every raster? Data: I have a large study area represented by approximately 50 raster layers, one per variable (e.g. elevation, slope, vegetation, etc...) my study area is divided into 10 sub-region Methods: I need to draw random samples (using `sampleRandom()` or `sample()` depending on object) from each region of each raster; therefore approximately 500 samples. Issue: I can do this relatively easy by using the `extract()` function for small datasets, but it will take forever with these data sets. I am hoping there is a way to sample from within the area of the sub-regions without having to use `extract()` Below is code (some thankfully borrowed from J. Evans) that shows the basics of what I am doing.   This works fine for small raster data sets, but trying to use `extract()` on large raster data sets fails. Thank you for your time.               require(raster)          r1 <- raster(ncol=10, nrow=10, xmn=0, ymn =0, xmx=100, ymx=100)     r1[] <- runif(ncell(r1),0,100)     r2 <- raster(ncol=10, nrow=10, xmn=0, ymn =0, xmx=100, ymx=100)     r2[] <- runif(ncell(r1),100,200)     r3 <- raster(ncol=10, nrow=10, xmn=0, ymn =0, xmx=100, ymx=100)     r3[] <- runif(ncell(r1),200,300)          region1 <- rbind(c(0,0), c(50,0), c(50,50), c(0,50), c(0,0))     region2 <- rbind(c(50,0), c(100,0), c(100,50), c(50,50), c(50,0))     region3 <- rbind(c(0,50), c(50,50), c(50,100), c(0,100), c(0,50))     region4 <- rbind(c(50,50), c(100,50), c(100,100), c(50,100), c(50,50))          polys <- SpatialPolygons(list(Polygons(list(Polygon(region1)), "region1"),Polygons(list(Polygon(region2)), "region2"), Polygons(list(Polygon(region3)), "region3"), Polygons(list(Polygon(region4)), "region4")))          plot(r1)     plot(polys, add=TRUE)           r_stack <- stack(r1, r2, r3)     names(r_stack) <- c("r1", "r2", "r3")          for(i in 1:length(polys)){         region <- polys[i]         for(k in 1:length(names(r_stack))){             raster_region <- extract(r_stack[[k]], region)             sample_values <- sample(unlist(raster_region), 10)             print(sample_values)             ### add samples to some matrix...         }     }     ### Do analysis on values...


title: Summarize polygon area using another polygon

text: I have a shapefile representing public/private lands in Montana. This is a multipart shapefile where I've simplified attributes so that #1==private and #0==public. Next I have a shapefile of counties for the state. How could I easily summarize the total area of public vs private by county? Do I need to rasterize and go about it that way? I have an ArcEditor and Spatial Analyst license. thanks in advance, J


title: How to check Topology in QGIS?

text: I have two Shapefiles with Polygons ( City and Forerst) Is there a Plugin who is able check the topology (QGIS)? Like "Knots do not overlap" or something like that? I dont want that they overlap. Does anybody know a Plugin? Would be great!


title: How to determine vector between two lat/lon points

text: For all the math geniuses out there...given two lat/lon points (point A and point B), what is the simplest way to find the vector needed to reach point B from point A. For example, the vector needed to reach 36, -85 from 36, -82 would be 270°. I am not entirely sure if "vector" is the technically correct term so if there is a more accurate term, comment and let me know. The fact that the earth is a sphere is throwing me off...I feel as if I could tackle this problem if given a plane but the sphere complicates it too much for me. A chunk of code would be ideal but a simple algorithm that I can translate into code myself is perfectly fine also.


title: Is there an ArcGIS to PostGIS function cross-reference table?

text: I was looking back at the presentation called "PostGIS Tips for Power Users" by Paul Ramsey and got to thinking... Has anyone ever seen an ArcGIS Geoprocessing Tool to PostGIS function cross-reference table? In particular, pages 21-27. For ArcGIS users attempting to transition to PostGIS, I would think this could be very helpful. I know of the following resources, but not aware of anything to help users learn the right terminology in PostGIS from what they already know in ArcGIS.   * PostGIS 2.0 Cheatsheet   * PostGIS Manual For example, I was thinking something like this: ![enter image description here](http://i.stack.imgur.com/w8F04.jpg)


title: How to remove clouds in High resolution Imagery(.5 mGeoeye)

text: i have imagery with high resolution .5m, large parts covered with clouds,so it's obscure many urban areas. thank you.


title: Coordinate Conversion in PostGIS

text: I have some problem with the representation of the coordinates. I need to convert the coordinates from "EPSG: 32633" to "EPSG: 4326", and then create geometries in PostGIS. For the conversion (from java) I use GeoTools.               double xc = 586073.93799999997;     double yc = 4644834.9100000001;     Coordinate c = new Coordinate(xc, yc);     Point p = gf.createPoint(c);     CoordinateReferenceSystem utmCrs = CRS.decode("EPSG:32633");     CoordinateReferenceSystem tCrs = CRS.decode("EPSG:4326");     MathTransform mathTransform = CRS.findMathTransform(utmCrs, tCrs, false);     Point p1 = (Point) JTS.transform(p, mathTransform);     System.out.println(p1.getCoordinate());      with the output: > (40.271915834122524, 17.83405351767443, NaN) However, this representation is wrong for PostGIS.               create table coordinate(     id int  primary key,     tipo varchar(15)     );     SELECT AddGeometryColumn('coordinate','luogo','4326','POINT',2);     INSERT INTO coordinate(id, tipo, luogo) VALUES(1,'calcare', ST_GeomFromText('POINT(40.271915834122524, 17.83405351767443)', 4326));      returns > ERROR: parse error - invalid geometry HINT: "POINT(40.271915834122524, " > <\-- parse error at position 26 within geometr how can I solve my problem? ... the wkt of GeoTools is the same as that contained in spatial_ref_sys:               crs = CRS.decode("EPSG:4326");     wkt = crs.toWKT();     System.out.println("\n\nwkt for EPSG:4326");          wkt for EPSG:4326     GEOGCS["WGS 84",        DATUM["World Geodetic System 1984",          SPHEROID["WGS 84", 6378137.0, 298.257223563, AUTHORITY["EPSG","7030"]],          AUTHORITY["EPSG","6326"]],        PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]],        UNIT["degree", 0.017453292519943295],        AXIS["Geodetic latitude", NORTH],        AXIS["Geodetic longitude", EAST],        AUTHORITY["EPSG","4326"]]


title: General strategies for debugging model builder

text: When coding a standard program in a text editor there are many tools for debugging problems (syntax errors, debuggers, logging information to the console, unit tests, etc.) What are some ways to do these things in model builder? Is it possible? Should I just use the Python API instead? Thanks!


title: Road-construction (alignment) tools

text: Are there any tools available in any GIS environment for route alignment? To be specific, the tool should take a DEM and road-construction parameters -- distance and gradient, roadfill, drainage crossings, starting point and end point -- as input? I have used path distance in ArcGIS but it is very basic. Are there any tools which have been developed completely for this specific task in GIS environment. Or is there any better option for this task?


title: How to Build Query Table from two OLE DB tables using ArcPy?

text: When I run the following code to create a subset of a table, the error says "An invalid SQL statement was used" [Query]. The error shows in line 39. Is there any idea you can figure this out?               import arcpy     # MakeQueryTableOLEDB.py     # Description: Create a query table from two OLE DB tables using a limited set of     # fields and establishing an equal join.          try:         # Local variables...         tableList = ["Database Connections/OLE DB Connection.odc/ORADBA.DLY_AVG_DATA"]              fieldList = [["ORADBA.DLY_AVG_DATA.OBJECTID", 'OBJECTID'],["ORADBA.DLY_AVG_DATA.STATION_ID", 'STATION_ID'],\                     ["ORADBA.DLY_AVG_DATA.SENSOR_NUM", 'SENSOR_NUM'],["ORADBA.DLY_AVG_DATA.DMONTH", 'DMONTH'],\                     [ "ORADBA.DLY_AVG_DATA.DDAY", 'DDAY'], ["ORADBA.DLY_AVG_DATA.NO_OF_VALUES", 'NO_OF_VALUES'],\                     ["ORADBA.DLY_AVG_DATA.DLY_AVG", 'DLY_AVG']]         whereClause = "ORADBA.DLY_AVG_DATA.STATION_ID = 'ADM'"         keyField = "ORADBA.DLY_AVG_DATA.OBJECTID"         lyrName = "Query"         #QueryTable = "C:/RG/Projects/CriticalDataNetwork/QueryTableTest.dbf"         # Make Query Table...         arcpy.MakeQueryTable_management(tableList,  lyrName, "USE_KEY_FIELDS", keyField, fieldList, whereClause)              # Print the total rows         print arcpy.GetCount_management(lyrName)              # Print the fields         Fields = arcpy.ListFields(lyrName)         for field in Fields:             print Field.name                   # Save as a dBASE file         arcpy.CopyRows_management(lyrName, "C:/RG/Projects/CriticalDataNetwork/QueryTest.dbf")          except Exception, e:         # If an error occurred, print line number and error message         import traceback, sys         tb = sys.exc_info()[2]         print "Line %i" % tb.tb_lineno         print e.message


title: Clipping several raster layers to the same exact size QGIS

text: I'm new to QGIS, sorry if this is a newbie question... I need to clip about 40 different raster layers from one vector layer (that has the boundaries of 4 neighboring countries) but I need them all to have the same exact size, which didn't happen when I tried clipping the rasters one by one... How can I do this? I thought about using the first cut raster to clip all others but I can't clip a raster with another raster... or can I? Thanks for any help


title: Editing PostGIS topology in a web browser

text: I am trying to conceptualize the best process for editing topology from PostGIS in a web browser. Ideally what I would like to accomplish is the loading of a topology (imagine the lower 48 states of the US), adjusting the borders, and updating a PostGIS table with the changes. I am fairly comfortable with the delivery, representation, and uploading of the geometry, but I am a bit stuck on the editing part. Although it seems like Open Layers 3 supports topology editing (despite this answer), it does not seem to enforce any rules (for example, it allows the user to drag a boundary to overlap another polygon), nor does it seem like Open Layers supports exporting vector layers as a GeoJSON, though I'm sure it's possible to write that functionality. Likewise, Leaflet offers the excellent Leaflet.draw plugin along with GeoJSON export functionality, but does not offer the ability to edit a topology. Unless I'm missing something, it seems like my best bet at this point is to write a tool that allows for the topological editing of polygons, most likely for Leaflet because I am more comfortable with it. I imagine the process would be something like:   * Click on a boundary   * Find all features that have a geometry at this location   * Create a vertex at this location in all found feature geometries   * While dragging, update all features that share this vertex   * When done dragging, make sure that no segments of the altered geometries intersect other features Ideally it would function very similarly to the Open Layers example above, but would also enforce some rules and make exporting easy. Does this sound like a reasonable process? Am I missing anything? Do you have any specific suggestions for any part of the process, i.e. any existing tools that accomplish these pieces of functionality? Thank you!


title: True long int fields in a File Geodatabase

text: I have been wondering about the 'long ints' in FGDBs. They are 32 bit, which is just a regular int in most places. Is it possible to actually have a field be a 64 bit int? I have some bigint's (64 bit) in a PostGIS database that I am translating into a FGDB with ogr2ogr and they all max out in long ints (2,147,483,647) in the FGDB. I know I can cast them as text and it will get the numbers in there, but that makes it harder to deal with them as numbers. I played with the doubles a little bit, but they start to go into scientific notation once they are too large, which is also not useful.


title: Pythonaddins SaveDialog filter

text: If you look in the **ArcGIS 10.1 help** , the Pythonaddins module SaveDialog shows the following:               SaveDialog({title}, {name_text}, {starting_location})      But if you open the python window in ArcMap, you get a little something extra:               SaveDialog({title}, {name_text}, {starting_location}, {filter}, {filter_label})      When you use the SaveDialog, it doesn't behave very well, and has a dropdown for "Save as type" that I cannot figure out how to access. It would seem that, although the help doesn't contain instructions for doing so, they have built in a means to control the dialog so you can filter it for only certain types. I've tried many alternatives, but cannot figure out how to get this dialog to let me filter. Does anyone have an idea of how this would be implemented (since it's completely undocumented)?


title: Upgraded ESRI .NET ADF application (9.3 to 10.1) intermittently missing service issue

text: I also asked this question on the ESRI forum but I'm not sure if anyone will see it through the tumbleweeds over there. I had an application based on the .NET Web ADF originally built in 9.3. The server is being upgraded to 10.1 so the application needs to move with it (re- architecting is not an option). The automated upgrade utility for moving ADF applications to 10.1 was not the smoothest process and I have since fixed a couple of very obvious bugs in the ESRI JavaScript. However there is one problem I can't figure out: Intermittently, after a full page reload, the map fails to render. The application calls ESRI.ArcGIS.ADF.Web.UI.WebControls.MapHandler.ashx requesting a URL for the new map image. Sometimes that handler returns               {"code":400,"type":"error","errorMessage":"Resource <service name> not found"}      The application then requests an image from ESRI.ArcGIS.ADF.Web.MimeImage.ashx passing an ImgID parameter that seems to be randomly made-up - I can't see it in any previous requests and it obviously isn't provided by MapHandler.ashx. The server then responds with a 404. If I reload the application again the same requests are sent with exactly the same parameters (in the example I saw the post-reload request even passed the same 't' parameter which I assume should represent time) and a successful response is returned. Can anyone tell me why my service is intermittently reported as unavailable? ArcGIS Server is not logging any errors or warnings. Based on this and how quickly the request is re-issued and successfully answered I don't think the service is actually having trouble. Comparing a working and a failing request in Fiddler gives me no clues as the same headers and cookies seem to be sent. **EDIT** The only reason I can imagine is if the server control 'MapResourceManager' or something associated with it isn't fully initialised by the time the first request is sent, but as such timing issues are the responsibility of the framework I would be surprised if this was the problem. However I've seen widespread use of setTimeout in the ADF's JavaScript and I'm really hoping this issue isn't the result of some flawed assumption about timing.


title: How to change tooltip of an esriControl?

text: Is there a way to change an esriControl's tooltip? I tried using the ICommand interface but the tooltip propertry is read-only. If I have to build a whole new ICommand and/or ITool is there anyway I can extend an esriControl so I do not have to build a new control from scratch?


title: km2 or sqkm - which is more correct?

text: What is more correct, km2 or sqkm? I'd prefer to use km with a superscript 2, but it's not always possible to format text that way so I was just wondering if there is a convention that should be used when you can't use superscript formatting.


title: How to place elevation numbers on contours with uphill orientation and position?

text: I want to orient my elevation labels uphill. Is that possible? I'm currently using the labeling tool from QGIS (not that one in layer settings). I want the labels curved and uphill oriented (and I would like to define how regularly the numbers are displayed). It's all for a print-map.


title: Looking for a QGIS Architecture Diagram

text: I would like to know the architecture of QGIS: relationships between QGIS and other open sources, and what elements are implemented in the QGIS. Usually, There is a diagram showing the architecture of a software such as arcobjects. Even it is simple like Geotools's architecture, it is useful to understand internal structure of a software. Thanks a lot.


title: Splitting DXF layers into Shapfile layers in QGIS?

text: When i load a dxf into QGIS all the layers are present in the attribute table. Is it possible to split this "entities" layer in QGIS into its component parts and save each component part as an individual shapefile. I am at present selecting each layer and right clicking and saving as .shp. But is there a fast way because this way is time consuming and boring! Cheers.


title: Creating Digital Elevation Model (DEM) from Google Earth Image?

text: I am very new to GIS and given the task to create DEM's of agricultural area. I have studied the advantages of DEM in agriculture now i just want to know how i should start to develop the DEM's ? How i can create the Digital Elevation Model (DEM) from Google Earth Images ??


title: How do I add "padding" to features when applying d.out.file to map display on GRASS?

text: Thing is this. I've done up a script to automate some spatial interpolation and eventually output the interpolated raster and contour lines as an image of the map display. I'm doing this using `d.out.file`. The issue now is, the features (mainly an island + the contours) are rendered right to the edge of the output image. I've tried using `d.zoom` to zoom out prior to the `d.out.file` command but to no avail. Does any GRASS expert here know how I can place a "padding" so that my features don't end up right on the edge of the image? Post-processing won't cut it either as there are grids too. Also, it has to be done programmatically and not via the GUI. **edit:** one more thing, I've got a raster mask in place so I'm not sure if that is causing the issue. Thanks in advance!


title: How to place points along a line in a specific offset using python / arcpy?

text: I have a street shapefile with information about the street name, house number from / to on the left side and house number from / to on the right side. My aim is to create a point shapefile that visualizes all houses along the streets. A first prototype is already implemented. With a Search-cursor I am reading out the total number of houses by building the range by using the paramaters **house number from / to on the left side** and **house number from / to on the right side**. The total amount is divided by the length of the street so that I have a value that show me in which distance I have to place the points along the line to map all houses. The result is that all house numbers are placed like a string of pearls. Is there a way to place the numbers with an specific offset along the line? The script is written in python using arcpy and ArcGIS 10.1 SP1. Thanks!


title: How can I calculate the optimal routes while touching all available roads in a network?

text: For preperation of the winter season we want to calculate the most optimal routes for sprinkling salt on the roads. The analysis knows the following criteria: \- vehicles start and stop at a single loading point \- all available roads need to be sprinkled with salt \- one route can take no longer than a certin time (lets assume 2 hours) \- because of the limited load of salt per vehicle, the distance of a route is limited to the amount of salt available. (lets assume 10 km) Network analyst of ArcGIS (10.0) assumes you have a start- and an endpoint to calculate the route. However, in this case it is not about calculating the fastest route from origin to destination, but about the most optimal routes to cover as much road distance as possible within a limited time frame. Together with a collegue we're thinking about solutions, but we would appreciate it if you could give feedback on this interesting challenge. Now we are thinking about calculating midpoints for every road section and use those as destinations to calculate the route. If someone has better suggestions, please share! Cheers, Mark


title: What GIS analysis techniques were used in determining Gale Crater as the landing site for the Curiosity Rover?

text: What specific GIS analysis was used in determining Gale Crater as the landing site for the Curiosity Rover? I'm assuming factors such as slope, elevation, ruggedness, etc were used as well as potential for finding water (what determines this). I'm looking for a journal article (if it exists) or some technical document outlining the processes used for choosing Gale crater over other potential landing sites. Any help would be appreciated. This is just for my own curiosity (pun intended).


title: Is there a quicker way of checking if a feature exists with OGR

text: I'm writing lot's of features to a FileGDB using OGR. The slow point in my script is checking that the feature isn't already in the Feature Class, based on the ID. Currently I'm using a function and:               crfLayer.SetAttributeFilter("ID = " + str(crfID))     Count = crfLayer.GetFeatureCount()     if Count == 0:         return 0     else:         return 1      I've tried looping creating a list from the file and checking if the ID is in it, however this is slower. Does anybody have any alternatives to the code above?


title: Is there a way to Collect Lat-Long in a Google Docs Form?

text: I'm helping out a small non profit. They are collecting some simple data using a Google Docs Form. The need to collect exact latlong (upto 5 meters accuracy is fine). Ideally the end user would click on a map and that lat-long would be submitted along with the rest of the form. I have seen the mapbox blog talk about this here. But that is using the mapbox API. I was wondering if there exists a ready made template to collect the point on top of Google Maps. I don't want to spend half a day converting the mapbox sample to Google Maps API, and then find out that something already existed. So Is anyone aware of template where a point was collected from Google Maps in Google Forms?


title: Is it possible to create OGC wms service based on the arcgis mxd file and data in java

text: We are trying to create a OGC wms using J2EE, and we know there are some opensource solution, like geoserver-wms. However most of them will use the `SLD` to render the data to an image. Which does not meet our requirement. Then I wonder if we can read the `mxd` file in java and generate data according to the wms parameters? * * * Update: We have arcgis engine license. Since we are trying to create a OGC WMS servcie, so we will have to generate an image and return it to user.


title: ArcObjects Geoprocessing Failures

text: I am working on launching the GenerateNearTable using the ArcObjects Assembly method. I have the advanced license (10.1), so the tool is available. The code to launch geoprocessing is:                  Private Sub generateNearTable(ByRef FeatureLayer)              'Grab the feature class for the GP tool         Dim lcount As Long = 10         Dim closest As Boolean = False              'Setup the geoprocessor         Dim gp As ESRI.ArcGIS.Geoprocessor.Geoprocessor = New ESRI.ArcGIS.Geoprocessor.Geoprocessor         Dim neartable As GenerateNearTable = New GenerateNearTable()                   neartable.in_features = FeatureLayer         neartable.near_features = FeatureLayer         neartable.out_table = "C:\Users\jlaura\Desktop\test4" ' Hard code to test.         neartable.closest = closest         neartable.closest_count = lcount 'This should be a parameter              Dim addtomap As Boolean = gp.AddOutputsToMap         gp.SetEnvironmentValue("workspace", "C:\Users\jlaura\Desktop\test")         gp.TemporaryMapLayers = False         gp.AddOutputsToMap = False              Try             gp.Execute(neartable, Nothing)             returnmessages(gp)         Catch ex As Exception             returnmessages(gp)         End Try               End Sub      This is causing an error 99999. The results are being returned to the results tab, so I am able to double click on the failed tool and run it successfully via ArcMap. New to ArcObjects, so I'm definitely missing something obvious.   1. Do I need to check out a license like I would with ArcPy?   2. Am I parameterizing in a strange way?   3. Do I need to pass the dataset (IDataset) instead of the featureclass? #This failed when I tried it, but perhaps I missed something... As an extension (that I can open as a new question if people think it is too disconnected) - to programatically set the output table, can I generate an empty, 'fieldless' standalone table and let the tool handle writing. Updated Code:               Private Sub generateNearTable(ByVal distlayer, ByVal knn, ByVal DistanceTableOut)              Dim FeatureLayer = GetFLayerByName(distlayer)         Dim NearFeatureLayer = GetFLayerByName(distlayer) 'Fails if I use this.              'Setup Params         Dim lcount As Long = knn         Dim closest As Boolean = False              'Setup the geoprocessor         Dim gp As New ESRI.ArcGIS.Geoprocessor.Geoprocessor         Dim genneartable As GenerateNearTable = New GenerateNearTable()              'Setup the Gp Environment         Dim addtomap As Boolean = gp.AddOutputsToMap         gp.TemporaryMapLayers = False         gp.AddOutputsToMap = True 'Add the derived table to the map doc.         gp.OverwriteOutput = True              'Set Params         genneartable.in_features = FeatureLayer         genneartable.near_features = "C:\Users\jlaura\Desktop\Zunil\Zunil_secondaries_proj.shp"         genneartable.out_table = DistanceTableOut         genneartable.closest = closest         genneartable.closest_count = lcount 'This should be a parameter              'Execute         Try             gp.Execute(genneartable, Nothing)             If gp.MessageCount > 0 Then                 For count As Integer = 0 To gp.MessageCount - 1                     MsgBox(gp.GetMessage(count), MsgBoxStyle.OkCancel, "Working")                 Next             End If         Catch ex As Exception              End Try               End Sub      EDIT: I am doing something wrong with the feature layer inputs. By hard coding the path to the shapefiles I am able to successfully run the tool. Back to the drawing board, re: featurelayer access. EDIT2: The tools runs properly if the featurelayer is passed for the `in_features` and hard coded for the near_features. In this case should I be reaccessing the feature layer from the input?