ArcGIS Online - My Hosted Services

I am looking for a solution to a problem I have using ArcGIS Online. I am currently creating 5 Tiled Map Packages which I tile on my ArcGIS for Desktop then Publish to my ArcGIS Online Account to save using any credits. Once the Tiled Map Package is Published I can then email my colleague the URL (REST) which he then feeds into his ArcGIS Server. He has asked me if I can group the Services together into one URL so that he can display the 5 Published Services as a Group instead of individual Layers. Is this possible? Any help would be greatly appreciated!! Thanks Ryan


Simplifying adjacent polygons using QGIS "Simplify Geometries" tool

I have a big Shape file with 200,000 polygons but the file is too big for the application that will use it. I would like to generalise these polygons and create much smaller file. I have tried to do this in QGIS using "Simplify Geometries" tool. I have tried this on a very small portion of my Shape file and I am happy with the results but have small problem which I do not know how to solve. I have appended the image to show the problem. WHITE boundaries are boundaries before simplification. RED boundaries are after simplification. The problem areas are circled with orange colour. As you can see "simplified boundaries" contain some areas that after simplification do not belong to any of the previous polygons. How do I simplify this Shape file and keep common boundaries between polygons as a single line and eliminate areas that have been created after simplification? Thanks in advance for any help. ![Simplified Polygons in QGIS](http://i.stack.imgur.com/S6fPQ.png) ADDITIONAL INFORMATION: After following suggested instructions   1. Convert polygons to lines    2. Simplify lines   3. Convert lines back to polygons  I get the following results. It is slightly better than before but still I get areas outside original polygons (see white areas in the picture below). Andrew suggested that I associate each line with adjacent polygon but I do not know how to do it in QGIS. Thanks in advance for any help. ![enter image description here](http://i.stack.imgur.com/AKY9s.png) David Hi guys, I am sorry to bother you with this issue but I am still trying to complete simplifying polygons process. I have followed instructions suggested in this post and when I thought I got acceptable results one more problem appeared. Below is the summary of what I did with images of results. All steps have been performed using PROCESSING toolbox in QGIS.   1. First, I converted polygons to lines   2. Second, I cleaned lines using v.clean with break option. This resulted in CLEANED VECTOR LAYER and ERROR LAYER. Error Layer (green dots on the provided image) has green dot on every line intersection. I do not know what this means. Any input on the error layer would be appreciated.   3. Third, I simplified lines using "SIMPLIFY GEOMETRY" tool (result in red lines in the images). Comparing simplified red lines and original orange lines I am happy with the result.   4. Fourth, I wanted to convert lines back to polygons using "Polygonize" tool in "PROCESSING" and got results that I did not expect (second image with a shaded dark olive colour polygons). It looks like polygonizer connected error layer intersections and created completely different set of polygons ignoring simplified lines. Could someone explain what went wrong here? Why red simplified lines were not converted properly to simplified polygons? As you can see I am not an GIS expert but need to do this myself. I do not know if I am too far from the solution but it is frustrating that I can see simplified lines that I found acceptable for what I'd like to do but cannot convert them to polygons. Thanks in advance for any help! Here are two images with results of last two steps of the process: ![Simplified Lines and Original Lines](http://i.stack.imgur.com/xj0LN.jpg) ![Results of Plygonized simplified lines](http://i.stack.imgur.com/oBTCb.jpg)


How to make polygon volume in TIN surface?

I want to measure surface volume of an area in different elevation layer. For measuring the surface volume I have to make a Polygon Volume as well. I have problem with making this polygon over the TIN layer and define the different elevation for it ? Like the link below: http://resources.arcgis.com/en/help/main/10.1/index.html#//00q900000037000000 for example I have a TIN surface which min contour is 25 m and max contour is 140 , then I want to measure all the surface above 25 m , 35 m and ... ! I just wonder how to make this polygon on top?


OSM - Extract Streets from OSM Data (With/Without Python)

I'm currently in the works of creating a pythonic way of extracting street names from an OSM dump for the Province of Quebec. Originally, I wanted to do this in Python to make it easier to integrate the process in future code, but I'm failing miserably, and after weeks of on and off research without getting anywhere, so now I'll honestly take what I can get. What's the best way to extract all the street names from a dump such as these? http://download.geofabrik.de/north-america/canada/quebec.html Here's my original code if anyone can make it work. At the moment this code maxs my 16GB of ram and 32GB of swap in about a minute before it evens begins going through each element.               from lxml import etree     import xmltodict, sys, gc     from pymongo import MongoClient          #Ultimate fix of life.     reload(sys)     sys.setdefaultencoding("utf-8")          client = MongoClient()     db = client.re     streetsDB = db.streets          hwTypes = ['motorway', 'trunk', 'primary', 'secondary', 'tertiary', 'pedestrian', 'unclassified', 'service']          #Enable Garbadge Collection     gc.enable()          def fast_iter(context, func):         print 'Begun processing'         placement = 1         # http://www.ibm.com/developerworks/xml/library/x-hiperfparse/         # Author: Liza Daly         for event, elem in context:             placement += 1             print placement             func(elem)             elem.clear()             while elem.getprevious() is not None:                 del elem.getparent()[0]         del context          def process_element(elem):         data = etree.tostring(elem)         data = xmltodict.parse(data)         keys = data['way'].keys()         if 'tag' in keys:             if isinstance(data['way']['tag'], dict):                         if data['way']['tag']['@k'] == 'highway':                             if data['way']['tag']['@v'] in hwTypes:                                 streetsDB.insert(data)             else:                 for y in data['way']['tag']:                     if y['@k'] == 'highway':                         if y['@v'] in hwTypes:                             streetsDB.insert(data)                             break         del data         del keys         gc.collect()          context = etree.iterparse('quebec-latest.osm', tag='way' )     fast_iter(context,process_element)


How to create the custom tool in ArcGIS desktop with using of python script from model bulder?

I created some model builder for some customization to reduce time in our routine process. I exported my model builder to python script by availability of export options. It doesn't work if i bring to tool bar. any one suggest me clearly how to access the script in tool bar with some standard clear steps to follow? as well as whether i need to do any modification on my exported python script?


VBA functions use in ArcGIS Field Calculator

Would this be ok to use them for the field calculator in ArcGIS 9.3 ? See here MS Excel VBA Functions Or should I fix them in MS Excel and save as a file and then use the join with the shapefile or file geodatabase. Will this work ? Will ArcGIS read it in table ? Please advice as I am looking for information on how to use field calculator. Thanks


How to use nested loops with cursors?

I use the following nested loops to read the first point feature of a dataset and the rest points in that dataset. Then read the second point feature of a dataset and the rest points in that dataset and so on. The problem with these nested loops is that it reads the first point and the rest points only once. Why this happens?               Dim pFCursor As IFeatureCursor     Set pFCursor = pFeatureClass.Search(Nothing, True)          Dim pFCursor2 As IFeatureCursor     Set pFCursor2 = pFeatureClass.Search(Nothing, True)               Dim pFeature As IFeature     Dim pFeature2 As IFeature          Dim pPointA As IPoint     Dim pPointB As IPoint               Set pFeature = pFCursor.NextFeature     Do Until pFeature Is Nothing          Set pPointA = pFeature.Shape          MsgBox pFeature.Value(intPosFID) & "F1"             Set pFeature2 = pFCursor2.NextFeature        Do Until pFeature2 Is Nothing                        Set pPointB = pFeature2.Shape                   MsgBox pFeature2.Value(intPosFID) & "F2"                  Set pFeature2 = pFCursor2.NextFeature        Loop                Set pFeature = pFCursor.NextFeature      Loop


How to create varying opacity within a layer in Mapserver?

Suppose you have a mapfile containing a background layer and a layer with a lot of overlapping linestrings. Now, for the overlapping linestrings, it might be good to visualize their density, thus at places where a lot of lines do overlap they might be drawn more "solid" than at places where only a few lines are present. Using QGIS, for example, I can do this by just turning the layer opacity to some small value (say five percent). In Mapserver, however, setting the opacity of the layer, the opacity will just make the whole layer see through, but will not add up for each single element. How can I archive the desired effect?


What type/model of pointing device do you use for GIS?

I've been using basic mice for a long time but have been thinking of moving either to an ergonomic mouse or trackball. Trackballs seem like they may be better for me, but it seems like manufacturers have neglected the trackball market for a long time and there aren't that many good options available. Something like the Evoluent vertical mouse looks interesting as well. So what exotic pointing device do you use, and why do you like it or not like it? * * * I ended up going the trackball route and chose the Logitech Cordless Optical Trackman: ![Optical Trackman](http://i.stack.imgur.com/rq9Gr.png) It was the closest thing I could find to the physical form factor of the Microsoft Intellimouse Trackball Explorer, which was the best reviewed trackball I found in my search, but is no longer manufactured (though for a few hundred bucks you can get old ones on Ebay): ![Old MS trackball](http://i.stack.imgur.com/NCgEK.png) I would have preferred a USB wired version but the only one Logitech makes in this form is cordless, with a huge wireless USB dongle. Actually, this trackball doesn't even show up on Logitech's site anymore, I wonder if they are discontinuing it. It's available at Office Depot and Staples as I write this, I've had mine for several months now. The **reasons I chose the trackball** , and this one in particular, were:   1. Similar physical layout to the Microsoft model I mentioned. That MS trackball has legions of rabid fans online and they will spend hundreds of dollars to replace their dead Intellimouse Trackball with a new one. I didn't find any currently-available trackball with as dedicated a fan base.   2. It requires no more desk space than the unit itself occupies, unlike a mouse which needs a mousepad. This also means that if I'm digitizing something big, I never run out of mousepad and have to pick up a mouse, move it, and put it back down again, which can be a pain and cause errors.   3. Very little force is required to get the ball/pointer moving, much less than getting a mouse to move on a mousepad. It makes fine pointer movements a lot easier.   4. This particular one has the ball mounted on the top, unlike the Logitech jvangeld has, or the M570, another Logitech trackball. Both of those require you to use your thumb to move the ball, which seemed to me like it would get more tiring over time and also that finer movements would be harder than with the ball on top (where I can use 3 fingers to move it).   5. It has 5 extra buttons beyond Left/Middle/Right mouse buttons, those can be programmed to do whatever using the included software (or AutoHotKey). **A few Cons** I've noticed:   1. The small pegs that keep the ball centered in its socket need to be cleaned regularly. I've been doing it only when I noticed the wheel not rolling as freely but it should probably be done once a week at least. It seems like those pegs are going to wear out long before anything else does.   2. The mouse wheel's scrolling action is a bit stiff, so is the left mouse button.   3. It's apparently an old model and the USB wireless dongle is very big: several feet of cord with a dongle about 3" square on the end of it. It's not travel-friendly.


How can I find out how to convert this point to lat lon?

I have received a shapefile with measurement points to be analyzed, but no matter how I try, I can't find how to convert to lat lon. I suppose it would be good to know the projection but when I look into the file (with R and summary()) it says "Is projected: NA" and "proj4string : [NA]". One measurement-point I know is close to the border Peru-Chile has x,y = 354521, 7997417.8. How can I find out how to convert this point (and many others) to lat lon? I tried to follow the instructions of similar questions here, but to no avail. Thanks, /Chris


How to render GeoJSON tiles from Tilestream for polymaps?

Have been using polymaps successfully till now. Need to generate GeoJSON tiles for Canada Provinces. Found Tilestream and got it working with Tilemill. Is there a way I can make Tilestream hand out GeoJSON tiles for polymaps without Wax? http://my_tilestream_server:8888/provinces/{Z}/{X}/{Y}.json ## EDIT Any help here would be greatly appreciated. Still trying to make either Tilestream / Tilemill or Tilestache work correctly. Does anyone know where I can find a good shapefile for Canada and its Provinces (Canadian Fips code would be a bonus)? Only ones I have found are not accurate at zoom levels greater than 8. 1 - Tilestream / Tilemill are great, but need a way to dynamically add additional data to the Hover. Can't seem to find a way to do this without including it in the dbf. 2 - Installed Tilestache, but having problems converting data from shapefile to sql. I think I am messing missing something in the conversion to sql. Or is there a data source I can import that already has this info?


ArcGIS - How do I group corner-touching contiguous polygons?

Using ArcGIS 10, I have a Feature Class which I converted from a Raster (GRID) (Raster To Polygon Tool).   * There are approx. 180,000 polygons.    * The attribute table contains a [grid_code] field which I need to preserve.   * If one or more polygons are contiguous with corners touching AND they have the same [grid_code] value, then I'd like to group them (a dissolve basically) into one polygon (record).   * Then I need to calculate the grouped polygon's area. Using the Select By Location Tool, I'm able to successfully select the polygons that I need to group, but how do I then convert them into ONE polygon, while preserving the [grid_code] ?? Thanks.


How to check if a shapefile is empty with arcpy

I have a bunch of shapefiles in a folder. But some of them are empty, i.e. there are no rows within the shapefile. Is there a simply to find those empty shapefile out?


Too many features selected with the ArcMap Select by Polygon tool

Here is my problem: First of all you have to know that I have ArcMap 10.1 SP1. My features are stored in Oracle 11g database in SDO_GEOMETRY format. I have imported a layer of accidents location in ArcMap. So when I try to select all features included in a polygon (using the Select by Polygon), then too many features are selected. I mean some feature outside my polygon are part of selected feature. So I draw a triangle on my map and I expect that all features included into the triangle will be selected. The result I obtain is all features included in the rectangle envelope of the triangle that are selected, which is obvious wrong. Is someone can help me on this one? You have to know that I was in SDE_GEOMETRY before and I wasn't had the problem when I tried to select by polygon. Thank you. ![Here is my polygon selection \(triangle\)](http://i.stack.imgur.com/dsvi4.png) ![This is the wrong result of the polygon selection](http://i.stack.imgur.com/QUuzh.png) ![The rectangle envelope that seems to be used by ArcMap](http://i.stack.imgur.com/nWO6E.png)


Is AMEOBA ArcToolbox compatible with ArcGIS for Desktop 10.2?

Is the **AMEOBA ArcToolbox** compatible with ArcGIS 10.2? My error (when trying to run one of its tools) is noted below. ![enter image description here](http://i.stack.imgur.com/hvfcG.jpg) (click to enlarge)


How to display lines with Arrows on a FeatureLayer?

I've setup a feature service with a symbology. One symbol is a line with an arrow at the end. The description for such a symbol is not part of the FeatureServer response data. So it is not supported? I only see this (not sure why the width is 6):               UniqueValueInfos:     Value: Arrow      Label: Arrow      Description:      Symbol:     Style: esriSLSSolid      Color: [0, 0, 0, 255]      Width: 6      If not supported, any tips on how I might proceed? Perhaps a termination icon attribute and then some custom gfx code? Anything already exist for this?


sql query to get origins and destinations

I have a table of cities with an ID and geometry (200K rows) and a table of travels, with only two ID: idfrom, idto (7K rows) I'm wondering how could I select the geometries for each one in order to make a line later. What I'm finding is something like:               idfrom, geom_from, idto, geom_to      So, I can use the ST_makeline postgis function. I've tried with this SQL expression:               select r1.idfrom, city.geom as from_geom, r2.idto,ity.geom as to_geom     from routes as r1, routes as r2, city     where r1.idfrom=city.id     and r2.idto=city.id      But it seems that it does not work, It returns the same values in 'from' and in 'to' columns. I've tried more queries, but still have not found the solution Any clue or suggestion? Thanks!


How to add Attachment to an object in ArcGIS Online?

I created attachments for some objects of a theme (as a File Geodatabase Feature Class) in ArcGIS 10.0: images (JPG) and doc files. Then moved the same shp theme in my web-map of ArcGIS Online, and want to add in any way these attachments for these objects there as well. So, my main goal - find a way to see this attached jpg and doc files when I choose concrete object in the web-map of ArcGIS Online. Also I’d like to find opportunity to place in the web-map of ArcGIS Online - URL in attribute table for any concrete object. Could you advise me please – how it’s possible?


How to make a predetermined Polygon (Circle) From 4 Coordinates, The Radius and Center Longitude and Latitude in Openlayers

I know how to make a predetermined Polygon in Openlayers (see below). Please have a look at http://studio.netnrg.co.za/dev/rnd/blake.loizides/openlayers                     var feature2 = new OpenLayers.Feature.Vector(           OpenLayers.Geometry.fromWKT("POLYGON((-120.828125 -50.3515625, -80.1875 -80.0078125, -40.40625 -20.4140625, -120.828125 -50.3515625))"                 )             );           vectors.addFeatures([feature2]);      I have made a Edit Circle (40 points in Openlayers) (My Client wanted to see a circle on the map for zones) The problem I have is that I can only store the following details and not the 40 points or coordinates or Vertices to make a circle in Openlayers. I can only **STORE 4 COORDINATES** In the Database I store the Bounds Top Bottom Left and Right (bottom: -29.054153382197, left: 23.782779686138, right: 23.79286499053, top: -29.044068077805) I will replace the coordinates with the ones in my database. The Radius (525.627982303084) And the Center Longitude and Latitude (lon=23.787822338334,lat=-29.049110730001) Now I made some calculations below. My Question is how do I now recreate a (circle) with the information i have above that is stored in my database? I know how to get data from my database and display it on the map, What I don't know is how to make a circle if I cant store the 40 points and only have the Radius the Bounds and the Center Lon Lat to make a circle . Please if anyone can help me I would be greatly appreciative. I only want a Circle for the look and for no other reason. ![My Calculations](http://i.stack.imgur.com/P9GvP.jpg) ![Circle from Openlayers](http://i.stack.imgur.com/EGy4K.png)


How to merge different classes/clusters in one single raster in arcgis 9.3?

I've created a signature file with 8 clusters and performed an unsupervised classification. Later I have grouped values (spatial) classes together to form 4 classes I needed (in symbology) but that just "colors" them the same, they're still treated as separate categories. I need to merge them to have only 4 groups/categories. What to do?


Opening .shp file under Add Vector Layer

When I try to open streets.ship file from the exercise data file downloaded for working in the Quantum GIS Training Manual, I get an error message: Invalid Data Source C:/Users/new user/Documents/streets.ship is not a valid or recognized data source. What am I doing wrong?


Qgis & processing doesn't links with Grass, r or saga

After having formatting my computer (W7-64bits), I have re-installed QGIS 2.0.1-64 bits on it . But I can not use GRASS, SAGA or R via Toolbox (or Processing). I have checked the path in the Options, and I have tried to uninstall and re-install it again, but it doesn’t work yet. In the beginning of the re-install, I have a message telling that I QGIS is already installed, but it didn't remain any folder or files link with Qgis. I checked as well in the option and plugin setting, but everything seem to be fine (plugin activate, plugin repositories connected...). Did I install something else than QGIS? did I miss a step, etc.? Many thanks for help!


Linking Google Sheets with PostGIS

I am curious if there is a way to link a PostGIS table with a Google Sheets? Both tables (Post and Google) have a unique ID that can be used as a linkage. I am hoping to be able to build a view that will display information from both tables into a single view. Is this something that can be done, if so, any pointers would be appreciated. Thanks...


Creating many origin-destination routes with pgRouting

I'm using a PostGIS / pgRouting database setup with a road network. I'm trying to create an origin-destination table with the complete route information, so not just the sum of the length, like in this post. I'm trying to get something like this as output table:               Origin | Destination | vertex_id | cost     1        2             415         14.2     1        2             508         4.3     1        3             919         2.4     1        3             1024        6.8      The following code for calculating a single route works fine, but how do I automate this for a lot of origin-destination pairs?               SELECT * FROM shortest_path('     SELECT gid AS id,     start_id::int4 AS source,     end_id::int4 AS target,     length::float8 AS cost     FROM network',     1, -- start_id     2, -- end_id     false, false);      I've tried something with subqueries, like this where I'm getting the start_ids from table `nodes`:               SELECT id,     (SELECT * FROM shortest_path('     SELECT gid AS id,     start_id::int4 AS source,     end_id::int4 AS target,     length::float8 AS cost     FROM network',     id::int4, -- start_id     450969, -- end_id     false, false))     FROM nodes;      But this doesn't work because the subquery returns multiple rows (as expected): `ERROR: subquery must return only one column`. Any help?


How do I resolve error DataGridView/Label1 is not declared?

' set the active display Dim disp As ESRI.ArcGISExplorer.Mapping.MapDisplay = ESRI.ArcGISExplorer.Application.Application.ActiveMapDisplay                   ' define the item that selected by user         Dim SelItem As ESRI.ArcGISExplorer.Application.SelectedItemsCollection = ESRI.ArcGISExplorer.Application.Application.SelectedItems              ' define the name of selected item         Dim SIName As String              ' define a featurelayer for selected item         Dim FLitem As FeatureLayer              ' get the name of user selected item         SIName = SelItem.Item(0).ToString()              'display the name on the label         Label1.Text = SIName              ' find the featurelayer that has the name which equals to the one user selected         FLitem = disp.Map.FindByName(SIName)              ' create a new bindtable and fill it with featurelayer's table         Dim FLbind As TableBindingAdapter         FLbind = New TableBindingAdapter(FLitem.Table)         FLbind.Fill()              ' using the datagridview to show the bindtable which is the selected rows         DataGridView1.DataSource = FLbind


How to use overview map with atlas plugin in QGIS 2

I'm using atlas plugin to batch generate maps. I insert an overview map but it doesn't get updated over iteration and generation of maps.Is it possible? If it's not, is there any alternative?


ArcPy replaceDataSource

I'm attempting to replace a datasource in an outline map frame using ArcPy. The rest of the script seems to run OK, but after the script has run the datasource the datasource of that layer shows an error. A section of my code is below:               for df in arcpy.mapping.ListDataFrames(mxd):         if df.name == 'Overview' or df.name == 'Keyplan':             df.scale = 500000         for lyr in arcpy.mapping.ListLayers(mxd,"",df):             if lyr.datasetName == 'blah':                 arcpy.mapping.RemoveLayer(df,lyr)                 lyr.replaceDataSource(r'\\directory','FILEGDB_WORKSPACE','FCName')                 arcpy.mapping.AddLayer(df,lyr,"TOP")          arcpy.RefreshTOC()     arcpy.RefreshActiveView()     fullpath = os.path.join(fileDir,'ALTERED_' + fileName)     mxd.saveACopy(fullpath)     del mxd      I have had a look at http://forums.arcgis.com/threads/23803-issues-updating- data-source-via-arcpy.mapping?p=78411&viewfull=1 and used some of the tricks there, but I can't seem to make it work. Picture of layer datasource before: ![enter image description here](http://i.stack.imgur.com/oZVDt.png) and after ![enter image description here](http://i.stack.imgur.com/Xbzoj.png) You'll see it has the correct layer name, but no actual layer. Any ideas what could be causing this. I thought at one point relative paths, but apparently they do not matter when using UNC paths, as I am. This appears to work when I do it in ArcGIS itself using the "CURRENT" map, so pretty sure it's not a pathname error. Any ideas?


Estimating sunlight hours from slope, aspect, and latitude

I'd like to estimate annual sun hours (or insolation, or solar radiation) for 4 points, 1 each on the four sides of the summit of a hypothetical flat-sided mountain. I have slope and aspect for each point, and latitude of the mountain, and for the purposes of the question, am not considering other topographical features or weather. I know that there are tools in ArcGIS and GRASS to do this. However, it seemed to me that looking at annual hours (rather than having to deal with seasonality), and ignoring other topography should make this problem simpler. In fact, would it even be necessary to consider latitude, if I'm after annual totals? So, is this something that it's possible to estimate with a formula? Or must I open up ArcGIS and load a DEM to do it?


Setting time animation by number of frames

I'm trying to use the ArcGIS 10 animation system to produce a time series of points by year. There are only 65 years in the dataset, and I want to make it so there is only one frame per year. I have tried setting it to play a certain number of frames and frame duration, or by total duration, but when it exports, it alters the frame rate and adds a number of duplicate frames when I export to AVI or as sequential images. Is there a way to export only a single frame for each year? The interval and units are set to 1 year.


Setting up MapServer WMS and use it with OpenLayers

I have set up MapServer WMS on the localhost by the following map file:               MAP       NAME "Prvi_WMS"       SIZE 600 300              # 600 by 300 pixel image output       PROJECTION            "init=epsg:4326"       END           #EXTENT -180 -90 180 90         WEB         IMAGEPATH "/ms4w/tmp/ms_tmp/"         IMAGEURL "/ms_tmp/"         METADATA           "wms_title" "My Global Map WMS Server"           "wms_onlineresource" "http://localhost/cgi-bin/mapserv.exe?mode=map&map=c:/ms4w/Apache/htdocs/MapFile06_wms.map&"           "wms_enable_request" "*"           "wms_srs" "EPSG:4326"         END       END            LAYER              NAME countries         TYPE POLYGON         STATUS DEFAULT         DATA "/ms4w/data/countries_simpl.shp"         MINSCALE 1000         MAXSCALE 1000000000              PROJECTION            "init=epsg:4326"         END           EXTENT -180 -90 180 90              CLASSITEM 'NAME'         CLASS           NAME 'Bulgaria'           EXPRESSION 'Bulgaria'           OUTLINECOLOR 100 100 100           COLOR 255 255 150         END              CLASS           NAME 'All Countries'           EXPRESSION ("[NAME]" ne 'Bulgaria')           STYLE             OUTLINECOLOR 100 100 100           END         END         METADATA           "wms_title" "World Boundaries"           "wms_onlineresource" "http://localhost/cgi-bin/mapserv.exe?"            "wms_srs" "EPSG:4326"           "wms_enable_request" "*"         END       END     END      then, I checked GetCapabilities and GetMap requests and the results seem fine. However, I have problems getting the map inside the OpenLayers. Here is how I define my wms layer:                             <title>Karta</title>         <link rel="stylesheet" href="openlayers/theme/default/style.css" type="text/css">
Primer prekrivanja slojev in izbire podlag

           The wms layer shown in OpenLayers is completely white. What am I missing? Thank you!


Is there a difference between the Planchon & Darboux and Wang and Liu depression filling algorithms? Other than speed?

Can anyone tell me, based on actual analytical experience, if there is a difference between these two depression filling algorithms, other than the speed at which they process and fill depressions (sinks) in DEMs? A fast, simple and versatile algorithm to fill the depressions of digital elevation models _Olivier Planchon, Frederic Darboux_ and An efficient method for identifying and filling surface depressions in digital elevation models for hydrologic analysis and modelling _Wang and Liu_ Thanks.


Purely client-side geocoding options for a webmap

I am creating a simple webmap that should be able to:   1. Display external maptiles (any of the following: Cloudmade, Mapquest, Google).   2. Overlay additional static layer in Geojson or KML format, that can show a popup box on a mouse-click.   3. Use a search box for geocoding that is limited only to towns and cities. And then zoom in to the geocoded location at a specific zoom level (z8-12 in webmap speak). I'm looking at several webmapping framework like Openlayers, Leaflet, Polymaps. So far, I was able to do 1 and 2. I'm having difficulty integrating the third feature. There are several geocoding APIs available (OSM's Nominatim, Geonames and Google's geocoder). However, I can't find any that limits geocoding to towns/cities. A major challenge is to do this purely on the browser/javascript, no server side processing since the webserver is only limited to static htmls.


Event listeners after each plotting of a point of a Polygon Vector in OpenLayers

I am quite new to OpenLayers. I have a polygonLayer for users to draw polygon vectors. The plotting feature is working fine.               polygonLayer = new OpenLayers.Layer.Vector("Polygon Layer");     map.addLayers([ polygonLayer ]);     var polygonControl = new OpenLayers.Control.DrawFeature(polygonLayer, OpenLayers.Handler.Polygon);     map.addControl(polygonControl);     polygonControl.activate();     var editControl = new OpenLayers.Control.ModifyFeature(polygonLayer);     map.addControl(editControl);     editControl.activate(RESIZE);      What I want to do is to have event listeners after each plotting of a point. Supposedly, I want to show popups to assist the user in between each point of plotting. How could I achieve it?


How can we improve the OpenLayers documentation?

I have difficulty understanding the documentation. Is it just me? Are there alternative sources of information?


What does "Size/Width" stand for in the symbol selector?

Can anyone explain how to interpret a symbol size of 10 or let's say 100? What kind of units are these? I ask because recently I had to deal with a reference scale and I wasn't really able to explain what symbol size means.


OpenLayers plugin for QGis problem with maps loading

I am new in GIS. I use QGis 2.0.1 with OpenLayers plugin downloaded from official repository. When I click to Add OpenStreetMap layer the map is not showing and load process is neverending. What is the problem and how can I solve it ? Thanks for your answers.


How to convert orthophotomap to .ecw file?

I am looking for a method for convert orthophotomap (.adf files) to ERDAS .ecw file. Have you got any ideas ? Or maybe you have got any other idea how i can reduce a size of orthophotomap .adf files? Now i have got 16 GB file I want to compress it to something like 3-4GB- is it even possible? I read that .ecw is highly compressed for raster formats...


QGIS connecting to ArcGIS world imagery - not loading correctly

Have connected QGIS to ArcGIS World Imagery using python console as described on http://www.northrivergeographic.com/qgis-adding-an-arcgis-rest-service. QGIS Seems to connect and make a raster layer the image is not correct. I'm very new to QGIS so would appreciate any ideas on how to resolve this problem. Thanks in advance ![Problem with World Imagery](http://i.stack.imgur.com/sX84q.png)


Mosaic tool not licensed?

I am very new to arcgis in general, and am having trouble stiching several rasters together. I know that the mosaic tool can do this, but when I try to create a new mosaic dataset I get the error message "you do not have the necessary license to execute the selected tool". The spatial anaylist extension is turned ON, and I have checked it several times, and it is otherwise working. Is there something else besides the spatial analyist I need? I couldn't figure out how to post on Esri's forum, and hoped someone here could help me. I am getting desperate! Any help is appreciated. Thank you!


Create features from ArcGIS Topology errors

I have imported field users shapefiles into an ESRI geodatabase so that I can create some topology rulesets that identifies the errors as needed. Unfortunately the field users do not have a license that can read gdb topology rules (or obviously fix them). I would like to be able to generate a new set of shapefiles that capture where all of the errors and exceptions are at so the users can go in an fix the obvious ones themselves (saving me a bunch of work!). This can also be a quick way of organizing the errors based on size (area or length) so that you deal with the biggest issues first. Does anyone have any suggestions on a way to generate features from ArcGIS Topology errors & exceptions? I have an ArcInfo license available so any tool within the package is at my disposal. While I am not an ArcObjects programmer, I am very familiar with Python scripts or ModelBuilder to develop such a tool/process. **UPDATE** : I did just find this ArcGIS Idea and developer sample, but I was hoping for something that I would be able to put into a toolbox, not a custom VB module.


Dimensional analysis/conversion problem: Map scale to image resolution (dpi)

Given a map scale of 1:17,200 and a raster whose cell size is 6 inches, what resolution (dots per inch) should I choose when exporting the map to an image in order to maintain the raster's native resolution (1 cell = 1 dot). Need a quick refresher on the conversion formulas.


It is possible to edit selected features with Python?

I'm trying to edit some fields of the selected features in ArcMap. It's possible to do this with Python (ArcPy)?


IUniqueValueRenderer with three fields

I have found some ArcObjects code to get unique values on one field for a IUniqueValueRenderer but I can't conceptually wrap my head around how to do that with 3 fields. Would this be .Net code or is there a property/method in the IUniqueValueRenderer? FYI, my goal is to get a list of layer files with python (these layer files have symbology from a master symbology layer which generally has more unique value combinations than what the target layer file has..meaning that some of the values in the IUniqueValueRenderer needs to be removed with the RemoveValue) then call an this arcobjects exe that will, once passed the path to a layer, finds it's unique values (from 3 fields) then check the renderer and start to remove the ones that don't apply. Thanks.


Why is my .prj file not being read by QGIS?

First off, I am new to GIS. I am using qgis (1.8.0), and have been given a shapefile containing .shp, .prj, .shx, .dbf, however qgis will not load the .prj file stating an invalid data source. Altering the project `properties>CRS to WGS84 Epsg:4326` gives correct latitude for the project but incorrect longitude. The .prj file reads: > > PROJCS["British_National_Grid",GEOGCS["GCS_OSGB_1936",DATUM["D_OSGB_1936",SPHEROID["Airy_1830",6377563.396,299.3249646]],PRIMEM["Greenwich",0.0],UNIT["Degree",0.0174532925199433]],PROJECTION["Transverse_Mercator"],PARAMETER["False_Easting",400000.0],PARAMETER["False_Northing",-100000.0],PARAMETER["Central_Meridian",-2.0],PARAMETER["Scale_Factor",0.999601272],PARAMETER["Latitude_Of_Origin",49.0],UNIT["Meter",1.0]]


Can I create OGC Services via Geoserver automatically?

I am new to GeoServer and wonder how I can create OGC Services via Geoserver automatically. I have to publish my data (Shapefiles) with a Java program as OGC Services. How can I do that? BTW: I checked out the source code of Geoserver, but I don't know how to package as right jar file to use.


ArcGIS - VBA - How can I create perpendicular polylines of a specified length that retain the attributes of the original polylines

i, I am using ArcGIS 9.2 and am desperately trying to create perpendicular polylines to existing polylines that I have. I would want the perpendicular lines to cross the original lines at their midpoint and I want the perpendicular lines to maintain the attributes of the original lines and finally I'd like to be able to control how long I want the perpendicular lines to be (eg. 50 metres) A script I found and which I've pasted below, takes me a good part of the way there in that it creates perpendicular lines at the midpoints but it doesn't retain the attributes or allow me to set the length. Can anyone suggest what I need to change to bring over the attributes and let me set the length of lines that are drawn. Thank you very much for any help that you can give me. I'm new at this so I apologize if I haven't passted the code in the propoer formatting                   Option Explicit     Sub Test()         Dim pEditor As IEditor         Set pEditor = Application.FindExtensionByName("ESRI Object Editor")         If pEditor.EditState <> esriStateEditing Then             MsgBox "start editing first"             Exit Sub         End If          Dim pEL As IEditLayers     Set pEL = pEditor     If pEL.CurrentLayer.FeatureClass.ShapeType <> esriGeometryPolyline Then         MsgBox "target layer must be polyline"         Exit Sub     End If          If pEditor.Map.SelectionCount = 0 Then         MsgBox "nothing selected"         Exit Sub     End If          Dim pEnumFeat As IEnumFeature     Set pEnumFeat = pEditor.Map.FeatureSelection     pEnumFeat.Reset     pEditor.StartOperation     Dim pFeat As IFeature     Set pFeat = pEnumFeat.Next     Do Until pFeat Is Nothing         If pFeat.Shape.GeometryType = esriGeometryPolyline Then             Dim pNewFeat As IFeature             Set pNewFeat = pEL.CurrentLayer.FeatureClass.CreateFeature             Set pNewFeat.Shape = MakePerp(pFeat.Shape, _                                           GetLength(pFeat.Shape))             pNewFeat.Store         End If         Set pFeat = pEnumFeat.Next     Loop     pEditor.StopOperation "makeperps"     Dim pAV As IActiveView     Set pAV = pEditor.Map     pAV.Refresh     End Sub          Function GetLength(pPolyline As IPolyline) As Double         GetLength = pPolyline.Length     End Function          Function MakePerp(pBaseLine As IPolyline, ByVal dLen As Double) As IPolyline          Dim pPointColl As IPointCollection     Set pPointColl = New Polyline     Dim pPerpLine As ILine     Set pPerpLine = New Line          pBaseLine.QueryNormal esriNoExtension, 0.5, True, dLen / 2#, pPerpLine     pPointColl.AddPoint pPerpLine.ToPoint     pBaseLine.QueryNormal esriNoExtension, 0.5, True, dLen / -2#, pPerpLine     pPointColl.AddPoint pPerpLine.ToPoint     Set MakePerp = pPointColl     End Function


Opensource Tools for Optimizing Server Performance

I use a server with a basic Microsoft OS primarily for data storage and as a central repository for edited vector data. These data are not intended to be used as mapping or GIS products available publicly on the web. Additionally, I have always used a standard file structure for organizing data on the server, where files exist within folders. Given this information, I am interested in optimizing my server performance with opensource database management tools such as PostgreSQL and PostGIS. I do not have a strong IT background and am looking for a starting point. So far, I have downloaded PostgreSQL/PostGIS, although I find the process of loading data cumbersome and find myself asking why go through the hassle?. Questions:   * Where should I begin in bringing my GIS server and data management practices up to par?   * What opensource/commonly available software can be utilized to optimize a server designed for editing vector data and storing GIS data?


How to copy geodatabase on ftp site to local geodatabase using Python?

There is a geodatabase sitting on an ftp site that I would like to download with a Python script. Right now I'm thinking one way to do this is to copy the ftp geodatabase to a geodatabase on my computer. Below is the script I've started. Does anyone know how I can alter this script so that I obtain ftp gdb? Thank you UPDATE: Below is my final, working code as based on the answer @om_hennners provided. Thanks!               import arcpy, os, sys     from arcpy import env     arcpy.env.overwriteOutput = True     from ftplib import FTP          directory = "/group/geodb" #location of gdb on ftp     folder = "D:\\temp\\"      out_gdb = "data.gdb"     out_path = folder + os.sep + out_gdb     copy_gdb = "hydro.gdb" # This is the gdb I would like to copy from the ftp  site     ftp = FTP("10.4.2.22")     ftp.login("user", "pass")          ftp.cwd(os.path.join(directory, copy_gdb))     print "Changed to " + os.path.join(directory, copy_gdb)          filenames = ftp.nlst()     print filenames          print "starting to write"     for f in filenames:         with open(os.path.join(out_path, f), 'wb') as local_file:         ftp.retrbinary('RETR '+ f, local_file.write)                     ftp.close()     print "closed ftp connection"


Differences between ExtJS and Dojo Framework

I have been using OpenLayers and the ESRI's JavaScript API and found that both using different set of JavaScript libraries have some implications on the quality of application that they deliver. Dojo has three dominions, including dijit and dojox but there is no such aspect in ExtJS. What I have come to believe is that, this is the most prominent feature that distinguishes Dojo from ExtJS. This provides the opportunity to change the look and feel of Dojo based apps more easily as compared to ExtJS. The look and feel of ExtJS can be very easily noticed. I am still at a pretty early stages of learning and so might have a wrong notion, but I would definitely appreciate if you could suggest as to what makes ExtJS framework different from Dojo Framework? Also, I keep stumbling on these terminologies of Framework and Libraries, could you also clarify, when are they Framework and when are they Library or do both mean the same?


shape file coordinate system

I added a shape file of Iran to Qgis. 'Tehran' has 35.6833° N, 51.4167° E coordiantes, however Qgis shows Tehran with 4772860,31296767. How i can view coordinates of a city or location in latitude and longitude form in qgis? Thanks.


Vendor who sells Zip+4 from Semaphore corporation

I need a vendor who sells zip4 CD from Semaphore corporation. I would like to input the address data and zip and then mapped/output that correct address and Zip+4.


GDAL RasterizeLayer Doesn't Burn All Polygons to Raster

I am trying to burn a shapefile to a raster using GDAL's RasterizeLayer. I pre-create an area of interest raster from a different shapefile, given a specific pixel size. This AOI then serves as a base for all following rasterizations (same number of collumns and rows, same projection and geotransform). The issue occurs, however, when I go to burn the shapes to their own raster, based on the same pixel size and projections. The link below (don't have enough rep to post the image), shows the original shapefile in tan, and the dark pink where RasterizeLayer has burned data. The light pink is nodata values for the dark pink raster data. The grey is the AOI based on which the shapefile burn was completed. Given the extents of the shapefile polygons, I would expect to see burn values in the bottom two corners, as well as the two pixels underneath the data that does show. Obviously, however, this is not the case. Image for Problem- Finished Raster Burns As follows is the code that I used to generate these. All of the shapes were created using QGIS, and were all created in the same projection. (It should be noted that the gridding in the picture shown was just to give an idea of the pixel size I was using.)               from osgeo import ogr     from osgeo import gdal          aoi_uri = 'AOI_Raster.tif'     aoi_raster = gdal.Open(aoi_uri)          def new_raster_from_base(base, outputURI, format, nodata, datatype):              cols = base.RasterXSize         rows = base.RasterYSize         projection = base.GetProjection()         geotransform = base.GetGeoTransform()         bands = base.RasterCount              driver = gdal.GetDriverByName(format)              new_raster = driver.Create(str(outputURI), cols, rows, bands, datatype)         new_raster.SetProjection(projection)         new_raster.SetGeoTransform(geotransform)              for i in range(bands):             new_raster.GetRasterBand(i + 1).SetNoDataValue(nodata)             new_raster.GetRasterBand(i + 1).Fill(nodata)              return new_raster          shape_uri = 'activity_3.shp'     shape_datasource = ogr.Open(shape_uri)     shape_layer = shape_datasource.GetLayer()          raster_out = 'new_raster.tif'          raster_dataset = new_raster_from_base(aoi_raster, raster_out, 'GTiff',                                     -1, gdal.GDT_Int32)     band = raster_dataset.GetRasterBand(1)     nodata = band.GetNoDataValue()          band.Fill(nodata)          gdal.RasterizeLayer(raster_dataset, [1], shape_layer, burn_values=[1])      My question therefore, is: Is this a bug in GDAL? Or is RasterizeLayer burning data based on something other than just the presence or lack of a polygon within a specified pixel area? Edit: The files that I was using can be found here.


How to access existing graphs in arcGIS 10

I am trying to access an existing graph in arcGIS version 10. The student workbook I am using was written for version 9.3. Version 10 does not have "tools" option in the toolbar. The instructions guide students to the 'Tools' option.   I will be teaching students to use arcGIS using 'Mapping our World' student workbooks (written for version 9.3) using version 10.


ArcSDE for PostgreSQL installation

I am trying to install ArcSDE 10 for PostgreSQL in CentOS. I tried doing that but am stuck with same errors repeatedly and am thinking of giving one more shot from the scratch. I just want to make sure the versions/releases before I start over again. I am thinking of using CentOS 5.4 (for i686 architecture), PostgreSQL 8.3.8/8.4.1, PostGIS 1.4.0 and ArcSDE 10. Please let me know if I have to change any version/releases of any software. Thanks!


four color theorem for QGIS 2.0 or Leaflet?

I recently looked at some other posts relating to this topic but as far as I can tell there is currently not a working plugin or script for QGIS 2.0 to achieve a four color theorem for styling polygon vector data (it doesn't necessarily have to be 4 colors, just so that each neighboring polygon is a different color, though keeping it at a minimum number of 4-6 colors would be elegant). I will ultimately be displaying my data using Leaflet JS, so if there is a plug-in for Leaflet to accomplish the same that would work as well. If anyone has a working script or plugin to share I would appreciate it, thanks!


LayerSwitcher problem in FireFox: layer after turn off and on are not displayed

I have a problem displaying layers using the LayerSwitcher. The problem only occurs on FireFox and not in the IE.   On some applications, layers which I turn off and again on, are not displayed until I do a pan or zoom. Is this a problem with Firefox or is something wrong with the website?


How to change icon size toolbar in Arcgis Viewer for Silverlight?

i have created a new web Map using Arcgis Viewer for Silverlight,and i want to change icon size toolbar from small to large size, and how to disable base map in my app, because i want to add imagery instead of base map. ![Free Web Proxy](http://www.freeimagehosting.net/newuploads/krf44.png)


Export .mxd to .dxf

Converting a .mxd to .dxf is something we need to do, and we do not only need the layers but also the layout of a map on our .mxd to be transfered to .dxf. I've done some research on the internet but it seems like an very difficult feat. In the **ArcGIS Discussion Forums** , some people mention the possiblity of going from .mxd to .ai (Adoble Illustrator) followed by a conversion from .ai to .dxf, have any one tried it ? Others indicate that they are working with .mxd to .pdf then .pdf to .dxf, I gave this a try but the layout is not well converted. Is there any other way we might proceed to have better results ?


How to add Opencarto to a Geotools project?

I've installed geotools through maven (and m2ecplise) but i'm not able to add opencarto library to my project. (As you can guess, I'm new to java) I tried by adding opencarto's main folder as library through the Project panel, but still, maven is not proposing me opencarto as dependency when I try to find it. Any idea of what i'm doing wrong? (this is probably a question for you Julien. I sent you a mail through sourceforge. I'm not sure you got it ... anyway, it's a good thing to make noise around opencarto since it seems very interesting!) Thank you


Are there any open source tools to detect and fill sinks on a DEM?

Are there any open source or inexpensive tools to detect and fill sinks on a DEM? ESRI Spatial Analyst is just out of my price range.


Creating a .shp with OSM and import it to POSTGIS

![enter image description here](http://i.stack.imgur.com/PZT9u.png) Any suggestions ? I can't launch osm2po I don't know why, when I click on the jar file, nothing happen :p


Conversion of Latitude/Longitude to different map projection

I have two m*n matrices coming from a WRF-ARW simulation: one exclusively with longitude values and one with only latitude values. For both matrices I know the originatin projection (lambert conic conformal, with true latitudes and parameter set by me). I would like to convert them to a regular plate carree projection (the one employed by Google Earth, to be precise): is it possible to do so? If yes, which would be the best instrument and how, exactly? I'm currently having at hand R (where I've and may have both QGis and MATLAB (but without any toolbox) Thanks for any help!


Load and display a layer using PYQGIS

I'm having a go at using pyqgis from within the QGIS python console to load a layer from SQL Server and display it. I feel like I'm part way there but I'm having a real issue getting the layer to display. I've loaded the layer, it appears in the layers dock but nothing appears in the map window (even after selecting zoom to layer extent and checking styles). What's also weird is that if I open attributes there is data in there and I can save the layer as a shapefile, reopen it and it appears fine. Is there something I'm missing? Here's my code:               from PyQt4.QtCore import QFileInfo,QSettings     from qgis.core import QgsRasterLayer, QgsCoordinateReferenceSystem     from qgis.gui import *     from PyQt4.QtCore import *          s = QSettings()     oldValidation = str(s.value( "/Projections/defaultBehaviour", "useGlobal" ))     s.setValue( "/Projections/defaultBehaviour", "useGlobal" )          uri = "MSSQL:server=My_Server;database=My_Database;tables=dbo.My_View;trusted_connection=yes"     vlayer = QgsVectorLayer(uri, "my_table", "ogr")          vlayer.setCrs( QgsCoordinateReferenceSystem(4326, QgsCoordinateReferenceSystem.EpsgCrsId) )          s.setValue( "/Projections/defaultBehaviour", oldValidation )          QgsMapLayerRegistry.instance().addMapLayer(vlayer)


FME Data migration for Feature class with Attachments

I want to migrate data from one point feature class to another point feature class in same Database schema using FME job. I can easily do this task without using any transformers. But now both point feature classes (reader and writer) has attachment tables. Attachment table has created by using arc catalog. The feature class and its attachment table have some relations. Can anyone please let me know how i can migrate data for attachment tables while I do data migration between feature classes?


How to do the loop when using GRASS module r.water.outlet?

I have seen this post:Watershed Analysis with GRASS, and I have a .txt file which contains coordinates of many points, and I want to use r.water.outlet to get basin area of each point. I type the following code in the post:               while read X Y; do     .... ;     done < cross_pts.txt      and in GRASS command line, it shows: ![enter image description here](http://i.stack.imgur.com/XDRvK.png) I already import the DEM file and successfully run r.watershed to get the drainage file in GRASS, but I do not know how to follow the guide in the post. I am new to using command line and know nothing about bash script. How can I loop through all different coordinates? I am using GRASS 6.4.3 under windows8.


Raster: mosaic first, or project first?

I have a thousand raster tiles in un-projected decimal degrees from which to build a seamless mosaic in a metres-based coordinate system. The final image will be about 70,000px by 50,000px. Which is the better approach, thinking in terms of processing (time) efficiency and quality of result:   * mosaic into a single image and then transform to the final projection?   * or transform each tile and then mosaic?   * (or don't mosaic at all and use a virtual raster table in the final projection?) I can use either or any of ArcGIS Desktop, GDAL Utilities, QGIS and a have 64bit OS with 12GB RAM.


Using ModelBuilder, how to define directory while using a Python script?

I have a model (shown below) in which at the very end of it I run a Python script ('SumFieldInsertNew'). One issue I'm having is that I'm creating this model for an end-user that will have little to no knowledge of most things computer-wise. Ergo, I'm trying to make it as hassle-free and flexible as possible. My question is, I want the user to be able to name the file the script will run itself on, but am unsure of how to do something like this. So essentially, I want to make 'Empty FC' a parameter in which the user can decide what they want their file to be named. Then have the Python script directory code (where it's defining the feature class), know that whatever the user named it is conditional to it being able to run its processes. I appreciate any help and suggestions, thank you! ![enter image description here](http://i.stack.imgur.com/oNhpE.jpg) My model runs as follows:   1. It establishes a connection to an Oracle dB   2. It creates a Query Layer based on the pre-filled 'Query' parameter   3. It creates an empty feature class (Empty FC) into a pre-filled GdB parameter   4. It then exectures the 'Copy Features' tool on the Query Layer and outputs the copy into the Empty FC   5. This newly populated feature class is then given a new field, called 'SUM_Quantity_Solid'   6. In which lastly, this new field is populated via the Python script 'SUM_FIELD_INSERT_NEW', which (as seen below), SUMs up a field and then inserts the values. > >     > import arcpy >     > >     > # Define the feature class >         fc = r'User_Defined_Path' >     > >     > # find the unique 'SEGMENT_LENGTH' values >         Slist = list() for row in arcpy.da.SearchCursor(fc, > 'SEGMENT_LENGTH'): >     >     # if the value isn't in the list then add it to the list >     >     if not row[0] in Slist: >     >         Slist.append(row[0]) >     > >     > for Value in Slist: >     >     # definition query to limit the rows in the cursor >     >     DefQ = 'SEGMENT_LENGTH = ' + str(Value) >     > >     >     # Use a generator expression to populate a list from the > 'QUANTITY_SOLID' field >     >     b = sum(row[0] for row in arcpy.da.SearchCursor(fc, > 'QUANTITY_SOLID'),DefQ) >     > >     >     with arcpy.da.UpdateCursor(fc, ['QUANTITY_SOLID_SUM'],DefQ) as > cursor: >     >         for row in cursor: >     >             row[0] = b >     >             cursor.updateRow(row) >   As referenced, I would like to be able to have > fc = r'User_Defined_Path' be what the user initially defined the 'EMPTY FC' filename parameter as.


Analyse effect of groups on raster maps

I work with a mechanistic model that uses climate data to simulate some specific output for different time steps and spatial locations. Now I want to investigate the effect of different climate data sources on the model output. To do this, I varied the climate data and ran the model with otherwise constant settings; the groups are:   * Climate model that produces the climate data   * Emission scenarios with different assumption underlying the climate data from the same climate model   * Spatial resolution of the climate data   * Downscaling method to arrive at the specific spatial resolution I have produced now a bunch of rasters with my model driven by the varying climate data. Now I would like to extract the effects of the different groups of climate data on the model output and produce raster maps showing at which spatial locations there is the most variation in model output. I was thinking about applying a principal component analysis to avoid pairwise comparisons of the simulations. Would a PCA be an appropiate method or is there a better way of comparing the different groups? And what would be the next step after doing the PCA and calculating the loadings?


R: Download a large DEM, change projection, and adjust to smaller scale

This is a process that takes just a few seconds in GIS software. My attempt to do it in R uses a large amount of memory then fails. Is there something wrong in my code, or is this just something R cannot do? I have read R can work inside Grass, can I use a Grass function from inside R?               library(raster)          # I have many environmental rasters in this format     new_r <- raster(ncol=615, nrow=626, xmn=-156.2, xmx=-154.8, ymn=18.89, ymx=20.30)     res(new_r) <- 0.00225     projection(new_r) <- "+proj=longlat +ellps=GRS80 +datum=NAD83 +no_defs +towgs84=0,0,0"          R> new_r ### not too big with a few hundred cells per side     class       : RasterLayer      dimensions  : 627, 622, 1  (nrow, ncol, nlayers)     ncell       : 389994      resolution  : 0.00225, 0.00225  (x, y)     projection  : +proj=longlat +ellps=GRS80 +datum=NAD83 +no_defs +towgs84=0,0,0      extent      : -156.2, -154.8, 18.89, 20.3  (xmin, xmax, ymin, ymax)     values      : none          # I get the DEM at much higher resolution (zipfile is 182Mb)     zipurl <- "ftp://soest.hawaii.edu/coastal/webftp/Hawaii/dem/Hawaii_DEM.zip"     DEMzip <- download.file(zipurl, destfile = "DEMzip")     unzip("DEMzip", exdir = "HIDEM")     HIDEM <- raster("HIDEM/hawaii_dem")          R> HIDEM ### 10m resolution, file is way too big     class       : RasterLayer      dimensions  : 15067, 13136, 1  (nrow, ncol, nlayers)     ncell       : 197920112      resolution  : 10, 10  (x, y)     projection  : +proj=utm +zone=5 +ellps=GRS80 +datum=NAD83 +units=m +no_defs +towgs84=0,0,0      extent      : 179066, 310426, 2093087, 2243757  (xmin, xmax, ymin, ymax)     values      : HIDEM/hawaii_dem      min value   : 0      max value   : 4200           # the following line fails (after a long time)     new_HIDEM <- projectRaster(HIDEM, new_r)


Is it possible to convert ArcGIS 10 styles to version 9.3?

I have a style file (*.style) generated in ArcGIS 10.0. And I need to use symbols from this file on machine with ArcGIS 9.3. Is there some way to convert style file from format 10.0 to 9.3?


Can GeoServer 2.0.2 access ArcSDE 10.0 databases?

It looks like ArcSDE 10.0 databases are not yet supported... At least documentation on the geoserver-2.0.2-arcsde-plugin doesn't appear to address it specifically, and I can't get it to work. Is that right? Thanks!


What product does "EDN ArcGIS Enterprise Advanced" key unlock?

I have a key to unlock something called "EDN ArcGIS Enterprise Advanced" I tried to unlock ArcEditor with the key - seemed like the obvious choice, but that didn't work. The key starts with first three letters "ECP". It seems like it unlocks a Server, but I thought the server license cost $30K and I don't think we paid for that license, but it is available to me in the EDN portal customers.esri.com


Openlayers 2.1x Load WMS Layers From PHP

I have been looking to see if it's possible to add WMS layers afrom a PHP script. I'd like to do this to move a lot of the code from the client side to the server. Currently I load the layers in the initialize function with:               /*Declare layers + properties:*/     var layer1= new OpenLayers.Layer.WMS("Layer1", "/geoserver/wms", {layers: 'Workspace:Layer1', styles: '', srs: 'EPSG:27700', format: format, tiled: 'true',transparent: 'TRUE'}, {tileSize: new OpenLayers.Size(540,640), buffer: 0, displayOutsideMaxExtent: true,isBaseLayer: false,visibility : false,transitionEffect: 'none' } );     var layer2= new OpenLayers.Layer.WMS("Layer2", "/geoserver/wms", {layers: 'Workspace:Layer2', styles: '', srs: 'EPSG:27700', format: format, tiled: 'true',transparent: 'TRUE'}, {tileSize: new OpenLayers.Size(540,640), buffer: 0, displayOutsideMaxExtent: true,isBaseLayer: false,visibility : false,transitionEffect: 'none' } );          /*Add layers to layer array*/     layerArray.push({key : 'layer1', lyr : layer1});     layerArray.push({key : 'layer2', lyr : layer2});          /*Add layer to the map*/     for(var n =0; n<layerArray.length; n++)     {         obj = layerArray[n];         map.addLayer(obj.lyr);     }      Is there a way to move the layer information and layer array population to PHP so the only code in the initialize function is the code to add the layers to the map? Thanks for any advice, ** Edit ** Thank you for the advice. More details as requested: I have been asked to produce a site that contains an Openlayers 2.1x map that utilises vector drawing, select, measure, geolocate, cql filters and all mapping layers are stored in Geoserver 2.33 and are accessed using WMS. Geoserver is run on the webserver. I have completed all of this using javascript to populate the layer names (as indicated above) and sketchSymbolizer properties etc. I also use javascript to create any calculations and cql filters. The client who has asked for the site has expressed concerns with the ease of taking the javascript and replicating the site, plus he feels large calculation in javascript are affecting the performance of older smart phones. I have moved calculations to php and got them running smoothly (if not a little slower for fast machines, but adding a loader and everyone is happy again), but there is still a want to move all of the openlayers initialize items into php as well. I can not get openlayers to initialize, I have used the code below to create an array, but feel I'm missing the point completely.               <!--Demo code held in html for proof of concept before moving to an external file-->     <?php          $al = array(     array(     "key" => "layer1",     "lyr" => "new OpenLayers.Layer.WMS("Layer1", "/geoserver/wms", {layers: 'Workspace:Layer1', styles: '', srs: 'EPSG:27700', format: format, tiled: 'true',transparent: 'TRUE'}, {tileSize: new OpenLayers.Size(540,640), buffer: 0, displayOutsideMaxExtent: true,isBaseLayer: false,visibility : false,transitionEffect: 'none' } )"     ),      array(      "key" => "layer2",      "lyr" => "new OpenLayers.Layer.WMS("Layer2", "/geoserver/wms", {layers: 'Workspace:Layer2', styles: '', srs: 'EPSG:27700', format: format, tiled: 'true',transparent: 'TRUE'}, {tileSize: new OpenLayers.Size(540,640), buffer: 0, displayOutsideMaxExtent: true,isBaseLayer: false,visibility : false,transitionEffect: 'none' } )"     )     );     $arrayObject = new ArrayObject($al);     ?>     <script type="text/javascript">         var layerArray= <?php echo json_encode($arrayObject); ?>;     </script>      This produces the error "annot read property 'wrapDateLine' of null". In light of "why" do I want to, my hope was to hide variables used for creating the layers and other elements, but if there is a more effective / better way, I'm very open to trying anything. Kind regards,


unable to create or edit any layer on geoserver

ve been using geoserver for some time now with no problems untill now. After adding some layers and working on my layouts now I am unable to add or modify layers. I can add stores and create workspaces but if I try to add a new layer, after clicking the save button the server just reloads in the same "Edit Layer" page and no changes are saved nor a layer is added. Tryed to reboot, uninstall, install windows service, WAR, different version and got the same result. Then I tried to restor my geoserver data dir from my backup and layers are back but I can not modify nor add new ones!! Geoserver logs show nothing unusual, no errors, no warnings no nothing!! May it be some windows configuration that went bad? I´m on windows server 2008 Please Help me!!


How do you use GIS to detect percentage of Land Use Change

I am a student in GIS I have a dissertation about land use change. I have questions about using ArcMap 10 for this aim, how can using GIS to detect the percentage of land use change. For example we want to identify land use change in Sheffield city between 1990 until 2000 Land use of Urban %30 1990 also it is changed to %40 in 2000 how can detect this result I mean (%30 and %40).


Compressions to support when Implementing GeoTiff Reader?

I'm writing a GeoTiff reader for GeoTrellis for Google Summer of Code. I have read in all tags and geokeys and now I'm working with the decoding process. Currently the Tiff 6.0 spec recommends that baseline readers support the following decompressions:   * Uncompressed (duh)   * PackBits   * Huffman And the extension recommends the support of:   * LZW   * Jpeg (Old spec) I have implemented all the baseline compression types and LZW but as I understand it the compression type 6 (Jpeg old spec) is deprecated and replaced by compression type 7 (new Jpeg). Except these compression types mentioned Deflate compression (compression type 8) is mentioned in an extension document by Adobe. As I understand it nobody uses the old Jpeg compression, is there any **common** other compressions I should support? That isn't bound for a single reader/writer software brand? And I should implement the compression type 7 and 8? The goal with the reader is to support as many GeoTiffs as possible. For the interested the reader is written in the Scala programming language.


Displaying coordinate references in GeoTools?

I'm using this library to finish my final project and I have to say I'm no familiar with the maps worlds. My question is if geotools has a method to show automatically the coordinates references in the map, and how use it? That would be really helpful because I spend a lot of time trying to find it without success. I use the quickstart example to show a shapefile with the worlds map. The user can select a rectangular zone of interes and then the application draw the Sea Surface Temperature. The map generated needs the references of where are positioned because we export it to a png image. One example of what I am looking for are the blue lines in the map below. ![enter image description here](http://i.stack.imgur.com/YxBLB.jpg)


Loading OSM data into QGIS

I have been trying to use OSM maps and my data is 60 GB around. I have imported it to QGIS as a vector layer but it seems the data is so big that the rendering is taking forever. My system has Intel Dual Core and 2GB of RAM. So can anyone please tell me how i can work with my data without any rendering problems. please help me work around this problem.


Overpass XML query with multiple polygons or bounding boxes

I am having trouble creating an OpenStreeMap Overpass XML combined polygon/bounding box query to return what I want. Can someone provide a simple example or correct the following XML? I want to construct a query to retrieve data from multiple polygons or bounding boxes. Specifically, I want to retrieve all data within a small polygon and only certain data types, major roads, in a larger containing bounding box and return the combined result. This XML shows conceptually what I want to do but is returning all data for the larger area rather than all data for the small area and major roads for the larger area:               <osm-script timeout="20" element-limit="1073741824">         <union into="_">             <!-- Smaller area, all data -->             <polygon-query into="_" bounds="54.12 -2.779  54.143 -2.779  54.143 -2.761 54.12 -2.761" />             <recurse from="_" into="_" type="up"/>             <recurse from="_" into="_" type="down"/>         </union>               <union into="_">             <!-- Larger area, major roads only -->             <bbox-query s="54.12" n="54.143" w="-2.779" e="-2.744" />             <has-kv k="highway" v="motorway"/>                <recurse from="_" into="_" type="up"/>             <recurse from="_" into="_" type="down"/>         </union>         <print from="_" limit="" order="id"/>     </osm-script>


Gdal creating raster rotated wrong way

I have been trying to read in a stack of images for a timeseries and created an output raster that contains the max pixel values of the timeseries. The arrays were created using numpy               img_stack = np.empty((image_width, image_length, len(image_files)), np.dtype('f'))           final_values = np.empty((image_length, image_width), np.dtype('f'))      I imported the image using the following loop:               for i, fname in enumerate(image_files):         img = gdal.Open(os.path.join(fn + ("\{0}".format(fname)))).ReadAsArray()         img_stack[:, :, i] = img      then I read the values into a 2D list using the following nested loop:               for x in range(image_length):         for y in range(image_width):             final_values[x][y] = 0             for z in range(len(image_files)):                 if (final_values[x][y] < img_stack[y, 0, z]):                     final_values[x][y] = (img_stack[y, 0, z])                 else:                     next      then use               output_raster = dr.Create(final, image_width, image_length, 1, gdal.GDT_Float32)     output_raster.GetRasterBand(1).WriteArray(final_values)      the output contains the correct pixel values, but is rotated 90 degrees to the left. I believe the issue is due to GDAL reading the image as an array first by width and then by length, but I am not sure. Is there an issues with my code or is there a function I could use to rotate it to the proper orientation? Thanks


Getting fail to store Polygon Shape Object in Geometry Datatype with hibernate spatial API

Tried to store Polygon shape in database with HibernateSpatial API. I have followed this tutorial to store values. I am using Mysql as my database server. So Made some couple of changes for my scenario. Getting the following ERROR to store it in the database. **: _:ERROR:_ :** **ERROR [org.hibernate.engine.jdbc.spi.SqlExceptionHelper] (http-/192.168.5.244:8080-1) Data truncation: Cannot get geometry object from data you send to the GEOMETRY field** **ERROR [com.service.HibernateSaveService] (http-/192.168.5.244:8080-1) Exception :_: could not execute statement** My code is below. **ENTITY CLASS :_:**               package com.entity;          // Generated Mar 1, 2014 12:22:56 PM by Hibernate Tools 3.4.0.CR1          import static javax.persistence.GenerationType.IDENTITY;          import javax.persistence.Column;     import javax.persistence.Entity;     import javax.persistence.GeneratedValue;     import javax.persistence.Id;     import javax.persistence.Table;     import javax.persistence.UniqueConstraint;          import org.hibernate.annotations.Type;          import com.vividsolutions.jts.geom.Polygon;          /**      * HibernateSpatial generated by hbm2java      */     @Entity     @Table(name = "HIBERNATE_SPATIAL", catalog = "geocodedb", uniqueConstraints = @UniqueConstraint(columnNames = "GEO_REGION_NAME"))     public class HibernateSpatialEntity implements java.io.Serializable {              /**          *           */         private static final long serialVersionUID = 1L;         private Integer prdGeoRegionId;         private String geoRegionName;         private String geoRegionDescription;              @Type(type="org.hibernate.spatial.GeometryType")         private Polygon geoRegionDefinition;              @Id         @GeneratedValue(strategy = IDENTITY)         @Column(name = "PRD_GEO_REGION_ID", unique = true, nullable = false)         public Integer getPrdGeoRegionId() {             return this.prdGeoRegionId;         }              public void setPrdGeoRegionId(Integer prdGeoRegionId) {             this.prdGeoRegionId = prdGeoRegionId;         }              @Column(name = "GEO_REGION_NAME", unique = true, nullable = false, length = 50)         public String getGeoRegionName() {             return this.geoRegionName;         }              public void setGeoRegionName(String geoRegionName) {             this.geoRegionName = geoRegionName;         }              @Column(name = "GEO_REGION_DESCRIPTION", length = 70)         public String getGeoRegionDescription() {             return this.geoRegionDescription;         }              public void setGeoRegionDescription(String geoRegionDescription) {             this.geoRegionDescription = geoRegionDescription;         }              @Column(name = "GEO_REGION_DEFINITION")         public Polygon getGeoRegionDefinition() {             return this.geoRegionDefinition;         }              public void setGeoRegionDefinition(Polygon geoRegionDefinition) {             this.geoRegionDefinition = geoRegionDefinition;         }          }      **PERSIST CLASS :_:**               package com.service;          import java.util.Random;     import javax.ejb.Stateless;     import javax.ejb.TransactionManagement;     import javax.ejb.TransactionManagementType;     import javax.persistence.EntityManager;     import javax.persistence.Persistence;     import javax.ws.rs.Consumes;     import javax.ws.rs.POST;     import javax.ws.rs.Path;     import javax.ws.rs.Produces;          import org.hibernate.Session;     import org.slf4j.Logger;     import org.slf4j.LoggerFactory;          import com.solartis.entity.HibernateSpatialEntity;     import com.vividsolutions.jts.geom.Geometry;     import com.vividsolutions.jts.geom.Polygon;     import com.vividsolutions.jts.io.ParseException;     import com.vividsolutions.jts.io.WKTReader;          @Stateless     @TransactionManagement(TransactionManagementType.BEAN)     @Path("/")     public class HibernateSaveService {         /**          * mLogger - attribute used for to write the log          */         private static final Logger mLogger = LoggerFactory                 .getLogger(HibernateSaveService.class);              @POST         @Path("/saveHibernate")         @Consumes("*/*")         @Produces("*/*")         public String saveShapes(String jsonRequest) {             String statusResult =   null;             String methodName = "saveShapes";             StringBuilder regionName    =   null;             Session session =   null;             Random random   =   null;             HibernateSpatialEntity hibernateEntity  =   null;             Polygon poly    =   null;             String singleShape  =   "Polygon((-50.123456789 -50.123456789,50.98765432175 -50.753951654852,50 50,-50 50,-50.123456789 -50.123456789))";                  mLogger.debug("Entered inside the Method :_: " + methodName );             Geometry geom   =   wktToGeometry(singleShape);             try {                 EntityManager entity = Persistence.createEntityManagerFactory(                         "spatialPersistenceUnit").createEntityManager();                 session = (Session) entity.getDelegate();                 if (!geom.getGeometryType().equals("Polygon")) {                     throw new RuntimeException("Geometry must be a Polygon. Got a " + geom.getGeometryType());                 }                 else {                     //poly.                     mLogger.info("GeoType :_: " + geom.getGeometryType());                     hibernateEntity =   new HibernateSpatialEntity();                     random  =   new Random();                     hibernateEntity.setGeoRegionDefinition((Polygon)geom);                     hibernateEntity.setGeoRegionDescription("ABCDEFGHIJKLMNOPQRSTUVWXYZ");                     regionName  =   new StringBuilder();                     regionName.append("RegionName :_: ").append(random.nextInt(10000));                     hibernateEntity.setGeoRegionName(regionName.toString());                     mLogger.info("StoredName :_: " + regionName.toString());                     session.beginTransaction();                     session.save(hibernateEntity);                 }                  } catch(Exception exception) {                 mLogger.error("Exception :_: " + exception.getMessage());             } finally {                 session.clear();                 if ((session != null) && (session.isOpen())) {                     session.close();                 }             }             mLogger.debug("Exited from the Method :_: " + methodName);             return statusResult;         }              private Geometry wktToGeometry(String wktPoint) {             WKTReader fromText = new WKTReader();             Geometry geom = null;             try {                 geom = fromText.read(wktPoint);             } catch (ParseException e) {                 throw new RuntimeException("Not a WKT string:" + wktPoint);             }             return geom;         }     }      **DATABASE SCHEMA**               +------------------------+--------------+------+-----+---------+----------------+     | Field                  | Type         | Null | Key | Default | Extra          |     +------------------------+--------------+------+-----+---------+----------------+     | PRD_GEO_REGION_ID      | int(11)      | NO   | PRI | NULL    | auto_increment |     | OWNER_ID               | decimal(3,0) | NO   |     | NULL    |                |     | GEO_REGION_NAME        | varchar(50)  | NO   | UNI | NULL    |                |     | GEO_REGION_DESCRIPTION | varchar(70)  | YES  |     | NULL    |                |     | GEO_REGION_DEFINITION  | geometry     | YES  |     | NULL    |                |     +------------------------+--------------+------+-----+---------+----------------+      Regards, ArunRaj.


Error: can't display dash lines (broken lines) when zooming PageLayout

I am working with PageLayoutControl, its mapframe include elements: lines, texts. And I am getting an error when zooming: If I zoom the pagelayout to 100% ( smaller than 110%), it can display "Dash Lines". But, if zoom to 120% (or larger), all "dash lines" become "solid lines". Please see attached screenshots:![zoom to 100%](http://i.stack.imgur.com/QHbSs.png) ![zoom to 120%](http://i.stack.imgur.com/Jfjsh.png) And another problem: when I print the pagelayout, all dash lines become solid lines. Are there anyone know what's wrong ? Thanks and regards, Tai


GeoExplorer: Setting more/custom zoom levels

I have BoundlessGeo's OpenGeo Suite (Geoserver,PostGIS,GeoWebCache,GeoExplorer) installed, and love the web functionality of GeoExplorer, but need my maps to zoom in a few steps further. So far the only way I can think would be to edit the GeoExplorer.js file (which appears a min/concat of Ext and OpenLayers) to change the default zoom levels, but that seems like a bad way to go. Any ideas? Is this functionality that should be added to GeoExplorer? Perhaps I could add zoom levels after the map is loaded? Cheers


How to perform Random Forest land classification?

This is a follow-up to a previous post: Machine Learning Algorithms for Land Classification. It seems that the Random Forest (RF) classification method is gaining much momentum in the remote sensing world. I am particularly interested in RF due to many of its strengths:   * A nonparametric approach suited to remote sensing data   * High reported classification accuracy   * Variable importance is reported Given these strengths, I would like to perform Random Forest land classification using high resolution 4 band imagery. There is a lot of material and research touting the advantages of Random Forest, yet very little information exists on how to actually perform the classification analysis. I am familiar with RF regression using R and would prefer to use this environment to run the RF classification algorithm. **How do I collect, process and input training data (i.e. based on high resolution CIR aerial imagery) into theRandom Forest algorithm using R? Any step-wise advice on how to produce a classified land cover raster would be greatly appreciated.**


Can I create a Styled Map for only one section of Google Maps?

I have been reading everything I can to find out how this can be done. It looks like the styled map applies only to the whole map, and I can't find a way to only apply it to a certain polygon region.


Random Dips/Pits in Elevation surface in ArcGlobe 10.2

I have an elevation raster which I load into ArcGlobe as an elevation source and appears to have large pits and dips in it when I place floating features such as roads and imagery. The odd thing is that the phenomenon does not show up in ArcScene and if I import the raster into a FGDB the pits disappear. ![enter image description here](http://i.stack.imgur.com/RgkpH.png) It's seriously looking like Apple Maps and any help would be nice.


How can I display GeoEye1 imagery in colour in Opticks?

Anyone familiar with the use of the open-source remote sensing software Opticks? Im new to opticks, I load GeoEye1 imagery but it displayed as greyscale. It supposed to be displayed as bands combination true colour or false colour composites.


Adding layers using ArcGISDynamicMapServiceLayer

I am right now adding individual layers using `FeatureLayer` but it takes quite some time to load as the map is for US Counties (3076 count). I tried to load it using `ArcGISDynamicMapServiceLayer` but I was not successful. It seems that `ArcGISDynamicMapServiceLayer` can load mapserver but not the layers in the mapserver. I used the following code but it did not work.               var ftrLyr = new esri.layers.ArcGISDynamicMapServiceLayer(url, {         visibleLayers: [1]     });      Is it not possible to load layers in a mapserver using `DynamicMapService`?


How can I create a portable script for coordinate conversions to run from USB?

I would like to be able to convert Lat long values in wgs84 to the swedish RT90, the values are in CSV files under headings X and Y. I would like to have this portable, an eg runned in a script in python idle or something similar. I have a python script today on USB that uses the pyshp module to convert csv files to shapefiles, but some of the csv are in wgs84 which would be nice to have converted in the same process. Can gdal be used, and can it be made portable?


How can I measure area from geographic coordinates?

If I have polygons in geographic coordinates (WGS84), how do I measure the total area each takes on the surface of the earth, taking into account the curvature of the earth?


How to convert PostGIS table to Shapefile in Python?

I want to convert a PostGIS table into a shapefile (without using pgsql2shp). In order to create a geometry in the shapefile I have to give the Xmin,Ymin and Xmax,Ymax, and the geometry which I have in my PostGIS table is an irregularly shaped one (I can get the exterior using bounding box but that will include some extra area more than my area of ineterest). Is there any method by which I can get the task done? I want to do the thing programmatically and using Python.


Seeking shapefiles of business locations

Does anyone know where to find shapefiles of businesses such as grocery stores or gas stations? Thanks!!


Moving "last name" from last to first in the Attribute Table

How do you get _e.g._ "John D. Doe" to "Doe, John D." in the AT using VB Scripting without going individually going through each line? I am using ArcGIS 9.3. I should note that not all the names are the same; they are different. Is it possible to select those in the column with different names and just reverse the "first, last name" to "last name, first"?


Tilecache positions my gdal_translated map 20km south of where I want it

I have created a geotiff file through gdal_translate:               gdal_translate -a_ullr 547812.2384341886 6834734.336439242 547967.5011478929 6834622.35081697 -a_srs EPSG:900913 pdz3.jpg pdz.tif      When I load this tif file in qgis, I see my map on the right location. Also, when I ask the gdalinfo from this file I get the following information:               Driver: GTiff/GeoTIFF     Files: pdz.tif     Size is 1024, 724     Coordinate System is:     PROJCS["Google Maps Global Mercator",         GEOGCS["WGS 84",             DATUM["WGS_1984",                 SPHEROID["WGS 84",6378137,298.257223563,                     AUTHORITY["EPSG","7030"]],                 AUTHORITY["EPSG","6326"]],             PRIMEM["Greenwich",0],             UNIT["degree",0.0174532925199433],             AUTHORITY["EPSG","4326"]],         PROJECTION["Mercator_1SP"],         PARAMETER["central_meridian",0],         PARAMETER["scale_factor",1],         PARAMETER["false_easting",0],         PARAMETER["false_northing",0],         UNIT["metre",1,             AUTHORITY["EPSG","9001"]]]     Origin = (547812.238434188649990,6834734.336439241655171)     Pixel Size = (0.151623743851815,-0.154676273856209)     Metadata:       AREA_OR_POINT=Area     Image Structure Metadata:       INTERLEAVE=PIXEL     Corner Coordinates:     Upper Left  (  547812.238, 6834734.336) (  4d55'15.89"E, 52d22'37.58"N)     Lower Left  (  547812.238, 6834622.351) (  4d55'15.89"E, 52d22'35.36"N)     Upper Right (  547967.501, 6834734.336) (  4d55'20.91"E, 52d22'37.58"N)     Lower Right (  547967.501, 6834622.351) (  4d55'20.91"E, 52d22'35.36"N)     Center      (  547889.870, 6834678.344) (  4d55'18.40"E, 52d22'36.47"N)     Band 1 Block=1024x2 Type=Byte, ColorInterp=Red     Band 2 Block=1024x2 Type=Byte, ColorInterp=Green     Band 3 Block=1024x2 Type=Byte, ColorInterp=Blue      I tested the 4d55'15.89"E, 52d22'37.58"N values in google map, and they also point to the proper location. However, when I want this geotiff tiled - through tilecache - the map is placed about 20km south of the point I want it to appear (and appears in qgis). When I georeference the map in qgis self, and send it to tilecache I get the map on the right location. Calling gdalinfo on the qgis referenced tif I get roughly the same information (a couple of meters difference, due to mismatch clicks ;-)) The tilecache configuration file of the map is               [cache]     type=GoogleDisk     base=/home/gworx/data/tilecache          [jvv-chan-5-pdz]     type=GDAL     file=/home/gworx/data/layer-files/1320850804986-pdz.tif     bbox=547812.238434,6834734.33644,547967.501148,6834622.35082     spherical_mercator=true     srs=EPSG:900913     metaTile=true     extension=jpeg     mime_type=image/jpeg     tms_type=google      The bbox values are in line with the values I've set when calling gdal_translate. I am pretty much clueless at this moment on how to get this to work. If anyone can point me in the right direction on what I could be doing wrong, I would greatly appreciate it!


Creating .shp from .dbf in ArcMap?

How do I create a shapefile using the latitude and longitude information in my dBASE table? I am using ArcMap 10.1


ArcObjects: How to get the midpoint of a set of polylines?

How can I find the midpoint (along the curve) of a set of polylines using ArcObjects?


how do you do a buffer analysis of a geographical location?

I want to use a buffer analysis to identify the geodemographics of Incinerators. I have downloaded the Output Area Classification (OAC) map and also know the locations of all the Incinerators I want to analyse from an online map. I want to know how i can perform a buffer analysis from the geographical point in which the Incinerators are located based on the OAC map.


How to add CloudMade map to ArcGIS Android?

I want to add a CloudMade map to Esri's Android HelloWorld application. I'm following the documentation, but when I run the app, LogCat throws an EsriService Exception and the app closes. Below is the relevant code. Does anyone have an idea on what I need to change to display a CloudMade basemap? I posted this question on ArcGIS Forum a few days ago, but no response there yet. I'm working with ArcGIS Runtime SDK for Android 10.2.2. Thanks, Patty               public class HelloWorld extends Activity {     MapView mMapView;           // Called when the activity is first created.      @Override      public void onCreate(Bundle savedInstanceState) {      super.onCreate(savedInstanceState);      setContentView(R.layout.main);      mMapView = (MapView)findViewById(R.id.map);      String[] tileServerURL = {"http://b.tile.cloudmade.com/{KEY}/1/256"};      OpenStreetMapLayer tileLayer = new OpenStreetMapLayer(tileServerURL, 1, 16, "CloudMade", true);      tileLayer.setKey("myKey");            mMapView.addLayer(tileLayer); //when this line is commented out the app runs, but only with an Esri base map


Creating rectangle over a buffer circle with specified angles/rotation?

I have a layer that is a circle buffer around a point. What I want to do is create a rectangle around the circle, but have every rectangle be perfectly lined up N/S (0degrees of angle/rotation). I have tried using Minimum Bounding Geometry, but the rectangles are all rotated a little different each time instead of all being the same. Is there any way to do this? Thank you for the help.


How I get SRID from Geometry field

what is the function for getting SRID from Geometry field. I have inserted a shapefile into SQLServer table with Geometry field defined. I know the projection of shapefile is Swreff99_1800 (EPSG:3011). Now my question is How I get this SRID 3011 by using some OGC method of geometry field?


How do GPS receivers map out directions?

Although I've never really dug into it, I've always assumed mapping out directions in a GPS unit is done with A* or Fibonacci-heaps at a higher level and Dijkstra's algorithm for smaller graphs. Are there any other algorithms being used that aren't immediately obvious?


Location data for Nuclear Power stations in France

As it says on the title, does anyone know a place I can look at and/or download location data for Nuclear power installations in France? I believe it is well known data, but I cannot find it anywhere. It's a very innocuous request, based on choosing where to buy a new pad in France and is nothing sinister. Thanks in advance.


WFS Layer does not refresh or give new data when map moved in openlayers

I am trying to add a WFS Layer with a BBOX strategy. When I change the extent I want GeoServer to return only the features in that extent. Here is the code I am using to add layer.               wfs = new OpenLayers.Layer.Vector("Fields_WFS", {         strategies : [new OpenLayers.Strategy.BBOX()],         protocol : new OpenLayers.Protocol.WFS({         version : "1.1.0",         url : GEOSERVER_HOST + ":" + GEOSERVER_PORT + "/geoserver/wfs",         featurePrefix : "rpid",         featureType : "FIELDS_SUBSET_BNG_OSGB36_1",          featureNS : "<URI>",         srsName : "EPSG:27700" //default         }),     });      The map is set to a maxextent initially. And I have restricted WFS features to max 100 in geoserver settings. My problem here is it only shows the same features every time I load and when I move or zoom +- it does not refresh WFS features. Thanks, JP


Reading Coded Domains from ArcSDE 10 RDBMS using SQL/ADO.net

I need to read ArcGIS Domain Coded Value pairs (value/description) from SQL Server / Oracle in SQL (i.e. not using ArcObjects). This is an XML field. I can get the domain definition thus select gi.definition from dbo.gdb_items gi, dbo.gdb_itemtypes git where git.Name='Coded Value Domain' and gi.Name='MAJORTOWN' (where MAJORTOWN is my domain name)...just looking for the next, XML-parsing part. I want an IDictionary of the value,description. Anyone done this?


How Trusted is ZipInfo.com's data?

I'm not in GIS, and I'm working on a nutty redball express project that involves matching some zip+4 data to congressional districts for the upcoming 113th congress. We have found very few sources of this data as it hasn't been officially published by the census bureau. One of the ones we did find was ZipInfo.com. The problem is I'm not sure if we can trust their data -- I'm not horribly familiar with the field, and my google fu leads to no smoking guns either for or against. So has anyone had experience with their data? Is it worth the price of admission?