Visualization of a big graph in 3D from postgresql/postgis

We have a graph in 3D in a PostgreSQL/PostGIS database. It contains something like 50'000 edges and 50'000 vertices, with x-y coordinates and floors. We would like to visualize this graph with straight lines as edges in 3D and to add on top of this few more time-dependent data related to the graph. What would be the best technology to do this? Anyone with similar experience? We would prefer this technology to be:   * open source   * as easy as possible   * but still able to manage this graph and time-dependent data Would Processing be able to manage this? Or would it be better to go for OpenGL? Is there any other technology that would be more appropriate for this? And how long would it take, if you've any experience, to visualize the graph with each technology? I've hear that processing already reduces the frame rate with 10'000 links. Also, it seems interesting to use a bit of both. OpenGL seems nice since you can program it in any language (c, c++, objective-c, c#, java, etc.) and results could already be obtained in one week.


Problems with Combining Layers in QGIS

I am trying to add the political boundaries baselayer from http://www.cec.org/Page.asp?PageID=924&ContentID=2336 and then on top add points where I gathered coordinates from Google Earth (using a delimited text layer). I can see both layers in QGIS but my points are not in the right location and I have to zoom in way too close to be able to see them. It seems to be an issue of scale or projection or something. Im new to QGIS and GIS in general, what things can I try to try and make these layers fit better.


How to handle add feature event in the FeatureLayer/Editor

I use a FeatureLayer created from the featureCollection. I also use esri.dijit.editing.Editor to add a feature to the map. Every time feature has been added to this layer I want to send it to my remote service using rest API. I tried to use FeatureLayer's onEditsComplete event:               <!doctype html>                     <title>...</title>         <link rel="stylesheet" type="text/css" href="http://serverapi.arcgisonline.com/jsapi/arcgis/2.2/js/dojo/dijit/themes/claro/claro.css">

         <button dojoType="dijit.form.Button" onclick="initEditor();">             Set up editor         </button>                  But it fires before user has entered information in the dialog created when user clicks on the map. So it seems useless. I tried onAddAttachmentComplete event but it never fires. Do you know how to specify that mentioned action shall be done after user entered text in the editor's dialog?


How can I switch line direction in QGIS?

Something that is hopefully pretty easy. I would like to be able to swap the direction of a polyline/line in QGIS. I made a custom tool to do this in MapInfo a couple of years ago, however I can't seem to find anything for QGIS. Does anyone know of a tool to do this? _If one doesn't exist then I'm happy to have a go at making it, just didn't want to if one has already been made._


Loading Shape file in flexviewer 2.2

I want to load Shape file in Flexviewer 2.2 Any sample would help a lot. Thanks


Automatic updates to OpenLayers code

I have my map code in openlayers JS. My input data is in a text file and I am able to render the points on the map. I have to give the user a right to create a layer by entering the data in a form. Writing a program to insert data into a text file is trivial. But I am stuck at how to include this layers code in the JS programmatically so that after giving the layers details in the form if the user reloads the page then the new layer should be visible, Any solutions??


How to clip a layer file

I have a world imagery layer that I want to clip to a single county. I would just save as a JPEG and edit using other software, but this uses a lot of processing power if I want a high resolution. I would like to clip the image before exporting to JPEG. Is this possible?


How to reproject to other projections than WGS84 or Webmercator?

I have an OpenLayers map with a vector layer to add points on. I'm using Google Maps as the base layer. After the point is added, I need to store the data in my database as lat/long and x/y. Can I transform the projection of the x/y points or convert the lat/long to x/y? The x/y that is on the point is in spherical mercator so it doesn't work for me. Update: I'm not having any trouble converting the x/y to lat/long. I want to convert the x/y that I have to a different projection. I'm not sure that's the right terminology, but I need to convert the x/y values. Update: The SRID I need is 3734. When I drop the point, the x/y is like: x: -9058948.6891836 y: 5035248.6338434 But what I need looks more like: x: 2282175.42 y: 545126.58 I'm also using PostGIS so if this is possible to do with PostGIS I'm open to that also.


Geometry Calculation

Wanted to ask about accuracy of areas calculated by GIS compared to actual ones on ground. Because I have data of my country without areas for the districts within Uganda. Are the areas obtained using GIS accurate and do they reflect the actual areas on ground or in the field?. Thank You. Senoga Moses


FME Feature Type with Python

Here, the python for create a list of feature type but....               import sys     from pyfme import *     try:         log=FMELogfile("shapeToTAB.log")              log.log("Opening reader",FME_INFORM)         reader=FMEReader("GEODATABASE_FILE")         reader.open("C:\depot\export\TE-PMB.gdb")              GeomList = []              schemaFeature=FMEFeature()         log.log("Copying schema features")         while reader.readSchema(schemaFeature):             GeomList.append(schemaFeature.getFeatureType())             log.logFeature(schemaFeature)              reader.close()         print GeomList          except FMEException, err:         print "FMEException: %s" % err         sys.exit(1)      I would like **return** a list of feature type of my geodatabase. Is it possible ? Thanks


Interpreting Metes and Bound description?

![enter image description here](http://i.stack.imgur.com/FjMmR.jpg) This is the Metes/Bounds description - **5.3 acre parcel in the NW4, commencing at the Northest corner, thence South 789 ft to the place of beginning running thence: East 30 feet to the East side of the present County Road; thence S. 17 degrees - 05' E. 74 ft; thence S. 26-30 degrees E. 500 ft; thence S 1-40 degrees E 382 ft; thence S. 46 degrees W. 356 ft; thence West 30 feet to the west line of the section; thence North 1148 ft to the place of beginning.** * * * My question that I have a problem with is """thence South 789 ft to the place of beginning running thence: East 30 feet to the East side of the present County Road""" The parcel that I did with the metes/bound is in Yellow, slightly below the Borie text. The parcel is supposed to be in the white/green box to its left. I started my m/b off at "S 17-05 e / 74" and proceeded from there. It is supposed to be 5.3ac, but I am getting 4.7ac afterwards. when I try putting in W / 30ft and then N and 1148 ft via metes/bounds, there is no movemenet. What is it I should be doing or not doing? The ArcGIS 9.3desktop is not very helpful. How does the "thence South 789 ft to the place of beginning running thence: East 30 feet to the East side of the present County Road;" AND "thence West 30 feet to the west line of the section; thence North 1148 ft to the place of beginning" come into play? I'm using 9.3 ArcView if that helps.


loading CanVec directory full of .shp files

I'm totally new to QGIS, and trying to view a Canadian Topo map **as published** on an iMac - I don't want to change it. http://www.geogratis.ca/geogratis/en/option/select.do?id=5460AA9D-54CD-8349-C95E-1A4D03172FDF&datasets=031g05   1. the directory contains a lot of files - can all be loaded in one go so I get the whole map, not just pieces of it?   2. the CanVec attributes (colour, symbol etc.) aren't retrieved when loading a shape file with vector load, so that doesn't work to see it the way the map maker intended. I suspect these problems are related, in that there's a file somewhere that lets QGIS tie everything together, but I can't figure out what it is or how to tell QGIS how to use it.


Cannot See Georeferencer Results

I am working in QGIS 1.8.0 Lisbosa. Currently I am trying to georeference a building blueprint in a .TIF format to a Google Streets layer. After selecting four points (the corners of the building), I run the function (Thin Plate Spline - Cubic Spline) with a target SRS of UTM Zone 18N (EPSG: 3725). However, the results don't show up on the map. When I click zoom to layer, the screenview doesn't move. So either it isn't working, or the results are invisible. Any idea what I'm doing wrong here? I know I've Georeferenced similar blueprints before, but I'm not sure what I'm doing differently.


Labels for Very Small Features Not Displaying

I'm labeling features in a map book (MapLogic) in ArcMap using Maplex labeling. I have the "Never Remove..." box checked within Maplex for the layers, however, some of my features that have a very small area (e.g. 0.005 sqft) are missing labels - they're just not placing. All labels are dynamic. When I try to create annotations for those layers (tried selecting just the features whose labels weren't placing, and all the features within the map), no annotations were created for those unplaced labels (even when I selected 'create annotation for unplaced..') Is this a bug that others have experienced, or is there some tweek within Maplex settings that I need to do? Here're a few screenshots of one of the layer's Maplex label settings that won't place - the feature's area is <0.005 sqft. **EDIT** : This is in ArcGIS 10.1 ![screenshot1](http://i.stack.imgur.com/17d2T.jpg) ![screenshot2](http://i.stack.imgur.com/oST4J.jpg) ![screenshot3](http://i.stack.imgur.com/KAQ05.jpg) ![screenshot4](http://i.stack.imgur.com/zSZc7.jpg)


Is there a tutorial explaining the basics of WMS?

I'm building some software that needs to download satellite map images to display behind some survey charts. I understand that, to download these images, I need access to a Web Map Service. The only problem is, I can't find any tutorials explaining how to use one. Does anyone have any links to some basic tutorials on WMS? In addition, can anyone point me in the direction of some useful WMS services? Price/subscriptions shouldn't be too big of a problem. Edit: The program is written in C++, if that has any affect on anything. It will need relatively high detail imagery, enough to make out roads and houses. I don't have any experience with XML, which it seems most (all?) WMS use.


3d borehole view for ArcScene

I am trying to make a 3D Representation on the well from the following data. I cnat put it in a table work from there and use python. Just not sure the math. Anyone Help I do have the starting coordinates for the Borehole. LOG DEPTH The depth recorded on the field logs for the borehole.   TRUE DEPTH The true vertical depth corresponding to the above depth. Corrected from the start of the analysis.   HOLE TILT AND AZIMUTH The sampled borehole orientation.   AXIAL COORDINATES The coordinates North and East from the target origin.   POLAR COORDINATES The polar or radial coordinates of the borehole.   N.B. The reference point for all bearing angles on this listing is given at the top of each sheet. All interpretations are opinions based on inferences from electrical or other measurements and we cannot, and do not, guarantee the accuracy or correctness of any interpretations, and we shall not, except in the case of gross or wilful negligence on our part, be liable or responsible for any loss, costs, damages or expenses incurred or sustained by anyone resulting from any interpretation made by any of our officers, agents or employees. These interpretations are also subject to our general terms and conditions as set out in our current price schedule. All co-ordinates with respect to True North                   DEPT       TDEP          VTLT    VAZT        AXLN       AXLE     BBRG       BRAD     ~ASCII Log Data        1600.000    1600.00       0.2     100.8       0.00       0.01      101       0.01        1625.000    1625.00       0.3     141.9      -0.03       0.04      128       0.05        1650.000    1650.00       0.2     164.2      -0.05       0.05      138       0.07        1675.000    1675.00       0.2     341.4      -0.03       0.04      130       0.05        1700.000    1700.00       0.1     327.9      -0.02       0.03      122       0.03        1725.000    1725.00       0.2     284.2      -0.01       0.02      121       0.02        1750.000    1750.00       0.2     284.3       0.00      -0.04      275       0.04        1775.000    1775.00       0.2     323.7       0.06      -0.08      308       0.11        1800.000    1800.00       0.2     359.7       0.12      -0.08      326       0.15        1825.000    1825.00       0.2       0.6       0.20      -0.08      338       0.21        1850.000    1850.00       0.2     356.3       0.25      -0.08      342       0.27        1875.000    1875.00       0.2     355.8       0.33      -0.09      345       0.34        1900.000    1900.00       0.2     354.2       0.41      -0.10      347       0.42        1925.000    1925.00       0.3     350.7       0.50      -0.11      348       0.51        1950.000    1950.00       0.2     343.3       0.58      -0.13      347       0.60        1975.000    1975.00       0.3     341.9       0.68      -0.17      346       0.70        2000.000    2000.00       0.3     347.9       0.80      -0.19      347       0.83        2025.000    2025.00       0.3     353.7       0.93      -0.20      348       0.95        2050.000    2050.00       0.3     348.7       1.03      -0.22      348       1.06        2075.000    2075.00       0.3       7.6       1.13      -0.21      349       1.15        2100.000    2100.00       0.3     351.3       1.23      -0.23      349       1.25        2125.000    2125.00       0.3     344.1       1.35      -0.26      349       1.38        2150.000    2150.00       0.3     333.3       1.45      -0.31      348       1.48        2175.000    2175.00       0.3     334.4       1.54      -0.35      347       1.58        2200.000    2199.99       0.3     336.0       1.66      -0.41      346       1.71        2225.000    2224.99       0.4     340.1       1.81      -0.46      346       1.86        2250.000    2249.99       0.4     346.3       1.99      -0.50      346       2.05        2275.000    2274.99       0.4     333.4       2.16      -0.58      345       2.23        2300.000    2299.99       0.5     330.6       2.35      -0.70      343       2.45        2325.000    2324.99       0.6     336.0       2.58      -0.80      343       2.70        2350.000    2349.99       0.7     339.5       2.85      -0.90      343       2.98        2375.000    2374.99       0.6     337.2       3.09      -1.00      342       3.25        2400.000    2399.99       0.7     337.4       3.37      -1.12      342       3.55        2425.000    2424.98       0.7     340.2       3.65      -1.22      342       3.85        2450.000    2449.98       0.8     342.8       3.97      -1.32      342       4.18        2475.000    2474.98       0.8     343.2       4.29      -1.41      342       4.52        2500.000    2499.98       0.7     341.8       4.60      -1.51      342       4.84        2525.000    2524.98       0.7     340.8       4.90      -1.62      342       5.16        2550.000    2549.97       0.7     340.0       5.20      -1.73      342       5.48        2575.000    2574.97       0.8     337.1       5.51      -1.86      341       5.82        2600.000    2599.97       0.8     338.4       5.83      -1.99      341       6.16        2625.000    2624.97       0.8     338.4       6.14      -2.11      341       6.49        2650.000    2649.96       0.8     340.9       6.48      -2.23      341       6.85        2675.000    2674.96       0.8     335.6       6.79      -2.37      341       7.20        2700.000    2699.96       0.8     337.8       7.11      -2.50      341       7.53        2725.000    2724.96       0.8     336.7       7.42      -2.63      340       7.87        2750.000    2749.95       0.8     334.0       7.72      -2.78      340       8.21        2775.000    2774.95       0.9     335.9       8.08      -2.94      340       8.59        2800.000    2799.95       0.9     338.6       8.44      -3.08      340       8.99        2825.000    2824.95       0.8     333.7       8.77      -3.24      340       9.35        2850.000    2849.94       0.8     335.8       9.10      -3.39      340       9.72        2875.000    2874.94       0.8     334.1       9.43      -3.55      339      10.08        2900.000    2899.94       0.8     334.1       9.76      -3.71      339      10.44        2925.000    2924.93       0.8     331.3      10.07      -3.88      339      10.80        2950.000    2949.93       0.8     332.0      10.38      -4.05      339      11.14        2975.000    2974.93       0.8     328.4      10.67      -4.23      338      11.48        3000.000    2999.93       0.8     327.9      10.96      -4.41      338      11.81        3025.000    3024.93       0.8     328.1      11.25      -4.59      338      12.15        3050.000    3049.92       0.8     329.2      11.54      -4.76      338      12.48        3075.000    3074.92       0.8     328.7      11.83      -4.94      337      12.82        3100.000    3099.92       0.7     327.6      12.11      -5.11      337      13.14        3125.000    3124.92       0.8     328.3      12.39      -5.28      337      13.47        3150.000    3149.91       0.8     327.0      12.69      -5.48      337      13.82        3175.000    3174.91       0.8     330.4      12.99      -5.65      336      14.17        3200.000    3199.91       0.8     329.4      13.31      -5.84      336      14.54        3225.000    3224.91       0.8     329.7      13.62      -6.02      336      14.89        3250.000    3249.90       0.8     331.5      13.95      -6.20      336      15.26        3275.000    3274.90       0.9     333.3      14.28      -6.36      336      15.63        3300.000    3299.90       0.8     337.1      14.59      -6.50      336      15.97        3325.000    3324.90       0.8     339.6      14.91      -6.61      336      16.31        3350.000    3349.89       0.8     339.4      15.23      -6.73      336      16.65        3375.000    3374.89       0.7     342.4      15.53      -6.83      336      16.96        3400.000    3399.89       0.7     345.5      15.83      -6.91      336      17.27        3425.000    3424.89       0.8     345.4      16.15      -6.99      337      17.60        3450.000    3449.88       0.8     351.2      16.47      -7.04      337      17.92        3475.000    3474.88       0.8     351.2      16.82      -7.09      337      18.26        3500.000    3499.88       0.8     351.9      17.15      -7.14      337      18.58        3525.000    3524.88       0.8     352.5      17.52      -7.19      338      18.94        3550.000    3549.87       0.8     355.0      17.88      -7.22      338      19.28        3575.000    3574.87       0.8     353.9      18.25      -7.26      338      19.64        3600.000    3599.87       0.9     356.1      18.64      -7.28      339      20.01        3625.000    3624.87       0.9     355.7      19.03      -7.31      339      20.39        3650.000    3649.86       1.0     354.4      19.47      -7.36      339      20.81        3675.000    3674.86       1.0     353.8      19.90      -7.40      340      21.23        3700.000    3699.85       1.0     353.9      20.33      -7.45      340      21.65        3725.000    3724.85       1.0     355.1      20.75      -7.48      340      22.06        3750.000    3749.85       1.0     352.3      21.17      -7.54      340      22.48        3775.000    3774.84       0.9     348.8      21.55      -7.62      341      22.86        3800.000    3799.84       0.8     340.2      21.90      -7.74      341      23.23        3825.000    3824.84       0.8     340.5      22.24      -7.86      341      23.59        3850.000    3849.84       0.8     344.0      22.57      -7.96      341      23.93        3875.000    3874.83       0.8     346.3      22.91      -8.04      341      24.28        3900.000    3899.83       0.8     347.2      23.26      -8.12      341      24.63        3925.000    3924.83       0.8     349.5      23.60      -8.18      341      24.98        3950.000    3949.83       0.8     345.9      23.94      -8.27      341      25.33        3975.000    3974.82       0.8     348.7      24.27      -8.33      341      25.66        4000.000    3999.82       0.8     347.6      24.59      -8.40      341      25.99        4025.000    4024.82       0.8     349.6      24.94      -8.47      341      26.34        4050.000    4049.82       0.8     348.0      25.29      -8.54      341      26.69        4075.000    4074.81       0.9     346.7      25.66      -8.63      341      27.07        4100.000    4099.81       0.8     344.6      26.01      -8.72      341      27.43        4125.000    4124.81       0.9     343.6      26.36      -8.83      341      27.80        4150.000    4149.81       0.8     342.4      26.72      -8.94      342      28.17        4175.000    4174.80       0.9     341.2      27.08      -9.06      341      28.55        4200.000    4199.80       0.9     340.4      27.45      -9.19      341      28.95        4225.000    4224.80       0.9     341.3      27.84      -9.33      341      29.36        4250.000    4249.79       0.9     343.3      28.20      -9.43      342      29.74        4275.000    4274.79       0.8     343.3      28.55      -9.54      342      30.10        4300.000    4299.79       0.9     340.7      28.90      -9.66      342      30.48        4325.000    4324.78       0.9     337.4      29.27      -9.82      341      30.87        4350.000    4349.78       0.9     337.6      29.64      -9.97      341      31.27        4375.000    4374.78       0.9     336.7      30.02     -10.13      341      31.68        4400.000    4399.77       0.9     329.6      30.37     -10.33      341      32.08        4425.000    4424.77       1.0     331.8      30.73     -10.53      341      32.49        4450.000    4449.77       0.9     330.9      31.08     -10.72      341      32.88        4475.000    4474.77       0.9     332.7      31.42     -10.90      341      33.26        4500.000    4499.76       0.9     333.0      31.75     -11.07      341      33.63        4525.000    4524.76       0.8     331.9      32.08     -11.24      341      33.99        4550.000    4549.76       0.8     332.6      32.38     -11.40      341      34.33        4575.000    4574.75       0.7     333.6      32.66     -11.54      341      34.64        4600.000    4599.75       0.7     334.4      32.93     -11.67      340      34.94        4625.000    4624.75       0.6     333.6      33.16     -11.78      340      35.20        4650.000    4649.75       0.6     332.7      33.40     -11.90      340      35.46        4675.000    4674.75       0.6     333.2      33.64     -12.02      340      35.72        4700.000    4699.75       0.5     333.6      33.85     -12.13      340      35.96        4725.000    4724.75       0.6     333.5      34.08     -12.24      340      36.21        4750.000    4749.75       0.5     332.1      34.28     -12.35      340      36.44        4775.000    4774.74       0.5     329.7      34.47     -12.46      340      36.65        4800.000    4799.74       0.5     331.6      34.66     -12.56      340      36.87        4825.000    4824.74       0.5     333.0      34.87     -12.67      340      37.10        4850.000    4849.74       0.6     331.8      35.08     -12.78      340      37.33        4875.000    4874.74       0.6     331.6      35.29     -12.90      340      37.58        4900.000    4899.74       0.6     330.6      35.50     -13.02      340      37.81        4925.000    4924.74       0.6     330.7      35.73     -13.14      340      38.07        4950.000    4949.74       0.5     330.7      35.93     -13.26      340      38.30        4975.000    4974.74       0.5     328.3      36.12     -13.37      340      38.52        5000.000    4999.73       0.5     327.0      36.31     -13.49      340      38.73        5025.000    5024.73       0.5     328.8      36.49     -13.61      340      38.95        5050.000    5049.73       0.6     324.3      36.69     -13.75      339      39.18        5075.000    5074.73       0.6     326.0      36.89     -13.88      339      39.42        5100.000    5099.73       0.5     326.3      37.08     -14.01      339      39.64        5125.000    5124.73       0.5     324.9      37.27     -14.14      339      39.87        5150.000    5149.73       0.6     322.2      37.47     -14.29      339      40.10        5175.000    5174.73       0.6     320.6      37.66     -14.45      339      40.33        5200.000    5199.73       0.5     322.2      37.84     -14.60      339      40.56        5225.000    5224.72       0.5     321.5      38.02     -14.74      339      40.78        5250.000    5249.72       0.6     319.1      38.21     -14.90      339      41.01        5275.000    5274.72       0.7     320.2      38.43     -15.08      339      41.28        5300.000    5299.72       0.7     315.4      38.64     -15.30      338      41.56        5325.000    5324.72       0.8     319.4      38.91     -15.53      338      41.90        5350.000    5349.71       0.9     318.4      39.19     -15.78      338      42.25        5375.000    5374.71       1.0     327.2      39.56     -16.01      338      42.68        5400.000    5399.71       1.1     334.9      40.01     -16.22      338      43.17        5425.000    5424.70       1.4     328.5      40.52     -16.53      338      43.76        5450.000    5449.68       2.1     326.1      41.27     -17.04      338      44.65        5475.000    5474.65       3.0     325.9      42.32     -17.75      337      45.90        5500.000    5499.58       4.1     325.0      43.78     -18.77      337      47.64        5525.000    5524.47       5.4     324.4      45.68     -20.13      336      49.92        5550.000    5549.31       6.7     324.1      48.02     -21.83      336      52.74        5575.000    5574.04       8.3     323.6      50.93     -23.98      335      56.29        5600.000    5598.71       9.4     323.2      54.20     -26.42      334      60.30        5625.000    5623.32      10.1     322.3      57.64     -29.09      333      64.56        5650.000    5647.90      10.6     322.2      61.26     -31.89      333      69.07        5675.000    5672.44      11.0     322.1      65.01     -34.81      332      73.75        5700.000    5696.94      11.5     322.3      68.99     -37.89      331      78.71        5725.000    5721.44      11.5     322.0      72.90     -40.95      331      83.62        5750.000    5745.92      11.6     321.8      76.85     -44.06      330      88.58        5775.000    5770.42      11.6     322.0      80.79     -47.14      330      93.54        5800.000    5794.90      11.7     322.0      84.78     -50.26      329      98.56        5825.000    5819.36      11.8     321.6      88.84     -53.48      329     103.70        5850.000    5843.79      12.3     321.9      93.02     -56.75      329     108.97        5875.000    5868.21      12.4     321.7      97.21     -60.06      328     114.27        5900.000    5892.62      12.5     321.9     101.47     -63.40      328     119.64        5925.000    5917.02      12.6     321.9     105.73     -66.75      328     125.04        5950.000    5941.41      12.6     321.3     110.03     -70.19      327     130.51        5975.000    5965.78      12.9     321.1     114.37     -73.68      327     136.05        6000.000    5990.15      12.9     320.4     118.67     -77.24      327     141.59        6025.000    6014.44      13.8     320.9     123.27     -80.97      327     147.48        6050.000    6038.65      14.4     320.6     128.07     -84.91      326     153.66        6075.000    6062.77      15.2     319.9     133.14     -89.17      326     160.24        6100.000    6086.86      15.5     319.6     138.21     -93.49      326     166.86        6125.000    6110.93      15.7     319.0     143.31     -97.93      326     173.58        6150.000    6134.97      15.9     318.7     148.45    -102.44      325     180.37        6175.000    6158.97      16.3     318.9     153.72    -107.03      325     187.31        6200.000    6182.94      16.5     319.2     159.14    -111.71      325     194.43        6225.000    6206.93      16.3     319.3     164.44    -116.27      325     201.39        6250.000    6230.97      16.0     319.4     169.67    -120.74      325     208.24        6275.000    6255.01      15.9     319.5     174.86    -125.18      324     215.05        6300.000    6279.06      15.9     319.4     180.06    -129.62      324     221.86        6325.000    6303.12      15.7     319.2     185.21    -134.08      324     228.65        6350.000    6327.25      15.2     319.0     190.14    -138.36      324     235.15        6375.000    6351.43      14.7     319.5     194.95    -142.48      324     241.46        6400.000    6375.67      14.2     319.7     199.62    -146.44      324     247.58        6425.000    6399.95      13.8     319.8     204.16    -150.27      324     253.50        6450.000    6424.27      13.3     320.1     208.62    -154.00      324     259.30


Arc 10 field calculator scripts

I need scripts for the field calculator that:   1. extracts the first word from a column, i.e. "Jones School" - "Jones", and   2. A script that will capitalize this new field, "Jones" - "JONES". This is Arc 10.


GeoServer converting PostGIS timestamp with time zone

In a PostGIS view i have a field of data type "timestamp with time zone". In GeoServer this field is show with datatype "Timestamp" ![enter image description here](http://i.stack.imgur.com/WcxjJ.png) My problem is that i want the Timestamp to be show in my application exactly as it is in my Database as **UTM. +00**. What i get back from GeoServer via WFS request is: > 2010-08-20T16:56:15.095-04:00 What is stored in the PostGIS database is: > 2010-08-20 20:56:15.095+00 **How do i tell GeoServer not to convert the time so that i see the time exactly how it is in the database as UTM? Is it on Apache, System Time, GeoServer config?** Running GeoServer on latest stable Suse linux Ent.


How to divide a shapefile in QGIS

German original: ich beschäfte mich gerade mit QGIS. Ich habe über die Funktion mein Shapefile hinein geladen. Jetzt ist mein Problem: In dem Shapefile sind unterschiedliche Punkte, wie Punkte für einen Pfeil oder nur für eine Linie,... und diese sind alle in dem einen Shape. Ist es irgendwie möglich, dass man den Shape irgendwie auftrennt. Also das es statt einem dann angenommen drei gibt oder so? Ich hoffe, ihr wisst was ich meine. Könnt ihr mir weiter helfen? Translation: I am working currently with QGIS. I have loaded the shapefile into my function. Now my problem is: in the shapefile are different with different kinds of points (one for arrow, one for ...) and these are all in one shapefile. I would like to know how can I create a separate shape file for each kind of point. So there are three instead of one then accepted? Can you help me please?


How can I see the coordinate transformation parameters in QGIS?

How can I see (and edit if possible) the transformation method and parameters used by QGIS for "on-the-fly" transformation, and then for saving data in another CRS? Like in ArcGIS? Here is the similar question but answered differently, I don't want to create a new _coordinate system_ , I want to check the _transformation_ between existing in Qgis coordinate systems. In ArcGIS I can select between several transformation methods for the pair of CRS, and what about QGIS? I've changed the default WGS84 CRS to another local CRS. For example, if I set for the project EPSG 3763 (ETRS PT TM 06), described as `+proj=tmerc +lat_0=39.66825833333333 +lon_0=-8.133108333333334 +k=1 +x_0=0 +y_0=0 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs`, and add data in EPSG 20790 (Datum Lisboa) `+proj=tmerc +lat_0=39.66666666666666 +lon_0=1 +k=1 +x_0=200000 +y_0=300000 +ellps=intl +towgs84=-304.046,-60.576,103.64,0,0,0,0 +pm=lisbon +units=m +no_defs`, it is reprojected on-the-fly and looks more or less correctly. But if I cannot see is it really correct I want to know the transformation used. There are several reccomended methods of different accuracy, and I want to know which one was used (if it was one of them). (I am not fluent in QGIS, just learning). Thanks in advance!


Merging polygons for rendering with Geoserver

For my application, I need to merge the polygons of different records into one polygon, but only for the rendering (WMS Map). If it's not clear, here's the goal: a PostGIS table contains the buffers around linear structures (roads, railways,...). Each element of each the structure has its own buffer and I'd like to show on the map only the union of the buffers grouped by type (one for the rails, one for the roads,...). I've tried the following:   * having only one buffer object per kind in the table (multipolygon): very slow, since it requires to load the entire polygon each time it is accessed, and not very flexible when the source data is modified   * creating a view in the PostGIS table that groups every thing using the following query: SELECT min(restriction.id) AS gid, st_multi(st_union(restriction.geom)) AS geom, restriction.type FROM public.restriction GROUP BY restriction.type Also very slow, since the union is done each time the view is accessed   * Using Geoserver SQL View's feature. According to the following post, a where statement is added to the view to limit the request to the queried bounding box, but is even slower than the PostGIS view: https://getsatisfaction.com/opengeo/topics/postgis_parameterized_view_dynamic_bounding_box It is not possible to it with SLDs. Is there any way to improve the query or to achieve the desired result, but faster? I'm using the latest version of OpenGeo (3.1).


(Pgrouting) Finding the best path within some maximal length

Is it possible to find the best path (where cost is not length) with the length limit? For example, find the best path that is shorter than 10km? I currently use pgr_astar from Pgrouting to find the best path, so is there any function/trick in Pgrouting which is able to handle this?


how to produce xml file in SLD to add discrete colors and color legend

I have to produce a xml file, according to Styled Layer Descriptor (SLD) (for more informations you can see: http://www.opengeospatial.org/standards/sld). The file is required by GeoServer and it has the purpose to add a color legend. How would this xml file be written, according to standard SLD ? Can you indicate me some link with examples? Thank you for your attention. Best regards, B.O.


Grouping analysis of points to divide a polygon into n number of sub-regions

I'm attempting to create _n_ sub-regions from a polygon based on a cluster analysis of a bunch of overlapping raster layers (physical layers, eg: depth, currents, waves). Briefly, I can create a regular grid across the polygon, then extract physical attributes from the physical raster layers (eg: Gridspot or equivalent tool) then, run a cluster analysis restricted to _n_ number of clusters (in R or other stats package). I can then, identify each cluster-group, and plot them back in GIS (QGIS or ArcMap). I envisage however, that some cluster-groups will be dispersed (spatial outliers and not meaningful), whereas some will be clumped (worthy of being a sub-region). I could then manually draw around representative clumps to create _n_ sub-regions. Is there a tool like ArcMap 10.1 Grouping Analysis that can be run in QGIS? I'm only running 10.0. Are there suggestions of a better way to do cluster analysis of multiple raster layers, to create n sub-regions (bio-regions)?


Batch geocode with multiple services?

I've got a batch of addresses to geocode (on the order of 10,000), I've looked at various services for doing the geocoding, and each service seems to be able to locate a portion between 50% and 90% of the addresses with address-level specificity. Are there any advantages/disadvantages to getting the results from multiple services and trying to combine them?


Polymaps / Tilestache tiling grid mismatch

I have set up a basic PostGIS / Tilestache / gunicorn / nginx / Polymaps test environment, rendering my dataset from PostGIS to GeoJSON and shipping it out to be rendered in Polymaps. That works nicely, but for one small detail: Polymaps seems not to be requesting all tiles. If my dataset encompasses a world map, at zoom level 0, it requests 0/0/0.geojson, but not 0/1/0.geojson (the Eastern hemisphere). If I zoom further into the visible dataset, the right halves of it seem to disappear. With a more local dataset (Norway, display centered around lat/lon 60/10) nothing at all is displayed. One correct tile to load and display would for example be 2/4/0.geojson, which spans the whole dataset. Polymaps never loads any tile with an index above 3. If I load the tile directly, ignoring the tiling scheme, it is loaded and displayed correctly. Any clues to what I am missing? Polymaps code below:               var po = org.polymaps;          var map = po.map()     .container(document.getElementById("map").appendChild(po.svg("svg")))     .center({lat: 60, lon: 10})     .zoomRange([0, 20])     .zoom(10)     .add(po.interact())     .add(po.hash());      loading the tiled layer:               map.add(po.geoJson()      .url("/boundary/{Z}/{X}/{Y}.geojson")     .id("boundary")     .on("load", loadFeature))      loading the one tile directly (works):               map.add(po.geoJson()     .url("/test.geojson")     .id("boundary2")     .on("load", loadFeature))


Using OSM for routes?

Does anyone have any tutorial, link, video etc on how I can use Open Street Map for routes? I know google maps and saw that I can easily do that. But this application will have an internal module and we can not afford the license.


How to make a vectorlayer through of sentence SQL on PostGRES

How do I create a vectorlayer to through SQL on PostGiS on Python 2.7. I is working with QGIS 1.8 and I needed create a vectorLayer through SQL on PostGIS.


Computing distance of a georeferenced .jpg to each polygon (town) in a shapefile

I have a .jpg that is a highway that spans multiple towns, and a shapefile with boundaries of these towns. I have georeferenced the .jpg so that it is aligned with the shapefile of the towns. My ultimate goal is to compute the distance of each town--ideally from its centroid--to the nearest point along this highway. I am able to convert the georeferenced highway into a shapefile. I have converted the shpaefile of the towns to a set of points to ease the distance calculation. When I do the Near calculation under Toolbox, however, I do not get any distances computed. When I try to do a Near calculation, I get Null for the distance. When I try a Spatial Join, I get a -1. I cannot figure out what issue I can be having. I am using Arcmap in ArcGis 10.2.


How do I remove the orthoimage layer from the new USGS geospatial PDF topo maps

I would like to remove the orthoimage layer from a USGS topo geoPDF map. I have Adobe Acrobat Pro 9, but can't find a way to delete a layer. I've tried setting the layer to "Never export", then exporting the pdf, but the image is still there (and in the process I have somehow deleted the geo-referencing info.) Thanks


How can I create geometry in SQL2012 with custom SRID?

I want to create geometry using EPSG 2039 in mssql 2012. This srid doesn't exist in sys.spatial_reference_systems any idea how I can do it ?


How to convert shp to PostGIS dump?

**Given a file my_file.shp** such those from Gadm e.g. US administrative divisions, where the coordinate reference system is latitude/longitude with WGS84 datum. Input could actually be:   1. Shapefile,   2. ESRI personal geodatabase,    3. ESRI file geodatabase,    4. Google Earth .kmz,   5. R (SpatialPolygonsDataFrame) **How to convert this input into a WGS84 PostGIS dump ?**


Superimposition of label in CartoDB

I'm currently doing a carte with CartoDB where I use CSS in order to show label but my issue is that I have entities which are very close and their label are difficult to read because of their superimposition. I try some of the CSS properties' of CartoCSS (like text-min-padding or text- placement) but I failed. Here is my CSS code               /** category visualization */          #doc {        marker-opacity: 0.9;        marker-line-color: #FFF;        marker-line-width: 1;        marker-line-opacity: 1;        marker-placement: point;        marker-type: ellipse;        marker-width: 40;        marker-allow-overlap: true;     }          #doc[type="ONG"] {        marker-fill: #A6CEE3;     }     #doc[type="Think Tanks"] {        marker-fill: #1F78B4;     }     #doc[type="Politique"] {        marker-fill: #B2DF8A;     }          #doc::labels {       text-name: [name];       text-placement: point;       text-face-name: 'DejaVu Sans Book';       text-size: 20;       text-fill: #000;       text-halo-fill: #FFF;       text-halo-radius: 1;       text-dy: 0;       text-allow-overlap: true;       text-placement-type: simple;       text-label-position-tolerance: 0;     }      Thanks a lot. Have a nice day A.H


Don't cancel the click event on a SelectFeature

On my OpenLayers map I have a custom Control that has a ClickHandler. I also have a SelectFeature on one of my layers. I need the ClickHandler to fire even if a SelectFeature is clicked. I tried:   * Switch the SelectFeature to hover only.   * unregister the click event on the SelectFeature from the map.events (in 100 different ways, I think the best one being in the onSelect of the SelectFeature: map.events.unregister('click', this.handlers.feature, this.handlers.feature.click); Nothing helps. If I click 3 pixels to the right, next to the marker, the ClickHandler fires perfectly fine. How can I tell the SelectFeature to NOT cancel the click event? I dont mind wether it fires or not. A long as my ClickHandler also fires.


Apply JQuery in Django Forms

I have an application for web mapping. I had already a maps(land parcels) displayed in django templates. [see image below] ![enter image description here](http://i.stack.imgur.com/wTjTt.png) Now, what I want to achieve is to display a single parcel by searching PIN through the use of dropdown list in django. [see image] ![enter image description here](http://i.stack.imgur.com/Sr3Cu.png) I want to achieve like this: [see image below] ![enter image description here](http://i.stack.imgur.com/rOtIN.png) This is my database: [see image] ![enter image description here](http://i.stack.imgur.com/nfPAW.png) I read an article applying the JQuery is much preferred., Is anyone knows how to integrate JQuery in Django?? and how to achieve the image below.,?? thanks in advance guys,. Any help is very much appreciated..


ValueError using MapProxy decorate image middleware with WMS

I am successfully using Mapproxys decorate image middleware to apply a watermark to WMTS tiles. I am using the following code:               def annotate_img(image, service, layers, environ, query_extent, **kw):              # Required opacity of watermark         opacity = float(0.5)              # Get path to watermark png         wm_path = Image.open(os.path.join(os.path.dirname(__file__), 'watermark_256.png'))              # Grab tile and Convert to RGBA         img = image.as_image()         if img.mode != 'RGBA':             img = img.convert('RGBA')              # Create new RGBA image for watermark         wm = Image.new('RGBA', img.size, (0,0,0,0))         wm.paste(wm_path)              # Adjust opacity of watermark         alpha = wm.split()[3]         alpha_c = ImageEnhance.Brightness(alpha).enhance(opacity)         wm.putalpha(alpha_c)              # Create new composite image using the wm as a mask         compo = Image.composite(wm, img, wm)         final_img = Image.new('RGBA', img.size, (0,0,0,0))         final_img.paste(compo)              return ImageSource(final_img, image.image_opts)      As I said this is working fine for WMTS, however, when I try to view a WMS layer I get the following error in my log:               self.im.paste(im, box)     ValueError: images do not match


OpenLayers - How to construct OpenLayers.Bounds

I would like to create Bounds for a layer but I am not sure in what order I have to provide the min_x / max_x and min_y / max_y coordinates. is it correct to put it in that order: min_y, min_x, max_y, max_x               var bounds = new OpenLayers.Bounds.fromArray([{{ wms.min_y }}, {{ wms.min_x }}, {{ wms.max_y }}, {{ wms.max_x }}, ]);


Feature Layer (Polygon) Rendering issue in ArcGIS Server 10.1

I have a feature layer consisting of 96,000 features, but we are displaying a maximum of 13,000 features in our map application. There is an issue that for some zoom levels the features are visible and for other, intermediate zoom levels the features are not visible. I have deleted the feature class and restarted the service in case there is issue with the feature class, but that also did not help.Thank you in advance and any kind of help will be appreciated.


GeoNetwork 2.6 and Geoserver CRS codes

I am a relatively new convert to the GeoNetwork (on 2.6) and presently setting up a server in the UK but have been having issues with the crs code used in , I think this is where I need to make the changes but please correct me if wrong! What I am trying to do is to project the data in the "mapViewer" and "mapSearch" (within config-gui.xml) in a projection other than epsg:4326, which for my region is epsg:27700 (the UK), so that when a viewer looks at the map the data are projected in a manner familiar to the viewer. However, when I set the crs code to 27700 I find that the map elements of GN do not load. I would be very grateful if somebody could please assist as I have tried many different options and appear to be getting nowhere - are there some other files that I should be modifying? Many thanks in advance for any assistance with this. Mark


Is it possible to create add-ins with VS Express?

Last year, when version 10 was still beta, I contacted ESRI to recommend a migration path for VBA conversion. I was reassured that the free VS Express will be sufficient. Now in version 10 and with all VBA customization at "document level" since UI Controls cannot be added to the "Normal.mxt" template created a significant problem with supporting / upgrading VBA customizations. Seems that I finally need to migrate the tools so i set out to complete the Walkthrough example only to find that I don't seem to be able to complete it. The Walktrhough is obviously customized for those with full version of VS. I have a very limited knowledge of .Net and so i thought the Add-ins should be simple enough migration path. Is it even possible to create an Add-in with VS Express? Does anyone know of a good example using VS Express rather than the full version.


ArcGIS Desktop's registry value under "software" key[CLOSE]

> **Possible Duplicate:**   >  How can I programmatically get the path of “Python.exe” used by ArcMap I am trying to get some information from registry about ArcGIS. I am using ArcGIS Desktop 10.0 and For this the registry looks like: ![registry key for arcgis](http://i.stack.imgur.com/gIXWX.png) My questions are:   1. What does it mean when the "Python10.0" subkey has a "True" value?   2. When "Python10.0" subkey has a "False" value?   3. Does the registry key values same for every version of ArcGIS? * * * **[Close]** I have edited How can I programmatically get the path of “Python.exe” used by ArcMap question. Thanks @blah238 and @whuber♦ to get me the wrong. Now, this is a close question. Please help me for previous one. * * * ### Closely related At How can I programmatically get the path of "Python.exe" used by ArcMap I asked how to obtain the path to `Python.exe` programmatically within ArcMap. Of the several solutions offered, one of them suggested using the registry, prompting the present question.


Anyone used Informatica Powercenter and FME together?

FME apparently integrates with Informatica Powercenter. I use Informatica quite a bit to move data around, push and pull to ArcSDE, and to slurp mass quantities of XML into SQL Server. I'm particularly interested in extracting DWG attributes and features and pushing them to featureclasses in SDE and into tabular SQL Server databases as well. Does anyone have any experience with anything like this or just Informatica and FME together, in general? I'm currently extracting data from DWGs through ArcObjects, and it can be painful. So, I guess my specific question to anyone with experience would be just how well does FME integrate with Informatica?


OpenLayers get single layer from WMS multiple layers

Openlayers enables to declare multiple layers in one WMS layer class.                var nasa_wms = new OpenLayers.Layer.WMS(             "NASA Global Mosaic",             "http://wms.jpl.nasa.gov/wms.cgi",              {                layers: "modis, global_mosaic"             });      When I add this layer to map, the items are coming together on the map. But sometimes I want to use a single layer to apply CQRS filter to **"modis"** layer or " **global_mosaic** ". So I should use a single wms layer and add CQRS filter. Is this possible, or should I separate layers as two layers.


How to add north arrow to map using ArcGIS silverlight API?

I think it's a simple question but i don't answer it. I want to add North arrow in my webgis using silverlight API. How can I do it? Is there any sample code to add north arrow using silverlight api ?


how to display movements of a moving object in cartodb

I have few GPS tracks with lat long and timestamp. I want to present them on a map and to show how the movement took place using the time series location values (X,Y,T). I searched a lot but could not find any suitable map visualization technique for moving objects. Cartodb offers visualization on the web with an appealing presentation. I was wondering if it can be done in Cartodb? Any suggestion on how to visualize the movement behaviour on a map would be highly appreciated. I am looking to develop something like this- http://here.com/livingcities/ Click on the metro cities on the webpage (the link provided above) and you can see the visualization.


Using feature service with JavaScript API and ArcGIS for Server on localhost?

I installed ArcGIS 10.2 for Server on my localhost and created and published a feature service on it. Now I want to use the feature service I created on a simple html page using Javascript. How can I do that? I read the resources on arcgis.com but it just described how to do it with the maps on arcgis.com! I don't need a base map from arcgis.com. actually I want to use the web page offline(just on localhost).


Why does a cache service need to load first when a dynamic service is involved?

I notice if you have 2 services 1 cached and 1 dynamic you have to load the cache first for it to display. Does anyone know why?


50k Mexico Topo Maps Sources?

Anyone know of another source (preferably free) for these 50k Mexico Topo Maps other than www.omnimap and www.eastview? There is a free WMS that has them but only for UTM 15 and 16 here: http://antroposig.ciesas.edu.mx:8399/arcgis/services/MapotecaDigitalCIESAS01/MapServer/WMSServer The areas we need are in UTM 12,13,and 14. Carta topográfica 1: 50 000


How to set up QGIS print composer options from the python console?

If I have a print composer called "ABC", I would like to "tell" its composer that I would like to   1. update "Map 0" to simulate the "set to map canvas extent" button activation   2. update certain text labels   3. activate "export to image"  Thank you.


best free easy GPS software? (not necessarily 'real' GIS)

Our search and rescue group is looking for some advice on what free and user- friendly (no computer savvy required) software folks could use to interact with their own GPS. Doesn't have to be a 'real' GIS tool. The main features we're looking for are downloading, uploading, saving, viewing, printing, and organizing of tracks, waypoints, and routes. We'd also like to list some pro's and con's of each. Garmin BaseCamp is an obvious choice. Trimble Outdoors got another recommendation so far - I guess parts are free and parts are subscription-based? Any others? Thanks UPDATE: while this calls for subjectivity, which can be seen as out-of-place for this web site, there's plenty of precedent: some of the most popular long- standing posts on this site call for subjectivity (what makes a map beautiful; what free tools should every GIS user have; etc) - subjectivity is important too when used with judgment. So - thanks for your opinions, keep em coming. Smart phone apps too.


Can we compare a Service in Arc GIS server with normal webservice

In Arc GIS server there are many types of services are available (like feature service, Map Service, Geo-processing service etc.,). My question is in general **web service** whether it is ASP.NET web service or WCF service are like same or different with respect to map services in Arc GIS server.


Identifying shapes in TIGERweb WMS SVG

I'm using the US Census WMS service describe here. Using code similar to that found in this example, but specifying a data format of "image/svg+xml", I've been able to get the TIGERweb WMS service to give me SVG that is rendered in my map. However, each shape in the SVG just the actual path, with no identifying data. Is there any way to get TIGERweb to annotate the shapes it serves so I can tell what they are (and process them accordingly with JS, CSS, etc.)? Like, I would like to get an SVG where each path is the boundary of a census block, tract, etc., and where the each path element in the SVG has an attribute giving the block number, FIPS code, or whatever piece of info will allow me to identify what the shape represents. As it is, I can't see how one could use the TIGERweb SVG data for anything except a plain rendering of the boundaries of counties, tracts, or what have you. Without any identifying information, there's no way to spiff up the shapes with special styling to represent data associated with the census areas.


Offline SLIPPY Map (Tiles) Database for Leaflet

How can create an offline cache of a large number of SLIPPY tile (PNG, 256x256) inside a web browser and then display them from that offline cache. For example, can one use HTML5 IndexedDB as an storage location for offline map tiles?


OpenLayers label without using StyleMap

I need to show label of each feature based on attribute field value. I know that I can use ${attribute} notation for this purpose. But this approach works only if I make StyleMap and pass it to vector layer constructor. But in this case I can not modify style of each features separately. How I can make layer with labels and with not null value of style property of each feature? My current config:               var custom_style = {         fill: true,         fillColor: "#feb24c",         fillOpacity: 0.2,         strokeColor: '#3182bd',         strokeWidth: 1     };          var l = new OpenLayers.Layer.Vector("Coverage", {         strategies: [new OpenLayers.Strategy.BBOX()],         protocol: new OpenLayers.Protocol.HTTP({             url: "http://***:8088/collect",             format: new ol.Format.GeoJSON()         }),         style: custom_style     });      If I add `label: ${attribute}` to custom_style object, then text `${attribute}` will be shown on map, not value of attribute.


Using Geometry commands in MS SQL Server 2008

I am looking to find the most effect way to use SDE and ST_Geometry SQL commands. Can anyone tell me how/what commands I might use to perform a spatial query in MS Sql Server 2008, to Select property parcels that are contained by or cross the boundary of a flood area layer? I then have to export the result (containing, owner, and address details) to some format that can be used to do a mail out. I would also have to create a site specific map for each property affected. Showing the neighbourhood parcels, specific owners parcel, flood potential layer and roads. Any help would be much appreciated.


GPS tracking using smartphones

We're looking to launch a study that examines spatial behaviour on a university campus. Given the near ubiquity of smartphones among our study population, we're excited about the possibility of using GPS tracking via a smartphone app to garner microscale travel behaviour data that could be analysed in conjunction with retrospective interviews. We're keen to minimise privacy intrusion (and battery drain) by only tracking behaviour on campus. We'd also like to minimise responded burden by making the application operate silently "in the background" and easy to opt-out of at any time after giving consent. What apps are available to collect this data? Are there any ethics or privacy issues we should be aware of?


How to populate the geometry column in a PostGIS table?

After loading a polygon shapefile into PostGIS, the geometry column (geom) is blank for several features. In one case, only a single feature has the geometry calculated and all others are blank. I'm using the PostGIS SHP and DBF loader v.2.0. I am creating a spatial index, and using COPY rather than INSERT. I have loaded several other point datasets and they geom column is populated just fine... (Otherwise, things are working just fine for the polygons: when I create a new feature, the geometry is calculated / geom field populated...etc.) Have I done something wrong when importing the data? * * * ### Solution by @ThomasG77 (in comments) PgAdmin3 [doesn't] display features that are too complex but they are present. You can try `SELECT length(ST_astext(the_geom))` from your_table to confirm or you can find "true" empty geometry with `SELECT * from your_table WHERE the_geom IS NULL;`.


How to take point data from a network and input it as TO/FROM endpoint data in streets layer?

I am building a network of my city's roads in ArcMap 10.1. Right now I have three layers: 1) Streets_To-From (lines representing the southbound and westbound sides of every road and one-way streets in the same direction), 2) Streets_From-To (lines representing the northbound and eastbound sides of every road and one-way streets in the same direction), and 3) Nodes (points representing all of the intersections in my city). My goal here is to take the OBJECTID data of all the nodes and put them into the line data of each layer. For example, along with the data I already have for each road segment (speed limit, street name, etc), I want to have fields named "FromNode" and "ToNode" which will correspond with the nodes in the Nodes layer. Since there is a duplicate of every two-way road, a lot of the data will be the same except "FromNode" and "ToNode" will have opposite values. Is there a way to do this other than manual input? Thanks very much and please ask me questions if this doesn't make sense!


Enhance performance of a join query for mapnik using postgis

I use mapnik to render map with datasources from a postgis database. One query is a join on two table, one with a geometry (a polygon) and the other one with some fields I want to use for rendering. Unfortunately, this query is very slow... Here is the slow query:               SELECT table1.id, table1.geometry, table2.fk_table1, table2.field FROM table1     INNER JOIN table2 ON (table1.id = table2.fk_table1)     WHERE table1.geometry && SetSRID('BOX3D(x1 y1, x2 y2)'::box3d,900913);      Here is the query plan I got from an EXPLAIN:               Hash Join  (cost=118896.94..262210.34 rows=27942 width=1099)         Hash Cond: (table1.id = table2.fk_table1)         ->  Seq Scan on table2  (cost=0.00..73740.89 rows=4786689 width=16)         ->  Hash  (cost=113876.63..113876.63 rows=33865 width=1083)             ->  Bitmap Heap Scan on table1  (cost=1112.15..113876.63 rows=33865 width=1083)      This query can be faster if the bbox is first applied on table1.geometry. However, mapnik does not help me there. Indeed, mapnik builds a query from a table/subquery by adding a WHERE clause on it                 SELECT AsBinary("geometry") AS geom       FROM subquery        WHERE "geometry" && SetSRID('BOX3D(x1 y1,x2 y2)'::box3d,900913)      With the following subquery:               ( SELECT table1.id, table1.geometry, table2.fk_table1, table2.field FROM table1     INNER JOIN table2 ON (table1.id = table2.fk_table1) ) as subquery      Basically, I need to enhance the speed of this query knowing that mapnik will add a where clause on the subquery. Thx for your help.


How to save a Leaflet or OpenLayers map for opening later?

Can a webmap state be saved to a text file or something for restoring later? Like if I set up a map with a particular zoom and choose my layers and draw a polygon or two, can it be saved to a file and reopened in the same state later on? Thanks


Overview of GPS tracking system technologies

Background: I've been asked by a friend to estimate the cost of setting up a GPS tracking system for a small fleet of buses, so we can do real time reporting of bus locations via a webmap. In my research, I saw three kinds of active GPS trackers:   1. Radio   2. SMS   3. Satellite Some major considerations are   1. minimizing long term costs, and    2. making the system as reliable as reasonably possible.  Secondary concerns are minimizing up front cost, latency, difficulty of installation/setup, level of support required after the initial setup is complete. The buses run within a 5 mile radius, so radio is feasible. With an eye towards reducing long term costs, radio seems good: it would be the only way to run the system without a middleman. SMS requires cell phone towers, and satellites require a space program. I am concerned that a DIY solution with radios might prove to be less reliable, though. But if we bought SMS transmitters that provided the data through a third party API, they could go out of business, or have unacceptable down time or latency. * * * **My question(s):**   * What would each of these three solutions look like?    * What hardware is involved?    * What are the trade-offs between them?    * What problems have you encountered?    * Have you used a system that _Just Works_?


Remove stale/obsolete geometry columns from the geometry_columns table

I understand that the function `Populate_Geometry_Columns()` inserts records in the `geometry_columns` table for geometry columns that are not yet listed there. However, an invocation of `Probe_Geometry_Columns()` reports:               select probe_geometry_columns();                   probe_geometry_columns                    ---------------------------------------------------      probed:1086 inserted:0 conflicts:1086 stale:53612     (1 row)      This is after a call to `Populate_Geometry_Columns()`. Is there an easy/documented way how to get rid of the conflicts and the stale entries? This has been posted a while ago to dba.sx with no replies: http://dba.stackexchange.com/q/14975/6934.


PostGIS Select Union of Polygons

I have a Postgres 8.4 PostGIS 2.0 enabled database and I've loaded it up with a bunch of shapefiles. I'm trying to get the database to do on-the-fly unions of geometries. Here is the query I'm working with:               SELECT ST_AsGeoJSON(ST_Union(the_geom)) FROM (SELECT geom FROM forecast_zones WHERE      state = 'ND' AND zone IN ('002','011','021','035','046','050')) as the_geom;      The Subquery to pull the geometries of the individual zones works fine:               SELECT geom FROM forecast_zones WHERE state = 'ND' AND zone IN      ('002','011','021','035','046','050');      But whenever I try to pass that to ST_Union I get undefined function errors. **EDIT:** Here is the query I ended up using that works:               SELECT ST_AsGeoJSON(ST_Union(forecast_zones.geom)) AS geom FROM forecast_zones     WHERE forecast_zones.state = 'ND' AND forecast_zones.zone IN     ('002','011','021','035','046','050')     GROUP BY forecast_zones.state


Reference materials for learning Google API V3 with PHP and MySQL (tracking)

I am new to this Google API. Now i wanted to learn Google API with PHP and MySQL. Actually i wanted to create an application that will track x,y coordinates coordinates on real time bases for android mobile. There will be many Android mobiles which will send data on real time bases and using php and mySQL i wanted to display it on Google map and refresh display and data after certain time interval. So with respect to this tracking system part can any one suggest me any references material and web site.


Import polylines and polygons into QGIS from CSV

I'm assuming this is a simple question, though my efforts have not yet found the answer! I want to import a large number of polylines and/ or polygons into one one or more QGIS layers in order to plot the data. Each will only have simple attributes (e.g. name). I need to know:   1. The required format of the CSV file to enable step 2.   2. The command or plug in required to create a polyline or polygon layer in QGIS. For information, the input CSV files are written by bespoke Excel VBA or by FORTRAN codes, so I have full control of these output formats. Essentially I'm looking for the polyline and polygon equivalent of the delimited text plugin. Thanks!


How do you emulate a right click on a tablet/mobile device using OpenLayers?

I'm trying to make my OpenLayers mapping application work for mobile devices and one feature I haven't figured out how to emulate is the 'handleRightClicks' Control Click option. Anyone know how to make this work with tablets/mobile devices? Thanks!


How do I add TIFF files to one layer/raster catalog?

We have several TIFF images (already tagged with position information) we want to load as a single layer into the map. It's the first time I dig into this raster/topology thing, but at least I managed to create a layer with one image as an `IRasterDataset`:               Type factoryType = Type.GetTypeFromProgID("esriDataSourcesRaster.RasterWorkspaceFactory");     IWorkspaceFactory wsFactory = (IWorkspaceFactory)Activator.CreateInstance(factoryType);     IRasterWorkspace2 rasterWS = (IRasterWorkspace2)wsFactory.OpenFromFile(@"c:\sourcefiles\");     IRasterDataset rds = rasterWS.OpenRasterDataset(@"example.tif");         IRasterLayer rasterLayer = new RasterLayerClass();     rasterLayer.CreateFromDataset(rds);     _map.AddLayer(rasterLayer);      This works well but I need to add all TIFF files to this raster layer. So I assume that `IRasterCatalog` could suit my needs. But my attempt to use it fails already on creation of the catalog object (taken from sample):               Geoprocessor geoProcessor = new Geoprocessor();         CreateRasterCatalog createRasterCatalog = new ESRI.ArcGIS.DataManagementTools.CreateRasterCatalog();          //Set parameters     createRasterCatalog.out_path = outPath;     createRasterCatalog.out_name = catalogName;     createRasterCatalog.raster_spatial_reference = map.SpatialReference;     createRasterCatalog.spatial_reference = = map.SpatialReference;              //Execute the tool to create a raster catalog     geoProcessor.Execute(createRasterCatalog, null);      On Exceute() I get a COMException, which isn't very helpfull. I tried to set the SpatialReference to `IGPCoordinateSystem` objects. I made sure that the output path exists. Still the COMException raises. Has anyone a clue for me? Is this the right approach at all?


Is there a QGIS mailing list that only announces new releases?

Is there a mailing list for qgis that announces only new releases / new binary packages? Thanks.


ArcGIS Tracking Server

I want to create a vehicle tracking application. I have ArcGIS Desktop and ArcGIS Server. Can I create application for real time with these two software? I saw in ESRI site another server called Tracking Server. Does ArcGIS Server include Tracking Server or the Tracking Server is another product? Thanks


Change order of added FeatureGraphicsLayer in GraphicsDataSet

I am adding two FeatureGraphicsLayer to a GraphicDataSet like so:                graphicsFunctionality.GraphicsDataSet.Tables.Add(featGraphicsLayer); ` It looks like it puts the last one in on the very top of all the other FeatureGraphicsLayers. I need to add one to the bottom of the GraphicsDataSet. Is this possible? Is there a ways to change the order of FeatureGraphicsLayers once they are in the GraphicsDataSet?


Can I edit WMS layer or KML layer?

I am using Openlayers and until now i have been using OSM as background and i got points from my db (mySQL) and add them to vector layer, then i could modify them and save it back to the server. now i am not using OSM anymore, I got data as shp file, and this is the background of the map, I am using GeoServer to show the files. I want to make an option to modify this layers (WMS). I cant insert everything to the database because its too much and the map will work slowly if i will add the data from the database. So i need a way to modify the vectors that i am getting from GeoServer and to save the changes, or if its not possible, maybe to do it with KML layer(but this is a bad option). there is a way to do it? I haven't find anywhere....


How to correct a wrong projection of shapefile?

i have a shape file wit UTM zone 16. This projection is wrong. It must be UTM zone 32. How can i correct the projection? When i convert (f-tool) to UTM zone 32, all the vector are rotated. The property of the layer shows funny numbers that is not correct:               xMin,yMin -4417494.80;6867783.98 : xMax,yMax -4416554.49;6869279.82      Correct property shold be:               xMin,yMin 3410596.14;5317812.48 : xMax,yMax 3411876.35;5318694.50


How to convert from ITM ref to GPS ref ?

How can I convert from ITM ref (Israel new map) to GPS ref and from GPS ref to ITM ref? I need any C# sample code or formula for this. For example: 32.086156,34.769239 = 178364,665948 Or : 29.548783,34.952316 = 194997,384546


Specify world file to use in ArcGIS

I have a BMP file where I know the coords of the corners, and hence have a generated a World File according to ESRI's standards. It has the the right filename, but when I add the image in ArcMap, it complains there's no spatial information. Things I've tried:   * Turned on the "Use world file to define the coordinates of a raster" option   * Different extensions for the world file - BMPW, BMW, BPW   * Trying the help... the directions it gives are for a _link_ file, which has to be imported manually for each raster (I'll have hundreds to these rasters) So the question is: how do I get ArcGIS to use the world file?


Cadastral Fabric 9.3

I am using ArcView without any extensions. Is there an alternate program out there which uses Cadastral Editing (Parcel Editing) or something similar? I downloaded the ArcEditor 10.0 trial version but the parcel editing does not function as I don't know if its due to it being trial. Thanks.


What is a good resource for making choropleth maps using excel data?

I would like to make a choropleth map of Malaysian regions, colour-coding them by population density. Is there a place where I can find templates of such a map?


How do I build a custom form for entering attributes in QGIS

I've done this many times in ArcGIS and ArcPad but I don't even know where to start building a custom form for data entry in QGIS. I would like to do this for field operations--QGIS would make a good field tool with a little customization


OID for the last added feature

How can i get the Object ID for the last feature created ? i use `m_editEvents_OnSketchFinished()` event to run some code, but i can't get the OID of the most recent added feature.


Qgis : QGraphicsTextItem's position not updated?

I'm using      `QgisInterface->mapCanvas()->scene()->addItem(QGraphicsTextItem)`      to add a QGraphicsTextItem to the map canvas. The problem is that when I use the pan tool to move the display the added QGraphicsTextItem is not moving !   How can I make it update its position whenever the pan tool is used on the map ?


How to display a particular area?

Hy, I'm developing a plugin in python for Qgis. I would like the user to select a particular area in a list (let say a country list) and zoom in the map on the corresponding area. Actually I'm doing it 'ugly':I create a layer (from a wkt), zoom to this particular layer, and delete this layer... Is there a way to do it properly ?


cannot submit address to Google geocoder on iPad: ArcGIS JS webmap app

I am working on a web map app for desktop and iPads. The app uses a Google geocoder to locate addresses or intersections, added with the standard code in 'init': geocoder = new google.maps.Geocoder(); Everything about the geocoder works fine on my desktop using IE9, Chrome, or Firefox 20.0. But on an iPad (I am testing with one using v 5.1), I cannot submit the address I type in the text box...The 'Go(enter)' button does nothing, and no clicks or presses on the screen do anything. Oddly, I can enter and successfully search for addresses using the same app on an iPad mini, by pressing 'Go' (v 6.1.3) ... Has anyone encountered this problem, or does anyone have any suggestions? I am brand new at testing for iOS! I can add the geocoder code, but there is kind of a lot of it..I don't think the error is in the geocoder setup, but I am not sure... // J


Highlight paths between two points

I'm new to qgis and very sorry if this is a simple question to ask. IS there a way to highlight the available paths (roads) between two selected points in a road map(in qgis) ? I want to find the possible routes between two points in map, which less than or equal to a specified distance. For example, highlighting the roads(paths) between point A and point B which are less than 5km (if available). And the two points are selected by the user on map, using mouse clicks.


Arcgis 10.0 desktop: how to clear All transformations from map

I have a custom extension which adds some custom transformations to the map. When turned off, the extension is supposed to remove all the added transformations.   However, when calling `IGeoTransformationOperationSet.RemoveAll()`, nothing happens. The transformations are removed from the set but nothing changes in the map - the transformed features stay in the same place and do not move back to their original position. I'm using Arcgis desktop 10.0 and my code is written in c#. How can I turn off/on the custom transformations? What am I missing? Thanks.


r.regression.multi - ERROR: Not enough valid cells available

I am running a few analyses with r.regression.multi in Grass 7. For a few rasters that I try to correlate, the following error comes up, however: "ERROR: Not enough valid cells available". I cannot see how the rasters I am trying to correlate in theses cases differ from the rasters where the analysis works. They match in dimension and location. Any ideas where the problem could lay? Addition: Here a screen shot of another situation: the regression is only calculated for a certain area. ![enter image description here](http://i.stack.imgur.com/vuDm9.jpg)


Copy .prj file into QGIS projection folder

I want a custom projection used in ArcGIS be available in QGIS. I have the .prj file copied out of a ArcGIS folder (containing different projections). Now, I want this file to be added into QGIS. Is there a folder where I can copy this file into? The custom projection is:               PROJCS["ETRS_1989_UTM_Zone_33N7",     GEOGCS["GCS_ETRS_1989",     DATUM["D_ETRS_1989",     SPHEROID["GRS_1980",6378137.0,298.257222101]],     PRIMEM["Greenwich",0.0],     UNIT["Degree",0.0174532925199433]],     PROJECTION["Transverse_Mercator"],     PARAMETER["False_Easting",3500000.0],     PARAMETER["False_Northing",0.0],     PARAMETER["Central_Meridian",15.0],     PARAMETER["Scale_Factor",0.9996],     PARAMETER["Latitude_Of_Origin",0.0],     UNIT["Meter",1.0]]


Is there an Attribute Table Feature limit?

I am copying a table with 11,000+ features but only about 700 show up on the QGIS attribute table. Is there a limit? how can I avoid this problem? Thank You in advance!


How to add nodes into graph

I'm new to qgis and in here I want to find a path between two selected points on the map(Roads-vector layer). The points are selected by the user, using mouse clicks. So here I used the astar algorithm to find path between two points. *********************************astar.py************************************   import heapq class AStar(object):   def **init** (self, graphAstar):   self.graphAstar = graphAstar               def heuristic(self, node, start, end):         raise NotImplementedError          def search(self, start, end):     openset = set()         closedset = set()         current = start         openHeap = []         openset.add(current)         openHeap.append((0,current))         while openset:         temp = heapq.heappop(openHeap)         current = temp[1]             if current == end:                 path = []                 while current.parent:                     path.append(current)                     current = current.parent                 path.append(current)                 return path[::-1]             openset.remove(current)             closedset.add(current)             for node in self.graphAstar[current]:                 if node in closedset:                     continue                 if node in openset:                     new_g = current.gg + current.move_cost(node)                     if node.gg > new_g:                         node.gg = new_g                         node.parent = current                 else:                     node.gg = current.gg + current.move_cost(node)                     node.H = self.heuristic(node, start, end)                     node.parent = current                     openset.add(node)             heapq.heappush(openHeap, (node.H,node))         return None               class AStarNode(object):         def __init__(self):             self.gg = 0             self.H = 0             self.parent = None              def move_cost(self, other):             raise NotImplementedError      *****************************astar_grid.py*******************************               from astar import AStar, AStarNode     from math import sqrt          class AStarGrid(AStar):         def heuristic(self, node, start, end):             return sqrt((end.x - node.x)**2 + (end.y - node.y)**2)          class AStarGridNode(AStarNode):         def __init__(self, x, y):             self.x, self.y = x, y             super(AStarGridNode, self).__init__()              def move_cost(self, other):             diagonal = abs(self.x - other.x) == 1 and abs(self.y - other.y) == 1             return 14 if diagonal else 10      and in the main code, the following method is used to create graph from vector layer. ****************************plugin.py************************************               def make_graph(self, mapinfo):         nodes = [[AStarGridNode(x, y) for y in range(mapinfo['height'])] for x in range(mapinfo['width'])]         graphAstar = {}         for x, y in product(range(mapinfo['width']), range(mapinfo['height'])):             node = nodes[x][y]             graphAstar[node] = []             for i, j in product([-1, 0, 1], [-1, 0, 1]):                     if not (0 <= x + i < mapinfo['width']): continue                     if not (0 <= y + j < mapinfo['height']): continue                     graphAstar[nodes[x][y]].append(nodes[x+i][y+j])         return graphAstar, nodes      And I called that method in FindRoutes method..               def findRoutes(self):     # QMessageBox.information( self.iface.mainWindow(),"Info", "in findRoutes function" )     # vl = qgis.utils.iface.mapCanvas().currentLayer()     vl=self.canvas.currentLayer()     director = QgsLineVectorLayerDirector( vl, -1, '', '', '', 3 )     properter = QgsDistanceArcProperter()     director.addProperter( properter )     # crs = qgis.utils.iface.mapCanvas().mapRenderer().destinationCrs()     crs = self.canvas.mapRenderer().destinationCrs()     builder = QgsGraphBuilder( crs )          global x1     global y1     global x2     global y2     pStart = QgsPoint( x1, y1 )     pStop = QgsPoint( x2, y2 )     graphAstar, nodes = self.make_graph({ "width": 8, "height": 8 })     paths = AStarGrid(graphAstar)          start, end = ??     path = paths.search(start, end)      My question is, how to pass the start and end coordinates to the function above? Because passing them just as coordinates (start, end = pStart, pStop) does not work. How do add them to the graph created as nodes? Or is there any easy way to do it? Please help me to to find a solution to this problem. Thank You


Cancel event in Python Gp tools with arcpy

In GP-tools written in C# (ArcObjects) we can listen for the Cancel-event and react on that. I wonder if there is anything similar in Python GP-tools applying arcpy. I googled that issue and the only thing I found was the Cancel-method on the Result-Object (so that you can ask the state of the result of a GP-tool, i.e. if it was canceled). I'm aware that you can cancel GP-tools running in the background in the Result Window, but this Cancel-event is catched by the system and therefore the developer or the GP-tool cannot react on that. Thank You in advance!


How do I access MPS Atlas from python in ArcGIS 10?

I'm having trouble finding documentation on accessing MPS Atlas features using python in ArcGIS 10. I need to be able to modify dataframe manager properties, swap map sheets, etc. Is it part of arcpy.mapping? From what I can see, there is no documentation of it on the ESRI Resource Center.


free GIS software that can orthorectify?

Would someone care to recommend a free and/or open source GIS application that can orthorectify map images?


Online/cloud storage options for GIS file-based data

I am currently looking into online/cloud-based storage solutions in general and also for my GIS data which consists of shapefiles, rasters, file geodatabases, and the like (no relational databases). Here are my requirements:   1. Unlimited storage (need several hundred gigs to start with)   2. Either unlimited file size or can handle files many GBs in size   3. Just need it for one Windows 7 box right now, but Mac support needed as well   4. Less than $100/yr - most are around $60 looks like   5. An Android app would be nice, not to back up the phone, but to access my backed up files from the phone   6. Backs up external hard drives   7. Can seed the backup with a hard drive and then get the restore sent to me on a HDD as well From what I have looked at so far, it looks like CrashPlan might be my best bet, except no mobile apps available. Mozy used to be unlimited, by apparently changed and now is not. Carbonite looks promising and has mobile apps, but doesn't backup removable hard drives and I can't find anything about if they let you seed your backup with a removable drive. Anyone currently using any of these services? Know of another service that might fit the bill? One service/provider per answer please so they can be voted on properly. Thanks.


ArcGIS10 Debugging

First ever post to this forum, I have a C# extension (legacy from another developer) which I am trying to debug in the ArcGIS 10 Environ (64bit machine/Windows7). I have built the ext in "Debug" mode and have a .dll, and a .tlb file available and ready to go in the "Debug" folder of my project. My User Account Control Settings are set on the lowest setting possible, I am wide open. This is a virtual machine. When I run debug, and then try to "Add from File" I get a message of "No Objects Added" from this project. Frustratingly enough, I have another extension which I wrote and I am able to fire it up and debug it. Try as I might I have not found any difference between the project set ups for either of these two extensions. Does anyone have any suggestions?


Provincial Base Mapping Data Sources in Canada

I'm very familiar with the Ontario Base Mapping Data, and am looking for similar data for other Canadian Provinces. Perhaps we can start a list of resources here? I'm looking for something at a larger scale than the standard NTS data. Ontario Base Mapping - 1:10,000/1:20,000 (depending on latitude) - available from Land Information Ontario (via an Unrestricted Use Licence) or from GeographyNetwork.ca


Is there a way to clip an image with the basic ArcView package (9.3)?

Is there a way to trim a geo-referenced air photo to an irregularly shaped boundary? I tried the clip tool shown in my attachment, but it didn't work. It looks as though the output extent field requires a rectangularly shaped polygon as the clipping feature. Can I even do this in ArcView? I've never had a request to do anything like this before. ![enter image description here](http://i.stack.imgur.com/Eac8R.jpg) * * * I do use Photoshop and other similar software packages. I could just go in to an environment like this & start painting everything outside of my boundary in white, but then I would have to come back in and geo-reference the image again, and it just strikes me that there must be a faster and easier way (even though I only have ArcView 9.3).


Get all Nodes from a List of adjacent LineStrings in PostGIS

I use PostGis and pgRouting. I have a table with geometries of the adjacent edges returned by a shortest path query. I want to concatenate these edges and retrieve a list of all nodes these edges consist of. How do I do this in SQL? I think the solution would be a mix of LineMerge and GeomUnion, but is there a smart way to do it? Thanks.


How to position contour labels in Maplex above the line?

I would like to utilize the **"contour placement"** options of **"Uphill alignment"** and **"Laddering"** while positioning the label above the contour not over it. This does not seem to be an option under the contour placement where it's forcing me to place the contour over the line. (Which makes no sense in ArcGIS since there is no easy way to mask that portion of the line efficiently. Possible but very cumbersome and not dynamic) I tried the other placement options in Maplex but then "Uphill Alingment" and "Laddering" is not longer an option. Does anyone has a good way of doing this? ![enter image description here](http://i.stack.imgur.com/m5AFc.jpg)


export openlayers as .shp file in Qgis

Is it possible to save an openlayers map as a shape or tab file in Qgis. I would like to have a part of an area from google maps as a shape file for later use?


Styles from Geomedia (mdb with style library) to be used in qgis?

We are trying to move some projects from Geomedia Professional to QGIS. Since we have large style libraries (stored as mdb) we are looking for a way to convert the styles to a format that QGIS can read. There is no export functionality in geomedia. Export in Access from mdb to xml failed when importing in stylemanager in QGIS. Can anyone help? (Windows 7; Geomedia Professional 13; QGIS 1.8.0)


Geocode address once in google maps, save in database, read geocoded result from there on

Is it possible to geocode an address, save the result in a database, then read that result from there on? I'm using a Google map that is geocoding hundreds of addresses which obviously hit the geocode limit. Once geocoded, the address is very unlikely to change again so it makes sense to save the geocode result in a database. I'm using Wordpress as a platform so if possible I'd like to save the geocoded result as a field in the database that relates to the particular company. **THE CLIENT HAS SPECIFICALLY ASKED FOR GOOGLE TO BE USED THROUGHOUT THE PROJECT THUS RULING OUT YAHOO OR ANYTHING ELSE.** UPDATE: Everyone's saying about what I'm asking is going against their TOS but isn't their very own example doing the same thing - developers.google.com/maps/articles/phpsqlgeocode?


Does SHAPE_Length take into account elevation values?

Does the `Shape_Length` attribute take into account the elevation? For instance, if I have 3 points: `A`, `B`, and `C`               A.X = 10     A.Y = 10     A.Z = 10     B.X = 10     B.Y = 10     B.Z = 50     C.X = 0     C.Y = 0     C.Z = 10      Then I create two lines: `L1`, and `L2`.               L1.fromPoint = A     L1.toPoint = C     L2.fromPoint = B     L2.toPoint = C      What I am observing in the results is that the `SHAPE_Length` attribute is the same for both `L1` and `L2`. I have verified that all the points and all the polylines are ZAware. So, am I wrong to assume the `SHAPE_Length` attribute takes into account the Z of the shapes? **UPDATE:** I have no problem storing the Z Values or manipulating them. I am using IPolyline. ~~**UPDATE 2:** Apparently if you have 3D Analyst the SHAPE_Length does include the Z value. If not, it is strictly XY. IF that is the case is there a way in .NET ArcObjects to update the SHAPE_Length field to include the Z value in the calculation? Is there anything I need to do or watch out for when manually updating this attribute?~~ **UPDATE 3:** From what I’ve been hearing, SHAPE_Length is 2D only and isn’t editable.


Workflow to convert a DWG CAD file to a shapefile

I'm trying to develop a workflow where I can convert a `.dwg` file to `.shp`. I don't have a good understanding of which tools I should be using. My preference is to use software that runs on `OSX`, is open-source (or free). My current workflow is as follows:   1. Convert `.dwg` to `.dxf` using Draftsight (http://www.3ds.com/products-services/draftsight/overview/ )   2. Load `.dxf` into `Sketchup`   3. Scale/rotate imagery in `Sketchup` and geolocate   4. View in `Google Earth` Once I reach step 4, I'm unsure how best to proceed. As far as I can tell, `Google Earth` is actually rendering a `Collada` model, and not loading a `kml` file that I can access, so I don't know how to extract the underlying file. I thought about using `ogr2ogr` (http://www.gdal.org/ogr2ogr.html ) but I'm not sure if I can interact with a `CAD` file this way (as I need to specify the projection and location). I also have `QGIS` and `PostGIS` available. I also considered `FME` but would prefer not to have to have a Windows OS just for this conversion process. I've also read through the following related question: Ways to transform from DWG/DXF to GIS (Experiences are welcome)


How to obtain the maximum/minimum longitude and latitude for a polygon?

I have the MSA/CMSA shape file which contains 274 metropolitan areas. How do I have to obtain the maximum/minimum longitude and latitude **for an area**? So, it's like drawing a box right outside and getting the 4 corner points in the lon/lat format. Can we do this using Python to iterate the process for all the cities and to export the coordinates to a text file as follows?:               msa     min_lon   max_lon   min_lat  max_lat     6162   -76.XXXX  -74.XXXX   38.XXXX  40.XXXX     9094   -97.XXXX  -96.XXXX   37.XXXX  38.XXXX     ....     ....      Thanks!


Problem mixing the google maps and WMS layer

I have been trying to use both google maps and WMS layers as an option in the layerswitching panel. But when I switch the base layer from WMS to Google Maps layer and then back to WMS, _the bounds of the map do not come right and the map is not displayed properly. Even the vector overlay features which display correctly at 1st in the WMS, get crowded into small area on the google maps. I think there is some problem with the re-instantiation of the maps variable during layer switch, due to which the map bounds are not correct._ But still on the whole, please tell me how to solve this problem as it would be great to be able to use any map service in a single HTML page.


Can't import arcpy in IDLE, pyscripter, or even ArcGIS10.1's Python window

After using arcpy and python for months with ArcGIS 10.1 without any problems I'm encountering something I can't work around or even find help with online. Something as simple as "import arcpy" is failing for no reason I can think of:               >>> import arcpy     Traceback (most recent call last):       File "<interactive input>", line 1, in <module>       File "C:\Program Files\ArcGIS\Desktop10.1\arcpy\arcpy\__init__.py", line 24, in <module>         from arcpy.toolbox import *       File "C:\Program Files\ArcGIS\Desktop10.1\arcpy\arcpy\toolbox.py", line 342, in <module>         from management import Graph, GraphTemplate       File "C:\Program Files\ArcGIS\Desktop10.1\arcpy\arcpy\management.py", line 153     SyntaxError: Non-ASCII character '\xe2' in file C:\Program Files\ArcGIS\Desktop10.1\arcpy\arcpy\management.py on line 154, but no encoding declared; see http://www.python.org/peps/pep-0263.html for details      I've looked around and I've seen other people with similar symptomes but nothing with the same last message, the `Non-ASCII character '\xe2' in file C:\Program Files\ArcGIS\Desktop10.1\arcpy\arcpy\management.py on line 154, but no encoding declared;` part seems different. I haven't done anything like installing new python extensions or versions or any other type of software. The only thing remotely related was that I had an FME process working in the background while I was testing my code. I didn't modify any pre-existing ArcGIS python scripts. I have tried using IDLE installed by ArcGIS 10.1, pyscripter, and the Python window that can be opened inside of ArcGIS. I tried rebooting on the off chance it was some fluke. The only thing I can think of to try next is a fresh instal, and if that doesn't work maybe going as extreme as wiping the PC's OS and starting over from scratch. Any help would be greatly appreciated.


How do you use 'reclass' in GRASS?

It could be user error, but for some reason, the reclass tool is just hanging and hanging... QGIS 1.8 Compiled against GDAL/OGR 1.9.2 Ubuntu 12.04.1   * Running the tool from the Sextante toolbox   * Input is a 10' cellsize raster ~10mi x 10mi in area   * grass region cellsize = 10   * Reclass file is in order from what I can gather from the help When run: - There is no process operating in the system monitor - Dialogue hangs and cursor shows progress spinny thing... Not sure what else to say? Does the GRASS tool run from this dialogue on a .tif? Or does it have to be run from the GRASS tools specifically? Thanks!


"SDE error 6" when publishing GeoServer layer

I have successfully created an ArcSDE datastore in GeoServer 2.2.4. When I try to publish a feature class in the SDE database, the "Compute from data" link returns nothing, and the "Compute from native bounds" link produces the following error.               Error computing the geographic bounds:[SDE error 6][Error desc=UNKNOWN ERROR CODE. ][Extended desc= ]      The Native SRS says UNKNOWN in the text box but has the correct system to the right "GCS_WGS_1984...". Because I do not have either of the bounding boxes filled in, obviously I can not save the layer. I have tried feature classes with other coordinate systems, deleting and re- creating the datastore and creating whole new feature classes. Extensive searching on the internet has not turned up any kind of solution. Any help would be great! Thanks!


edit 10.1 feature service

has anyone edited a 10.1 feat. service using desktop 10? I can view the service and add it to the map, but when I go to edit features/create local copy for editing, I get an "error creating the..." I perform the same exact steps on desktop 10.1 and it works fine. I have a bug report submitted, just wondering if anyone else experienced this. thanks!


Create a shapefile from geometry

Please excuse my lingo, I am still learning some aspects of geo spatial and arcpy. I have a python script that takes a shapefile. The shapefile consists of a polyline. I have traversed the polyline, done some interpolation and retrieved some points. Long story short, I'm trying to generate polygons along the polyline based on intervals of the M values. These polygons have four points (my geometries). How can I take these geometries, or points, that are in my python script and output them into a shapefile. I know I'll have to pass my points to a function that produces a shapefile, but that is my dilemma. How exaclty do I create the function that creates my shapefile for the geometries?


FME - Possible to return a value from workspace?

I'm calling my workspace from the command line (either via Python or .NET to a system/os call). I'm parsing DWGs and pushing attributes to SQL Server in the workspace. Either the data validates in the workspace or it doesn't. I need to return a value of 'pass/fail' to my calling application from my workspace based on whether or not the data validated within the workspace. Is this possible?


How to create a table with one polygon and one point?

using postgresql and postgis, i would like to create a table with one polygon and one point. then use st_contains to check if the point is inside the polygon or not. so first, I should create a table :               CREATE TABLE my_schema.my_table (id serial);      then add               SELECT AddGeometryColumn( 'my_table', 'geo_polygon', 4326, 'POLYGON', 2);     UPDATE my_table       SET geo_polygon=GeomFromText('POLYGON((0 0,4 0,4 4,0 4,0 0))',4326)      WHERE id='p1';               SELECT AddGeometryColumn( 'my_table', 'geo_point', 4326, 'POINT', 2);     UPDATE my_table       SET geo_point=GeomFromText('POINT((2 2))',4326)      WHERE id='p2';      Finally, to check containing i should use :               SELECT ST_Contains(geo_polygon, geo_ponts, 4326))          FROM my_table         WHERE ???      I am not sure about the condition here. do I have to use the ids ? and how ?