title: Cannot Import Feature Class into Personal Geodatabase in ArcGIS 10.1

text: My objective is really simple. I have a shapefile and would like to convert it to a Feature Class. My goal is to use HEM tools, which I think requires me to convert my shapefile of interest to a feature class. I have a geodatabase target where I have been able to create **new** feature classes and feature datasets just fine. Importing a feature class is really problematic though. I get different errors depending upon the method used. Right clicking on the geodatabase gives me the "Cannot Create Output Error" implying that the output folder cannot be found....but I can see it, and the tool can see it as well or it would let me run. Right-clicking on the shapefile and exporting to a geodatabase gives me the generic 99999 error. I've tried changing path names to eliminate spacing, have tried writing to different locations to avoid permission issues (their shouldn't be any since I am an administrator). Nothing works. I'm just missing whatever it is that I need to change. Using the Feature Class to Feature Class tool works if I write to a folder...but then my end product is another shapefile, not a feature class. I new that would happen but just did it to test out writing to my Desktop Folder. Oddly, I've tried writing to geodatabases that I have written to before and I can't import feature classes to them either. Strange since they are populated with feature classes that I have imported in the past. I did just upgrade from 10.0 to 10.1...not sure if that has anything to do with it. Thanks for your help, Vin D'Angelo


title: Is there a compelling case to learn/use Python Toolboxes (new at ArcGIS 10.1) over Python Script Tools?

text: I've written a few Python Toolboxes (which are new at ArcGIS 10.1), but am yet to decide whether/when I should write them rather than Python Script Tools in a standard toolbox. I thought the **Online Help** might enlighten me when it prefaces some dot points with: > Once created, tools in a Python toolbox provide many advantages However, the five advantages listed all seem to be over not being able to use Python to write tools, and none seem to specify an advantage of Python Toolboxes over Python Script Tools. The two advantages that I can think of are:   * I can now write a "pure" Python tool in a single Python script without having to hook it up to a separately authored dialog with its Tool Validation looking like it was tacked on but I'm happy to be pragmatic rather than pure in this regard   * I could now use code (Python or any language capable of writing text files) to automate the writing of Python toolboxes but I am yet to come across a requirement to do this Am I overlooking the compelling case that led Esri to provide the Python Toolbox capability and, if so, what is it?


title: IMG and geoTIFF raster extraction

text: everyone. I'm very new to GIS technologies. Somehow, i need to extract raster(RGB data) data from geoTIFF and IMG files. What is the fastest way doing it?


title: ArcGIS Server Error - ClientCommunicatorAdmin restart

text: I am seeing some strange behavior each night with our ArcGIS Sever installation where it will throw the following error and become unresponsive.               Sep 01, 2014 4:03:06 AM ClientCommunicatorAdmin restart     WARNING: Failed to restart: java.rmi.ConnectException: Connection refused to host: 192.168.30.14; nested exception is:       java.net.ConnectException: Connection timed out: connect     Sep 01, 2014 4:03:27 AM RMIConnector RMIClientCommunicatorAdmin-doStop     WARNING: Failed to call the method close():java.rmi.ConnectException: Connection refused to host: 192.168.30.14; nested exception is:       java.net.ConnectException: Connection timed out: connect     Sep 01, 2014 4:03:27 AM ClientCommunicatorAdmin Checker-run     WARNING: Failed to check connection: java.net.ConnectException: Connection timed out: connect     Sep 01, 2014 4:03:27 AM ClientCommunicatorAdmin Checker-run     WARNING: stopping      I am still able to navigate to the the manager page, however, each service is listed as "stopping" or "stopped". Restarting the service for ArcGIS Server seems to resolve the problem until approximately the same time on the next day We have tried re-installing but the problem is still reoccurring.. Has anyone run into this problem before and, if so, what did you do to resolve it?


title: Is it possible to fix a vector layer after I badly georeferenced a map?

text: Sorry in advance if I missed another question on this topic. I searched for everything I could think of, but I didn't find anything with a solution to my problem. I am quite new to GIS tools and recently georeferenced a map of Germany. I didn't notice that I did a poor job on it though and although in some parts it is correct other's are way off. The problem I have now is that I created a vector layer with line objects based on the wrong georeferencing I did on the map. I know I can fix the original georeferencing quite easily, but is there a way though to fix the vector layer more or less automatically? Is there any way I can somehow adjust the position of the objects in my vector layer according to the new georeferencing of that map? I would very much like to avoid having to manually create all those line objects again, but don't know how I could possibly fix it otherwise. Thanks a lot to anyone who can help me with this, it's much appreciated! :) Edit: I forgot to mention, I am using Quantum GIS 1.8, but if it was possible with any other Freeware I'd be open to try that


title: Spatial adjustment for vector shp files in QGIS

text: How to do spatial adjustment for vector shp files in QGIS? I have a set of 4 cadastral maps of adjacent locations. I have converted them form DWG to DXF and then to vector shapefiles. They are not geo-referenced/ geo-rectified? I want them to be linked to a co-ordinate system (geo- referenced/ spatial adjustment) so that I can join the four parts together. I also have a shapefile of gps points corresponding to one of the four points. Can you please tell me how to do this in QGIS or GRASS, not ArcMap? Thanks, New_QGIS_user


title: How can I transform Real-time Local position system coordinates into GPS/WGS84?

text: I have a local coordinate system built with Ubisense sensors and tags (it is a real-time locating system used to automatically identify and track the location of objects or people). The local coordinate system is simply a cartesian coordinate system, which I can define the origin and axies at will except the rule that x-y panel must be horizontal, thus the z axis is vertical. After setting up the system, I can get coordinates of tags in the contained area where the sensors can get signals emitted from the tags and locate the tags. Now, I have the problem to convert the (x,y,z) coordinates into global system coordinates. I made a model myself which took the origin of local coordinate system as a reference point and calculate the difference of other points in lat and lon to make the transformation. However, I guess there must be more standard and accurate way to deal with such problem. Since I am a noob in this field, I did google something and read some pages. However, there is rarely content relevant with this topic(much I found is from some standard to another, however is not a standard coordinate system, it is cartesian but not geocentric). Anyone can give any hints or suggestions for this problem? Your help or attention to this question will be much appreciated.


title: ESRI 3D generate file format to ShapeZ

text: Can GDAL (ogr2ogr) convert ESRI 3D generate file format (class PointZ, LineZ and PolylineZ) to 3D Shapefiles? The format description can be found here: 3D Formats information. The PolylineZ format is similar to the GEOMETRY=AS_WKT, but not completely.               # GEOMETRY=AS_WKT     WKT,ID,     "POLYGON ((4377044.79 5585587.57 245.08,4377041.42 5585583.35 244.9,4377044.47 5585579.64 244.9,4377044.79 5585587.57 245.08))",78      ESRI 3D generate file format:               # PolylineZ format      78     4377044.790000 5585587.570000 245.080000     4377041.420000 5585583.350000 244.900000     4377044.470000 5585579.640000 244.900000     END      The Transformation of WKT to ShapeZ works.               ogr2ogr  -f "ESRI Shapefile"   export/poly3d.shp   ./import/wtk_poly.csv     Warning 1: Value 'POLYGON ((4377044.79 5585587.57 245.08,4377041.42 5585583.35 244.9,4377044.47 5585579.64 244.9,4377044.79 5585587.57 245.08))' of field WKT has been truncated to 80 characters.     This warning will not be emitted any more for that layer.      Any suggestions to convert the PolylineZ Format?


title: Problem with shooting star with pgRouting: doesn't respect the one ways

text: I used osm2po to create a table for pgRouting. I notice a problem: for a route with dijkstra or a* the route is correct and respect the one ways. But with shooting*, it doesn't respect the one ways. The route is from 11 baulacre geneve (46.2159082, 6.1394522) to coq inde cologny (46.2068106, 6.1919172). Here are the queries: dijkstra/a*:               SELECT topo.id, topo.osm_name, topo.geom_way     FROM Shortest_Path('SELECT * FROM ch_2po_4pgr', 263951, 271357, true, true)     AS path, ch_2po_4pgr AS topo     WHERE path.edge_id = topo.id;      shooting*:               SELECT topo.id, topo.osm_name, topo.geom_way     FROM Shortest_Path_shooting_star(        'SELECT *, geom_way as rule, reverse_cost as to_cost FROM ch_2po_4pgr', 321089, 309056, true, true)     AS path, ch_2po_4pgr AS topo     WHERE path.edge_id = topo.id;      Did you notice the same problem with shooting* in others situations? Here is a picture. The orange line is for a*. The red line is for shooting*. Since with shooting* with use edge_id instead vertex_id, maybe the algo starts at the beginning of the way (starting point is in pink), but the problem with the one ways is in the pink circle: the red line uses the "rue des Alpes" which is a one way. ![enter image description here](http://i.stack.imgur.com/d7jZN.png) Thanks to all for your help!


title: Can ArcMap automatically highlight an object on the map its label?

text: I was wondering if there is a way to highlight a feature (e.g. line) only by selecting its label on the map. Can ESRI autmatically highlight this object or do I need to create some "connection" between line and its label?


title: 'Iterate Field Values' for Multiple Fields in ModelBuilder

text: I have a model that iterates a process for each unique field value of [Type] field. The model aggregates polygons of a like Type, and at the end of the model a new field is created and populated with each unique [Type] for the aggregation using "%Value%" from the output of the iterator. These iterations are then appended and contained within a newly created feature class. I would like to retain this attribution for each set of polygons(original [Type] value), but I would like to include another set of values to the iteration. That is, instead of iterating through each unique [Type] value, I would like to iterate through [Type] AND [Grouping]. The [Grouping] field is a simple Short field, containing only two unique values. Ideally, the aggregation would aggregate polygon clusters with like [Type] AND [Grouping] values, yet attribute the polygon with a single [Type] value, as well as a field with a [Grouping] value. The problem that I am running into is that if i concatenate the two fields in before running the iterator based on the concatenated field, the value output is concatenated, and my "%Value%" output no longer is correct for populating the field at the end of the model. Is there a way that I can use an iterator to select [Type] AND [Grouping], yet call on the [Type] value of the selection to populate the field at the end of the model? I can post images of the model if it would help. I'm using ArcGIS 10. Thanks, Eric


title: How to find nearest point for each point in current layer?

text: I have a point layer and want to know the distance of nearest point of current layer for each point in the layer. I found some tools in ArcToolbox about it , but they need two layers for input and they can find the distance of nearest object between that two layers. I tried to use them and i added my point layer in that tools ( Same layer for both inputs ) but they couldn't work properly because they found same point in second layer and gave me zero distance for every points. How can i find distance of nearest point for each point in my layer ? It is better for me to find an answer in ArcToolbox . Something such as Tools or Models. Thanks alot


title: Repairing broken links for new paths of layers to shapefiles using ArcPy?

text: I have an .mxd with broken link layers and each layer has a different new path. I know the python script for repairing the broken links of the layers, but I want to put all the new paths of the layers and then let python do and finish all the process. The script that I have does it one at a time. I know if all the shapefiles are in one path then it's done by performing it once. but what if they al have different paths? To clarify, let's say I have 3 layers with broken links and each one has a different new path now. I'd like to have a GUI with 3 inputs for the new paths altogether or let script do the process one by one and each time asking for the next path to input.


title: Deep merge data from multiple layers

text: Let me explain what I mean by "Deep Merge" and what I am trying to do. Suppose that I have 3 layers (Layer1, Layer2 and Layer3) with the same field names (field_a, field_b and field_c). For any record, each layer has a value only in one of fields. For example, One row in Layer1 has a non-empty value only in field_a. One row in Layer2 only in field_b. One row in Layer3 only in field_c. The tricky part is that, different rows in the same layer may not have a non-empty value in the same field. My goal is to merge all 3 layers into one target layer, which contains all non-empty values from all 3 layers. In my mind, it involves some serious loops or joins. So I am wondering if anybody knows any existing tool (open source or commercial) to do that. Thank you very much!


title: cannot reproject postgis in mapserver

text: I am trying to set up a wms server using point data in a postgis (v 1.5.1, pgsql v 8.4) database and mapserver (v 1.5.1) (I know these are not the most current versions of anything, but my server is running debian stable, if neccesary, I'll upgrade them). The data is a global data set, so the points are stored in lat/lon as postgis geography data type. I manage to display the data (using either arcmap 10 or Qgis 1.8 or directly in a web server) in lat,lon, but I would like to be able to serve the data in other projections as well. I have tried to set up LAEA north pole (EPSG 3408) as well. In the GetCapabilities, I get a part like this:               <Layer>         <Name>xxxx</Name>          <Title>yyyy</Title>          <Abstract>zzzz</Abstract>          <CRS>EPSG:4326</CRS>          <CRS>EPSG:3408</CRS>          <EX_GeographicBoundingBox>            <westBoundLongitude>-50</westBoundLongitude>             <eastBoundLongitude>90</eastBoundLongitude>             <southBoundLatitude>-50</southBoundLatitude>             <northBoundLatitude>90</northBoundLatitude>          </EX_GeographicBoundingBox>         <BoundingBox CRS="EPSG:4326" minx="-50" miny="-50" maxx="90" maxy="90" />      </Layer>      I.e my two projections are listed, but I only get the bounding box for the original projection. SHould I also get the other one? As far as I can see from the instructions I can find, I just need to list the projections I want to use in wms_srs in the METADATA section - in the layer, in the MAP or both places (I've put it both places, but I've also tried it in just the layer or MAP, with the same result)


title: Is it possible to create a point within a Polygon (not a centroid)?

text: Is it possible (without writing Python code - sorry I'm not a developer) to create a point layer from a polygon layer but ensure the point is within the polygon? The ftools centroid tool is cool but some centroids are created outside the polygon (i.e. irregular polygon). In ArcGIS there is a "Feature to Point" tool similar to the ftools 'polygon centroids' tool but has a checkbox to force the point to be within the polygon.


title: How to import R-generated .gri/.grd raster files in GRASS?

text: I received data in .grd + .gri format, where the first one is a small text file and the second one (having the same name) is a binary file. I am unable to open them in GRASS/QGIS. Any hint ? Missing library ? I tried to use Google but unfortunately "r" word is too small, and it ends to refer always to the r.* commands in GRASS... This is the content of the .grd file:               [general]      creator=R package 'raster'      created= 2013-12-14 15:44:30      [georeference]      nrows= 134      ncols= 143      xmin= 56000      ymin= 1613000      xmax= 1200000      ymax= 2685000      projection= NA      [data]      datatype= FLT8S      byteorder= little      nbands= 5      bandorder= BIL      minvalue= 0:0:0:0:0      maxvalue= 0.945:1:1:1:1      nodatavalue= -1.7e+308      [legend]      legendtype=       values=       color=       [description]      layername= Value:Twin:Psum:pH:CN      history=        Thank you...


title: ArcGIS Python Field Calculation

text: Ok so i've got a problem with the field calculator in ArcGIS: I'm already pretty far with it since i got some help :P but it still doesnt work. What i've got: -A complete table of several states of the world with their specific life expectancies. -A world map (not yet joined) What i have to do: classify the states in 5 groups by their life expectancies in an extra field in the field calculator. What i did: I generated the following code block:               def reclass(a):         if (a <= 50):             return 1         elif (a <=60):             return 2         elif (a <=70):             return 3         elif (a <=75):             return 4         elif (a >75):             return 5         else             return 999      =               reclass(!LIFE_EXP00!)      (I'm not sure if u can see it in this post but i was respecting the indentation. I added a new field called "Field_1" and i took LIFE_EXP00 from the already given data in the table. But still ArcGIS tells me after every single time i try it, that there was a failure during processing. What am I doing wrong? Please help me.


title: Missing Georeferencer

text: I installed Qgis Master 64 bits on win7 64 16gb of ram and I am unable to locate the georeferencer. All other Gdal application are there, Gdal is active, but no georeferencer. I tried yesterday, tried also reinstall and also I anstalled again today. No georeferencer at all. Any Ideas?


title: Creating MySQL database for storing imported KML?

text: I need to create a MySQL database for storing coordinates from a KML file with the structure as follows:               <Placemark>             <name>someplace</name>             <styleUrl>#m_ylw-pushpin</styleUrl>             <Polygon>                 <tessellate>1</tessellate>                 <outerBoundaryIs>                     <LinearRing>                         <coordinates>                             73.13693971060766,19.25310233594502,0 73.13693833446244,19.25310233595927,0 73.13693403840304,19.2530984545477,0 73.13693244215831,19.25309329030865,0 73.13693090444444,19.25308942538022,0 73.13692794507664,19.25308428308517,0 73.13692494004152,19.25307787269911,0 73.13692467634385,19.25307148170985,0 73.1369258304161,19.25306638278946,0 73.13692698179315,19.2530612962149,0 73.1369322686994,19.25305748929283,0 73.13693618566978,19.25305368927265,0 73.13694020972429,19.25305242409976,0 73.1369455322031,19.25304989609698,0 73.13695090854537,19.25304863322264,0 73.13695628164153,19.25304737111269,0 73.13696294560927,19.25304484915161,0 73.13696837570761,19.25304484914826,0 73.13697652081103,19.2530448491438,0 73.13698202045055,19.25304610973631,0 73.13698895342023,19.2530486332095,0 73.13699582091516,19.25304989604919,0 73.13700004851012,19.25305242405945,0 73.13700428114488,19.25305495514516,0 73.137008596689,19.25305875746043,0 73.13701292005057,19.25306256664898,0 73.13701733179597,19.25306765632393,0 73.1370191023439,19.25307403574487,0 73.13702208277293,19.25307787264039,0 73.13702249452535,19.25308428301952,0 73.1370229075506,19.25309071284612,0 73.13702323885482,19.25309587081664,0 73.13702227801461,19.25310233591246,0 73.13702260952107,19.25310752220441,0 73.1370201827993,19.25311272114856,0 73.1370177501589,19.2531179327989,0 73.13701531155583,19.25312315721282,0 73.1370114025442,19.25312708392549,0 73.13700602252482,19.25312970574124,0 73.1370006359243,19.25313233078209,0 73.13699509258757,19.25313233079353,0 73.13698816344083,19.25313233078499,0 73.13698116356548,19.25313101786714,0 73.13697555292208,19.25312970576528,0 73.13696994578734,19.2531283944655,0 73.13696434201735,19.25312708395525,0 73.13695999711467,19.25312315723566,0 73.13695697863325,19.25311793285053,0 73.13695264975624,19.25311402288048,0 73.13694832879452,19.2531101201258,0 73.13694395681273,19.25310492752692,0 73.13693971060766,19.25310233594502,0                          </coordinates>                     </LinearRing>                 </outerBoundaryIs>             </Polygon>         </Placemark>         <Placemark>             <name>place2</name>             <styleUrl>#m_ylw-pushpin</styleUrl>             <Polygon>                 <tessellate>1</tessellate>                 <outerBoundaryIs>                     <LinearRing>                         <coordinates>                             73.13723765697174,19.25290222725077,0 73.13723989677642,19.25289857411988,0 73.13724050975087,19.25289250009896,0 73.13724140951585,19.2528888642842,0 73.13724216325605,19.25288402638158,0 73.13724690239067,19.2528791999123,0 73.1372491215834,19.25287558746117,0 73.13725413788218,19.25287318265183,0 73.13725944636737,19.25287318265056,0 73.13726740917318,19.25287318261767,0 73.1372740448365,19.25287318261646,0 73.13727818083784,19.25287438462484,0 73.13728161429077,19.25288040541489,0 73.1372835696157,19.25288523471525,0 73.13728286450872,19.25289007539774,0 73.13727964512475,19.25289614231522,0 73.1372761061091,19.25289979111273,0 73.13727389776194,19.25290344626298,0 73.13727065467811,19.25290955268617,0 73.13726828301977,19.25291200031773,0 73.13726337940838,19.25291567720809,0 73.13725846706093,19.25291936058711,0 73.13725324778412,19.25292058984915,0 73.13724877730209,19.25291690427552,0 73.13724550831267,19.25291200035216,0 73.13724224708763,19.25290710800144,0 73.13724152171221,19.25290100884397,0 73.13723765697174,19.25290222725077,0                          </coordinates>                     </LinearRing>                 </outerBoundaryIs>             </Polygon>         </Placemark>      How do I store the above in MySQL database? I need to store coordinates of polygons. say FE: in kml file, i have highlighted 2 areas. I may required to calculate area(sq. mt) or distance from point a to point b. add custom placemarks with data generated by database


title: Test scripts for google maps using selenium tool?

text: I want to test Google maps using selenium but since selenium works based on object identification properties whereas Google maps does not have DOM elements.Anybody here who can help me out to perform zoom in zoom out operation on maps?


title: NYC Subway Line Station Order?

text: I've currently got data on the NYC subway stations (including name, position, and how many people entered and existed during a timeframe). Also, I have the data on which stations have which train lines. However, I don't have the order of these stations on the line. Is there anywhere where this data is already formatted in some way that can be easily used programmatically? I could put them in order in some csv file on my own, but I'd rather avoid that if I can. Thank you much!


title: Can osm way records in SQL spatial be merged to one geographic record if they have same Id?

text: I've imported England shp file of roads to mssql using OSM2MSSQL and I've noticed that the same road has couple of records with same Id but different geographic (probably it was segmented). Is there a way to convert these records to one record with entire geographic data? I want to have unique Id in my streets table instead of having two tables, one with unique Id and one with geographic data.


title: Different typesets for labeling within one project

text: I'm searching a way to show labels with a western typeset, as well as a differnt one. In my case thats western (English/German/...) PLUS Russian, which I don't speak at all. But I'd like to be able to display the native labels - for example names of cities. My idea was to use a second layer and join both, but I cant find a way :-( Not that I cant find a way to join, but have cyrillic displayed :-/ QGIS doesn't seam to save/recall if the typeset is changed. Layer->Settings->Labeling->Font->Typeset Please "suggest" a menu path - I'm not a native speaker of either, nor do I use the English version of QGIS.


title: Does the Calculation of Spatial Weights Matrix in ArcGIS take into account 3D?

text: This question is about Spatial Weights Matrix in ArcGIS. Does it take into account 3D or not? I want to calculate Spatial Weights Matrix with SPACE_TIME_WINDOW. Anton


title: Is there an equivalent of QgsMapToolIdentify in QGIS 1.8?

text: I'm using quantum gis 1.8 Lisbon, and I need to select a point that is on the map by mouse, to know their attributes. But I can not use the QgsMapToolIdentify function, that would show me another function that would help me to activate the layer by mouse. This code are thus not working               myLayer iface.activeLayer = ()      myLayer.selectedFeatures objects = ()      objets.attributeMap ()      object = objects [0]      object.attributeMap attributes = ()      attributes [0]. toString ()      thank you * * * We introduce points in the field of view of qgis, I plan to get information from one point when I click the mouse on that point. There is QgsMapToolIdentify function, but I'm using quantum gis 1.8 Lisbon and this function does not exist. So I have to call another function, but do not know which. thank you


title: Mapinfo problem importing mid/mif

text: I am trying to import a mid/mif into MapInfo. I follow the usual routine of Table --> Import and select the correct mid/mif. I get the following error: (f:\FILENAME.mif:8) Found [(] while searching for [An identifier] I have searched online and found no mention of this error. Any help?


title: How do I stop event propagation with rightclick on Leaflet marker?

text: I didn't find a rightclick event in Leaflet markers so I started listening to the mousedown event. In the method that resolves this event I check               if (event.originalEvent.button == 2) {         do_something();         // now i want to stop propagation of the rightclick         return;     }      This does work in that do_something() is called, but the event propagates to the map and a contextmenu is being shown. How do I prevent the event to propagate to the map? Does this work in all major (i.e. IE) browsers?


title: Getting buffer wizard-like speed in Python (vs using MultipleRingBuffer_analysis)

text: I've got a script that takes a single input polygon feature, throws 10 buffers around it at tenths of a specified input distance, and then symbolises the output based on an existing layer file. Pretty simple. However, the arcpy.MultipleRingBuffer_analysis operation is amazingly slow. It takes upwards of two minutes to generate the buffers, even for very basic polygon inputs - the same result can be had in about two seconds by using the buffer wizard tool. Problem is, the buffer wizard can't be accessed through arcpy. So obviously it's _possible_ to quickly generate multiple ring buffers - does anyone have any insight as to how the buffer wizard tool is doing it, and how that might be replicated in Python?


title: Select layer by attribute using SQL statement LEN() in ArcMap

text: as the title already suggests... I want to query the length of a string column and thought to use: LEN([ColumnName]) > 5 However, that doesn't work at all. I'm querying a table in a fGDB in ArcMap 10.1, hence the notation of [ColumnName]. Can someone help me? The error message I'm getting from ArcGIS is that I'm using an invalid SQL statement and that does not really help. Cheers, Tom


title: How to re-calculate area, perim, p2a,width and length in QGIS?

text: I cut a part out of one feature and added it to another feature in a polygon shapefile. My problem now is that both features has the same values in the area-, perim-, p2a-, width- and length-columns? The features are very much different in size, how do I get QGIS to sort of "re-calculate" the areas etc? Greatful for any help!


title: Turn off field programmatically in ArcMap via VB.Net

text: I'm creating an addin button in Arcmap via VB.net that when clicked, runs the identify tool. But I want only selected fields to appear in the identify dialog. My idea is to hide the unnecessary fields and then show them again once the identify tool is deactivated. Is it possible to programmatically turn off/on a field?


title: Remove non-alphanumeric characters from a string in Mapbasic

text: I have a user-entered string which I am using as a search term. I would like to display the results in a table named `"SearchResults_"+SearchString`. However the string could feasibly contain characters that can't be used for table names - spaces would be the most common and there are plenty of other strange characters in place names. **I would like to strip all special, non-alphanumeric characters from a string.** As far as I can see there aren't any built-in functions that would achieve this. I can loop through the characters easily enough but comparing each one to a white-list is proving harder than I expected. I could do it with a epic if/else statement but that's just dirty! VB.net has IsLetterOrDigit() which would be nice. PHP has ctype_alnum(). In Python I'd just check it against a white-list string. Any ideas?


title: Free tool to convert .lpk to .shp file

text: I have a .lpk file from ESRI that I would like to convert to a shapefile. How can I do this, using only free software?


title: Practical way of managing polygons and unions of them

text: I have a table in PostGIS with a quite fine partition of a region, which I'd like to dissolve at two higher levels by means of two foreign keys (e.g. municipality and state). The partition is a real one, i.e. the areas don't overlap, they just touch each other. I thought it'd be nice to do that in views, to follow the "don't insert into the database something which can be calculated from other existing data". However, a simple               SELECT id_municipality, ST_Union(the_geom)     FROM areas     GROUP BY id_municipality      takes around 22s on a core i7 for "just" around 36000 areas with a POLYGON geometry field. I indexed the geometry field with a GIST index, but with no luck. I am now thinking that the only practical way would be to insert a geometry column in the municipalities table and calculate it beforehand. Unfortunately, I need to share my results to less technical people, so installing experimental packages (like topological extensions) is not practical.


title: Create a Layer in GeoServer using REST

text: As part of a possible solution to this question i want create a layer in GeoServer via REST that supplies a SQL statement along with a userid. Reading I think i would create a feature type and publish it? I need to do this all in REST. Does anyone have any code examples?


title: How can you map counties to congressional districts?

text: Given county A, I need to know what congressional districts are within it.


title: Using IDocumentEvents_Event from a .NET 4 assembly within ArcMap

text: Are there a known issues when working with events from .NET 4 COM assemblies within ArcMap? For example, when using IDocumentEvents_Event from a .NET 4 assembly in ArcObjects, the event appears to work properly, but disables the standard right-click context menu behavior in ArcMap. Changing nothing in the assembly except the target to .NET 3.5 or .NET 2.0 will "restore" the default behavior. Is this a COM issue, or an ArcObjects issue? If the latter, is there a potential workaround?


title: topology rules for structuring topology data to get line feature as topoparent in topology graph

text: i Created topology graph of a polygon feature and i copied geometry of topology edges into line feature class,and topology nodes into point feature class i created topology layer containing featureclasses as polygon,line,point . using IEnumtopoparent i am getting FID of individual points ,lines. i want to get polygon arclist,for getting this using topologyGRAPH.getparentedges(polygon feature class,OID) i am accessing individual edges of polygon. using topoparent.m_FID iam getting fid of parents. but iam not getting line feature as topoparent to get line feature as topoparent for tpology edges what topology rules should I applY when creating topology layer pls helpe me.....


title: PostGIS template corrupted

text: Somehow my postGIS template got overwritten by one with data from one of the databases I previously created. I've tried everything I read on goggle but could not delete or replace the template with the correct one... I am hoping that someone here would know the best way to fix this problem without having to do a complete a database software. Thanks Chris


title: Geotools throws "Cannot encode multilinestring (yet)" while Shape import on Linux

text: I've encountered a strange platform-dependent SHP import problem. Everything works fine on Windows, shape files are converted fine to database. There is although small little warning: > [[ACTIVE] ExecuteThread: '9' for queue: 'weblogic.kernel.Default (self- > tuning)'] WARN geotools.jdbc - No mapping for > com.vividsolutions.jts.geom.MultiLineString Obviously SHP contains Multilinestring. So why the same code on Windows works, but on Linux we've got exception: > UnsupportedOperationException("Cannot encode MultiLineString (yet)") Would appreciate if someone could help. Cheers!


title: PostGIS import csv encoding problem

text: I want to import csv into PostGIS database, I receive error massage:               alpha=# COPY project FROM '/home/alpha/database/project.csv' DELIMITERS ',' CSV;     ERROR:  invalid byte sequence for encoding "UTF8": 0xc1da     CONTEXT:  COPY project, line 1      so I find my csv's coding is `ISO-8859-1` and set `encoding = ISO_8859_1`               #SET client_encoding = ISO_8859_1;      COPY worked but I use               SELECT * FROM project      it's give me a lot of garbled words! so how to make it right??


title: Where can I get shapefiles for Indian Constituencies according to the latest delimitation?

text: The GIS data officially provided by the government is only available as PDF maps. Some people have shared shapefiles for the previous delimitation but these are outdated. Parliamentary constituency data can be also be viewed here but no source has been provided. Would anyone know a way to get GIS data for the assembly and parliamentary constituencies?


title: Points warped when exporting to shape file (position off in Qgis when imported)

text: I have just purchased a used Trimble GeoXH and I have done some surveys. I did the differential corrections in Trimble Pathfinder Office 4.10. I have exported the corrected and uncorrected files and have imported them into QGIS 2.0.1 and both file's points are distorted - horribly. I assume this is a projection problem on export. When I look at the points in Pathfinder Office, they are extremely accurate - even uncorrected. It seems I am doing something wrong on the export. I have terrasync pro 3.30 and have set the Datum to NAD 1983 (Canada) in TerraSync on the gps - I am in Western Canada. I have the project CRS set in QGIS to NAD83 UTM Zone11N EPSG:26911 (this is where the data was captured - BC / Western Canada) and I have "ON the fly CRS Transformation" turned on. I had it set to EPSG:3857 as stated in another post for use with the google hybrid plugin, but that is still warped. I have imported data prepared from a government source into QGIS, and everything displays properly. I would appreciate some advice on how to correct this issue. I tried exporting to a ESRI shape file using the EPSG:26911 projection file I downloaded from spatialreference.org. If I am using a faulty projection file or there is a better export format for me to use, I would appreciate some guidance.


title: Is there a ArcGIS python script development API?

text: I'm a somewhat experienced software engineer new to the GIS world. I started working on a GIS project that supports ArcGIS 9.3 and 10.0 and looking through existing code it looks like the main interface library is `arcgisscripting`. However, a google search of "arcgisscripting API" or any other variant I can think of only turns up blog posts, forum threads, or broken links. Is there a core API or reference for `arcgisscripting` or writing Python scripts for ArcGIS? Am I doing it wrong and need to try something else?


title: Align and georeference GIS layers exported from CAD

text: I have two separate CAD files of the same locale. The first one has surface data like street curbs, building footprints, etc. It is not georeferenced, but when I add it to other items in a projected (NAD 1983 FIPS US feet Connecticut) data frame, it aligns perfectly. The other CAD file contains subsurface utilities and a surface basemap as a reference. When I try adding them, they show up very very small not too far from the properly projected data. When I try to project subsurface layers, they don't show up at all. Trying to georeference manually gives me an error after I manually attempt to enlarge it using georef toolbar tools, even if I delete the spatial index. I have ArcGIS 10.1 SP1 and a trial of CAD 2014. Screenshot displays properties boxes, the surface data on the right, subsurface on the left. ![enter image description here](http://i.stack.imgur.com/CMT0X.jpg)


title: Flex Viewer: How do I return a subset of the input data using the geoprocessing widget?

text: I have a set of point data (representing crashes) that has been published as a feature service and which I have pulled into a Flex (3.5) application. I want the users to be able to submit a query and create a subset of the original data - for example, to access all the crashes that occurred on a Tuesday, or which were fatal, etc. To do this I have created a Geoprocessing widget based on a very simple model (Get expression, Select by attribute; copy features) that has been published as a GP service. However, I haven't been able to get the GP widget to work properly. Once run, it adds the output 'Query' layer to the layer list, but there are no features within that layer (ie, nothing on the map). Geoprocessing Widget config file:               <?xml version="1.0"?>     <configuration label="Select by Attribute">         <taskurl>http://ags-dev2:6080/arcgis/rest/services/CAS_Crashes/SelbyAttr/GPServer/SelbyAttr</taskurl>         <shareresults>true</shareresults>         <inputparams>             <param name="Expression"                    label="Expression"                    required="false"                    type="string"                    visible="true">             </param>         </inputparams>         <outputparams>             <param name="Output"                    label="Query"                    type="featurerecordset"                    visible="true"                    geometrytype="point">                 <simplerenderer>                     <simplemarkersymbol color="0x2D9C49" width="4"/>                 </simplerenderer>             </param>         </outputparams>     </configuration>      Does anyone have any ideas about what I might need to change to bring the features through to the output layer? I'm using ArcGIS Desktop 10.2, ArcGIS Server 10.2, Flex 3.5. Thanks!


title: Total station problem

text: **I have total station data in local coordinates in excel format. Now i want to transform it into global coordinate (in utm) to correctly overlay on the actual place. What is the method?**


title: relative path for image elements in ArcMap

text: Is there a way to set relative path for image elements (like a company logo) in ArcMap projects? The classic "Store relative path" option in the Map Properties doesn't seem to be working. Thanks!


title: Cascading GeoWebCache Servers

text: I am trying to address some latency issues accessing US based servers from Asia We are using GeoWebCache in front of two GeoServer instances, based in the USA. One thought that occurred to me was to deploy additional GeoWebCache Servers, but instead of having them hit a GeoServer instance, locally or remotely, have them front another GeoWebCache instance (Cascaded) located in the USA. The thought behind this is that the tiles will already be cached on the "Master" GeoWebCache which simply serves the requested tiles to the GeoWebCache instance requesting it in Asia. This should eliminate the need to have additional Geoserver Instances to service Tile requests Has anyone done this before or have any thoughts as to whether it will work as I expect?


title: Fastest way to test if raster contains data within a polygonal area

text: I am writing a python script for a geoprocessing tool in arcmap. I have a polygon in WGS 1984 Web Mercator (Auxiliary Sphere) and a raster grid, raster1, in WGS 1984. I reproject the polygon to create a mesh of points that are contained within the boundary of the polygon and laid out in a uniform grid. Now, the goal is to check raster1 at each point in the mesh and detrimine if the corrosponding cell has a value at that xy. If it does remove the point from the mesh. If it doesn't, keep the point. In the end I want to have a fragment of the mesh of points which contains all points for which the raster had no value. My first attempt was to loop through the points and use arcpy.GetCellValue_management() to see what the value was at each point. But the GetCellValue_management fuction seems to take a painfully slow 2-3 seconds (I have a lot of points). Anyone know why this takes so long? An identify on the same raster in arcmap is almost instant! So now I'm thinking about using Clip_management to clip out the intersection of the polygon and the raster and then RasterToPolygon to get a polygon representing that intersection. Then finally using Erase_analysis to erase from the point mesh those points in the intersection. What do you guys think? Is plan two good? Whats the fastest way of doing this?


title: VB Script for selecting multiple records

text: > **Possible Duplicate:**   >  "If then" VB Script code (or Python) equivalent to SQL's IN ('x','y',z') > expression in ArcGIS Field Calculator question? going to ask this again (hopefully more specifically) Please see screen shot linked below for hopefully better explanation. VB Script does not recognize SQL code "in (x,y,z)" for selecting multiple records. Have searched high and low but still cannot locate the VBScript version for this and am starting to feel that it does not exist which cannot be possible. Please help.![enter image description here](http://i.stack.imgur.com/TRCTy.jpg) ![enter image description here](http://i.stack.imgur.com/clwP8.jpg)


title: Address locator not using modified .cls file

text: I have made a customization to the US Streets locator style for our application. We've changed "SR" to map to "STATE ROAD" instead of "HC" in the us_addr.cls file in the C:\Program Files (x86)\ArcGIS directory. When I run StanEdit, it maps "SR" correctly. When I do a test in ArcMap using the Geocoding Standardize Addresses tool, picking the the US Streets with City, State, Zip and AltName, it works correctly also. I've created a locator in ArcCatalog using the US Streets with City, State, Zip and AltName style and registered it on my local server (which is the same machine I did the above testing on). When I run the batch geocoding tool using that locator or use the Find Address tool, it maps "SR" to "HC". So I'm thinking that I'm not editing the correct us_addr.cls file. But I have no idea where the locator is pulling from. Is there any documentation anywhere (besides the Geocoding Rule Base Developer Guide - been all thru that) that explains how to modify the locator syles? We are running ArcGIS 9.3.1 SP 2.


title: Ruby bindings to GDAL

text: I'm looking for a guide on how to get up and running with ruby and gdal on Mac. I've seen some examples on how to use gdal in ruby but never found any tutorials on where to download and install it. Any ideas highly appreciated!


title: TRI, TPI, or Roughness

text: I am working on a site selection projection where I will be unable to perform any in situ observations. A key component of the selection process is the ruggedness of the underlying terrain. I plan on using GDALdem to create the derived ruggedness maps. I have read both the Wilson, et al. 2007 paper and the GDALdem documentation. While both offer insight into the algorithms used to generate the indices, only the Wilson paper offers an assessment of index suitability. It is essential, because I can not perform in situ observations, that the index selected over estimate the surface roughness. In your experience, which index have you utilized, for what application, and why?


title: Redistribute linestring vertices to specified distance

text: I use a groundwater model program called GMS, and it has a very useful function "Redistribute vertices" (more info). I use this function all the time to both simplify and provide _regular_ spacings between points (vertices) in linestrings. See this image from GMS where the black lines are original stream with irregular spacings, and the grey lines have redistributed vertices using 50.0 m ![from GMS](http://i.stack.imgur.com/XMRRx.png) I haven't seen the equivalent function anywhere else in the GIS-processing world (Esri, JTS, PostGIS, etc.). Redistribute vertices requires a specified spacing distance, whereas simplify-type functions use a tolerance and don't use a specified spacing. Does anyone know if there is something equivalent in either ArcGIS or open source GIS projects?


title: How do you add a search box to the map when creating a CartoDB visualization at runtime?

text: I know that it is possible to add the CartoDB search box when you link to a visualization such as this. However, I have a map which is created at runtime, using the same method as here. Does anybody know, when a visualization is created at runtime, if there is any way to add the CartoDB search box?


title: Why is my north arrow pixelated when being imported as emf?

text: I've tried to import my own North Arrow with the Style Manager (ArcGIS 10.1) both as .emf and .bmp but the North Arrow is pixelated although the resolution of the picture should not cause this problem (i've saved the Arrow in a DinA3 file) and although the vetor format of the picture. Has anyone an idea how to solve this problem?


title: How to use paging in a WFS query?

text: I am running a `GetFeature` request against a WFS server which does not support to download all data at once. Can I use `ogc:PropertyIsGreaterThanOrEqualTo` to split the dataset into chunks at download them step by step? I wonder how I can construct the actual URL to include the paging filter. Here is what I tried:               http://example.com/wfs.aspx?request=GetFeature&service=WFS&version=1.1.0 \     &typeName=example:example&maxFeatures=50000 \     &FILTER=<ogc:Filter xmlns:ogc="http://www.opengis.net/ogc"> \     <ogc:PropertyIsGreaterThanOrEqualTo><ogc:PropertyName>OBJECTID</ogc:PropertyName> \     <ogc:Literal>50000</ogc:Literal></ogc:PropertyIsGreaterThanOrEqualTo></ogc:Filter>      This is not a valid URI though... That's why uri-encoded the filter part:               http://example.com/wfs.aspx?request=GetFeature&service=WFS&version=1.1.0 \     &typeName=example:example&maxFeatures=50000 \         &FILTER=%3Cogc%3AFilter%20xmlns%3Aogc%3D%22http%3A%2F%2Fwww.opengis.net \     %2Fogc%22%3E%3Cogc%3APropertyIsGreaterThanOrEqualTo%3E%3Cogc%3APropertyName%3EOBJECTID \     %3C%2Fogc%3APropertyName%3E%3Cogc%3ALiteral%3E50000%3C%2Fogc%3ALiteral \     %3E%3C%2Fogc%3APropertyIsGreaterThanOrEqualTo%3E%3C%2Fogc%3AFilter%3E      This actually works. Is there a "nicer" way to do this?


title: projection change to create buffer in QGIS

text: I am about to become desperate... I like to create buffer around several points - if possible - with a declaration of the size in km of these buffers. Just to show you what I am working with: ![enter image description here](http://i.stack.imgur.com/iMoTL.jpg) My problem seems to be the projection of my data. Its in the WGS84 projection and my background map is in the RGF93 / Lambert-93 projection (I dont know why that works together). With the WGS84 projection the distances will be shown in degrees, so i tried to change the projection of the background map, the projection of the data or both together but I dont find any other projection where they fit together. Do you have any suggestions how to create these buffers with a clearly shown distance?


title: MSSQL server 2008 dissolve sql query

text: I have a set of spatial objects in a MSSQL server 2008 R2 trying to dissolve adjoining objects where one of the attributes are equal. The table below is somewhat representative to my original table:               id     year     GEOM     01     1999     0x79F...     02     2013     0x79F...     03     1988     0x79F...     04     2013     0x79F...     05     1999     0x79F...     06     2013     0x79F...      Many of the objects above are neighbours so I would like to aggregate the spatial objects so that I in the end of the day end up with a smaller number of table entries/objects where possible. I am not looking for myltipolygons. I have looked around and a propper dissolve functionality like the one available in for example ArcGIS desktop. MSSQL 2008 R2 does not seem to sport this functionality directly. Is there a way to send a query to the database which does this job for me?


title: How to aggregate values multiple from multiple points along a route as an attribute of it?

text: I have a route network with a "Name" field that identify every route. On the other hand I have a shape file of points. Each point corresponds to the number of people using a route at a time t. I would like to apply those counts to the routes but I do not know which tool to use, knowing that some routes have more than one count along them. I have interpolated using the inverse distance weighing but it is giving values over the entire areas instead of just along the trails. I also tried to use calibrate routes but it does not seem to be working :s Anyone has an idea? Thanks a lot


How to get a kml from stored in a blob file of oracle with arcpy? i try get data or kml files contains in a database oracle in a field of a table, using cursos options of arcpy.               import arcpy, os     from arcpy import env     source =r'C:\connectionFiles\file.sde\tableoracle'     cursor = arcpy.SearchCursor(source)      ... the conection work corretly, but the cursor's no work weel for this type of information .... any idea, THANKS FOR YOUR ANSWERS :)


title: Is it possible to remove (Actions) and (Derived) from the QGIS Identify dialog?

text: When identifying a layer in QGIS 2.0.1 Dufour, the Identify Results dialog shows (Actions) and (Derived) keywords: ![enter image description here](http://i.stack.imgur.com/u3AlB.png) Is it possible to remove these? Thanks


title: Suggestions on how to create an Owner Index for a Plat Book?

text: I've been tasked with creating a county wide plat book for publication this june. The cartography portion is fairly straightforward as I am using ESRI's Data Driven pages for ArcInfo 10 to produce the files by PLSS township/range (also my index layer) as pdf's than will be cleaning them up in Adobe Illustrator. I think I will find trouble when I need to create the owner index table for the back of the plat book as I am fairly incompetent with MS Access and was wondering if there is an efficient way to do this right from my .mxd as the information I need will be in the attribute table of my parcel layer. I need to figure out how to format the table as follows: 5 columns/page with 3 headers/column listing the information that will look something like this one column example of: Last/First Name, Section No., Pg No. Adams,Ryan.........14 25 Hood, Patterson....19 32 Tweedy, Jeff..........22 18 I found this post here that gives me some suggestions, the link to pyfpdf here being the one of most interest to me. Has anyone used pyfpdf to create an owner index, or any index for that matter, along the lines of what I spelled out?


title: Removing duplicates from large LAS files

text: Are there any publicly and free (not necessarily open-source) applications that can remove duplicate points from LiDAR clouds? I would like something on the lines of lasduplicate -unique_xyz, which seems to work fine, but LasTools are no longer free and have certain input limits above which the output is tainted by noise, missing information and it also seems to be incomplete. Also, I've tried SAGA GIS's "remove duplicate points", but this option is almost impossible to work with for large sets (more than a couple of million and SAGA GIS hangs). _LATER EDIT:_ a Python script able to filter out such an input is provided in my answer to this question. You should force a garbage collect from time to time if you run into memory problems.


title: Why is it important to ensure that the shapefiles have a spatial reference system when created?

text: Why is it important to ensure that shapefiles have a spatial reference system when created?


title: Union shapes which is not touching

text: I want to union two shapes which is not touching each other using Arcobjects. I do not want to perform small buffer on two shapes and union them . This is because when we apply buffer round buffer 0.002 it creates lot of additional vertices. This is very much custom requirement. Please help me how to proceed on this. Thanks in advance


title: queried layer is not updating after drawing new feature in original shapefile

text: I have a project in qgis 1.8 where I added the same shapefile 3 times. On 2 of them I made a query. The third one is for edditing purposes. When I change the value of one of the features (in the third shapefile) it is immediately rendered in the corresponding queried shapefile. But, when I draw a new feature, nothing changes at all. But when I close and reopen my project, the new feature is drawn in the corresponding queried layer. Is there a way to update a queried layer with new drawn features, without restarting your project? thanks


title: Creating a new map frame: Error HRESULT E_FAIL has been returned from a call to a COM component

text: I am trying to create a new map frame in the PageLayout Control. But I got the error: **Error HRESULT E_FAIL has been returned from a call to a COM component** at the method AddElement. Please see below codes:                                   IEnvelope envelope = new EnvelopeClass();                         envelope.PutCoords(Upper_LeftX, Upper_LeftY - _Height, Upper_LeftX + _Width,                             Upper_LeftY);                              IMap map = new MapClass();                              IMapFrame mapFrame = new MapFrameClass();                         mapFrame.Map = map;                         IElement mapElement = mapFrame as IElement;                         mapElement.Geometry = envelope;                              pageControl.ActiveView.GraphicsContainer.AddElement(mapElement  , 0); // throw error                         pageControl.ActiveView.FocusMap = mapFrame.Map;                         pageControl.ActiveView.PartialRefresh(esriViewDrawPhase.esriViewGraphics, null, null);      Do you know why ? Thanks and regards, Tai


title: Live update (the database & Map) possible in ArcGIS Server 10.1?

text: I am using sql server data base contains lat/lon values and using this information I have created the shape files and published on server.. When I get the update say 100 records then I need to recreate the shape file and merge it to (previous one) like this I am doing daily for last few days.. It is difficult to do daily for all records.. So is there any workaround?? I am using ArcGIS Desktop 10.1 and ArcGIS Server 10.1+ MS sql database 2008 r2 Any help will be great!! thanks in advance : )


title: Calculate Topographic wetness index in SAGA

text: i am currently working on a project which involves the quantification of the wetness in certain areas. However it produces quite different results from what i expected. I'll describe the steps i have taken. To do so i used the ASTER global DEMv2 and clipped it to the parts to be examined. I have to clip the raster, because using the original DEM results in various errors saying "memory allocation failed". ![ASTER DEM](http://i.stack.imgur.com/482Hz.jpg) Then i calculated the slope and Catchement area via the SAGA function Catchment Area (Mass-Flux-Method, parameter: flow-spit method area cell). I also tried the Parallel and recursive method, however the results are mostly the same. Then i calculated the Topographic Wetness Index, but the results look very mixed (somewhere between a rainbow and random fuzzy bits) and values are always negative. ![example of the resulting TWI](http://i.stack.imgur.com/aQv2z.png) I expected that the resulting raster should look like and display index values like the following image in the pdf PDF, page 14. Has someone here experience with the quantification of topographic wetness in an area and could provide a tutorial how to do it right? I would appreciate using SAGA or QGIS, but i also have ArcGis 9 at Hand. By the way: If i just use the calculate the SAGA wetness index with the clipped elevation data above, SAGA always keeps crashing (i believe because of the large no-data areas, because smaller smarts work)


title: How to autoclamp to ground automatically in Google Earth using a dynamic KML created from ArcGIS Server 10

text: It loads as an absolute everytime so terrain must be turned off to see the lines. Any idea how to get it so they default to clamped to ground? Using newest version of GE.


title: Reset local memory in ArcMap 10

text: I often am unable to delete files due to a lock or need to reset the data frame coordinate system. Very often starting a new session of Arc will allow certain geoprocessing or data management tasks to run effectively. What solutions exist to reset the local memory usage in ArcMap 10 other than closing and starting a new session? For the Python programmers out there, is there a way to add code to scripts that clears local memory before or after running geoprocessing tasks? Thanks for the advice.


title: Best method to aggregate 1m DEM mosaic to 10m DEM

text: I have a mosaic of a 1m LiDar derived DEM. I need to output a subset of the data as a 10m DEM. I am currently using the aggregate tool in ARCGIS 10 to produce a mean value for each new 10m pixel. Any advice on whether this is the best technique for such as task? Is the mean value the best approach with this type of data?


title: How can I transform layout page units to map units?

text: I am trying to build an add-in that will create a polygon from the current map extent. This is not a problem in Data View where I can easily get the envelope of an iMxDocument:iActiveView in data view. I need to do the same in Layout view, preferably without switching back and forth between layout and data views. When in Layout the iActiveView:Extent:Envelope returns an envelope coordinates of the map position with respect to the layout. (position of the map frame on paper so to speak) How do I transform these into actual map coordinates? I tried using iDisplayTransoformation:ToMapPoint but the input needs to be a screen point and the coordinates I have are of the map frame with respect to layout, not screen coordinates. In other words i need to transform page units to map units whereas iDisplayTransformation will transform only screen nits to map units or screen units to page units. Thanks


title: ArcGIS Search searching only in my own toolboxes

text: when I want to search for some tool and I write e.g. Merge into Search window (ctrl+F), ArcGIS (ArcMap) finds only tools named Merge or similar that are located in my own toolboxes. It doesn´t find normal ESRI installed Merge located somewhere in... management toolbox I think. Why? It just searches through my own toolboxes and not ESRI installed toolboxes. If I delete search index in options, it is ok then (it finds all Merge tools - default and also mine) but just for a while. Then something unnoticable happens and search is corrupted again. Does anybody know how I can fix this? Thanks jonlew Win8 64bit, ArcGIS for Desktop - Basic 10.1 SP1


title: Splitting polygon into different regions using ArcMap?

text: Is there any way that I can divide the below image into different areas? Just like a Google Map, where you can see the country borders. I am using ArcMap 10.0. ![enter image description here](http://i.stack.imgur.com/aqG9m.jpg)


title: WMS GetFeatureInfo and Google maps api V3

text: Is there a way to request GetFeatureInfo from a WMS service with the Google maps api v3?


title: update cursor to calculate field

text: This is my first post, so any help is greatly appreciated. I'm trying to use an update cursor to calculate a field. Basically, I'm trying to avoid division by 0 by looking in row[0] and see if it's 0. My code thus far is:               import arcpy     fc="C:/Temp/crimeGrid.shp"     fields=('MEAN_GRIDC', 'SUM_Join_C', 'crimeRate')     with arcpy.da.UpdateCursor(fc,fields) as cursor:         for row in cursor:             if (row[0] > 0):                 row[2] = row[1]/row[0]             else:                 row[2] = 0             cursor.updateRow(row)      Any ideas or help on what I'm missing? I can only assume (yes I know what assuming does) that it's the calculation aspect of my code. Thanks!!! **UPDATE:** Sorry for the ridiculous question. I continually thought it was in the division of rows, so I looked into my fields a little closer. Turns out, Row[0] and Row[1] were Double and Row[2] was Long. Fixed that aspect and the code works as is.


title: GDAL Python bindings not picking up OGR PostgreSQL driver

text: I have a windows server (Microsoft Windows Server 2003) with:   * Python 2.7,    * PostgreSQL 9.04-1,    * GDAL 1.8.1 (installed from the windows installer (gdal-18-1500-core.msi) hosted at http://vbkto.dyndns.org/sdk/PackageList.aspx?file=release-1500-gdal-1-8-mapserver-6-0.zip), and   * GDAL Python bindings installed from the windows installer (GDAL-1.8.1.win32-py2.7.msi) at the same url as above. If I do               ogrinfo PG:"dbname=postgis_tpl user=postgres password=xxxxxx"      connection with the database is successful:               INFO: Open of `PG:dbname=postgis_tpl user=postgres password=xxxxxx'       using driver `PostgreSQL' successful.     1: world_merc (Multi Polygon)      However, from the Python shell, if I do:               >>> from osgeo import ogr     >>> print ogr.Open("PG:dbname=postgis_tpl user=postgres password=xxxxxx")      I get:               None      and:               >>> print ogr.GetDriverByName("postgresql")      gives me:               None      and if I list all the drivers available via the Python bindings:               >>> for i in range(0,ogr.GetDriverCount()):     ...     print ogr.GetDriver(i).name      the following drivers are listed (PostgreSQL is missing):               ESRI Shapefile, MapInfo File, UK .NTF, SDTS, TIGER, S57, DGN, VRT, REC, Memory, BNA, CSV, NAS, GML, GPX, KML, GeoJSON, GMT, SQLite, ODBC, PGeo, MSSQLSpatial, PCIDSK, XPlane, AVCBin, AVCE00, DXF, Geoconcept, GeoRSS, GPSTrackMaker, VFK, PGDump, GPSBabel, SUA, OpenAir, PDS, WFS, HTF, AeronavFAA      My environment variables are (I think) set correctly:               GDAL_DRIVER_PATH = C:\Program Files\GDAL\gdalplugins      (There's nothing that looks like the PostgreSQL driver in the gdal-plugins folder anyway).               GDAL_DATA = C:\Program Files\GDAL\gdal-data      and               C:\Program Files\GDAL      is in my PATH (GDAL is installed at C:\Program Files\GDAL). So what am I missing? Why can't Python see the PostgreSQL driver? Do I need to install its bindings separately?


title: Application of GIS in passive and agent defence

text: I am working on passive defense modeling, using GIS. Modeling passive defense can be better implemented using GIS and I've come up with my own methods. It would be far better for me if I know there are any previous works on this area. I want to be informed about any application of GIS in passive and agent defense that are publicly available. For example, where to locate the military bases, so that we have the minimum casualties in the case of external attack.


title: Difference between regionalization and image segmentation

text: Regionalization is a special form of clustering where spatial units are grouped together, based on a set of defined criteria (cluster homogeneity) and a set of contiguity or adjacency constraints. Isn't this exactly the same thing that image segmentation tries to accomplish: finding contiguous and homogeneous clusters in image data? What is the difference?


title: Sextante Plugin error in QGis

text: I am an user of QGis 1.9.0. When I go to Fetch plugin option Sextante plugin is broken seems what is the problem?Even I agin downloaded the plugin in Qgis official resportary but it displays me the like below:![enter image description here](http://i.stack.imgur.com/cbIRZ.jpg)


title: Changing CSV layer attribute value from text to number in QGIS?

text: I read in a CSV file for linking to a .shp file. All the variables that are just numbers came in as text. How do I change the text variables to numeric in QGIS?


title: How to prevent QGis (MMQGIS plugin) from interpreting numbers as strings

text: I'm using the MMQGis plugin to "Join Attributes From CSV" into a shape file and am having trouble getting QGis to interpret some numbers as numbers (rather than strings which is what its doing). This is becoming a real problem because the values do not behave like they should when trying to order them, select above a range, etc. I've seen two somewhat similar questions to this but the answer to this question instructs the user to import the csv as a layer. But the CSV data I'm trying to import is not geocoded and I'm using values (things like zip code, etc) to join in the CSVs, so I don't think that question's answers might help here. Thanks


title: String to Integer in QGIS

text: I have a csv (below) of values I would like to represent in a choropleth map. I've managed to join them to a shapefie via fips code successfully in this format, but the last four columns do not show up under the option for graduated symbolization. (I also changed them from percentage to decimals in excel first) In the QGIS Properties field tab, I notice these were all in string. How do I convert these attribute columns to real or integer form so I can map them? Do I do this in QGIS or Excel? http://www.nytimes.com/packages/text/newsgraphics/2009/1128-foodstamps- data.csv


title: Why is QGIS reading in my CSV numeric data as string data, and how do I change it?

text: QGIS keeps reading in my numeric data from .csv files as string data. This means I have to figure out how to read it in with the proper data type, or convert it back to numeric inside of QGIS, neither of which I know how to do. Why is it doing this, and how do I fix it?


title: Trying to run raster calculator to build an expression for ls factor

text: I am trying to run a raster calculator, but it failed. I believe that ArcGIS 10 does not understand when I use the variable resolution. What do I need to put in its place? I took the equation from this sheet USLE from North Carolina State University and checked on 2.3 for ArcGIS 8.1. Here is the code I used: Power (("flowaccNW") * resolution/22.1, 0.6) * Power(Sin("slope_degnw") * 0.01745) / 0.09, 1.3 I am using ArcGIS 10 SP4 using the Map Algebra raster calculator. Thanks


title: Working with WKT data in XML File in SQL Server

text: Anyone have experience with importing XML files into SQL server 2008 or 2012 that have a WKT element that includes the geometry for each feature? I can get this to easily work using FME software but wondering if this can't be done directly using SQL Server. In FME I used the GeometryReplacer transformer. These XML files have all of the WKT data in a single XML element for each feature.


title: ArcSDE 10 FeatureLayer is read only, FeatureTypes property is null, not sure why

text: I have a WPF FeatureLayer backed by an ArcSDE feature class that I hope to use for editing. After the layer is done initializing its readOnly property is true, I think that this points to some problem with the ArcSDE setup but I can't figure out what it could be. I've verified that the connection is working, inspecting the FeatureLayer shows that it has downloaded the remote table name. The post installation setup for ArcSDE had no trouble, same with creating some empty feature classes, adding them to an mxd file and using it to host a feature service. Is there any other setup I may have missed? Thanks!


title: Can't get correct number of populated places for a place

text: I am trying to get a number of populated places for example for Andorra since it is a small number. Here is a geonames site with 10 places: http://www.geonames.org/AD/largest-cities-in-andorra.html But if I do this query: http://api.geonames.org/search?featureCode=PPL&username=genadinik&country=AD&style=full&maxRows=100 Only 4 out of the results have population. Am I doing something wrong? Is there a more accurate way of getting populated cities there? Thanks!


title: In the Symbology Properties of an ArcGIS layer, why is the count sometimes not shown?

text: Sometimes the count column does populate for me, but most times it does not seem to. I haven't been able to figure out a rhyme or reason to it. There is no definition query or anything out of the ordinary. Is this a bug or is there something you can do to make sure that field is populated? I am running ArcView 10 on Windows 7 x64. ![screencapture of symbology window](http://i.stack.imgur.com/Se5qv.jpg)


title: How to join fields which must still be created in ArcGIS Model Builder?

text: I am in the process of creating a **generic** model in ArcGIS Model Builder. At one point in my model, I have to do a join between an existing layer (entry parameter) and a table issued from the "Zonal Statistics as table" tool. I get an error when I launch the model from the Toolbox, as the field to do the join doesn't exist. When I create a "Zone field", "Input join field" and "Output join fields" as parameters, I can get the fields to be defined by the user. However, he has to write the name of the "Zone field". Is it possible to select that value in a list (like the fields from my existing layer)? If I am doing things wrong, how can I get the model to say which fields to do the join on? Regards, PS: One layer is a .shp that is defined as a parameter. The other is a layer produced by a raster calculator in the geodatabase then reclassed and to which a "Zonal Statistics as table" was applied. Version 1 that returns an error: ![Here is a view](http://i.stack.imgur.com/slQ3M.jpg) Version 2 that works, but doesn`t enable the user to pick the field the join needs to be on in a list: ![enter image description here](http://i.stack.imgur.com/O35RZ.jpg)


title: What is the default unit for QGIS scales and how do I set it to KM?

text: I want to render a specific layer only between 0-12km but have no idea how to set it in QGIS because the scales dimensions are high numbers. My project and layers CRS is EPSG:2039 Thank you .


title: Can Compound Topographic Index (CTI) / Topographic Wetness Index (TWI) produce negative values?

text: I am calculating Topographic Wetness Index for a study area using SRTM DEM (90 and 30m). I've seen several versions of TWI (CTI). I have been using: TWI = ln(As +1 / tanB) where As = catchment area, B = slope (degrees/radians??) When I run this function, the output contains negative values (using slope in degrees and radians). Are negative values acceptable? Or does this indicate an error in the analysis?? I understand the output should be scaled to 1-10, but negative values seems counter intuitive. Could someone explain what negative values could mean (due to flat areas, no data values?) On a side not, the calculations are being run on on a very large area, which overaps several watersheds. Is it appropriate to run TWI over such a large extent?


title: Detect connections along shapefile in QGIS?

text: I have a set of polylines and a set of points that lie along or very near to the polylines, and I need a way to measure which points are connected to each other along the polylines, and what the distances between each point are. Having the results in a neat grid would be nice, but I'll take what I can get. I have the vague idea that I need to use GRASS or SEXTANTE, but I don't know where to go from there (then again, I didn't know what GIS was a few months ago). Also, I have access to ArcGIS, but I'd rather use QGIS if possible. To clarify the data a bit, my polylines are waterways in the Talladega National Forest, and the points are beaver ponds that myself and my research professor sampled. The points are either on or very close to a waterway, close enough that Road Graph can use them. I also have an area estimate polygon for each point, if direct connection is needed. I know that some of the points lie on different waterways from the others, but I need an accurate and repeatable way to find those connections.


title: how to obtain elevation points for a road using contour?

text: I have contour lines, & roads centerlines feature classes, I need to get points with elevation in the starts, intersections, & the ends of the roads, how can I get this information? I used the "Intersect" approximity analysis tool, but I think it is not accurate enough. is there a way in 3D analyst or Spatial analyst? Thanks In the case that roads overlap each other on different elevation, is there better way to still find the elevation for each street centerline segment (from node, and to node)? I have both contour line and DEM for the area. But the roads overlap a lot. Thanks


title: What spatial analysis did the US Supreme Court review in the Perry v. Perez opinion?

text: The Supreme Court recently issued an opinion on Perry v. Perez which seems like it must have involved evidence based on spatial analysis. The L.A. Times reports: > The 9-0 decision set aside a new map of congressional districts drawn by a > special federal court in San Antonio that gave Latinos and Democrats a good > chance to win three or possibly four new seats in the House of > Representatives. Does anyone know what type of spatial analysis - if any - was used to create the evidence presented? **Update:** Thanks to mkennedy for interpreting this, and also for pointing me in the direction of _preclearance_. I see that the Federal Register Vol 76 No. 27 Part III goes into what I was wondering about under the " _Retrogressive_ Effects". Fortunately, this is very googleable.


title: How to show "First Vertical Derivative of the Magnetic Field" survey data on a map?

text: I will be given a magnetic survey raster and I am being asked to display this on a map as a "First Vertical Derivative of the Magnetic Field". I've done this before with a different older dataset but I basically tried to match the raster colour ramp to a an example of a hard copy map. I am not a geoscientist or geologist but **I would like to understand (in a nutshell) what "First Vertical Derivative of the Magnetic Field" is and how to correctly symbolize it on a map.** I've seen some examples on the web and it seems fairly straight forward but, is there a predefined universal way of displaying this type of data?


title: What I'm doing wrong with the command r.out.ascii in GRASS GIS?

text: I have raster in .bil format which I want to save in GRASS GIS with r.out.ascii. I get this result in .txt file:               north: 1N     south: 0     east: 1E     west: 0     rows: 1     cols: 1     55537       I suppose it is something wrong, because of dimensions of this .txt file. I didn't change any options in r.out.ascii window. Thank you for your help.


title: Monte Carlo simulation using QGIS and pgRouting on optimal sidewalk construction

text: I am new to spatial analysis and would appreciate some general direction on a project I am attempting, outlined below (I am starting from scratch). **GOAL:** To find the _best_ locations to install 2000 feet of sidewalks in my hometown in order to connect the most households to the Central Business District (CBD), where "connected" means within 1.2 miles walking of the CBD. I have shapefiles showing existing structures (households), roads, and sidewalks (already installed). Here is my proposed solution / thought-process:   1. Convert the in-place sidewalk network into a database of nodes that are connected by weights (i.e. distances). Is there a way to directly do this in QGIS (or other program) by clicking on all intersections?   2. Calculate the number of households that are within 1.2 miles walking of the Central Business District (e.g. a lat-long point or polygon) using the routing capabilities of pgRouting or something else. This will be the base case "household access" value.   3. Using the road layer as a guide, randomly place an additional 2000 feet (say, in 10 foot segments) of sidewalks onto the sidewalk layer. This is the equivalent of constructing a bunch of new sidewalks arbitrarily.   4. Re-calculate the nodes and weights using the new pedestrian network as in (1), and then re-calculate the number of households that are now within 1.2 miles of the CBD as in (2). It should increase with the additional sidewalks. Save the locations of the additional sidewalks and the associated "household access" value to a file (e.g. spreadsheet).   5. Repeat steps (3) and (4) 10000 times, similar to a Monte Carlo simulation. Using the 10000 sets of data points, choose the sidewalk placement locations that maximize the number of households within 1.2 miles of the CBD. * * * **Does this thought process sound realistic? Does anyone have any suggestions?** -- I would like to accomplish this using some combination of QGIS and R, however I am open to learning PostGIS and/or Python (or anything else) to achieve the goal. Any thoughts or direction will be helpful. Many thanks!


title: How to programatically generate my own wms file?

text: I'm new to OpenLayers and as far as I understood layers are retrieved from an wms file. I know this is not precisely a development forum but I was just wondering if there is a way to fool this framework and make it read a dynamically generated wms with a programming language. I mean, instead of giving it the URL of the wms, make it load a PHP with certain parameters and make the map display the data I want. Would it be the same? Is there a way to build a wms file? Where can I find a good tutorial to do this? My project involves a map like this one: http://maps.nrel.gov/biomass. Look how different colors are painted over the map depending on the product I choose on the left. This is exactly what I'm trying to achieve. I have the data in a relational database, not in a wms that's why I'm asking if there's a way I can "make" my own wms or make OpenLayers framework read my data from it to display it on the map. Thanks


title: Accessing non-spatial PostgreSQL tables using OGR Python bindings

text: When connecting to PostgreSQL/PostGIS using the OGR Python bindings is it possible to get non-spatial table names (i.e. OGR layers) from the connection? Currently, doing:               conn = ogr.Open("PG: #params")     for layer in conn:       print layer.GetName()      Will only print the names of spatial tables, however I can access non-spatial tables by (for example) directly specifying the table name, e.g.:               layer = conn.GetLayerByName('non_spatial_table_name')      The OGR docs refer to the "PG_LIST_ALL_TABLES=YES" parameter, but I can't find how this links to API (let alone the Python SWIG bindings!). Any pointers appreciated. Thanks, Tom


title: How to extract data for a specific area from OSM?

text: I need to extract all streets available on a particular area of a City. I'm able to extract all streets name using this query:               select distinct P2.*     from      planet_osm_polygon P1,     planet_osm_line P2     where P2.name is not null AND P1.osm_id = -80277 and ST_Intersects(P2.way, P1.way)      order by P2.name      (in this casa -324211 is the city of Toronto). Now, I need to extract the streets name of only a particular subarea of the City. I can found the information about the subarea on planet_osm_point table but, how can I extract only the line or rows are in one particular subarea? With this I can extract all subarea of Toronto:               select  P2.osm_id, P2.boundary, P2.name     from      planet_osm_polygon P1,     planet_osm_point P2     where P2.name is not null AND P1.osm_id = -324211 and ST_Intersects(P2.way, P1.way)  and P2.place = 'suburb'     order by P2.name      is there a way to do that? --- UPDATE --- One solution is, starting from a point, create a square or round polygon and extracts all streets that intersect this new polygon, something like that:               select asText(way) from planet_osm_point where osm_id = 349354234      and, when I use the results here:               select  P2.*     from      planet_osm_line P2     where P2.name is not null AND ST_Intersects(P2.way, (ST_Buffer('POINT(-8840304.91218709 5429889.22642508)', 50) ) )      I get this error message: Operation on mixed SRID geometries On a blog I read how to update my SRID table geometry column but, is it the right way to fix that? thanks for your help! Any suggestion (not only a solution) will be accepted and very helpful! Andrea PS: I imported data using osm2pgsql but, if I understand well (I'm newbie here) I can found more information using osmosis, correct?


title: How to scrape addresses from a webpage for geocoding purposes?

text: I am trying to get addresses from a webpage for mapping purposes. I know I could copy and paste the addresses into an Excel file and geocode them but I was wondering if there was any faster way to get the address from a webpage and create a point location file. Thanks in advance.


title: How to group connected linestrings in PostGIS?

text: I have a table of streets I've selected based on a set of attributes (let's say it's `speed_limit < 25`). There are groups of streets that are locally contiguous; I'd like to group these sets of connected linestrings into GeometryCollections. In the image below, there would be two GeometryCollections: one with the red lines and one with the blue lines. ![enter image description here](http://i.stack.imgur.com/WNlxX.png) I tried running a couple of "dissolve, deaggregate" queries along the lines of:               SELECT (ST_Dump(st_union)).geom     FROM          (SELECT ST_Union(geom) FROM roads) sq      With everything I've tried, I either end up with a single feature (`ST_Union`) or my original geometry (`ST_Dump` of `ST_Union`). Maybe it's possible to do this with some kind of `WITH RECURSIVE` magic?


title: How to populate a text field with double quotes using field calculator

text: How would I populate a text field with double quotes using the Field Calculator. For example, I want to populate a field as 24". Is there a way to populate fields using the field calculator with this value? If I use "24"", I get an error because of the two double quotes. Thanks in advance!


title: Consume geojson web service in ArcMap

text: I need to consume data from a REST endpoint that returns features in geojson format within a model in ArcGIS 10.2. Is there a way to access a URL instead of opening a file? If there is no way to access a URL directly, would this require a custom python script to access the URL and save it to a local file, or is there a better way?


What defines a Class in the Heat map plugin? As in: "Black holes like that are usually caused by raster values exceeding the range of the color ramp. Try adding another **class** with a higher associated value. That should fix it." This quote comes from an earlier anwered question: Why is there a black hole in the Heatmap? I'm getting all sorts of clunky results and I'm not sure whether I need to find the right settings, alter my data, or both. I'm working on identifying all important parameters: obviously Radius, and cell size, their relation, extents (I'm foggy on), projection (I'm still not positive how it affects the heat map as long as it sit on top of the data. And it does.), projection units,... I'd like to know better specifically what affects the curve of the histogram. Also is the **color ramp** synonymous with the histogram curve? In addition to lots of other issues, I too have black holes. But I have to say it is awesome fun ferreting all this out. And Thank you, Underdark for pointing the way. I'm starting to find some light with the raster histogram.


title: Method for finding circle center points that cover an extent

text: If I have a geographic extent, is there any way to find all of the center points of all of the 5km circles that would be needed to cover the entire extent? I'm working with an API that will return results that occur within a 5km radius of a given latitude and longitude. I would like to come up with a way to get all of the items occurring within the extent.


title: Simplest way to share map and data with ArcGIS 10.x and 9.x users

text: What is the simplest way to share a map and data with both ArcGIS 10.0 and 9.3 users, portably? This question could be rephrased as _"Instead of waiting for Esri to (maybe)implement backwards compatible map packages, what can we do?"_ Assume that one is starting with project resembling the following, all in v10:               D:\workspace:         MyCoolMap.mxd:             default.gdb:                 Transportation\roads                 Transportation\trails                 Hydro\Lakes                 Hydro\Rivers                 some_random_points                 a_linked_table           a_shapefile.shp           a_linked_table_for_fun.csv           {some_coded_value_domains}      I've tried saving the mxd as 9.3, and then create a map package from that, but get an _"error compressing map package"_. Using the Feature Class To Geodatabase (Multiple) tool is unsatisfactory because it can't handle Feature Datasets. Another failed approach was to export the v10.gdb to binary xml and import it into a v93.gdb, _("Import xml data failed. The domain field type does not match that of the field it is being assigned to [feature class name]")_ Got any other approaches that might bear fruit?


title: How to create a ClassBreakRenderer without data coming from Feature Layer?

text: I would like to create a ClassBreakRenderer to represent integer values. The symbols I am using are defined as PictureMarkerSymbols with JSON and they are just a filled blue circle with different sizes that were created with ArcMap 10.1 and published into ArcGIS Server. This is just to introduce the thing. The problem is that in all the examples I see on the ArcGIS Javascript API about Rendering, the examples **do always** get the data from a **particular column** of a Feature Layer, but **my data do NOT come** from this data source: **they come from an array of integers in client side.** Just the symbol comes from the server-side, but not the integer value. **Is there any way of "tricking" the constructor** of ClassBreakRenderer **in order to accept an array** _(or collection, dictionary, set, whatever)_ of integers to be classified in symbols? There **should be "life"** a part of Feature Layers, isn't it? Thanks for your help! * * * **EDIT** The thing is that **my data source is an array** of integers like this:               var myArray = [15, 37, 47, 12, 125, 87, 36, 66, 10, 4, 18, 19, 23]      AND NOT:               var featureLayer = new esri.layers.FeatureLayer("http://sampleserver1.arcgisonline.com/ArcGIS/rest/services/Demographics/ESRI_Census_USA/MapServer/3", {               mode: esri.layers.FeatureLayer.MODE_SNAPSHOT,               outFields: ["*"],      And in principle, the ClassBreakRenderer receives a column of that Feature Layer as:               var renderer = new esri.renderer.ClassBreaksRenderer(symbol, "POP07_SQMI");      But **I do NOT have that column "POP07_SQMI" I have an array called "myArray"** and my question is that if it is possible to use the array instead of the column name that I don't have. Because world should not be limited to the existence of feature layers to retrieve data from them.... :S


title: Are there any tricks to using a schema other than public in PostGIS?

text: I'm currently setting up a fresh install of PostGIS 2.0.2 and PostgreSQL 9.1.6 on Ubuntu. I've recently come across some information indicating that using the public schema to store all data is not a good idea.   For this reason, I've set up a schema called data and made myself the owner. Here are my questions:   1. Besides setting the owner, do I need to pay attention to anything on the Privileges tab when creating this new schema (through pgAdmin III)?   2. Could I get the same benefits by storing my data in the public schema and dumping all data into a separate schema before doing a backup/restore (this would save a few keystrokes when using ogr2ogr)?   3. Will I run into any trouble by not having the default PostGIS tables and views in my new data schema (they are in the public schema within the same database)? Thank you!


title: Does tile server renders every request?

text: I am trying to understand how tile server handle requests. I red about Mapnik which "creates" the map using OSM, Does "creates" means that every request is rendered and returned? If yes what is the performance hit (global map)? Or it means that Mapnik Creates all tiles needed on initialization and then serves the right one upon request? if yes what is the total storage needed for all tiles(for global map)? I am trying to find a way to build my own tile server (global map) while rendering it only once.


title: Grass won't import shapefiles

text: I have a problem involving the import vector data command, specifically importing shape files. I installed QGIS on a pc at work and there I could import shape files into GRASS GIS 6.4.3RC2 (which is installed along with QGIS) and then export them as .DXF but it took a lot of time to do that so I wanted to finish on my pc at home which is much better than those at the office. I downloaded QGIS again at home and made the same steps as at work but now GRASS wont import the same shapefiles... This is what the program tells me:               Traceback (most recent call last):       File "C:\PROGRA~2\QUANTU~1\apps\Python27\lib\site-     packages\wx-2.8-msw-unicode\wx\lib\filebrowsebutton.py",     line 140, in OnChanged          self.changeCallback(evt)       File "C:\PROGRA~2\QUANTU~1\apps\grass\grass-6.4.3RC2\etc\w     xpython\gui_core\gselect.py", line 1601, in OnSetDsn          self._reloadLayers()       File "C:\PROGRA~2\QUANTU~1\apps\grass\grass-6.4.3RC2\etc\w     xpython\gui_core\gselect.py", line 1620, in _reloadLayers          dsn = dsn)       File "C:\PROGRA~2\QUANTU~1\apps\grass\grass-6.4.3RC2\etc\w     xpython\core\gcmd.py", line 660, in RunCommand          ps = grass.start_command(GetRealCmd(prog), flags, overwrite,     quiet, verbose, **kwargs)       File "C:\PROGRA~2\QUANTU~1\apps\grass\grass-6.4.3RC2\etc\p     ython\grass\script\core.py", line 177, in start_command          return Popen(args, **popts)       File "C:\PROGRA~2\QUANTU~1\apps\grass\grass-6.4.3RC2\etc\p     ython\grass\script\core.py", line 56, in __init__          startupinfo, creationflags)       File     "C:\PROGRA~2\QUANTU~1\apps\Python27\lib\subprocess.py", line     679, in __init__          errread, errwrite)       File     "C:\PROGRA~2\QUANTU~1\apps\Python27\lib\subprocess.py", line     867, in _execute_child          args = '{} /c "{}"'.format (comspec, args)     UnicodeEncodeError     :     'ascii' codec can't encode character u'\xfc' in position 43:     ordinal not in range(128)      can anyone help me here?! Thank you for your help


title: Field within a Field - Donut Parcel?

text: Firstly, apologies as I'm a newbie and have only been studying GIS for a few months. I am digitising field boundaries and have a field within a field, and am wondering what the best way of representing this in a single layer is, so that the area of the internal field is subtracted from the area of the main field. I have read an explanation about donut parcels but am not sure if this is the right solution for what I want to do? I am using ArcMap. Thanks in advance for your suggestions.


title: Understanding the QGIS API documentation to write python code

text: I am not sure how to use the API documentation to code in Python. For example, I want to determine whether a layer is a point, line or polygon layer and found the following in the API documentation: ![enter image description here](http://i.stack.imgur.com/qjx5J.png) I have selected the lakes shapefile from the QGIS alaska dataset. In the QGis Python console, I have typed:               from PyQt4.QtCore import *     from PyQt4.QtGui import *     from qgis.core import *     canvas = self.iface.mapCanvas()     cLayer = canvas.currentLayer()     type(cLayer)      The returned value is: `<class 'qgis.core.QgsVectorLayer'>`. How do I use the above information in the API documentation to (1) determine the type of the selected layer and (2) return the value of the listed enumerate types (e.g. polygon)? Whilst I am sure someone can tell me the answer, I would really like to know how to work out the answer using the API documentation and any other available documentation.


title: Count number of attributes for unique field values

text: Running ArcGIS 10.1, I have two feature classes in a geodatabase. One is a point FC for polling stations, the other is a polygon FC for states. The polling station FC contains over 1000 points and has a states field indicating which state the polling station is located in. The states FC has an attribute for each state and has a blank field for the number of polling stations in each state. I need a python script that will get the count of how many polling stations are located in each state from the polling station FC, and then write those numbers to the column in the states FC - obviously to their corresponding state. I've tried a few times to start this, but I can't seem to figure out the best way to do this. Any help is much appreciated!


title: Why can't I save csv as a shapefile

text: I am trying to follow the instructions given here, Tutorial: Making Heatmaps using QGIS and GRASS , to build a heat map with my own data. I taken the first step and imported a csv data table into the layer bar, but when I attempt to save it as a shapefile, QGIS tells me that is not a valid or recognized data source. What am I doing wrong?


title: Identify polygon overlaps within single shapefile

text: Using either QGIS (1.9) or ArcInfo (9.3.1) I'm trying to identify unique overlaps of prehistoric building footprints from an archaeological site recorded within a single polygon shapefile. While the Overlay>Interset and Geoprocessing>Intersect tools in ArcInfo and QGIS, respectively, work to some degree both return multiple results for each unique overlap. An appropriate parallel for my problem I suppose can be illustrated by a simple 3 unit Venn diagramme in which there are 3 unique intersections or overlaps, yet when running an Intersect routine (particularly that in ArcInfo) I often get up to 12 unique polygons returned - 2 for each shared 'leaf-shaped' intersection, and a further 6 for the centre overlap (each intersection being counted twice). I'm trying to produce a shapefile illustrating those unique intersections. Ideally I would love for the attribute table of this produced shp to contain details of each intersection i.e. Feature_X intersects with Feature_Y, but I fear I may be asking for the moon on a stick. Many thanks in advance for any and all advice, Marc Edit: Following the excellent suggestion of @PolyGeo, I have attached a rough sketch below to hopefully better illustrate my current problems and what I'm hoping to achieve. Forgive the hand-drawn sketch, I'm away from my main PC. ![A hand-drawn sketch of my current problem re: intersecting polygons](http://i.stack.imgur.com/oOAav.jpg)


title: How to convert a Visio diagram to a GIS layer

text: I have a office floorplan that is a standard Visio drawing file (*.VSD) that I would like to convert into a GIS feature class (shapefile, GDB, PostGIS, or any other applicable format, it doesn't matter much). A couple of options occurred to me, but I am wondering if anyone know of any other tools/techniques for converting a Visio file to a feature layer? My ideas:   * Simply load the VSD file into a Desktop GIS software and heads-up digitize the features (walls, doors, cubical's, network ports, etc)   * Convert the Visio to Raster, then using ArcScan to help automate the generation of feature linework. Another idea I found in the question: Transforming a polygon to ESRI shapefile? was to:   * "export Visio Shapes (file, save as) to a .dxf (AutoCAD exchange file)", which seems like it might be the easiest method I have found so far. What additional tools, techniques, ideas, or suggestions might I be overlooking? My preference would be to use ArcGIS 10.0 tools to accomplish this task, but am open to other possibilities.