Help getting my flex service publicly accessible

My organization is trying to get our web maps up and running. I have created a simple web map service on our in house ArcGIS server, using the Arcgis flexviewer that is visible locally, works great, and is visible on arcgis.com. We have a 2 server set up where Server A is our internal server that everyone in the company stores their data on. this is the machine we have our install of ArcGIS 10.1 for server on. It is not public, and for internal staff use only.Server B on the other hand is off site and set up by a company that hosts our website to specifically host our web maps. It is in a DMZ and set up to get through our company firewall. It is just a simple web server running IIS 7 and has the ArcGIS Web Adaptor installed .They, like me, have no experience in setting something like this up so I'm learning as much as I can along the way.. We have our flex map up and running [here] which is great. What is not so great is the error that is popping up when the flex map is opened. It looks like my data, which consists of a a few poly lines, some points and some labels, are not making it over to Server B from Server A. I guess my real question is what do I need to do to get the data from Server A show up in my flex app on server B. If anyone can shed some light on this that would be great.


Openlayer plugin Google Hybrid date

Comment connaître la date de capture d'une image Google Hybrid sur QGIS? How to know when capturing an image on Google Hybrid QGIS? **French Speakers - Please help fix the above machine translation**


NDS file format conversion to an executable GIS format

Could we open .nds files in any of the GIS mapping environment. it stands for Navigation Data Standards, having the navigation features stored in it. There are 2 types of nds file formats. One is executable one (in sqlite interface) and the other non executable (because of its binary format). I want to customize these nds files to view the data stored in it. I require a software technology to open these files and view the stored data on any of the mapping software which is compatible.


How to build a query tool for spatialite data in QGIS?

I am trying to build some tool for a collegue of mine to query data from a spatialite database while clicking on vector features. My first guess is, that I could do this with an "action" I define for the layer. The example actions provided for layers unfortunately do not give me any hint how to perform, and more detailed documentation on actions is quite absent. I need a query, cause each feature in layer has many corresponding entries in the database table (one-to-many). While performing the action, I would like to pop up simple table which lists the corresponding data sharing the same id as the feature. In the example action, i can query an attribute of the layer with this:               echo "[% "id" %]"      which then looks like this: ![enter image description here](http://i.stack.imgur.com/MxJ3W.jpg) Ok, that's a start. But what i need is a table which lists all items in a spatialite table with the same id-attribute as the feature. My data structure:               Layer attribute:      id          Data table:      id | name | year      The layer itself is stored in the spatialite database too, if thats of any help. Can anyone give me a hint on how to proceed or direct me to some other web ressources on this topic? Thanx Bernd


Automaticallly select and replace all occurences of a font face

I have around 50 mxd file that is done some times earlier. The text font used in labeling was not licensed. Now I can Change the text in layout element. But it has also some text in map view too. But arcpy.mapping.ListLayoutElements can not find the text document. Is there any solution in arcpy?


OpenLayers: how to display the character degrees in the label

I would like to display temperature in label. So how can I write DEGREE SIGN? to be display as ° I was trying write: & # 176; or & deg; or \ & deg; and nothing. Can anybody knows? I greet.


Field Type Date invalid when run in Python, but not in Model

Using ArcGIS Desktop Advanced 10.1 and PyScripter When running the arcpy.FeatureClassToFeatureClass_conversion, I get the error that the field type of Date is not valid. However, when running this same procedure in the model where the python was exported from, it runs successfully. I have verified that the field type is in fact Date for the field and that the field length is sufficient.   This seeems to be an issue after moving to 10.1 Error Message:               ExecuteError: ERROR 000278: Field Mapping error: Field type invalid for output field xxx and input field xxx      Any ideas?


FME: Can I add multiple feature types to the SQL Server Spatial reader?

The title's it, really. With .TAB file reader, you can add multiple layers to the reader and use the 'merge feature type' check box with an '*' to read all the layers listed in the source files, which can then be separated out again on featuretype name. With the SQL Server Spatial reader, you can specify many layers in the database selection, but it seems that the 'merge feature type' option doesn't work. Is this a bug, or expected behavior, or am I missing something? EDIT: I have found that I can do this sometimes... so the above may be incorrect. I'll need to test some more.


How to Create a Vector Field from Separate u and v Data Files?

I am using Quantum GIS and trying to create a vector field using data from two separate netcdfs. One contains the u data, one contains the v data. Can this be done, and if so, how? Thanks. Jason


Get the junction id for the roads in a network dataset

I have a network dataset with roads and junctions. I am wondering what class I can use in arcpy to extract FromJunctionId and the ToJunctionId for each road segment in the network dataset. Thanks.


Lost geometry when joining spatial and non-spatial tables

I have two tables I am joining: 1 spatial (Table 1) and one not spatial (Table 2). Table 1 looks like this: | ID | geometry | Table 2: | ID | ... | When I perform the query: > CREATE TABLE AS SELECT * FROM t1, t2 WHERE t1.ID = t2.ID; I get a table with the geometry column included. But unfortunately it is no longer recognized as a spatial table (spatialite_gui at least). Also, this mean I can't do what I really want to do, which is dump the table as shapefile. So my question is: Why is the new table not recognised as a spatial table with a geometry column?


How to write labels with varying significant figures/decimal places?

I have results that range from ex. 113,009 down to 0.00495 - Can someone recommend a way to write/construct a label expression for the following in ArcMap: • If less than 1, use 2 decimal places (or three if absolutely necessary to show the first non-zero value) • From 1 to 10, use 1 decimal place • If 10 or greater, use no decimal places Or a way to quickly generate a new text column/field which adheres to the above?


Why does DataGridView data not update for ArcMap add-in?

I have created a dockwindow（name:FrmWellConfig) with a datagridview(datagridview1) control in a new ArcMap add-in project. I then bind the data:               datagridview1.datasouce=wellconfig.tables[wellcolunms];      Why does the datagridview not update? Code behind tool：               FrmWellConfig frm = new FrmWellConfig(SettingWorkspace.AppBasePath, wellConfigFileName, frmSave.MapType, SettingWorkspace.workspace1.CnnStr_InfoDB);      Code behind FrmWellconfig(wellConfig1 is a dataset)：                /// <summary>     /// Designer class of the dockable window add-in. It contains user interfaces that     /// make up the dockable window.     /// </summary>     public partial class FrmWellConfig : UserControl     {         private string cnnStr_OLEDB;         private bool nFlag = false;           public static string appBasePath;         private string wellConfigFileName = "";         sylMapType mapType;                   public sylMapType MapType         {             get { return mapType; }         }         public WellConfig WellConfigData         {             get { return wellConfig1; }         }              public FrmWellConfig(object hook)         {             InitializeComponent();             this.Hook = hook;             this.dataGridView1.DataError += delegate(object sender, DataGridViewDataErrorEventArgs e) { };         }                   public FrmWellConfig(string AppBasePath, string WellConfigFileName, sylMapType _mapType, string CnnStr_OLEDB)         {                  InitializeComponent();             nFlag = true;                    frmType = sylChildFormType.DialogForm;             cnnStr_OLEDB = CnnStr_OLEDB;             appBasePath = AppBasePath;                mapType = _mapType;             dataGridView1.Columns[0].Visible = false;             wellConfigFileName = WellConfigFileName;             wellConfig1.ReadXml(wellConfigFileName);         }              /// <summary>         /// Host object of the dockable window         /// </summary>         private object Hook         {             get;             set;         }              /// <summary>         /// Implementation class of the dockable window add-in. It is responsible for          /// creating and disposing the user interface class of the dockable window.         /// </summary>         public class AddinImpl : ESRI.ArcGIS.Desktop.AddIns.DockableWindow         {             private FrmWellConfig m_windowUI;                  public AddinImpl()             {             }                  protected override IntPtr OnCreateChild()             {                 m_windowUI = new FrmWellConfig(this.Hook);                 return m_windowUI.Handle;             }                  protected override void Dispose(bool disposing)             {                 if (m_windowUI != null)                     m_windowUI.Dispose(disposing);                      base.Dispose(disposing);             }              }


Is it possible to export boundaries from Google Map marker?

I primarily have a doubt that is it possible to export a marked boundaries from google map marker ?if possible kindly tell how to to do this?


Add Surface Information in 3d-analyst

One of our guys is using the "Functional Surface->Add Surface Information" feature in 3D-analyst tool in ArcMap 10.1 to add elevation data to an XY list. It's the only feature we're using from the 3D-analyst. Our trial period has ended and we're trying to figure out if we need to buy the whole 3D-analyst tool for $2500 just for this feature, or if there's a cheaper way to do this task since we don't use anything else in 3D-analyst. Since it's a pretty simple task, I might just develop a tool myself to do this using the USGS web service given the cost of 3D-analyst. Is there a cheaper / easier way to do this?


How to create and re-use my own map template in Arcmap 10.1?

I am working on a project which requires 5 different analyses on the same data. I created five map documents (.mxd) and completed my analyses. Now, the problem is that I want to export my results in map format but I want them to be uniform in display i.e. having a same map template. Is there a way to save one of the maps as a template and then use it for creating other 4 maps? I tried saving one document in 'Map Templates' folder in ArcGis folder but a message box appears which says that I don't have permission to save in that location and I must contact administrator for that, even though I am the only user and administrator of the computer. I don't want to manually add all map elements for 5 different maps, thats too tedious!


How to reopen a plugin-window in QGIS?

I am currently taking my first steps in playing around with programming a plugin in python 2.7 for QGIS 1.8. While my plugin does what it should, I would like to make it more userfriendly. At the moment, the user is expected to choose an input- and an output-file, which are then processed. In case the user forgets to choose one (or both) of the files, the plugin is displaying an information window - and then quits completely after the user hits "OK".               if src_filepath =="" and dst_filepath =="":         QMessageBox.information(None, "DEBUG:","Please specify files!")`      However, I would rather like to reload the main plugin-window with the "choose file"-dialog. I suspect that the solution is rather easy, but I was not able to invoke this window, so any help is appreciated. The interface code derives more or less directly from the "Plugin Builder" in QGIS.


Force dimensionality with shp2pgsql

After invoking: `shp2pgsql -W "latin1" theshape.shp > thesql.sql` My console responds with:               Shapefile type: PolygonZ     Postgis type: MULTIPOLYGON[2]      Note the shapefile type is `PolygonZ`. The result is, I have to manually edit the `thesql.sql` file and wrap the geometry column in `ST_Force_2D()` before executing the SQL, otherwise, Postgres complains I'm trying to import a 3D shape into a 2D column. How can I force dimensionality to 2 with `shp2pgsql`? Alternatively, can I force it in QGIS?


Upload layers in Geoexplorer

I have OGS 4.02, OS is Windows Server 2008 R2 Standard.I have questions about Geoexplorer. I know that anonymous users can not save and publish map and also they can not upload layers. If I create users which are in Geoserver admins (access mode is admin for specific workspace) or editors (access mode is write for specific workspace), why they can not upload their layers (zipped shapefiles) in their workspace, creating maybe even new stores on the fly. Why is that option only avaliable for one person, user with admin role? Is it something with folder security in Windows, is it necessary to create also users and specific folders for them. I know that after succesfull upload, two folders are created: one in C:\ProgramData\Boundless\OpenGeo\geoserver\uploads (tmp prefix) and another in workspace where data are uploaded. So who can only write in that folders through geoexplorer and geoserver port 8080? Is it necessary to create postgis store to enable upload for some power users (editors or small admins)?


TableToTable gives unexpected character after line continuation character?

I'm working on updating a python script from arcmap 9.3.1 to 10.2.1 and im getting an error for "unexpected character after line continuation character" and im unable to diagnose the character.               arcpy.TableToTable_conversion(sde_parcel_base_temp, sde, "parcel_base", "", "pin \"PIN\" true true false 10 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",PIN,-1,-1;pwd \"PWD\" true true false 13 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",PWD,-1,-1;add_num \"Address Number\" true true false 10 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",add_num,-1,-1;add_prefix \"Address Prefix\" true true false 1 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",add_prefix,-1,-1;add_strt \"Address Street\" true true false 50 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",add_strt,-1,-1;add_type \"Address Type\" true true false 10 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",add_type,-1,-1;add_sdir \"Address Suffix\" true true false 2 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",add_sdir,-1,-1;add_munic \"Address Municipality\" true true false 50 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",prop_city,-1,-1;add_zip \"Address ZIP\" true true false 10 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",prop_zip,-1,-1;add_full \"Full Street Address\" true true false 50 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",prop_street,-1,-1;ownr_name1 \"Owner Name 1\" true true false 50 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",owner1,0,49;ownr_name2 \"Owner Name 2\" true true false 50 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",owner2,0,49;ownr_addr \"Owner Address\" true true false 50 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",own_street,-1,-1;ownr_city \"Owner City\" true true false 50 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",own_city,-1,-1;ownr_state \"Owner State\" true true false 2 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",own_state,-1,-1;ownr_zip \"Owner Zip\" true true false 10 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",own_zip,-1,-1;acres "\Legal Acreage\" true true false 50 Double 0 0 ,First,#,"+ sde_parcel_base_temp +",legal_acreage,-1,-1;prop_code \"Property Class Code\" true true false 10 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",property_class,-1,-1;nghbr_code \"Neighborhood Code\" true true false 10 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",neighborhood,-1,-1;lgl1 \"Legal Line 1\" true true false 80 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",lrsn,-1,-1;lgl2 \"Legal Line 2\" true true false 80 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",lrsn,-1,-1;lgl3 \"Legal Line 3\" true true false 80 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",lrsn,-1,-1;lgl4 \"Legal Line 4\" true true false 80 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",lrsn,-1,-1;lgl5 \"Legal Line 5\" true true false 80 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",lrsn,-1,-1;lgl6 \"Legal Line 6\" true true false 80 Text 0 0 ,First,#,"+ sde_parcel_base_temp +",lrsn,-1,-1;lrsn \"Proval LRSN\" true true false 50 Long 0 0 ,First,#,"+ sde_parcel_base_temp +",lrsn,-1,-1;updated \"Last Update\" true true false 50 Date 0 0 ,First,#,"+ sde_parcel_base_temp +",last_update,-1,-1", "")      Above is the part of the script giving me problems. Anybody with a trained eye able to point out my problem for me?


Create a subset of a dataset with a specific grid size in ArcGIS using Data Management Tool

I am wondering how you can create a subset of a dataset that is a specific grid size (in pixels) using the Clip tool within the Data Management Tool box for raster data. I was trying to find a hint where to set the grid size in the tool's Environment Settings, but the result, after running the tool, did not apply the new grid size. Please does anyone handled this before? Thanks.


SQL Command to convert pgRouting result to GeoJSON

My goal is to produce a Web UI for pgRouting using this tutorial. I have already imported my data from an OSM pBF to PostGIS using osm2po. Using the pgRoutingLayer plugin for QGIS, I was able to figure out the SQL code to create routes using driving time as the cost.               SELECT ways.*, route.cost AS route_cost FROM ways JOIN (         SELECT * FROM shortest_path(         'SELECT gid AS id,                  source::int4 AS source,                  target::int4 AS target,                  cost::float8 AS cost,                 reverse_cost::float8 AS reverse_cost FROM ways',                 62,19169, True, True))          AS route ON ways.gid= route.edge_id;      My next step is to convert this into GeoJSON as done in Part 1 of the tutorial. This is what seemed logical:               SELECT * INTO export     FROM ways.*, route.cost AS route_cost FROM ways JOIN (         SELECT * FROM shortest_path(         'SELECT gid AS id,                  source::int4 AS source,                  target::int4 AS target,                  cost::float8 AS cost,                 reverse_cost::float8 AS reverse_cost FROM ways',                 62,19169, True, True))          AS route ON ways.gid= route.edge_id;     SELECT ST_AsGeoJSON(ST_UNION(the_geom)) AS geom_union     FROM export;          ERROR:  syntax error at or near ","     LINE 2:  FROM ways.*, route.cost AS route_cost FROM ways JOIN (                         ^      but I can't get it to work. Honestly, I'm not very well acquainted with SQL, can someone point out my mistake? edit: This seems to work, but does it do the same thing?               CREATE TEMPORARY TABLE route_result AS      SELECT ways.*, route.cost AS route_cost FROM ways JOIN (         SELECT * FROM shortest_path(         'SELECT gid AS id,                  source::int4 AS source,                  target::int4 AS target,                  cost::float8 AS cost,                 reverse_cost::float8 AS reverse_cost FROM ways',                 62,19169, True, True))          AS route ON ways.gid= route.edge_id;      SELECT ST_AsGeoJSON(ST_UNION(the_geom)) AS geom_union      FROM route_result;


How would I get number of tiles loaded, and total number loading?

So new to OpenLayers 3, and I'm trying to create a "loading bar" to indicate how much of a layer has loaded. I'm thinking of making it based on the number of tiles loaded so far divided by the total number of tiles that need to be loaded to render that whole layer from the WMS server-- but I've been looking around the API and I don't know how to get either of those numbers in OpenLayers 3. Does anyone know anything about this, or alternatively does anyone know a different approach I could use to create a layer loading bar? Thanks in advance! :]


Running a batch script in GRASS command line - v.dissolve and v.db problems

I have been writing a number of scripts for automating isochrone production in GRASS just using the command line. I create files with isochrone lines in them and then a buffer around each of the lines. I then want to dissolve the buffers, so use the command: v.dissolve input=AOS_buffer_XXX output=AOS_buffmerge_XXX Each time this completes it's first file, the script stops. I then want to put an attribute into each line that reflects the filename and I thought I could do this with the command: v.db.addcol map=AOS_buffmerge_XXX columns=GTW_ref INT [--verbose then v.db.update map=AOS_buffmerge_XXX column=GTW_ref value=XXX [--verbose] However when I run this on the several hundred shapes, it ends after the first execution - it seems to allow me to do one, but not significant numbers. Can anyone help me out with why these three functions in GRASS command line won't script? (I'm using GRASS 6.4.2 - would an upgrade help?) Example of the script (I have put annotation where the script just stops - no error message, it just stops):               v.in.ogr dsn=C:\0Data\FK_openspace\network\cleanedfalkirknetwork.shp output=fk_net009 --verbose -–overwrite               v.in.ogr dsn=C:\0Data\FK_openspace\OpenSpaceAndGates\Calculations\AmenityOpenSpace.shp output=AOS_gates1 where="Reference = 1" --verbose --overwrite     v.in.ogr dsn=C:\0Data\FK_openspace\OpenSpaceAndGates\Calculations\AmenityOpenSpace.shp output=AOS_gates2 where="Reference = 2" --verbose --overwrite     v.in.ogr dsn=C:\0Data\FK_openspace\OpenSpaceAndGates\Calculations\AmenityOpenSpace.shp output=AOS_gates3 where="Reference = 3" --verbose --overwrite     v.in.ogr dsn=C:\0Data\FK_openspace\OpenSpaceAndGates\Calculations\AmenityOpenSpace.shp output=AOS_gates4 where="Reference = 4" --verbose --overwrite          v.net input=fk_net009 points=AOS_gates1 output=AOS_fixednet1 operation=connect alayer=1 nlayer=2 thresh=500 --verbose --overwrite     v.net input=fk_net009 points=AOS_gates2 output=AOS_fixednet2 operation=connect alayer=1 nlayer=2 thresh=500 --verbose --overwrite     v.net input=fk_net009 points=AOS_gates3 output=AOS_fixednet3 operation=connect alayer=1 nlayer=2 thresh=500 --verbose --overwrite     v.net input=fk_net009 points=AOS_gates4 output=AOS_fixednet4 operation=connect alayer=1 nlayer=2 thresh=500 --verbose --overwrite          v.net.iso input=AOS_fixednet1 output=AOS_output1 ccats=0-99999999 costs=400 --verbose --overwrite     v.net.iso input=AOS_fixednet2 output=AOS_output2 ccats=0-99999999 costs=400 --verbose --overwrite     v.net.iso input=AOS_fixednet3 output=AOS_output3 ccats=0-99999999 costs=400 --verbose --overwrite     v.net.iso input=AOS_fixednet4 output=AOS_output4 ccats=0-99999999 costs=400 --verbose --overwrite          v.out.ogr input=AOS_output1 type=line dsn=C:\0Data\FK_openspace\Outputs\advanced_outputs\AOS layer=1 format=ESRI_Shapefile --verbose      v.out.ogr input=AOS_output2 type=line dsn=C:\0Data\FK_openspace\Outputs\advanced_outputs\AOS layer=1 format=ESRI_Shapefile --verbose      v.out.ogr input=AOS_output3 type=line dsn=C:\0Data\FK_openspace\Outputs\advanced_outputs\AOS layer=1 format=ESRI_Shapefile --verbose      v.out.ogr input=AOS_output4 type=line dsn=C:\0Data\FK_openspace\Outputs\advanced_outputs\AOS layer=1 format=ESRI_Shapefile --verbose           v.in.ogr dsn=C:\0Data\FK_openspace\Outputs\advanced_outputs\AOS\AOS_output1.shp output=AOS_line_1 where="cat = 1" --verbose --overwrite     v.in.ogr dsn=C:\0Data\FK_openspace\Outputs\advanced_outputs\AOS\AOS_output2.shp output=AOS_line_2 where="cat = 1" --verbose --overwrite     v.in.ogr dsn=C:\0Data\FK_openspace\Outputs\advanced_outputs\AOS\AOS_output3.shp output=AOS_line_3 where="cat = 1" --verbose --overwrite     v.in.ogr dsn=C:\0Data\FK_openspace\Outputs\advanced_outputs\AOS\AOS_output4.shp output=AOS_line_4 where="cat = 1" --verbose --overwrite          v.buffer input=AOS_line_1 output=AOS_buffer_1 type=line distance=35 --overwrite --verbose     v.buffer input=AOS_line_2 output=AOS_buffer_2 type=line distance=35 --overwrite --verbose     v.buffer input=AOS_line_3 output=AOS_buffer_3 type=line distance=35 --overwrite --verbose     v.buffer input=AOS_line_4 output=AOS_buffer_4 type=line distance=35 --overwrite --verbose          v.dissolve input=AOS_buffer_1 output=AOS_buffmerge_1 --verbose     v.dissolve input=AOS_buffer_2 output=AOS_buffmerge_2 --verbose     v.dissolve input=AOS_buffer_3 output=AOS_buffmerge_3 --verbose     v.dissolve input=AOS_buffer_4 output=AOS_buffmerge_4 --verbose          v.db.addcol map=AOS_buffmerge_1  columns=GTW_ref INT --verbose      v.db.addcol map=AOS_buffmerge_2  columns=GTW_ref INT --verbose     v.db.addcol map=AOS_buffmerge_3  columns=GTW_ref INT --verbose     v.db.addcol map=AOS_buffmerge_4  columns=GTW_ref INT --verbose          v.db.update map=AOS_buffmerge_1 column=GTW_ref value=1 --verbose     v.db.update map=AOS_buffmerge_2 column=GTW_ref value=2 --verbose     v.db.update map=AOS_buffmerge_3 column=GTW_ref value=3 --verbose     v.db.update map=AOS_buffmerge_4 column=GTW_ref value=4 --verbose          v.out.ogr input=AOS_buffmerge_1 type=area dsn=C:\0Data\FK_openspace\Outputs\advanced_outputs\AOS\final layer=1 format=ESRI_Shapefile --verbose --overwrite     v.out.ogr input=AOS_buffmerge_2 type=area dsn=C:\0Data\FK_openspace\Outputs\advanced_outputs\AOS\final layer=1 format=ESRI_Shapefile --verbose --overwrite     v.out.ogr input=AOS_buffmerge_3 type=area dsn=C:\0Data\FK_openspace\Outputs\advanced_outputs\AOS\final layer=1 format=ESRI_Shapefile --verbose --overwrite     v.out.ogr input=AOS_buffmerge_4 type=area dsn=C:\0Data\FK_openspace\Outputs\advanced_outputs\AOS\final layer=1 format=ESRI_Shapefile --verbose --overwrite


OpenLayers.Format.GML parser problem

I'm trying to parse WFS reponse using `OpenLayers.Format.GML` parser. However, an exception thrown in my try-catch block is `TypeError: a is undefined`. Any pointers as to what I may be doing wrong will be appreciated.               protocol = new OpenLayers.Protocol.WFS({         version: "1.1.0",         url: WFS_HOST,         featureType: "pl",         featureNS: "ws",         srsName: "EPSG:900913"     });          wfsLayer = new OpenLayers.Layer.Vector("PE Vector", {         strategies: [new OpenLayers.Strategy.BBOX()],         protocol: protocol     });          var _Callback = function(resp) {         console.log(resp);         console.log(resp.error);              try {             var gmlParser = new OpenLayers.Format.GML();             gmlParser.extractAttributes = true;             var features = gmlParser.read(resp.responseText);                  if(features){                 for(var i = 0; i < features.length; i++) {                     var feature = features[i];                                         // More operations                 }             }         } catch(e) {             alert("Error: " + e);         }     };          response = protocol.read({         maxFeatures: 100,         callback: _Callback     });      My environment: GeoServer 2.1.3, OpenLayers 2.11, PostGIS 1.5.


How to create a hole in a polygon in QGIS?

I'm trying to use the tool "delete ring" in order to create a hole in a polygon. The geometry of the hole is another feature of the same shapefile I cannot manage to delete the part of the polygone. Does someone have a tutorial or some explanation?


ECW on QGIS 1.8 (Ubuntu 11.10)

I'm trying to deal with an ECW file with QuantumGis 1.8 Lisboa, on Ubuntu 11.10. Obviously it tells me that .ecw is not a supported kind of file on QuantumGis. I have tried to install all the libraries but it still doesn't work. What sould I do? Thanks, Jacopo


Using render_list to render only map extract

I am trying to render an area where I have added an extract from the planet file with render_list. The command I am using is               render_list -x 159 -X 179 -y 160 -Y 410 -n 2 --socket=/var/run/renderd/renderd.sock --min-zoom=10 --max-zoom=10      to render zoom level 10. I already rendered 0-9 zoom levels for the whole planet, but want to render only the areas where I imported maps for the bigger zoom levels. Running the command above gives message "Starting 2 rendering threads" and stays for a long time without doing anything. Any ideas why?


Postgis nested query trouble

I am having difficulty writing a nested postgis query. I have a single table from which I need to create a view in PostgreSQL. I have created the view with the following command               create view rrop08.rrop08_lws_f AS     select * from rrop08.rrop08_gk3_polygon      where planznr='04.01' OR planznr ='04.02';      The problem here is that the table contains another field called "bez". I need to carry out the first query but also with bez=0; so that none of the 04.01 or 04.02 rows have a bez value other than 0. I have tried this...               create view rrop08.rrop08_lws_f AS     select * from rrop08.rrop08_gk3_polygon      where planznr='04.01' OR planznr ='04.02' AND bez=0;      but this retrieves all rows with bez=0 and not just the 04.01 or 04.02 ones. I need to first select 04.01 or 04.02 and then on this query get those with bez=0...ie.nested. I´m a bit lost on how to do this. THanks for any help, Rob


How can i resolve problem with bounds for export?

How could I set up bounds in TileMill if i don't want to upload the country around? I simply want to upload my created island in `TileMill` to `Mapbox`. Is there a way to do it? I am using OSM Bright for the design of my map.


Calculate Field python code has "Expected end of statement" error

I am creating a line shapefile from a csv in a Python (2.7) script. I am assigning an ID field from the csv to a variable and then later using it in the **Calculate Field** tool. The ID field originally looks like this '011947-OCT-084' which ArcGIS(10.1) does not seem to like, so I removed the '-'s and when I use that variable in the tool I get an error: > ExecuteError: ERROR 999999: Error executing function.   >  Expected end of statement   >  Failed to execute (CalculateField). Here is my code ("temp" is the shapefile temporarily stored in in_memory):               with open(inCSV, 'rb') as input:         with open(outCSV, 'wb') as output:             csvwriter = csv.DictWriter(output, headers, delimiter = ',')             csvwriter.writeheader()             for row in csv.DictReader(input):                 uniqueID = row['Common_id']                 newID = uniqueID.replace("-", "")                 arcpy.AddField_management("temp", "uniqueID", "TEXT")                 arcpy.CalculateField_management("temp", "uniqueID", newID)                 row['Common_id'] = newID                 csvwriter.writerow(row)      I'm not entirely sure what the cause behind this is, as I've used variables in the Calculate Tool field before without a problem. Thoughts?


IL East 1201 coordinate system

What is the difference between NAD CORS96 (Illinois East 1201) and the "normal NAD1983? We have data that was GPS'ed in CORS96, but data was imported into ArcMap using the "normal" NAD1983? What effect, if any, will it have on data? For example, will it affect the location and the area calculations of the data?


Does arcpy have a method for accessing related records within a cursor?

I have a line feature class that represent sampling segments. I have a data table of field notes in the same Geodatabase that is connected to the feature class with a Relationship Class. Is there a built in method for arcpy to access related records? The work around would be to add another cursor for the second table, but you'd think that if you take the time to relate your data there would be a quick and easy method for accessing related records programatically. Software: ArcGIS 10.1


Interpolation iteration using ArcGIS- arcpy

My current project requires interpolation of state level weather data. Approximately 5 years of daily data needs to be interpolated. I still need to figure out the best interpolation method. However, for this exercise I chose IDW. I found a script in ArcGIS forum and tried to use it. However, I am not sure why the script is not working. Any help would be greatly appreciated.                #-------------------------------------------------------------------------------         # Name:        Inverse Distance Method Interpolation         # Purpose:         #         # Author:               #         # Created:     19/09/2013         # Copyright:            # Licence:              #-------------------------------------------------------------------------------         # ---------------------------------------------------------------------------         # rainfall_idw2012.py         # Created on: Tue Apr 08 2008 10:47:39 AM         # ---------------------------------------------------------------------------                   # Import system modules         import arcpy         import sys         import os         import arcgisscripting         import string              arcpy.env.overwriteOutput =  True              # Create the Geoprocessor object         gp = arcgisscripting.create()              # Check out any necessary licenses         gp.checkoutextension("Spatial")              # Load required toolboxes...         gp.AddToolbox(r"C:\Program Files (x86)\ArcGIS\Desktop10.2\ArcToolbox\Toolboxes\Spatial Analyst Tools.tbx")              #Output_variance_of_prediction_raster = ""         modelPath = r"D:\Geodatabase\Rainfall\r2012"                   # Setting workspace         gp.workspace = modelPath                   # Input Variables         stations = r"D:\Geodatabase\MyGeodatabase.gdb\mr2012"         studyArea = r"D:\Geodatabase\MyGeodatabase.gdb\boundary"                   try:             # Get list of dau=ily events for looping             daily_events = gp.ListFields(r"D:\Geodatabase\MyGeodatabase.gdb\mr2012", "r2012*")                  daily_events.reset()                  # Get the first reading and start loop                  dailyevent = daily_events.Next()                  while dailyevent:                      rain_amount = dailyevent.Name                      # Build output data layer name from the rainfall:                      r_level = modelPath + _rain_amount                      extract = modelPath + _rain_amount                 dailyevent = daily_events.Next()                      gp.AddMessage ("Current Operation: Preparing Interpolating using IDW")                      # Process: IDW interpolation                      gp.IDW_sa(stations, rain_amount, r_level, "1000", "2" "VARIABLE 12", "")                                # Process: Extract by Mask...                 gp.ExtractByMask_sa(r_level, studyArea, extract)              except:             gp.AddMessage("IDW failed " + gp.GetMessages())             print "IDW failed ", gp.GetMessages()


Add only unique features from another shapefile

I've got a problem that I am sure at least one of you have dealt with. I'm new to a company that makes use of a rather large database of point data. For whatever reason, the woman who was here before me made a huge amount of duplicate points - the attributes may be a little different but the coordinates are the same. I've found that I can remove the points either in QGIS using the mmqgis plugin, or in ArcMap using the "delete identical" tool. However, when I use either method, I have noticed that there are some points that are removed completely. For most, it leaves one point but I have found a few that are completely gone. Do any of you know of a method to add unique features to an attribute table? I'm thinking that I can compare the original table to the "cleaned" one and add the missing points back. Thank you so much!


EXCEPTION - bad allocation-printing

I have a street map, from Durango, Mexico. And I want to print it. But Qgis show the next error " EXCEPTION bad allocation". I´m going to print it in 1.50 x 1.06 mts sheet. What can I do? Greetings Roberto Guerrero Durango, Dgo., México.


ArcObjects 10/10.1 Get rid of Invalid layer symbol in TOC

I've written an add-in (VS2010, C#) that finds specific invalid feature layers in ArcMap and fixes them. After the layer is fixed, the red invalid symbol still appears in the TOC. An activeview.refresh doesn't take care of it. How does one programatically get that layer entry to appear as valid without manual intervention?


Google Maps v3 with Openlayers - First Zoom Level

I use Google Maps v3 in my application and I noticed a strange behavior, When I refresh the web page, the map zoomed out to display the whole map then get back to the zoom level I specify in my code. This happen fast but every time I refresh the page. How can I fix it? **EDIT** Part of mu code:               var geographic = new OpenLayers.Projection("EPSG:4326");     var mercator = new OpenLayers.Projection("EPSG:3857");     maxExtent: new OpenLayers.Bounds(-20037508, -20037508, 20037508, 20037508.34).transform(geographic,mercator);          map = new OpenLayers.Map('map',{             projection: new OpenLayers.Projection('EPSG:3857'),              displayProjection: new OpenLayers.Projection("EPSG:4326")      });           gmap = new OpenLayers.Layer.Google("Google",{ numZoomLevels: 20});        map.addLayers([gmap]);


How to add text-only labels on Leaflet map with no icon

I'm looking for a way to show text on Leaflet map using markers or anything else without showing any icon. I want to show text only and be able to style and rotate it... Any suggestion? Thanks a lot Etienne


How to get an attribute of the selected feature and pass it into a URL?

As a total Novice with Python, I thought 2 weeks would be enough time to learn how to build some tools for our GIS... I was wrong and I'm running up on a deadline. I am attempting to get a value from a selected feature in geodatabase, and pass that value into a URL link / query after the equals sign. Opening the link to the image viewer works, if I don't try to pass in the field value. After 9 days of reading training manuals my eyes are raw. How can I fix this...               import arcpy     import pythonaddins     import webbrowser     from threading import Thread          def GetParameterAsText (str.Unit_No):     "This prints a passed through string into this function"          import str     row.getValue(Unit_No)          def OpenBrowserURL ():         url = 'http://fcweb/caviewer#unit='         webbrowser.open(url,new=2)          class CAVButton(object):     """Implementation for PythonButton_addin.CAVbutton (Button)"""     def __init__(self):         self.enabled = True         self.checked = False     def onClick(self):             t = Thread(target=OpenBrowserURL)             t.start()             t.join()


How can I set the settings of an Address Locator within a Model or arcpy?

I am creating a model that processes a large road dataset, then creates an address locator from the roads. The Create Address Locator tool does not allow me to set options like Spelling Sensitivity or Minimum Candidate Score. **Is it possible to set these options in Model Builder (or at least, within an arcpy Python script)?** I'm using ArcGIS 10, but I'm also interested in whether it's possible in 10.1.


Click to Hover for WMSGetFeatureInfo

I am using the following code and wanted to know how I can change it from click event to hover or mousemove event. I know there is a eventListner but I am unable to locate the click event where I must change it to mousemove event.               info = new OpenLayers.Control.WMSGetFeatureInfo({                 url: 'http://localhost:8080/geoserver/iPlant/wms',                  title: 'Identify features by clicking',                 queryVisible: true,                 eventListeners: {                     getfeatureinfo: function(event) {                         map.addPopup(new OpenLayers.Popup.FramedCloud(                             "chicken",                              map.getLonLatFromPixel(event.xy),                             null,                             event.text,                             null,                             true                         ));                     }                 }             });             map.addControl(info);             info.activate();      **Update 2.0** I made the changes and also added the PopUp, but only the click function is being recognized.               info = {                 click: new OpenLayers.Control.WMSGetFeatureInfo({                 url : 'http://niles.iplantcollaborative.org:8080/geoserver/wms',                 title : 'Identify features by clicking',                 queryVisible : true,                 eventListeners : {                 getfeatureinfo : function(event) {                 deletePopUp();                 popup = new OpenLayers.Popup.FramedCloud("popInfo", map.getLonLatFromPixel(event.xy), null, event.text, null, true);                 map.addPopup(popup);}}                 }),                 hover: new OpenLayers.Control.WMSGetFeatureInfo({                 url : 'http://niles.iplantcollaborative.org:8080/geoserver/wms',                 title : 'Identify features by clicking',                 queryVisible : true,                         hover : true,                 eventListeners : {                 getfeatureinfo : function(event) {                 deletePopUp();                 popup = new OpenLayers.Popup.FramedCloud("popInfo", map.getLonLatFromPixel(event.xy), null, event.text, null, true);                 map.addPopup(popup);                 }}                 })                 };                 map.addControl(info.hover);                 info.hover.activate();


Find the intersection of two rasters

I have two raster images, Boolean classified (1 or 2). They overlap in given areas. I need a raster that is only the overlapping sections. Using Arcinfo 10. Answer can only use raster methods, no converting to polygons, points, etc...![Just for reference here's a screen shot. ](http://i.stack.imgur.com/9JseK.jpg)


OSM Country import - missing border-near lakes, rivers

I'm importing a country (e.g Switzerland) from the cloudmade OSM branch and I'm failing to get the bordercrossing lakes on display. The same issue came up, when I tried to import the geofabrik .dpf branch this morning. Somehow it makes sense and I don't think it's a bug, but I struggle to get the border crossing stuff to work. How to solve this issue without importing the european-western fileset? Thanks! Martin


Seeking a GIS for indoor mapping and direction finding

I've seen a few options for indoor GIS platforms such as Google's indoor mapping, implementations of OpenLayers, Open Street Map, and many other plausible solutions. There's even a post on here asking a very similar question to mine but all the answers involved how to provide locations to users. For my use I don't care to map where the user is but rather simply provide a usable map of our building and in my ideal world I would be able to provide static directions (think Google Maps without location for our building). Is there anything out there that can help us achieve this? I'm willing to learn as long as there's documentation.


proj4 cookbook?

Is there a proj4 "cookbook"? I discovered proj4 from @MerseyViking excellent answer to: Create Mercator map with arbitrary center/orientation? but most of the documentation appears to be fairly high level. EDIT: Details on what I'm trying to do:   * Rotate the globe so that any point I choose is at 0,0 (where the prime meridian crosses the equator) or the North Pole or the South Pole, my choice.    * Apply any projection with any "zoomlevel" I choose to the result.    * "Zoomlevel" meaning that the projection includes only points close to the point I chose earlier.  The linked question shows me how to do this for the Mercator projection. How to generalize? I realize I could write code to do 3D rigid rotations (for step 1 above), but can proj/cs2cs do this for me?


How to import data with SRID 4269 to GEOGRAPHY type?

In PostGIS, we are importing various US census bureau shapefiles (with SRID = 4269) using **shp2pgsql** loader. I have medium-sized database with several tables containing fields of type geography (lat/lon points, polygons, multipolygons, etc) that span locations across the US on which we do various distance calculations, “within” polygon, etc type of geographic functions as well as display data using Google Maps Api. We currently have no geometry type fields (only SRID used in database is 4326 via geography type). Because they cover a fairly large area and to keep things simple and consistent (avoid projection details etc), I would like to store the census shapes (ie boundaries for states, metros, counties,block groups..) in geography data type instead of geometry as well (unless there is a reason NOT to do this or someone can explain why we should store as geometry instead). My question is, if I use shp2pgsql loader with –G option to create a MULTIPOLYGON geography field for these boundaries (and thus use a 4326 SRID) does the loader “convert” the data from a 4269 projection to 4326? Is a conversion really necessary? Do I need to first load into temp table with geometry type then perform update to transform/cast into geography type and then store that result? Thank you for any advice or insight related to this.


How to programatically populate a shapefile field with polygon areas in OGR?

I'm trying to programatically populate a shp field with polygon areas:               import sys     import ogr          ds = ogr.Open( 'tttttttttt.shp', update = 1 )     if ds is None:         print "Open failed./n"         sys.exit( 1 )          lyr = ds.GetLayerByName( "tttttttttt" )     lyr.ResetReading()          field_defn = ogr.FieldDefn( "Area", ogr.OFTReal )     lyr.CreateField(field_defn)          for i in lyr:         # feat = lyr.GetFeature(i)          geom = i.GetGeometryRef()         area = geom.GetArea()         print 'Area =', area         lyr.SetFeature(i)         i.SetField( "Area", area )          ds = None      But I get an empty values in the field... Thanks


Automating attribute completion of new polygons

Good Afternoon I am currently outlining property boundaries and then copying and pasting in the necessary attributes on a daily basis. I would like to speed up this process and remove the chance of error by introducing some automation. From my database I am able to export (to an excel file), all of the records for the coming day. I would like some kind of drop down box or search function which allows me to select the correct record from the excel file for the property boundary which I have just created. I would like to achieve this using a python script. However, I have searched online for related problems or tutorials and cant find any solutions or similar problems. I would like some pointer about where to start with this and how to approach it. Are there any specific arcgis or python tools/modules which can help me achieve this. Many thanks


Install both postgis and postgresql-9.1-postgis packages?

I've installed the PostgreSQL package (ver. 9.1) on a Debian 7 server using the command "apt-get install postgresql postgresql-contrib". I now want to install PostGIS as well. Should I install both the "postgis" package and the "postgresql-9.1-postgis" package or just the latter? It's not clear from the descriptions on their Debian wiki pages what the differences are between the two packages. I ask because I've seen some blog posts where the author has installed both packages but yet on the OpenStreetMap PostGIS wiki page, they just install the postgresql-9.1-postgis package. Yet this Debian page seems to imply that you should install both packages as postgis contains common PostGIS files while postgresql-9.1-postgis contains files that specifically support PostgreSQL 9.1 Thanks.


Create custom projection using either QGIS (or ArcGIS)?

I have searched for an answer to this but don't see anything relating to what I am trying to do. What I have are a series of Geographic Coordinates, Latitude, Longitude, Ellipsoid heights. What I want are a local system which is translated from NAD83/UTM. I know the parameters of the shift which works with UTM Zone 9N. My parameters are: ****Origin North: 5984910.533m, Origin East: 519681.595m**** ****Trans North: 4973.397m, Trans East: 2068.553m**** ****Rotation: 10°12’53”, Scale: 1.0003933608**** ****Trans Height: ‐1.907m, Slope North: 0.00001, Slope East: 0.00004**** This should result in coordinates similar to N = 5000.000m, E = 1500.000m, Elev = 0.000m Is there any good help for the EPSG notation? Perhaps I need to steal from the UTM projection and combine with another? Thanks, Sean Here is some data in Geographic(WGS84/NAD83) and Local ground **Geographic Coordinates Latitude, Longitude, Ellipsoid Height(meters) (Deg.MMSSsssss)** ie (54.00595628 = 54 degrees 00 minutes 59.5628 seconds) 1 54.00595628 -128.4231334 30.48 MON 6 54.0059556 -128.4154351 14.607 MON 9 54.00451717 -128.4150654 13.147 MON 19 54.00019706 -128.4139788 9.238 MON 20 54.00133493 -128.415171 13.803 MON **Local ground system Northing, Easting, Elevation (meters)** 1 5562.562 1570.005 28.547 CONC PILLAR MON 6 5445.703 2233.204 12.713 CONC PILLAR MON 9 4996.025 2222.510 11.250 CONC PILLAR MON 19 3646.242 2186.193 7.329 CONC PILLAR MON 20 4030.357 2033.239 11.848 CONC PILLAR MON


Alternatives to using Arcpy

I seem to use ESRI's Arcpy site package for virtually all of my python geoprocessing. To ESRI's credit, these are an incredible suite of tools that can help accomplish a great deal. However, I would also like to create geoprocessing scripts outside of the ESRI Arcpy domain. For example, if I want to clip a raster to a polygon, I would start with the following script from ESRI:               # Import system modules     import arcpy     from arcpy import env     from arcpy.sa import *          # Set environment settings     env.workspace = "C:/sapyexamples/data"          # Set local variables     inRaster = "elevation"     inMaskData = "mask.shp"          # Check out the ArcGIS Spatial Analyst extension license     arcpy.CheckOutExtension("Spatial")          # Execute ExtractByMask     outExtractByMask = ExtractByMask(inRaster, inMaskData)          # Save the output      outExtractByMask.save("C:/sapyexamples/output/extractmask")        I'm not sure how I would accomplish the same task programmatically without Arcpy. My questions for the serious programmers out there: **What collection of python tools do you use to accomplish tasks that ESRI users would accomplish with the Arcpy site package? Where do I begin?**


Can Geomedia be automated using VBA

Can I automate Geomedia using VBA? The Geomedia training manual "Fundamentals of Geomedia Development Part 1" talks about automating using visual basic but I do not have Visual Basic.


Are there only two methods of interpolation in QGIS?

I did recently a map of maximum temperature for specific area in Quantum GIS 1.8.0 and its interpolation. However, I was wondering if there are only two method of interpolation, apart from ‘Triangulated Irregular Network (TIN)’ or ‘Inverse Distance Weighted (IDW)’. Maybe is there any extra plugin for that?


Find the next vertices at the intersect line

I'm Trying to calculate the angle between two intersected lines at geometric network (simple edge)in order to check if the lines are close to be perpendicular. I already have an Array of the intersection vertices and a function that calculate the angle using Pythagorean theorem. How I can get the next vertices at the intersect line, in order to calculate the angle between this vertices and the intersect point using the function ?


SQL to get points within certain distance

I got about 2 million points in a spatialite database would like to get which points are withina radius of 30 km from another point. All data is in WGS84 but I need distance in km. Also it would be handy to summarize by some field I tried spatial query in QGIS but crashed due to the size


rasters as independent variables in logistic regression: must they have the same values classification?

I am planning to build a predictive model for land productivity using arcGIS 10.1. and logistic regression. In my model, the Independent Variables would be, among others, distance from rivers, distance from spring-line, distance from the coastline. I am wondering if it is sounder/useful/essential to use the same value classification for all the three aforementioned rasters, or if it would not pose any problem in analytical perspective to have rasters with different classification.


Creating contiguous / adjacent ellipses

I've only just started using drawing tools in ArcGIS, and have the following question I haven't found the answer to online. Is there a way of taking a drawn ellipse (at a slight angle) and adding say 50-100 of the same ellipse, arranging them as efficiently as possible so that there's as little blank space as possible around them (but no overlaps)? I have tried doing this by hand but I can't get the snapping tools to work properly (plus it's pretty time consuming), so ideally looking for a way to get this done quickly and introducing little/no human error into the process.


ODBC connection issue on 64bit architecture

I'm migrating an addin to a 64bit machine using Arcmap 9.3.1. Upon execution, it returns the following error - "[Microsoft][ODBC Driver Manager] The specified DSN contains an architecture mismatch between the Driver and Application". I'm aware through sleuthing that there exists two versions of the odbcad32.exe; one in the System32 folder, and one in sysWOW64. Oracle 11g connections are already set up such that one should be able to configure either driver in the control panel, but upon attempting to configure the necessary driver using the 32 bit version, the very same error occurs. Only the 64 bit tests correctly, but the software uses 32. My question is then - is it possible to point to odbcad32 in the system32 folder, and otherwise, why is it that the ODBC driver might be returning the architecture mismatch error upon configuration in the control panel? I also tried rebuilding the tool source (VB .NET if that matters) in x64 but this did nothing.


Convert GPS to cartesian and scale to a bounding box

I've little knowledge about geographic systems and i'm gonna try to explain my problem the best i can. I've gps coordinates captured from a mobile device WGS84 store in a file, each file stores a path the user did while capturing the coordinates. What i want, is to get the "shape" of that path, i don´t care in this use case, about the geographical location of the path. I need to draw the "shape" of the path on a screen with a specific width and height and it doesn't need to be very accurate. I was thinking about converting the lat/long to ECEF coordinates that gives me x/y/z and i know how to it, but know i need to translate the coordinates to screen coordinates where the origin(0,0) is the top left corner but i don´t know how to do that. Can you give some help on this ? Is there a better way to achieve this ? Thanks.


Is there a On Map Move event in OpenLayers?

In my OpenLayers web app I have registered an `onmousemove` event which performs certain functions. What I am really after though is a way to register and handle the user panning or zooming the map. Thanks in advance


Lidar and Opengeo suite

Opengeo suite has provided a tutorial on how to manage LiDAR data using postgresql. I was trying to follow the tutorial in this link (http://workshops.boundlessgeo.com/tutorial-lidar/). I am new to linux and using ubuntu 12.04. I have successfully installed Opengeo suite 4.1. I started off with the tutorial smoothly , installed point cloud extension and pdal and downloaded the example data. I only got stuck at the section "Loading LiDAR into the database". In this section, I followed the tutorial and implemented following codes successfully, `CREATE EXTENSION postgis; CREATE EXTENSION pointcloud; CREATE EXTENSION pointcloud_postgis;` And also the following pdal commands run smoothly , `pdal info --input 20090429_42122c8225_ld_p23.laz --schema` `pdal info --input 20090429_42122c8225_ld_p23.laz --metadata --xml` Now I faced the issue while implementing the following command,               pdal pipeline laz2pg.xml      From the terminal I got the error,               mzaman@ubuntu:~/Downloads$ pdal pipeline laz2pg.xml      `Caught PDAL exception: unable to connect to database` This is to mention that, I have created a database in Postgresql called 'LiDAR' using pgAdmin. I kept another terminal open and run the following command using users 'postgres' and 'mzaman' mzaman@ubuntu:~/Downloads$ sudo -u mzaman psql LiDAR psql (9.3.4) Type "help" for help. LiDAR=# I also tried with the user 'postgres',               entermzaman@ubuntu:~/Downloads$ sudo -u postgres psql LiDAR      But the first terminal where I run pdal command, still shows the error `Caught PDAL exception: unable to connect to database. My quetion is, how to bypass this error `Caught PDAL exception: unable to connect to database`? Thanks `


ArcGIS 10 Time series to kml

Hey all I am having trouble with the KML export. My time series will not work properly in google earth, it has time start in 0001 when it should start in the year 2001. My incidents time is coded `YYYYMMDDHHSS` do you think if I cut it back to `YYYYMMDD` that it would work? It is a geocoded series of terrorist attacks in Afghanistan. Works fine in arcmap 10 and arcglobe 10.


Loading an older google satellite image with OpenLayers plugin in QGIS

I am using the OpenLayers plugin in QGIS. I want to use the google satellite layer, but there is a cloud on the specific image. Is it possible to load one of the older google satellite images? Like loading the older satellite images on google earth.


Where to download QGIS win7 x64 installer?

When I try to print big map (1200mm x 1200mm x 300dpi) on pretty strong computer (i7 975 12GB RAM) withj Windows 7 I obtain error Error - 'Image too big. Creation of image with bla bla pixels failed' - same on QGIS 1.7.4 and 1.8.0 versions. Can anybody link x64 build of QGIS please - I'd like to check printing on that version? I don't wanna install programming tools and see same error after...


What is difference between Map event onExtentChange and Map.extent property?

I am a newbie and I am looking to persist the map extent across sessions so that it can be loaded instead of the default extent. I believe what I need to do is to store the extent in a cookie after every map extent change. Then when the user reloads the page, set the map extent to the extent last saved. I came across the 'onExtentChange' event as well as the read only property map.extent and map.SetExtent. What is the difference between using the event or using the map.Extent method? Can someone please explain?


Is there any automated process to download POI from Wikimapia?

Is there any automated process to download point of intrest from Wikimapia. I want to download some features such as Schools


How to edit feature in sqlserver geodatabase

I am using arcmap 10.2.1 and arcobject 10.2.1,I just want to conncet to a geodatabase and edit,add,delete a feature from a feature class using arcobject,I tried to connect to database and it works fine but I cant start edit and add a feature to it,here is my code                 public void open_Workspace()         {                   Type factoryType = Type.GetTypeFromProgID("esriDataSourcesGDB.SqlWorkspaceFactory");         IWorkspaceFactory workspaceFactory = (IWorkspaceFactory)Activator.CreateInstance(factoryType);         IWorkspaceFactory2 workspaceFactory2 = workspaceFactory as IWorkspaceFactory2;         // Build a connection string.         String[] connectionProps =      {         "dbclient=SQLServer", "serverinstance=CRAZYVIRUS-HP",              "database=land_gis", "authentication_mode=OSA"      };         String connString = String.Join(";", connectionProps);              IWorkspace workspace = workspaceFactory2.OpenFromString(connString, 0);         ESRI.ArcGIS.Geodatabase.IFeatureWorkspace featureWorkspace = (ESRI.ArcGIS.Geodatabase.IFeatureWorkspace)workspace; // Explict Cast         ESRI.ArcGIS.Geodatabase.IFeatureClass featureClass = featureWorkspace.OpenFeatureClass("DBO.test");              ESRI.ArcGIS.Carto.IFeatureLayer featureLayer = new ESRI.ArcGIS.Carto.FeatureLayerClass();                   featureLayer.FeatureClass = featureClass;                        featureLayer.Name = featureClass.AliasName;         featureLayer.Visible = true;                        IWorkspaceEdit workspaceEdit = (IWorkspaceEdit)workspace;         IWorkspaceEdit2 workspaceEdit2 = (IWorkspaceEdit2)workspaceEdit;         IMultiuserWorkspaceEdit muWorkspaceEdit = (IMultiuserWorkspaceEdit)workspace;              Make sure that non-versioned editing is supported. If not, throw an exception.        if (!muWorkspaceEdit.SupportsMultiuserEditSessionMode            (esriMultiuserEditSessionMode.esriMESMNonVersioned))        {            throw new ArgumentException(                "The workspace does not support non-versioned editing.");        }              // Start a non-versioned edit session.        muWorkspaceEdit.StartMultiuserEditing        (esriMultiuserEditSessionMode.esriMESMNonVersioned);              featureClass.CreateFeature();                   // Stop the edit session. The saveEdits parameter indicates the edit session         // will be committed.         workspaceEdit.StopEditing(true);          }      when I run above code I got this error `Unable to cast COM object of type 'System.__ComObject' to interface type 'ESRI.ArcGIS.Geodatabase.IWorkspaceEdit2'. This operation failed because the QueryInterface call on the COM component for the interface with IID '{199A1F9D-435A-4118-9B3F-4E0B8F984AA0}' failed due to the following error: No such interface supported (Exception from HRESULT: 0x80004002 (E_NOINTERFACE)).` I tried to enable versionning or disable it but there isnot any chance to make it work..Can you help me find the soluton into this error? Is there any thing wrong with my connection?


Combine aerial imagery pictures to create a layer

We are looking for software or a web service that we can put in pictures we took using our weather balloon to create a single image, that we can use as a layer for our GIS project. Preferably free, or free for education. EDIT: We use a ESRI system


How to get a legend for a WMS layer and search in it?

How to bring legend and search for WMS layer in QGis? I can see the WMS layer in Qgis but can't make any legend for identifying feature types and searching in WMS layer?


Bubble Plots in R not displaying

I was attempting to follow a tutorial to create a bubble plot of school district proficiency rates. ## Data shapefile of school districts (Census2010.shp) https://drive.google.com/folderview?id=0ByHQd1ckVHrSYVB4QmNtUWktTGM&usp=sharing District level data: https://docs.google.com/file/d/0ByHQd1ckVHrSdkdlbUsyREtrZms/edit?usp=sharing ## Code reading in data               dist <- read.csv('/Users/Path/to/District.csv', header = T)     state <- readShapePoly('/Users/Path/to/Census2010')      merge shapefile and district data               GA <- state     GA@data <-         merge(GA@data, dist, by.x="NCESID", by.y="NCESID", all.x=T, sort=F )      creating bubble plot               plotvar <- GA@data$CCRPI     nclr <- 8     plotclr <- brewer.pal(nclr,"PuRd")     max.symbol.size=12     min.symbol.size=1     class <- classIntervals(plotvar, nclr, style="quantile") #Some school districts did not report data     colcode <- findColours(class, plotclr)     symbol.size <- ((plotvar-min(plotvar))/         (max(plotvar)-min(plotvar))*(max.symbol.size - min.symbol.size) + min.symbol.size)     plot(GA)     GA.cntr <- coordinates(GA)     points(GA.cntr, pch=16, col = colcode, cex = symbol.size)     points(GA.cntr, cex = symbol.size)      Though the resulting map is just hollow school district outlines. What am I missing here?


Which labelling tool/engine in QGIS will be supported in the long term?

Why does QGIS appear to have two labelling tools? (Note this is not new and old symbology, but that...) there are two different GUIs, one in the toolbar (which gives the Data Defined Settings option) and then in the Layer Properties a tab which leads to Advanced | Date Defined Placement + others. I can turn both on and have two indenpendent lables. I ask, because what will be supported in the long term? ciao


Satellite imagery options for small/medium sized countries?

I need to create land cover classifications for a couple countries in central/eastern Asia. The largest country is 50,000km^2 and the smallest is 10,000km^2. What are my options for recent imagery? At this stage resolution is not my primary concern, so Landsat would be acceptable, but the scan line errors in Landsat 7 prevent me from collecting useful imagery from 2003 onward. Likewise, I am willing to pay for imagery, but very high resolution data may require too many images to feasibly process. My goal is to create a land cover classification that can be compared with a previous classification (~1990) created from a composite of Landsat imagery and aerial photos. **EDIT:** It appears that Landsat 8 data is now available through Glovis - imagery goes back to about April 2013. Finding cloud-free imagery is a bit difficult given how few images there are.


Can GeoTools parse GML CRS definitions?

A kind of branching of this question. Soon I will need to handle **GML** _Coordinate Reference Sysyem_ definitions in a Java servlet. Links:   * **GML OGC** reference.    * Concrete example. I know that GeoTools library can parse GML, and I can see in http://docs.geotools.org/latest/userguide/library/xml/geometry.html some good documentation regarding `Geometry` features. I'd like to know whether I will be able to parse such definitions, and maybe translate them to WKT. Thanks.


Duplicate items in the legend

I have recently created a web map that uses the legend snippet from the ESRI JavaScript API website. https://developers.arcgis.com/en/jav...et_legend.html It works great for the most part, but has one issue, there are duplicate entries in the legend. The rest services that I am using are all symbolized properly without duplicates. The method I'm using to add the services to the map is through a feature layer. ![enter image description here](http://i.stack.imgur.com/xJSr5.png) This is the exact code in my .js file.               dojo.connect(myMap, 'onLayersAddResult', function (results) {         var layerInfo = dojo.map(results, function (layer, index) {             return { layer: layer.layer, title: layer.layer.name };         });              if (layerInfo.length > 0) {             var legendDijit = new esri.dijit.Legend({                 map: myMap,                 layerInfos: layerInfo},                 "legendDiv");                 legendDijit.startup();             }         });          //Layers are added to map     myMap.addLayers([featureLayer1, DMSL_StudyArea]);      Any information would help. Thanks


How to create proper PostGIS points out of GRIB2 file?

I'm trying to load NCEP NOAA weather data from GRIB2 format into PostGIS and then visualize it using QGIS etc. I'm using wgrib2 library to convert the GRIB2 data to csv first.               wgrib2 -match "(:TMP:surface:anl:|:TMP:2 m above ground:anl:)" gfs.t00z.sfluxgrbf00.grib2 -csv fcast.csv      Then I copy CSV to PostGIS and create a geometry column:               COPY tmp (begindate, enddate, name, name2, value, lat, lon)      FROM     '/Users/nl/Work/ncep/fcast.csv'     WITH DELIMITER ',' CSV;     SELECT addGeometryColumn ('public', 'tmp', 'the_geom', 4326, 'POINT', 2);     UPDATE tmp SET lon = lon-360 WHERE lon > 180;     UPDATE tmp SET the_geom=ST_SetSRID(ST_MakePoint(lat,lon), 4326);      When looking at the data in QGIS, I see a strange result: ![enter image description here](http://i.stack.imgur.com/bkVja.jpg) When I expect to see the regular grid - something like this: ![enter image description here](http://i.stack.imgur.com/UKcYC.jpg) I see, that someone has done it already here: How to process GRIB2 Data for Visualization? Can anyone give some details?


Distance of points in polygon to nearest polygon edge

I have a polygon shapefile with about a dozen polygons (Area.shp). I have an equation that calculates a value for points (or small areas) as a function of distance to the edge of these polygons. I would ultimately like to apply this function to all points within each of these polygons and sum them (within each polygon) to get a single value for the polygon. It seems like one of those things that shouldn't be too difficult, but I can't seem to get a working solution. I had an idea of converting the shapefile to raster then raster to points then taking the distance of each point to the nearest polygon edge, then apply the equation. It seemed like a reasonable, if clunk solution but unfortunately, I get all sorts of errors during various conversions using ArcMap 10. Does anyone have any suggestions on how to do this? I am happy to use ArcGIS, QGIS, GRASS, Python, or export to R or any other solution. I am a novice in all of these systems (except R) but am willing to use whatever works. [UPDATE] Additional Information: When I was thinking of points, it was really as a way to approximate numerical integration. I have an equation where the number of animals within a particular patch of habitat (polygon) is a function of distance to the edge of the habitat (nearest edge of the polygon). The equation is a 4pt logistic function               N(i) = C + A/(1 + e^-(x-D)/B)      where all values are known (previously estimated) constants except x which is distance to the edge. N(i) is the number of animals in a small area, roughly 2 m^2. I'm interested in how the size and shape of habitat patches (polygons) affects the total N (sum of N(i) within a polygon) in a patch. I have a variety of patches identified from GIS layers as a single shape file with about a dozen polygons. For the purpose of this study, I can assume that no polygons overlap.


How to make Properties > Change Symbol > Symbol Selector > Edit Symbol dialog loading faster?

Is there any way to tweak ArcMap to make the loading of the Symbol Selector and Symbol Editor dialogs faster? I have a lot of symbols in my style files and I want to keep it that way since it is essential for my workflow. I understand that style files are access mdb databases so i compacted them but that didn't help although it did reduce the size of the databases. There is a lot of lag (5 - 6 seconds) between clicking the "Change Symbol" in the "Properties dialog" and another 5 - 6 seconds between clicking the "Edit Symbol" button in the "Symbol Selector" and the Editor dialog coming up. Another place where things get quite slow is editing custom symbols, especially complex symbols composed of many layers - the properties in the Symbol Property editor take a while to refresh and it refreshes after every single change. I am planning on going through my symbols library and purge any that are unused but i have doubts it will help much. This is a common issue even in with a new install of ArcGIS with only the dafault style file. Does anyone know a way to fix this? ![enter image description here](http://i.stack.imgur.com/JHWsG.jpg) ![enter image description here](http://i.stack.imgur.com/5Mare.jpg)


Trigger to fill attributes (data, client ip) in table for each new objects, PostGIS

In PostgreSQL/PostGIS database stores shapefiles. Client users software - QGIS. Please help me write PostgreSQL/PostGIS trigger, which for each newly added/updated feature in layer added 'time create' and 'client ip address' in attribute table. Attribute table contain next fields: 'gid' - primary key, 'updated' - data create/update feature, 'nameuser' - client ip address. I write next trigger:               CREATE OR REPLACE FUNCTION process_audit() RETURNS TRIGGER AS $audit_roads$     BEGIN          IF (TG_OP = 'INSERT') THEN      INSERT INTO roads SELECT now(), inet_client_addr(), NEW.*;      RETURN NEW;           ELSIF (TG_OP = 'UPDATE') THEN      INSERT INTO roads SELECT now(), inet_client_addr(), NEW.*;      RETURN NEW;           END IF;           RETURN NULL;       END;     $audit_roads$ LANGUAGE plpgsql;               CREATE TRIGGER audit_roads     AFTER INSERT OR UPDATE OR DELETE ON roads FOR EACH ROW EXECUTE PROCEDURE process_audit();      but when I add new feature in QGIS and save them, I receive message: "INSERT has more expressions than target columns" Tracking Edit History using Triggers \- there is example of Tracking Edit History using Triggers, but history saved in nyc_streets_history table (not in nyc_streets) Thanks!


Comparing ShapeFile Data in C#

I have made a List containing all the original data (.shp files), and I have another List containing the modified data. I basically want to check all the modified file data against the original... So here lies the problem, I next a loop within a loop to check if the original exists within the modified List... and it always seems to find a match? Is my approach right? This is the sample code I have for PointZ type shapes                foreach (PointZ origrecord in origPointZList)             {                 Boolean foundIt = false;                      foreach (PointZ modrecord in origPointZList)                 {                     if ((origrecord.X == modrecord.X)                         && (origrecord.Y == modrecord.Y)                         && (origrecord.Z == modrecord.Z)                         && (origrecord.M == modrecord.M))                     {                         foundIt = true;                         break;                          }                     else                     {                         foundIt = false;                     }                      }                 if (!foundIt)                 {                     result += "Record not found";                 }


Start postgresql+postgis without installation on Windows

I have used Windows7x64 and want start PostgreSQL+POSTGIS without installation process. I copy PostgresSQl\9.1 folder from other laptop to c:\opt. I need create .bat file with enviroment variables and add pg_ctl.exe with any parameters. is it possible start Postgresql server without creating Windows services and using any regedit manipulations? Sorry if my question is not correct compiled


Running Erdas command in MS-DOS

I have generated the following command line in ERDAS, and want to run it in MS-DOS. Getting an error- " ** _system cannot find the path specified"_**. I have checked the path by opening it in windows explorer, it navigates to the corresponding file, have tried multiple options, by changing the location or going through .bat approach, still it won't work! C:/ProgramFiles/Intergraph/ERDAS IMAGINE 2014/bin/x64URelease/modeler.exe @"C:/Users/WRD/AppData/Local/Temp/ImagineCommandLine_modeler_3780003780"


clip polygon with the extents of contour polylines

I have a polygon shape that represents the boundary of Fennoscandinavia and a polyline shapefile representing two contour intervals for a trait. I would like to create a new polygon representing the portion of the Fennoscandinavia boundary polygon between the two polylines using ArcMap 10. There does not appear to be any tools in the 'data management', 'conversion' or 'Spatial Analyst' toolboxes that work. Any suggestions? ![Fennoscandinavia with contour intervals for trait frequency](http://i.stack.imgur.com/tUwXS.png)


ArcGIS Server 10.1 and SQL Server 2012 SP1

We have just built and new ArcGIS/ArcSDE system with SQL Server 2012 SP1. We are experiencing very long periods (up to 2 minutes) just to connect to the database in Desktop. Desktop does connect and we can publish map services, but everything is painfully slow. We have ruled out network issues, and wonder if we are missing something. Has anyone experienced this?


Python console not working at all

My python console won't work at all. I can't open it. My QGIS version is 2.4. OS is win7 64. I've reinstalled QGIS a few times but I always get the same result. Strange. Something went wrong with the source code I guess? Here's the error :               Traceback (most recent call last):   File "", line 2, in    File "C:/OSGEO4~1/apps/qgis/./python\console\console.py", line 43, in      show_console          _console = PythonConsole( parent )   File "C:/OSGEO4~1/apps/qgis/./python\console\console.py", line 75, in      __init__          self.console = PythonConsoleWidget(self)   File "C:/OSGEO4~1/apps/qgis/./python\console\console.py", line 101, in      __init__          self.shellOut = ShellOutputScintilla(self)   File "C:/OSGEO4~1/apps/qgis/./python\console\console_output.py", line 103,      in __init__          self.insertInitText()   File "C:/OSGEO4~1/apps/qgis/./python\console\console_output.py", line 145,      in insertInitText          socket.gethostname()) UnicodeDecodeError: 'ascii' codec can't decode byte 0xe9 in position 9: ordinal not in range(128) Version de      Python :      2.7.5 (default, May 15 2013, 22:44:16) [MSC v.1500 64 bit (AMD64)] Version de QGIS :      2.4.0-Chugiak 'Chugiak', 8fdd08a Chemin vers Python : ['C:/OSGEO4~1/apps/qgis/./python/plugins\\processing',


Preserving formatting of Scales & Legends in Templates

My QGIS templates don't seem to preserve the formatting of my scales or legends. Is there a way to save these so they don't have to be redone each time? (I'm doing a series of maps that will be the same scale and have the same legends.) Thanks TVZ


Update fields automatically after adding new geometry

I have two layers, one for "States" (Polygon) and the other for "Important_Cities_location" (Point). "States" feature class is completed with all attributes. I want to update the state name field automatically in "Important_Cities_Location" when I place them in the respective "States" spatially. Can anyone help?


How can I determine the layer name in an edit session?

I have a problem with an editor extension in an Add-In for ArcGIS 10. The editor extension listens on different events occuring on shapefile layers via registering listeners:               mEditEvents.OnChangeFeature += OnChangeFeature;     mEditEvents.OnCreateFeature += OnCreateFeature;     mEditEvents.OnDeleteFeature += OnDeleteFeature;      All these callbacks get are a parameter of type IObject, which can be cast to IFeature. With this IFeature I can get the feature class and its alias name (IFeature.Class.AliasName) which seems to be the name of the shapefile associated with this feature, but I want the name of the layer associated with this feature which can be different from the shapefile's name. Is there any way how to get this layer name? Thanks in advance! Update: I found a solution (see my comment below) - thanks for the help!


Calculating Speed of Android device in Navigation mode using ArcGIS Runtime SDK?

I am working on a Navigation App and I am using ArcGIS runtime sdk for this. It looks pretty straight forward and I was able to create a simple app. Now I want to calculate the speed of the device while navigating. I was wondering if that is possible or not as I did not find anything in this regard in ArcGIS documentation.


Symbol Level in ArcGIS 9: 'Broken' look of symbols, how to resolve

I am using symbol levels in ArcMap to style roads so that I have a banded feel. My problem is that at intersections and sharp turns, the cartography gives a broken effect in the road. How can this be resolved. I have tweaked and searched without success. . ![Pic](http://i.stack.imgur.com/xYZiP.jpg) Unwanted effect pointers in blue. UPDATE EDIT: After further investigation I had an effect closer to what I wanted but with some side effect. From left to right; First image shows (i)the settings in ArcMap (Symbol Property Editor) , (ii) Resulting effect solving the 'craked effect' but undesirable capping at the end and (iii) the desired effect at the end of a road. .![enter image description here](http://i.stack.imgur.com/L6wMS.jpg)


How can we display an animated output for watershed management in openlayers

Here i'm using Tomcat(7.0)+GeoServer(2.1.3)+Openlayers(2.11). How can we display an output like Flow of water for the **forest watershed management** application. I have to show the **water flow** Based on the slope range range of the area, rain fall(qty in mm), duration of the rain, and so on. For this right now I have only 2 methods to generate this output: 1)1 is to fill the color of layer by percentage wise. 2)2nd is to render the layer by percentage wise. can Anyone guide me how it is possible in Openlayers? All help is greatly appreciated... Best Regards, _**RamRamesh_**..:)


Raster to Vector seems not to work

> **Possible Duplicate:**   >  Qgis polygonize problem - DLL load fail When I'm trying to invoke the Polygonize Module in the Raster Tab I get the following error message (I did give bothe the source raster and the polygone output layer): Traceback (most recent call last): File "C:\PROGRA~1\Quantum GIS Lisboa\bin\gdal_polygonize.py", line 34, in import gdal, ogr, osr File "C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27\lib\site-packages\gdal.py", line 2, in from osgeo.gdal import deprecation_warn File "C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27\lib\site-packages\osgeo__init__.py", line 21, in _gdal = swig_import_helper() File "C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27\lib\site-packages\osgeo__init__.py", line 17, in swig_import_helper _mod = imp.load_module('_gdal', fp, pathname, description) ImportError: DLL load failed: Die angegebene Prozedur wurde nicht gefunden. Could anybody give me a hint on what's the problem, thx, Ludwig


How can I fix DLL loadfailed with Polygonize in QGIS 1.8?

> **Possible Duplicate:**   >  Qgis polygonize problem - DLL load fail I'm getting the following error message when I try to use the polygonize tool in QGIS 1.8 (windows 7):               Traceback (most recent call last):      File "C:\PROGRA~2\Quantum GIS Lisboa\bin\gdal_polygonize.py", line 34, in      import gdal, ogr, osr      File "C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages\gdal.py", line 2, in      from osgeo.gdal import deprecation_warn      File "C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages\osgeo\__init__.py", line 21, in      _gdal = swig_import_helper()      File "C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages\osgeo\__init__.py", line 17, in swig_import_helper      _mod = imp.load_module('_gdal', fp, pathname, description)      ImportError: DLL load failed: No se encontró el proceso especificado.       I can perform this operation with the same raster file on QGIS 1.7 on this very same computer. (This is for a class I'm giving and all the students have installed 1.8).


DLL load error with Gdal polygonize in QGIS 1.8 (windows 7)

> **Possible Duplicate:**   >  Qgis polygonize problem - DLL load fail [I had already posted this question, but was closed and redirected to a supposedly duplicate question. I'm re-posting since I believe this is note exactly the same case.] I'm getting the following error message when I try to use the polygonize tool in QGIS 1.8 (windows 7): Traceback (most recent call last): File "C:\PROGRA~2\Quantum GIS Lisboa\bin\gdal_polygonize.py", line 34, in import gdal, ogr, osr File "C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages\gdal.py", line 2, in from osgeo.gdal import deprecation_warn File "C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages\osgeo__init__.py", line 21, in _gdal = swig_import_helper() File "C:\PROGRA~2\Quantum GIS Lisboa\apps\Python27\lib\site-packages\osgeo__init__.py", line 17, in swig_import_helper _mod = imp.load_module('_gdal', fp, pathname, description) ImportError: DLL load failed: No se encontró el proceso especificado. I have already tried or checked the following with no success:   * Making sure there is no PYTHONPATH variable   * adding "C:\PROGRA~2\Quantum GIS Lisboa\bin\", "C:\PROGRA~2\Quantum GIS Lisboa\bin" and "C:\PROGRA~2\Quantum GIS Lisboa\" to the beginning of PATH variable   * Checking that GDAL was not installed separately   * Installed Visual Basic 2010 Express (couldn't find Visual Studio 2010 as suggested in the other post)   * "Path to the GDAL executables" and "Path to the GDAL Python modules" are blank in he GDAL setting menu within QGIS. Any ideas what this could be? I wonder why would it be the path variable if all other Gdal are working fine. Wouldn't it have a DLL load error for all or other Gdal modules?


Polygonize Raster to Vector Error

I have a problem when I am trying to Polygonize (Raster To Vector). I have a Corine map with 44 different land uses. I run a query and I save the selected areas with Bands 0 and 1. Now I am trying to Polygonize the map and to save the selected areas as vector data. It gives me the following ERROR:               Traceback (most recent call last):      File "C:\PROGRA~1\Quantum GIS Lisboa\bin\gdal_polygonize.py", line 34, in      import gdal, ogr, osr      File "C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27\lib\site-packages\gdal.py", line 2, in      from osgeo.gdal import deprecation_warn      File "C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27\lib\site-packages\osgeo\__init__.py", line 21, in      _gdal = swig_import_helper()      File "C:\PROGRA~1\Quantum GIS Lisboa\apps\Python27\lib\site-packages\osgeo\__init__.py", line 17, in swig_import_helper      _mod = imp.load_module('_gdal', fp, pathname, description)      ImportError: DLL load failed: The specified procedure could not be found.


How to add a google/OSM/bing etc style map layer to Quantum GIS

I am just getting started using QGIS. I am aiming to create a heat map of Census tracts in Georgia based on population. I have all of the data - the population for each census tract in Georgia. I've been able to create the heat map successfully by importing a census tract shape file and following this tutorial: http://qgis.spatialthoughts.com/2012/03/using-tabular-data-in-qgis.html However I am stuck on how to add a nice looking map layer underneath which will show common geographical features - highways, cities, water courses etc. I have tried importing OSM data from cloudmade but there is far too much detail. The map is supposed to be fairly high level but the shape files from there show every single road. I am trying to replace MS MapPoint (no 2010 census tracts). The functionality I am trying to get would allow for different levels of detail depending on how far I zoom. Generally just showing major roads / cities etc. I don't mind if this is automatic or manual. Could anyone point me in the right direction? Any help would be much appreciated. Thanks in advance.


How to convert Lat & Long from Onemap.sg to Google Maps

I have a set of address for Singapore, I would like to take the Lat & Long from onemap.sg (SYV21: X, Y / Easting, Northing) and then store it in my databse , convert into (WGS84: Longitude, Latitude) used by google maps and then display on Google Maps. I have looked at this -- http://sgonemap.codeplex.com/ but could not get much out of it. It will be really nice if someone can help me on this.


How to find intersection points between two circles?

We have two points (centers of two circles) and their radius in meters, those radius make the circle. We need to find intersection points. For example we have lat1 = 55.685025, lng1 = 21.118995, r1 = 150 and lat2 = 55.682393, lng2 = 21.121387, r2 = 250. Below you can find our current formula:               // Find a and h.     $a = ($circle_1_r * $circle_1_r - $circle_2_r * $circle_2_r + $distance * $distance) / (2 * $distance);     $h = sqrt($circle_1_r * $circle_1_r - $a * $a);          // Find P2.     $circle_3_x = $circle_1_x + $a * ($circle_2_x - $circle_1_x) / $distance;     $circle_3_y = $circle_1_y + $a * ($circle_2_y - $circle_1_y) / $distance;          // Get the points P3.     $intersection_1 = $this->newLatLngPoint(         ($circle_3_x + $h * ($circle_2_y - $circle_1_y) / $distance),         ($circle_3_y - $h * ($circle_2_x - $circle_1_x) / $distance)     );          $intersection_2 = $this->newLatLngPoint(         ($circle_3_x - $h * ($circle_2_y - $circle_1_y) / $distance),         ($circle_3_y + $h * ($circle_2_x - $circle_1_x) / $distance)     );      We find such intersection points (yellow markers), however those locations doesn't match in real world. ![enter image description here](http://i.stack.imgur.com/4VUif.png) Someone, can help to find where the problem is and how to sort it ? P.S. Does the altitude (Height above mean sea level) affect the final result? I don't use it but may be should?


Problem with public key, when i use mapquest directions service?

I recently registered in mapquest service and get api key, then i tried to use make request to directions service, but get following response:               Ext.data.JsonP.callback12({         "route": {             "routeError": {                 "message": "",                 "errorCode": 3             }         },         "info": {             "copyright": {                 "text": "© 2013 MapQuest, Inc.",                 "imageUrl": "http://api.mqcdn.com/res/mqlogo.gif",                 "imageAltText": "© 2013 MapQuest, Inc."             },             "statuscode": 403,             "messages": ["This is not a valid key. Please check that you have entered this correctly. If you do not have a key, you can obtain a free key by registering at http://developer.mapquest.com."]         }     });


Can a new point automatically be populated with data from the containing polygon

I need some help to create a script. Imagine that I'm using 3 shapes: 1 point and 2 polygons. When I create a new point I want it to assume some attributes of the polygons in which that point is inserted. Any idea?


How to set my own CRS in QGIS

I want to transform WGS to MRT (GCS for Malaysia) using QGIS. Is "save as" to a new CRS the correct method? With ArcGIS, I can create transformation method to do transformation but I'm not sure in QGIS. If it is possible, how? And what is the string for xyz rotation, xyz translation, scale factor?


Problems in gdalwarp

I was trying to use gdalwarp inside a qgis 2.0 model via processing but I think that the implementation is not working properly. So I would like to know if anyone has encountered this problem? ![enter image description here](http://i.stack.imgur.com/whjU8.png)


How to change windows can't find D:\Program

Today, I setup QGIS. When I run msys, it appear 'windows can't find D:\Program'. Why? Thank you!


UnboundLocalError: local variable referenced before assignment

I have following simple function to get percent values for different cover types from a raster. It gives me following error: UnboundLocalError: local variable 'a' referenced before assignment which isn't clear to me. Any suggestions, please.               def percent(ras):         temp = arcpy.SearchCursor(ras,"","","Value;Count","")         for row in temp:             theval = row.getValue("Value")             thecount = row.getValue("Count")             if row.getValue("Value") == 1:                 a = thecount             if row.getValue("Value") == 2:                 b = thecount             if row.getValue("Value") == 3:                 c = thecount             if row.getValue("Value") == 4:                 d = thecount         return a,b,c,d


How to extrude polygon features from a TIN surface in ArcScene?

How do I Extrude polygon feature class from a TIN surface using ArcScene 10.1?


I'm using osm2pgsql to import an extract of OpenStreetMap's data so I can geocode an address to a lat/lng. But I'm new to OSM and mapping so I'm having trouble figuring out what data is relevant and what is unnecessary. What is the bare minimum data I would need to import from an OSM data file in order to geocode an address? Thanks!


Can I Combine Areas in Quantum GIS?

I am working on a provincial annual health bulletin. I have a .shp file that has GIS data about the smallest divisions called "GN Area". Two or three GN Areas form one MOH area. I want to get the MOH area details displayed in the map. How can I do it? I am using QGIS. Is there any way where GN Area data exported to excel and add the relevant MOH area and import it back to Quantum GIS and display MOH Areas rather than GN ares? If that is not possible with Quantum GIS, is there any other method to do it?


ArcGIS JavaScript API layerDrawingOptions seem restricted to 1 option, and it doesn't have any effect

I have two problems with ArcGIS Server's (10.0) dynamic layers when calling Layer.setLayerDrawingOptions from v3.2 of the JavaScript API. 1) It doesn't work. I set the layerDrawingOptions array for a dynamic layer and the resulting image is rendered the same as if I hadn't changed anything. I have verified that the layer ID is correct, and also tried several approaches to defining the LayerDrawingOptions object. When I inspect the request that goes to my map server it contains the following `dynamicLayers` property:               [{"id":0,"source":{"type":"mapLayer","mapLayerId":0},"definitionExpression":"ASSET_ID IS NULL"},{"id":1,"source":{"type":"mapLayer","mapLayerId":1},"definitionExpression":"ASSET_ID IN (27)","drawingInfo":{"renderer":{"type":"simple","symbol":{"color":[255,0,0,255],"outline":{"color":[0,0,0,255],"width":1,"type":"esriSLS","style":"esriSLSSolid"},"type":"esriSFS","style":"esriSFSSolid"}},"transparency":1,"scaleSymbols":false,"showLabels":false}},{"id":2,"source":{"type":"mapLayer","mapLayerId":2},"definitionExpression":"ASSET_ID IS NULL"},{"id":3,"source":{"type":"mapLayer","mapLayerId":3},"definitionExpression":"ASSET_ID IS NULL"}]      The feature with ASSET_ID = 27 is rendered but in its default style, however `"color":[255,0,0,255]` says it should be solid red. The map service referenced by this layer only contains a single polygon layer and that is ID 1. I can't see any way that I am accidentally attempting to style the wrong layer. 2) If I define multiple layerDrawingOptions in the array, one for each layer within the map service, I am forced to provide a proxy URL (even though the resulting request is nowhere near the 2000 character limit). When what should be the map request goes to the map server via the proxy, with multiple layerDrawingOptions, the map server returns a JSON string that describes an image (height, width, bbox etc.) but has an empty href (""). Although there are no errors in my AGS log I assume this means I'm not allowed to set multiple layerDrawingOptions in a single layer. Is this actually true? So if anyone can provide thoughts or experience on 1) why layerDrawingOptions seem to have no effect, or 2) why I am seemingly not allowed to set more than one option for a map service, I'd really love to hear it/them.


How to get the number of objects at a pixel position?

Is it possible to get the total number of objects (Lines,Points) at a given pixel position in the map?


Arcpy remove invalid character

I've built a model that takes the name of the Excel worksheet containing xy data, and uses that name as the name of the output shapefile. ( Worksheet "massachusetts_sightings" in an xlsx becomes "massachusetts_sightings$.shp"). I do this using Parse Path. Notice that a "$" is added. If I then feed this shapefile, massachusetts_sightings$.shp into a clip tool (say, to clip out points in the data that are not in my study area), I get an error for an invalid character in the clip tool. My suspicion is the $. I know I could remove this with Calculate Value...but I can't figure out the code block. My parse path output is called Value: ![Codeblock](http://i.stack.imgur.com/BFZRy.png) Can someone help troubleshoot? I haven't used Python in years-               # Import arcpy module     import arcpy          # Load required toolboxes     arcpy.ImportToolbox("Model Functions")          # Set Geoprocessing environments     arcpy.env.scratchWorkspace = "X:\\Tools\\XY\\point_upload\\Scratch"     arcpy.env.workspace = "X:\\Tools\\XY\\point_upload\\Data"          # Script arguments     XY_Table = arcpy.GetParameterAsText(0)          Define_Coordinate_System = arcpy.GetParameterAsText(1)          Projection = arcpy.GetParameterAsText(2)          Geographic_Transformation = arcpy.GetParameterAsText(3)          Study_Area = arcpy.GetParameterAsText(4)          Projected_and_Clipped_Data = arcpy.GetParameterAsText(5)     if Projected_and_Clipped_Data == '#' or not Projected_and_Clipped_Data:         Projected_and_Clipped_Data = "X:\\Tools\\XY\\point_upload\\Data\\%fixed%_Clip.shp" # provide a default value if unspecified          # Local variables:     xy_event = XY_Table     v_defined = xy_event     defined_data = v_defined     v_projected = defined_data     Value = XY_Table          # Process: Parse Path     arcpy.ParsePath_mb(XY_Table, "NAME")          # Process: Make XY Event Layer     tempEnvironment0 = gp.scratchWorkspace     arcpy.env.scratchWorkspace = "X:\\Tools\\XY\\point_upload\\Scratch"     tempEnvironment1 = gp.workspace     arcpy.env.workspace = "X:\\Tools\\XY\\point_upload\\Data"     arcpy.MakeXYEventLayer_management(XY_Table, "", "", xy_event, "", "")     arcpy.env.scratchWorkspace = tempEnvironment0     arcpy.env.workspace = tempEnvironment1          # Process: Copy Features     tempEnvironment0 = gp.scratchWorkspace     arcpy.env.scratchWorkspace = "X:\\Tools\\XY\\point_upload\\Scratch"     arcpy.CopyFeatures_management(xy_event, v_defined, "", "0", "0", "0")     arcpy.env.scratchWorkspace = tempEnvironment0          # Process: Define Projection     arcpy.DefineProjection_management(v_defined, Define_Coordinate_System)          # Process: Project     tempEnvironment0 = gp.scratchWorkspace     arcpy.env.scratchWorkspace = "X:\\Tools\\XY\\point_upload\\Scratch"     tempEnvironment1 = gp.workspace     arcpy.env.workspace = "X:\\Tools\\XY\\point_upload\\Data"     arcpy.Project_management(defined_data, v_projected, Projection, Geographic_Transformation, "")     arcpy.env.scratchWorkspace = tempEnvironment0     arcpy.env.workspace = tempEnvironment1          # Process: Clip     arcpy.Clip_analysis(v_projected, Study_Area, Projected_and_Clipped_Data, "")          # Process: Calculate Value     arcpy.CalculateValue_management("\"%Value%\".strip(\"$\")", "\\n", "String")