How to show labels in arcglobe for extruded features

I am trying to visualize lithology data in ArcGlobe, i.e, the representation of different types of soil data under the earth. This I have done properly by creating Z aware featureclass and extrusion properties. I have attached a picture, a representation in Globe. In this diagram each color in the 3D bars is one feature representing the type of soil. Here I want to show the labels from the attribute table having the depth information of each feature. But when I am applying extrusion automatically, the Show Labels menu gets disabled. Is there any way to add labels programmatically (vertically next to each feature) like 3D text elements etc? Give me some ideas please. Thank you.... ![screen shot of 3D representation](http://i.stack.imgur.com/sYdba.jpg)


IEditSketch and IEditSketch2

What is the difference between IEditSketch and IEditSketch2? Also why there are same interfaces with the same name only with a difference "2" at the end of the interface's name ?


Defining a scale number from the Helmert transformation

I know this is more bit of a cartographic issue, but maybe somebody would be able to help. I have an old map georeferenced to the openstreetmap with 50 identical points. I want to get the scale of the old map. I am using Helmert 4-parameter tranformation (similarity) with the equations:   _X=ax-by+X0_   _Y=ay+bx+Y0_   where _x_ and _y_ are the old map coordinates in meters (same as coordinates used from openstreetmap), _X0_ and _Y0_ are translations and then there is one scaling and rotation. Scaling parameter _m_ could be obtained as:   _m=sqrt(a^2+b^2)_   Rotation is:   _alpha=arctan(b/a)_ The resolution of the old map is 400 dpi (0.0635 mm per pixel). I have all the parameters counted now with the least squares method.   Does anybody know how to get the scale of the old map from these parameters? I want the scale in a common form 1 : M (1 : 600 000).   _EDIT_ : Further on I will use iteratively reweighted least squares and I will be omitting those identical points with the biggest spatial error (I am working with very old maps which have mistakes). I will do this to make the definition of scale more accurate. That's why I need to use the transformation parameters for the map scale definition and not only to compare the distances.   Thanks a lot.


pivot table output

In an attempt to summarise landclasses, the columns of my pivot table output are a bit mixed up with the field info regarding the three input fields that I've selected. As shown in snapshot 3 the description field have the info of the grid_slo field, the grid_slo field have the info of the grid_alt field, and grid-alt showing zero. The original descriptions that should show under the descriptio field has been lost. Numerous attempts to fix this problem has been unsuccessful. Can anyone help please? ![attribute table](http://i.stack.imgur.com/3ta5U.png) ![pivot table input](http://i.stack.imgur.com/gueFu.png) ![snapshot showing the output, with the description field having the info of the grid_slo field, the grid_slo field having the info of the grid_alt field, and grid-alt showing zero. The original descriptions that should show under the descriptio field has been lost](http://i.stack.imgur.com/3nwD5.png)


Missing date column from interactivity callback

I can't seem to get "date" columns to return for popups/infoWindows. I can retrieve other data types without an issue. I have the interactivity set correctly from what I can tell, but the featureClick callback only returns one column. My code relies on v2 of the CartoDB API. Here's my call to cartoDB.createLayer: cartodb.createLayer( map, { type : "cartodb", options : layerData } ) where layerData is               interactivity: "name,admin_date,size"     query: "SELECT size, name, admin_date, the_geom_webmercator FROM book1"     table_name: "book1"     tile_style: "#book1{marker-width:5;marker-opacity: 0.9;marker-line-color: white;marker-fill:#e6e6e6;marker-line-width: 1;marker-line-opacity: 0.9;marker-placement: point;marker-type: ellipse;marker-allow-overlap: true;} #book1{[size>=0]{marker-fill:#b8cde9;}[size>=18229430]{marker-fill:#6693cc;}[size>=22061458]{marker-fill:#2f5f9d;}}"     user_name: "c2es"      The featureClick callback only return the first and last columns. Any ideas what's wrong?


How to close spaces within a polygon in QGIS 2.2?

I have a shapefile displaying the borders of several city units. One of its polygons, however, had free space consisting of four vertices in it. I created another polygon inside the free space and then used `Dissolve` to close the free space. The result was the following layer which still contains some errors (see the close-ups below). What is the best way to resolve this in QGIS 2.2? Using some kind of `v.clean`? Thank you for your help. ![enter image description here](http://i.imgur.com/lqqUSoe.jpg) ![enter image description here](http://i.imgur.com/GtecT6e.jpg) ![enter image description here](http://i.imgur.com/tKzyup1.jpg) ![enter image description here](https://i.imgur.com/dkLf9h0.png) ![enter image description here](https://i.imgur.com/FxdGrAO.png) ![enter image description here](https://i.imgur.com/tT12UPE.png) Edit: The entire data seems to be corrupted. Using `Dissolve` on the entire shape results in what the fourth screenshot shows. Edit 2: Added fifth screenshot. Edit 3: Added sixth screenshot.


QGIS export ESRI Shapefile to SpatiaLite format issue

I am trying to export a shapefile of 3221 mixed POLYGON and MULTIPOLYGON geometry elements and every time, I get this same error: ![enter image description here](http://i.stack.imgur.com/Bv6K2.png) Does anyone have any idea what's going on? And no, this is not an exact duplicate of http://stackoverflow.com/questions/15203403 as I am working solely within QGIS. Any help is greatly appreciated. Thank you! - Brandon


Is there a way to get centroid of geographic object in MSSQL?

I have a table of OSM geographic roads in MSSQL 2012. I want to create a new column with the centroid of each road (So I can set center on OpenLayers map for each road). How can I do it ?


How to get the best accuracy of coordinates for web GIS- application?

I would like to get advice from you about one thing. I would like to develop my own web GIS application. It would be about some sport places in my hometown; it would be good practice for my education - I'm still a student. First, I need to get coordinates of the features, then I think to put it in AutoCad, then export to .shp file, then work with it in QGIS and use Postgres and PostGis for databases and than use Mapserver and so on. But I need help with the first step, getting the coordinates. I don't have any high accurate instrument like some good GPS, so is there any chance to get coordinates of the places by Google Earth? How accurate is that and is it even possible to do it like that? I'm looking forward to someone's answer and please don't get mad if I said something stupid because as I said, I'm in the beginning of getting to know all this GIS web development and things like that. Thanks! * * * Thank you all for your answers, it all helped me in some way. The only thing that I'm not sure is accuracy of osm and bing maps. I was editing some stuff in JOSM recently and I've realized that bing maps don't really fits in the right place if you compare it to gpx data. (there is at least 2m difference) I'm not still quite sure is it accurate enough for me. I would like to achieve 1m maximum level of accuracy.. :/


v.dissolve to multipolygon

Using v.dissolve in GRASS, I'd like to get a multipolyon as a result. However, I end up with distinct polygons with the same 'dissolve value.' For example; lets assume I have two polygon features which are not touching each other spatially. Both have a 'DISS" value of 1. If I run v.dissolve on this DISS value, I get the same thing in the output as the input (i.e. two seperate polygons with a DISS value of 1). What I would like to be able to get is an output that contains one single multipolygon. How do I go about this?


Will QGIS fulfill my requirements? moveable, rotatable, draggable, complex, interactive icons

I'm trying to compile a list of mapping software that will fit my requirements for a desktop application I want to build. I'm trying to build something similar to MissionPlanner: https://code.google.com/p/ardupilot- mega/wiki/Mission Here is a quick moc-up I did for my application: http://i.imgur.com/52DaJtQ.png Requirements:   * I want to be flexible with what I can draw on the map: lines, rectangles, polygons    * I want to be flexible with icons: they should be able to move, rotate   * I want flexibility with icon interaction: for example, being able to easily drag a waypoint would be great   * I want to be able to use my own data for maps: for example an aeronautical chart made from a folder of .tiff images   * I want to sell my application commercially   * I don't have a preference on programming language yet Is QGIS right for me?


How to create a catchment area using LSOA and point file data

I have to create a catchment area (based on LSOA) for the 19 underground stations within Camden borough and it has to be in a polygon vector file. I know I have to somehow group the SOAs but don’t know which geoprocessing tool to use.


Arcpy CalculateField using formatting syntax

I have 2 feature classes (FC) that I am merging together and then performing some field calculations on. The 1st FC has 100k+ records in it, and the second FC is a blank FC that I'm using as a template. I need to merge the two FC's together and then field calc the values from FC1 into FC2.   The problem with FC2 is that the fields are required to often be of a different type. For example, in FC1 field "PID" is a double and in FC2 "PIN" is a string. I need to calculate the values of "PID" into "PIN" with a prefix of "139-". This pattern is repeated several times with different fields. My simple code example:               data = r"C:\Temp\Test.gdb\PAR_PROPINFO_SCO_VW2"         template = r"C:\Temp\Test.gdb\Template"         out = r"C:\Temp\Test.gdb\Scott_CO"         arcpy.Merge_management([data,template],out)         qry= '"139-"' + "+'PID'"         arcpy.CalculateField_management(out,"PIN",qry)      In this particular instance, qry threw a syntax error. If I just use:               arcpy.CalculateField_management(out,"PIN","PID")      I get a type mismatch error. Since I'm on Arc 10.0, I was trying to avoid a cursor with the amount of records involved. This process will run daily and the # of records will only increase. I was thinking I could just get away with a brute force merge and calc method. However, I run into type mismatches. I thought field mappings, but that is beyond my current level.


Unusual results for geoprocessing speed tests

I have been observing unusual performance with a Python geoprocessing script. The (attached) script performs the following actions:   1. Use a search cursor to look up UTM zone corresponding to polygon features   2. Create spatial reference object based on search cursor results   3. Convert .csv to feature layer and then to a point feature class I've noticed markedly different processing times based on how the script is run:   * **32-bit** processing using **IDLE** = 203 seconds   * **32-bit** processing foreground **script tool** = 91 seconds   * **64-bit** processing background **script tool** = 206 seconds _Why would this script perform so differently given the above conditions?_ I certainly would not expect the 32-bit script tool run in the foreground to be 2X as fast as the other methods. * * *               import arcpy, os, time          ###IDLE Parameters     ##fc = r'C:\path\to\polygon\fc\with\utm\zones\and\features'     ##outws = r'C:\out\location'     ##arcpy.env.workspace = r'C:\workspace'          ####################     ## Script tool parameters     fc = arcpy.GetParameterAsText(0)    # Feature class     outws = arcpy.GetParameterAsText(1) # Folder     arcpy.env.workspace = arcpy.GetParameterAsText(2)   # Workspace     ####################          # Tables are .csv     tables = arcpy.ListTables()          start = time.clock()          # Look up which UTM zone .csv features are in     for t in tables:         quad = t[7:17]         print quad         whereClause = """ "QUADID" LIKE '%s' """ % quad         with arcpy.da.SearchCursor(fc, ("QUADID","ZONE"), whereClause) as cursor:             for row in cursor:                 if row[0] == quad:                     utmZone = row[1]                     if utmZone == 10:                         sr = arcpy.SpatialReference(26910)  # NAD_1983_UTM_Zone_10N                     elif utmZone == 11:                         sr = arcpy.SpatialReference(26911)  # NAD_1983_UTM_Zone_11N                     elif utmZone == 12:                         sr = arcpy.SpatialReference(26912)  # NAD_1983_UTM_Zone_12N                     elif utmZone == 13:                         sr = arcpy.SpatialReference(26913)   # NAD_1983_UTM_Zone_13N                     else:                         print "The UTM Zone is outside 10-13"                 else:                     pass              # Convert .csv to feature class         try:             outLayer = "in_memory"             # Now with the sr defined, create the XY Event Layer             arcpy.MakeXYEventLayer_management(t, "x", "y", outLayer, sr, "z")             arcpy.FeatureClassToFeatureClass_conversion(outLayer, outws, t[7:17])             arcpy.Delete_management("in_memory")             end = time.clock()             print "In_memory method finished in %s seconds" % (end - start)              except:             # Print any error messages             print arcpy.GetMessages(2)          print "Processing complete"


how to open a layer's table in a map document

Trying to create a map document programmatically. Added a featureclass as a layer. When I open the map document in ArcMap, I want the table for that layer to be open. Found `TableView` (doc) and `AddTableView` (doc), but these seem to be only for standalone tables. I also looked at the properties for Layer (doc) and for DataFrame (doc) and didn't see methods that looked right. What I'm shooting for is to double-click on the mxd and have something like image below appear. ![enter image description here](http://i.stack.imgur.com/jaEzR.png) Any ideas? Thx


How to create a multi-colored raster image based on long,latt,color information?

I want to create a multi-colored raster layer to overlay on my map in OpenLayers I have the following information in my DataBase:               D_Cell_ID          Long            Lat         ColorCode     055113             121.4698        40.8417     FF6600     057113             121.4698        40.8325     FF6600     .     .     .     063125             121.0344        40.4658     FFCC00     063127             121.0344        40.4567     FFCC00      I would like to have the output something similar like this:Example ![raster overlay](http://i.stack.imgur.com/I3Kye.jpg) Hence I need to write some code to create a raster image based on long,latt,color information. Do anyone know how to do this? or any library tool which can import into Java environment to do this? Method and direction is needed. Thank you!


How to specify fields to export from QGIS?

I'm using QGIS 2.0.1 and am connected to a MS SQL Server spatially enabled database (not ESRI). I am able to connect to the layer and load it into my TOC just fine. I have hidden a number of fields that are not important for visibility sake. Now I want to export that data to another data type (we'll use shapefile for this example). When I right-click the layer, choose "Save As" and define the output file, I get the following error: > Export to vector file failed. Error: unsupported type for field GlobalID _**Note 1** : If I check "Skip attribute creation" in the dialog box, the shapefile is created with no errors._ _**Note 2** : ArcGIS does not have this problem, it just converts the GUID field to "TEXT" type._ Looking at the table, one of the fields that I have hidden is a GUID of SQL type <uniqueidentifier>, which is in all likelyhood the culprit. In an attempt to resolve this, I made the field "hidden". I would have assumed that because the field was hidden, QGIS would not try to export this field, but that turned out to not be the case. **Can you define the fields you want to export from QGIS?** Is it possible for QGIS to export just the "Visible" fields? Or is there another alternative way of doing this from QGIS (maybe define the fields via an ogr2ogr query)? This seems like a shortcoming of QGIS if I cannot define which fields to export. **UPDATE:** I did a test using a sample shapefile to see if I could limit the output fields. The following ogr2ogr command line worked and eliminated all fields except for "name":               ogr2ogr -f "ESRI Shapefile" test.shp city.shp -sql "select name from city"      but using QGIS and adding:               -sql "select name from city"      to the "layer" section of the "Save As" dialog did not result in a shapefile with just the single field. So I don't really see how I can manipulate the file output to show just a select few fields. Any help is appreciated.


draw multipart polygon with one surrounding border

As you can see on the image below, the waterbody is divided by the border line because it consists of several polygons. ![enter image description here](http://i.stack.imgur.com/73vLm.png) I could merge the polygon, but I dont want to change the shapefile. Is there a option in QGIS to draw te waterbody just with a border sourrounding its "absolute extend"?


How to filter big datasets in Geoserver WMS?

Is there a good way to filter datasets with approx 50,000 features in Geoserver? I'd have to select between 100 and 200 features based on ids. These ids are different from the FeatureIds. Couldn't find a "IN" operator similar to SQL in OGC Filter Encoding Standard. Combining multiple PropertyIsEqualTo operators performs really poorly.


Cannot run the Arc samples

When I try to compile I get Error 1 Cannot register assembly "C:\Users\Patrick\Desktop\DynamicBiking\CSharp\bin\Debug\DynamicBiking.dll" - access denied. Please make sure you're running the application as administrator. Access to the registry key 'HKEY_CLASSES_ROOT\DynamicBiking.BikingTrackModeCmd' is denied. DynamicBiking2010 Trying to compile this one: http://help.arcgis.com/en/sdk/10.0/arcobjects_net/SampleZips/101efef5-9071-49ac-9ce9-c5a5fe5cc896CSharp.zip I have administrator priveleges


using shapely: translating between Polygons and MultiPolygons

[EDIT: the solution to this was simply to use OGR to read shapefiles. See geographika's example.] In an ESRI shapefile, there is no distinction between Polygons and MultiPolygons. Furthermore, there is no explicit distinction between interior holes and exterior rings (besides the "handedness" of a given polygon). So after reading a shapefile, I have a list of coordinate sequences describing rings, but without some more intensive processing, I cannot distinguish which of these rings are exterior rings, interior holes, or additional polygons. It appears that for shapely's Polygon and MultiPolygon constructors, there must be a clear distinction between exterior and interior rings, so how should I move from an unclear list of rings to an ordered set of separated polygons, with clearly designated interior and exterior rings? To summarize: if I have a list of polygon rings, but I don't know which rings are holes in the interior or are separate polygons, how should I best sort them into separate polygons with designated interior holes? I'm looking for a simple algorithmic solution that I can implement in python, can be used to process hundreds of polygons in ~a minute or less, and I'm doing this in order to perform a large number of intersections.


Counting Events Per Polygon Feature

I would like to calculate the total number of events within each feature of a polygon using a dbf file and a polygon shapefile. Is there an easy way to automate this so I can then produce charts per feature area? Thanks, Haskett


Out of range error when converting from geometry to geography in PostGIS

I am attempting to convert a geometry column in the postgres database to a geography column. I created a new column like this:               alter table scene add column footprint_new geography(POLYGON, 4326);      and I'm trying to populate the new column like so:               update scene set footprint_new = footprint::geography;      and I get the following error:               ERROR:  Coordinate values are out of range [-180 -90, 180 90] for GEOGRAPHY type      Now, it isn't too surprising that there are some out of range values in this database. I want to skip those out of range geometries during the update and identify them later so they can be fixed. How can I do that?


How to use more than one CPU?

I'm creating a real time traffic application. I read the AVL car data from a web service and insert the records in a PostgreSQL table Then a trigger calculates the closest link and if there is a previous position for the same car I calculate the route to get the distance travelled in that time frame, then I get the average speed for that route. The thing is I get around 3000 points each 5 min, so I need to process ~600 records every min or 10 each second. I'm very close to those numbers because I filter the road network before calling the routing function and the routes are very small. I understand PostgreSQL doesn't support multithread in one connection but support multiple connections so could use the other CPU cores. I can make a different table for each State (25), but how can I partition the problem so each trigger creates a new connection and use more CPU cores instead of only one?


Using primary keys with ArcSDE (Oracle)?

Is it possible, more advisable to assign a primary key to the OBJECTID column in a featureclass in SDE. I am using Oracle and some of my flat tables (junction tables) need to SEE a primary key in the featureclass table in order to model a M-M relationship correctly. By default ,no primary keys are made by Oracle when you insert and SDE layer (which makes sense as you are doing this all through ESRI tools anyway) although I know SDE creates an index in SDE, but I did not want to cause any kind of corruption if I start meddling with things on the Oracle side - is there any danger of this? Before submitting this question I also did a little more reading around the subject, and I seem to be finding that it is bad practice to use OBJECTID as a form of primmary key in any case - so would I be better off creating my own primary key field and then let Oracle maintain this with triggers and sequnces? If I do this, effectively behind SDE's back, will I upset SDE and all it's metadata etc? **UPDATE** I have found this nugget of info from Neil Clemmons - what do you think? > Using the OID may not be a good choice. For instance, let's say you have a > road centerline feature class related to a maintenance table using the OID > field. If you decide to merge two centerlines into the same feature, ArcMap > does this by deleting the two original features and creating a new one. > Depending on how you have the relate setup, the records in the related table > could be deleted. At the very least, the new feature would not be related to > any records in the table because it has a new OID. If you based the relate > on another unique field, the new feature would retain the attributes of one > of the original features and the relationship with the table would be > maintained. and more Neil Clemmons advice... > You shouldn't use the OID as a foreign key. This field is created and > handled by the ArcGIS software and is present so that a unique identifier is > always available for the ArcObjects functions that need one. You have no > control over the value of this field. If you need such a field for use as a > foreign key, you should create one and populate it yourself. This way, the > value will stay with a feature throughout the lifetime of that feature.


ArcObjects ITextSymbol.GetTextSize in Python revisited: usage with StartDrawing and FinishDrawing methods

I am using ArcGIS 10.0/ArcObjects and Python 2.6. I am trying to get the usage of the ITextSymbol.GetTextSize method correct. The syntax, as established in an earlier question, is:               dblx, dbly = pTextSymbol.GetTextSize(pDisplay.hDC, pDisplay.DisplayTransformation, txt)      and the outparameters are returned just fine. My problem is in figuring out how to use this method so that it does not screw up ArcMap's screen cache. ESRI documentation indicates that it should be surrounded by calls to IDisplay.StartDrawing and IDisplay.FinishDrawing so that the input hDC will be correct. (Indeed, calling GetTextSize without the start and finish drawing results in an invalid hDC.) The syntax that I am using for StartDrawing is as follows:               pDisplay.StartDrawing(0, esriDisplay.esriNoScreenCache)      I have tried both 0 (as discussed in documentation) and pDipsplay.hDC for the first parameter of StartDrawing, but the ArcMap screen still gets stuck in one place. (Panning just causes a redraw from cache right back to the same place.) Either input for the first parameter results in valid outparams but a stuck screen. The problem could also be with the esirNoScreenCache constant, I suppose, but I don't see how. I am not trying to draw anything to the screen at this point. I'm just operating on a TextSymbol, not a TextElement or the like. pDisplay and txt are both valid, in-process objects, so the problem, I think, is with the hDC or the screen cache. I'm at a loss to figure out which. Any insights are greatly welcomed.


Spatial Classification

I am not sure what to call my problem. From this image, I have a set of classified polygons (colored). I would like a way to classify the white polygons in the image. I don’t have any attributes to base the classification on other than proximity to a classified polygon. I have tried a naïve approach by using the near tool in Arc. This works fine enough but it becomes sensitive to outliers. For example, a set of unclassified polygons may be surrounded by polygons of category "A" but if there is a single “island” polygon of category "B" nearby, they’ll be classified the same as the island polygon. I am looking for a way to rank the likelihood that an unclassified polygon belongs to a particular category based on proximity to arbitrary k-polygons. For example, if there is a polygon of category "B" immediately adjacent to an unclassified polygon but the next five closest polygons are of category "A", I'd like the system to rank the unclassified polygon as category "A". I am wondering if there’s anything in R that would solve my problem. Really, any off the shelf method would work (python, qgis, postgis, arc). ![enter image description here](http://i.stack.imgur.com/ZBg9e.jpg)


Unable to Install External QGIS Plugins

I'm using: QGIS version 1.8.0-Lisboa on Windows XP, installed using OSGeo4W I am behind a proxy but I have configured it in Options > Network. When I go to Plugins > Fetch Python Plugins.. All of the the Repositories are "Status: Connected" however on the Plugins tab I can only see the 6 installed plugins (Repository: only locally available) - See images below: ![Repositories ](http://i.stack.imgur.com/l1qgS.png) ![Plugins](http://i.stack.imgur.com/Y67kx.png) Any ideas why I cannot see any plugins from the "connected" repositories? Many thanks. Update - As suggested by vinayan I can install plugins by manually downloading them and putting them in the relevant folder and using the QGIS Plugin Manager. However still no joy using the Fetch Python Plugins.


issue with WEB ADF based Web application

I have a WEB ADF based application that I am having some issues with. When I start app in debug mode, it shows me a blank page but I can right click and choose the option view in browser and it works fine and loads properly. I am trying to figure out what I need to do to load this page propely in a debug mode. Appreciate any ideas. Thanks Jay


Does anyone know a good way to handle z ordering when using GraphicTracker in ArcEngine?

From the GraphicTracker Documentation: > The first graphic added will appear underneath subsequent graphics if they > have overlapping display areas. This is not true when using ArcEngine 10.0 with polygon fill symbols on picture marker symbols - sometimes it works (1st time for sure) and sometimes it does not. Has anyone figured out a proper way to manage z ordering of symbols when using ArcEngine and using GraphicTracker?


Polygon shapefile with null values

Does anyone know where I can find a sample shapefile that has null entries in it (preferably a polygon shapefile)? I'm not sure how to create one with null entries. Thanks, Jeff


PostGIS - Intersects feature to split into linestring

I would like to know if it's possible to intersect feature datasets (multipolygon, polygon, multilinestring and linestring) to get a result as a set of linestrings being split at each intersection point. ![enter image description here](http://i.stack.imgur.com/kR4n5.png) Output should give:   * 4 linestrings for the square    * 3 linestrings for the line   * 2 linestrings for the polygon


Lost qgis_core.dll

I´ve been working with QGIS 2.4 since two months ago. After doing some changes in my computer, I think I've lost something important for running QGIS 2.4 The message is: > qgis.exe - System error. > > The program can´t start because qgis-core.dll is missing from your computer. > Try reinstalling the program to fix the problem I reinstalled the 2.4 32 and 64, the 2.2 version and the problem is still there. I can't open QGIS. It shows the first window and crashes. The QGIS browser works fine.... This is the path of my computer: > C:\Program > Files\Gimp-2.8\bin;C:\Python34\;C:\Python34\Scripts;C:\Python34\Lib\site- > packages;C:\Program Files (x86)\Intel\iCLS Client\;C:\Program > Files\Intel\iCLS > Client\;%SystemRoot%\system32;%SystemRoot%;%SystemRoot%\System32\Wbem;%SYSTEMROOT%\System32\WindowsPowerShell\v1.0\;C:\Program > Files\Hewlett-Packard\SimplePass\;C:\Program Files (x86)\Windows > Live\Shared;C:\Program Files (x86)\ATI Technologies\ATI.ACE\Core- > Static;C:\Program Files\Intel\Intel(R) Management Engine > Components\DAL;C:\Program Files\Intel\Intel(R) Management Engine > Components\IPT;C:\Program Files (x86)\Intel\Intel(R) Management Engine > Components\DAL;C:\Program Files (x86)\Intel\Intel(R) Management Engine > Components\IPT and the python path is: > C:\Python34\Lib;C:\Python34\DLLs;C:\Python34\Lib\lib- > tk;C:\Python34\Scripts;C:\Python34\Lib\site-packages;C:\Python34\projects I don´t know what to do to fix the problem and I can´t find anything on the web.


How to measure correct distances in a project with OpenLayers background?

I have 3 lat/long points as follows:               A = {35.257246,139.721256}     B = {35.256541,139.721491}     C = {35.256326,139.7203271}      * * * 1) Measurement result by Movable Type Scripts: By using the distance calculation method at Movable Type Scripts, the distance results are:               distance(A,B)= 80.88745771063131(m)     distance(B,C)= 101.58773310547808(m)      I believe these results are correct. * * * 2) Measurement result by QGIS measure line tool: The QGIS desktop 2.0.1 and the Projection CRS is set to `[+proj=longlat +datum=WGS84 +no_defs]`. I use the Plugins/OpenLayers plugin/Add Google Satelite and import 3 lat/long points into QGIS. The results are:               distance(A,B) = 100.374(m)     distance(B,C) = 135.114(m)      I need to generate a map from a large set of lat/long points as above and Google Satellite background on QGIS. Could you show me how to solve this problem?


Error during visualization of many layers on Chrome (ExtJS, GeoExt, OpenLayers)

I've noticed that on the latest version of Chrome (36.0.01985) my application crashes when a lot of layers are displayed (and/or the browser is maximized and the map is zoomed out). In particular, the panels are not shown at all and the map is only partially rendered. No matter the kind of layer (raster, vector), the problem is always present when there is a lot of data to render. The problem doesn't occur on other browsers. I've also tested it also on an older version of Chrome (26) and interestingly, there was no error. I'm actually using ExtJS 3.4, GeoExt 1.1 and OpenLayers 2.13.1 (I know they are oldies, but till now they worked very well). I'm also wondering if an upgrade of these tools could solve the problem or only make me feel younger. Tips?


Python Syntax Error when running Python Snippet of successful geoprocess

I am having an issue with a simple python snippet. Here's the background: I wanted to include an additional step in a nightly update script that replaces characters in a field that simply do not play nice with a web application that allows users to query the data. These characters included angle brackets ("<" and ">") which I was trying to replace with "Less Than" and "Greater Than". So when I was testing the Field Calculation, I finally got it to work the way I wanted to, and as I have done with so many other snippets, I simply copied the Python Snippet of the geoprocess, and pasted into my script. Below is the snippet:               arcpy.CalculateField_management("path/to/my/layer","field","RemoveBad( !field!)","PYTHON_9.3","def RemoveBad(x):/n   x = x.replace(\"<=\", \" Less Than or Equal To \" )/n   x = x.replace(\">=\", \" Greater Than or Equal To \")/n   x = x.replace(\"<\", \" Less Than \")/n   x = x.replace(\">\", \" Greater Than \")/n   x = x.replace(\"=\", \" Equal To \")/n   x = x.strip()/n   return x")      Now in my script, I do have a log that reports geoprocessing errors, so when I run the script, I am getting the following:               Executing: CalculateField path\to\my\layer field RemoveBad( !field!) PYTHON_9.3 "def RemoveBad(x):/n   x = x.replace('<=', ' Less Than or Equal To ' )/n   x = x.replace('>=', ' Greater Than or Equal To ')/n   x = x.replace('<', ' Less Than ')/n   x = x.replace('>', ' Greater Than ')/n   x = x.replace('=', ' Equal To ')/n   x = x.strip()/n   return x"     Start Time: Mon Aug 25 09:35:32 2014     Failed to execute. Parameters are not valid.     ERROR 000989: Python syntax error: Parsing error SyntaxError: invalid syntax (line 1)     Failed to execute (CalculateField).     Failed at Mon Aug 25 09:35:32 2014 (Elapsed Time: 0.00 seconds)      So I started combing through the script but I cannot find the issue, which leads me to believe it has something to do with the angle brackets. To reiterate, this python snippet is from a successful Field Calculator geoprocess. Any assistance in pointing out where the syntax error is occurring would be appreciated. **UPDATE** : Per Alex's feedback I created the following variable:               RemoveBad = """def RemoveBad(x):        x = x.replace(\"<=\", \" Less Than or Equal To \" )        x = x.replace(\">=\", \" Greater Than or Equal To \")        x = x.replace(\"<\", \" Less Than \")        x = x.replace(\">\", \" Greater Than \")        x = x.replace(\"=\", \" Equal To \")        x = x.strip()        return x"""      Then I updated the snippet:                arcpy.CalculateField_management("path/to/my/layer","field","RemoveBad( !field!)","PYTHON_9.3",RemoveBad)      The error is now gone and script is working as expected.


Using auto increment code with counter

I was wondering if it's possible to use the standard auto increment code to work on groups of records:               FeatureClass = "newFC"      recs = collections.defaultdict(int) autoIncrementEXP = "autoIncrement()"      expCalc = """def autoIncrement(!Group!): recs[!Group!] += 1 return recs[!Group!]"""     arcpy.CalculateField_management(FeatureClass, "Number",autoIncrementEXP,"PYTHON",expCalc)      So instead of having the auto increment portion of the script stop when it reaches the end of the records, have it create sequential numbers using field value groups if that makes sense. You would start with something like the `Group` field below and end with the `Number` field containing the increment values:               Group   Number     A       1     A       2     A       3     B       1     B       2     C       1     C       2     C       3     C       4     C       5      Thanks for any insight you can provide.


Can ArcGIS Desktop 10 Add-in Use IronPython?

Coming from the same background as the original poster in this thread, I wonder if anyone has had any success creating an Arc10 Add-In using IronPython. I have been able to access ArcObjects via IronPython using the info in the linked thread, assorted blog postings, and a few posts on the ESRI forums. These neophyte attempts have all been run from the command line and I would to be able to gain access through either a toolbar or a button. One thought that I had was to create the add-in using the .Net 3.5 template, add a python file, and then call a function in the python file from the button's C# file. It appears though that add-ins are not yet supported in .Net4? Here is the code from the .cs file which controls button functionality.               using System;     using System.Collections.Generic;     using System.Text;     using System.IO;          //I need to then add references to IronPython     using IronPython;     using IronPython.Hosting;     using Microsoft.Scripting;     using Microsoft.Scripting.Hosting;          namespace PyButtonTest     {         public class ZoomToLayer : ESRI.ArcGIS.Desktop.AddIns.Button         {             public ZoomToLayer()             {             }                  protected override void OnClick()             {                 var ipy = Python.CreateRuntime();                 dynamic test = ipy.UseFile("ZoomToLayer.py");                 test.ZoomToActiveLayerInTOC(); //Call this function in the py file                      ArcMap.Application.CurrentTool = null;             }             protected override void OnUpdate()             {                 Enabled = ArcMap.Application != null;             }         }          }      As far as I can tell the Add-in format takes care of getting the application hook. I should then be able to, in the python file, after referencing and importing all of the necessary libraries get the map document with (I think?)               mxDoc = ArcMapUI.IMxDocument.ActiveView      From here I am stumped. Google has not turned up anything helpful. Has anyone else tried this? Got something to work? Do I need to wait until add-ins are supported in .Net4 to go this route?


QGIS Contour plugin: Creating complex isolines leads to isolated areas

I use the contour plugin (v1.1) within QGIS 2.4 to create contour lines. I have streets with different travel costs. I observe problems for streets with low travel costs if they are surrounded by streets with high travel costs. The low cost street is "piercing" into the low speed area (cf. picture). The points are created as the travel cost markers from Dijkstra algorithm using pgrouting 2.0.5. The Dijkstra algo. calculates correct values for every street. The problem I observe is now, that due to low speed areas left and right to the highway, the isolines are connecting both areas over the highway (yellow). Based on values, the correct behavior would be the one shown in green. **Can someone recommend me a way of preventing such results?** Before I used GDAL TIN interpolation and contour mapping, which leads to better results if it comes to complex networks. On the other hand it adds transformation errors (rasterizing) and extra steps to the calculations. ![Section of my map showing two 60 sec. isolines without connection. Green: The likely correct isoline. Yellow: The wrong connections for the isolines.](http://i.stack.imgur.com/iyVhw.png) ![Original result](http://i.stack.imgur.com/x0yki.png)


Counts the number of points in a polygon in R

I have two classes sharing the same CRS (Latitute and Longitude):   1. `bolognaQuartieriMap`: a `SpatialPolygonDataFrame` containing data of a city boroughts.   2. `crashPoints`: a `SpatialPointsDataFrame` containing data of accidents. They are well plotted using:               plot(bolognaQuartieriMap)     title("Crash per quartiere")     plot(crashPoints, col="red",add=TRUE)      What I need is to get the number of points (`crashPoints`) in each polygon that constitute `bolognaQuartieriMap`. I was suggested to use `over()` but I did not succeed.


OpenLayers: Switch layer off when switching another layer on

As the title says, I want to turn layers off when turning another layer on. I wrote a function that works. However, I have to click the layer twice in the LayerSwitcher to turn the layer on that I want to see. The first click turns the other layers off but only a second click turns the actual layer tho show on. Can somebody tell me where I am wrong in my code?                   function getVLayers(showLayer){             for (var i=0; i < map.layers.length; i += 1){                 if (map.layers[i].visibility){                     if (map.layers[i].name == showLayer && map.layers[i].name == "OpenStreetMap"){                          } else {                         map.layers[i].setVisibility(false);                     }                 }             }         };         map.events.register("changelayer", this, function(e){             getVLayers(e.layer.name);         });


QGIS GRASS Missing Dependancy

I want to perform a view shed analysis in QGIS using GRASS 7 r.viewshed. When I double click on r.viewshed in GRASS GIS 7 commands in the processing toolbox of QGIS I get the following message: > _Missing dependency.This algorithm cannot be run :-( This algorithm requires > GRASS GIS 7 to be run. Unfortunately, it seems that GRASS GIS 7 is not > installed in your system, or it is not correctly configured to be used from > QGIS Click here to know more about how to install and configure GRASS GIS 7 > to be used with QGIS_ GRASS commands seem to work such as r.los but the GRASS 7 commands don't. I am using a MAC OS X Mavericks with the latest version of QGIS 2.4.0 Chugiak.


raster2pgsql does nothing

I try to load rasters through raster2pgsql. I try to call it through pgAdmin command line, but when I call any `raster2pgsql` command (even `raster2pgsql -G`), it does nothing, I just get another line to write more commands. No output, no time doing something. Non-gdal psql commands work as expected. `\raster2pgsql` leads to "invalid command" error. I am new to the command line and psql, so I'm stuck here. I have PostgreSQL 9.3.4, PostGIS 2.1.3 and pgAdmin III 1.18.1 on WIndows XP 32 bit. I made sure that `raster2pgsql.exe` is in its place in `bin` folder, but I don't know of other ways to validate whether it is installed properly.


How to include data from other tables in the info-window in CartoDB

I have 3 layers(table_1,table_2,table_3) of data. In the info-window, which was created for the second layer. I would like to display one column of data from table_1 . When I used `{{table_1.name}}`to display the data, I am getting the error , `"Error in line 40: Cannot read property 'table_1' of undefined` Is there any way to display data from one table in another layer of data?


Formatting Cells in ArcMap

I am trying to join tables together and need the PLSS to line up correctly. In my table with geometries I have "CO060020N0650W0SN330" in a field. I want to make a new column and make it say 2N65W33 so I can join a table. Is there a formula I can use in the Field Calculator to get this result?


When styling with data defined properties or Labeling, Is there a way to replace an ID representing field with the actual values depicted by the ID?

I am currently using data defined properties and labeling to present pipes based on diameter. I now have a dictionary table that has id's for every combination of diameter and material and I'm using it with the value relation widget. I must find a way to access the value represented by the key in the diameter || material field. An ugly way that solves this would be to use:               CASE        WHEN id_material_diam in (id1, id2) THEN result_formula_containing_diameter        WHEN id_material_diam diameter in (id3...) THEN ...     END      Is there an elegant way that I can access the actual values for which the ID is the key? Thank you PS. QGIS 2.2 Valmiera


GDAL equivalent of ArcGIS Mosaic

The word 'mosaic' is often used for quite different things within the worlds of ArcGIS and GDAL. With GDAL and it's downstream applications like Qgis it is common to see mosaic used in the context of _"assemble all these tiled images so they appear as one big image"_ , with gdalbuildvrt and `gdalmerge` put forth as a convenient means of getting there. ArcGIS has the same starting point, _"assemble all these tiled images so they appear as one big image"_ , and then goes on to add _"while doing so [blend, min, max, mean, sum, ...] the overlapping pixels"_ and a few other things. To put in uncharitably, in this regard ArcGIS is smart and GDAL is stupid. However it's pretty evident GDAL is _not_ stupid, so my question is: in the GDAL and/or Qgis tool chain what is the equivalent workflow to ArcGIS mosaic? _(note, I'm not talking about the Mosaic Dataset, which is a different beast altogether)_


R - Dealing with missing data in SpatialPolygonDataFrames for Moran.test

I'm trying to use Moran.test on a SpatialPolygonDataFrame consisting of 7194 elements in R. I know that there is around 150 polygons with NA values. First I generate a spatial weights matrix: `W_cont_msoa <- poly2nb(all.data, queen=T) W_cont_msoa_mat <- nb2listw(W_cont_msoa, style="W", zero.policy=TRUE)` Then when I carry out Moran.test the code below returns:               > moran.test(residuals, listw=W_cont_msoa_mat, zero.policy=T)     Error in moran.test(residuals, listw = W_cont_msoa_mat, zero.policy = T) :      objects of different length      I've also tried to this but it hasn't worked:               > moran.test(residuals, listw=W_cont_msoa_mat, zero.policy=T, na.rm=TRUE)     Error in moran.test(residuals, listw = W_cont_msoa_mat, zero.policy = T,  :      unused argument(s) (na.rm = TRUE)      I've tried to go back to the original data and remove missing data using complete cases but I received the following:               > all.data[complete.cases(all.data)]     Error in all.data[complete.cases(all.data)] :      error in evaluating the argument 'i' in selecting a method for function '[': Error in       complete.cases(all.data) : invalid 'type' (S4) of argument      The structure of the data is as follows:               `Formal class 'SpatialPolygonsDataFrame' [package "sp"] with 5 slots     ..@ data       :'data.frame':   7194 obs. of  23 variables:     .. ..$ MSOA01CD   : Factor w/ 7194 levels "E02000001","E02000002",..: 1 2 3 4 5 6 7 8    9 10 ...        .. ..$ MSOA01NM   : Factor w/ 7194 levels "Adur 001","Adur 002",..: 1468 125 126 127     128 129 130 131 132 133 ...          .. ..$ averagesp  : num [1:7194] 11.6 20.9 20.3 21.4 19.1 ...        .. ..$ nga        : int [1:7194] 41 112 177 97 120 NA 142 161 122 113 ...        .. ..$ mediansp   : num [1:7194] 11.5 13.1 12.6 13.5 10 ...          .. ..$ mbps2      : int [1:7194] 125 74 104 48 83 NA 64 67 76 25 ...         .. ..$ msoa.total : int [1:7194] 2372 114 191 97 126 NA 157 208 130 119 ...          .. ..$ p.nga      : num [1:7194] 1.73 98.25 92.67 100 95.24 ...          .. ..$ p.mbps2    : num [1:7194] 5.27 64.91 54.45 49.48 65.87 ...        .. ..$ residents  : Factor w/ 4057 levels "10,003","10,004",..: 2222 1748 27 1261       3216 NA 3380 403 3093 3256 ...        .. ..$ area       : Factor w/ 6908 levels "1,002.23","1,002.58",..: 3707 2869 2840      3304 839 NA 1899 1327 6811 748 ...          .. ..$ density    : num [1:7194] 25.5 31.3 46.9 24.8 72.1 NA 50.6 81.5 87.4 76.8 ...          .. ..$ micro.bc   : num [1:7194] 36 32 83 33 39 NA 74 60 44 32 ...           .. ..$ small.bc   : num [1:7194] 49 14 64 3 14 NA 50 14 14 73 ...        .. ..$ medium.bc  : num [1:7194] 850 0 15 0 0 NA 10 0 0 0 ...        .. ..$ large.bc   : num [1:7194] 245 0 0 0 0 NA 0 0 0 0 ...           .. ..$ primary    : num [1:7194] 60 0 0 0 0 NA 0 0 0 0 ...          .. ..$ secondary  : num [1:7194] 135 0 20 0 0 NA 0 0 5 5 ...         .. ..$ tertiary   : num [1:7194] 16600 95 335 100 145 NA 300 200 160 80 ...          .. ..$ kibs       : num [1:7194] 11805 20 80 20 30 ...           .. ..$ edu        : num [1:7194] 2205 2681 811 2658 116 ...          .. ..$ p.non.white: num [1:7194] 21.4 35 45.4 19 33.7 NA 32.8 30.4 38.8 30.8 ...        .. ..$ claimants  : num [1:7194] 117 245 316 254 295 ...         ..@ polygons   :List of 7194         **Can anyone suggest how I overcome this problem please?** Thanks, Ed


ArcGIS Publish query layer as Feature Service fails

I have a query layer that's simply a DB view call in PostgreSQL, the view is as such: ` CREATE OR REPLACE VIEW daily_query AS SELECT b.id, b.first_name, b.last_name, c.assignment_date::date AS assignment_date, d.case_id, d.case_sequence, d.status, a.address_complete, a.the_geom_point FROM universe a, staff b, enumerator_assignments c, enumerator_case_assignments d WHERE b.id = c.enumerator_id AND c.objectid = d.enumerator_assignments_id AND a.id = d.case_id;` The query layer is simply a call to the view `select * from daily_query` The layer renders fine in ArcMap - but then I try to publish the layer as a Service to ArcGIS Server, I want this to be a FeatureService so I enable "Feature Access" Capability and I only allow "Query" operation. and I don't allow "Geometry Updates". When I run "Analyze" it returns with "00134: Layer's data source is not supported" The Service Analyzes and Publishes fine if I don't enable "Feature Access" on it. I'm running 10.2.2 Any ideas?


How do I change the targeted .NET 4.0 Framework to version 3.5 in an existing VS 2010 EXPRESS project?

I am upgrading from ArcGIS desktop 10.0 to 10.1 and I have a bunch of addins developed in VS 2008 Express using ArcObjects 10.0. I would like to upgrade everything **to 10.1 using VS 2010 Express**. I understand that for ArcGIS Desktop **.NET Framework 3.5 Service Pack 1 is still required.** I converted my first existing addin to 10.1 in VS 2010 Express but it seems to be targeting .Net 4.0. I only found instructions on changing the target version in full version of VS but I can't find this feature in the Express edition. Can someone shed some light on this? EDIT: I am adding a screenshot from Project -> Project Properties -> Application to show where I would expect to change the Target Framework but I don't see where it could be changed. ![enter image description here](http://i.stack.imgur.com/C9nwd.jpg) **EDIT 2:** I found an option to change Target Framework (All Configurations) under **Project -> Project Properties -> Compile -> Advanced Compile Options -> Target Framework (All Configurations)** Can someone confirm that his is where I change the target framework in the 2010 Express edition?


MapInfo License Server - Can It Serve Multiple Versions of MapInfo

I have a MapInfo License Server Utility that serves 30+ MI v10.5 concurrent licenses across my organisation. v11.5 has now landed on my desk (with new serial numbers and access codes), so can I activate my v11.5 licenses within the same instance of License Server utility, whilst retaining the v10.5 licenses, in order to gradually migrate users onto v11.5 over the coming months?


How to calculate the shortest path between two points characterized by latitude and longitude with pgrouting?

I have to calculate the **length** of the shortest path between A (lat1, long1) and B (lat2, long2). But all the solutions I've found do not get a correct result because they only calculate the distance between the nearest vertex of A to the nearest vertex of B.


How can I select features by map clicks

I'm developing a flex web application using the ArcGIS flex API 3.1. I have a map with a polyline feature layer. I'd like to be able to select polylines from the map by clicking on the features on the map. Could you guide me to a sample? Regards.


Algorithm for finding irrregular polygon centroid (label point)

I need to find a centroid (or label point) for irregularly shaped polygons in Google Maps. I'm showing InfoWindows for parcels and need a place to anchor the the InfoWindow that is guaranteed to be on the surface. See images below. ![alt text](http://i.stack.imgur.com/sGcEX.png) ![alt text](http://i.stack.imgur.com/pQDX0.png) In reality I don't need anything Google Maps specific, just looking for an idea of how to automatically find this point. My first idea was to find the "false" centroid by taking the average lat and lngs and the randomly placing points out from there until I find one that intersects the polygon. I already have the point-in-polygon code. This just seems awfully "hacky" to me. I should note that I don't have access to any of the server side code outputting the geometry so I can't do anything like ST_PointOnSurface(the_geom).


How to display only one pie chart for polygon features in QGIS?

I'm working with Qgis on a project and I have a small problem. I have to do some demographic research and therefore I have a .shp file from a given country with 10 districts for example. Now, I know how to insert pie charts but the problem is that I only want to insert one pie chart for one district because this area is of importance to me for example. So I added a column in the attribute table with data for the district that I want to display and left all the other values 'NULL'. However, if I insert pie charts, they appear for all the districts, even though I've not inserted any values in the attribute table for the other 9 districts that are of no interest to me. Does anybody know how I can only display one pie chart for 1 feature of the polygon instead of having 1 good pie chart and 9 empty on the map?


QGIS Cross Tab of Towns and Routes

Using QGIS, I have a layer with 100 towns. I have another layer with 200 cycle routes, each with a different name. I want to create (automatically natch) a table that will help the people in each town, for example, get a list of all of the routes that pass within 1Km of their town centre.


Understanding web mapping tools

Before I start, this is going to be a general question but with specific cases. I've been trying to understand GIS tools and applications to achieve the following:   1. Build my own basemap   2. Publish the basemap online   3. Add my own data layers to the online basemap   4. Perform analysis (on-the-fly)    5. Add interactivity to the map So, In my first attempt, I built a spatially enabled database on PostGIS, queried this using PHP, overlayed it on google maps using their API and added interactivity using their api itself. I carried out analysis like area, denisty etc. on the fly using postGIS functions. For my next attempt, I downloaded TileMill and rendered some really beautiful maps and added zoom level functions to it. I also exported my MBtile files. When I came across trying to publish this online, I hit quite a few names, but got lost and confused:   1. TileStache   2. TileCache   3. TileStream   4. Mapserver   5. Mapnik   6. GeoServer   7. Mapbox So, I understand these are mapservers. But what exactly does a mapserver do? Why is it used? Why can't i server my tiles from a normal web Apache server? And for MBtiles, how many servers are out there which lets you publish it directly? I can also export my basemap as a png/vsvg format. Can these formats be used to make a basemap through a mapserver? When researching interactivity I came across:   1. Leaflet   2. ModestMaps   3. OpenLayers   4. Mapstraction I understand these are Javascript Mapping API's. But on what base map can I use these? Can it be used on Google maps? Can I use them on my own base maps too? I basically want to understand the differences between these tools and which would be the best software stack to build a fully customized, interactive map with my own basemap. Sorry if my question sounds to general. I will definitely be here to clarify my questions


Change Vector projection in Openlayers

I have added vector layer to map                   vactorTest = new OpenLayers.Layer.Vector("Editable Features", {     strategies: [new OpenLayers.Strategy.BBOX(), saveStrategy],     projection: new OpenLayers.Projection("EPSG:900913"),     protocol: new OpenLayers.Protocol.WFS({         srsName: "EPSG:900913",     url: "http://localhost:8080/geoserver/wfs",     featureType: "ROADS",     featureNS: "http://localhost:8080/MIA",     featurePrefix: "MIA",     geometryName: "GEOM",     version: "1.1.0"     })     });          map.addLayer(vactorTest);      And than added polygon draw control               var drawPolygon = new OpenLayers.Control.DrawFeature(     vactorTest, OpenLayers.Handler.Polygon,     {     title: "Draw Polygon",     displayClass: "olControlDrawFeaturePolygon",     multi: true     }     );      with drawn polygon I want to filter layer               testLayer.mergeNewParams({         "CQL_FILTER": mergeFilter     });      I use within filter: WITHIN(GEOM, POLYGON((4693721.49445 672633.89194, 4560244.10812 322255.75282, 4703732.29843 422363.79257, 4693721.49445 672633.89194))) target layer is in Georgia, and if I draw polygon near Italy than filter works.... How can I change vector projection like my layer to get right result?


Is there a way to figure out what files an MXD points to?

Let's say I load up two layers, and then save an .mxd file. Now I want to get the path names of my layers from the .mxd file. Is there a way to do this in ArcEngine? I'm using ArcEngine 10 with VS2010 & C#


Problem circle with zoom level

when I zoom level 16 and over I have two circle segments, I don't know if it's normal ( _projection_ ) or a bug. Have you an idea ? Openlayers : **2.13.1**   Problem With Chrome, IE and Firefox formula circle               function CircleMax(radius)     {     var pointsCircle = [];           var side = 160;      var angle;     var projectedLonLat, geometryPoint;     var points = [];         for (var i = 0; i < side; i++)     {         angle = (i * 360 / side) + side;         projectedLonLat = OpenLayers.Util.destinationVincenty(centerCircle, angle, radius);         projectedLonLat.transform(new OpenLayers.Projection("EPSG:4326"), new OpenLayers.Projection("EPSG:900913"));                     geometryPoint = new OpenLayers.Geometry.Point(projectedLonLat.lon, projectedLonLat.lat);         pointsCircle.push(geometryPoint);     }     var layerCircleStyle = OpenLayers.Util.extend({}, OpenLayers.Feature.Vector.style['default']);           var style_poly = OpenLayers.Util.extend({}, layerCircleStyle);       style_poly.strokeColor = "red";     style_poly.fillColor = "transparent";     var linearRing = new OpenLayers.Geometry.LinearRing(pointsCircle);     var polygonFeature = new OpenLayers.Feature.Vector(new OpenLayers.Geometry.Polygon([linearRing]),{},style_poly);     circleMaxLayer.addFeatures([polygonFeature]);      } Zoom level 6 ![Zoom le 6 ](http://i.stack.imgur.com/fKj0y.png) Zoom level 16 and over ![image description here](http://i.stack.imgur.com/2dVqA.png)


Using multiple basemaps with masking/clipping?

Is there a way to limit a basemap to within a boundary in arcmap 10? I am working on a project within San Joaquin county in California. I want the areas surrounding the county to be displayed with a light grey basemap with minimal details. Within the county itself, however, I want to use a more detailed basemap that has the information I am trying to portray. Is there a way to add a second basemap but have it be limited to the county boundaries? Thank you so much for any help you can give me.


how to solve Arcgis API for Java Script - Asp proxy error

I am using ArcGis 10.1 for ArcGIS Api for Javascript. I have faced the problem of {"code":400,"message":"Unable to complete operation.","details":[]}. while running the application in IIS and my Proxy.config file is look like this:               <serverUrl url="http://localhost:6080/arcgis/rest/services/" --------------------myservices                matchAll="true"></serverUrl>      How can i rectify this error..


Null values coerced to zero (0) when exporting to csv or txt

I extracted values from several rasters, and in some cases, the value is null. However, when I export the attribute table to a csv or a txt file (from ArcMap), the null values are converted to 0. Is there any way to prevent this from happening? These zeros are a real downer, as throw off my analysis, which is done in R.


convert XML to PostgreSQL using MATLAB

In a weeks time I should be able to transfer the GIS data provided in the XML format to PostgreSQL database. I have prepared the PostgreSQL database (empty) and looking at Matlab scripts. (I am on ubuntu platform). Matlab is our existing platform.The work is passed on to me by another person who used Mapping tools on MATLAB to export the data to viewers like ArcGIS/Google Earth.This was created for a prototype. (However, in the final product I will be exporting PostgreSQL to QGIS to view the data). I am pretty new to the GIS data/table.Could you suggest me the steps in-detail to have the data in PostgreSQL?


which framework is good for closed intranet

My requirement is simple. By latitude/longitude, draw icon on the map, draw a line between two icons.I may need put a customer picture over the map. It is not difficult by google map javascript or bing map javascript. The problem is that the system is running in a closed intranet environment. No internet access. Which framework shall I choose? Does it provide map? The system only focus on a specific area, such like a national park, it does not need the map for the whole world. Thanks!!!


How to apply Log Transform to Flow Accumulation?

I have run a flow accumulation for a DEM that I'm working with. I'm trying to apply a log transform because I've heard it could give a better visual representation. How would I go about doing this process.... Can this be done without Python?


AVG and GROUP BY in ArcMap

I have a Road Network shp file.(attributes include Names and Volumes ) Each road in my network consists of many links with the same name and different volumes. I want to display Road Names and volumes as labels using sql query or python function to group by names and have volume Average. ** In normal SQL query I would type:               SELECT ROAD_NM, AVG(VOL) as VOL     FROM TABLE     GROUP BY ROAD_NM;


How can I convert a DEM into areas of "flatness", with an arbitrary fill/cut to area ratio?

I want to find areas of relative "flatness" in a DEM. How can I find flat areas? I want to find a neighbourhood of contiguous cells that have an elevation variance that is within a certain tolerance threshold. I can use the ESRI suite of tools, or any open source GIS.


Find maximum and minimum value of points within a polygon buffer?

I have a series of 500m buffer zones along a polyline. I also have a series of points following the line, each point represents a 10m x 10m cell. I calculated the values of a field in each point known as "F18". I want to find the average, maximum and minimum values in F18 for all points within each buffer zone and add fields to the attribute table of the buffer for each. ArcGIS 10.1, arcinfo license.


Custom scalebar units in ArcMap

Is it possible to create a scalebar in ArcMap that shows distances in custom units? I understand how to do this manually, but it would be great if I could define a custom unit of measurement (let's say "1 swimming pool") and use that to dynamically generate a scalebar.


Django QGIS API, PyQt4 integrated

I am developing a web page through Django, where I want to insert images as static files. These images are produced with QGIS commands, which exist in QGIS API, PyQt4, and other API. So, my question is: Is it possible to incorporate the QGIOS libraries in Django? How do I do that? Thanks!


Eclipse export as runnable jar file error

I am building a destop application , which takes user input for personal geodatabase and finds feature class in it . It works fine till i run it in eclipse , but when i export it as a runnable jar file from eclipse Juno export option , it export it with warning and when i run the jar file from command promt i get following error ---               java.lang.IllegalArgumentException: URI is not hierarchical         at java.io.File.<init>(Unknown Source)         at com.esri.arcgis.interop.NativeLoader.getArcGISHome(Unknown Source)         at com.esri.arcgis.interop.NativeLoader.loadLibrary(Unknown Source)         at com.esri.arcgis.system.EngineInitializer.initializeEngine(Unknown Source)         at default.FeatureClassFinder.main(FeatureClassFinder.java:69)         at default.BrowsePath.actionPerformed(BrowsePath.java:139)         at javax.swing.AbstractButton.fireActionPerformed(Unknown Source)         at javax.swing.AbstractButton$Handler.actionPerformed(Unknown Source)         at javax.swing.DefaultButtonModel.fireActionPerformed(Unknown Source)         at javax.swing.DefaultButtonModel.setPressed(Unknown Source)         at javax.swing.plaf.basic.BasicButtonListener.mouseReleased(Unknown Source)         at java.awt.Component.processMouseEvent(Unknown Source)         at javax.swing.JComponent.processMouseEvent(Unknown Source)         at java.awt.Component.processEvent(Unknown Source)         at java.awt.Container.processEvent(Unknown Source)         at java.awt.Component.dispatchEventImpl(Unknown Source)         at java.awt.Container.dispatchEventImpl(Unknown Source)         at java.awt.Component.dispatchEvent(Unknown Source)         at java.awt.LightweightDispatcher.retargetMouseEvent(Unknown Source)         at java.awt.LightweightDispatcher.processMouseEvent(Unknown Source)         at java.awt.LightweightDispatcher.dispatchEvent(Unknown Source)         at java.awt.Container.dispatchEventImpl(Unknown Source)         at java.awt.Window.dispatchEventImpl(Unknown Source)         at java.awt.Component.dispatchEvent(Unknown Source)         at java.awt.EventQueue.dispatchEventImpl(Unknown Source)         at java.awt.EventQueue.access$200(Unknown Source)         at java.awt.EventQueue$3.run(Unknown Source)         at java.awt.EventQueue$3.run(Unknown Source)         at java.security.AccessController.doPrivileged(Native Method)         at java.security.ProtectionDomain$1.doIntersectionPrivilege(Unknown Source)         at java.security.ProtectionDomain$1.doIntersectionPrivilege(Unknown Source)         at java.awt.EventQueue$4.run(Unknown Source)         at java.awt.EventQueue$4.run(Unknown Source)         at java.security.AccessController.doPrivileged(Native Method)         at java.security.ProtectionDomain$1.doIntersectionPrivilege(Unknown Source)         at java.awt.EventQueue.dispatchEvent(Unknown Source)         at java.awt.EventDispatchThread.pumpOneEventForFilters(Unknown Source)         at java.awt.EventDispatchThread.pumpEventsForFilter(Unknown Source)         at java.awt.EventDispatchThread.pumpEventsForHierarchy(Unknown Source)         at java.awt.EventDispatchThread.pumpEvents(Unknown Source)         at java.awt.EventDispatchThread.pumpEvents(Unknown Source)         at java.awt.EventDispatchThread.run(Unknown Source)     Error: null     Sample failed.  Exiting...     java.lang.NullPointerException         at com.esri.arcgis.interop.NativeLoader.getArcGISHome(Unknown Source)         at com.esri.arcgis.interop.NativeLoader.loadLibrary(Unknown Source)         at com.esri.arcgis.system.EngineInitializer.initializeEngine(Unknown Source)         at default.FeatureClassFinder.main(FeatureClassFinder.java:69)         at default.BrowsePath.actionPerformed(BrowsePath.java:139)         at javax.swing.AbstractButton.fireActionPerformed(Unknown Source)         at javax.swing.AbstractButton$Handler.actionPerformed(Unknown Source)         at javax.swing.DefaultButtonModel.fireActionPerformed(Unknown Source)         at javax.swing.DefaultButtonModel.setPressed(Unknown Source)         at javax.swing.plaf.basic.BasicButtonListener.mouseReleased(Unknown Source)         at java.awt.Component.processMouseEvent(Unknown Source)         at javax.swing.JComponent.processMouseEvent(Unknown Source)         at java.awt.Component.processEvent(Unknown Source)         at java.awt.Container.processEvent(Unknown Source)         at java.awt.Component.dispatchEventImpl(Unknown Source)         at java.awt.Container.dispatchEventImpl(Unknown Source)         at java.awt.Component.dispatchEvent(Unknown Source)         at java.awt.LightweightDispatcher.retargetMouseEvent(Unknown Source)         at java.awt.LightweightDispatcher.processMouseEvent(Unknown Source)         at java.awt.LightweightDispatcher.dispatchEvent(Unknown Source)         at java.awt.Container.dispatchEventImpl(Unknown Source)         at java.awt.Window.dispatchEventImpl(Unknown Source)         at java.awt.Component.dispatchEvent(Unknown Source)         at java.awt.EventQueue.dispatchEventImpl(Unknown Source)         at java.awt.EventQueue.access$200(Unknown Source)         at java.awt.EventQueue$3.run(Unknown Source)         at java.awt.EventQueue$3.run(Unknown Source)         at java.security.AccessController.doPrivileged(Native Method)         at java.security.ProtectionDomain$1.doIntersectionPrivilege(Unknown Source)         at java.security.ProtectionDomain$1.doIntersectionPrivilege(Unknown Source)         at java.awt.EventQueue$4.run(Unknown Source)         at java.awt.EventQueue$4.run(Unknown Source)         at java.security.AccessController.doPrivileged(Native Method)         at java.security.ProtectionDomain$1.doIntersectionPrivilege(Unknown Source)         at java.awt.EventQueue.dispatchEvent(Unknown Source)         at java.awt.EventDispatchThread.pumpOneEventForFilters(Unknown Source)         at java.awt.EventDispatchThread.pumpEventsForFilter(Unknown Source)         at java.awt.EventDispatchThread.pumpEventsForHierarchy(Unknown Source)         at java.awt.EventDispatchThread.pumpEvents(Unknown Source)         at java.awt.EventDispatchThread.pumpEvents(Unknown Source)         at java.awt.EventDispatchThread.run(Unknown Source)      I don't have any idea whats wrong neing done. Please provide me the coorect way to export my application. Thanks in advance////


examples for Text3dElement

i am working in arcgis10. in java. can you all provide me with links to some sample program of `Text3dElement`. til now what i have implemented is given below                              Text3DElement td=new Text3DElement();                       td.setText("run away");                    IElement e=new Text3DElement();                    e=(IElement)td;                    e.setGeometry(point);                          iGlbCnt=new GlobeGraphicsLayer();                   iGlbCnt.addElement(e,0);      ` i was originally using `TextElement`, but it was not working properly on 3d globe. hence, i am trying to work with `Text3dElement` but unable to find any sort of example. thnks for ur help.


How to install PostGIS and pgrouting on Windows if Stack Builder fails to connect?

I know this has probably been covered and that there are plenty of guides out there, but for the life of me I cannot get pgrouting installed. All guides/how-tos are different or assume prior knowledge. I'm a QGIS newbie and please assume I know nothing! **How do install pgrouting? And what do I need first and how do I install whatever is needed? Can some helpful soul offer a step-by-step guide?** I've been trying to do this for days. I can get Postgres installed but the Stack Builder will not connect. I am thoroughly confused. I'm using Windows 7 Enterprise 64bit Service Pack 1, QGIS Lisboa. I'm now having trouble getting PostgreSQL to install.               Loading additional SQL modules...      Executing cscript //NoLogo "C:\Program Files (x86)\PostgreSQL\9.2\installer\server\loadmodules.vbs" "postgres" " _*_ *" "C:\Program Files (x86)\PostgreSQL\9.2" "C:\Program Files (x86)\PostgreSQL\9.2\data" 5432 Script exit code: 2 Script output: Installing the adminpack module in the postgres database... Executing 'C:\Users\ucesga1\AppData\Local\Temp\rad87BB2.bat'... psql: could not connect to server: Connection refused (0x0000274D/10061) Is the server running on host "localhost" (::1) and accepting TCP/IP connections on port 5432? could not connect to server: Connection refused (0x0000274D/10061) Is the server running on host "localhost" (127.0.0.1) and accepting TCP/IP connections on port 5432? Failed to install the 'adminpack' module in the 'postgres' database loadmodules.vbs ran to completion Script stderr: Program ended with an error exit code Error running cscript //NoLogo "C:\Program Files (x86)\PostgreSQL\9.2\installer\server\loadmodules.vbs" "postgres" " _*_ *" "C:\Program Files (x86)\PostgreSQL\9.2" "C:\Program Files (x86)\PostgreSQL\9.2\data" 5432 : Program ended with an error exit code ![enter image description here](http://i.stack.imgur.com/hCSC3.jpg)


WMS: point "pointer" symbol with offset - getFeatureInfo search area slight mismatch

I'd like to ask what ways there might be to adjust a point's getFeatureInfo search area so that it matches a "pointer" (i.e slightly offset) symbol. We are running GeoServer 2.2.0 (part of an Opengeo Suite) and the symbology(SLD) for some point layers is various forms of "pointers", that is to say that the center of the icon is not at the center of the point, rather it points to the point, like an arrow or a badge with a small pointer at the bottom or a map pin. (SVG graphics.) There appears to be two ways one can get the point symbols offset: I. Modify the SVG graphic to contain a transparent area in one end, see exhibit A. This method provides an empty area for getFeatureInfo corresponding to the transparent area in the SVG symbol as in geoserver's docs on point slds - one can however then click to the side of the visible symbol and still get a hit, which is unwanted. II. Use the offset function in the SLD for the point layer. The effect this has is that while the "pointer", or pin symbol is properly offset the search area remains centered on the point. See B. Many thanks for any suggestions.


Split all OSM roads within boundingbox in 20m segments and save to new table

I have a PostGIS Database with all streets from osm within a given boundingbox. Now I want to split all roads from my table planet_osm_roads into segments of max 20 meters length and save those segments to a new postgis table. The geometry of the roads is save in "way geometry(linestring,900913)". Can anyone give me a hint on how to accomplish this? Cheers, Martin


Arcpad Message Box and Keyboard Button Issue

I'm using ArcPad 8 along with ArcPad Studio to create a custom application for use on Windows Mobile 6 devices. There are data entry forms and on those forms are buttons that, when clicked, present the user some information in a message box. For a button's onclick event, I specify an event script (using JScript) such as:               Application.MessageBox("Some information for the user...");      That all works fine. However, when the user dismisses the message box, the button for toggling the input keyboard/transcriber/recognizer that is located in the lower right corner of the application window immediately disappears, which prevents the user from showing and hiding the keyboard as they attempt to enter data in other various controls on the form. Has anyone else ever had this problem and know of any solutions? Thanks!


Topological relationship between two polygon feature layers

I have two polygon feature layers: one representing land parcels and one representing roads. Some land parcels (polygons) coincide with a part of a road (another polygon feature). How can I check which land parcels coincide in any of their boundaries with a road? Of course help will be great if a code in VBA and ArcObjects is already available. I note also that I have considered the topology rules defined in http://webhelp.esri.com/arcgisdesktop/9.3/index.cfm?TopicName=Topology_rules and also and the ITopologyRule interface but I could not be sure which rule in both cases to use for solving my problem. Thanks. Demetris


Any freely available vector data in cylindrical projection?

I'm trying to develop an application in PHP to plot various points in Europe from a database, and serve them dynamically on a map. I think this should be relatively straightforward when using a map in the cylindrical projection - but I just cannot find one. I only require country borders, and pretty low level of detail. The smallest section of the map to be displayed at one time is likely to be the UK, and the largest output size is probably around 1280 x 1024px. Does anyone know of a freely available (or paid, if free isn't an option) vector? Many thanks, Jack


Arcmap: save all layers in Table of Contents to layer files

Using python, how can one save all layers in the Arcmap (v10) table of contents pane to layer files on disk?


Clearing map cache - ArcGIS Engine

Is there a programmatic way to clear the cache of a map in ArcGIS Engine 10.0? It is apparent that there is some caching scheme operating under the covers, and I'd like to clear all of its cached data. To be clear it isn't _leaking_ memory, I just don't want it to be caching anything. This is a fairly vanilla (*.mxd). Just a raster dataset, no basemap layers or anything fancy going on. ![graph of memory](http://i.stack.imgur.com/CFpby.png)


QGIS Mosaicing with null values in images

I have a set of landuse landcover data which needs to be mosaiced. Whenever i try using QGIS Raster>miscellenous>merge tool. giving no datavalue as '0', it returns either **(A)**                ' Traceback (most recent call last):      File "C:\PROGRA~2\QGISVA~1\bin\gdal_merge.py", line 509, in      sys.exit(main())      File "C:\PROGRA~2\QGISVA~1\bin\gdal_merge.py", line 495, in main      fi.copy_into( t_fh, band, band, nodata )      File "C:\PROGRA~2\QGISVA~1\bin\gdal_merge.py", line 241, in copy_into      nodata_arg )      File "C:\PROGRA~2\QGISVA~1\bin\gdal_merge.py", line 61, in raster_copy      nodata )      File "C:\PROGRA~2\QGISVA~1\bin\gdal_merge.py", line 101, in raster_copy_with_nodata      to_write = Numeric.choose( nodata_test, (data_src, data_dst) )      File "C:\PROGRA~2\QGISVA~1\apps\Python27\lib\site-packages\numpy\core\fromnumeric.py", line 296, in choose      return choose(choices, out=out, mode=mode)      MemoryError       or **(B)** ' The process crashed sometime after starting successfully' I am using QGIS 2.2 on Windows 8 PC. I have a set of 30 images in .tiff format which needs to be mosaiced. I also tried doing it with 2-3 images but the same problem persists.


Saving in arc ascii grid format

This might be a really simple question but I can't figure it out! I need to save a DEM in Arc ascii grid format so that I can load it in to Matlab. If I explain a bit about what I am doing maybe it will help: I have a 3D plane produced from the CFtool in matlab and I need to work out the volume between this plane and my DEM. I was going to import into matlab and try working it out. I have tried the other way around, i.e. exporting my plane from matlab but this doesn't seem to work. Anybody know how to save the DEM (in tif at the mo) as arc ascii grid format? Or have any other ideas on what to to work out my volume. Working in Arcmap10 and Matlab R2012a


Trace Downstream task: which ArcObject interface and method?

Does anyone know what ArcObjects interface and method the Trace Downstream trace task on Arcmap's Utility Network Analyst toolbar calls?


OpenLayers plugin fails to install

When I try to install OpenLayers plugin for QGIS 2.2.0 Valmiera I get the error message the plugin is broken; python said: 'QPyNULLVARIANT'object has no attribute '_ _getitem_ ' Does anyone have any ideas; i am a complete beginner, so I hope for a simple fix. Or, where else can I get base maps. I am using mac OS X 10.9.2 Thanks.


Problems installing r.roughness

I am trying to install r.roughness in GRASS 6.4.3. It starts compiling but at the end I get this error:               ERROR: Module did not install properly      I cannot identify the problem.   Can somebody that has come across the same error help me?


Identify corrupt shapefiles or rasters via Python or ArcObjects

I'm interested in a way to identify rasters (e.g., ESRI GRID format) or shapefiles that are corrupt that would prevent them from being used in an analysis tool. I have a large number in structured folder tree and would like to be able to script something that will walk the file tree and be able to ID rasters that, for example if you tried to view it in ArcCatalog it would not allow it to be displayed and accordingly would make any geoprocessing tool fail. Same deal for shapefiles or feature classes. My preference would be somehow to do this in Python (and I wouldn't be opposed to a non-ESRI python module or other) but I'm open for anything. I guess you could hack something together by trying to validate geometry or capitalizing on some aspect of the file design, but what would make the most sense? Or is there another way to do it that is specifically designed for such a purpose? My goal is to just be able to run a simple Quality Control mechanism on some large week-long preprocessing files to make sure they're all sound before I start doing the actual batch processing. thanks, Tom


Changing the line color of a marker box using OpenLayers

I'm trying to figure out how to change the line color of a marker box in openlayers. I've been learning from the following example: ![enter image description here](http://i.stack.imgur.com/KYWHt.png)               <!DOCTYPE html>                            <meta http-equiv="Content-Type" content="text/html; charset=utf-8">             <meta name="viewport" content="width=device-width, initial-scale=1.0, maximum-scale=1.0, user-scalable=0">             <meta name="apple-mobile-web-app-capable" content="yes">             <title>OpenLayers Boxes Vector Example</title>             <link rel="stylesheet" href="http://openlayers.org/dev/theme/default/style.css" type="text/css">             <link rel="stylesheet" href="http://openlayers.org/dev/examples/style.css" type="text/css">
Boxes Example Vector

                 box, vector, annotation

                 Demonstrate marker and box type annotations on a map.

                         What I tried doing was replacing the `boxes` line with the following:                           //var boxes  = new OpenLayers.Layer.Vector( "Boxes" );                 var styleMapBoxes = new OpenLayers.StyleMap({borderColor: '#008000'});             var boxes  = new OpenLayers.Layer.Boxes( 'Boxes', {styleMap: styleMapBoxes});        but it did not work. I've read through the example here, but still do not understand how I can edit the box style. According to the documentation `borderColor` is a valid attribute.


Algorithm to match trajectory

I've data in the format (lattitude,longitude,timestamp). What I want to do is between timestamp1 and timestamp2 I want to find out which people might have taken the same routes. Is there some kind of trajectory matching algorithm? Let's say we have 100 users each has a set of points of (lattitude,longitude,timestamp) I want to know which set of users might have matching trajectory. Matching trajectory would be the same route taken, hence in a given set of timestamps the lattitude and longitude should be same or close enough as well as the timestamp should be same or close enough. Close enough can be about 30 seconds for timestamp while for space let it be like 200 metres. I can do this via a brute force approach and I'm looking for better solutions. Posted it on SO: http://stackoverflow.com/questions/16725427/design-an- algortihm-to-match-trajectory I hope I would get more thoughts here.


Is there an industry-standard official mapping of GLONASS satellites to "PRN" values?

In the U.S. GPS system, satellites are uniquely identified by the PRN number broadcast by each satellite. This ID is usually visible on most GPS receiver platforms that provide access to individual GPS satellite signal information to identify the broadcasting satellite (e.g., accessible on Android via the GPSStatus.Listener and the GpsSatellite.getPrn() method (http://goo.gl/QLse3)). However, the Russian GLONASS system does not currently have the concept of a PRN number broadcast by every satellite, as many of the legacy satellites are an FDMA system, not a CDMA system such as the U.S. So, in GLONASS these satellites broadcast the same code on a different frequency. However, an ID is still available for GLONASS satellites for many GPS receiver positioning platforms that support dual GPS/GLONASS positioning, often via the same access point as the PRN value for U.S. GPS satellites (e.g., on Android via the GpsSatellite.getPrn() method, with seemlingly standardized values across devices I've seen from 65 to 88. See the first screenshot listed in the GPSTest app on Google Play (https://play.google.com/store/apps/details?id=com.android.gpstest) for an example list of U.S. and GLONASS satellites obtained from the GPSStatus.Listener on a Samsung Galaxy S3). I haven't yet been able to find an official industry-standard source for how these ID values for GLONASS satellites are determined. However, proprietary documentation for particular products have shown GLONASS satellites numbered from 65 to 88 [1] and 38 to 61 [2]. Is there any industry-standard official mapping between GLONASS satellites and "PRN" values?


Unable to reference GPS waypoints correctly in ArcGIS

I'm working on a project. Waypoints were taken in the field with a Garmin Dakota 20. I've tried two different ways of putting them into ArcGIS. One way is through converting the GPX file to a shapefile via DNRGPS. This results in the coordinate system being wrong (even if changing it in the program) and the waypoints end up at the wrong place. The other is with the data interoperability tool which results in the waypoinst missing spatial reference. Even after adding such in through the catalogue this doesnt move the Waypoints anywhere. What am I doing wrong? I want them matched to the current data frame with its maps (Sweref99) **EDIT:** After editing once again through the ArcCatalog i tried changing it to Geographic Coordinate System instead of the projected (still Sweref99) and now, instead of them being nowhere, I get them at the same spot I did through dNRGPS (quite a bit off, but at least on the map).


ArcObjects Background Process (Thread) Progress Bar?

Using a toolbar, written in ArcObjects, it is possible to push processing to a separate thread so that the GUI does not lock. Is it possible to have the thread update a progress bar in the same position (lower right - arcmap) that a geoprocessing tool does? I must be missing some obvious documentation to gain access to that progress bar update location.


Why does install of add-on r.clump2 for GRASS 7 give error?

I am trying to install the add-on r.clump2 for grass7. What is wrong and how do I fix it?. Thanks The error code im receiving is the following:               > (Fri Jul  5 15:13:42 2013)                                                 > g.extension -s extension=r.clump2     > svnurl=http://svn.osgeo.org/grass/grass-addons/grass7 Fetching     > <r.clump2> from GRASS-Addons SVN (be patient)... Compiling...     > flag.c:1:23: fatal error: grass/gis.h: No such file or directory      > #include <grass/gis.h>     >                        ^ compilation terminated. make: *** [OBJ.x86_64-apple-darwin10.8.0/flag.o] Error 1


Parallel Computing and Raster Operations in ArcObjects gives cryptic error

This is a problem I have been struggling with for awhile. The code below takes 1-3 input feature layers (pFLayer) that were stored into a list called pFeatureLayerCollection, casts them into IGeoDatasets and executes EucDistanceFull in parallel to get distance rasters from them. The output IGeoDatasets are added to a list. After the Parallel.ForEach ends, if there is more than 1 output raster, I attempt to get the localminimum of the rasters using the method below. It runs fine if I give it one input raster. It also runs fine if I run it with multiple rasters sequentially (not in parallel) However, if I try to execute it with more than one raster in parallel, it fails with this error every time: **Attempted to read or write protected memory. This is often an indication that other memory is corrupt. (Esri.ArcGIS.SpatialAnalyst at ESRI.ArcGIS.SpatialAnalyst.RasterLocalOpClass.LocalStatistics(IGeoDataset LocalData, esriGeoAnalysisStatisticsEnum Type)**                            ESRI.ArcGIS.Geometry.IEnvelope pExt; //create a new envelope to hold the feature class' size information                 pExt = new EnvelopeClass(); //instantiate the envelope as an object                 pExt.XMin = **[Some given known double value]**;                 pExt.YMin = **[Some given known double value]**;                 pExt.XMax = **[Some given known double value]**;                 pExt.YMax = **[Some given known double value]**;                      string StrOutRasPath = **[Some given valid filepath]**;                      List<IGeoDataset> pMinimumGeodatasetCollection = new List<IGeoDataset>(); //store output distance rasters                           Parallel.ForEach(pFeatureLayerCollection, pFLayer =>                 {                     ESRI.ArcGIS.DataSourcesRaster.IRasterWorkspace pTempRWS;                     ESRI.ArcGIS.Geodatabase.IGeoDataset pOutputRaster = null;                     ESRI.ArcGIS.GeoAnalyst.IRasterAnalysisEnvironment pEnv = null; //Create a raster analysis environment. This sets up an area to do raster operations in                     ESRI.ArcGIS.SpatialAnalyst.IDistanceOp pDistanceOP; //Create raster DistanceOp object                     pDistanceOP = new RasterDistanceOpClass(); //define the newly created object as a new raster distance operation object                          string g_tempFolder = null;                     lock (strOutRasPath) //threadsafe creation of folders from method                     {                         g_tempFolder = CreateTempFolder(strOutRasPath);                     }                               lock (pWSF)                     {                         pTempRWS = (IRasterWorkspace)pWSF.OpenFromFile(g_tempFolder, 0);                     }                                         ESRI.ArcGIS.Geodatabase.IWorkspaceFactory pWSFact; //create general workspace factory (this is the directory you will be working in)                     //pWSFact = new ShapefileWorkspaceFactoryClass(); //set general workspace as a shapefile workspace. WARNING: USING THIS CAUSES A STRANGE COM REWRAPPING ERROR. IT IS AN ESRI BUG                          Type t = Type.GetTypeFromProgID("esriDataSourcesFile.ShapefileWorkspaceFactory"); //DO NOT REMOVE THIS LINE OR THE FOLLOWING TWO LINES. Fixes strange casting error from COM rewrapping. (ESRI BUGFIX)                     System.Object obj2 = Activator.CreateInstance(t);                          pWSFact = (ShapefileWorkspaceFactory)obj2;                          ESRI.ArcGIS.Geodatabase.IFeatureWorkspace pFeatureWorkspace; //create a feature workspace to handle features                     pFeatureWorkspace = (IFeatureWorkspace)pWSFact.OpenFromFile(g_tempFolder, 0); //set feature workspace as the shapefile workspace defined earlier, and specify a file location                                    ESRI.ArcGIS.Geodatabase.IGeoDataset pTempMinRaster;                                  ESRI.ArcGIS.Geodatabase.IGeoDataset pTempGeodataset;                     ESRI.ArcGIS.Carto.IFeatureLayer pTempFeatureLayer;                          ////////////                     //Loop through all features.                     ////////////                          ESRI.ArcGIS.GeoAnalyst.IConversionOp pConvOp;                     pConvOp = new RasterConversionOpClass();                          ESRI.ArcGIS.Carto.IFeatureSelection pFeatureSelection;                     pFeatureSelection = (IFeatureSelection)pFLayer; // ***IMPORTANT*** THIS SETS THE CURRENT FEATURE LAYER IN COMBOBOX TO THE FEATURE SELECTION TO BE USED FOR CONVERSION TO RASTER                      ESRI.ArcGIS.Geodatabase.ISelectionSet pSelectionSet;                     pSelectionSet = pFeatureSelection.SelectionSet; //Defines a selection set of features. The pFeatureLayer information is stored here.                          ESRI.ArcGIS.Geodatabase.IFeatureCursor pFCursor2;                     ESRI.ArcGIS.Geodatabase.IFeature pFeature2;                     ESRI.ArcGIS.Geodatabase.IQueryFilter pQueryFilter;                          pQueryFilter = new QueryFilterClass();                     pFCursor2 = pFLayer.Search(null, false);                     pSelectionSet = pFeatureSelection.SelectionSet;                     pFeature2 = pFCursor2.NextFeature();                          pEnv = (IRasterAnalysisEnvironment)pDistanceOP;                     pEnv.OutWorkspace = (IWorkspace)pTempRWS;                     pEnv.SetCellSize(esriRasterEnvSettingEnum.esriRasterEnvValue, cellsize);                          lock (pExt) //threadsafe grab of pExt from outside Parallel.Foreach                     {                         pEnv.SetExtent(esriRasterEnvSettingEnum.esriRasterEnvValue, pExt);                     }                                        ESRI.ArcGIS.Geodatabase.IGeoDataset temporaryGeodataset = (IGeoDataset)pFLayer.FeatureClass;                              pOutputRaster = pDistanceOP.EucDistanceFull(temporaryGeodataset, true, false, false);                                    lock (pMinimumGeodatasetCollection) //threadsafe adding to the list of output geodatasets                     {                         pMinimumGeodatasetCollection.Add(pOutputRaster);                     }                          pFeatureSelection.Clear();                      });                      //LIST CREATION COMPLETE. NOW IF TWO OR MORE RASTERS WERE CREATED, GET THE MINIMUM USING LOCALMINIMUM METHOD                      if (pMinimumGeodatasetCollection.Count == 1) //if just one input feature class                 {                     pMinimumRaster = pMinimumGeodatasetCollection[0];                     pFeatureLayer = pFeatureLayerCollection[0];                      }                      else if (pMinimumGeodatasetCollection.Count == 2) //if two input feature classes, get local minimum                 {                          pMinimumRaster = LocalMinimum2Rasters(pMinimumGeodatasetCollection[0], pMinimumGeodatasetCollection[1], strOutRasPath + "\\tmpVoro0004");                      }      There is the bulk of the code, and below is the implemented LocalMinimum method. LocalMinimum3Rasters does the exact same as LocalMinimum2Rasters, but acommodates for one more raster. The final line of this code is where the program is failing.                   public IGeoDataset LocalMinimum3Rasters(IGeoDataset pRaster1, IGeoDataset pRaster2, IGeoDataset pRaster3, string strOutPath)         {                  ESRI.ArcGIS.DataSourcesRaster.IRasterBandCollection pRBCollTmp;             ESRI.ArcGIS.DataSourcesRaster.IRasterBand pRBand;                  ESRI.ArcGIS.DataSourcesRaster.IRasterBandCollection pRBColl;             pRBColl = new RasterClass();             pRBCollTmp = (IRasterBandCollection)pRaster1;             pRBand = pRBCollTmp.Item(0);             pRBColl.AppendBand(pRBand);             pRBCollTmp = (IRasterBandCollection)pRaster2;             pRBand = pRBCollTmp.Item(0);             pRBColl.AppendBand(pRBand);             pRBCollTmp = (IRasterBandCollection)pRaster3;             pRBand = pRBCollTmp.Item(0);             pRBColl.AppendBand(pRBand);                  //Create a RasterLocalOp operator             ESRI.ArcGIS.SpatialAnalyst.ILocalOp pLocalOp = new RasterLocalOpClass();                  ESRI.ArcGIS.Geodatabase.IWorkspace pWs;             ESRI.ArcGIS.Geodatabase.IWorkspaceFactory pWSF = new RasterWorkspaceFactoryClass();             pWs = pWSF.OpenFromFile(strOutPath, 0);                  return pLocalOp.LocalStatistics((IGeoDataset)pRBColl, esriGeoAnalysisStatisticsEnum.esriGeoAnalysisStatsMinimum);              }      What could be causing this error? It seems to be directly tied to parallelism - I've noted that all operations work perfectly fine if I use a regular foreach instead of Parallel.ForEach. Extra information: The program creates the two distance rasters every time - If I go to the temporary folders I created in this code after the program crashes, I will find the two temporary rasters wholly intact and perfectly created.


How to fill selected feature(s) with transparent red fill using ArcGIS Server 10?

I'm writing a Server Object Extension (SOE) for ArcGIS Server 10, and I'd like to change the symbology of features selected from my feature layer such that the features are shaded with a transparent red fill. Here's the code I have so far. This code results in selected features being outlined in red, and filled with a solid red fill, where it is supposed to be a transparent red fill.                       featureClassToQuery = (IFeatureClass)mapServerDataAccess.GetDisplayDataSource(mapName, layerIndex);                  parcelSelectionQueryFilter = new QueryFilter()             {                 WhereClause = string.Format("APN='{0}'", apn)             };                  featureSelection = (IFeatureSelection)featureLayer;             featureSelection.SelectionSymbol = (ESRI.ArcGIS.Display.ISymbol)new ESRI.ArcGIS.Display.SimpleFillSymbol()             {                 Color = new ESRI.ArcGIS.Display.RgbColor() { Blue = 0, Green = 0, Red = 255, Transparency = 50 },                 Outline = (ESRI.ArcGIS.Display.ILineSymbol)new ESRI.ArcGIS.Display.SimpleLineSymbol()                 {                     Color = new ESRI.ArcGIS.Display.RgbColor() { Blue = 0, Green = 0, Red = 0 },                     Style = ESRI.ArcGIS.Display.esriSimpleLineStyle.esriSLSSolid,                     Width = 2.00                 },                 Style = ESRI.ArcGIS.Display.esriSimpleFillStyle.esriSFSSolid             };             featureSelection.SetSelectionSymbol = true;             featureSelection.SelectFeatures(parcelSelectionQueryFilter, esriSelectionResultEnum.esriSelectionResultNew, false);