Adding ArcMap standard toolbar to the form in C# .NET WinForms project

I have Windows Forms project with ArcObjects framework. I have already added axMapControl1 to the form. ![enter image description here](http://i.stack.imgur.com/zbdMe.png) The question is How can I add a standard toolbar to the form and link it with axMapControl1 ? Like this: ![enter image description here](http://i.stack.imgur.com/oyUaK.png) Is there any standard toolbar control in the ArcObjects or should make it by myself?


ArcGIS tool for indicating changes in land use

I have a couple of shapefiles in which there is some information about land use (polygons) such as water, forests, soils, buildings etc. Each of the shapefiles represents a certain year of the same area of interest. I want to indicate the area which was the biggest change over the years. Can you suggest me a procedure to do it in ArcGIS 10 ?


Retrieve WMS service from GeoWebCache with OpenLayers

Ok guys, here we go ! For my needs, I basically added a **standalone version** of **GeoWebCache**. After some configuration, the **GeoWebCache** retrieve the data from **GeoServer** using Capabilities. Now, in my **OpenLayers** code, I would like to question **GeoWebCache** to retrieve my data, and not **GeoServer** (if it is possible). For the moment, I have something like this (working) :                   var fond_carto = new OpenLayers.Layer.WMS(             "Fond cartographique 2",              "http://179.172.277.77:8080/geoserver/wms",             {                 layers: 'departement',                 TILED: 'true'             }      I tried stuff like that (not working) :                   var fond_carto = new OpenLayers.Layer.WMS(             "Fond cartographique 2",              "http://179.172.277.77:8080/geowebcache/sevice/wms",             {                 layers: 'departement',                 TILED: 'true'             }       Is there a way to accomplish that ? After some reading, I assume I'm a bit lost with this standalone version of **GeoWebCache**. It is written it works as a proxy between client ( **OpenLayers** ) and server ( **GeoServer** ). But I read that the TILED: 'true' option means that **GeoServer** will send the layers for tiling to **GeoWebCache**. As I desactivated the integrated version of **GeoWebCache** in **GeoServer** , and as my WMS request are still tiled, I assume this is true. But then it's not really a proxy. Could it works as follow: _openlayers => standalone geowebcache => geoserver => standalone geowebcache => openlayers_ Or it has to work that way : _openlayers => geoserver => standalone geowebcache => geoserver => openlayers_ **Thanks for any tips or any clarifications :-) !** # EDIT After some more configuration in geowebcache.xml, I still have the error but I can add some details to help. This is the URL of the request that is not displaying (pink tiles) :               http://179.172.277.77:8088/geowebcache/sevice/wms?LAYERS=departement&TILED=true&SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&STYLES=&FORMAT=image%2Fjpeg&SRS=EPSG%3A4326&BBOX=2.8125,48.515625,3.515625,49.21875&WIDTH=256&HEIGHT=256:8080/geowebcache/sevice/wms?LAYERS=departement&TILED=true&SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&STYLES=&FORMAT=image%2Fjpeg&SRS=EPSG%3A4326&BBOX=2.8125,48.515625,3.515625,49.21875&WIDTH=256&HEIGHT=256 :      and when I try to copy/paste this url, I got a GWC error :               400: Thread 35 Unknown layer departement. Check the logfiles, it may not have loaded properly.      I don't get why I got this error. Here's my geowebcache.xml description of the layer :               <wmsLayer>       <name>departement</name>       <mimeFormats>         <string>image/gif</string>         <string>image/jpeg</string>         <string>image/png</string>         <string>image/png8</string>       </mimeFormats>       <gridSubsets>         <gridSubset>           <gridSetName>EPSG:4326</gridSetName>         </gridSubset>       </gridSubsets>       <wmsUrl>         <string>http://179.172.277.77:8080/geoserver/wms</string>       </wmsUrl>       <wmsLayers>aaz:departement</wmsLayers>     </wmsLayer>


How to get Database name for layer in ArcObjects?

I am writing an ArcGIS ArcMap Add-In and need to determine what database a given layer from. How can I determine this information? I attempted to use the FeatureLayer to get a FeatureClass, and then get the FeatureDataset from that. From FeatureDataset I can get the name and parse out the first component that seems to be the DB name. However this approach seems really hacky and does not work if the layer is not part of a FeatureDataset. When I look at the layer's properties, it clearly shows the Database (see the highlighted line in the screenshot) so I think I can get the DB from a layer but I don't know what interface to use. ![enter image description here](http://i.stack.imgur.com/bpHL1.gif)


ArcEngine: Drawing graphics returned from BackgroundWorker

My application is drawing graphics objects on the map(as simple as adding IElement to GraphicsContainer). When I create graphics elements and add them to GraphicsContainer in BackgroundWorker DoWork method, ActiveView Refresh method does not show any added graphics. Are there any triks to take into account when drawing on the map in different thread? (if at all possible) Any help will be greatly appreciated.


Measure area of a polygon in specific units using ArcPy

I am trying to find the area of a particular polygon (`OriginalPolygon`, below), and I can do that... But, I'd like it in acres regardless of the original units of the polygon.(Data is coming into this script in a range of projections, but we're always interested in "acreage.") If I was calculating geometry to store in a field, then I could use `feat.area@acres` \-- but I don't need this stored as an attribute, I just need it as a temporary variable to tell the user. This question implies that unless I'm actually calculating geometry, the units of the shape area measurement are tied to the projection and therefore fixed. So: Is there any way to calculate geometry in such a way that the value is only a numerical output that I can display as the script runs, or am I stuck with making an attribute field, calculating geometry, reading the field value to tell the user, and then deleting the field?               rows = arcpy.SearchCursor(OriginalPolygon)     shapeName = arcpy.Describe(OriginalPolygon).shapeFieldName     for row in rows:         feat = row.getValue(shapeName)         PolygonArea = feat.area


How to locate a point on google maps/earth using cadastral maps in Kenya

I'm an architect, so not too well versed with all this datum/transformation talk. My little problem: i want to be able to enter a position using UTM or deg/sec/min coordinates into google earth and fly there. the cadastral maps we have come with eastings and northings and a general reference that reads: coordinates referred to Equator and 37* E (37deg E), all in meters station:A +X (Easting):116379.70 -Y(Northing) : 6309.10 any help out there??


Unmodified OpenGeo Suite Geoexplorer proxy fails with HTTP 500, regardless of servlet

We're running an Opengeo Suite (3.0.1 currently) and an unknown amount of time ago the proxy redirects to other data stores GeoExplorer does* started failing with HTTP 500 (Internal server error). *               var app = new GeoExplorer.Composer({             authStatus: {{status}},             proxy: "../proxy/?url=",      applied to every store with an absolute URL(doesn't matter whether valid or not), like the default               suite: {url: "http://v2.suite.opengeo.org/geoserver/ows      so that the first call becomes `GET /geoexplorer/proxy/?url=http://v2.suite.opengeo.org/geoserver/ows?_dc=4734342332&SERVICE=WMS`. In this environment a newly installed and unmodified Suite does this, and an unmodified GeoExplorer webapp does the same thing in Tomcat. The specific error, same from jetty as Tomcat, is " _Error: No valid JSGI response: [object Object] (ringo/jsgi/connector.js#88)_ ". The environment is a Windows Server 2008. In other environments (on other machines) the Suite has no such issues, so I'm wondering whether it has to do with some network setting or other in Windows? Further, I'm not completely sure what this proxy is for and have been unable to find any documentation for it. Any hints would be very welcome :) Thanks for any hints. Edit. Further research has revealed that this error message (HTTP500, ringo/jsgi/connector.js:88) appears if the machine hosting the servlet has no internet connection. It has been observed on a different machine, running Win 7. I wonder if the servlets in the affected Win 2008 server environment are somehow prevented from reaching the internet, and if this can be detected in any logs / audit files of Win 2008 server R2?


How do I convert MULTIPOLYS into LINESTRINGS in PostGIS

I have a table with a geom column as a geometry(MULTIPOLYGON) I would like to create an additional column of type geometry(LINESTRING) How can I do this using a PostGIS Query? I need to be able to use the linestring in Geoserver to style labels along the line.


Area calculation of an object without coordinates

I have a photo of an object with a scale and I want to calculate its area. I guess I have to create a polygon of the area to be calculated to use the 'field calculator'. But how do I tell the program what a centimeter (cm) is?


Getting compilation error on esriLicenseProductCode.esriLicenseProductCodeArcView

I am getting compilation error on esriLicenseProductCode.esriLicenseProductCodeArcView   1. I have installed Arcgis 10.1 and authorized ArcEditor. I am not able to authorized ArcView.   2. I have not installed Runtime Engine. Am I missing anything any suggestions.


Represent full globe in Oracle Spatial

In SQL Server 2012 I can do this:               -- Construct a new FullGlobe object (a WGS84 ellipsoid)     DECLARE @Earth geography = 'FULLGLOBE'          -- Calculate the area of the earth     SELECT PlanetArea = @Earth.STArea()      Is there a similar construct in Oracle Spatial? Alternatively what would be the best way to represent the entire earth in an 8307/4326 SRID? Is there anything in the WKT standard that speaks to this?


How to install Quantum GIS Lisboa and scipy

I'm trying to install some plugins of the QGIS but an error occurs. The program asks the scipy program. I installed Python, the Numpy and scipy. Also installed the program R. However, QGIS continues asking for scipy. Could you help me?


Convert Gauss-Kruger to UTM

How to covert Gauss-Kruger coordinates to UTM coordinates with GeoTools java library? As I have figured out the location (zone, country) is important. I need the calculation for Slovenia.


Join is deleting fields in arcpy

I'm working on a script to automate the conversion of CAD .dwgs of rooms to feature classes, I'm nearing completion and am now working on the second portion of the script (I've seperated it into two modules as the first part takes about an hour to run and the second part should execute MUCH faster). Part of the script involves adding two fields to one layer, and then joining it to another layer. The problem is, when I execute "arcpy.AddJoin_management" the fields I added aren't in the attribute table. Whats puzzeling is the fields show up when I crate a shapefile from the layer when I perform "arcpy.CopyFeatures_management" before I execute the join.               import arcpy     from arcpy import env     import os     import math          # Set the MTolerance to 0.02     arcpy.env.MTolerance = 0.02                         # Set overwrite option          arcpy.env.overwriteOutput = True          #Set Workspace          arcpy.env.workspace = "in_memory"          #input dataset from part 2     converted_CAD = r"C:\Users\JOE.AWESOME\test_script_output\test.shp"          #create a minimum boundary geometry     arcpy.MinimumBoundingGeometry_management (converted_CAD, 'bounding_geom', 'RECTANGLE_BY_AREA','', '', 'MBG_FIELDS')                         #create layers for an attribute join     arcpy.MakeFeatureLayer_management (converted_CAD, 'CAD_layer')     arcpy.MakeFeatureLayer_management ('bounding_geom', 'bounding_geom_layer')          #Add a field to the table for the difference between the Minimum Bounding Geometry and the CAD_layer area          arcpy.AddField_management ('CAD_layer', 'area_diffe', 'FLOAT')          arcpy.AddField_management ('CAD_layer', 'area_sqM', 'FLOAT')                    #join the two layers for analysis          arcpy.AddJoin_management ('CAD_layer', 'FID', 'bounding_geom_layer', 'ORIG_FID')          #calculate the area of the polygons in the CAD_layer          arcpy.CalculateField_management ('CAD_layer', 'area_sqM', '!shape.area@squaremeters!', 'PYTHON_9.3')          arcpy.CopyFeatures_management ('CAD_layer', r"C:\Users\JOE.AWESOME\test_script_output\join_test.shp")      I've tried changing them to a FLOAT from a DOUBLE, but that didn't change anything. I'm stuck, maybe it's something little. In the words of Michael Bolton, "I always miss some mundane detail."


Exporting selected features to a new shape file

I have a street layer and I only need to export a subset of features to a new shape file where the "TestDone" attribute field is not '0' (Zero). My current code ( which is pasted below ) exports the total street layer. I can use a query filter to get the non zero values out and assign it to an ICursor. But then not sure how to conenter code herevert it to a IfeatureClass. ( Not really sure whether the ICursor and then trying to convert to IFeatureClass is the correct way to go even ) How can I just get features that have a non zero value in the "TestDone" attribute field and then convert that to a feature layer as the second argument in the ExportFeatureClass method ( please see below code) needs a IFeatureClass. Any help with some code guidance will be greatly appreciated. _**_*CODE _*_ ****               private void btnExport_Click(object sender, EventArgs e)     {     SaveFileDialog saveFileDialog1 = new SaveFileDialog();     saveFileDialog1.Filter = "Shape file|*.shp";     saveFileDialog1.DefaultExt = "shp";     saveFileDialog1.AddExtension = true;     saveFileDialog1.InitialDirectory = @"C:/";     saveFileDialog1.Title = "Check data";     saveFileDialog1.ShowDialog();     if (saveFileDialog1.FileName != "")     {     string strFile = System .IO.Path.GetFileName ( saveFileDialog1.FileName) ;     IMap pMap = default(IMap);     IWorkspaceFactory pWorkspaceFactory = default(IWorkspaceFactory);     IWorkspace pWorkspace = default(IWorkspace);     IFileName pFileName = default(IFileName);     IFeatureLayer2 pFeatureLayer = default(IFeatureLayer2);     IFeatureClass pFeatureClass = default(IFeatureClass);     IPropertySet pPropertySet = default(IPropertySet);     IGeometryDef pGeometryDef = default(IGeometryDef);     IGeometryDefEdit pGeometryDefEdit = default(IGeometryDefEdit);     string deletefile = null;          pMap = clsFunctions.pMxDocument.FocusMap;     pFileName = new FileNameClass();     pFileName.Path = System .IO.Path.GetDirectoryName (saveFileDialog1.FileName);     pPropertySet = new PropertySet();     pPropertySet.SetProperty("DATABASE", pFileName.Path);     pWorkspaceFactory = new ESRI.ArcGIS.DataSourcesFile.ShapefileWorkspaceFact ory();     pWorkspace = pWorkspaceFactory.Open(pPropertySet, 0);          IGeographicCoordinateSystem pGCS = default(IGeographicCoordinateSystem);     SpatialReferenceEnvironment pSpatialReferenceEnv = default(SpatialReferenceEnvironment);     ISpatialReference2 pSpatialReference = default(ISpatialReference2);     pSpatialReferenceEnv = new SpatialReferenceEnvironment();     pGCS = pSpatialReferenceEnv.CreateGeographicCoordinateSys tem(Convert.ToInt16(esriSRGeoCSType.esriSRGeoCS_WG S1984));     pSpatialReference = (ISpatialReference2)pGCS;          ILayer pLayer;     int intI = 0;              for (intI = 0; intI <= pMap.LayerCount - 1; intI++)         {         pLayer = pMap.get_Layer(intI);         if (pLayer.Name == cmbStreets.Text)         {         pFeatureLayer = (IFeatureLayer2)pLayer;         }         }          pFeatureClass = pFeatureLayer.FeatureClass;     pGeometryDef = new GeometryDef();     pGeometryDefEdit = (IGeometryDefEdit)pGeometryDef;     pGeometryDefEdit.GeometryType_2 = pFeatureClass.ShapeType;     pGeometryDefEdit.SpatialReference_2 = pSpatialReference;          ExportFeatureClass(strFile, pFeatureClass, pWorkspace, pSpatialReference, pGeometryDef);     }                    public void ExportFeatureClass(string strLayerName, IFeatureClass pFeatureClass, IWorkspace pOutWorkspace, ISpatialReference pSpatialReference, IGeometryDef pGeometryDef)     {     IFeatureClassName pInFeatureClassName = default(IFeatureClassName);     IDataset pDataset = default(IDataset);     IWorkspaceName pOutWorkspaceName = default(IWorkspaceName);     IFeatureClassName pOutFeatureClassName = default(IFeatureClassName);     IDatasetName pDatasetName = default(IDatasetName);     IExportOperation pExportOperation = default(IExportOperation);     IWorkspace pInWorkspace = default(IWorkspace);          pDataset = (IDataset)pFeatureClass;     pInFeatureClassName = (IFeatureClassName)pDataset.FullName;     pInWorkspace = pDataset.Workspace;     pDataset = (IDataset)pOutWorkspace;     pOutWorkspaceName = (IWorkspaceName)pDataset.FullName;          pOutFeatureClassName = new FeatureClassNameClass();     pDatasetName = (IDatasetName)pOutFeatureClassName;     pDatasetName.Name = strLayerName;     pDatasetName.WorkspaceName = pOutWorkspaceName;          pExportOperation = new ExportOperation();     pExportOperation.ExportFeatureClass((IDatasetName) pInFeatureClassName, null, null, pGeometryDef, pOutFeatureClassName, 0);          //Cleanup     pInFeatureClassName = null;     pDataset = null;     pOutWorkspaceName = null;     pOutFeatureClassName = null;     pDatasetName = null;     pExportOperation = null;     pInWorkspace = null;     }


Why does arcpy append from file geodatabase to SDE fail?

I'm trying to append a feature class from a file geodatabase to an Arcsde featureclass using this line of code: `arcpy.Append_management("parcel_copy", sde_path+"parcel" ,"NO_TEST")` This line works under Python console in ArcMap , but when I execute it under Pyscripter IDE, all the code works fine except the append statement, and it does not return an error. * * * Thanks for you responses, i have tried to isolate the `append` statements in a little script, and it worked, but i steel don't understand why it has not worked. i wonder if i have not deleted a featurelayer refering to a featureclass in SDE, that locks this featureclass, so i cannot append to it. any explanation is welcome


Looking for gis server for pictometry that can do bing-like "birds eye" tilted view and rotation

I need an open source gis server that can present these pictometry images, I guess they are also called 2.5D. The camera on a plane takes 5 shots, one orthogonal, other 4 in all directions. North, west... 2.5D images are the main focus here, user would have to be able to rotate the tilted view in 4 directions. Additionally I would like to have those standard gis server features like vector layers, queries, labels, redlining... So far I have decent experience with mapguide, open source and the the autodesk one, and I can't recall those features. No rotation or tilting at all, I might be wrong. I have very basic xp with geoserver and mapserver. I'm looking for advice in which direction to go with research. I was thinking geo-server is the best bet, but what interface... Also, I cant really say that I've been working with any kind of 3d desktop gis. I would take some advice on that too. Any suggestions? Thx.


GRASS r.out.mat memory issues

we have built a model with pythn/GRASS that calls the r.out.mat module. The area is quite large 25,000x25,000 cells. The machine we are running it on is Windows 7 machine with 16GB of RAM. After a while the script aborts complaining about an out of memory issues. Now, I would not have expected to have such issues with so much RAM even though the manual page for r.out.mat does say that the module loads the whole map into RAM so it requires a huge amount of RAM. Any idea on how we can optimise the script avoid such problems? If this is a limitation of the module, I was wondering whether r.out.mat could be used to read the map area in patches, dump out the results in single files and then some other module patches all the r.out.mat results together. Any suggestions on how we can overcome this issue? Thanks Antonio


Arcpad 10...using sql with VBScript

I am trying to get a count of related records (Tableform) in order to validate Editform entry (using VBS). Any help getting myCount? I Know there are some syntax issues here, but is this even possible in Arcpad axf? Seems simple enough, no brainer using VBA but something like this:               Sub Validate()         myCount = (Select count(*) from tblRelatedRecs where ID = Pcontrol.value Group by ID            If myCount > 5 then             Application.MessageBox ("Huma, you erred!")             ThisEvent.Result = False         End If     End Sub


Zonal statistics returns wrong values in QGIS

I used zonal statistics for a simple task: to transmit raster pixel values to vectormap which contains squares as big as the pixel size is. The problem is that the values of the outcome are not correct, for example: pixel value on the rastermap is 11, but when I look the same square on the vectormap, it is 76. What can cause the problem? Thank you for answering!


Problems with other applications after installing QGIS 2.0.1

Since the official QGIS installation 2.0.1 (Virtual Server) various core applications from other departments in the municipality do not work always properly. Does anyone know the problem? Thank you for your help Best regards


Viewer Function

Is there a Viewer function/display in QGis as it is in ArcGIS. Or is there an other possibility to limit the access (just visualisation, no analysis). Thank you for your answer.


Snippets applicable for ArcObjects 10.2.1

I am trying to create ArcGIS desktop add-in using Visual Studio Professional 2012. But when I use ArcGIS Snippets Finder to insert code, I found the applicable versions as follows: ' Applicable ArcGIS Product Versions: ' 9.2 ' 9.3 ' 9.3.1 ' 10.0 But I am using ArcGIS 10.2.1, and the ArcObjects obtained from ArcGIS 10.2.1 Desktop installation CD. I want to know where can get the snippets applicable for ArcGIS 10.2.1 (or above). And what version of Visual Studio can support ArcObjects for 10.2.1?


Socket Server to receive GPS location data

I am working on a web application which receives data sent from trackers via TCP Packets containing location information. Up until this point we have been working with dummy data within the mysql database but we are now ready to start building the server which will receive data from the remote devices. As we are just exploring our options at the moment, I thought I would see does anyone have any suggestions on the best route to take when setting up a server to receive TCP packets of data and then store them in a remote database? Our general requirements are as follows:   * Hosted on Amazon EC2   * Allow multiple concurrent connections to receive data (potentially in the thousands)   * Store data in remote MySQL database Any guidance would be greatly appreciated!


ImageLab Program Problem

I am new to GeoTools. I was trying the samples in http://docs.geotools.org/latest/userguide/tutorial/raster/image.html.   I am not able to run the ImageLab sample code. I know it would be some lack in understanding or a simple mistake. I am getting the following exception: Exception in thread "main" java.lang.UnsupportedOperationException: Trying to get a reader from an unknown format. at org.geotools.coverage.grid.io.UnknownFormat.getReader(UnknownFormat.java:62) at org.geotools.demo.ImageLab.displayLayers(ImageLab.java:93) at org.geotools.demo.ImageLab.getLayersAndDisplay(ImageLab.java:77) at org.geotools.demo.ImageLab.main(ImageLab.java:51) What actually needs to be done? Also, is it possible to use normal jpeg image or any other way around? In the GeoTools website the mentioned bluemarble.tif is not available and I have been using normal jpg image. In another website while seraching for more info on this, I read that it uses a world file. What is a world file and how it is related to this and how do I use it? Please help, I have been stuck with this for a long time. Thanks in advance.


Naming layers in DXF using OGR2OGR

We have been using ogr2ogr to export data from PostGIS into DXF. Using the default behaviour of the export, all data are put into layer '0'. We would put data into a layer name from the PG table name. It is not really clear if the correct solution is manipulating the header.dxf under the GDAL library or if this can be achieved by using some of the command line options from ogr2ogr. Thanks for any input. Rgs. Ulrik


Postgres Log Error

I am using Geoserver + GeoWebCache + Postgres for an application. On the application side there is no bugs but on postgres I am getting this log. Please explain and help me to solve this > 2014-03-14 10:14:50 IST ERROR: relation "public.gt_pk_metadata" does not > exist at character 15 2014-03-14 10:14:50 IST STATEMENT: SELECT * FROM > public.GT_PK_METADATA 2014-03-14 10:14:50 IST ERROR: relation > "public.gt_pk_metadata" does not exist at character 15 2014-03-14 10:14:50 > IST STATEMENT: SELECT * FROM public.GT_PK_METADATA 2014-03-14 10:14:50 IST > ERROR: relation "public.gt_pk_metadata" does not exist at character 15 > 2014-03-14 10:14:50 IST STATEMENT: SELECT * FROM public.GT_PK_METADATA > 2014-03-14 10:14:50 IST ERROR: relation "public.gt_pk_metadata" does not > exist at character 15 2014-03-14 10:14:50 IST STATEMENT: SELECT * FROM > public.GT_PK_METADATA 2014-03-14 10:14:50 IST ERROR: relation > "public.gt_pk_metadata" does not exist at character 15 2014-03-14 10:14:50 > IST STATEMENT: SELECT * FROM public.GT_PK_METADATA 2014-03-14 10:14:50 IST > ERROR: relation "public.gt_pk_metadata" does not exist at character 15 > 2014-03-14 10:14:50 IST STATEMENT: SELECT * FROM public.GT_PK_METADATA > 2014-03-14 10:14:50 IST ERROR: relation "public.gt_pk_metadata" does not > exist at character 15 2014-03-14 10:14:50 IST STATEMENT: SELECT * FROM > public.GT_PK_METADATA 2014-03-14 10:14:50 IST ERROR: relation > "public.gt_pk_metadata" does not exist at character 15 2014-03-14 10:14:50 > IST STATEMENT: SELECT * FROM public.GT_PK_METADATA 2014-03-14 10:14:50 IST > ERROR: relation "public.gt_pk_metadata" does not exist at character 15 > 2014-03-14 10:14:50 IST STATEMENT: SELECT * FROM public.GT_PK_METADATA > 2014-03-14 10:14:50 IST ERROR: relation "public.gt_pk_metadata" does not > exist at character 15 2014-03-14 10:14:50 IST STATEMENT: SELECT * FROM > public.GT_PK_METADATA 2014-03-14 10:14:50 IST ERROR: relation > "public.gt_pk_metadata" does not exist at character 15 2014-03-14 10:14:50 > IST STATEMENT: SELECT * FROM public.GT_PK_METADATA 2014-03-14 10:14:50 IST > ERROR: relation "public.gt_pk_metadata" does not exist at character 15 > 2014-03-14 10:14:50 IST STATEMENT: SELECT * FROM public.GT_PK_METADATA > 2014-03-14 10:14:50 IST ERROR: relation "public.gt_pk_metadata" does not > exist at character 15 2014-03-14 10:14:50 IST STATEMENT: SELECT * FROM > public.GT_PK_METADATA 2014-03-14 10:14:50 IST ERROR: relation > "public.gt_pk_metadata" does not exist at character 15 2014-03-14 10:14:50 > IST STATEMENT: SELECT * FROM public.GT_PK_METADATA 2014-03-14 10:14:50 IST > ERROR: relation "public.gt_pk_metadata" does not exist at character 15 > 2014-03-14 10:14:50 IST STATEMENT: SELECT * FROM public.GT_PK_METADATA > 2014-03-14 10:14:50 IST ERROR: relation "public.gt_pk_metadata" does not > exist at character 15 2014-03-14 10:14:50 IST STATEMENT: SELECT * FROM > public.GT_PK_METADATA


Does the ArcGIS JavaScript API let me build / package my application using Dojo?

I'm currently getting to grips with the new (3.1) ArcGIS JavaScript API and how it bundles Dojo and relies on Dojo's Asynchronous Module Definition (AMD) approach. I'm familiar with AMD and recently completed a project using Require, Backbone, and OpenLayers in combination - overall a very pleasant experience. However I'm worried by   * how tightly-coupled Dojo is with the ArcGIS JavaScript API,   * that trying to use Require instead of Dojo AMD throws an exception, and   * that Dojo AMD doesn't seem to have the same user-base and momentum as Require, making it difficult to research / resolve problems without hitting the forums My main concern is that I will write an ArcGIS JavaScript API application and not be able to build / package it into a single / small number of scripts. This will leave me with an application making 100+ HTTP requests at startup for all of the modules it depends on. I've seen this article covering 3.0 but it doesn't mention the build process and I believe 3.0 is quite different. I've also seen this piece on building / packaging an application that uses the ArcGIS JavaScript API but the comments suggest this is pre-Dojo 1.7 AMD. Have ESRI really only half implemented the AMD approach, and in a way that will actually make my application slower? Why can I not find ESRI documentation on building / packaging a JavaScript API-based application? Any thoughts or experience here would be much appreciated. **EDIT** Based on the answers so far it seems that no, ESRI hasn't provided any information on building your Dojo application after inflating the number of scripts involved. It seems the only useful information available is from non-official sources and slightly outdated (enough to cause problems). This leads to me wonder why more people aren't vocal about the API's new improvements / deficiencies   * am I being too critical / demanding in expecting some more information here?   * is no one using the latest version of the API?   * should I use an older version of the API that doesn't depend on AMD?


Is it possible to set the output of a spatial etl tool in model builder?

I am using the interoperability extension and building spatial ETL tools, but I can't seem to specify the output location of the ETL tool so it's impossible to use this dynamically. Has anyone figured out how to do this?


Keyboard shortcut for GPS data capture in QGIS

I would like to setup keyboard shortcuts for GPS data capture with QGIS. Using GPS Information dialog box, I could use Mouse-Left-Click to "Add feature" and "Add track point", but I would like to assign keyboard shortcuts for them so that I can use a programmable Mouse. Settings --> Configure Shortcuts lists only few commands Is there a way I can call the API behind that and assign to a shortcut key? I appreciate your time. Thanks, Thiru


What multiplier would you use for calculating estimated population?

Given a number of address points in a given area, what multiplier would you use for calculating estimated population? I have seen 2.5 in a few places, but I was wondering what the generally accepted multiplier would be.


ArcGIS/arcpy: Can a variable be set to reference an in_memory workspace object?

This is an offshoot question related to a question I posted earlier. I have the following arcpy code that inputs a Python list of polygon geometry objects to the Union_analysis tool, which then writes the result to the in_memory workspace. The result is then input to the Dissolve_management tool, which in turn writes the result to memory. This all occurs near the end of a function, and the final goal here is to return the single geometry from the Dissolve to the calling function, where it can be passed along. Hence my question: **How can (can I?) I set the Dissolve output to a variable so that I can return it as a geometry object?** The first thing that will be done with it is to read its geometry into a point array to use as the inner ring for another polygon.                   if len(lstHSGeoms) > 1:             #Union geometries into one geometry             arcpy.Union_analysis(lstHSGeoms, "in_memory\\lstHSUnionGeom", "ONLY_FID")             arcpy.Dissolve_management("in_memory\\lstHSUnionGeom","in_memory\\lstHSOneGeom", None, None, "SINGLE_PART")             hsOneGeom = lstHSOneGeom      The last line is the problem here. No surprise, arcpy complains that the global name `lstHSOneGeom` is not defined. Putting in the whole in_memory path as the assignment value, i.e. `"in_memory\\lstHSOneGeom"` just assigns the simple string to the `hsOneGeom` variable. There is no reference to the in_memory object. I originally coded this block as the ArcGIS help files would have you do, which is to create empty arcpy Geometry objects and to set them as the outputs of the tools. When then used as inputs, they were rejected as invalid. At least I could keep a reference to the geometry objects in a variable, but the tools just wouldn't recognize them as valid input. So, I went with in_memory workspace update that you see in the snippet above. (Thanks to @dmahr for this solution in the previous question referenced above.) I now don't get errors from the tool, but I do get an error from the last line. Any suggestions? BTW, I should mention the vitals: ArcGIS/ArcView 10.0, Windows 7, 8 GB RAM.


Disable Certain MapTips in Graphics Layer for ArcGIS Silverlight

I am using the ESRI Silverlight API to allow users to use the Identify functionality on a map. The graphics layer on this map is set to use Map Tips, and a pin will be applied to the graphics layer when a user clicks on the map. However, in specific cases, I want the map tips to be disabled. I am using the MouseEnter event on the graphics layer to catch each graphic as the user rolls over them. Once inside, I can determine whether or not I want the Map Tip to show up based on the graphics' attributes. Is there a way for me to disable the Map Tips inside of the MouseEnter event? I don't want to disable them for the entire graphics layer. Thanks, Seth


Getting points within a range from a postgres database?

ok so this is query i'm making which returns no result:               SELECT          building.gid     FROM          building     WHERE         ST_Dwithin(building.geom, ST_GeomFromText('POINT(long lat)',4326), 400)      If I run this in QGIS, I get a syntax error, but if I run this on pgAdmin, i get no error, but simply 0 rows.... my long and lat are in degrees. what could possibly go wrong?


How to convert ESRI SHP to KML with correct coordinate system?

I am working on this SHP file, containing geographical district of Hong Kong. I would like to export this to KML and export it to Google Map for overlay. After I converted it to KML in QGIS (just save as kml format) and opened it in Google Map, it seems the coordinate is out of the entire map world region. Is there any conversion I need to do from ESRI SHP coordinate system to Google Map compatible coordinate system? Thank you for any pointer!


How to use Google Maps API V3 to diplay circle with Openlayers

I found bellow code to display Circle by click event with Radius 60.0m.(google maps API v3). http://stackoverflow.com/questions/15111403/google-maps- api-v3-change-circle-radius-by-click-event I want to display a circle (within a limit of 3000m, 5000m ) and the popup window to appear attribute details of buffer zone. (school name and the school type & etc.,). Please help me to properly place a this JavaScript code into my HTML page.


How to realize two overlaying layers? One with Popups, the other one with Info-div!

I'm trying to display two overlaying layers (geojson). Layer1 should give me information by hovering over it (works fine based on the Leaflet Choropleth Map Example. Layer2 is situated within Layer1 and should give me information by clicking which does not work. It seems as if Layer1 is in front and disables the interaction with the other one. I already tried the `bringToFront()` and `bringToBack()` methods to handle the problem. To checkout what I mean look here: http://jsfiddle.net/uFPwB/2/


Are there any resources or standards when writing a GIS data specification?

I have been searching for sources of GIS data specifications to assist with documenting a number of datasets that I'm currently working on. I'm used to looking at Ordnance Survey specifications, such as the one located here: http://www.ordnancesurvey.co.uk/docs/technical-specifications/os-mastermap- water-layer-technical-specification.pdf. I'm guessing that a good specification will be one which is independent of applications and formats, and describes the data in the detail which is necessary for its purpose and intended audience. Are there any resources to assist with the creation of GIS data specifications? Is there an industry standard way to do them so that they can be easily understood by both users and GIS developers? I'd also be interested to know what kind of documentation people do generally create for datasets which are to be widely used within their organisations. Many thanks.


Where is the network message log tab in QGIS on Windows7?

I want to find out what web requests is QGIS exactly sending to my web server when it connects to a WMS and WFS layer. In my Linux installation of QGIS 2.0 I can see three tabs in the Log Messages window: General, Plugins, Network. In my Windows installation of QGIS 2.0 I can only see two tabs in the Log Messages window: General, Plugins. How do I activate the "Network" log messages tab in QGIS 2.0 on Windows 7?


Is it possible to use .bundle and .bundlex files in QGIS?

I am working with a colleague that has given me UK OS data and Aerial image files to work with in .bundle and .bundlex form. This isn't a huge amount of data as we are primarily working to the 1:50,000 scale to show what is going on in a forest environment. I need to open this and so far am getting nowhere. Initial searches reveal nothing on QGIS. I did find something about opening it with ESRI GIS Desktop, but I am not on that system. If this can't be done with QGIS, I have an ArcGIS 3.2, on an old machine at home. Is that a possibility? Although I think I will run into speed and memory problems on that four year old laptop. All the files are on my Macbook running QGIS 1.8. (QGIS 2.0 not supported with OS 10.9.1) OK... Nothing like Saturday night for chasing down the meaning of arcane GIS terms. The compact cache groups many tiles together in one large file called a bundle. A bundle contains up to 16,384 tiles. The result is a cache with dozens or hundreds of files, instead of thousands or millions. If you look at a compact cache on disk, you can see the bundle files with the .bundle extension. You'll also see some corresponding index files with the extension .bundlx. It seems all ESRI GIS product read .bundle and .bundlex files? http://resources.arcgis.com/EN/HELP/MAIN/10.1/index.html#//0154000002Z3000000.HTM * * * FINAL UPDATE: Curious to note that having posted this here, and talked to a half dozen others I get no answers. The person who gave me the files is hitting a wall, doesn't know what it is. It works fine on her system, but we can not get anything out that will open in QGIS. We started to export AI files but as the files are scalable we need to printout dozens and dozens of sections on her system and she has no time to do that. (What happened to MrSid? it was much easier.) Upshot question here.... Is this -really- such an oddity that nobody knows what it is or how to come at it? The slightest guidance would be most helpful. I am now pressing up against a deadline. \--> Meanwhile I've unpacked a Six year old PC and spent two days cleaning it out and downloaeding the 60 day eval version of Esri, to have a go with that this afternoon. I don't have high hopes about this either to be honest. But we shall see. Nothing is easy at the moment, but if it were easy...


Looking for examples of bayesian modelling and smoothing in a GIS

I am looking to assemble clear examples of Bayesian modelling / smoothing in a GIS. I have yet to see GIS being used to conduct Bayesian style analytical models and would like to see what is currently out there. I am interested in both the application of Bayesian modelling and its current state of inclusion in traditional and/or non-proprietary (opensource) software. At the moment I am very interested in Bayesian modelling for transportation planning, specifically transit ridership and modal split/share estimations.


Striping Artifact in Aster DEM Mosaic

I'm working with SRTM and ASTER data for a university project. My goal is to delineate watersheds and evaluation of these with ArcGis 10 and ArcHydro. But a striping error occurs with the flow direction or hillshade as seen in the attached picture. I have set advanced arcmap setting to Cubic and also projected the data with cubic resampling. ![](http://i.imgur.com/3UVnrEC.jpg) This Strange gdaldem hillshade artifacts is a smiliar question, but for QGis :( Can I get rid of this? And if not, does this alter the following steps after flow direction?


Is there a way to keep labels and symbols from displaying inside an extent indicator?

I've been working on a utilities map and have been asked to turn off the labels and symbols within the extent indicator on the larger scale map and show them only in the detailed map. Is there a way to do this without converting labels to annotation and making a copy of and editing each individual shapefile specifically for this map?


MapInfo Professional 11.5 sample data

would you be able to point me in the direction of finding a copy of the sample data that comes on the CD's for learning MapInfo Professional 11.5?


Two boundary sets intersect - some are genuine diffferences, others slivers resulting from poor digitising. How to differentiate?

**Problem:** I have two boundary sets. Intersecting reveals that there are genuine differences between them - as expected. However, there are also a host of slivers due to their being digitised by different agencies. These mismatches may have statistical implications. I have tried several approaches to clearly differentiate one mismatch type from the other. **Current thinking:**   1. I get the MBR of each polygon in boundary set 1 and get the ratio of the area of the MBR to the area of the actual polygon. I do the same for boundary set 2.   2. I get the MBR of each intersection (these have been disaggregated to singlepart) and get the ratio of the area of each intersection's MBR to the area of the actual intersection.   3. The distribution of the ratio for slivers is likely to be materially different to the distribution of the ratio for genuine differences.  **QUESTIONS** Do you agree that this seems like a sensible approach? How easy is it to obtain the MBR (minXmaxXminYmaxY) for a polygon in ArcGIS 10.1 Advanced? Either point-and-click or Python suggestions welcome. If my thinking is flawed, please point out how, and suggest an alternative.


Atmospheric correction for Landsat 8 in ERDAS?

I am attempting to run some classification on Landsat 8 imagery, however, I need to perform atmospheric correction on the images in order to reveal true spectral signatures. I am using ERDAS, and the autonomous atmospheric correction tool does not support the Landsat 8 sensor. Unfortunately, I do not have access to ENVI which supports landsat 8, or a licence to ATCOR, which seems to be a ERDAS add on for atmospheric correction. Any Ideas? Thanks


How to repair Geometric Network in Versioned ArcSDE Geodatabase?

Our Versioned ArcSDE workflow has a feature dataset with Geometric Network in the parent (Default), which users derive in their own versions. It is now observed that in some versions, the geometric network is no longer valid and does not abide connectivity rules.   * How can we repair connectivity/enable Geometric Network again in these versions?    * Or should we do this on the default dataset, rather than in user versions?    * What are the tools which may help us to retrieve the Geometric Network ?


How to map and color localities (provinces) with an external CSV using Quantum GIS

I use a Swiss shape file map (called PLZO_PLZ) available on http://www.cadastre.ch/internet/cadastre/fr/home/products/plz/data.html (ZIP file Shape LV03) PLZO_PLZ contains a column called PLZ which correspond to a swiss locality. For instance number PLZ 1202 is in Geneva. By the way, PLZ can also be called NPA (same meaning). My aim is to get a colored map according to the NPA frequency coming from an external CSV file: http://pastebin.com/raw.php?i=iAjV88KQ For instance NPA not appearing in the CSV should be not colored (white) and NPA coming few time should appear in light blue and the more frequent the NPA is the more dark blue should be the map. According to my CSV NPA 1203 should be the darkest. Is it possible to obtain what I am looking for with Quantum GIS? (I tried to use "Joins" function in proprieties but it's not working) ![enter image description here](http://i.stack.imgur.com/1ygmV.png) **Edit:** I try to join like this : ![enter image description here](http://i.stack.imgur.com/HP9dV.png)


Practical Building Extraction Algorithm Implementation

I've been reading many papers over 3D building extraction using LiDAR data and Aerial images. I see that there are many different algorithms for doing such things, such as this and this. Those papers, as far as I can see, describe mostly abstract concepts and math. I am a software developer, so how could I practically implement said algorithms into a piece of software? Also, what algorithm is best for extracting building data with minimal vegetation interference?


Where to obtain free 15m ASTER image for New Mexico?

I am looking for a free Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) image that will let me to perform NDVI for the pilot crop health/soil moisture that can go as far back ten years or so. It needs to have:   * 15 m (meters)   * temporal   * Ten years back or so   * Perform NDVI This will be in a small area somewhere near Albuquerque, NM.


Mimimum distance between points (Postgres/PostGIS)

I have a table (tab_Cities with cities which own more than 50000 inhabitants. The table also contains lat, long and the_geom as columns. I can already see the points via QGIS. Now I want to write a postgresql query which gives me the minimum distance of the points but I only now how to get the distance of two points. But now I have around 60 points and I want to have the smallest distance of ALL points. How can I do this? I found some similar questions in other threads but they were using two tables. I am using one table only. Thank you for help! :-).


PyQgis question: How to change or set active layer?

Developing a qgis plugin, I found that if no layer is selected in the Layers panel I can iterate over all layers and do anything to them, but "activate" them. Could someone give me a hint please? Thank you!


Read layer symbology class color with ArcPy

I am using ArcPy within ArcGIS 10.1 and attempting to read symbology from a *.lyr file created by a raster dataset within ArcMap. My goal is to be able read out the RGB (color) values assigned to the class breaks on the dataset. It’s seems easy enough to get the lyr file into an object with:                   lyrFile = arcpy.mapping.Layer("C:/.../ReadLyrFileSymbology/conus_foresttype6.img.lyr")      However, after that I am unable to find the method or property to get symbology color (RGB) on a layer. Any pointers would be greatly appreciated.


Get Cell Value Webservice for ArcGIS Raster (LiDAR) without using ArcGIS

I want to expose the getCellValue function on Lidar Data, as a web service. But I don't have access to ArcGIS (server or desktop). What options do I have to fetch cell value, given lat lon, from a raster in ArcGIS format? It would be great if someone can point me to some open source code/library that I can use to create a web service. I am open to any other option except ArcGIS. I know that Geoserver has a raster rendering feature: can it also find the cell value for a lat-long input?


Convert Lambert Conformal Conic to WGS84

I am trying to convert my "Lambert Conformal Conic map with 2 standard parallels" to the basic/standard WGS84 projection that Google Maps use. I have 30 of these maps which I need to convert. They are all in `.tiff` format with a corresponding `.tfw` file each. What my final goal is:   1. Convert Maps to WGS84   2. Merge all 30 maps   3. Change them to tiles For now I am only struggling with the conversion of projection. Does anybody now of any ways/programs/techniques or anything (preferably free) to accomplish this? _I have heard/read about GDAL etc. but I do run Win Vista, not MAC._


Geowebcache "seed" doesn't seem to work

I'm using Geowebcache seed command (via web interface at `http://localhost:8080/geoserver/gwc/rest/seed/workspace:layerName`) and it looks like everything runs smoothly but actually no tiles are saved, at least where I'd expect them to be saved `(/data_dir/gwc)`. I have GeoServer installed on top of Tomcat 7 as WAR. I also checked the logs (enabling `GEOTOOLS_DEVELOPER` mode) and it doesn't give any errors. Actually if enabling `VERBOSE` logging mode I also get something like...               ...         2014-07-07 09:38:08,766 DEBUG [wms.map] - Writing png image ...     2014-07-07 09:38:08,771 DEBUG [wms.map] - Writing png image ... done!     2014-07-07 09:38:08,772 DEBUG [wms.map] - Writing png image ...     2014-07-07 09:38:08,776 DEBUG [wms.map] - Writing png image ... done!     2014-07-07 09:38:08,777 INFO [seed.SeedTask] - GWC Seeder Thread-1 completed (re)seeding layer workspace:layerName after 448 tiles and 45 seconds.     ...      ...and it look like everything went fine. But in the end I have no tiles. By the way, the tiles I'm trying to save to EPSG:900913 projection coming from EPSG:3045 PostGIS store (already done this operation successfully on another VM having GeoServer installed on top of Tomcat6 via `apt-get` instead than via WAR). Any suggestion on how to inspect/solve the issue? **Update:** Following @EvilGenius advice I was able to find the missing tiles (in fact they quite strangely were served by the Geowebcache to the frontend) in `/tmp/tomcat7-tomcat7-tmp/geowebcache/workspace_layerName`. Now the question becomes: how comes they are not copied/moved to the final location (`data_dir/gwc` folder)? again, how can I inspect that?


Unable to locate PostGIS in OpenGeo Suite using command line

I am working on this new ubuntu 12.04 server and I need to work on the Postgres database which is PostGIS enabled. The problem is I am not sure how to locate the database using command line. The server already has OpenGeo suite installed on it and so it definitely has PostGIS. It is already installed on the server. I just need to locate it. How do I locate it?


Annoying little line in color ramp legends - how do I remove it?

This was never a problem with ArcGIS 9 but it appears that in 10.0 Esri has added small horizontal lines as the default when including color ramps in legends. I want to get rid of these annoying lines, but I haven't figured out how to. Anyone else run into this problem?


How to choose selectable layers in ArcMap 10?

Maybe I'm missing something, but I cannot figure out how to toggle selectable layers in ArcMap 10. In 9.3, I remember being able to choose from a list which layers I wanted to be able to select features from in the map. Let's say I'm editing two different feature classes for a building - 'FireExtinguishers' and 'EvacuationRoutes'. While editing, I don't want to accidentally select the 'Rooms' feature class. I want to be able to turn off the selectability of 'Rooms', but the only way I know how to change this is to right-click a layer in the Table of Contents and choose Selection -> Make This The Only Selectable Layer. I would like to have finer control over what layers are and are not selectable. Does anyone know how to do this?


Building own OSM server. Searching/routing/rendering

I need to build my own OSM server to reach customers from other countries. I have the following requirements for the resulting system:   1. Search:   1.1. "Hierarchical" (full hierarchy from country to building/object - country > admin level 3 > admin level 4 > ... > settlement(city/town/...)[ > suburb] > street > building/object):     * It is possible to use "placex" table produced by Nominatim. What is your experience with it? Are there any alternatives? 1.2. Fuzzy search:     * Are there any working projects/examples for fuzzy search on OSM/GIS data?     * Is there any list of the frequently used words (for example: street, avenue, restaurant, hotel) for different languages/locales?      * Are there any use for such libraries as Sphinx or Apache Lucene etc.?  1.3. "Free" search (one parameter as "q" in Nominatim):     * The only available project that I've found is Nominatim. Are there any alternatives?   2. Routing (fastest route, possibly taking traffic jams into account):      * What is your "favourite" engine?     * OSRM seems to be the fastest one. How does it handle OSM diffs?   3. Rendering/Serving tiles:      * Do you run Mapnik after every OSM diff?


pythonaddin button return value

I have asked this question a while back but within arcgis environment without using addin. Now i am trying to do the same task to get pixel value when user click anywhere on screen, to use in later steps. Following pythonaddin button code snippet returns nothing. Any suggestions where is the catch?               class Val(object):     """Implementation for Pythonaddin_addin.genW (Button)"""     def __init__(self):         self.enabled = True         self.checked = False     def onClick(self):         geog = arcpy.GetParameterAsText(0)         ras = arcpy.PointToRaster_conversion(geog, "ID", "in_memory/ras", "MAXIMUM", "",30)         pass


Zoom to link button in Popup, killing the Map

I am using the traditional Popup of JavaScript API and it seems to not work properly in my application. When you click on the map, the Popup comes fine but when you click on the Zoom to link button, the map zooms in and then if you move the map, the map disappears. It only returns when you refresh the app. Not sure what is causing it?


Reproject JPEG to TIFF?

I want to import an elevation layer into GeoServer for later WMS access. This is the base layer I would like to import - http://en.wikipedia.org/wiki/File:Srtm_ramp2.world.21600x10800.jpg . What I'd like to do is reproject the layer to Mercator so it looks like this http://upload.wikimedia.org/wikipedia/commons/thumb/f/f4/Mercator_projection_SW.jpg/350px- Mercator_projection_SW.jpg , then somehow convert it to GeoTIFF so that I can import it into GeoServer. How should I proceed? I tried using gdal_warp to reproject my layer but it warps the layer completely leaving nothing behind. I had some success reprojecting the layer using G.Projector, but it doesn't work for bigger images and it doesn't georeference the image anyway (GeoServer reports no layers in GeoTIFF file).


gdal_translate wms

Help me, pls, vs gdal_translate and wms i created this xml conf file `<GDAL_WMS> <Service name="TMS"> <ServerUrl>http://213.168.24.219/tilecache/tilecache.cgi/1.0.0/</ServerUrl> <Layer>maorto</Layer> <Format>jpeg</Format> </Service> <DataWindow> <UpperLeftX>856520.000000</UpperLeftX> <UpperLeftY>6799520.000000</UpperLeftY> <LowerRightX>365000.000000</LowerRightX> <LowerRightY>6308000.000000</LowerRightY> <TileLevel>12</TileLevel> <TileCountX>2</TileCountX> <TileCountY>1</TileCountY> <YOrigin>top</YOrigin> </DataWindow> <Projection>EPSG:3301</Projection> <BlockSizeX>256</BlockSizeX> <BlockSizeY>256</BlockSizeY> <BandsCount>3</BandsCount> </GDAL_WMS>` But it doesnt work =( what did i do wrong? _edited to add error message_ :               c:\site>gdal_translate -of gtiff frmt_wms_Est_tms.xml 1.gtiff     Input file size is 2097152, 1048576 0     ERROR 4: `/vsimem/wms/000000001911DDF0/wmsresult.dat' not recognised as a supported file format.      ERROR 1: GDALWMS: Unable to open downloaded block.      ERROR 1: GDALWMS: ReadBlockFromFile (213.168.24.219/tilecache/tilecache.cgi/1.0.0/maorto/12/0/…) failed.      ERROR 1: frmt_wms_Est_tms.xml, band 1: IReadBlock failed at X offset 0, Y offset 0     ERROR 1: GetBlockRef failed at X block offset 0, Y block offset 0      * * * unactual question. this is more simple i thought.                               c:\site>gdalinfo http://213.168.24.219/tilecache/tilecache.cgi/1.0.0/orto                 Driver: WMS/OGC Web Map Service                 Files: none associated                 Size is 1048576, 1048576                 Coordinate System is:                 PROJCS["Estonian Coordinate System of 1997",                     GEOGCS["EST97",                         DATUM["Estonia_1997",                             SPHEROID["GRS 1980",6378137,298.257222101,                                 AUTHORITY["EPSG","7019"]],                             TOWGS84[0,0,0,0,0,0,0],                             AUTHORITY["EPSG","6180"]],                         PRIMEM["Greenwich",0,                             AUTHORITY["EPSG","8901"]],                         UNIT["degree",0.0174532925199433,                             AUTHORITY["EPSG","9122"]],                         AUTHORITY["EPSG","4180"]],                     PROJECTION["Lambert_Conformal_Conic_2SP"],                     PARAMETER["standard_parallel_1",59.33333333333334],                     PARAMETER["standard_parallel_2",58],                     PARAMETER["latitude_of_origin",57.51755393055556],                     PARAMETER["central_meridian",24],                     PARAMETER["false_easting",500000],                     PARAMETER["false_northing",6375000],                     UNIT["metre",1,                         AUTHORITY["EPSG","9001"]],                     AUTHORITY["EPSG","3301"]]                 Origin = (365000.000000000000000,6692000.000000000000000)                 Pixel Size = (0.366210937500000,-0.366210937500000)                 Image Structure Metadata:                   INTERLEAVE=PIXEL                 Corner Coordinates:                 Upper Left  (  365000.000, 6692000.000) ( 21d33'20.83"E, 60d20'27.98"N)                 Lower Left  (  365000.000, 6308000.000) ( 21d47' 6.90"E, 56d53'45.84"N)                 Upper Right (  749000.000, 6692000.000) ( 28d30'12.36"E, 60d17'17.35"N)                 Lower Right (  749000.000, 6308000.000) ( 28d 4'53.13"E, 56d50'52.99"N)                 Center      (  557000.000, 6500000.000) ( 24d58'53.39"E, 58d38' 9.97"N)                 Band 1 Block=256x256 Type=Byte, ColorInterp=Red                   Overviews: 524288x524288, 262144x262144, 131072x131072, 65536x65536, 32768x32768, 16384x16384, 8192x8192, 4096x4096, 2048x2048, 1024x1024, 512x512, 256x256                 Band 2 Block=256x256 Type=Byte, ColorInterp=Green                   Overviews: 524288x524288, 262144x262144, 131072x131072, 65536x65536, 32768x32768, 16384x16384, 8192x8192, 4096x4096, 2048x2048, 1024x1024, 512x512, 256x256                 Band 3 Block=256x256 Type=Byte, ColorInterp=Blue                   Overviews: 524288x524288, 262144x262144, 131072x131072, 65536x65536, 32768x32768, 16384x16384, 8192x8192, 4096x4096, 2048x2048, 1024x1024, 512x512, 256x256      and xml with this data works. but output image is very large. How can i download it by small pieces?


Setting up THREDDS catalogs for ocean model data

I've got a bunch of NetCDF files from ocean modeling runs. Some are 3-day forecasts that are run daily, some are hindcast runs where each hindcast is split into multiple files like `ocean_his_0001.nc`, `ocean_his_0002.nc`, etc, and some are just random netcdf files required for input (forcing files, grid files, etc). I would like to serve these data using the Unidata THREDDS Data Server, so my data can be made available via OPeNDAP, WCS and WMS, but the documentation is a bit overwhelming. Are there some simple suggestions for how to set up THREDDS catalogs to handle this diverse collection of model-related netcdf files?


MapInfo: How to get the population contained within 2 geographical tables

I using 3 different tables which are also geographical objects:   * population data; points on the map representing inhabitants (1 or more)   * region data; the map is divided in regions   * coverage data; it contains several coverage areas overlapping or not What I am searching for is: What is the percentage of people covered by a certain coverage area in a certain region? I tried several things but nothing worked, for example with SQL select I tried: Select Columns: Regions.name, CoverageArea.name, Sum(Population.inhabitants) from Tables: Regions, CoverageArea, Population where Condition: Regions.obj Contains Population.obj And CoverageArea.obj Contains Population.obj Group by Columns: 1, 2 Order by Columns: 1, 2 into Table Named: Selection "And I get Incorrect tables are joined. Invalid join condition in WHERE clause." Hope I was clear enough :s ** edit ** Few tests I did. I can't understand why the thrid SQL query doesn't work... ![Population per coverage area](http://i.stack.imgur.com/9Qe9x.jpg) ![Population per region](http://i.stack.imgur.com/WBJgs.jpg) ![Trying to get population per region and coverage area](http://i.stack.imgur.com/7oKtu.jpg) ![Error message](http://i.stack.imgur.com/HYb3V.jpg)


Edit shapefile using modelbuilder

I have got the folloving problem: I need to change some of the values of one attribute in my point shapefile. I want to do it in modelbuilder, but I dont know how. I have found a script which can do it in python. The script is as following:               import arcpy     ...      ... features = arcpy.UpdateCursor(r"C:\Users\peter.valent\Documents\ArcGIS\Data\Vah po LM\test_points.shp")     ... for feature in features:     ...         feature.precipitation = value     ...         features.updateRow(feature)     ... del feature,features      This is what I put to modelbuilder Expression:                myfcn(r"%file_name%","%value%")      Code Block:               import arcpy          def myfcn(path,value):     features = arcpy.UpdateCursor(path)     for feature in features:             feature.precipitation = value             features.updateRow(feature)     del feature,features


How do I snip/split polylines at chainages?

What is the best way to snip/split lines at chainages? E.g I have a csv with the uniq ID of the polyline item and a start and end chainage. I want to extract the lines between the chainages stated. I have access to MapInfo Professional, QGIS and ArcGIS for Desktop, and I am also OK at using MapBasic.


WMS Query does not return all records (Geoserver)

This is a tricky one - so bear with me... i'm a Geoserver newbie... Basically, i have two very similar queries that return species records for a given species in a given geograpic area. One query uses WFS (returns JSON data that I then parse and map on Google maps); the other uses WMS and returns a map with the records just plotted as points. Despite there being no obvious filtering on the WMS (apart from the species name and bounding box), the WMS only returns a fraction of the records. I'm running Geoserver 2.3.3 on Windows server 2008. All the data I am querying is in shapefiles. Here is the WMS query:                   <?xml version="1.0" encoding="UTF-8"?>     <ogc:GetMap xmlns:ogc="http://www.opengis.net/ows" xmlns:gml="http://www.opengis.net/gml" version="1.1.1" service="WMS">          <StyledLayerDescriptor version="1.0.0">          <NamedLayer>         <Name>wales_lrcs</Name>           <NamedStyle>             <Name>polygon_empty_red</Name>         </NamedStyle>     </NamedLayer>          <NamedLayer>         <Name>bis_4_other_species</Name>           <UserStyle xmlns="http://www.opengis.net/sld">           <FeatureTypeStyle>           <Rule>           <Filter xmlns="http://www.opengis.net/ogc">             <And>               <ogc:PropertyIsLike wildCard="*" singleChar="." escape="!">               <ogc:PropertyName>SPECIES</ogc:PropertyName>               <ogc:Literal>Amanita muscaria*</ogc:Literal>               </ogc:PropertyIsLike>             </And>           </Filter>           <PointSymbolizer>             <Geometry>               <PropertyName xmlns="http://www.opengis.net/ogc">the_geom</PropertyName>             </Geometry>             <Graphic>               <Mark>                 <WellKnownName>SQUARE</WellKnownName>                 <Fill>                   <CssParameter name="fill">#0000CC</CssParameter>                 </Fill>               </Mark>               <Size>4</Size>             </Graphic>           </PointSymbolizer>           </Rule>         </FeatureTypeStyle>         </UserStyle>     </NamedLayer>          </StyledLayerDescriptor>          <BoundingBox srsName="http://www.opengis.net/gml/srs/epsg.xml#4326">         <gml:coord>             <gml:X>-5.625</gml:X>             <gml:Y>51.255</gml:Y>         </gml:coord>         <gml:coord>             <gml:X>-2.395</gml:X>             <gml:Y>53.572</gml:Y>         </gml:coord>     </BoundingBox>          <Output>         <Format>image/png</Format>         <Transparent>false</Transparent>         <Size>             <Width>500</Width>             <Height>600</Height>         </Size>     </Output>          <Exceptions>application/vnd.ogc.se+xml</Exceptions>     </ogc:GetMap>      And here is the WFS query:               <wfs:GetFeature      service="WFS"      version="1.1.0"      outputFormat="JSON"           xmlns:wfs="http://www.opengis.net/wfs"      xmlns:ogc="http://www.opengis.net/ogc"      xmlns:gml="http://www.opengis.net/gml"      xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"      xsi:schemaLocation="http://www.opengis.net/wfs      http://schemas.opengis.net/wfs/1.1.0/wfs.xsd">     <wfs:Query typeName="bis_workspace:bis_4_other_species">     <ogc:Filter>         <ogc:And>             <ogc:DWithin>                 <ogc:PropertyName>the_geom</ogc:PropertyName>                 <gml:Envelope srsName="http://www.opengis.net/gml/srs/epsg.xml#4326">                     <gml:lowerCorner>-5.625 51.255</gml:lowerCorner>                     <gml:upperCorner>-2.395 53.572</gml:upperCorner>                 </gml:Envelope>                 <ogc:Distance units="meter">0</ogc:Distance>             </ogc:DWithin>             <ogc:PropertyIsLike wildCard="*" singleChar="." escape="!">                 <ogc:PropertyName>SPECIES</ogc:PropertyName>                 <ogc:Literal>Amanita muscaria*</ogc:Literal>             </ogc:PropertyIsLike>         </ogc:And>     </ogc:Filter>     </wfs:Query>     </wfs:GetFeature>      Map from the WMS query (pretty much only one record returned!): ![WMS](http://i.stack.imgur.com/Wnknk.png) And the WFS points plotted on Google Maps (returns 32 features): ![WFS](http://i.stack.imgur.com/3QM90.jpg) Any help would be very much appreciated! Steve * * * Here is the output from the Geoserver log with developer debugging turned on. Log is for the WMS query - I'm assuming that's where the issue is:                   2013-09-12 17:01:37,784 DEBUG [org.geoserver.ows.OWSHandlerMapping] - No handler mapping found for [/web/]     2013-09-12 17:01:37,784 DEBUG [org.geoserver.ows.OWSHandlerMapping] - No handler mapping found for [/web/]     2013-09-12 17:01:37,784 DEBUG [org.geoserver.ows.OWSHandlerMapping] - No handler mapping found for [/web/]     2013-09-12 17:01:37,785 DEBUG [org.geoserver] - Thread 18 locking in mode READ     2013-09-12 17:01:37,785 DEBUG [org.geoserver] - Thread 18 got the lock in mode READ     2013-09-12 17:01:37,790 DEBUG [org.geoserver] - Thread 18 releasing the lock in mode READ     2013-09-12 17:01:37,790 DEBUG [org.geoserver.security.filter.GeoServerSecurityContextPersistenceFilter$1] - SecurityContextHolder now cleared, as request processing completed     2013-09-12 17:01:37,790 DEBUG [org.geoserver.filters] - Compressing output for mimetype: text/xml; charset=utf-8     2013-09-12 17:01:41,540 DEBUG [org.geoserver.security] - Start searching for expired authentication tokens     2013-09-12 17:01:41,540 DEBUG [org.geoserver.security] - Number of expired authentication tokens found: 0     2013-09-12 17:01:41,540 DEBUG [org.geoserver.security] - End searching for expired authentication tokens     2013-09-12 17:01:43,728 DEBUG [org.geoserver.ows.OWSHandlerMapping] - No handler mapping found for [/web/]     2013-09-12 17:01:43,729 DEBUG [org.geoserver.ows.OWSHandlerMapping] - No handler mapping found for [/web/]     2013-09-12 17:01:43,729 DEBUG [org.geoserver.ows.OWSHandlerMapping] - No handler mapping found for [/web/]     2013-09-12 17:01:43,730 DEBUG [org.geoserver] - Thread 18 locking in mode READ     2013-09-12 17:01:43,730 DEBUG [org.geoserver] - Thread 18 got the lock in mode READ     2013-09-12 17:01:43,738 DEBUG [org.geoserver] - Thread 18 releasing the lock in mode READ     2013-09-12 17:01:43,739 DEBUG [org.geoserver.security.filter.GeoServerSecurityContextPersistenceFilter$1] - SecurityContextHolder now cleared, as request processing completed     2013-09-12 17:01:43,739 DEBUG [org.geoserver.filters] - Compressing output for mimetype: text/xml; charset=utf-8     2013-09-12 17:01:43,949 DEBUG [org.geoserver.ows.OWSHandlerMapping] - No handler mapping found for [/web/]     2013-09-12 17:01:43,949 DEBUG [org.geoserver.ows.OWSHandlerMapping] - No handler mapping found for [/web/]     2013-09-12 17:01:43,949 DEBUG [org.geoserver.ows.OWSHandlerMapping] - No handler mapping found for [/web/]     2013-09-12 17:01:43,950 DEBUG [org.geoserver] - Thread 18 locking in mode READ     2013-09-12 17:01:43,950 DEBUG [org.geoserver] - Thread 18 got the lock in mode READ     2013-09-12 17:01:43,952 DEBUG [org.geoserver] - Thread 18 releasing the lock in mode READ     2013-09-12 17:01:43,952 DEBUG [org.geoserver.security.filter.GeoServerSecurityContextPersistenceFilter$1] - SecurityContextHolder now cleared, as request processing completed     2013-09-12 17:01:43,952 DEBUG [org.geoserver.filters] - Compressing output for mimetype: text/html; charset=utf-8     2013-09-12 17:01:44,037 DEBUG [org.geoserver.ows.OWSHandlerMapping] - No handler mapping found for [/TestWfsPost]     2013-09-12 17:01:44,037 DEBUG [org.geoserver.ows.OWSHandlerMapping] - No handler mapping found for [/TestWfsPost]     2013-09-12 17:01:44,037 DEBUG [org.geoserver.ows.OWSHandlerMapping] - No handler mapping found for [/TestWfsPost]     2013-09-12 17:01:44,037 DEBUG [org.geoserver.ows.OWSHandlerMapping] - Mapping [/TestWfsPost] to HandlerExecutionChain with handler [org.springframework.web.servlet.mvc.ServletWrappingController@ae2984] and 1 interceptor     2013-09-12 17:01:44,041 DEBUG [org.geoserver.ows.OWSHandlerMapping] - No handler mapping found for [/wms]     2013-09-12 17:01:44,041 DEBUG [org.geoserver.ows.OWSHandlerMapping] - No handler mapping found for [/wms]     2013-09-12 17:01:44,041 DEBUG [org.geoserver.ows.OWSHandlerMapping] - No handler mapping found for [/wms]     2013-09-12 17:01:44,041 DEBUG [org.geoserver.ows.OWSHandlerMapping] - No handler mapping found for [/wms]     2013-09-12 17:01:44,041 DEBUG [org.geoserver.ows.OWSHandlerMapping] - Mapping [/wms] to HandlerExecutionChain with handler [org.geoserver.ows.Dispatcher@1958a34] and 1 interceptor     2013-09-12 17:01:44,042 DEBUG [org.vfny.geoserver.requests] - First 4 bytes of XML doc are : 3C ('<') 3F ('?') 78 ('x') 6D ('m')     2013-09-12 17:01:44,042 DEBUG [org.vfny.geoserver.requests] - Charset detection phase 1. Inferred encoding: UTF-8     2013-09-12 17:01:44,042 DEBUG [org.vfny.geoserver.requests] - Charset detection phase 2. Charset in XML declaration is `UTF-8`.     2013-09-12 17:01:44,042 DEBUG [org.vfny.geoserver.requests] - Trying to create reader basing on existing charset information: `UTF-8`.     2013-09-12 17:01:44,042 DEBUG [org.geoserver.ows] - Raw XML request starts with: <?xml version="1.0" encoding="UTF-8"?>     <ogc:GetMap xmlns:ogc="http://www.opengis.net/ows" xmlns:gml="http://www.opengis.net/gml" version="1.1.1" service="WMS">     <StyledLayerDescriptor version="1.0.0">     <NamedLayer>         <Name>wales_lrcs</Name>           <NamedStyle>             <Name>polygon_empty_red</Name>         </NamedStyle>     </NamedLayer>     <NamedLayer>         <Name>bis_4_other_species</Name>           <UserStyle xmlns="http://www.opengis.net/sld">           <FeatureTypeStyle>           <Rule>           <Filter xmlns="http://www.opengis.net/ogc">             <And>               <ogc:PropertyIsLike wildCard="*" singleChar="." escape="!">               <ogc:PropertyName>SPECIES</ogc:PropertyName>               <ogc:Literal>Amanita muscaria*</ogc:Literal>               </ogc:PropertyIsLike>             </And>           </Filter>           <PointSymbolizer>             <Geometry>               <PropertyName xmlns="http://www.opengis.net/ogc">the_geom</PropertyName>             </Geometry>             <Graphic>               <Mark>                 <WellKnownName>SQUARE</W...     2013-09-12 17:01:44,045 INFO [org.geoserver.wms] -      Request: getServiceInfo     2013-09-12 17:01:44,064 DEBUG [org.geoserver.wms] - setting up map     2013-09-12 17:01:44,085 DEBUG [org.geoserver.wms.map] - setting up 500x600 image     2013-09-12 17:01:55,360 INFO [org.geoserver.wms] -      Request: getMap         Filters = null         Tiled = false         Palette = null         SRS = EPSG:4326         FeatureVersion = null         Styles = [StyleImpl[ name=polygon_empty_red], StyleImpl[ name=Default Styler]]         Layers = [org.geoserver.wms.MapLayerInfo@2c1e2b26, org.geoserver.wms.MapLayerInfo@c1e22fea]         MaxFeatures = null         Bbox = Env[-5.625 : -2.395, 51.255 : 53.572]         Env = {}         Angle = 0.0         CQLFilter = null         Elevation = []         FeatureId = null         RemoteOwsType = null         RemoteOwsURL = null         FormatOptions = {}         StartIndex = null         ViewParams = null         Crs = GEOGCS["WGS 84",        DATUM["World Geodetic System 1984",          SPHEROID["WGS 84", 6378137.0, 298.257223563, AUTHORITY["EPSG","7030"]],          AUTHORITY["EPSG","6326"]],        PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]],        UNIT["degree", 0.017453292519943295],        AXIS["Geodetic longitude", EAST],        AXIS["Geodetic latitude", NORTH],        AUTHORITY["EPSG","4326"]]         Transparent = false         ValidateSchema = false         SldVersion = null         BgColor = java.awt.Color[r=255,g=255,b=255]         SldBody = null         Sld = null         TilesOrigin = null         Exceptions = application/vnd.ogc.se+xml         Height = 600         Width = 500         Filter = null         Time = []         Format = image/png         Buffer = 0         Request = GetMap         RawKvp = null         Get = false         BaseUrl = http://localhost:8080/geoserver/         RequestCharset = null         Version = 1.1.1     2013-09-12 17:01:55,361 DEBUG [org.geoserver.wms.map] - Writing png image ...     2013-09-12 17:01:55,475 DEBUG [org.geoserver.wms.map] - Writing png image ... done!     2013-09-12 17:01:55,475 DEBUG [org.geoserver.filters] - Not compressing output for mimetype: image/png     2013-09-12 17:01:55,476 DEBUG [org.geoserver.security.filter.GeoServerSecurityContextPersistenceFilter$1] - SecurityContextHolder now cleared, as request processing completed     2013-09-12 17:01:55,483 DEBUG [org.geoserver.security.filter.GeoServerSecurityContextPersistenceFilter$1] - SecurityContextHolder now cleared, as request processing completed


level of spatial aggregation of polygons

I am looking for a method that allows me to obtain the spatial aggregation of polygons for each cell of a grid with 10kmx10km resolution. In each cell (10kmx10km), I want to be able to consider that big polygons that share a line segment as well as big polygons within certain distances (e.g., 1km), represent areas with high aggregation. Smaller polygons with the same spatial distribution represent a medium aggregation cell. While small disperse polygons represent areas with low aggregation. What is the best way to do this (in ArcGIS, QGIS or R)?


Metric/parameters to compare length of polyline

What would be a good metric to compare the length of two polylines? A very naive approach would be to define a maximum delta which the two lines can differ. I'm looking for more approaches. For example one that might also take a function over the length of the polyline as a parameter. Does anybody has algorithms/functions to share and has even evaluated parameters that work quite good?


Adding information about adjacent features to the data

I am working with ArcGIS 10.0. As a simple example let's say I have a map of the United States. I want to each row of the data which states have adjacent boundaries with that state - so that in the row for South Carolina there would be new variables adjacent_1 and adjacent_2 that list North Carolina and Georgia (or anything similar to this) Any thoughts? I know I can select on "features that touch the boundary of the selected feature" but then I'm not sure how to add the information to the data. Thanks in advance


How to export dxf of window view from postGis MasterMap layer

I have managed to get qgis to show my Mastermap data which is stored in a postGis database. It's national coverage so seems to take ages to load but that's another question. I basically want to export the map that I can see (or at least one layer of it) as a DXF. When I try to do this (layer save as) it starts to generate a temp version of the dxf file but then just hangs and I have to use task manager to close it and start all over again. It's as if it's trying to read the whole database. The postGis database was created using Edonica. Some of the field names are different to those you get if you load a gml file directly in to qgis so this may be the problem. Or do I have to make a selection of the area I want? Thanks for any help Graham


GRASS Add Vector Point

I need to do a viewshed analysis using GRASS. I've got the DEMs loaded; and I simply want to input several (2-4) points (like 41.12345N, 122.34567W from GPS). What is the simplest, most straightforward, way of accomplishing this?


Store elements of Imap into BLOB fields using IPersistStream

I have an IMap (include elements: text, line) and I want to store it into a field of a table in SQL, so I can load it again later. I used below codes, but I got an the error :               public class Test      {       IStream stream;       IMap myMap;       void Save()       {         IPersistStream perStream;         IMemoryBlobStream memoryStream = new MemoryBlobStreamClass();         stream = memoryStream as IStream;         perStream = (IPersistStream ) myMap ;         perStream.Save(stream, 0);         }       void Load()       {          IMemoryBlobStream memoryStream = new MemoryBlobStreamClass();          IPersistStream pstm = new MapClass();            pstm.Load(stream); **// throw error: Error HRESULT E_FAIL has been returned from a call to a COM component.**               myMap = (IMap)pstm;            }     }      **Are there anyone can help ?** Thanks and regards, * * * I found out one way to save elements to an binary file and can load it again later as below codes:               void Save (AxPageLayoutControl pageControl)     {                     IPersistStream perStream;                     IMemoryBlobStream memoryStream = new MemoryBlobStreamClass();                          IGraphicsContainer graphicsContainer = pageControl.ActiveView.GraphicsContainer;                     graphicsContainer.Reset();                     IElement pElement = graphicsContainer.Next();                     while (pElement != null)                     {                         if (pElement is IMapFrame)                         {                             perStream = (IPersistStream)pElement ;                             perStream.Save(memoryStream, 0);                                  memoryStream.SaveToFile("C:/image.bmp");                              break;                              }                         pElement = graphicsContainer.Next();                     }     }          void Load (AxPageLayoutControl pageControl)     {                 IPersistStream perStream;                 IMemoryBlobStream memoryStream = new MemoryBlobStreamClass();                 IGraphicsContainer graphicsContainer = pageControl.ActiveView.GraphicsContainer;                 graphicsContainer.Reset();                 IElement pElement = graphicsContainer.Next();                 while (pElement != null)                 {                     if (pElement is IMapFrame)                     {                              memoryStream.LoadFromFile("C:/image.bmp");                         perStream = (IPersistStream)pElement;                         perStream.Load(memoryStream);                         IMapFrame frame = new MapFrameClass();                         frame   = (IMapFrame)perStream;                         pElement = (IElement)frame;                         break;                          }                                     pElement = graphicsContainer.Next();                 }                 pageControl.ActiveView.PartialRefresh(esriViewDrawPhase.esriViewGraphics, null, null);     }      But, I am finding a way to convert IPersistStream or IMemoryBlobStream to byte[] format, so I can store in SQL Server and load later. Are there anyone know this ? Thanks and regards,


Transform bookmarks to data driven pages in ArcGIS?

We have this large project in ArcGIS with 30 maps with different extents using Bookmarks to pan between these maps/pages. Now we would like to use Data driven pages instead. The question is, is it possible to convert bookmarks to data driven pages or do we have to create a lot of polygons covering each map extent?


How do you get the coordinates from a click or drag event in the Google Maps API V3?

I have made a Google Version 3 Geocoder , I want to be able to pick up the coordinates of the marker when it is dragged or clicked. Below is my code:               <!DOCTYPE html>               <meta name="viewport" content="initial-scale=1.0, user-scalable=no"/>     <meta http-equiv="content-type" content="text/html; charset=UTF-8"/>     <title>Google Maps JavaScript API v3 Example: Geocoding Simple</title>     <link href="http://code.google.com/apis/maps/documentation/javascript/examples/default.css" rel="stylesheet" type="text/css" />
     <input id="address" type="textbox" value="Sydney, NSW">          <input type="button" value="Geocode" onclick="codeAddress()">

                I have tried to use the following code to do this but it does not seem to work.                      // Javascript//            google.maps.event.addListener(marker, 'dragend', function(evt){            document.getElementById('current').innerHTML = '
Marker dropped: Current Lat: ' + evt.latLng.lat().toFixed(3) + ' Current Lng: ' + evt.latLng.lng().toFixed(3) + '
';            });                 google.maps.event.addListener(marker, 'dragstart', function(evt){            document.getElementById('current').innerHTML = '
Currently dragging marker...
';            });           map.setCenter(marker.position);      marker.setMap(map);           //HTML//

Nothing yet...


arcgis python script to tool issue

I have a script to investigate mxd's and save them as lyr symbology/kmz and sld which works fine when run from the .py but does not run when used as a tool from arctoolbox. I have tried arcpy.AddMessage statements to see what's wrong but it just does not work. Please help... SCRIPT               # Author: George Corea, Atherton Tablelands GIS     # georgec@atgis.com.au; info@atgis.com.au          # Licence:Creative Commons          import arcpy, string, datetime, shutil, os     import arcpy.mapping as MAP          #Read input parameters from script tool     MXDList = string.split(arcpy.GetParameterAsText(0), ";")     ProjectPath = arcpy.GetParameterAsText(1)          #MXDList=r'P:\2012\183_TownPlanning_Symbology\Working\TP_Biodiversity.mxd'     #ProjectPath=r'P:\2012\183_TownPlanning_Symbology\Working'          count=0     #Loop through each MXD and print     #for MXDPath in MXDList:     #MXDFile=r'P:\2012\189_Townplanning_scheme\Working\TRC_PlanningScheme_LayoutMaster_v10.mxd'          count=count+1     arcpy.AddMessage('starting...'+str(MXDList))     ##try:     for MXDFile in MXDList:              mxd=arcpy.mapping.MapDocument(MXDFile)         outPath = ProjectPath+'\\'+mxd.filePath[mxd.filePath.rfind('\\')+1:mxd.filePath.rfind('.')]         arcpy.AddMessage('working on 1...'+str(mxd)+str(outPath))         try:             os.mkdir(outPath)         except:             pass         dfList = arcpy.mapping.ListDataFrames(mxd)         arcpy.AddMessage('working on 2...'+str(dfList)+str(outPath))         for df in dfList:         #msd = outPath+'.msd'         #arcpy.mapping.ConvertToMSD(mxd, msd, df, "NORMAL", "NORMAL")         #arcpy.AddMessage(str(count)+'...' +str(outPath)+'...'+'\n')             lyrList=arcpy.mapping.ListLayers(mxd, "", df)             outPath = ProjectPath+'\\'+mxd.filePath[mxd.filePath.rfind('\\')+1:mxd.filePath.rfind('.')]+'\\'+str(df.name)             #os.mkdir(outPath)             arcpy.AddMessage('working on 3...'+str(lyrList)+str(outPath))             for lyrFile in lyrList:                 arcpy.AddMessage(str(lyrFile))                 if lyrFile.isFeatureLayer == True:                     if lyrFile.visible == True:                         arcpy.AddMessage(str(lyrFile)+' is visible')                         #print str(lyrFile)+' is visible...exporting dataset'                         outFileN=str(arcpy.ValidateTableName(lyrFile.longName[lyrFile.longName.rfind('\\')+1:]))                         try:                             arcpy.FeatureClassToFeatureClass_conversion(lyrFile.dataSource, outPath, outFileN)                             arcpy.SaveToLayerFile_management(lyrFile,outPath+'\\'+outFileN+'.lyr', "ABSOLUTE")                             arcpy.LayerToKML_conversion(outPath+'\\'+outFileN+'.lyr', outPath+'\\'+outFileN+'.kmz')                         except:                             errorm=arcpy.GetMessages()                             arcpy.AddMessage('Error...'+str(errorm)+' ...continuing')                           arcpy.AddMessage(str(lyrFile)+' is visible. Not Processing')                 #print str(lyrFile)+' is not visible'         ##except:         ##    errorM=arcpy.GetMessages()         ##    arcpy.AddMessage(str(count)+str(errorM)+'\n continuing...\n')         ##    print errorM         ##    #MAP.PrintMap(MXD, printer)         #Remove variable reference to file         del mxd, msd      ERROR ![enter image description here](http://i.stack.imgur.com/b0vQJ.jpg) and that's it. The text 'working 2...' from the arcpy.AddMessages doesn't work. When run in the Python Interpreter window you get               >>> import os     >>> os.getcwd()     'Q:\\scripts\\py'     >>> MXDFile=r'P:\2011\Job_031_TownPlanning_SeriesProduction\Working\mxd\Nov14_Landslide\18_TownPlanning_B&L_Overlay_Ver4a_Chillagoe.mxd'     >>> mxd=arcpy.mapping.MapDocument(MXDFile)     >>> ProjectPath=r'Q:\scripts\py\junk'     >>> outPath = ProjectPath+'\\'+mxd.filePath[mxd.filePath.rfind('\\')+1:mxd.filePath.rfind('.')]     >>> dfList = arcpy.mapping.ListDataFrames(mxd)     >>> for df in dfList:     ...     lyrList=arcpy.mapping.ListLayers(mxd, "", df)     ...     print lyrList     ...          [<map group layer u'Residences'>, <map layer u'Residential Structure'>, <map layer u'Residential Structure'>, <map group layer u'Base_Top'>, <map layer u'Regional Activity Centre'>, <map layer u'Homesteads'>, <map layer u'State Controlled Road'>, <map layer u'Local Road'>, <map layer u'Local Road'>, <map layer u'Wetland'>, <map layer u'Waterbody'>, <map layer u'Protected Area'>, <map layer u'TRC Boundary'>, <map layer u'Property Boundary'>, <map layer u'Watercourse'>, <map layer u'Tri_RiversTRCClip'>, <map layer u'Rural_DrainageTRCClip'>, <map group layer u'Hazard'>, <map group layer u'NewLandslide'>, <map layer u'100k_trc_tiles_TVM'>, <map layer u'BULLOCK_CREEK'>, <map layer u'BULIMBA'>, <map layer u'BLACKDOWN'>, <map layer u'BELLEVUE'>, <map layer u'MAYTOWN'>, <map layer u'M_SURPRISE'>, <map layer u'LYNDBROOK'>, <map layer u'KIRRAMA'>, <map layer u'KINGFISH'>, <map layer u'HIGHBURY'>, <map layer u'CASHMERE'>, <map layer u'CAIRNS'>, <map layer u'MUNGANA'>, <map layer u'MOUNT'>, <map layer u'MOSSMAN'>, <map layer u'GALLOWAY'>, <map layer u'CONJUBOY'>, <map layer u'CHILLAGOE'>, <map layer u'BARTLE_FRERE'>, <map layer u'ATHERTON'>, <map layer u'ABINGDON_DOWNS'>, <map layer u'WALSH'>, <map layer u'WALLABADAH'>, <map layer u'VALLEY_LAGOONS'>, <map layer u'TULLY'>, <map layer u'STRATHMORE'>, <map layer u'ST_RONANS'>, <map layer u'SOUTH'>, <map layer u'RUMULA'>, <map layer u'RAVENSHOE'>, <map layer u'Landslide Zone'>, <map layer u'High Bushfire Hazard'>, <map layer u'High Bushfire Hazard Buffer '>, <map layer u'Medium Bushfire Hazard'>, <map layer u'Medium Bushfire Hazard Buffer'>, <map group layer u'Base_Bottom'>, <map layer u'Protected Area'>]     [<map group layer u'Residences'>, <map layer u'Residential Structure'>, <map layer u'Residential Structure'>, <map group layer u'Base_Top'>, <map layer u'Regional Activity Centre'>, <map layer u'Homesteads'>, <map layer u'State Controlled Road'>, <map layer u'Local Road'>, <map layer u'Local Road'>, <map layer u'Wetland'>, <map layer u'Waterbody'>, <map layer u'Protected Area'>, <map layer u'TRC Boundary'>, <map layer u'Property Boundary'>, <map layer u'River'>, <map layer u'Tri_RiversTRCClip'>, <map layer u'Rural_DrainageTRCClip'>, <map group layer u'Hazard'>, <map group layer u'NewLandslide'>, <map layer u'100k_trc_tiles_TVM'>, <map layer u'BULLOCK_CREEK'>, <map layer u'BULIMBA'>, <map layer u'BLACKDOWN'>, <map layer u'BELLEVUE'>, <map layer u'MAYTOWN'>, <map layer u'M_SURPRISE'>, <map layer u'LYNDBROOK'>, <map layer u'KIRRAMA'>, <map layer u'KINGFISH'>, <map layer u'HIGHBURY'>, <map layer u'CASHMERE'>, <map layer u'CAIRNS'>, <map layer u'MUNGANA'>, <map layer u'MOUNT'>, <map layer u'MOSSMAN'>, <map layer u'GALLOWAY'>, <map layer u'CONJUBOY'>, <map layer u'CHILLAGOE'>, <map layer u'BARTLE_FRERE'>, <map layer u'ATHERTON'>, <map layer u'ABINGDON_DOWNS'>, <map layer u'WALSH'>, <map layer u'WALLABADAH'>, <map layer u'VALLEY_LAGOONS'>, <map layer u'TULLY'>, <map layer u'STRATHMORE'>, <map layer u'ST_RONANS'>, <map layer u'SOUTH'>, <map layer u'RUMULA'>, <map layer u'RAVENSHOE'>, <map layer u'Landslide Zone'>, <map layer u'High Bushfire Hazard'>, <map layer u'High Bushfire Hazard Buffer '>, <map layer u'Medium Bushfire Hazard'>, <map layer u'Medium Bushfire Hazard Buffer'>, <map group layer u'Base_Bottom'>, <map layer u'Protected Area'>]     [<map layer u'Regional Activity Centre'>, <map layer u'State Controlled Road'>, <map layer u'Frame_Extents2'>, <map layer u'TRC Boundary'>]     >>>       TOOL PARAMS ![enter image description here](http://i.stack.imgur.com/MR1Us.jpg) **=== RESPONSE to Answers/Comments in post ===**               Messages     Executing: InvestigateMXDs P:\2012\096_CCRC_2_Biodiversity\Working\BiodiversityOverlay_70k_Series3b.mxd D:\junk     Start Time: Mon Aug 06 09:01:22 2012     Running script InvestigateMXDs...     Now processing ...       Map: P:\2012\096_CCRC_2_Biodiversity\Working\BiodiversityOverlay_70k_Series3b.mxd      So basically it hangs on the df list. There is no further error message even if you try to capture it. Test on simple mxd is below. Note that this works fine when run as a python script directly. ![enter image description here](http://i.stack.imgur.com/D9yNi.jpg)               Messages     Executing: InvestigateMXDs P:\2012\183_TownPlanning_Symbology\Working\test.mxd P:\2012\183_TownPlanning_Symbology\Working\junk     Start Time: Tue Aug 07 09:11:35 2012     Running script InvestigateMXDs...     Now processing ...       Map: P:\2012\183_TownPlanning_Symbology\Working\test.mxd      When run as a python we get the following in less than 30s > > >     starting...['test.mxd']     working on 1...<geoprocessing Map object object at 0x029471E0>P:\2012\183_TownPlanning_Symbology\Working\test          working on 2...[<DataFrame object at 0xee12950[0xee22fc0]>]P:\2012\183_TownPlanning_Symbology\Working\test     working on 3...in: P:\2012\183_TownPlanning_Symbology\Working\test\Layers     TRC_Boundary is visible PROCESSING...     TRC_Boundary is not visible. Not Processing     Map_Extents is visible PROCESSING...     Map_Extents is not visible. Not Processing      and directory structure ![enter image description here](http://i.stack.imgur.com/D9yNi.jpg) **##Issue with modified script from Polygeo** Same problems. ![enter image description here](http://i.stack.imgur.com/id0jH.jpg) **=== END Response ===**


new administrator - Advice to choise a best web sharing GIS DATA with authorization information level

i'm a new "GIS 's needer" and i have just Pleasantly surprised to find your community ( thanks google!!!! ) I need for my office's project a appication who permit web sharing(only reading)some gps point on the satellite's map of my country(like google maps) but also this solution will permit to manage users ...i want just say that there 'll be differents users to read gps map and i want lanage the level of authorization to see information of each point e.g :when administrator"i" update the map with GPS coordonate of "A" point with some information to represent a location where somebody do a business trip I would have wanted that my boss have permission to see all information of the business trip ,but that my office colleague hav just authorization to see the name of man who do business trip... please sorry for my english i live in french's country and i hope your advice to do my best for this project thank you !!!!


Render data on top of Google Maps

I have an interpolated temperature data samples that covers a map evenly. The data is in the form of (lat,lon,temp). I want to render that data on top of google maps. I tried google heat maps and the rendering is inaccurate. I assume one way to do it is to have a layer stored in geoserver. I don't want to use this option. I would prefer if the data are included with the webpage and rendered at runtime. Kindly advice what technique to follow. yours sincerely **Update** t seems that is is possible to render polygons in google maps. You could set the opacity, and color. What I did is interpolating temperature at different degrees, and then creating a large number of polygons that get rendered on top of the map. The color of each polygon is changed based on the approximated temperature value. I got the following results: ![Rendering large number of polygons on top of google maps](http://i.stack.imgur.com/6AOJ5.jpg) Now I have the following isse: File size is too large, about 1MB is needed to store the code of java script. How can I send the java script code in compact format to the client? If anyone knows how to compact the data then please let me know. Also let me know if there is an alternative way to render polygons on top of google maps. thank you.


Popup is stationary does not move with the location

I am using the GeoExt example and used the same code in my application. This is my code The problem I am facing is that the Popup box in my application on `http://128.196.142.94/geo/test/test.html` is stationary. Where ever you may click, the popup simply does not move. Not sure where I am going wrong?


Create polygon layer from polyline layer

Our systems expect polygon shapefiles, representing a collection of non- intersecting and geometrically valid (multi)polygons. Specifically, farm fields. We often receive our representations of farms and their fields as polyline shapefiles containing the field boundary polylines, typically with some small geocoding errors like dangling ends. Is there a reliable way to transform the polyline layer into a polygon layer, where the polygons perimeters are calculated as minimal loops around polylines in the shapefile? I understand I'll probably have to tidy up the polylines first. I am primarily interested in solutions involving Quantum GIS or other free tools, since my company doesn't buy any of the professional packages.


Find Duplicate Feature Classes in ESRI Geodatabase

I have inherited a SDE instance that has 1000+ feature classes within which are around 100 Feature Datasets. I have found through manual exploration that there are at least some Feature Classes that have been duplicated across the different Feature Datasets. Though the names are different everything else is the same. Is there a way to find and possibly catalogue these duplicate Feature Classes without searching manually? Metadata is very poor so this neither a solution nor a way to even to narrow down my search. SDE and Desktop are version 10.0


how to use arcsde need from basics

i need to access tables which are in oracle database by arccatalogue. i am new to this arcmap software, so plz suggest me detailed steps to access database tables from acrmap. thank you.


Is there a way to view KML natively in ArcMap?

Answering this question: How do I import/export KML to ArcGIS Desktop? raised another one. There are loads of tools to import/export kml/kmz. Are there tools available for native viewing of KML in ArcMap? I'm curious if anyone has built a custom featurelayer(s) that reads kml, similar to the way ZigGIS uses PostGIS as a data source. Or maybe something like inmemoryworkspacefactory, which probably doesn't count as native, but would have some advantages: http://geobabble.wordpress.com/2007/07/10/consuming-georss-in-arcmap-with- inmemoryworkspacefactory/