title: Why does data imported from AutoCAD appear different in ArcMap 10.1?

text: I am importing information from AutoCAD into Arcmap 10.1. I found after some issues that I need to give the incoming information geographic location, which works fine. The problem is that all the data appears the same. I have two maps with different points of collections, but when I bring over both maps, they give me the same data. Pic1: is the data that I am pulling over from AutoCad   ![enter image description here](http://i.stack.imgur.com/Ol26j.png) Pic2: is what appears in Arcmap 10.1.   ![enter image description here](http://i.stack.imgur.com/vTB8G.png) If you look at the bottom, you will notice that they are both suppose to be the same, but they are not. The title is different, and also when I bring it over to ArcMap 10.1, it has alot more points of data in the different areas.   I have bringing over alot of different maps, is it possible that my maps are reading data from a different map that I brought over before I was looking at these?


title: Exporting basemap & Layers as an image?

text: I need a open-source server-side component for a web-application that can export basemaps + layers to an image format (saved on a db). Being slightly new to GIS, I've managed to find a few tools that do this on the client-side.   * leaflet's plugin Leaflet.print does this. I'm using leaflet to display CloudMade tiles on the browser. it has npm packages, so I'm guessing it can run server-side as well.   * OpenLayers 3 (beta) has this as an example but I've not tested it out.   * MapFish Print    * GeoServer which uses MapFish Print again. Is there any component that could take in map bounds & a features as geoJson and output image file on the server-side ? It would not be possible to use ArcGIS or proprietary tools


title: Sqlite, Python 2.7 and Spatialite

text: I am trying to spatially enable my sqlite database:               import sqlite3     db = r'c:\temp\mydata.sqlite'     conn = sqlite3.connect(db)     conn.enable_load_extension(True)     cur = conn.cursor()     cur.execute(r"SELECT load_extension("<path>\libspatialite-4.dll");")     conn.commit()     conn.close()     del conn      Python just stops working when I run this code. It doesn't crash in windows 7 x64, it just closes. It's very strange. No error is raised. When I run the commands directly from python.exe in the command propmpt, I see python crashing: ![spatialite crashing python](http://i.stack.imgur.com/XCybZ.png) Am I doing anything wrong? Thank you


title: Specifics of ESRI's use of FlexNet Publisher TrustedStorage at ArcGIS 10

text: So at ArcGIS 10, ESRI adopted Flexera's **FlexNet Publisher Trusted Storage** module for secure exchange of license file "fulfillment records" and transfer of product activations. Question what network traffic occurs and what corresponding firewall ports should be open to allow the traffic to pass? Also, what conflicts will arrise with other applications running the same FlexNet Publisher licensing tools-- Autodesk, SolidWorks, Adobe, etc. and how best to make things coexist? The FLexNet Publisher License Administration Guide doesn't address those specifics. And when looking at the ArcGIS 10 installations you'll now see a **FLEXNetLicensingService.exe** running and calls to **libFNP.dll** on license server and client. That in addition to the "old" FlexNet Publisher **lmgrd.exe** and ESRI's FlexNet **ARCGIS.EXE** Vendor daemon on the license server. How does it all behave? And how do we fix it when collisions or version conflicts occur? A very scary thing is that a number of Apps using FNP are placing configurations and executable files into the common directory:               C:\Program Files\Common Files\Macrovision Shared\FLEXnet Publisher      Where vendors can do what they want and update their FlexNet components at their convenience, without regard for release level of other applications present on the system. Looking at the FLEXERA site I don't have a sense that they make much of an effort to support the end users. License exchanges over VPN has been especially troublesome for end users.


title: Programmatically alter AliasName for fields defined in layer file (not underlying source)

text: Is there a way to alter the AliasName on a layer file only, and not the underlying shapefile or feature class? I have a few hundred layers that define joins between hundreds of tables with census data sitting in an MS Access database, and a shapefile with census tracts. I like to keep it this way to not have to store the feature geography for each of those data files. Each layer file maps out hundreds of variables with intuitive names such as B19325_089. All said, there are several thousand. I would like to change the NameAlias property to have meaningful aliases/variable labels based on a lookup value, but, can it be done/defined in the layer file and not in the actual geography file, as all the layer files create the join to the same shapefile? I got off to a decent start on feature classes in a file GDB using Mark Cederholm's and Matt Wilkie's suggestions, but am, as noted, curious about whether it is possible to have the alias definition reside in the layer files. Initially, I thought it very simple and doable with the geoprocessor object as ESRI lists the AliasName property as read/write, but that appears to not be possible, and arcobjects to be necessary.


title: Problems installing TopoJSON

text: I experienced problems with installing TopoJSON.               C:\Program Files (x86)\MSBuild\Microsoft.Cpp\v4.0\V120\Microsoft.Cpp.Platform.targets(64,5): error MSB8020: The build tools for Visual Studio 2010 (Platform oolset = 'v100') cannot be found. To build using the v100 build tools, please install Visual Studio 2010 build tools.  Alternatively, you may upgrade to the current Visual Studio tools by selecting the Project menu or right-click the solution, and then selecting "Upgrade Solution...". [C:\Users\Sigita\AppData\Roaming\npm\node_modules\topojson\node_modules\d3\node_modules\jsdom\node_modules\contextify\build\contextify.vcxproj]     gyp ERR! build error     gyp ERR! stack Error: `C:\Program Files (x86)\MSBuild\12.0\bin\msbuild.exe` failed with exit code: 1     gyp ERR! stack     at ChildProcess.onExit (C:\Program Files\nodejs\node_modules\npm\node_modules\node-gyp\lib\build.js:267:23)     gyp ERR! stack     at ChildProcess.EventEmitter.emit (events.js:98:17)     gyp ERR! stack     at Process.ChildProcess._handle.onexit (child_process.js:797:12)     gyp ERR! System Windows_NT 6.2.9200     gyp ERR! command "node" "C:\\Program Files\\nodejs\\node_modules\\npm\\node_modules\\node-gyp\\bin\\node-gyp.js" "rebuild"     gyp ERR! cwd C:\Users\Sigita\AppData\Roaming\npm\node_modules\topojson\node_modules\d3\node_modules\jsdom\node_modules\contextify     gyp ERR! node -v v0.10.26     gyp ERR! node-gyp -v v0.12.2     gyp ERR! not ok     ...     npm ERR! contextify@0.1.8 install: `node-gyp rebuild`     npm ERR! Exit status 1     npm ERR!     npm ERR! Failed at the contextify@0.1.8 install script.     npm ERR! This is most likely a problem with the contextify package,     npm ERR! not with npm itself.     npm ERR! Tell the author that this fails on your system:     npm ERR!     node-gyp rebuild     npm ERR! You can get their info via:     npm ERR!     npm owner ls contextify     npm ERR! There is likely additional logging output above.          npm ERR! System Windows_NT 6.2.9200     npm ERR! command "C:\\Program Files\\nodejs\\\\node.exe" "C:\\Program Files\\nodejs\\node_modules\\npm\\bin\\npm-cli.js" "install" "-g" "topojson"     npm ERR! cwd C:\Users\Sigita     npm ERR! node -v v0.10.26     npm ERR! npm -v 1.4.3     npm ERR! code ELIFECYCLE     npm ERR!     npm ERR! Additional logging details can be found in:     npm ERR!     C:\Users\Sigita\npm-debug.log     npm ERR! not ok code 0      I'm running it on Windows 8.1. with all updates installed, npm version 1.4.3., gyp 0.13.0 and node.js 0.10.26.


title: How to show QGIS plugin dialog always on top?

text: Is it possible to define for a pyqgis plugin, that the plugin dialog is always shown on top, after starting it from the toolbar?


title: How can I define the layer CRS for a memory layer in pyqgis

text: how can I define the CRS for a memory layer from the qgis project settings                   vpoi = QgsVectorLayer("Point", "vectpoi", "memory")         feature = QgsFeature()         feature.setGeometry(QgsGeometry.fromPoint(QgsPoint(10,10)))         provider = vpoly.dataProvider()         ...      and how can I define the CRS from another loaded layer? thanks a lot


title: Auto resize map to 100% after content pane is removed in ArcGIS JSAPI

text: Been at this for hours and I can't seem to figure it out. I'm using AG-JSAPI 3.6 I have a map layout with a bordercontainer and two contentpanes, One is the map div and the other is the left div (which holds the legend, and basemaps and such as child elements). Anyway, I have a checkbox element on the footer that when clicked, fires off a function that (should ideally) make the left div disappear and the map div expand to the width of the browser.               function toggleDiv() {         // If statement to check if the checkbox is checked.     if (dijit.byId("ToggleLeft").checked) {         //Sets the left div's display CSS property to inline and then setting the         //width property        dojo.setStyle(dijit.byId("leftDiv").domNode, 'display', 'inline');        dojo.byId("leftDiv").style.width = '285px';          }     else     {         // Sets the display css property to "none" which makes the element disappear         dojo.setStyle(dijit.byId("leftDiv").domNode, 'display', 'none');          };      } The thing is, when I resize the browser, the map automatically resizes to the width of the browser. However, I want to do this as soon as the leftDiv element disappears. What js code would I have to use to get this to happen? Any help would be appreciated. Thanks


title: How to iterate through maps & add new layer to data frame and legend of each using ArcPy?

text: Newbie for Python and wondering if the following is doable. Problem is like this: I have two sets of data:   1. one mxd file for each of multiple cities (all stored in the same folder) containing links to multiple layers of data.   2. additional dataset (multiple fields; point shapefile) with data for all cities stored in a file geodatabase (*.gdb) My goal is to have a python script to:   1. open all mxd file iteratively   2. open the additional dataset and select the data based on a query (City name) and add to the mxd file. 2) Label the added data in the map by selecting desired field 2) add the legend of the added data at the bottom of the original legend. Help with scripting would be highly appreciated.


title: gp.CreateFeatureClass(..) and gp.CopyFeatures(..) isn't copying all fields of feature class

text: How do I get the ArcGIS 9.3 geoprocessor to copy all of the fields of a feature class? Here is a description of what I'm trying to do: I am modifying an existing script created by a developer not on the program anymore. The script's purpose is to extract a given feature set into a zipped .gdb. Therefore, the script takes a feature set as input. In the FeatureSet that I am attempting to export (which is a feature set representing major lakes), the fields are DEPTH, NAME, OBJECTID, Shape_Area, Shape_Length, SQKM, SQMI, and SURF_ELEV. The .gdb resulting from this script, however, only contains OBJECTID, Shape ( _where'd that come from??_ ), Shape_Length, and Shape_Area. I am exporting the feature set via               gp.CopyFeatures(featureSet, scratchGdb + os.path.sep + lyrName)      The .gdb is successfully created, but only _some_ fields are there! Why am I losing some of the fields in this feature class when exporting it to a .gdb??? Through using Komodo IDE, I was able to remote debug the script, and I saw that the input feature set only had the four fields that ended up in the exported .gdb. Why the heck is this? Does this have to do with setting a schema for the feature set? Please help! Thanks!


title: Automate converting labels to annotation in ArcMap at multiple scales?

text: I need to be able to convert labels to Annotation feature classes at multiple scale levels for both address and street data. Similar in concept to generating tiled layers for a web service, I have 5 to 10 pre-defined extents that I need to export annotation for. The resulting annotation will be used in an ArcEngine application that is restricted to those zoom levels. When done manually, I turn on labels for the layers I want labeled, right-click each layer and choose 'Convert Labels to Annotation, save them to an annotation feature class, then repeat for each scale level. Does anyone know of a way of accomplishing this programmatically? All of my code so far has been in Python, but if necessary I'm open to other languages.


title: Valid ranges for bounding box in WMS 1.1.1

text: I'm connecting to a WMS server and found that the longitudes LatLonBoundingBox specified for some of the values are slightly outside of the +/-180 range (such as 180.05). Are these values considered valid (I couldn't find anything in the schema that seemed to enforcer a range)? If they are valid, should they be treated as rounding errors (maybe the values were auto generated) and clamped to +/-180, or should they be normalized in the +/-180 range, so 180.05 might become -179.95? Also note that in this particular case, the parent layer's bounds are listed as -90/90 and -180/180. Thanks.


title: Why do features in QGIS drawn on same datum as Google imagery from OpenLayers move?

text: I have just started using QGIS. I've loaded the "openlayers" plugin and added a google hybrid layer, which is based off the WGS84 datum to my understanding. Now I create a shapefile and add a new line vector layer based on the same datum. I draw a line and then zoom out. Now, the line has moved with regards to the google map (i'm drawing a road and now the road no longer lines up with the original road I'm using for a guide). My end goal is to make a georeferenced shapefile with some roads. Is my methodology flawed or is their a problem with my plugin (assuming the former of course!)?


title: Digitizing with openlayer plugin

text: I am facing problem when working with Google earth imagery. when I am digitizing at 1:1000 it just shift the shapefiles when i even pan the image. checked for CRS and everything but it doesn't work for me.


title: Digitized vectors don't match while zooming

text: I'm trying to digitize in QGIS with Google Maps as the background. Digitizing with constant zoom works fine. But if the zoom is changed, the digitized vectors don't match the Google Maps background image. Thanks for your help! Regards Hans


title: How to get Google satellite images (via the plugin) to behave well in QGIS with a map created with a different CRS?

text: How to get Google satellite images (via the plugin) to behave well in QGIS with a map created with a difference CRS? I have a problem using the Google Image layer plug in - at some scales, it won’t match the rest of the project (but does at others). The project started with Tiff files imported from GB Ordinance Survey and set the CRS to EPSG:27700 - OSGB 1936 / British National Grid (which I think these tiles default to). I then added a couple of shapefile layers (set to the same CRS). One has point data, the other is a polygon which defines a site boundary. I then added (via the plug in) the google satellite image layer. Whenever I add this layer to any project, this always defaults to a view where the whole world is repeated several times across the screen. If I zoom in by one click from there, it seems to find itself ok and the site boundary matches the OS tile and the satellite image. However on this project (and a previous one), when I zoom in too far, suddenly the site boundary doesn't match the same area on the google image (about 1:1000) - it is quite a jump. I can zoom back out any amount and the polygon and google image match fine. When the polygon and google image don’t match up, if I bring the OS tile above the Google image, I can see that the polygon and OS tile have remained lined- up. I know Google use a different CRS - WGS 84 / Pseudo Mercator, which I assume is the root of the problem, but there must be a fix. I’ve tried changing the google image layer’s CRS to match the rest of the project with the obvious result that it goes even further off (“sorry we don’t have image for here”) and right clicking on the OS layer and selecting “set project CRS from layer” doesn't help. Any guidance appreciated - Thanks!


title: QGIS Google satellite layer moving each time it's turned on

text: I have imported location data in a text layer, CRS = WGS84 and have put in Google Satellite and Google Street (WGS84/pseudo mercator) in but each time I switch on a Google layer it moves relative to the location points - 121 metres last time I measured. The movement appears to be random in direction and I believe the imagery also distorts. The project CRS is WGS84 EPSG:4326 (in the bottom right corner). Can anyone help this first time user please?


title: Created polygon appears to move when zooming in (qgis)?

text: I am using QGIS 2.20 on a Windows 8.1 64-bit laptop. I have a shapefile that includes a polygon that I created from digitizing a Google Satellite image from OpenLayers Plugin. Here is my problem: when I zoom in or out, the polygon MOVES (or at least appears to). Here is an example: http://imgur.com/0MOReoo then, ALL I did was zoom in using my mouse wheel, and the polygon layer was NOT in editing mode, and I got this: http://imgur.com/YKxXQRJ Any ideas? Is the polygon moving, or is the Google layer moving?


title: Create shapes using google streets layer of openlayers plugin

text: I am using openlayers plugin to use as the base layer of streets google. I have attached the picture that you can see that the more I give to zoom overlay layer moves. When I used another shape as I skirted base with it snapping, but I know that with google streets layer can not do it. My question is: how can I work around this situation automatically? Today I'm doing the editing always a given zoom level to know that this level until the overlap is correct.


title: Attribute Table does not honor subtypes - Flex

text: The domain has 8 coded values. I added a subtype so it would filter by subtype. The attribute inspector honors the subtypes. As you can see on the pictures, the subtype "Forestry" has two coded domain values. However the attribute table, as you can see on the other picture does not honor the subtypes and shows all the coded values for the domain. Does anybody know why this happens and if there is a fix? Thank you so much for any help! ![Correct filtering of domains by subtype](http://i.stack.imgur.com/9x0W4.jpg) ![Subtype not being honored in table](http://i.stack.imgur.com/hNEvw.jpg)


title: Cannot describe coverage. WCS

text: I have Map Server 6.0.3 installed and running in my local machine. The same for Apache HTTP Server 2.2.22 I have created a .map file using NOTEPAD++ for a WCS service that I want to access using either QGIS and ArcMap as CLIENTS. 1.- When I use the URL:               http://localhost:80/cgi-bin/mapserv.exe?map=/ms4w/apps/raster.map      QGIS CAN list the service (see it), fine, but when I try to add the layer into the canvas it reports:               WCS provider: Cannot get test dataset.     Raster layer: Provider is not valid (provider: wcs, URI: cache=PreferNetwork&     crs=EPSG:4326&format=GTiff&identifier=gebco&url=http://localhost:80/cgi-bin     /mapserv.exe?map%3D/ms4w/tmp/raster.map      If I try adding a GDAL DRIVER (xml file), see below, as a RASTER file:               <WCS_GDAL>       <ServiceURL>http://localhost:80/cgi-bin/mapserv.exe?map=/ms4w/app/raster.map</ServiceURL>       <version>1.0.0</version>       <CoverageName>gebco</CoverageName>       <supportedCRSs>EPSG:4326</supportedCRSs>       <PreferredFormat>gtiff</PreferredFormat>     </WCS_GDAL>      I get the following error:               msLoadMap(): Regular expression error. MS_DEFAULT_MAPFILE_PATTERN validation failed.     msEvalRegex(): Regular expression error. String failed expression test.       2.- ArcMap only list (but again not displaying it) the image by adding just:               http://localhost:80/cgi-bin/mapserv.exe?map=/ms4w/apps/raster.map      So, the same problem as using QGIS.. 3.- The only thing which seems to be working so far is adding the following URL into the browser:               http://localhost:80/cgi-bin/mapserv.exe?map=/ms4w/apps/raster.map     &SERVICE=wcs     &VERSION=1.0.0     &REQUEST=GetCoverage     &coverage=gebco     &CRS=EPSG:4326     &BBOX=32.4434100,-28.0,55.0017433,0.0     &WIDTH=2707&HEIGHT=3360     &FORMAT=GTiff      After that a pop-up window ask you to download a tiFF into your computer. But again using either ArcMap or QGIS as clients I am not able to display the raster using MapServer. My MAP FILE:               MAP         NAME WCS_server         STATUS ON         DEBUG ON       SIZE 2707 3360       PROJECTION         "init=epsg:4326"       END       EXTENT 32.4434100 -27.9999999974 55.0017433333 2.63999311301e-009       UNITS DD       SHAPEPATH "/ms4w/apps"       CONFIG "MS_ERRORFILE" "/ms4w/tmp/ms_error.txt"       MAXSIZE 3361            OUTPUTFORMAT         NAME GEOTIFF_16         DRIVER "GDAL/GTiff"         MIMETYPE "image/tiff"         IMAGEMODE INT16         EXTENSION "tif"       END            WEB         IMAGEPATH "C:/ms4w/tmp/ms_tmp/"          IMAGEURL "/ms_tmp/"         METADATA           "wcs_label" "GMap WCS Demo Server"                               "wcs_onlineresource" "http://localhost:80/cgi-bin/mapserv.exe?"            "wcs_enable_request" "*"                                   END       END            LAYER         NAME gebco         METADATA           "wcs_label" "gebco_blue"           "wcs_rangeset_name" "Range 1"           "wcs_rangeset_label" "My Label"         END         TYPE RASTER         STATUS ON         DATA gebco.tif         PROJECTION           "init=epsg:4326"         END       END     END      Checking the ms_error.txt for errors, it reports the following when accessing the service using QGIS:               [Wed Feb 26 11:41:55 2014].129000 CGI Request 1 on process 4936     [Wed Feb 26 11:41:55 2014].129000 msWCSParseRequest20(): WCS server error. Parameter     'Version' is already set.     [Wed Feb 26 11:41:55 2014].129000 msWCSDispatch20(): Parse error occurred.     [Wed Feb 26 11:41:55 2014].209000 CGI Request 1 on process 1900     [Wed Feb 26 11:41:55 2014].209000 msWCSDispatch20(): version and service are not     compliant with WCS 2.0.0     [Wed Feb 26 11:41:55 2014].209000 msWCSParseRequest(): request is KVP.      And when trying to add the layer to the map (QGIS):               [Wed Feb 26 11:54:05 2014].220000 CGI Request 1 on process 4472     [Wed Feb 26 11:54:05 2014].220000 msWCSParseRequest20(): WCS server error. Parameter     'Version' is already set.     [Wed Feb 26 11:54:05 2014].220000 msWCSDispatch20(): Parse error occurred.     [Wed Feb 26 11:54:08 2014].904000 CGI Request 1 on process 5892     [Wed Feb 26 11:54:08 2014].904000 msWCSParseRequest20(): WCS server error. Parameter     'Version' is already set.     [Wed Feb 26 11:54:08 2014].904000 msWCSDispatch20(): Parse error occurred.     [Wed Feb 26 11:54:08 2014].994000 CGI Request 1 on process 6072     [Wed Feb 26 11:54:08 2014].994000 msWCSDispatch20(): version and service are not     compliant with WCS 2.0.0     [Wed Feb 26 11:54:08 2014].994000 msWCSParseRequest(): request is KVP.     [Wed Feb 26 11:54:09 2014].84000 CGI Request 1 on process 5944     [Wed Feb 26 11:54:09 2014].84000 msWCSParseRequest20(): WCS server error. Parameter     'Version' is already set.     [Wed Feb 26 11:54:09 2014].84000 msWCSDispatch20(): Parse error occurred.      Now, when connecting the service via ArcGIS:               [Wed Feb 26 11:57:32 2014].802000 CGI Request 1 on process 5832     [Wed Feb 26 11:57:32 2014].802000 msWCSDispatch20(): version and service are not     compliant with WCS 2.0.0     [Wed Feb 26 11:57:32 2014].802000 msWCSParseRequest(): request is KVP.     [Wed Feb 26 11:57:32 2014].822000 Format 'png24' ignored since mimetype 'image/png'     duplicates another outputFormatObj.     [Wed Feb 26 11:57:32 2014].822000 Format 'GTiff' ignored since mimetype 'image/tiff'     duplicates another outputFormatObj.     [Wed Feb 26 11:57:32 2014].932000 CGI Request 1 on process 3540     [Wed Feb 26 11:57:32 2014].932000 msWCSDispatch20(): version and service are not     compliant with WCS 2.0.0     [Wed Feb 26 11:57:32 2014].932000 msWCSParseRequest(): request is KVP.     [Wed Feb 26 11:57:32 2014].932000 msWCSParseRequest(): Whole String: gebco     [Wed Feb 26 11:57:32 2014].932000 msLoadProjectionString(): Projection library error.     Unrecognised OGC CRS def 'urn:ogc:def:crs:OGC::imageCRS'.      And finally, when trying to display the image in ArcGIS:               [Wed Feb 26 12:00:57 2014].16000 CGI Request 1 on process 4632     [Wed Feb 26 12:00:57 2014].16000 msWCSDispatch20(): version and service are not     compliant with WCS 2.0.0     [Wed Feb 26 12:00:57 2014].16000 msWCSParseRequest(): request is KVP.     [Wed Feb 26 12:00:57 2014].16000 msWCSParseRequest(): Whole String: gebco     [Wed Feb 26 12:00:57 2014].16000 msLoadProjectionString(): Projection library error.     Unrecognised OGC CRS def 'urn:ogc:def:crs:OGC::imageCRS'.     [Wed Feb 26 12:00:58 2014].166000 CGI Request 1 on process 4440     [Wed Feb 26 12:00:58 2014].166000 msWCSDispatch20(): version and service are not     compliant with WCS 2.0.0     [Wed Feb 26 12:00:58 2014].166000 msWCSParseRequest(): request is KVP.     [Wed Feb 26 12:00:58 2014].166000 msWCSParseRequest(): Whole String: gebco     [Wed Feb 26 12:00:58 2014].166000 msWCSGetCoverageBands11(): WCS server error.     RangeSubset field name malformed, expected 'Range 1', got RangeSubset=Range      So, there are many things to have a look at, and obviously the errors in the ms_error.txt file gives us a clue, how to solve them?. But for instance, why it says:               version and service are not compliant with WCS 2.0.0?      When I am using VERSION=1.0.0 in the requests?


title: Adding Multipoint Features to ArcGIS Online

text: I was attempting to upload a zipped multipoint shapefile to AGO when I received this error:"Unable to import this shapefile. (Unable to get property 'intersects' of undefined or null reference)". I was unsure why I got the error but I noticed that AGO states at this time it does not support multipoint shapefiles. So is there some other format besides a .shp that you can upload a multipoint feature like a gdb feature class or an SDE layer or is it a no go all together.


title: Is there a tool to track edits that were made?

text: Is there a tool/report that keeps track of what edits were made that I could run prior to or after reconciling my versions? This would help me to know if there are any changes I need to make to any maps or databases that are disconnected from our SDE. I'm using ArcGIS 10.1


title: Extract values using buffers - Equation

text: I have run a linear regression and obtained an equation that describes the link between a group size and the area used by this group. Now, I would like to enter this equation into the GIS to obtain buffers of different sizes depending on the group size. For instance, a group of 50 individuals may end up with a buffer of 100m diameter. I was going to add a field in the attribute table "buffer" and enter the equation. After that, I was wondering if I could use my different buffers to extract values from rasters in a simple way? I mean that I would like to extract the average value but based on my buffer size and without having to make a hundred shapes... I hope that I have been clear Thank you so much for your time,


title: Topology rules for a 3d network

text: I am currently writing a spec to which people may end up using to guide the digitization of a large 3d network (exterior pedestrian network + building interiors). In 2d network drawing, it seems to be best practice to stick to a set of topology rules. The most prominent of these is enforcing no line intersection. If two links cross with grade separation (e.g. bridge/tunnel) then this is specified with data either on links or nodes (the link approach seems easier to me as is used in OS ITN for example). The advantage of this is that you know all intersections are errors, and it defends against a common problem - hard to detect visually - of intersecting lines when you intended to split them. In a fully 3d network, the likely problems are different:   1. accidental line intersection may still happen when the user intends to break lines, as 3d networks can have a lot of coplanar sections (e.g. floors of buildings).    2. depending on the view position when editing, many lines may appear to be joined when they are not, and vice versa.   3. depending on the view position when editing, lines may appear to connect to the wrong line, for example in plan view a line directly above/below the intended line _**Does anyone have experience of or recommendations for handling these problems?_** Here are my thoughts so far.   * Creating nodes (as in a point data set) to mark all joins between links, is a possibility and would help with problem 2. However, it would probably be quite time consuming, and fail to fix problem 3.   * A different approach may be to specify that all link endpoints must join to other links unless explicitly marked as dead ends (perhaps with point data). This fixes problem 2 with less work, but not 3.   * No intersections could be a requirement within a 3d tolerance (though ArcGIS seems to falsely report intersections where they don't exist in 3d) - this would fix problem 1   * Within buildings we could enforce that links remain level (do not change in z) unless explicitly marked as doing so - this would mostly fix problem 3 It is perhaps a good idea to split the dataset into parts: one for the exterior network and one for each building interior. In that case you could enforce different rules outdoors and indoors:   * Outdoors any apparent intersection in plan view is an error, and links must be noded to show intent. (The links will still be 3d to represent the exact shape of foot bridges etc). Dead ends are an error unless explicitly marked.   * Indoors, apparent intersections from plan view are allowed, but any links that change in height are an error unless explicitly marked as doing so. 3d intersections are an error within a certain 3d tolerance. Dead ends are an error unless explicitly marked.   * It may sometimes be easier to draw buildings floor by floor, in which case, apparent intersection in plan view is an error, height changes are not allowed and dead ends are an error unless explicitly marked.   * Apparent dead ends will appear where outdoor, indoor and floor networks join, so joining points between the different datasets will need to be explicitly marked in a similar way to dead ends. Any comments on these policies or alternative ideas would be welcome!


title: Arcgis rest server date bug - One day earlier - Flex

text: I am having that arcgis rest server bug where the attribute table and the attribute inspector show the dates as one day earlier. I found a possible solution on the web but there is an error in the following line because the attribute inspector has no "form" property. Is there any way to get around this? Thanks error in this line var formItems:* = editor.attributeInspector.form.getChildren();               / Place the following code below the function populateEditor      private var dateFormatter:DateFormatter = new DateFormatter(); private const newDateFormat:String = “DD/MM/YYYY”; private function attributeInspector_showFeatureHandler(event:AttributeInspectorEvent):void { formatDateFieldDisplay(event); } private function attributeInspector_updateFeatureHandler(event:AttributeInspectorEvent):void { formatDateFieldDisplay(event); } private function formatDateFieldDisplay(event:AttributeInspectorEvent):void { if(!(event.featureLayer && event.feature && event.feature.attributes)) return; // This ensures that we only validate when we are working at a field level (not an AttributeInspector level) if(event.field) { if(event.field.type == Field.TYPE_DATE) { if(event.newValue) // If it comes exclusively from an update … format the display { // Change the date format as desired event.target.formatString = newDateFormat; } } } else if(event.type.toLowerCase() == “showfeature”) { var formItems:* = editor.attributeInspector.form.getChildren();                 // Loop items and trigger update event       for (var i:Number = 0; i < formItems.length; i++)       {          if(formItems[i].data && formItems[i].data is Field)          {             var field:Field = formItems[i].data;             if(field.type == Field.TYPE_DATE &&                field.name in event.feature.attributes &&                    event.feature.attributes[field.name])             {                var event:AttributeInspectorEvent = new AttributeInspectorEvent(“updateFeature”,                                                                  true,                                                                  event.featureLayer,                                                                  event.feature,                                                                  field,                                                                  event.feature.attributes[field.name],                                                                  event.feature.attributes[field.name]);                var formItem:* = formItems[i].getChildren();                if(formItem && formItem.length > 0)                   formItem[0].dispatchEvent(event);             }          }       }      } } // This function should be in a separate util file to be reused if necessary by other classes. // This function comes with the ArcGIS Viewer for Flex 2.4. // If Esri updates this function so should the developer. private function msToDate(ms:Number, dateFormat:String, useUTC:Boolean):String { var date:Date = new Date(ms); if (date.milliseconds == 999) // Workaround for REST bug date.milliseconds++; if (useUTC) date.minutes += date.timezoneOffset; if (dateFormat) { dateFormatter.formatString = dateFormat; var result:String = dateFormatter.format(date);                 if (result)          return result;       else          return dateFormatter.error;      } else { return date.toLocaleString(); } } Source: http://gisengineering.blogspot.co.nz/2011_09_01_archive.html


title: Clip data frame using "Outline of Selected graphic(s)"

text: I did a subtraction of a data frame in ArcGIS using the following steps. I prepared a layout in ArcGIS. There is a lot of empty space at the top of layout page. So I want to remove this space. I used a graphic operation to remove it.   1. Created a graphic rectangle that the same size as the Data Frame   2. Created another rectangle covering empty space   3. Selected rectangle 2, then rectangle 1, and chose--> Graphic Opperations --> Subtract.    4. In the last step, I wanted to clip the subtracted image by "Outline of Selected graphic". But 'Outline of selected graphic' is not activated. I can't click it. I'm wondered why this is not activated? Thanks a lot.


title: Best practice for serving data to OpenLayers client

text: We have some polyline and point data in SHP format and would like to serve that data to an OpenLayers client. First we thought about converting the data to GeoJSON and send it to the client as a whole package, but then we saw that the data is quite large (>500MB) and downloading such an amount as a whole is not a good idea. So we came up with the idea to dynamically serve the data in a tiled format, using WMS. We're not really experienced when it comes to do stuff like this, but at the same time, serving some SHPs as WMS shouldn't be too difficult, right? So what is the easiest and smoothest way to do that? We've heard of GeoServer and MapServer, which one is easier to deploy (We're on CentOS 6). And is the caching of tiles included with serving WMS or does it have to be configured separately?


title: Problems with gxp Print plugin via opengeosuite sdk

text: I am trying to import the print plugin within an application built using the Opengeo Suite sdk and I have run into issues as is described in the following thread https://getsatisfaction.com/opengeo/topics/doubts_about_the_use_of_print_from_openlayers_ready_gxp After having followed all the instructions provided, I still cannot get the print plugin the work. The print button appears but I get the message               Uncaught TypeError: undefined is not a function           This seems to be on the following line:      var scaleLine = new OpenLayers.Control.ScaleLine({ ...       I have also tried using an alternate approach wherein I have a handler for the print button as is shown below:                   var app = new gxp.Viewer({         portalConfig: {             layout: "border",             region: "center",         ...         ...         ...         // layer sources         sources: {             local: {                 ptype: "gxp_wmscsource",                 url: "/geoserver/wms",                 version: "1.1.1"             },             osm: {                 ptype: "gxp_osmsource"             },             geonetwork: {                 ptype: "gxp_cataloguesource",                 url: "/geonetwork/srv/en/csw",                 title: "Geonetwork"             },             ol: {                 ptype: "gxp_olsource"             }         },              // map and layers         map: {             id: "mymap", // id needed to reference map in portalConfig above             title: "Map",             //projection: "EPSG:900913",             //center: [-10764594.758211, 4523072.3184791],             projection: "EPSG:4326",             center: [-97, 38],             zoom: 3,             layers: [{                 source: "ol",                 type: "OpenLayers.Layer.WMS",                 args: ["World map", "http://vmap0.tiles.osgeo.org/wms/vmap0", {layers: 'basic'}],                 group: "background"             }, {                 source: "ol",                 type: "OpenLayers.Layer",                 args: ["Blank"],                 visibility: false,                 group: "background"                      //source: "osm",                 //name: "mapnik",                 //group: "background"             }//, {                 //source: "local",                 //name: "usa:states",                 //selected: true             //}//, {                 //source: "google",                 //name: "ROADMAP",                 //group: "background"             //}             ],             bbar: [{                 text: "Print...",                 handler: showPrintWindow             }],             items: [{                 xtype: "gx_zoomslider",                 vertical: true,                 height: 100             }]         }          });          function showPrintWindow() {         var printWindow = new Ext.Window({             title: "Print",             modal: true,             border: false,             resizable: false,             width: 360,             autoHeight: true,             items: new GeoExt.ux.PrintPreview({                 autoHeight: true,                 printMapPanel: {                     // limit scales to those that can be previewed                     limitScales: true,                     // no zooming on the map                     map: {controls: [                         new OpenLayers.Control.Navigation({                             zoomBoxEnabled: false,                             zoomWheelEnabled: false                         }),                         new OpenLayers.Control.PanPanel()                     ]}                 },                 printProvider: {                     // using get for remote service access without same origin                     // restriction. For async requests, we would set method to "POST".                     method: "GET",                     //method: "POST",                          // capabilities from script tag in Printing.html.                     capabilities: printCapabilities,                     listeners: {                         "print": function() {printWindow.close();}                     }                 },                 includeLegend: true,                 mapTitle: "PrintMapPanel Demo",                 sourceMap: app.map             })         }).show().center();     }      From within the showPrintWindow function the sourceMap is set to app.map. However, when I click on the print button now, I get the error message:               Uncaught TypeError: Object #<Object> has no method 'getProjection' PrintMapPanel.js:176     GeoExt.PrintMapPanel.Ext.extend.initComponent PrintMapPanel.js:176     Ext.Component ext-all.js:7     L ext-base.js:7     L ext-base.js:7     L ext-base.js:7     L ext-base.js:7     L ext-base.js:7     GeoExt.ux.PrintPreview.Ext.extend.initComponent PrintPreview.js:156     Ext.Component ext-all.js:7     L ext-base.js:7     L ext-base.js:7     L ext-base.js:7     showPrintWindow app.js:185     Ext.Button.Ext.extend.onClick ext-all.js:7     I      I am not entirely sure how I can get this working and any help is greatly appreciated.


title: How to remove empty classes from legend in Arcgis

text: I have a classification image with only 12 classes. I'm preparing layout using this image. when I add legend into layout page, a lot of empty classes loaded together with these 12 classes. Totally 255 classes added to layout page. Then I tried to delete them in Symbology tab of data Properties window. How to remove these empty classes from legend. Is there any way to do this.   Thank you all


title: Maximum zoom levels for tilestache rendering mbtiles

text: I have been given an mbtiles file which contains tiles for zoom levels between 16 and 20. When I preview the mbtiles on tilestache, I am not able to zoom in beyond the 18th level. If I request the map at zoom level 19, tiles of zoom level 18 are being rendered. I have viewed the mbtiles files and sure enough there is text "minzoom16" and "maxzoom20". There is another mbtiles file with tiles from zoom levels 11 to 16. It is rendering tiles requested for zoom level 18 too; beyond zoom level that no tiles. This makes me wonder if there is any internal setting in the tilestache code that limits the server to rendering tiles up to a certain zoom level. I have searched the source, but couldn't find anything in particular; not that I am a python programmer, but still. So, does tilestache limit the zoom levels to a certain number? If yes, how can it be changed? Thank you. EDIT:               "bholakpur":     {        "provider": {"name": "mbtiles", "tileset": "/home/***/****/tiles/bholakpur.mbtiles"},        "preview": {"zoom": 20, "ext": "png"}     }      EDIT 2: Bing layer is the base layer for the map and as far as my tests go Bing has only 18 zoom levels. May be that is the reason the other layers are not requesting tiles beyond the 18th zoom level.


title: Refreshing the GeoCache and Boundaries in OpenGeo suite

text: Currently if I add a new site i can query the layer and see where the point should be however it does not always display in Geoexplorer or OpenLayers preview. However sometimes it does show in Openlayers but not Geoexplorer. I think this is down to the GeoCache, and the Bounding Box around the data. I clear my internet history each time too. Edit: I have tried deleting the Geoserver cache, reseeding the layer at zoom levels 0-15. The odd thing is at zoom level 2, all points appear and then when I zoom in a progressively loose points on the map. However the data is still visible, if I query the layer? Any suggestions Thanks


title: IMap.FeatureSelection get the selected feature field value always return -1

text: I am following the sample code of ArcObjects (vb.net) of IMap::FeatureSelection, I want to get 1 field value of all selected features. unfortunately I always get the value -1, any idea on it? The following is my code:               Public Sub LoopThruSelection2()       Dim pMxDoc As IMxDocument       Dim pMap As IMap       Dim pEnumFeature As IEnumFeature       Dim pEnumFeatureSetup As IEnumFeatureSetup       Dim pFeature As IFeature            pMxDoc = My.ArcMap.Document       pMap = pMxDoc.FocusMap       pEnumFeature = pMap.FeatureSelection       pEnumFeatureSetup = pEnumFeature       pEnumFeatureSetup.AllFields = True       pFeature = pEnumFeature.Next       Do While (Not pFeature Is Nothing)         MsgBox(pFeature.Value(pFeature.Fields.FindField("testdb.DBO.testLayer1.RoadName")))       pFeature = pEnumFeature.Next       Loop     End Sub


title: Is there a way to render lines with a smooth (bezier) curves in QGIS?

text: I am moving from creating maps with Illustrator to making them with QGIS and actual data. One of the things I haven't been able to re-create in QGIS are the smooth curved lines I can get in Illustrator. I am not talking about Join and Cap Styles, but how an entire line is rendered. I am looking for an answer that doesn't include exporting as SVG to Illustrator and finishing the map there. Also, I realize they could be considered an inaccurate representation but, for the most part, these maps are for giving riders an idea of where they are and don't necessarily have to be an exact representation. Here is an example of what I mean: ![enter image description here](http://i.stack.imgur.com/y2RW7.png)


title: Matplotlib downloads to override folder

text: Just completed install of QGIS (2.0.1 Dufour) on a Macbook Pro running on Mavericks via Kyngchaos download page. However, when downloading Matplotlib for the install, while everything ended up the correct python site package (2.7), the matplotlib wound up in a folder within the 2.7 site package called matplotlib-override. Does anyone have an explanation? When in the 'override folder QGIS could not find Matplotlib, but once removed and bumped up a level to the site package (2.7) folder QGIS installed cleanly (still running it through its paces for functionality).


title: How to process MODIS ndvi data using QGIS/GRASS/R?

text: I have MODIS ndvi data (MODIS_Grid_16 days NDVI) from the iberian peninsula (3 tiles), from February of 2000 to August of 2012 (really big data). I want to know how to deal with this kind of files.   My final goal is to extract the ndvi value to points with a associated date, i.e., I have sequencial dated points and I want to know the value of ndvi to that date point.) How do I (automatically) merge the 3 tiles that compose the iberian peninsula for each 16-days (2000-2012)?   (I think if I use R to identify each filename and use gdal merge, it might work, but maybe there's a simplier way) How do I create a time series from it and extract the values associated?


title: Why can't I see all records in my personal GDB in QGIS?

text: I am running QGIS 1.8. I have a layer from an ESRI Personal GDB. There are 431 features in the layer. When you load the layer into QGIS, and look at the Metadata properties, it shows the number of features as 431. ![layer metadata](http://i.stack.imgur.com/6nqnP.png) When I open the attribute table, the table only shows 180 features. ![Layer attribute table](http://i.stack.imgur.com/ueB1j.png) This is a layer that had been loaded in straight from the Personal GDB. No filters or queries had been added to it, so all features should be showing. I first came across this issue because another coworker was trying to access the same layer using QGIS 1.7. He was having the same issue. A couple items of note:   1. This layer was existing in the PGDB, and then using ArcGIS, I copied and replaced a new version from a File GDB into the PGDB.   2. I am able to load other layers from the same PGDB into QGIS where all the features are displayed.   3. The PGDB is a modified MS Access database. If you open the database with MS Access, then open the table for this layer, it contains all 431 records with attributes. Here is a screenshot of the layer in QGIS: ![QGIS Screenshot](http://i.stack.imgur.com/m7uDt.png) Here is a screenshot of the same layer in ArcGIS: ![ArcMap Screenshot](http://i.stack.imgur.com/8ZqHg.png) If anyone has some insight into why QGIS is able to see the actual number of features in the layer, but is only displaying a certain number, I would appreciate it.


title: Web mapping software for a huge raster time series?

text: I'm an ArcGIS Desktop user who is stepping into the world of web mapping for the first time. Following advice I found here on GIS.SE, I started with the free tutorials from OpenGeo. However, I began to realize that most of the demos and examples were geared towards vector data. My main project involves the display of a 300 frame time series of 5000 x 5000 pixel rasters. On my desktop computer these are stored in a single 5000x5000x300 16-bit integer BIP binary, which totals around 8GB. **My goal is to be able to click on a cell (of a single raster in the time series) and have a graph pop up showing the values of that pixel in the 300-element time series.** The raster data used for each time series graph must be stored losslessly, though the overlay maps that are clicked on can be lossy caches. Is there anything that may be better suited for this project (for a novice web GIS developer) than OpenGeo? Or should I just keep going with this setup? For reference, I have programming experience in Python, Java, and PHP. I don't have much experience with SQL. This is an open-ended project so I have plenty of time to learn new languages if need be. I already have a web server, though I'll probably have to switch hosts because GoDaddy doesn't seem to support PostgreSQL without getting a VPS. Thanks for your help! **EDIT: (Jan 13)** I'm still seeking information on exactly how best to store a 3 dimensional 16-bit integer BIP raster and be able to efficiently query a single z-axis "column" of data. I don't want to convert it to a 32-bit format (because that would double its file size from its current 16-bit form).


title: What's the best way to merge various road data into a single layer?

text: My company has two different road layers:               1) County roads     2) Old CAD roads      Each with their own pros and cons. I'm attempting to merge these into a single road layer that will be maintained. The County has pretty good road data as far as location. It is also M-enabled for routing, which is nice. But it's missing road names. Our CAD roads are dated and haven't been updated in a long time, but all the roads do have names. My plan is to use the County roads as the base layer and try to join the CAD road names to it. But I first need to get them to line up correctly. I believe the Integrate tool will work, but I haven't used it yet. But then I need join the CAD road names to the County roads and keep the output M-enabled. Any tips?


title: Proxy.jsp Internal Server Error

text: I am using proxy.jsp file to get WMS GetFeatureInfo popup, I have dropped proxy.jsp file in /usr/lib/cgi-bin and added allowedHosts. I changed OpenLayers.ProxyHost = "/cgi-bin/proxy.jsp?url=";   if I clicked on any feature it is displaying like **Internal Server Error The server encountered an internal error or misconfiguration and was unable to complete your request. Please contact the server administrator, webmaster@localhost and inform them of the time the error occurred, and anything you might have done that may have caused the error. More information about this error may be available in the server error log. Apache/2.2.22 (Ubuntu) Server at 192.168.1.1 Port 80** on the popup. Error I am getting in error log. > [Wed Jul 30 10:57:00 2014] [error] [client 127.0.0.1] (8)Exec format error: > exec of '/usr/lib/cgi-bin/proxy.jsp' failed, referer: > http://xxx.xxx.x.x/rjy/GISShops.html [Wed Jul 30 10:57:00 2014] [error] > [client 127.0.0.1] Premature end of script headers: proxy.jsp, referer: > http://xxx.xxx.x.x/rjy/GISShops.html Can anyone please help me how to cure this error? Here is my proxy.jsp file                   <%@page session="false"%>     <%@page import="java.net.*,java.io.*" %>     <%@page trimDirectiveWhitespaces="true"%>     <%          /**     * This is a white list proxy that could be used the prevent an error due to     * JavaScript Same Origin Policy.     *     * CAUTION: It might break some sites and it's a security risk because     * people can use this proxy to browse the web and possibly do bad and/or illegal stuff     * with it.     * It can load any content type.     * This proxy implementation was inspired by the proxy.cgi script of OpenLayers     * {@link http://openlayers.org}     * To use this in OpenLayers you have to set OpenLayers.ProxyHost = "Url/To/This/Proxyfile/proxy.jsp?";     * within your JavaScript code <br>     * The base code of the proxy has been provided by SNIPPLR     * {@link http://snipplr.com/view/17987/jsp-proxy-for-javascript-applications/}     *     * @author terrestris GmbH & Co. KG     * @author Christian Mayer     * @author Marc Jansen     *     * @license BSD see license.txt     *     */     String[] allowedHosts = {     "www.openlayers.org", "openlayers.org",     "labs.metacarta.com", "world.freemap.in",     "prototype.openmnnd.org", "geo.openplans.org",     "sigma.openplans.org", "demo.opengeo.org",     "www.openstreetmap.org", "sample.azavea.com",      "vmap0.tiles.osgeo.org", "localhost", "localhost:8181",      "localhost:80","192.168.1.1",     "192.168.1.1:8181","xxx.xxx.xxx.xxx",      "xxx.xxx.xxx.xxx:8181"     };     HttpURLConnection con = null;     try {     String reqUrl = request.getQueryString();     String decodedUrl = "";     if (reqUrl != null) {     reqUrl = URLDecoder.decode(reqUrl, "UTF-8");     }     else {     response.setStatus(400);     out.println("ERROR 400: No target specified for proxy.");     }          // extract the host     String host = "";     host = reqUrl.split("\\/")[2];     boolean allowed = false;          // check if host (with port) is in white list     for (String surl : allowedHosts) {     if (host.equalsIgnoreCase(surl)) {     allowed = true;     break;     }     }          // do the proxy action (load requested ressource and transport it to client)     // if host is in white list     if(allowed) {     // replace the white spaces with plus in URL     reqUrl = reqUrl.replaceAll(" ", "+");          // call the requested ressource     URL url = new URL(reqUrl);     con = (HttpURLConnection)url.openConnection();     con.setDoOutput(true);     con.setRequestMethod(request.getMethod());     String reqContenType = request.getContentType();     if(reqContenType != null) {     con.setRequestProperty("Content-Type", reqContenType);     }     else {     con.setRequestProperty("Content-Type", "application/x-www-form-urlencoded");     }          int clength = request.getContentLength();     if(clength > 0) {     con.setDoInput(true);     byte[] idata = new byte[clength];     request.getInputStream().read(idata, 0, clength);     con.getOutputStream().write(idata, 0, clength);     }          // respond to client     response.setContentType(con.getContentType());          BufferedReader rd = new BufferedReader(new InputStreamReader(con.getInputStream()));     String line;     int i = 0;     while ((line = rd.readLine()) != null) {     out.println(line);     }     rd.close();     }     else {     // deny access via HTTP status code 502     response.setStatus(502);     out.println("ERROR 502: This proxy does not allow you to access that location.");     }          } catch(Exception e) {          // resond an internal server error with the stacktrace     // on exception     response.setStatus(500);     byte[] idata = new byte[5000];          if(con.getErrorStream() != null) {     con.getErrorStream().read(idata, 0, 5000);     }          out.println("ERROR 500: An internal server error occured. " + e.getMessage() + " " + new String(idata));     }     %>


title: QGIS - Loading data by SQL-window in the DB-Manager is not working?

text: I have problems to add a subset from a PostGIS-table to QGIS by using 'load as new layer' in the SQL-window of the DB-Manager (QGIS-Version 2.2.0-Valmiera). I defined both the 'column with integer values' and the 'geometry column' (see image 'SQL-window') but the subset is not added to the QGIS-window. What could be the problem? Thank's a lot for any advice. Phil ![enter image description here](http://i.stack.imgur.com/IPuj0.jpg)


title: Permanently Change ArcMap's "Maximum Sample Size" Setting?

text: I've lately been working on feature classes with a huge number of polygons (e.g. all of the Census blocks within the Chicago metro area). I've been needing to symbolize them by quantities a fair amount to display demographic information, but because there are so many polygons I receive the `Maximum sample size reached. Not all records are being used. Use this sample or change maximum sample size.` message every single time. I know that I can can manually change the sample size by pressing the `Sampling...` button in the Classification window, but I'm wondering if it's possible to change the default sample size somewhere, so that I will no longer have to go through this whole rigmarole. I haven't been able to find anything in Esri's documentation.


title: ArcGIS Server security and Windows domains - issues for non-IE browsers

text: I'm looking for some advice / experience in ArcGIS Server security with a Windows domain. Normally in this situation I would do some prototyping but I don't have access to a domain. I'm planning to create ArcGIS map services that are restricted by roles stored in a domain. These map services will be used by a map page in a custom .NET application. Both the application and ArcGIS Server will run on the same domain and potentially in the same IIS site (e.g. Default Web Site). The .NET application will immediately require Windows authentication, so I know that IE users will pass through seamlessly but other browsers will prompt for a login. It's not clear to me if non-IE users will also be prompted for a login when first accessing a restricted map service.   * **Q1** : In the above configuration (same domain, same IIS site) would non-IE users get prompted a second time to access a restricted map service (e.g. `http://domain_name/ArcGIS/.../service_name`), or would the credentials they provided for `http://domain_name/application/` be used for ArcGIS?   * **Q2** : Would this change if a different IIS site (but still the same domain) was used for ArcGIS Server and the .NET application? Another (perhaps simpler) option I'm considering is having my map page's server-side code generate a security token on load by making a request to the ArcGIS token service with expiry time calculated as # hours until 17:30. This token would be passed to the client in JavaScript. In this approach only a single account needs access to the map service. The security token contains > the user's name, expiration time and other information (source)   * **Q3** does this 'other information' include anything that restricts the token to the IP or user-agent of the requestor - anything that will stop my server creating it on behalf of my client? **EDIT** from further research it seems that a token can be explicitly bound to a given IP (provided when creating it) or to a referrer. In this case Q3 is answered 'no' Any thoughts on these three questions would be super welcome. Thanks


title: Can I sell maps produced from ASTER Global DEM V2 data?

text: I don't quite understand the term-of-use of ASTER Global DEM V2 data: > By selecting 'Submit', I agree to the following: > > -I agree to redistribute the ASTER GDEM _only_ to individuals within my > organization or project of intended use or in response to disasters in > support of the GEO Disaster Theme. > > -When presenting or publishing ASTER GDEM data, I agree to include 'ASTER > GDEM is a product of METI and NASA.' **Does it mean I _can't_ sell maps produced using this data?**


title: Looking for free sdk/lib/framework for UI for openlayers app

text: I'm looking for free UI lib/sdk/framework to make openlayers app that I can have the map in one panel , the layer control in other panel and maybe hovering tool bar (instead of putting buttons on the map). Bottom line , I'm looking to easily and fast make good looking and practical UI wise app. Thx!


title: How to create a geodatabase feature class using a geodatabase table schema with ArcPy?

text: Is there a way in ArcPy to create an empty polygon geodatabase feature class that will acquire the schema of a geodatabase table? I know that it's possible in ArcCatalog (New -> Feature Class and then import all the fields from a table) but I wonder if it can be written in a Python script. The CreateFeatureclass_management tool allows only Feature Layers as schema templates.


title: How to best remove rgb values outside of aerial image?

text: I have just converted RGB-TIF images to RGB-ECW for reasons of space and performance in ArcGIS. However, after the conversion, there are relics outside of the images which are not nodata values but have the values 0,1,0 or 1,1,0 or 1,2,1 etc. All very low values which are shown as black remnants after I set nodata values to transparent. (See attatchment) I have ArcInfo and FME at hand and would like to either change these rgb values to 0,0,0 or clip all values outside of the aerial extents. Could anyone offer a tip on the most efficient way to achieve this? Image remnants: ![Image remnants](http://i.stack.imgur.com/o1zUs.jpg)


title: How to apply categorized symbols to polygon layer with PyQgis?

text: I'm fetching a polygon layer from postgres like this:               uri = QgsDataSourceURI()     uri.setConnection("localhost", "5432", "mydb", "postgres", "xxx")     uri.setDataSource("public", "my_table", "geom")     lyr = QgsVectorLayer(uri.uri(), "Blocks", "postgres")     QgsMapLayerRegistry.instance().addMapLayer(lyr)      This layer has an attribute field named _blocks_ having 6 values as 'A-6', 'B-4', 'G-5', 'D-2', 'P-7', 'C-3'. I want to apply _green_ color to that polygon where blocks value is 'B-4' and want to apply red color to the rest of polygons in this layer. How to do this? I've seen code to apply graduated symbols using PyQgis but I want categorized symbols. Thanks!


title: ArcSDE Problems: Microsoft SQL Server 2008 R2 on Windows 2008 R2 Standard (64-Bit)

text: I´m trying to set up ArcSDE 10.2 on a Microsoft SQL Server 2008 R2 running on a Windows Server 2008 R2 Standard 64-Bit machine. I can **create** the database connection, but when I want to **enable** it for spatial data I keep on getting the following error message:               Error: The version of your database server is not supported     ERROR installing/upgrading ArcSDE, Error = -494      This is strange because the ESRI page on System Requirements states, that this combination should work. Any idea what I could do, to get it fixed? Thanks EDIT: I´m using Microsoft SQL Server 2008 R2 64-Bit, or at least that is what Windows tells me (see image). I´ve tried without SP, than with SP1 and then SP2. Neither option worked. ![Installed Eindows Programs](http://i.stack.imgur.com/hwpid.png)


title: how to add a custom SRID to Spatialite

text: I'd like to create a spatialite database and load it with shapefiles using the Mississippi Transverse Mercator (MSTM) projection, but it doesn't have an EPSG definition in the spatialite library of projections. When loaded into a QGIS session, it correctly projects the shapefiles, automatically building a custom proj4 code for it from the esri prj files that come along with the data. I checked that QGIS accurately reprojects this data, on the fly, to online WGS84 Orthophoto (Bing Maps). It does. I'm using qspatialite to load these shapefiles into a single newly created db.sqlite; I try to populate the mandatory SRID field with the custom projection code definition number that QGIS generates, for example SRID= "100003" etc, which is on the dropdown picklist to choose from in the qspatialite plugin . Qspatialite can't load the shapefiles with any custom codes stored in QGIS though; it issues an error and fails to load the data--I figure because that isn't in the Spatialite EPSG library of projections. To double check that the plugin installation works, I did load the shapefiles successfully using a standard EPSG number, but my shapefiles are really in a custom projection CRS. Is there a way I can copy/paste QGIS's custom proj4 definition to the spatialite EPGS library and create a brand new custom SRID definition, then load my data using that new SRID number in the mandatory data loader field? Or, would spatialite-gui be better than qspatialite? Thanks to anyone who can explain how. these are the MSTM projection parameters (not proj4) Scale Factor: 0.9998335 Central Meridian: -89 45' 00'' Central Latitude: 32 30' 00'' False Northing: 1,300,000 meters False Easting: 500,000 meters Projection Units: meters Spheroid: GRS80 Datum: NAD83 --EDIT OF POST 1/24/2013 BELOW-- I just wanted to add a workaround method for bulk-batch importing shapefiles into a spatialite database, using Qspatialite's bulk load OGR-file button. This would be useful if the IMPORT OGR FILE command won't accept a custom EPSG number that you need to use. For myself, I just entered a similar but commonly used EPSG number using that also uses the same DATUM. For example, my custom EPSG: 3814 apparently gets little use and hasn't yet been included in QGIS spatial reference database (Lisboa 1.8.0). I found that the "custom" coordinate system uses NAD83. Instead of entering the SRID of 3814 that causes an import failure, I just used a similar NAD83, in this case EPSG:26995, that's preloaded into QGIS. They import just fine like that. Of course the imported shapefiles will not be projected correctly in QGIS relative to a basemap. To fix that I created a custom coordinate system in QGIS using the +proj code that QGIS generates to correctly project the individual shapefiles that are now importing to a spatialite database, that also included a 3814.prj spatial reference file. THEN after creating that CUSTOM CRS, I was able to load the files from my spatialite database and individually set the projection under the of each layer to my custom spatial reference code that I saved in QGIS. Presto, everything aligns perfectly with my basemap while using spatialite layers. And, I was able to use Qspatialite to bulk-batch over a hundred shapefiles and then a couple of hundred dbf tables of additional information. One by one would have been time consuming.


title: How to programmatically fill a datetime field with datetime.datetime.now() function

text: I have been trying to programmatically fill a datetime field using values coming from **datetime.datetime.now()** function.               #I create the temp shp using a URI     URI = 'MultiPolygon?crs=epsg:3857&field=id:integer&field=x:integer&field=y:integer&field=z:integer&field=bx:integer&field=by:integer&field=bz:integer&field=rating:integer&field=timestamp:datetime&index=yes'     # temp shp     mem_layer = QgsVectorLayer(URI, 'temp', 'memory')          # when I try to fill out the timestamp field with a variable now     now = datetime.datetime.now()     # I get an empty or zero value field      It looks like I have to cast the variable first to a string to make it work but in that case I lose the capability to make calculation over the field (I do need to keep track of time to evaluate the "duration" of a specific action triggered by the user) I tried with the following datatypes with no luck: **date** **datetime** **real** **double** Which datatype should I use when creating the shp in order for it to be filled with that specific value? Thanks in advance, m.


title: using PostGIS to calculate distance of a route over a large area

text: I'm new to PostGIS and looking to measure the distance of the Pacific Crest Trail which spans western North America from Canada to Mexico. I have a data set containing the lat lon and dates from GPX data taken on the trail from a hiker. I also have a shapefile of the entire trail that's fairly high resolution. What I'd like to do is to extract the total mileage gained for each of the GPS locations from computing the length of the trail data and then intersecting the two. I think the steps to doing this as are follows:   1. process the GPS location data to get lat lon positions at weekly intervals with additional special events like summiting Mt Whitney, trail start, end, etc.(this part is finished)   2. snap the lat lon GPS points to the nearest PCT shapefile nodes.   3. split the PCT shapefile into segments based on an intersection with the GPS points.   4. convert the total mileage from meters to miles.   5. total the distance gained by each GPS point as the dates increase. Some questions I have are: What projection or method of measuring distance would be most accurate for such a large north - south distance? I'm thinking of transforming the data into a custom equidistant projection but am unsure about setting the parameters for yielding the most accurate distance measurement. Or does PostGIS have other functionality that would yield accurate distance measurements without having to transform the data? Does this sound like a reasonable approach? Would you have any suggestions or tips that might help? Any advice you might have would be appreciated. Thanks in advance


title: How to set ArcGIS for Desktop license to lower level for testing whether Python script still runs?

text: I have a Python script that will need to be given to a client. We are currently evaluating what license level the client will need in order to run the script for their software acquisition purposes. We all have ArcGIS 10.1 Advanced license in our company. I checked all the tools and extensions used in the script against the web help and found that they should work on the Standard license level. Is there a way I can test the script against the Standard license level, when we have no Standard license level in or company (only the higher Advanced level)?


title: External geojson not showing in OpenLayers

text: I'm new to OpenLayers. I tried some recipes found here to load an external geojson file as an overlay on a map I'm creating and it's not working. The map is loading as well as my base layer, but no overlay... Anyone could help?                     Both file, trains and gares.geojson are in the same folders as my html. This is an extract from gares.geojson (in case the problem's there):               {     "type": "FeatureCollection",          "features": [     { "type": "Feature",      "properties": { "NO_STATION": 101.0, "NOM": "Gare Lucien-L'Allier", "NoAbord08": 18, "Parking": 0 },      "geometry": { "type": "Point", "coordinates": [ -73.570806, 45.494761 ] } }     ,     { "type": "Feature",      "properties": { "NO_STATION": 214.0, "NOM": "Gare Delson", "NoAbord08": 41, "Parking": 1 },      "geometry": { "type": "Point", "coordinates": [ -73.540803, 45.366832 ] } }          ]     }      I literally copy paste those two recipes and only changed the data source. I don't understand what's wrong... Thanks!


title: How to write Python script for spatial interpolation?

text: I am new to Python. I have 5 rainfall stations in a basin of 1400km2, and I want to interpolate the daily rainfall series for 14 years. I spent a long time on this, but am unable to write a program in Python for it. Please guide me on how to solve this problem I use ArcGIS 10.


title: Wrong $area from field calculator (QGIS)

text: The area of the polygons in a layer, calculated by the field calculator with "$area" is wrong (much too big) for some polygons. If I draw the concerned polygon again, the area is calculated correctly. It looks like the polygons really have a wrong size, but the centroids are calculated correctly. edit: Thanks for your replies so far. Here is an example with some polygons. For some, the area is calculated right, for some not (as I said much too big). For example the polygon with the ID 669. If I redraw it, the area is calculated right. The problem is, I have >1400 polygons, and I don't have the time to check and redraw every single one. Maybe someone got an explanation/solution?


title: Using st_intersects caused mismatch when a marking (tagging) shape marked a polygon overlapping only on the edge

text: I have accidentally marked a landshape as being tagged "CU GO-M" because it seems to share an _edge_ with a tagshape for "CU GO-M" but nothing else. I used the st_intersects method to join these shapes erroneously. Which method should I use? This might be a dumb question because I think it is st_within, but I'm worried that might kick out some false negatives (rather than a false- positive, ie this case). To summarize, I think any interior overlap match should be marked as tagged. This case marks the tagging for only an exterior (edge) intersect. ![4816 Wendover has been tagged within this GO-M tagging district](http://i.stack.imgur.com/rGZE9.jpg)


title: Sources for Geo-referenced data for Tiger populations in India or other countries

text: I'm currently an Undergrad student in a GIS class, and I'm trying to conduct some research on tiger populations and urbanization. Does anyone know where I might be able to find data for tiger populations?


title: Export basemap and vector layers into a single image with OpenLayers 2.11?

text: Has anyone been able to export OpenLayers 2.11 map and features as a static image? This example is exactly what I'd like to do, but it's using OL3 alpha which is not an option for me... http://dev.openlayers.org/sandbox/camptocamp/canvas/openlayers/examples/exportMapCanvas.html Thanks! Vincent


title: What are the top best webmapping API clients?

text: > **Possible Duplicate:**   >  Comparison of JavaScript mapping libraries? I am looking for the top APIs that I can use to display maps from a number of map servers/sources. (like openlayers and openfscales) I already know OpenSlales and Openlayers. So are there any other popular APIs for other languages (C++,Java,... ) ?


title: Map stack preferences

text: > **Possible Duplicate:**   >  Comparison of JavaScript mapping libraries? I've spent most of my GIS life doing analysis and 3D stuff and after skirting the edges of web-mapping for a while (using OpenLayers and dabbling with Mapnik), I am now getting a bit more serious about it. There seem to be so many possible routes each with pros and cons so I'm trying to get a handle on what constitutes a robust map-stack (Python APIs are a must for me). What do people use/recommend and why? Many thanks.


title: What technology to use for an Interactive web map?

text: > **Possible Duplicate:**   >  Comparison of JavaScript mapping libraries? I need to build a custom interactive web map of my country with all it's regions (Peru). Which tool or technology would be the best way to go? The idea is to have a map that shows the population of each region with a color scale: the more people the region has, the darker it's painted. Also I need to display some information as text and pictures when I click each region as the maps zooms in. I've seen around and there are lots of "pre-developed" tools that can be customized but I need to build my own. Flash is already discarded. I was thinking of Java applet but is just an idea. Thank you for your answers. **System requirements: Apache+PHP+MySQL**


title: Determining correct Coordinate System DHDN Bessel 1841

text: I am working on a data set and coordinates are supposed to be in DHDN Bessel 1841. Have a look here: https://www.pegelonline.wsv.de/gast/stammdaten?pegelnr=2330013 I created a shapefile with the cordinates. But I need to define the projection. I am using qgis and when I checked DHDN options, there were a lot like DHDN/ Gauss-Kruger 2,3,4,5 or DHDN/ 3-degree gauss zone 1,2,3,4,5. I am confused which one to use. Some insight will be of great help. Thanks NS


title: Selecting additional features in ArcMap 10

text: I'm having trouble selecting features in ArcMap 10. I have a shapefile of points, and I am in editing mode. I want to drag a box to select some points, then drag a new box and have these points added to my selection. Holding down shift while creating the new box will invert the selection - so any previously unselected point are added to the selection, but any already selected points that are in the overlap are deselected. This makes selecting a large number of points difficult, as I must zoom in each time to ensure that there is no overlap. How can I simply add points to my selection?


title: How to force network analyst to use metero station?

text: I'm new to use multimodal network analyst .I use esri tutorial for multimodal network analyst.It was very interesting for me .What i need is that i want to force network analyst to use metero station for network analyst .for example in esri tutorial multimodal network i don't want to route from street to street like the below image and i'd like to use metero station.How can i force network analyst to use metero station ? thanks for advance![enter image description here](http://i.stack.imgur.com/2uX09.png)


title: Will MapInfo only read PostGIS data that has been loaded via EasyLoader?

text: NMZ says in this thread that MapInfo will only read data in PostGIS that has been loaded by MapInfo's EasyLoader. I wasn't sure if this was still the case so I've loaded data into postgis via ogr2ogr but indeed MapInfo 12 hits a problem:               ODBC Error: ODBC RC=-1, ODBC SQLState=42601, DBMS RC=7, DBMS Msg=ERROR: zero-length delimited identifier at or near """"; Error while executing the query.      In the postgresql log I can see MapInfo attempting to execute a select statement;               2014-02-16 10:36:36 UTC ERROR:  zero-length delimited identifier at or near """" at character 74     2014-02-16 10:36:36 UTC STATEMENT:  select "ogc_fid", "app_ref", "id", "sdw"."planning"."polygons_ogr_28777"."", public.ST_AsBinary("wkb_geometry") as "wkb_geometry" from "sdw"."planning"."polygons_ogr_28777"      That ."" on the polygons table looks pretty weird to me. Taking a look at the query MapInfo is asking of PostGIS on an easyloader table looks like this;               2014-02-16 10:36:21 UTC LOG:  statement: select "app_ref", "id", "MI_PRINX", "MI_STYLE", public.ST_AsBinary("SP_GEOMETRY") as "SP_GEOMETRY" from "sdw"."planning"."polygons_bng5365"      No ."" in this one. I'm wondering if MapInfo is particularly to MI_PRINX being the primary key. I've got to do some further digging but would appreciate any insights people have. I've pasted the fuller log file here - loading the two tables via ogr (polygons_ogr_28777) and easyloader (polygons_bng5365) and then accessing the pair via MapInfo. The numbers in the table names do not reflect any srid btw. Just numbers.. Update: Sorry, I forgot to mention that I've built a mapinfo.mapinfo_mapcatalog. There's a trigger function which populates this table when an insert is made to geometry_columns. (PG1.5). The entries for the two tables are;               "spatialtype";"tablename";"ownername";"spatialcolumn"     ;"db_x_ll";"db_y_ll";"db_x_ur";"db_y_ur"     ;"view_x_ll";"view_y_ll";"view_x_ur";"view_y_ur"     ;"coordinatesystem"     ;"symbol"     ;"xcolumnname";"ycolumnname"     ;"renditiontype";"renditioncolumn";"renditiontable";"number_rows"          19.2;"polygons_bng5365";"planning";"SP_GEOMETRY"     ;487932.174983;162961.282073;487976.574637;163000.880738     ;487932.174983;160574.96515;490597.036446;163000.880738     ;"Earth Projection 8, 79, ""m"", -2, 49, 0.9996012717, 400000, -100000 Bounds (-7845061.1011, -15524202.1641) (8645061.1011, 4470074.53373)"     ;""     ;"NO_COLUMN";"NO_COLUMN"     ;1;"MI_STYLE";"";          19.3;"polygons_ogr_28777";"planning";"wkb_geometry"     ;0;0;700000;1300000     ;0;0;700000;1300000     ;"Earth Projection 8, 79, ""m"", -2, 49, 0.9996012717, 400000, -100000 Bounds (-7845061.1011, -15524202.1641) (8645061.1011, 4470074.53373)"     ;"SYMBOL (35,0,12) PEN (1,2,0) BRUSH (2,16777215,16777215) "     ;"NO_COLUMN";"NO_COLUMN"     ;1;"";"";      The mapinfo.mapinfo_mapcatalog definition is as follows;               -- DROP TABLE mapinfo.mapinfo_mapcatalog;          CREATE TABLE mapinfo.mapinfo_mapcatalog     (       spatialtype double precision,       tablename character varying(32),       ownername character varying(32),       spatialcolumn character varying(32),       db_x_ll double precision,       db_y_ll double precision,       db_x_ur double precision,       db_y_ur double precision,       view_x_ll double precision,       view_y_ll double precision,       view_x_ur double precision,       view_y_ur double precision,       coordinatesystem character varying(254),       symbol character varying(254),       xcolumnname character varying(32),       ycolumnname character varying(32),       renditiontype integer,       renditioncolumn character varying(32),       renditiontable character varying(32),       number_rows integer     )     WITH (       OIDS=FALSE     );     ALTER TABLE mapinfo.mapinfo_mapcatalog OWNER TO postgres;     GRANT ALL ON TABLE mapinfo.mapinfo_mapcatalog TO postgres;     GRANT ALL ON TABLE mapinfo.mapinfo_mapcatalog TO public;          -- Index: mapinfo.mapcatalog_idx          -- DROP INDEX mapinfo.mapcatalog_idx;          CREATE UNIQUE INDEX mapcatalog_idx       ON mapinfo.mapinfo_mapcatalog       USING btree       (tablename, ownername);      The MapInfo manual also says to create a MapInfo user account which I have done though not sure of its worth.


title: Is GeoServer Image Pyramid upscaling?

text: I have an _Image Pyramid_ in _GeoServer 2.2 RC3_ (not the latests - I know, but I have other pyramids working here). The total size is around 1TB and file size is around 1-2GB for each GeoTIFF. When looking at the results using WMS, it all looks blurred. Going into the GeoServer log, I found that when I request an image, it selects the level which I think is causing it to upscale the result to fit the screen: Request:   Size in pixels: 1632x960   Size in meters: 1406x827   I have a _GSD = 1 meter/pixel_ at the level GeoServer selects (which is 3). Shouldn't it select level 2 with a _GSD = 0.5 meter/pixel_ to ensure GeoServer downsamples the result - thus loosing information instead of creating it?


title: How can I join a spreadsheet to a point layer?

text: I have a layer with point data for the grave positions in our cemetery. We are gathering data from each headstone. If I enter the data to a spreadsheet can I later associate the spreadsheet with the point data layer. Will I have to go back to the point data layer and add all the columns from the spreadsheet as attributes to make the point layer searchable by data from the headstones? Would a better solution be to use a PostgreSQL to make a database to associate with the point layer? Am I going about this the right way?


title: Problem with QGIS unable to access virtual drive

text: After a motherboard failure I have had to reinstall QGIS Dufour on a Windows 7 laptop. As I share data with a desktop machine I use LiveDrive briefcase as a virtual drive. With all versions of QGIS I have always been able to access LiveDrive which appears as a mapped drive M: However, on this laptop, having installed PostgreSQL, PostGIS and Osgeo4W shell to be able to load MasterMap data using the Astun Loader Astun Loader I am now unable to access drive M: from QGIS, GRASS or SAGA. I did have some access issuses and I had use 'run as adminstrator' to get round it even though my account has administrator rights. I have tried doing a complete reinstall of QGIS without success. I don't know whether this a QGIS issue or not. Any suggestions to fix this are welcome.


title: Polygonize(raster to vector)

text: I have been install QGIS 2.0.1.when I use Polygonize (Raster to vector)for convert jpg to shp I get error message such as: > C:\windows\123.exe is not recognized as an internal or external command, > operable program or batch file. Please help me to fix it. thanks


title: How do I get the area of a WGS84 polygon in square meters?

text: The following gives me the area of the polygon in square degrees:               SELECT ST_Area(ST_GeomFromText('POLYGON((871325.790874952 6105405.3261047,871418.748307692 6105359.72944624,871346.22022442 6105215.141258,871254.85408906 6105261.72007212,871325.790874952 6105405.3261047))',4326)) As sqm;      Result: 16586.2319335938 How do I get the area in square meters? I tried the following:               SELECT ST_Area(ST_Transform(ST_GeomFromText('POLYGON((871325.790874952 6105405.3261047,871418.748307692 6105359.72944624,871346.22022442 6105215.141258,871254.85408906 6105261.72007212,871325.790874952 6105405.3261047))',4326)),31467) As sqm;      But it seems like ST_Transform doesn't work. I get the following error message:               ERROR:  function st_transform(geometry) does not exist     LINE 1: SELECT ST_Area(ST_Transform(ST_GeomFromText('POLYGON((871325...                            ^     HINT:  No function matches the given name and argument types. You might need to add explicit type casts.     ********** Error **********      ERROR: function st_transform(geometry) does not exist     SQL Status:42883     Hint:No function matches the given name and argument types. You might need to add explicit type casts.      Missing Proj support shouldn't be the problem, with SELECT PostGIS_full_version() I get               "POSTGIS="1.5.3" GEOS="3.3.3-CAPI-1.7.4" PROJ="Rel. 4.8.0, 6 March 2012" LIBXML="2.8.0" USE_STATS"


title: Creating new lines along specifically selected existing lines using QGIS?

text: I have tried to find a solution for this everywhere, with no result... My problem is that I want to be able to manually draw a line that automatically snaps to a line (creating an exact copy of the other lines geometry) that's on another layer (using something like this: How to create a line with start and end points which are in another layer in QGIS?). See example in attached image. This is how it would look if just using the snapping option. I will loose a few percent of length in the new lines compared to the lines in the other layer, and will later have to compensate that loss of length in some calculation (something I'd rather not do). I need to manually select which lines that should be created, therefore I can't use any automatic method that creates lines for the whole map (if such should exist). Is there any way to alter the lines afterwards, to match the other layers' lines? Or perhaps a plugin tool that I haven't yet found, that can help me in the manual creation of these lines? ![Example of snapping line](http://i.stack.imgur.com/dXVvV.png)


title: Finding left and right side of line in arcgis

text: i have road layer as line and polygon layer as building,i want to know left and right side of road in relation with buildings is it possible in Qgis or Arcgis?![enter image description here](http://i.stack.imgur.com/z6axD.jpg)


title: Calling QgsProjectionSelector from QGIS Python plugin?

text: I am trying to create a dialog that allows the user to select a coordinate system from either the recently used projections or a general list of all available projections. This would basically offer the same functionality as the standard dialog that pops up when you load in a shapefile with no known CRS. Is there a way to call the "Coordinate Reference System Selector"-dialog from within a plugin? I know it is possible to query the recently used projections and all available projections, but it seems like a lot of overhead to recreate something that is already there in QGIS.


title: glm with raster data: species decline~ ndvitrend

text: I have species decline data (number of population per year, over a period of 23 years) out of ndvi data for these 23 years I calculated the trend, and mapped it. now i'd like to relate the negative ndvi trend to species decline to see if this negative ndvi trend can explain the species decline. However, my trend data are raster data and my species decline data are just numbers in a csv table. I'd like to perform sth like this: lm1<- lm(speciesdecline~ndvitrend) in R. But i dont know how to do this with 2 different kinds of data. Anyone an idea? thanks!


title: ArcGIS Toolbox, parameters for python script not as positional arguments

text: I have an ArcGIS Toolbox with a bunch of parameters. A python script is used to perform some raster operations, based on this parameters. I am using Python argparse in the script to get the parameters. I do this because the program was originally developed as a command line tool using ArcPy. With argparse I can show some useful help text and use command line switches where the position of the argument in the argument list does not matter. Of course argparse offers much more advantages. ArcGIS is passing the parameters as positional arguments. So I can't use parameters like -k or --doSomeThing=100 Using a checkbox for example, that is ticked, this will pass a "true" to my script but I would like to have a "real" command line switch like "-k" or something like that. Has anyone an idea if this is possible? Maybe using tool validation?


title: QGIS 2: UUID generator edit widget only populating attribute with UUID if through attribute pop-up window enabled

text: I am trying to get the UUID generator edit widget to populate a UUID attribute when the attribute pop-up window is suppressed. Any ideas? (Seems I can only the get the auto-populate to happen when the attribute pop-up window is not suppressed)


title: using postgis function in pyqgis

text: I'd use Postgis function in Pyqgis script. This code:               qsql="(select gid, id, St_Buffer(the_geom,200) as geom from archi)"     uri = QgsDataSourceURI()     uri.setConnection("ip", "port", "x", "x", "x")     uri.setDataSource("",qsql,  "the_geom","","gid")     layer = QgsVectorLayer(uri.uri(), "test_buffer", "postgres")     if not layer.isValid:          print 'error'     else:          QgsMapLayerRegistry.instance().addMapLayer(layer)      load Vector layer, with NO buffer. What's the way to run Postgis Function in PyQgis? Thanks f


title: ogr2ogr/python overwrite PostGIS table

text: Using ogr2ogr command line I can overwrite or append to an existing PostGIS table like:               ogr2ogr -f "PostgreSQL"-overwrite PG:"host=myserver port=54321 user=user dbname=mydbname password=mypassword" C:/Temp/my.shp      Next, I would like to incorporate this logic into a python script. I've found some useful sites, however I'm not finding any examples using these methods. Open Source RS/GIS Python Week 1 PostGIS OGR Examples Append a shapefile to a postgis table using the GDAL/OGR Python interface How do you overwrite or append to a PostGIS table using ogr2ogr commands in python?


title: Arc GIS 10.2 Relating Tables in an OLE DB

text: I have a large registry of about 260,000 records in microsoft access with data on land quarter sections. This data was to much for Arc GIS to handle directly, so I created a Microsoft Jet 4.0 OLE DB connection to the database and added the query table I greated in Access. Now I want to relate the table with the quarter section feature class (polygons of all the quarter sections). With this, I need to be able to select one or more quarter sections on the feature class and bring up all the related data from my access table. However, Arc GIS help states that "Because OLE DB connections do not have an ObjectID field, which is a field maintained by ArcGIS that guarantees a unique ID for each row in the table, you won't be able to perform certain operations in ArcMap. These include selecting the features in the layer on the map in any way, applying a definition query, creating full relates, or editing the attributes. See About ObjectID fields for more information about working with tables that lack Object IDs." So I can not relate a table from the OLE connection with Arc GIS due to lack of ObjectID. Does anyone know a work around or a way to create the ObjectID so I can select related records from my large access dataset based on selected quarter sections on my feature class. Thanks!


title: CartoCSS detect odd vs even for use in labeling conditional

text: I liked the part in the advanced label styling about randomly asigning a 0 or 1 to control label placement and then using as a conditional               #labels {       text-name: "[name]";       text-face-name: "OpenSans Regular";       text-placement-type: simple;       text-placements: "E,NE,SE,W,NW,SW";       [dir=1] { text-placements: "W,NW,SW,E,NE,SE"; }     }      Rather than create a new variable, could I just use my FID and detect whether or not it is odd/even? Is there a way to detect this property in the CartoCSS?


title: Using da.SearchCursor to get count of unique values in a field

text: Is it possible to use arcpy.da.SearchCursor() to get a count of unique features in a field without a dictionary? So, by this I mean if I have a field called TEXT which contains 10,000 rows populated with either A or B. Can I use the da.SearchCursor to quickly give me the count for the total amount of A in TEXT and B in the TEXT? I know it is possible to use Summary Statistics or the Get Count tool but these involve having to make a feature layer and then selecting features or kicking out a table. I looking to have a few steps as possible in my script :-)


title: Is there a way of splitting roads at intersections with contours in QGIS?

text: I'm a beginner with QGIS. I need to split a polylines showing public roads at each intersection with a table containing elevation data (contours) which is also in polyline format. I will also need to add two columns to the roads table to record minimum and maximum elevation for each new record. I would appreciate any advice on how to go about this.


title: Transform projection from EPSG:4326 into EPSG:900913

text: . i've got this problem about projection. .i just wanted to know what am i missing. . i already displayed map base layer using openlayers through               base_layer =  new OpenLayers.Layer.WMS( "OpenLayers WMS",                "http://vmap0.tiles.osgeo.org/wms/vmap0?", {layers: 'basic'} );      i also displayed my vector layer through               var layer = new OpenLayers.Layer.Vector("Butuan");         map.addLayer(layer);         var format = new OpenLayers.Format.GeoJSON();         var feat = format.read({{my_geojson|safe}});         layer.addFeatures(feat);      and my position is               map.setCenter(new OpenLayers.LonLat(125.54915811889646, 8.90855351100384),11);      which positions the view as i desired it to be. when i tried to change my base layer to               base_layer =  new OpenLayers.Layer.OSM("OpenStreetMap");      the center of the map is redirected to another location. .i guess that this problem pertains to the projection. i tried to change the center into               map.setCenter(new OpenLayers.LonLat(125.54915811889646, 8.90855351100384),11).transform(geographic, mercator);      which is defined as               var geographic = new OpenLayers.Projection("EPSG:4326");            var mercator = new OpenLayers.Projection("EPSG:900913");      no changes has been noticed. .the position of the center is still redirected to another location. . can anyone suggest solutions on how to solve this problem?? thnx in advance. .you are really appreciated. . by the way, i am using django framework. . vectorformats and simplejson to convert the retrieved spatial data into geojson from postgis. .


title: Why does QGIS 2.2.0 displaying triangles instead of rectangles?

text: When displaying the demographic datas from INSEE (http://www.insee.fr/fr/themes/detail.asp?reg_id=0&ref_id=donnees- carroyees&page=donnees-detaillees/donnees-carroyees/donnees- carroyees-200m.htm), the results under qgis 2.2.0 is quite surprising :   * instead of rectangles I just see triangles in the bottom right corner (so covering half the surfaces)   * the same shape file opened with MapInfo 10.0 seems correct with all the rectangles correctly displayed and filled. I can't see the reason of such a difference, it seems for me to be a bug of Qgis ..


title: Refreshing MSSQL Layer in QGIS

text: When working on a map in QGIS I will often set up a view in my SQL Server spatial database and create a MSSQL layer based on this. As part of my workflow I often make adjustments to the SQL view based on input from the client. Is there a quick way of refreshing the layer in QGIS so that it picks up on these changes? Currently I end up closing and re-opening my project to pull any changes through.


title: Identify geographic location to points

text: I am quite new in QGIS and I really need your help. I have a map of Kenya divided into 68 districts and I have coordinates of 30.000 schools. The points have one variable besides the latitude and longitude data called "distance to nearest school". First I imported the map. Then I imported the spreedsheet and transformed it to a SHP file like the map. So now I have two SHP files. The points have three variables: "Latitude", "Longitude" and "Distance to nearest school". I want to have a fourth variable called "District" and take the mean of the "Distance to nearest school" for each district.


title: How can I get an error map of an interpolation in ArcGIS?

text: I have about 10 point shapefiles containing a lot of points each. (I would prefer not to join them into one shapefile, it's pain). And I am using interpolation method Topo to Raster to create a DTM. This method can deal with many separate datasets. I'm using ArcGIS 9.3 (and have access to 10). I would like to get errors of the estimates, to check if they have normal distribution. How can I get the interpolation errors? Like the -e flag option in GRASS. Is it available only for kriging? Is it the only way, to compute them manually comparing each dataset to the resulting raster? Thank you for advices.


title: Why does my zoomToExtent zoom to 0,0?

text: In this map I have added a line of code to make it zoom to the vector layer:               map.zoomToExtent(vector_Layer.getDataExtent());      This zooms in to the maximum zoom at the point in the Gulf of Guinea where the equator meets the Greenwich meridian. It is not supposed to zoom here, it is should show Scotland. I've passed the same bit of code in the console in Chrome and that works fine.   Does anyone know why it is not working properly when the page loads and does anyone know what I can do to fix it? Thanks.


title: Installing raster resampling module in QGIS 1.8

text: On the website of Sourcepole, I have recently come across the case of resampling images in Quantum GIS. Because I am rather a beginner in GIS software I would like to ask you for information about how I can enable/install the resampling option (as shown here) on QuantumGIS 1.8., (perhaps it is obvious, but not for me). I would be very grateful for any information.