title: writeOGR column limit error: "Creating Name field failed"

text: I have a SpatialPointsDataFrame with 1057 elements and 4760 columns and want to export it to a shapefile. When I try to do this with writeOGR() I get an error:               > writeOGR(pts,dsn="/home/joana/Desktop/",layer="PT", driver="ESRI Shapefile")     Error in writeOGR(pts, dsn = "/home/joana/Desktop/", layer = "PT", driver = "ESRI Shapefile") :        Creating Name field failed      When I reduce the column number, everything works fine.But I really need all of the columns, is there any way to work around it? I'm using Ubuntu 12.04/R/QGIS


title: How create one single coverage polygon feature class from XY point data?

text: I'm using ArcMap 10.2, and I need to map areas covered by London postcodes (E1 0xx, E1 1xx, E1 2xx, and so on). What I have is the coordinates for every single postcode (310,000 in total) in an Excel sheet, so I can create point feature classes, which I then want to group as explained above. My question is, _how can I create a polygon shapefile that will include a group of coordinates within E1 0xx postcode for example. I want that polygon shapefile to cover the are of those postcodes specifically. In the end I need to end up with multiple polygon shapefiles that will represent different Postcodes districts (E1 0xx, E1 1xx...)._ Tried the Coverage tool, but it's not seem to be doing what I expected it to. Just installed XTools Pro, maybe that could help me. Help is much appreciated. Thanks


title: How to plot geo-data using matplotlib/python

text: I am trying to plot a polygone on python, using different libararies, but no one of these worked with me. I tried vincent ,Shapely , .. and now matplotlib. I don't wanna use softwares like qgis I used this example: http://geographika.co.uk/visualising-geojson- in-15-minutes               from matplotlib import pyplot     from descartes import PolygonPatch     import math     import urllib2     import simplejson          data = []     with open('file.json') as f:         for line in f:             data.append(json.loads(line))               def plotFeature(coordlist, myplot):         #create a polygon geojson-like feature         poly = {"type": "Polygon", "coordinates": coordlist}         patch = PolygonPatch(poly, fc='#6699cc', ec='#6699cc', alpha=0.5, zorder=2)         #plot it on the graph         myplot.add_patch(patch)               for coordlist in data['features'][0]['geometry']['coordinates']:         plotFeature(coordlist, myplot)          fig = pyplot.figure(1, figsize=(10, 4), dpi=180)     ax = fig.add_subplot(121)     ax.savefig('myplot.png')      External file "file.json"               var states = [{         "type": "Feature",         "properties": {"party": "Republican"},         "geometry":          {         "type": "Polygon",         "coordinates": [             [                 [                     -84.32281494140625,                     34.9895035675793                 ],                 [                     -84.29122924804688,                     35.21981940793435                 ],                 [                     -84.24041748046875,                     35.25459097465022                 ],                 [                     -84.22531127929688,                     35.266925688950074                 ],                 [                     -84.20745849609375,                     35.26580442886754                 ],                 [                     -84.19921875,                     35.24674063355999                 ],                 [                     -84.16213989257812,                     35.24113278166642                 ],                 [                     -84.12368774414062,                     35.24898366572645                 ],                 [                     -84.09072875976562,                     35.24898366572645                 ],                 [                     -84.08798217773438,                     35.264683153268116                 ],                 [                     -84.04266357421875,                     35.27701633139884                 ],                 [                     -84.03030395507812,                     35.291589484566124                 ],                 [                     -84.0234375,                     35.306160014550784                 ],                 [                     -84.03305053710936,                     35.32745068492882                 ],                 [                     -84.03579711914062,                     35.34313496028189                 ],                 [                     -84.03579711914062,                     35.348735749472546                 ],                 [                     -84.01657104492188,                     35.35545618392078                 ],                 [                     -84.01107788085938,                     35.37337460834958                 ],                 [                     -84.00970458984374,                     35.39128905521763                 ],                 [                     -84.01931762695312,                     35.41479572901859                 ],                 [                     -84.00283813476562,                     35.429344044107154                 ],                 [                     -83.93692016601562,                     35.47409160773029                 ],                 [                     -83.91220092773438,                     35.47632833265728                 ],                 [                     -83.88885498046875,                     35.504282143299655                 ],                 [                     -83.88473510742186,                     35.516578738902936                 ],                 [                     -83.8751220703125,                     35.52104976129943                 ],                 [                     -83.85314941406249,                     35.52104976129943                 ],                 [                     -83.82843017578125,                     35.52104976129943                 ],                 [                     -83.8092041015625,                     35.53446133418443                 ],                 [                     -83.80233764648438,                     35.54116627999813                 ],                 [                     -83.76800537109374,                     35.56239491058853                 ],                 [                     -83.7432861328125,                     35.56239491058853                 ],                 [                     -83.71994018554688,                     35.56239491058853                 ],                 [                     -83.67050170898438,                     35.569097520776054                 ],                 [                     -83.6334228515625,                     35.570214567965984                 ],                 [                     -83.61007690429688,                     35.576916524038616                 ],                 [                     -83.59634399414061,                     35.574682600980914                 ],                 [                     -83.5894775390625,                     35.55904339525896                 ],                 [                     -83.55239868164062,                     35.56574628576276                 ],                 [                     -83.49746704101562,                     35.563512051219696                 ],                 [                     -83.47000122070312,                     35.586968406786475                 ],                 [                     -83.4466552734375,                     35.60818490437746                 ],                 [                     -83.37936401367188,                     35.63609277863135                 ],                 [                     -83.35739135742188,                     35.65618041632016                 ],                 [                     -83.32305908203124,                     35.66622234103479                 ],                 [                     -83.3148193359375,                     35.65394870599763                 ],                 [                     -83.29971313476561,                     35.660643649881614                 ],                 [                     -83.28598022460938,                     35.67180064238771                 ],                 [                     -83.26126098632811,                     35.6907639509368                 ],                 [                     -83.25714111328125,                     35.69968630125201                 ],                 [                     -83.25576782226562,                     35.715298012125295                 ],                 [                     -83.23516845703125,                     35.72310272092263                 ],                 [                     -83.19808959960936,                     35.72756221127198                 ],                 [                     -83.16238403320312,                     35.753199435570316                 ],                 [                     -83.15826416015625,                     35.76322914549896                 ],                 [                     -83.10333251953125,                     35.76991491635478                 ],                 [                     -83.08685302734375,                     35.7843988251953                 ],                 [                     -83.0511474609375,                     35.787740890986576                 ],                 [                     -83.01681518554688,                     35.78328477203738                 ],                 [                     -83.001708984375,                     35.77882840327371                 ],                 [                     -82.96737670898438,                     35.793310688351724                 ],                 [                     -82.94540405273438,                     35.820040281161                 ],                 [                     -82.9193115234375,                     35.85121343450061                 ],                 [                     -82.9083251953125,                     35.86902116501695                 ],                 [                     -82.90557861328125,                     35.87792352995116                 ],                 [                     -82.91244506835938,                     35.92353244718235                 ],                 [                     -82.88360595703125,                     35.94688293218141                 ],                 [                     -82.85614013671875,                     35.951329861522666                 ],                 [                     -82.8424072265625,                     35.94243575255426                 ],                 [                     -82.825927734375,                     35.92464453144099                 ],                 [                     -82.80670166015625,                     35.927980690382704                 ],                 [                     -82.80532836914062,                     35.94243575255426                 ],                 [                     -82.77923583984375,                     35.97356075349624                 ],                 [                     -82.78060913085938,                     35.99245209055831                 ],                 [                     -82.76138305664062,                     36.00356252895066                 ],                 [                     -82.69546508789062,                     36.04465753921525                 ],                 [                     -82.64465332031249,                     36.060201412392914                 ],                 [                     -82.61306762695312,                     36.060201412392914                 ],                 [                     -82.60620117187499,                     36.033552893400376                 ],                 [                     -82.60620117187499,                     35.991340960635405                 ],                 [                     -82.60620117187499,                     35.97911749857497                 ],                 [                     -82.5787353515625,                     35.96133453736691                 ],                 [                     -82.5677490234375,                     35.951329861522666                 ],                 [                     -82.53067016601562,                     35.97244935753683                 ],                 [                     -82.46475219726562,                     36.006895355244666                 ],                 [                     -82.41668701171875,                     36.070192281208456                 ],                 [                     -82.37960815429686,                     36.10126686921446                 ],                 [                     -82.35488891601562,                     36.117908916563685                 ],                 [                     -82.34115600585936,                     36.113471382052175                 ],                 [                     -82.29583740234375,                     36.13343831245866                 ],                 [                     -82.26287841796874,                     36.13565654678543                 ],                 [                     -82.23403930664062,                     36.13565654678543                 ],                 [                     -82.2216796875,                     36.154509006695                 ],                 [                     -82.20382690429688,                     36.15561783381855                 ],                 [                     -82.19009399414062,                     36.144528857027744                 ],                 [                     -82.15438842773438,                     36.15007354140755                 ],                 [                     -82.14065551757812,                     36.134547437460064                 ],                 [                     -82.1337890625,                     36.116799556445024                 ],                 [                     -82.12142944335938,                     36.10570509327921                 ],                 [                     -82.08984375,                     36.10792411128649                 ],                 [                     -82.05276489257811,                     36.12678323326429                 ],                 [                     -82.03628540039062,                     36.12900165569652                 ],                 [                     -81.91268920898438,                     36.29409768373033                 ],                 [                     -81.89071655273438,                     36.30959215409138                 ],                 [                     -81.86325073242188,                     36.33504067209607                 ],                 [                     -81.83029174804688,                     36.34499652561904                 ],                 [                     -81.80145263671875,                     36.35605709240176                 ],                 [                     -81.77947998046874,                     36.34610265300638                 ],                 [                     -81.76162719726562,                     36.33835943134047                 ],                 [                     -81.73690795898438,                     36.33835943134047                 ],                 [                     -81.71905517578125,                     36.33835943134047                 ],                 [                     -81.70669555664062,                     36.33504067209607                 ],                 [                     -81.70669555664062,                     36.342784223707234                 ],                 [                     -81.72317504882812,                     36.357163062654365                 ],                 [                     -81.73278808593749,                     36.379279167407965                 ],                 [                     -81.73690795898438,                     36.40028364332352                 ],                 [                     -81.73690795898438,                     36.41354670392876                 ],                 [                     -81.72454833984374,                     36.423492513472326                 ],                 [                     -81.71768188476562,                     36.445589751779174                 ],                 [                     -81.69845581054688,                     36.47541104282962                 ],                 [                     -81.69845581054688,                     36.51073994146672                 ],                 [                     -81.705322265625,                     36.53060536411363                 ],                 [                     -81.69158935546875,                     36.55929085774001                 ],                 [                     -81.68060302734375,                     36.56480607840351                 ],                 [                     -81.68197631835938,                     36.58686302344181                 ],                 [                     -81.04202270507812,                     36.56370306576917                 ],                 [                     -80.74264526367186,                     36.561496993252575                 ],                 [                     -79.89120483398438,                     36.54053616262899                 ],                 [                     -78.68408203124999,                     36.53943280355122                 ],                 [                     -77.88345336914062,                     36.54053616262899                 ],                 [                     -76.91665649414062,                     36.54163950596125                 ],                 [                     -76.91665649414062,                     36.55046568575947                 ],                 [                     -76.31103515625,                     36.551568887374                 ],                 [                     -75.79605102539062,                     36.54936246839778                 ],                 [                     -75.6298828125,                     36.07574221562703                 ],                 [                     -75.4925537109375,                     35.82226734114509                 ],                 [                     -75.3936767578125,                     35.639441068973916                 ],                 [                     -75.41015624999999,                     35.43829554739668                 ],                 [                     -75.43212890625,                     35.263561862152095                 ],                 [                     -75.487060546875,                     35.18727767598896                 ],                 [                     -75.5914306640625,                     35.17380831799959                 ],                 [                     -75.9210205078125,                     35.04798673426734                 ],                 [                     -76.17919921875,                     34.867904962568744                 ],                 [                     -76.41540527343749,                     34.62868797377061                 ],                 [                     -76.4593505859375,                     34.57442951865274                 ],                 [                     -76.53076171875,                     34.53371242139567                 ],                 [                     -76.5911865234375,                     34.551811369170494                 ],                 [                     -76.651611328125,                     34.615126683462194                 ],                 [                     -76.761474609375,                     34.63320791137959                 ],                 [                     -77.069091796875,                     34.59704151614417                 ],                 [                     -77.376708984375,                     34.45674800347809                 ],                 [                     -77.5909423828125,                     34.3207552752374                 ],                 [                     -77.8326416015625,                     33.97980872872457                 ],                 [                     -77.9150390625,                     33.80197351806589                 ],                 [                     -77.9754638671875,                     33.73804486328907                 ],                 [                     -78.11279296875,                     33.8521697014074                 ],                 [                     -78.2830810546875,                     33.8521697014074                 ],                 [                     -78.4808349609375,                     33.815666308702774                 ],                 [                     -79.6728515625,                     34.8047829195724                 ],                 [                     -80.782470703125,                     34.836349990763864                 ],                 [                     -80.782470703125,                     34.91746688928252                 ],                 [                     -80.9307861328125,                     35.092945313732635                 ],                 [                     -81.0516357421875,                     35.02999636902566                 ],                 [                     -81.0516357421875,                     35.05248370662468                 ],                 [                     -81.0516357421875,                     35.137879119634185                 ],                 [                     -82.3150634765625,                     35.19625600786368                 ],                 [                     -82.3590087890625,                     35.19625600786368                 ],                 [                     -82.40295410156249,                     35.22318504970181                 ],                 [                     -82.4688720703125,                     35.16931803601131                 ],                 [                     -82.6885986328125,                     35.1154153142536                 ],                 [                     -82.781982421875,                     35.06147690849717                 ],                 [                     -83.1060791015625,                     35.003003395276714                 ],                 [                     -83.616943359375,                     34.99850370014629                 ],                 [                     -84.05639648437499,                     34.985003130171066                 ],                 [                     -84.22119140625,                     34.985003130171066                 ],                 [                     -84.32281494140625,                     34.9895035675793                 ]             ],             [                 [                     -75.69030761718749,                     35.74205383068037                 ],                 [                     -75.5914306640625,                     35.74205383068037                 ],                 [                     -75.5419921875,                     35.585851593232356                 ],                 [                     -75.56396484375,                     35.32633026307483                 ],                 [                     -75.69030761718749,                     35.285984736065735                 ],                 [                     -75.970458984375,                     35.16482750605027                 ],                 [                     -76.2066650390625,                     34.994003757575776                 ],                 [                     -76.300048828125,                     35.02999636902566                 ],                 [                     -76.409912109375,                     35.07946034047981                 ],                 [                     -76.5252685546875,                     35.10642805736423                 ],                 [                     -76.4208984375,                     35.25907654252574                 ],                 [                     -76.3385009765625,                     35.294952147406576                 ],                 [                     -76.0858154296875,                     35.29943548054543                 ],                 [                     -75.948486328125,                     35.44277092585766                 ],                 [                     -75.8660888671875,                     35.53669637839501                 ],                 [                     -75.772705078125,                     35.567980458012094                 ],                 [                     -75.706787109375,                     35.634976650677295                 ],                 [                     -75.706787109375,                     35.74205383068037                 ],                 [                     -75.69030761718749,                     35.74205383068037                 ]             ]         ]     } } ]       But it doesn't seem to work out. I need to plot these coordinates with a map on background. Any example, library, tutorial to recommend ?


title: calculate distances between different components in a (cable-) network

text: I have to calculate all distances between every different components in a cable network (for example the distance between an amplifier and a modem) and have to show the results in a diagram (frequency distribution) i have components like amplifier as a point layer, and the cable as a line- layer I use ArcGIS 9.3 any suggestions? Thx in advance!


title: Customizing Popups in Openlayers

text:
Is it possible to customize the shape of popups in OpenLayers? For example, I am creating a pop-up using the following approach:               var popup = new OpenLayers.Popup.FramedCloud(                     "popup",                     OpenLayers.LonLat.fromString(feature.geometry.toShortString()),                     null,                     '
'+feature.data.Name + '
test.url' + '
',                     null,                     true                 );      This creates a popup with the following appearance: ![enter image description here](http://i.stack.imgur.com/GWd16.png) Is it possible to change from the rounded corners and a curved 'connection' to a more square pop-up?


title: How to get specific coordinate system to get arcGIS point shapefile in r?

text: I have a Point shapefile from arcGIS 10. I want to get it in r and also want to get it a specific predefined coordinate systen for my data. I can read the shapefile in R but summary in r Shows projectltion info like this `proj4string : [NA]' and if i want to transform the reference systen it Shows something like this: Error in spTransform(xSP, CRSobj, ...) : No transformation possible from NA reference System What can I do to deal with this Problem. Thanks a lot:-)


title: loading .shp exercise data with windows 8 interface error

text: i am just starting the training manual. i have QGIS installed and the exercise data downloaded and unzipped. when i try to use the add vector dialog, it asks me for a folder and i browse to the exercise data folder. when i hit enter i get an error that says the folder name is not valid. the name is the unzipped name and i know the .shp file is there when i look with windows explorer.


title: Characters displaying correctly in FGDB with ArcMap 10.1 now displayed incorrectly with 10.2 -- font coding issue?

text: I work with data from a variety of European countries, and am having trouble getting the non-standard characters (e.g. Turkish alphabet) to display correctly in my geodatabase tables. Sure, i expected this, but when I was running ArcGIS 10.1, they all displayed just fine! I tried the registry edit workaround suggested on the ESRI website, to no avail. Anybody have any better suggestions for reverting all of my '?' characters to their correct letters? And no, i have no idea which character set they came from originally, nor which coding system they should be using.. no idea at all.


title: Count number of unique raster pixels in polygon using ArcMap?

text: I have a polygon shapefile and I want to generate a table of a count for each pixel value in each polygon. My main goal is to find out what pixel is most abundant in each polygon. I am using ArcMap 10.


title: Applying version control to ArcGIS Models

text: Version control is an indispensible tool for software development, allowing one to reliably and cleanly step back in time to the last time X did it's job exactly right, or to see what changed between then and now -- typically used when trying to figure out why X is no longer working exactly right. However all the tools I know about for this work only on plain text files. Toolboxes, and thus their models, are binary. Does anyone have a workable method to bring versioning to them? UPDATE: _Versioning_ in this context is different from _backup_. There are any number of simple methods to create snapshots of files for a particular date/time -- windows backup, `xcopy /s d:\foobar\ x:\foobar_%date%`, `zip stuff_%date%.zip stuff\\*`, .... Applying a tool such as git, fossil, mercurial, subversion, or ... to a binary file is one step better than using xcopy or zip in so much as one can add a commit message, _"Model foobar %date% now overwrites previous results only if Baz doesn't exist"_ , but is still anemic compared to what that same toolset can do applied to text files: e.g. show me exactly what was changed between last year and today. ![WinMerge screenshot](http://i.stack.imgur.com/hTZOV.png)


title: Homebrew OGR with FGDB support on OSX 10.9

text: This answer worked beautifully in 10.8, but isn't working immediately in 10.9. The link to download the FileGDB API still works, but there are a bunch of options available I don't remember from before. Do we want the Mac version, or the Linux one? Version 1.2 or 1.3? I've tried the Mac version 1.3 so far.


title: ArcSDE - Registering table with Geography spatial column fails with error -138 (SQL Server 2008)

text: I have a table of point features I am trying to register with ArcSDE (9.3.1 on SQL Server 2008) using the following command:               sdelayer -o register -l ThirdPartyWellsSurface,Shape -e np -C OBJECTID,SDE -t GEOGRAPHY -R 1047 -i sde:sqlserver:hou-sqldev10 -s hou-sqldev10 -D sde -u user -p pass      Where 1047 is a valid SRID (in the SDE_spatial_references table) for this data. It is failing with the error:               Error: The given coordinate references are incompatible (-138).     Error: Cannot Create Layer.      Looks like this used to be a bug back in 9.3 that got fixed. I checked and I do not have any null instances in the table. Googling and the ESRI sites didn't turn up much. Running a describe on a similar dataset with SRID of 1047 yields:               C:\WINDOWS\system32>sdelayer -o describe_long -l thirdpartywellssurfacetest,shap     e -i sde:sqlserver:hou-sqldev10 -s hou-sqldev10 -D sde -u user -p pass               ArcSDE 9.3.1  for SQL Server Build 1632 Thu Feb 26 12:05:37  2009     Layer    Administration Utility     -----------------------------------------------------     Layer Description ....: <None>     Database              : SDE     Table Owner ..........: SWNGIS     Table Name ...........: THIRDPARTYWELLSSURFACETEST     Spatial Column .......: SHAPE     Layer Id .............: 121155     SRID .................: 1047     Minimum Shape Id .....: 1     Offset ...............:       falsex:       -400.000000       falsey:       -400.000000     System Units .........: 1145324610.666666     Z Offset..............:          0.000000     Z Units ..............:          1.000000     Measure Offset .......: <None>     Measure Units ........: <None>     XY Cluster Tolerance .:          0.000000055879     Spatial Index ........:       parameter:    SPIDX_GRID,GRID0=0.0178906151191854,FULL       exist:        Yes       array form:   0.0178906151191854,0,0     Layer Envelope .......:       minx:      -124.52074,        miny:        24.45000       maxx:       -72.80923,        maxy:        48.99887     Entities .............: npc     Layer Type ...........: SDE     Creation Date ........: 07/28/10 16:57:27     I/O Mode .............: NORMAL     Autolocking ..........: Enabled     Precision.............: High     User Privileges ......: SELECT, UPDATE, INSERT, DELETE     Coordinate System ....: GEOGCS["GCS_North_American_1927",DATUM["D_North_American     _1927",SPHEROID["Clarke_1866",6378206.4,294.9786982]],PRIMEM["Greenwich",0.0],UN     IT["Degree",0.0174532925199433]]      And I know all of my points lie within that layer envelope: ![alt text](http://i.stack.imgur.com/37W1C.png) Any ideas on what "The given coordinate references are incompatible (-138)" could be telling me here?


title: WGS1984 polygon labels incorrectly located in web mercator projection

text: I have polygons for states/provinces (wgs1984). When I label these polygons in ArcGIS (web mercator) some of the labels occur in random locations. For example, texas might show up in southern canada, while ontario and quebec will show up in virginia. Any suggestions?


title: Problem Inserting feature with Geotools WFS

text: With this tutorial http://docs.codehaus.org/display/GEOTOOLS/Data+Writing When I try to insert a feature               try {     DataStore geoToolsClient = getGeoToolsClient();          Transaction transaction = new DefaultTransaction();     try {          FeatureWriter<SimpleFeatureType, SimpleFeature> featureWriter = geoToolsClient.getFeatureWriterAppend(typeName, transaction);                      SimpleFeature simpleFeature = featureWriter.next(); // The problem is here!                   simpleFeature.setAttributes(featureToInsert.getAttributes());                        featureWriter.write();                        transaction.commit();                      } catch (Exception ex) {                     transaction.rollback();                 throw ex;                    }        This bug occurs at featureWriter.next(); statement which is supposed to provide an empty SimpleFeature. Then, we have to fill and save it. This attempt to create an empty SimpleFeature by giving a value to NULL nillable attributes and default non nillable attributes. We have a column of SDO_GEOMETRY type is not nillable in the bind table. When the "next ()" attempts to create a default value for this attribute, the following method is called to create a default instance: org.geotools.data.DataUtilities # defaultValue (Class Type) But then the attribute is "com.vividsolutions.jts.geom.Geometry" type and defaultValue method does not know how to create a default instance of this type. This exception is thrown : throw new IllegalArgumentException(type + " is not supported by this method"); And I’m not able to insert a spatial feature in the database Thanks for any ideas! Jeremie


title: Add-in; Failure to get code to execute when toolbar button pressed

text: I had written several VBA routines for ArcGIS up through v9.3 and have since upgraded to 10. Since VBA will be a thing of the past in the next version, I downloaded and installed VB Express 2008 and successfully (I think) migrated the VBA to VB.Net 2008. I saved out each of the routines to their respective *.vb files and ran the Add-In wizard for a new project, adding each routine as a button added to a toolbar. Everything compiles and I get my toolbar to show up in ArcMap, but when I press any of the buttons, they grey out and nothing executes. Similarly, the breakpoints I have inserted go hollow with a yellow triangle and when I hover over the breakpoint, it states that "The breakpoint will not currently be hit. No symbols have been loaded for this document." I think that somehow I am missing the magic of how to tie the XML stuff (which I understand controls the toolbar and buttons) to the VB.Net code (which perform the actions). Is there any VB code that links the code back to the XML so it functions? Any insight is greatly appreciated. Thanks, Dale Bridgford P.S (added later) To provide a bit more background, I am enclosing a portion of the code. This is typical of what I have cludged together from various posts. All variations that I have tried are equally inoperative. Again, any insight as to where I have gone astray would be greatly appreciated. Thanks, Dale Bridgford ` Imports ESRI.ArcGIS.Geometry Namespace ERMDTools               Public Class btnAddHyperlink         Inherits ESRI.ArcGIS.Desktop.AddIns.Button         Public Sub New()              End Sub              Protected Overrides Sub OnClick()             Dim pFeature As IFeature, pLayer As ILayer             Dim pMap As IMap             Dim pSRE As SpatialReferenceEnvironment             Dim pQF As IQueryFilter             Dim pMSR As ISpatialReference             Dim pGeoDataset As IGeoDataset             Dim pMxDoc As ESRI.ArcGIS.ArcMapUI.IMxDocument             pMxDoc = Nothing             pMxDoc = m_App.Document             pMap = pMxDoc.FocusMap             pMSR = pMap.SpatialReference `


title: ArcGIS 9.3. How do I increase Picture Marker Symbol performance?

text: I have a layer with 400 point features that display a picture marker symbol on the map. We expect this number to increase a great deal. Unfortunately it's already noticeably slow. This is being rendered through a GIS portal written in C#, but it's also not that great when viewing through ArcMap. I'm pretty sure it's a performance thing with the PictureMarkerSymbol. Is there a way to boost it's performance? If not, are there other options for displaying custom icons for map points?


title: Error ORA-00920: invalid relational operator

text: I am running this SQL in Oracle Spatial in postgis worked. In the Oracle error `[Err] ORA-00920: invalid relational operator` I need to select all geometries that touch on my streets layer, ie the corners.               SELECT DISTINCT         (b.cdidelog) AS Cruzamento     FROM         "GIVVI_EIXO_LOGR" A,         (             SELECT                 *             FROM                 "GIVVI_EIXO_LOGR"         ) b     WHERE         ST_Touches (A.geom, b.geom)     AND A.cdlog = 7775026     AND A.nmidelog <> b.nmidelog


title: Roles for uploading layers in geoexplorer

text: What kind of users I have to create for uploading layers? I tried with users that are not admins and they have specific role with write access on workspace, but it works only once or twice. What I need is following: many "public" users with anonymous login and access to public workspace (stores and layers) and several power users which can see all in public workspace but they have their common private workspace for writing. Private workspace is for uploading their shapefiles and geotifs, and that workspace first users don't see. So, is it possible to have non-admin users for that feature in geoexplorer? OS is Windows 2008 Server, I don't know is it necessary to give some permissions on tmp\upload and private folders?


title: "qgis-bin.exe has stopped working" whenever I try to open layer properties

text: I am using QGIS Desktop 2.2.0 on Windows 8. Whenever I try to open up any layer by double-clicking on it, a window pops up saying that qgis-bin.exe has stopped working. I have tried using different shapefiles, but the problem occurs every time. I also tried re-installing QGIS multiple times, and I even tried installing an older version of QGIS (2.0.1), but none of these solutions worked. I have run out of things to try and I haven't been able to find anyone who has had a similar problem. Please let me know if you have any ideas. Thanks in advance!![enter image description here](http://i.stack.imgur.com/DR3sP.png)


title: How to interact layers with a nested loop?

text: I've dealing with the code below with such a long time. I have to show in a raster image the result of calculating the cost of the electricity for a solar plant. The thing is that I have to interact a layer "DNI" (direct normal irradiation) with the formula, and the value of the formula changes with a factor for each year (costReductionFactorPlant) during the 30 years of the life plant. The factors are shown below. When I run the code, I receive the following error:               Runtime error      Traceback (most recent call last):       File "<string>", line 56, in <module>       File "c:\program files (x86)\arcgis\desktop10.1\arcpy\arcpy\sa\Functions.py", line 4049, in Times         in_raster_or_constant2)       File "c:\program files (x86)\arcgis\desktop10.1\arcpy\arcpy\sa\Utils.py", line 47, in swapper         result = wrapper(*args, **kwargs)       File "c:\program files (x86)\arcgis\desktop10.1\arcpy\arcpy\sa\Functions.py", line 4046, in wrapper         return _wrapLocalFunctionRaster(u"Times_sa", ["Times", in_raster_or_constant1, in_raster_or_constant2])     RuntimeError: ERROR 000732: Input Raster: Dataset D:_ArcGIS\Python\LCOE_test does not exist or is not supported     >>>       I'm afraid that there is something wrong when calling the layer "DNI". But no clue. Maybe there is also an error with the loop... I've tryed to look at it also, for me looks good, but I'm not an expert. I'm running ArcGis 10.1. All suggestions are more than welcome. CODE:               import arcpy, os, sys     from arcpy import env     from arcpy.sa import *          # Set Workspace     env.workspace = "D:\02_ArcGIS\Python\LCOE_test"          # Check out the ArcGIS Spatial Analyst extension license     arcpy.CheckOutExtension("Spatial")          # Variables definition     capacity = 100000     actualOutput = .94 * capacity     storageHours = 16     loadFactor = (.57 * storageHours + 7.1)/24     costReductionFactorPlant = [0.4790, 0.4533, 0.4297, 0.4080, 0.3880, 0.3696, 0.3526, 0.3368, 0.3222, 0.3087, 0.2961, 0.2844, 0.2735, 0.2633, 0.2538, 0.2449, 0.2366, 0.2288, 0.2215, 0.2146, 0.2081, 0.2021, 0.1964, 0.1910, 0.1859, 0.1811, 0.1765, 0.1723, 0.1682, 0.1643, 0.1607]     efficiency = [0.1822, 0.1827, 0.1831, 0.1834, 0.1837, 0.1840, 0.1841, 0.1843, 0.1844, 0.1845, 0.1846, 0.1847, 0.1847, 0.1848, 0.1848, 0.1849, 0.1849, 0.1849, 0.1849, 0.1849, 0.1849, 0.1850, 0.1850, 0.1850, 0.1850, 0.1850, 0.1850, 0.1850, 0.1850, 0.1850, 0.1850]     years = 30     i = .1     ir = .01     lcoePlantTotal = 0     productionPlant = loadFactor * actualOutput * 8760     DNI =  "D:\02_ArcGIS\Python\LCOE_test" #layer#          Count=len(costReductionFactorPlant)+ 1     counter=1          while counter < Count:         for i in range(counter,(counter+1)):                      # Set the cell size environment using a keyword.                 arcpy.env.cellSize = "MAXOF"                      # Set the extent environment using a keyword                 arcpy.env.extent = "MAXOF"                      #intermediate formula plant                      # size = (capacity * 8760 * loadFactor)/(DNI * efficiency[counter-1])                      size1 = (capacity * 8760 * loadFactor)                 timesConstant = efficiency[counter-1]                 size2 = Times(3, timesConstant)                 print size2                 size = Divide(size1, size2)                      # Calculate power block & storage                               powerBlock = capacity * 1484.9918 * costReductionFactorPlant[counter-1]                 storage = storageHours * capacity * 39.45759 * costReductionFactorPlant[counter-1]                      # Calculate mirrowField                  # mirrorField = size * 155.6975 * costReductionFactorPlant[cont-1]                 timesConstant2 = costReductionFactorPlant[counter-1]                 inconstant = Times(155.6975, timesConstant2)                 mirrorField = Times(size, inconstant)                      # calculate PowerTower                 powerTower = ((capacity * 121.03883) + (capacity * 128.50378)) * costReductionFactorPlant[counter-1]                      # Calculate investmentPlant                 investmentPlant = powerBlock + storage + mirrorField + powerTower                      # Calculate omPlant                 omPlant = Times(investmentPlant, 0.048)                      #formula: lcoePlant = ((investmentPlant * (((i * (1 + i)^counter)/((1 + i)^counter-1))+ ir))+omPlant)/productionPlant*100                 # Formula Part 1 = (investmentPlant * (((i * (1 + i)^counter) Power symbol = **                 InConstant3= i * ((1 + i)**counter)                 FormulaPart1 = Times(investmentPlant, InConstant3)                      # Formula Part 2 = ((1 + i)^counter-1))+ ir))+omPlant)/productionPlant*100                 FormulaPart2 =(((1 + i)**(counter-1)+ ir)+omPlant)/productionPlant*100                      # Formula                 lcoePlant = Divide(FormulaPart1, FormulaPart2)                      lcoePlantTotal = lcoePlantTotal + lcoePlant                      lcoePlantTotal.save = "D:\\02_ArcGIS\\Python" + str(lcoePlantTotal) + ".img"          print "lcoePlanttotal calculated"


title: Changin features' pointRadius based on the current zoom level

text: I need to have my features change their `pointRadius` based on the current `Zoom` level. My questions are:   1. Would it be possible to use `Rule Based Styling` techniques to overcome the problem?   2. If yes, it will work on `map zoom changed` as well?   3. If no, what are the alternative solutions and what you offer instead? (except redrawing features)


title: What does 'DBMS table not found' message mean?

text: I have a plain postgressql table (version 8.2). I am trying to connect to it to it from ArcGIS desktop 10.0, and use as a query layer, as given in this document: http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#/Connecting_to_a_database/00s500000037000000/ I can see the tables, but I get the error message 'DBMS table not found' ![enter image description here](http://i.stack.imgur.com/CUFQd.png) What exactly is the issue? How do I connect to it?


title: How do I convert a polygon into a point?

text: I want to convert polygons into points. I know how it's done in ArcGis and I think this may sound trivial but it's a serious issue :/. I'll appreciate any hints :)


I ST_Dwithin the function returns TRUE if the geometry A is completely inside B. I happen to have many polygons (representing regions) and need to know which streets are within each of these polygons. But it has some streets that hath part in a region and part in another. Which function is most appropriate to select postgis geometries that have any part inside another?


title: qgis web client installation

text: I recently joined github and I found this very interesting application. I was trying to test it and install the qgis/qgis-web-client application. I got the following error by following the installation guide: > filippo@filippo-Aspire-5742G:/etc/apache2/sites-available$ sudo a2ensite > qgis-web-client.conf ERROR: Site qgis-web-client does not exist! what's the problem of that? I went through the full installation then i tried to see the index.html and I got at the following link the right web page index.html but the link to the map "hello world" doesn't seem to work. Also when I tried GetCapabilities but I got a server not found answer. I guess I got a problem with the cgi-bin folder and its content. Then in 1. Check the symlink in cgi-bin is correct what I suppose to have been linked in here? And finally in 2. Check the QGIS libs are in your /etc/ld.so.conf path, how do I check this step? Many thanks in advance for the help


title: Do you know how to become Navteq (Here Maps) reseller?

text: Before we begin, I'd like to state that I know that Navteq has been acquired by Nokia and Nokia has been acquired by Microsoft. However, Navteq was out of the deal. When I go through pages of Here Maps (old Navteq), I cannot find any information regarding the distribution of raw data other than APIs. However, on the other hand, there are plenty of companies that claims to be Navteq data provider and they are able to provide sample datasets and license those datasets for commercial usage. As Nokia removed all pages related to Navteq and redirected the visitors to here.com I wonder if someone knows how to become Here Maps Premium Reseller?


title: How to convert AutoCAD file to Shapefile for map from ArcGIS API 3.2 for JavaScript?

text: I have an AutoCAD file that I want to convert to Shapfile to appear in a map by API for JavaScript in .net. What do I do?


title: Suggestions for vector format for polygons

text: I am using two polygon tables in postgis (A and B ). I need to make a vector file that incorporates information of the overlapped polygons B intersecting each polygon on (A). I can make a single shapefile for each polygon of A with the resulting intersections and other info of B. Is there a way in which I could "merge" all these resulting polygons (each with its own attributes) into a single file ? I was thinking that maybe using an array column but I also need other info related to the intersection like area and percentage of coverage.


title: How to upgrade from QGIS 1.4 to 1.8 on Linux Mint 12?

text: I've been trying to upgrade from the Enceladus 1.4 version of QGIS to 1.8 on a Thinkpad T42 running on Mint 12. When I try to upgrade from within QGIS help menu I receive the following error message: "Unable to communicate with QGIS Version server Unknown network socket error" From the official web site, if I use the command line to install QGIS using               sudo apt-get update     sudo apt-get install qgis     sudo apt-get install python-qgis     sudo apt-get install qgis-plugin-grass      I still get the 1.4 version. Any thoughts? Thanks.


title: How to quickly match many high-scoring geocoding ties

text: I am geocoding Portland Metropolitan Area in ArcMap 10.2 based on Street Address, City, State and Zip. I get decent matches (about 93%). However, I get thousands of high-scoring ties. For example 2 or 3 matches scoring 100%. These tied values are all "separate" points (e.g. two different points within a property polygon, and on point on the street facing the polygon). Is there a way to either:   1. non-manually match at random based on any one of the top scoring ties, or   2. take the average lat and average long of the top scoring ties? I am a noob, so I will happily entertain ways to improve this question based on constructive criticism.


title: What methods exist for geohash of polygon?

text: GeoHash is an excellent method to code the latitude and longitude of a point into a string. It's concise, due to the use of numbers + letters, and the precision can be tuned adding characters at the end. I'm looking for an equivalent method to code polygons, not points. I saw two ways to do it:   * geohash all polygon points, then concatenate them, and then compress the resulting string using classical text compression algorithms.   * reference all geohashes included into the polygon. It's the RecursivePrefixTree method - see http://www.opensourceconnections.com/2014/04/11/indexing-polygons-in-lucene-with-accuracy/. I would like the code could be manipulated by humans, ie. manually copied, not only by machines as an index. Edit: see also the discussion (in French) at http://seenthis.net/messages/269838.


title: Leaflet highlight feature by ID

text: I have a map application where:   1. I load polygons as GeoJSON objects    2. I draw all polygons in a single layer   3. I highlight a polygon on mouse hover What I need to do is the following: on mouse hover, I need to highlight multiple polygons (I get the list of polygons somehow - irrelevant -). I have the IDs of all the polygons I need to highlight, but I can't see how can I access them though the API.


title: JTX/Workflow Manager dependencies to control workflow

text: When creating a new job in JTX, one of the required fields (under "Additional Information") tab we have is "Asset Type". The three choices for "Asset Type" are "Water", "Sewer" and "Water and Sewer". My question is, is it possible to setup JTX where depending on the Asset Type, certain step(s) in the workflow will be skipped. For example, if I choose Water as my asset type, I may want to skip a certain step and if I choose sewer as the asset type I may want to skip a different step.


title: how do i assign values to shapes to be displayed in my map

text: I apologize for the simplistic nature of this question, I'm VERY new to QGIS. I'm trying to use QGIS to display political information. I have a map of my state, with each town as it's own shape. I have an ID for each town/shape. What I want to do is assign a value to each town/shape, and have a scaled display of these values by town/shape. Can anyone help me with the process for assigning values to each town/shape? Thanks!!!


title: How to computational efficiently check if there are any overlapping polygons in a group of polygons?

text: I need to do some processing with groups of polygons. If any of the polygons in the group overlaps other polygon the result is unpredictable and don't throw any error. So, I need to detect if there is any overlapping polygon inside the group and which polygons are overlapping. The groups are made of about 15 polygons, so checking each polygon against the others would be costly. I wonder if there is an efficient way to do that. The code is in Python and can use either the OGR bindings or Shapely. A theoretical answer is also appreciated.


title: Displaying adjacent polygon lines beside each other instead of overlapping

text: I would like to display adjacent polygons with the line work beside each other, not one on top of the other for that shared boundary. I do not want to edit the polygon boundaries - only want to change the way they are displayed. Is that possible?   ![image shows the blue polygon boundary is under the green polygon. Instead of it looking like this, I would like to be able to see both polygon boundaries beside each other where the boundaries are shared](http://i.stack.imgur.com/A3XIu.jpg)


title: Recent Documents list in ArcMap is not being updated

text: On a colleague's machine, the list of Recent Documents isn't being populated when closing an MXD.   * she opens map1.mxd   * saves the document and closes ArcMap   * re-opens ArcMap map1.mxd isn't shown in the _> Existing Maps > Recent_ section on the welcome screen, or under the _> File_ menu. Where is the list of recent documents stored? (Interestingly, her welcome screen does list 2 documents, just not the most recently-accessed ones.) Of course this is all within a corporate system with network drives, locked- down permissions, etc.... ArcMap 10.2.1 on Windows 7.


title: Are there any Javascript map labeling engines available?

text: Does anyone know of a javascript library containing functions that would allow me to generate labels for points on the map so that the labels don't overlap one another, and have optional leader lines connecting the label to the point? I'm using Esri's javascript API, but would also be interested in finding a generic conflict detection/resolution algorithm that I could adapt.


title: Can an entry in Esri GeoPortal Server contain both ArcGIS Server and *.zip links?

text: Does anyone know whether a dataset in Esri GeoPortal Server can support both Preview (in ArcGIS Server), and Download (of a zip file) on the same dataset? Eg this dataset has a link to Preview in ArcGIS Server. Could this same dataset also contain a Download option, which linked to a zip file? If so, how would that be set up? ![enter image description here](http://i.stack.imgur.com/nYNVH.jpg) ( _Open_ takes you to the REST endpoint. _Preview, KML, NMF_ and _LYR_ are auto-generated by ArcGIS Server. Can I also include a Zip file of related data, which I would create manually?)


title: Polygon Between Contour Line

text: I have a contour generated from SRTM,   ![enter image description here](http://i.stack.imgur.com/CpUeX.png) How to fill the space between that contour with polygon (for all part) but it generates field in attribute table automaticaly, example in image above between contour 6 and 7 then the polygon generated should have custom field "nilai" with value "6-7" in attribute table. I've tried polyline to polygon but it's not working I should enter the value of polygon manually, there are so many to be edited one by one. **1st Method Using Spatial Analyst License (you can zoom in youre browser using "ctrl and +")**   1. Created Polygon (per 1m for 1m - 30m using feature to polygon tool) from contour created from dtm using **contour list** tool. http://i.stack.imgur.com/OdBXL.png ![enter image description here](http://i.stack.imgur.com/OdBXL.png)   2. This is my DTM raster http://i.stack.imgur.com/mez5C.png ![enter image description here](http://i.stack.imgur.com/mez5C.png)   3. Creating Zonal using **Zonal Statistic as Table** tool http://i.stack.imgur.com/ieJlz.png ![enter image description here](http://i.stack.imgur.com/ieJlz.png)   4. Zonal result attribute table http://i.stack.imgur.com/8w8iG.png ![enter image description here](http://i.stack.imgur.com/8w8iG.png)   5. Join polygon with zonal result http://i.stack.imgur.com/LDhA6.png ![enter image description here](http://i.stack.imgur.com/LDhA6.png)   6. Attribute Table of joined result http://i.stack.imgur.com/b8Ve5.png ![enter image description here](http://i.stack.imgur.com/b8Ve5.png)   7. Checking the result after i **round** it, from that result maximum and minimum value are not match with the contour line http://i.stack.imgur.com/94Gdb.png ![enter image description here](http://i.stack.imgur.com/94Gdb.png) **2nd Method Using Spatial Analyst License** I use this script provide by Michael                   >>> import sys, os, arcpy          InPolygons = sys.argv[1]     NearContour = sys.argv[2]     OutPolygons = sys.argv[3]          ElevationFieldName = "contour" # change this to match your field name          InFolder = os.path.dirname(InPolygons) # the folder the input is in     Working = InFolder + "\\Working.gdb"     WorkingContour = Working+"\\Contour"     WorkingPoly = Working+"\\Poly"          arcpy.CreateFileGDB_management(InFolder,"Working.gdb") # create a working space     arcpy.AddMessage("Copy 1")     arcpy.CopyFeatures_management(NearContour,WorkingContour) # need to do this to take charge of names     arcpy.AddMessage("Copy 2")     arcpy.CopyFeatures_management(InPolygons,WorkingPoly)           # remove 'bad' edges, tagged with elevation = -9999 to show they are a boundary     arcpy.MakeFeatureLayer_management(WorkingContour,"Remove",ElevationFieldName + " = -9999")     arcpy.DeleteFeatures_management("Remove")     arcpy.Delete_management("Remove") # relax, it's just removing the layer, just to keep it clean          # Create a near table     arcpy.AddMessage("Near table")     arcpy.GenerateNearTable_analysis(WorkingPoly,WorkingContour,Working+"\\NearTable","0",closest="ALL")     arcpy.AddField_management(Working+"\\NearTable","Elev","DOUBLE")          try:         arcpy.AddIndex_management(Working+"\\NearTable","NEAR_FID","NEAR_FID")     except:         arcpy.AddWarning("Unable to index")          # join and copy elevations     arcpy.AddMessage("Join & calc")     arcpy.MakeTableView_management(Working+"\\NearTable","View")     arcpy.AddJoin_management("View","NEAR_FID",WorkingContour,"OBJECTID")     arcpy.CalculateField_management("View","NearTable.Elev","!Contour."+ElevationFieldName+"!","PYTHON")     arcpy.Delete_management("View") # relax, it's just removing the table view, just to keep it clean          # remove join and summary statistics     arcpy.AddMessage("Make stats")     #arcpy.RemoveJoin_management("View")     arcpy.Statistics_analysis(Working+"\\NearTable",Working+"\\NearStats",[["Elev","MIN"],["Elev","Max"]],"IN_FID")          # copy the min and max to the working contours     arcpy.AddMessage("Copy min/max")     arcpy.AddField_management(WorkingPoly,"MinElev","DOUBLE")     arcpy.AddField_management(WorkingPoly,"MaxElev","DOUBLE")     arcpy.MakeFeatureLayer_management(WorkingPoly,"Layer")     arcpy.AddJoin_management("Layer","OBJECTID",Working+"\\NearStats","IN_FID")     arcpy.CalculateField_management("Layer","Poly.MinElev","!NearStats.Min_Elev!","PYTHON")     arcpy.CalculateField_management("Layer","Poly.MaxElev","!NearStats.Max_Elev!","PYTHON")     arcpy.RemoveJoin_management("Layer")     arcpy.Delete_management("Layer") # relax, it's just removing the layer, just to keep it clean          # check for field     arcpy.AddMessage("Checking for field")     Fields = arcpy.ListFields(WorkingPoly,"nilai")     if not Fields:         arcpy.AddMessage("-Adding field")         arcpy.AddField_management(WorkingPoly,"nilai","TEXT",field_length="10")          # format the field with the min and max,     # this could be done with a field calculation     # but in this case I want a special case for isolations     arcpy.AddMessage("Cursoring")     with arcpy.da.UpdateCursor(WorkingPoly,("nilai","MinElev","MaxElev")) as UpCur:         for Feat in UpCur:             if (Feat[1]==Feat[2]):                 # values are the same = an isolation                 Feat[0] = str(int(Feat[1]))             else:                 Feat[0] = "%d - %d" % (Feat[1],Feat[2])             UpCur.updateRow(Feat)          arcpy.AddMessage("Copy out")     arcpy.CopyFeatures_management(WorkingPoly,OutPolygons)       _Processing :_ ![enter image description here](http://i.stack.imgur.com/TsKNi.png) _Result_ http://i.imgur.com/faaESEb.png ![enter image description here](http://i.stack.imgur.com/1SA9D.png) My Polygon and attribute table are shown in the printscreen. Please help, any sugestion using python or else. Thanks.


title: Export long lat from polyline as csv file in QGIS

text: I want to draw a line in QGIS and then, export the long-lat of each point into a Comma Separated Value (csv) file. I tried doing this by creating a new shapefile layer and then using the `Add Feature` tool to draw my line. I then tried to export my line via `Save A`s on the layer and exporting as CSV. Unfortunately, the only thing that gets saved to the CSV file is the id of the layer. I've been unable to figure out how to export the coordinates. Some help please.


title: How to filter date values using OGR 1.11.0?

text: With OGR 1.10.1 this works:               ogrinfo test2.shp -sql 'SELECT COUNT(*) FROM test2 WHERE DATE < "1970-01-01"'      Note that DATE is a column of type date (not string). (Here's test data if you want to try: http://ge.tt/6DBclJh/v/0) The same query in 1.11.0 causes               OGR[3] error 1: Type mismatch or improper type of arguments to < operator.      Has the syntax changed or is this really a bug?


title: How to make distance calculated using Latitude and Longitude more accurate

text: I'm coding a Web App whereby users can choose their city (start) and a destination city, and find out the approximate distance between both cities, using Latitude and Longitude of both cities to calcultate. Then i thought that since this will be done using the concept of a Straight Line on a sphere, it means that i'll be getting a result of a straight line which is the shortest, but in reality (driving route especially), will have hills, curves, turns. It might be practically impossible to have an accurate distance, but how can i make the app more accurate that just a straight line distance. I'll be most delighted by you opinions, suggestions, answers and comments. Thanks in advance.


title: Can Identify tool to return the raw values of rasters?

text: I have a raster in ArcGIS 10 SP4 set to bilinear resampling in its display properties. When using the identify tool it returns the resampled value instead of the raw value. Can we get the identify tool to natively return the raw value?


title: Geowebcache re-projection

text: I have one layer (raster geotiff image) published in geoserver with EPSG: 3011 projection. But I want this layer in EPSG 4326 projection for use it in google API web application. I have searched on web and get to know that Geowebcache can do this conversion by adding custom projection. I am not sure, if geowebcache will change the layer projection from EPSG: 3011 to EPSG: 4326. Is it possible in geowebcache? If yes, then, should i need to add custom projection with EPSG: 4326? and then assign that layer (layer with EPSG: 3011) to custom projection. need brief answer as i am new to geowebcache.


title: algorithm for displaying holes in a polygon by translating it into several hole-less polygons?

text: I'd like to plot polygons with holes on a framework that doesn't support holes. The approach to me seems to be to translate the single polygon into several ones adding up to the whole polygon, leaving out the holes. Is there a standard algorithm for that? Maybe even with code? Or is there another standard way to solve the problem of plotting polygons with holes?


title: How do you convert the radius decimals of a circle in Openlayers to Kilometers

text: > **Possible Duplicate:**   >  Unproject Radius with Openlayers with CRS:84 Projection I allow a user to draw a circle in Openlayers. I created an alert to show the radius of the circle. I now need to convert the radius to kilometers. For example 0.001788057579325097 (How many kilometers is this) How do I change this to kilometers? I think its just the decimals that is the problem.                     // The Drag function is now made and called              // Add Drag          drag = new OpenLayers.Control.DragFeature(polygonLayer, {          autoActivate: true,          onComplete: displayRadius          });                function displayRadius(circle) {          var area = circle.geometry.getArea();          var radius = 0.565352 * Math.sqrt(area);          alert(radius);           }


title: ARCPAD - VBscript to update values in attribute table

text: I am working on an automating tool for ArcPad for inventory issues on our mobile GIS machine (Trimble GeoXH - Geoexplorer 2008 series - Windows mobile 6.1 and Arcdap 10.0.3). I have developed the tool in VBscript. Short description what the tool must do: We have a point shapefile in Arcpad. Everytime when we manually add a point, it automaticaly counts the new ID by clicking on a button in the toolbar. The new ID is the last ID + 1. I succeeded already in looping through the attributes, finding the last ID and calculating a new ID. However, I can't add the new calculated ID for the added point in the attribute table. It's like I can't open the recordset in edit mode. Below you find my source code.   Do you know what I'm doing wrong? Thanks for helping! Kris Sourcecode Sub Test               Dim objDBF, strLaatstePaalID, NewPaalNR, NewPaalID, Nummer, PaalCode           Numbre = 0          Set objLyr = Map.Layers("Palen")     Set objDBF = objLyr.Records     'Look for record with the highest ID     'ID is text format like 'P_00450'     objDBF.MoveFirst     Do Until objDBF.EOF         PaalCode = objDBF.Fields("PAAL_ID").Value         If PaalCode <> "" Then 'the new added point (which is the last) has no value in the ID             If Int(Right(PaalCode,(Len(PaalCode)-2)))>Numbre Then                 'Make an integer of the text code                 Numbre = Int(Right(PaalCode,(Len(PaalCode)-2)))             End If         End If         objDBF.MoveNext     Loop          NewPaalNR = Numbre + 1          'Maak een nieuw text ID with the new numbre     If NewPaalNR < 10 Then         NewPaalID = "P_0000" & CStr(NewPaalNR)          ElseIf NewPaalNR < 100 Then         NewPaalID = "P_000" & CStr(NewPaalNR)          ElseIf NewPaalNR < 1000 Then         NewPaalID = "P_00" & CStr(NewPaalNR)          ElseIf NewPaalNR < 10000 Then         NewPaalID = "P_0" & CStr(NewPaalNR)          Else         NewPaalID = "P_" & CStr(NewPaalNR)          End If          MsgBox ("New Id: " & NewPaalID)          'go to the new added point     objDBF.MoveLast          if objDBF.Fields("PAAL_ID").Value ="" Then         'add the new value         objDBF.Fields("PAAL_ID").Value = NewPaalID         objDBF.Update         MsgBox("Paal_ID " & NewPaalID & " added in Palen.shp")     Else         MsgBox("No empty feature added")     End If            Set objDBF = Nothing      End Sub


title: Clipping a raster in R

text: I am building a map for the northeastern U.S.A. The map background needs to be either an altitude map or a mean annual temperature map. I have two rasters from Worldclim.org which give me these variables but I need to clip them to the extent of the states I am interested in. Any suggestions on how to do this. This is what I have so far:               #load libraries     library (sp)     library (rgdal)     library (raster)     library (maps)     library (mapproj)               #load data     state<- data (stateMapEnv)     elevation<-raster("alt.bil")     meantemp<-raster ("bio_1.asc")          #build the raw map     nestates<- c("maine", "vermont", "massachusetts", "new hampshire" ,"connecticut",       "rhode island","new york","pennsylvania", "new jersey",       "maryland", "delaware", "virginia", "west virginia")          map(database="state", regions = nestates, interior=T,  lwd=2)     map.axes()          #add site localities     sites<-read.csv("sites.csv", header=T)     lat<-sites$Latitude     lon<-sites$Longitude          map(database="state", regions = nestates, interior=T, lwd=2)     points (x=lon, y=lat, pch=17, cex=1.5, col="black")     map.axes()     library(maps)                                                                  #Add axes to  main map     map.scale(x=-73,y=38, relwidth=0.15, metric=T,  ratio=F)          #create an inset map           # Next, we create a new graphics space in the lower-right hand corner.  The numbers are proportional distances within the graphics window (xmin,xmax,ymin,ymax) on a scale of 0 to 1.       # "plt" is the key parameter to adjust         par(plt = c(0.1, 0.53, 0.57, 0.90), new = TRUE)            # I think this is the key command from http://www.stat.auckland.ac.nz/~paul/RGraphics/examples-map.R         plot.window(xlim=c(-127, -66),ylim=c(23,53))            # fill the box with white         polygon(c(0,360,360,0),c(0,0,90,90),col="white")            # draw the map         map(database="state", interior=T, add=TRUE, fill=FALSE)         map(database="state", regions=nestates, interior=TRUE, add=TRUE, fill=TRUE, col="grey")      The elevation and meantemp objects are the ones that need to be clipped to the area extent of the nestates object. Any input would help


title: How to create a circle vector layer with 12 sectors with Python and to save average value of underlying raster in each sector?

text: Hi I am writing a code with Python 2.7 using Arcpy 10.1 and I guess I am doing stuck with a silly problem that lies under my eyes. I have created a grid of points (with "Create fishnet" function) over a region where I have also raster data (for instance land cover or DEM). For each point I need to associate a circle in vector format (with center in each point) divided in "n" sectors of 1 km radius for example. The goal is to calculate the mean value of all pixels of the land cover or of the DEM situated in each sector and to save them in the attribute table of the circles. ![points grid](http://i.stack.imgur.com/VSpiB.jpg) I found this very helpful post (How to create a circle vector layer with 12 sectors with Python/pyqgis?) but I am somehow stuck where I have to save the data in shapefile instead of geojson and to go further with the estimate of the mean values for each sector. Here is part of the code take from the previous code:               # helper function to calculate point from relative polar coordinates (degrees)     def polar_point(origin_point, angle,  distance):         return [origin_point.x + math.sin(math.radians(angle)) * distance, origin_point.y +math.cos(math.radians(angle)) * distance]          features = []          for x in xrange(0,int(sectors)):         segment_vertices = []              # first the center and first point         segment_vertices.append(polar_point(center, x,y))         segment_vertices.append(polar_point(center, start + x*sector_width,radius))              # then the sector outline points         for z in xrange(1, steps_per_sector):             segment_vertices.append((polar_point(center, start + x * sector_width + z * step_angle_width,radius)))              # then again the center point to finish the polygon         segment_vertices.append(polar_point(center, start + x *  sector_width+sector_width,radius))         segment_vertices.append(polar_point(center, x,y))              # create feature         features.append(geojson.Feature(geometry=Polygon(segment_vertices)))          # prepare geojson feature collection     res = geojson.FeatureCollection(features = features)          # write to file     f = open('sector.shp', 'w')     f.write(geojson.dumps(res))     f.close() }      Any piece of advice out there please? Is any other approach maybe more practical/feasible? I was thinking also about "Euclidean direction" and "zonal Statistics", but it deals with many raster data and the calculation might be very long in case of large regions to analyze. Best regards.


title: How to handle adding a PostGIS layer in pyQgis without saving password?

text: I'm writing python plugin and trying to find the right way to add PostGIS layer on the map without storing user password. After adding layer on the map, I strip out the password from uri like this: `vlayer.dataProvider().setDataSourceUri(uri.removePassword(uri.uri()))` However, when user saves the project, full connection string stays including password string. I know that QGIS can popup user name/password dialog for me if I don't put these data on uri, but then I need to get somehow entered data, since I need to know user name and password anyway.


title: Plotting time-series data in ArcGIS

text: I'm attempting to produce interpolated plots of water quality time-series data from ArcSDE. This is the data format:               OID    WELL_ID  WELL_NAME  SAMPLE_ID  SAMPLE_DAT    PARAM   RESULT          90-------1---------WELL A-----19-24231----5/4/1966----CHLORIDE--------130          91-------1---------Well A-----19-24295------7/7/1966----CHLORIDE--------130          92-------2---------Well B-----19-24478----1/6/1967-----CHLORIDE---------122          93-------3---------Well C-----19-24597----5/5/1967-----CHLORIDE----------115      and onward for thousands of records. Each well has many records, and the sampling is erratic, some not being sampled for years. I've tried a few things, but I want to see if anyone else has encountered something similar and how they approached it.


title: Need help tuning gdalwarp to match source cell size and alignment

text: I was hoping gdalwarp could merge and crop MODIS tiff images such that the original cell size and alignment is maintained. I tried the command below, but cell size is off slightly and cell alignment is offset from 10 to 100 meters.               gdalwarp -cutline basins.shp -crop_to_cutline input1.tif input2.tif input3.tif outputmergeclip.tif      I tried setting `-tr` to the source cell size, `-te` to the extent of the shapefile snapped to the source grid cell alignment and `-tap`, but these didn't make any change in the result. Am I interpreting the options incorrectly? Am I on the wrong track with gdalwarp? I'm scripting this in Python and relatively new to both python(2.7) and gdal(1.9.2).


title: raster tile - find the correct tile for a given coordinate

text: I have a large raster map which I loaded into postgis with 1000x1000 meter tiles... the corresponding table has about 200 rows, each row representing a tile. I want to get the raster value for a specific coordinate, but I don't know how to do that most efficiently. I did something like               select ST_Value(rast, ST_Transform(ST_GeomFromText('POINT(lat lon)',4326),mysrid)) from rasterdata;      but it gives me a lot of               Attempting to get pixel value with out of range raster coordinates:...      messages because finally only one tile contains the coordinate that I'm searching for... Any help is highly appreciated...


title: Alias's in the Attribute Table's Column Preview/Column Filters

text: We've spent some time setting up alias's in Layers/Properties.../Fields (Fig. 1) for database fields displayed as the column headers in the Attribute Table (Fig. 2). Switching the layout in the Attribute Table from table view to form view and selecting between the Column Preview, the column names are now no longer aliased (Fig. 3). (This is also the case for the Column Filters, but is not shown here) Is there any way to get the filters to use the aliased names? If there is a PyQGIS solution I'd be happy to hear about it. ![Alias's](http://i.stack.imgur.com/SQhy7.png) ![AttribTable](http://i.stack.imgur.com/WRLJK.png) ![FormView](http://i.stack.imgur.com/NmNrg.png)


title: Using a tree menu beside openlayers

text: Is there an aplication to add tree menu for the map showed with openlayers? I mean show a map on my site exactly how is showing in google earth for example, my layers organized in folder/subfolder, with checkboxes.


title: Openstreet map loosing information

text: I succesfully created a shapefile from openstreetmap with the OpenStreetMap tool in QGIS 2.2. When I try to create a polyline layer, it says that there are 8000 values in the db file for the property highway, but after the import I have only 6000 in the attribute list. The problem is that I don't have all streets availaible for a given city and because of that I am unable to draw parcels where the streets are their borders. Is that problem normal or am I missing something? Please help!!!


title: ArcMap crashing in geoprocessing script with AreaSolarRadiation. Does a log file exist?

text: I am writing a greoprocessing script in python and running it in arcmap 10. The script works as expected until it gets to an area solar radiation calculation               AreaSolarRadiation(input_dem, "", "400", TimeMultipleDays(2011,1,31),7)      It seems every other time I try to run the geoprocessing tool in arcmap the application crashes when it gets to that solar radiation calculation. Sometimes the crash is followed by a prompt to send ESRI an error report. Sometimes I don't even get that. The dem I am calculating area solar for is pretty small so I don't believe the calculation is too intensive. The script never crashed arcmap before I started trying to do the solar calc. I am thinking if I could see a log file I might have a better idea of whats going on. Does anyone know if the arcmap application keeps a log file I could view that would shed light on why it crashed? This is running on windows server 2003.


title: Getting ISpatialReference from ArcObject mapcontrol in VBA

text: I am working with ArcObjects and ESRI for the first time, and I am trying to do what is probably a simple task. I have a map control in VBA that is associated with a mxd. I have a table with XY data, and XY data is non-spatial, but the XY is in the same coordinate system as the mxd. I can get the XY data into a ITable, but I cannot seem to get the ISpatialReference from the map control. I would like to create a symbolized layer from the non-spatial coordinates, and add it to the map control. Any samples in VBA of how to do this would be very appreciated. Thanks again for you time and patience.


title: Exporting a set of lat/long coordinates to SVG

text: Suppose I have a set of coordinates in latitude and longitude. What would be the easiest way to export those to an SVG file? I'd rather do this programatically (any language or library will do) rather than manually, but can revert to that if necessary.


title: Display full directory in ArcMap (10) recent file list in File Menu?

text: thought this would be a pretty basic question, but I haven't found anything here or on the ESRI forums on the topic. Is there any way to change the number of characters displayed for the file pathways listed in the File Menu list of recently used mxd files? The obvious solution is to create shorter file names, but sometimes that's not as easy you might think! Hover-over functionality for that list would be an ideal, long-term solution that I'd like to see, but that's neither here nor there... Thanks!


title: SELECT TOP query in ArcGIS

text: Is there a way to use a SELECT TOP (n) within ArcGIS desktop? Is appears that since the ArcGIS query builder is effectively building the SQL statement after the WHERE clause, that this may not be possible. Here is a typical SQL example of usage: > SELECT TOP 3 * FROM Layer


title: ArcGIS - Geoprocessing creating temporary files (e.g. "g_g193"). What produces these and how can they be cleaned up?

text: I have a couple of ArcGIS extensions that are running a lot of geoprocessing models with a mix of ArcObjects. I've noticed that the some, as of yet unidentified, geoprocessing steps are creating temporary raster files (e.g. g_g193, g_g983, g_g892). The funny thing is, these files are not being directed to my selected temporary workspace or the default Arc workspace. They are actually being sent to the 'bin' folder of my extension, which is obviously an issue since these temp files should not be stored with my included DLLs and EXE files. It's important to mention that my extension makes many calls to Arc geoprocessing, so it's difficult to location the source of these files. The interesting thing is, if I walk through all of my code using break points in .NET these files don't always get created. It's like Arc is not cleaning up these temp files if the processing is moving very quickly, but it will clean them up if given time (e.g. user steps through slowly). Regardless, I don't want these files anywhere near my 'bin' folder. Any insight would be greatly appreciated!


title: How to trigger a hover event for map graphics in ArcGIS Javascript API?

text: I want to trigger an event when the user hovers over a geometry that they drew. I know how to set up hover events for DOM objects, but not for graphic elements within the map window.


title: Export Raster Data to TIFF, BMP, or PNG File Using Python

text: I am writing a Python script to export a Raster Dataset to a TIFF, BMP, or PNG file, as I would manually by right-clicking on the layer, choosing Data | Export Data, and then make my choices in the Export Raster Data window. Can this be duplicated with ArcPy, and if so, how?


title: ArcGIS API Interpolation Functions

text: Does ArcGIS API provide interpolation functions, like if I want to interpolate from my own set of values? Edit:   I know that ArcGIS provides functions for interpolation like IDW, Spline and Kriging etc, but can I use those functions to interpolate against my own dataset?


title: What is the effective and legal way to get map tiles for the size of 2 or 3 states?

text: I would like to get map tiles possibly from Open Street Map (yes I am aware about tile usage policy and would not like to breach that) or other source for the purpose to serve it as offline maps to applications. Total area will be of the size of let’s say 2 or 3 states of US (300K square miles) and from **zoom level 1 to 12 levels** (possibly in future up to 17 levels), so I think it will be couple of GBs. Possibly it can be 2 ways, **generating those map tiles locally** by OSM data or **download tiles from some source with some utility**. **Or may be I have no idea** and there is something simple exist! What are the specific ways this can be achived quicker, easier and legal way?


title: How to aggregate data which meet certain conditions using QGIS?

text: I have a point shapefile with meteorological stations which contains three columns, which are: Date: The date on which the measurement has take place (it covers a range of ten years) Id: The meteorological's station unique id. Rainfall: The amount of rainfall observed at a specific station on a certain date. Rainfall on column date in not daily but decadal (rainfall is registered every ten days during the month) so for a month I should have three recordings in total for each year, in order to get the monthly total rainfall which is what I want to do. In column date the format is the following: **2012011** meaning that the rainfall for **2012** - **january** and the **first ten days** is accumulated. The problem is that not for all stations exist measurements for all three decads of the month. So, when I I get the total rainfall for a month, I want only the years that have complete decads to be taken into account. How could I not include into the sum the years on which there are decads missing from a certain month? Is there a select command that I could use or anything similar? I use QGIS 2.4


title: Is there a gallery of QGIS maps anywhere?

text: Specifically, I want to see the gorgeous map of Chugach State Park that you see when starting QGIS 2.4 "Chugiak". Is there a website out there where people share maps? I would like to share and see others' styles and ideas. Thanks!


title: From pdf to google maps overlay

text: First post, so be gentle. Context - I'm doing some research into healthcare provision and I'd like to visualise the data. I'd like to be able to put these PCT boundaries ![enter image description here](http://i.stack.imgur.com/3vohg.jpg) as an overlay to google maps and change their colours and so on to show data... I'm running osx and I'm fairly IT friendly, but I've got no idea where to start on this and google is sending me in all manner of different diagrams with language I'm struggling to get my head around. So I'm looking for any information people can give me on:   1. Recommendations of tools that run on OSX that might help   2. A general outline of the workflow I'm going to be going though   3. Places where I could put the finished files in an open-sourcy way so that other people might be able to make use of them?  Thanking you all kindly in advance... EDIT - should point out that I'm happy to trace the lines with a tablet if that's what it takes - doesn't have to be postcode perfect just as long as it's recognizable...


title: Polymaps Project + Geojson + D3

text: How are you doing ? I am Mehul Nahar working on a project which is similar to http://www.tableausoftware.com/learning/examples/gpa-analysis. I have used Polymaps platform, where I have used cloudmade API for the base map. My problem - I am unable to implement the Geojson file, in fact I am still unable to plot Lat&long. I was wondering if you guys could help me on this. If you can give me a right direction to work on and also if you have some sample code which you could share with will be real helpful. Look forward to hearing from you.


title: How to convert from UTM Zone 36N Mercator to WGS 84 degrees of lat and long?

text: I have acquired many coordinates from a points I extracted from polygons. The points are in WGS 84 UTM Zone 36N, and I want to use them on Google Maps. I found out that they use different systems (I am no GIS expert :S) "that uses a Mercator Projection appropriate for that zone, with Units of Meters", and I want to convert them to units of degrees of Latitude and Longitude. Example of points I have: 3390003.50135 3390003.50135


title: Dynamic Attribute Table

text: This seems silly to me, but I haven't hit this before. I have multiple (30+) point tables, 100 data points each, and want to start designing symbology and doing things like point distance. However, my data source is not finalized. (I have been asked to start now, but will receive updated coordinates in a week). Is there a table format (access table, gdb, etc) where I can calculate pt distance (multiple times per table), design a multi-part symbology based on attribute, throught the table (not making it a shp), and when the data is finalized, have the table auto-update? (and then, in theory, create a shp from that 'final' data that retains symbology, etc)? This must be possible- but i know i can't edit/replace a dbf that's part of a shp (without starting over), and it seems that edits can only be made in either ArcGIS OR in Access, but not both. thanks for ideas,


title: How to display the label with "Identify Features" in QGis?

text: I have been struggling for one hour with this simple problem : I have a raster layer. I would like to display the label and not only the band value when I use the "Identify Features" button. Is it possible? Thank you Best, Bap QGis 1.8.0-Lisboa ; Windows 7


title: Proper name for area defined by coordinate pair and radius

text: In defining geography objects within software, I have things such as a `City Area`, which is an area defined by the bounds of a city, or a `County Area`, defined by the county boundary. **What is the proper term for an area defined by a center coordinate pair with a radius?** I thought circle might be appropriate, but wasn't sure since that circle is bent along the earth's surface. Is there a proper name for this?


title: How to run a spatial query with a raster and vector layer in QGIS?

text: I have a soil map (vector-polygon layer) and a map of the crops farmed in the same area (raster layer, from the US Cropland Data Layer). I would like to calculate for a given soil type, the surface farmed per crop. How do I do that? The Zonal Statistics and ZonalStats plugins didn't solve my problem. Hope I am clear enough. Thanks


title: write multiple statistics analysis into the same output file with arcpy

text: How can I add multiple stats from the arcpy.Statistics_analysis to a single output file. Can I create a text file instead of a table ? Using one of the examples in help.arcgis, if I add other stat analysis, what is the best way to work with the output?               import arcpy     from arcpy import env     env.workspace = "C:/data/Habitat_Analysis.gdb"     arcpy.Statistics_analysis("futrds", "C:/output/output.gdb/stats1", [["Shape_Length", "SUM"]], "NM")     arcpy.Statistics_analysis("futrds", "C:/output/output.gdb/stats2", [["Shape_Length", "MEAN"]], "NM")     arcpy.Statistics_analysis("futrds", "C:/output/output.gdb/stats3", [["Shape_Length", "STD"]], "NM")


title: How accurate is approximating the Earth as a sphere?

text: What level of error do I encounter when approximating the earth as a sphere? Specifically, when dealing with the location of points and, for example, the great circle distances between them. Are there any studies on the average and worst case error compared to an ellipsoid? I'm wondering how much accuracy I'd be sacrificing if I go with a sphere for the sake of easier calculations. My particular scenario involves directly mapping WGS84 coordinates as if they were coordinates on a perfect sphere (with the mean radius defined by the IUGG) without any transformation.


title: Automatic calculation in one attribute table based on values in another attribute table using python

text: I am trying to fill the fields of one gdb table based on some statistics over the attributes of another shape file. As it was not possible in ARCgis directly, I am trying to use python coe for this case. I have table A where: Work_done Length 1 50 100 I would like to say that fill the legth field with sum of the Shape_Length in table B where "work_Done" field is equal another feild in table B. the following code is what I am thinking and I want to use it in the field_calculator but I got failed. I am not really pro in python therfore I would be really happy if sombody could guide me.                 import arcpy            from arcpy import env            # Set the current workspace        env.workspace = "C:/Python27/ArcGIS10.1"                 table = "Corriodr_stat.dbf"       ft = "Length"       fc = "Corridor_test.shp"       f1, f2 = "Corridor_S", "SHAPE_Leng"                 for row1 in sorted(arcpy.da.SearchCursor(fc, [f1, f2])):          for row2 in sorted(arcpy.da.UpdateCursor(fc, [ft])):             if fc.f1 <= 1:                table.row2 = getvalueField(sum row1.f2)             elif (fc.f1 > 1 and fc.f1 <= 60)                table.row2 = getvalueField(sum row1.f2)             elif (fc.f1 > 60 and fc.f1 <= 100):                table.row2 = getvalueField(sum row1.f2)


title: Get the specified feature by direction

text: I have two lines, while I want to get the line at the right hand according to a specified direction: ![enter image description here](http://i.stack.imgur.com/Zvw1m.png) The red line with arrow means the direction, while the line with blue color is the target to be selected. Is this possible?


title: removing spurious lines in raster before or after vector conversion?

text: I'm looking for a feature/function in either python opencv/ndimage/gdal/etc. or qgis/saga/grass that will help clean up a skeletonized image that I'd like to vectorize into a shapefile of lines/polylines. ![my messy skeleton](http://i.stack.imgur.com/lBqbJ.jpg) I spent considerable effort producing this network and I'd like to vectorize it and remove 'jaggies' (short lines to go nowhere) and also take the shortest path through the clumps of webby areas to produce as straight lines as possible. I realize this is a multi-step process and there may be possibilities in the python space or in gis. I know the steps I would do if I did it from scratch but it would take a while!


title: EM (electromagnetic) map to AWC (available water holding capacity) map

text: How do you use an equation from a regression analysis from Microsoft Excel to create a new available water holding capacity (AWC) map from an electromagnetic (EM) map. They are both to be in raster resolution. I have been told The equation describing the relationship can then be used to convert all EM values for grid cells from the map to AWC but I don't know how to do this in ArcMap.


title: Raster Tile Merge Performance of PostGIS 2.0

text: I'm attempting to intersect Australian Census boundary polygons with the 1" SRTM DEM (10.5bn points across Australia), and output the result as a single raster image on the fly. The DEM is stored in Postgres 9.2 x64 (with PostGIS 2.0) as 600x600 tiles with no overlap (i.e it's a mosaic) I've hit the known performance issue in PostGIS 2.0 using ST_Union to merge the intersecting full and partial raster tiles into a single record. For example: intersecting a 2,500 sq km polygon takes < 0.5s, which is excellent. But merging the resulting 18 rows into a single raster takes ~2 minutes. Is there any other function in PostGIS 2.0 to allow for tiles to be merged without having to do any comparison between the overlapping pixels (which don't exist here as it's a mosaiced dataset)? **Notes:**   * Plan B is to use gdal_merge in Python, but I'd prefer to keep this in the database for operational reasons.   * I've tried "ST_FasterUnion", thanks to this excellent post by Pierre Racine, but haven't experienced a significant improvement in performance   * I'm aware of ST_Union performance improvements coming in PostGIS 2.1, but don't have a spare dev environment to deploy and test.   * Postgres has already been optimised for my machine and has plenty of RAM and buffers to utilise.


title: How to find Closest Facility (euclidean distance)?

text: I have some facilities and incidents, and I want to find the closest facility for each incident. I've done it via Network Analyst (with a network dataset) but I'm wondering if you can do it with just Euclidean distance). I'm using ArcInfo 10.0.


title: Use subprocess.call with gdalwrap

text: I am now want to use subprocess.call with gdalwrap to clip raster from shapfile (in qgis 1.8.0) I use raster→projections→wrap(reproject)to generate command_line below               gdalwarp -t_srs "+proj=tmerc +lat_0=0 +lon_0=121 +k=0.9999 +x_0=250000 +y_0=0     +ellps=GRS80 +units=m +no_defs" -dstnodata 0 -q -cutline "C:/Program Files     (x86)/Quantum GIS Wroclaw/test/area.shp"     -dstalpha -of GTiff C:/Program Files     (x86)/Quantum GIS Wroclaw/test/wrapM_0.tif "C:/Program Files     (x86)/Quantum GIS Wroclaw/test/wrap"      and following is my code               import gdal     from subprocess import call          crs='"+proj=tmerc +lat_0=0 +lon_0=121 +k=0.9999 +x_0=250000 +y_0=0     +ellps=GRS80 +units=m +no_defs"'          mask='"C:/Program Files(x86)/Quantum GIS Wroclaw/test/area.shp"'          input='C:/Program Files(x86)/Quantum GIS Wroclaw/test/wrapM_0.tif'          output='"C:/Program Files(x86)/Quantum GIS Wroclaw/test/wrap"'          call(gdalwarp, -t_srs crs -dstnodata 0 -q -cutline msak -dstalpha -of GTiff     input output)      And python will give the error below                 File "<input>", line 1         call(gdalwarp, -t_srs crs -dstnodata 0 -q -cutline msak -dstalpha -of     GTiff input output)                                 ^     SyntaxError: invalid syntax      how to fix the error here??