title: OpenLayers WFS flip coordinates

text: I have a problem integrating WMS and WFS. Basically what i want is when user clicks on WMS then a clicked polygon is added to vector (and selected). The problem is that the WMS is in version 1.3.0 and WFS 1.1.0. OpenLayers automatically flips axis for WMS but when it calculates bbox for wfs request, it does not. I have the following code:                                  var map = new OpenLayers.Map({numZoomLevels: 20});                       var wms = new OpenLayers.Layer.WMS(                      "OpenLayers WMS",                      "http://sdi.geoportal.gov.pl/wms_dzkat/wmservice.aspx?",                      {layers: ['Dzialki','Numery_dzialek'], version:'1.3.0'}                  );                  select = new OpenLayers.Layer.Vector("Selection", {styleMap:                       new OpenLayers.Style(OpenLayers.Feature.Vector.style["select"])                  });              var protocol = new OpenLayers.Protocol.WFS({                      version: "1.1.0",                      //url:  "http://sdi.geoportal.gov.pl/WFS_DzKat/service.svc/post",                      url:'/GSIP/proxy/dopost/',                      featureType: "Dzialki",                      featureNS:'http://www.intergraph.com/geomedia/gml',                      featurePrefix:'gmgml',                      geometryName:"GEOMETRIA_DEG"                  });                       control = new OpenLayers.Control.GetFeature({                      protocol: protocol,                      box: true,                      //hover: true,                      multipleKey: "shiftKey",                      toggleKey: "ctrlKey"                  });               control.events.register("featureselected", this, function(e) {                      select.addFeatures([e.feature]);                  });                  control.events.register("featureunselected", this, function(e) {                      select.removeFeatures([e.feature]);                  });               map.addControl(control);                  control.activate();           And generated request after click:               <wfs:GetFeature xmlns:wfs="http://www.opengis.net/wfs" service="WFS" version="1.1.0" maxFeatures="10" xsi:schemaLocation="http://www.opengis.net/wfs http://schemas.opengis.net/wfs/1.1.0/wfs.xsd" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance">     <wfs:Query typeName="gmgml:Dzialki" srsName="EPSG:4326" xmlns:gmgml="http://www.intergraph.com/geomedia/gml">     <ogc:Filter xmlns:ogc="http://www.opengis.net/ogc">     <ogc:BBOX>     <ogc:PropertyName>GEOMETRIA_DEG</ogc:PropertyName>     <gml:Envelope xmlns:gml="http://www.opengis.net/gml" srsName="EPSG:4326">     <gml:lowerCorner>19.143386078712 52.371721634254</gml:lowerCorner>     <gml:upperCorner>19.143439722893 52.371775278435</gml:upperCorner>     </gml:Envelope>     </ogc:BBOX>     </ogc:Filter>     </wfs:Query>     </wfs:GetFeature>      The above request does not work (0 features returned) but when i flip coords so it looks like:               <gml:lowerCorner>52.371721634254 19.143386078712</gml:lowerCorner>         <gml:upperCorner>52.371775278435 19.143439722893</gml:upperCorner>      it works. How to force OpenLayers to build getFeature request with flipped coords? In WMS Layer there is attribute reverseAxisOrder but can't find it in Vector Layer


title: Moving data from columns in a table to rows by a key field?

text: I would like to shift data found in a column to a row by Key field. I have a table, within a geodatabase, set up similar to what's below.               KEY_    Value     123     Yellow     123     Blue     456     Red     456     Yellow     789     Green     246     Green     246     Orange      I want to create a table that moves the data in the Value field to rows based on the KEY_ field. The result will look like what's below.               KEY_    Value001    Value002     123     Yellow      Blue     456     Red         Yellow     789     Green     246     Green       Orange      Does anyone have an idea of how this can be done? I tried using the Transpose tool Esri offers but it didn't produce the result I'm looking for. I have basic understanding of Python.


title: How to get all OIDs from a table?

text: I have table which i need to get all the OID's from it. I tried this code :               List<int> OIDList = new List<int>();     IQueryFilter queryFilter = new QueryFilterClass();     queryFilter.SubFields = "OBJECTID";          for (int i = 1; i <= table.RowCount(queryFilter); i++) {         OIDList.Add(i);     }      But the problem is there is some values were skipped, e.g. the table contains 60 row but the OID for the last row is 100. in the code above, i can not catch any OID more than 60. How can i get all the OIDs?


title: GeoServer REST API Issue: PUT HTTP Method

text: I've the function below which creates or replaces an existing layer group:               public bool AddToLayerGroup(string layerGroupName, string layerName)             {                 try                 {                     // Retrive the layer group                     string cLayerGroup = GetLayerGroup(layerGroupName);                     // Add the new layer as a new layer node among the retrieved layer nodes                     XmlDocument newLayerGrpXml = AddLayerNodeToLayerGroupXml(layerName, cLayerGroup);                     // Convert the newly created layer group xml to string                     string newLayerGrpString = ConvertXmlToString(newLayerGrpXml);                     newLayerGrpString = newLayerGrpString.Replace("<?xml version=\"1.0\" encoding=\"utf-16\"?>", "");                          // Delete the existing layer group                     //DeleteLayerGroup(layerGroupName);                          // Now, do ur thing here                     string gUrl = GEOSERVER_HOST + "/rest/layergroups/" + layerGroupName + ".xml";                     WebRequest request = WebRequest.Create(gUrl);                          request.ContentType = "application/xml";                     request.Method = "PUT";                     request.Credentials = new NetworkCredential(GEOSERVER_USER, GEOSERVER_PASSWD);                          byte[] buffer = Encoding.GetEncoding("UTF-8").GetBytes(newLayerGrpString);                     Stream requestStream = request.GetRequestStream();                     requestStream.Write(buffer, 0, buffer.Length);                     requestStream.Close();                          WebResponse response = request.GetResponse();                          return true;                 }                 catch (Exception e)                 {                     Console.Write(e.Message);                     return false;                 }             }      This executes successfully. Tomcat logs also indicates successful execution of the PUT method. However, I can't see the created layer group using GeoServer web UI. Is there anything else I need to do after executing the above code? Trying to create another layer group with the same name fails with the error `layer group already exists`; a further proof that the layer group was created. But, why can't I see it using the GeoServer web UI. Any assistance will be highly appreciated.


title: How to access feature classes in file geodatabases with Python and GDAL?

text: I am trying to access a vector dataset in a ESRI File Geodatabase using Python + GDAL. I have successfully compiled GDAL with the file geodatabase API. The FileGDB driver is working correctly since entering               ogrinfo --formats      shows the FileGDB driver and entering               ogrinfo myfilegdb.gdb       gives me the correct information on the content of the database. However, I can't find out how to access the content itself in Python. For accessing a shapefile, I would write:               driver = ogr.GetDriverByName('ESRI Shapefile')     ds = driver.Open('shapefile.shp', 0)      When accessing a FileGDB feature class I would assume using the commands:               driver = ogr.GetDriverByName('FileGDB')     ds = driver.Open('myfilegdb.gdb/feature_class', 0)      but this doesn't seem to work since it can not identify/locate the data set. Does anyone know how to call individual feature classes from a ESRI FileGDB. I'm using Python 2.7, GDAL 1.9.1, filegdb api 1.2 on Ubuntu 12.04 x64. Thanks for any suggestions!


title: How to catch PostgreSQL/PostGIS errors in GDAL/OGR Python bindings?

text: Is it possible to catch PostgreSQL output error messages in OGR after a query has failed? For example               conn.ExecuteSQL('SELECT * FROM non_existing_table;')      In Postgres this returns an error message:               "Error: relation 'non_existing_table' does not exist."      In ogr we get the standard:               "Error: current transaction is aborted, commands ignored until end of transaction block"      It would be a bonus if we could catch/parse the Postgres error message. I suspect this is a bit of a long-shot given what I've read about SWIG/CPL, but worth asking to see if anyone's got ideas. Ideally this would be in Python, but I'll take C/C++ comments as well! Thanks.


title: Concurrent users using a DLL to access ArcObjects?

text: So currently I have some ASP.net code where users are updated GIS through a DLL file that uses ArcObjects to talk to the MapService and layer. My question is, will there be concurrent user issues? What if 10 people are submitting updates to the same items that need updating in ArcGIS map service? Will the DLL freeze up? Will ArcGIS service freeze up? Or does it handle concurrent users in a Queue-style via the ArcGIS service, the DLL, or IIS itself maybe realizes the DLL is in use?


title: How to get GeoJSON features using OpenLayers?

text: Can anybody explain how to get geojson features .Can somebody show way to get attributes from geojson text.               OpenLayers.ProxyHost = "/gt/proxy.cgi?url=";          var lon = 5;         var lat = 44;         var zoom = 5;         var map, layer;              //function init(){            var bounds = new OpenLayers.Bounds(                     -125051.196, 3967670.8507,                     1549596.1896, 4677667.0145                 );                 var options = {                     controls: [],                     maxExtent: bounds,                     maxResolution: 6541.59135,                     projection: "EPSG:23036",                     units: 'm'                 };         map = new OpenLayers.Map('map', options);         map.addControl(new OpenLayers.Control.PanZoomBar({                     position: new OpenLayers.Pixel(2,10)                 }));                      map.addControl(new OpenLayers.Control.Navigation());                 //map.addControl(new OpenLayers.Control.Scale($('scale')));                 map.addControl(new OpenLayers.Control.MousePosition({element: $('location')}));                  var geojson_format = new OpenLayers.Format.GeoJSON();             var vl2 = new OpenLayers.Layer.Vector("",             {isBaseLayer: true,              styleMap: new OpenLayers.StyleMap({'default':{                     strokeColor: "#00FF00",                     strokeOpacity: 1,                     strokeWidth: 1,                     fillColor: "#FF5500",                     fillOpacity: 0.5,                     label : "${ad}",                     fontSize: "8px",                     fontFamily: "Courier New, monospace",                      labelXOffset: "0.5",                     labelYOffset: "0.5"                 }})             });             map.addLayer(vl2);                 function handler(request) {                  vl2.addFeatures(geojson_format.read(request.responseText))             map.zoomToExtent(bounds);             }             var request = OpenLayers.Request.GET({             url: "myurljson&srsName=EPSG:23036",             params: {},             callback: handler         });      please show me some way


title: How to change position of columns in attribute table?

text: If I have an attribute table with columns A B C D, is it possible to move the columns to different positions, say, B C A D? When making a new column with the field calculator it would be useful to put it somewhere other than last position or to be able to move key columns close to it.


title: Looking for a developer-friendly web gis

text: I'm a recent Computer Science graduate who prides himself in his knowledge of software development and problem solving, but has no idea how gis systems work. (Read: **I'm a programmer** ) I've been charged with developing (ie. coding) a new **webmap** from scratch using software called GeoMedia Professional Web Maps. I am attempting to get support on developing the application, however it has been like trying to find gold in a pitch dark cave. I swing my axe wildly and sometimes get some gold dust (it's all trial and error, because of the lack of decent documentation). So barring that my attempts to use this software succeed, what would be a very **developer friendly** GIS software I could use? I would prefer anything that does **not** lock me into software (ie. Prefer not Microsoft), however if people have found one GIS in C#/asp.net happened to be the most intuitive, I would bite the bullet. (Preference: PHP, C++, or Python) [Btw I do have a GIS professional working with me on the project, but they know nothing about programming] # A little more background: The reason I'm thinking I may need to find new software, is that this GeoMedia stuff has all of their documented code as ASP/VBScript and only shows examples using access. I refuse to use ASP/VBScript (I already did some edits in it and I wanted to throw up as I was developing) and our data appears to be in oracle and in Geomedia's proprietary smart store DB. Surprisingly, Geomedia's documentation doesn't cover how to connect to there very own databases. I tried to search for books, or online help, or documentation and found nothing. When I called there tech support, they referred me to the local help files and the demo, which I already looked at, which again uses access. My only saving grace for this software is that it was written using COM objects, which means I can program in any language that has bindings to COM. So I've been writing in PHP, but still the COM objects returns errors that I'm not using them correctly, and the errors aren't descriptive enough to give me tips on troubleshooting the problems (The specific objects have issues, not using COM... I've gotten a bunch of the GeoMedia COM objects to work in PHP). EDIT: Although unimportant to the question, I found out COM exceptions could be tapped into to get better formatted information using exception. So that last point is a little less correct [still the errors could be more informative].


title: How to display scale on ArcGis-Engine map?

text: I embedded ArcEngine into a c#/wpf application and use it to display a map. Now for print my client needs to display the scale of the map like it is used on street maps (the small bar which shows how long 500m are, and/or the scale information like 1:25000). Prints are always on a set paper size like A4. Is there a simple way to achieve this in an ArcEngine application?


title: Metadata for Bing Maps Aerial in ArcGIS

text: Is there a way to dynamically pull up metadata on the Bing Maps Aerial web mapping service in ArcGIS? I am working at a 1:245 scale across various locations in OR and CA and need to know if Bing's coordinate system is accurate to within 0.5 meters. Additionally, I need to know the resolution. ESRI's description is pretty vague: > This web map contains the Bing Maps aerial imagery web mapping service, > which offers worldwide orthographic aerial and satellite imagery. Coverage > varies by region, with the most detailed coverage in the USA and United > Kingdom. Coverage in different areas within a country also varies in detail > based on the availability of imagery for that region. Bing Maps is > continuously adding imagery in new areas and updating coverage in areas of > existing coverage.


title: (New at Python Scripting) Creating a Loop so i can replace text in page layout in multiple mxds instead of one mxd at a time

text: I found some code on a website that will perform a search and replace on page layout text elements. Currently you can only do it one MXD at a time. I was hoping to change the code so I can do a search and replace on a folder with a bunch of MXDs in it. The script is below. I know I need a loop, just not sure how to do it; never programmed before.               import arcpy, string, os           #Read input parameters from script tool     Path = arcpy.GetParameterAsText(0)     oldText = arcpy.GetParameterAsText(1)     newText = arcpy.GetParameterAsText(2)     case = arcpy.GetParameter(3)     exact = arcpy.GetParameter(4)     outputMXD = arcpy.GetParameterAsText(5)          try:         #Referent the map document         mxd = arcpy.mapping.MapDocument(mxdPath)                       #Find all page layout text elements         for elm in arcpy.mapping.ListLayoutElements(mxd, "TEXT_ELEMENT"):                  if exact:                 if case:                     if oldText == elm.text:                         elmText = elm.text.replace(oldText, newText)                         elm.text = elmText                 else:                     if oldText.upper() == elm.text.upper():                         elmText = elm.text.upper().replace(oldText, newText)                         elm.text = elmText                else:                 if case:                     if oldText in elm.text:                         elmText = elm.text.replace(oldText, newText)                         elm.text = elmText                 else:                     if oldText.upper() in elm.text.upper():                         elmText = elm.text.upper().replace(oldText, newText)                         elm.text = elmText                           mxd.saveACopy(outputMXD)              del mxd          except Exception, e:         import traceback         map(arcpy.AddError, traceback.format_exc().split("\n"))         arcpy.AddError(str(e))


title: How do I get the count of NoData values in a raster?

text: I have created a raster in ArcGIS 10 with a large number of holes (NoData) in it. I want to fill the holes using something like the the following:               rOut = arcpy.sa.Con(arcpy.sa.IsNull(rGround),     arcpy.sa.FocalStatistics(rGround, arcpy.sa.NbrCircle(FillRadius,"CELL"), "MEAN"),     rGround)      It is desireable to fill the holes by averaging across a circle with a relatively small radius. But a small radius leaves larger holes unfilled, so I need to iterate this operation. To do that I would like to be able to get the number of NoData values in my raster, once it reaches zero I can stop filling. regards, Dave


title: OpenLayers behaves problematic after certain zoom level in Google Maps Satellite view

text: I have some polygon data in `.shp` format to show in Google Maps Satellite view. I can show them well until certain zoom level. ![enter image description here](http://i.stack.imgur.com/NVExJ.png) But after certain level (zoom level 17 or 18), `OpenLayers` layer and `Google Maps Satellite` layer displace each other like this. ![enter image description here](http://i.stack.imgur.com/rDwz6.jpg) We know `Google Maps Satellite` View changes it camera angle after certain zoom level and put it to 45 degrees. How can i resolve this with `OpenLayers`.


title: Is there a common interface for ST_Geometry SDO_Geomtry

text: I've worked a fair bit with Oracle Spatial and I've had some requests to transition to PostgreSQL Spatial. Is there a common interface for SDO_Geometry and ST_Geometry? I'm working in Java with EclipseLink ORM.


title: Converting JSON/GeoJSON files to Shapefiles

text: I have a quick question please about converting JSON/GeoJSON files. I'm trying to convert some JSON files into something I know how to work with (shapefile, R objects or PostGIS). I'm afraid I have no experience with JSON. I've tried converting them to shapefiles using ogr2ogr, readOGR in R, and FME. I've tried converting them and uploading them to a PostGIS database using FME. And I've tried converting them to R objects using the rjson package. Nothing works. ogrinfo can't read them and the couple of online JSON/GeoJSON validation tools I've tried say they are not valid JSON files? Would someone who understands JSON/GeoJSON mind taking a very quick look at one of the files please, and tell me if they're OK? If they're fine, that's great, I can go away and try and work it out. If not, I can save myself a little head scratching. I'm sure there's probably nothing wrong with them. Here's a link to one of the files in my Google Drive: Opendata Scotland JSON File


title: ArcGIS Geoprocessing Service fails on data loading

text: I have a Desktop Script tool that works fine when run on desktop, but does not work when the results are turned into a geoprocessing serice. Running ArcGIS 10.2 on Desktop and 10.2 on server. I have found that the geoprocessing service fails at the line of code where it tries to to load a data file I have saved to file. I have tried saving the data files using pickle, json, and also array files as .npy and .npz. It doesn't seem to make any difference. The geoprocessing service just will not load any array type file I provide. I have tried letting the geoprocessing packager save the file to the server, I have saved the files to a registered folder, and it doesn't seem to matter. Even text files will not load. Has anyone been successful at loading a numpy array, dictionary, or similar type file into a geoprocessing service on arcgis server? If so, how did you do it? I added error handling to my script, and here are the messages when I run the Geoprocessing Service:               Traceback Info:       File "E:\ARCGIS~2\DIRECT~1\ARCGIS~4\ARCGIS~2\Grace\HT5C33~1.GPS\EXTRAC~1\v101\HTAT_S~1\140307_HTAT_debug.py, line 74, in [module]      npzfile = numpy.load(array_package)            Error Info:       [Erno 2] No such file or directory: u'C:\\Users\\arcgis\\AppData\\Local\\Temp\\..\\cd\\tooldata\\package131123.npz      The E: drive is on the gis server machine and the .npz file is my data file I am trying to load. 140307_HTAT_debug.py is my source script. Apparently the npz file was copied over to the server, but not sure what that path means with the double dots, and why the server can't read it.


title: Using Python to connect to Multiple SDE Database Connections

text: I am writing a script that accesses multiple sde Database connections. I have set the script to include accessing the environment workspace:               import arcpy          arcpy.GetParameterAsText(0)     arcpy.GetParameterAsText(1)          arcpy.env.workspace = r"Database Connections\somedb.sde"          do some python stuff...          arcpy.env.workspace = r"Database Connections\someOtherdb.sde"      So, I open then tool, enter my parameters and run it. When it gets to setting the workspace, I get the regular popup window to prompt me to enter my username and password to the sde database. The problem I have it that the tool GUI runs on top of the popup. Any time I try to click on the popup window to enter my username and password, the GUI automatically pops on top again and I can't type anything. Is there any way for me to build the GUI so that I can enter my username and password aside from using the `arcpy.CreateArcSDEConnectionFile_management` function? My connection files already exist, so I don't need to create new ones. I just need to log in to it.


title: Why are Images used as label not rendered correctly by GeoServer?

text: When I try to render the label (png image) through GeoServer as an image, its quality is getting affected i.e. the image is not rendered sharply. Are there any settings that need to done in GeoServer or the SLD file to correct this?


title: For the non-IT savvy GIS person (not great w/ Linux & Apache) could someone suggest how to set up PostgreSQL, Geoserver, and ISS/Apache?

text: I've heard of GeoCommon as being an easy GIS solution, but I'd like to host the data locally for internal users. I know of a group called Open Geosuite - but the ongoing maintenance is something I'd like to avoid since I can do pretty well with maintenance and content management & web site publishing after the initial set-up. Does anyone know of an out-of-the-box solution or a group/company that can help with initial set up for these types of map service stacks? Thanks for any advice you can give for a noob!


title: Fantasy map, how do i slice large png and georefrence it the right way?

text: Im and IOS developer working on a map based ios game, and new to GIS stuff. I've read a lot of info and i started to get a hint about what i need to do in order to meet my goals. I have a png map i created, 10000 px width and height. I want to slice it to 256X256 px tiles, and add a pgw file to each slice to indicate it's location in the world. The map WILL grow over time as more players enters the game, so i need to be able to add more "space" dynamically. the png tiles are also dynamic and will be chaned during playgame (not the location, but the contents). i've successfully sliced the png with gdal2tiles but no pgw file is created, how can i do this? at the end i want to have set of directories, ex (tile_1), that contains: 1. the png file 2. the pgw file 3. some more files i needed internally again each tile must be in it's own directory. after i've done this, i want to serve the map with GEOSERVER, and show it on the iPhone/iPad with the mapbox ios framework. How can i slice and georefrenced the images? what projection system should i use? any more helpful data? (i'm really new so any details explanations will be welcome). Thanks alot in advanced. Or


title: Why does ST_Transform fail?

text: I got a postGIS error like this,               ERROR:  transform: couldn't project point (100.496 13.7118 0): latitude or longitude exceeded limits (-14)     SELECT ST_Length(ST_Transform(ST_GeomFromText(E'LINESTRING(100.495995129 13.7117836894,100.495962221169 13.7117761471941)',4326),27700));      This lat/lng comes from Bangok, Thailand and the 27700 is a british national grid. Do you think this casued the problem? May I mention that problems like this did not occur when I got lat/lng from Spain or USA. Any input would be much appreciated. Kind regards, A


title: Given a list of surveyor points, what is the algorithm used to simplify the points and retain the terrain characteristics?

text: Given a terrain and surveyor points sampled from from it or from its contours, is there any algorithm that one can use to simplify the points ( i.e., reduce the surveyor points number) and at the same time retaining the terrain characteristics? The reason I ask this is because the surveyors may take too much redundant points at flat terrain, and these points serve nothing but to slow down my volume calculation process, so I want them to get filtered out.


title: What are my options for creating a tile cache via a WMS server which requires http authentication?

text: I'd like to create a tile cache and offer the tiles via TMS/WMS. The source is a WMS server that is protected by basic http authentication, but no SSL. I'll also need to reproject the tiles. This is for my own personal use. I'd prefer something that can run on either Windows or Linux, and open source is a huge plus. From some searching it appears MapServer can access WMS servers requiring authentication, but doesn't cache. GeoServer can cache WMS, but can't use authentication. Lastly, TileCache can't reproject. What are my options?


title: What is a suitable Global/Universal Unique Identifier for a PostGIS database?

text: I have read that using OIDs as a primary key in a postgreSQL/PostGIS db is poor practice because there are instances where these can be reset. Sounds logical, but then what is a suitable alternative? I believe there is an option to use a "Universal Unique Identifer" UUID, but the large text and number value that spits out is horrible. Just a bit more background to my situation. I have all of my spatial tables created with a field called "gid" which is the primary key for that table and unique only to that table. I have an issue now because I want to relate my spatial tables (all with a "gid" field starting at 1 and incrementing) to one large table with the related information. Obviously for my relationship to work all of my spatial features need a unique identifier which differentiates them from one another. _EDITED_ Added this image as per Peters' comment. Peter this is the idea I have in my head, it may not be the best way to go about it or it might not even be good db design. I am interested in what you think. ![Conceptual diagram](http://i.stack.imgur.com/VjexJ.png) Any tips? Thanks Ando


title: How to migrate an extension from ArcGIS 9.3.1 to 10?

text: I was given an ArcGIS 9.3.1 extension built in Visual Studio (VS) that needs to be updated to work in ArcGIS 10. When I open the VS project on my computer (that has ArcGIS 10 installed on it) I get a few library reference errors. I attempted to remove the references then add them back in thinking that would resolve the error, however the errors are still present. Why are these references not being viewed by VS? Here are the errors: **Update** I changed the Target .NET Framework from 2.0 to 3.5 and added in the missing file and assemblies via the MSSetup>Primary output from HT.MSExtension(Active). This removed all the errors except error 4. Do I need to re-source ESRI.ArcGIS.ADF.dll to ESRI.ArcGIS.ADF.Local.dll? **Update2** I updated ADF to ADF.Local, however I'm still getting error 4. Any ideas? > Error 2 Unable to find source file 'C:\1_Projects\MS\Source\MS_ArcGIS- > Extension\MSSetup\stdole.dll' for assembly 'stdole.dll', located in > '[TARGETDIR]' C:\1_Projects\MS\Source\MS_ArcGIS- > Extension\MSSetup\MSSetup.vdproj MSSetup > > Error 1 Unable to find source file 'C:\1_Projects\MS\Source\MS_ArcGIS- > Extension\MSSetup\ESRI.ArcGIS.Editor.dll' for assembly > 'ESRI.ArcGIS.Editor.dll', located in '[TARGETDIR]' > C:\1_Projects\MS\Source\MS_ArcGIS-Extension\MSSetup\MSSetup.vdproj MSSetup > > Error 4 Could not load file or assembly 'file:///C:/Program Files > (x86)/ArcGIS/DeveloperKit10.0/DotNet/ESRI.ArcGIS.ADF.dll' or one of its > dependencies. An attempt was made to load a program with an incorrect > format. Line 195, position 5. C:\1_Projects\MS\Source\MS_ArcGIS- > Extension\HT.MSExtension\Forms\msManager.resx 195 5 HT.MSExtension > > Error 3 Unable to find source file 'C:\1_Projects\MS\Source\MS_ArcGIS- > Extension\MS.mxd' for file MS.mxd', located in '[TARGETDIR]', the file may > be absent or locked. C:\1_Projects\MS\Source\MS_ArcGIS- > Extension\MSSetup\MSSetup.vdproj MSSetup Thank you


title: How to calculate distances between several points along network?

text: I have several points (bicycle share hubs), and I want to calculate the distance between each of them using a network. I want a CSV/XLS that has three headers: hub_id1, hub_id2, distance. I've created a network dataset (after consultation), but how do I get the distance between each of the (blue) hubs pictured below? I know a dollop of Python, so if this requires some scripting, I'm comfortable trying to apply what I know. I'll be the first to say, though, that I've not scripted with spatial data before. FYI, both of these posts ask a similar question, but not using a network grid. Would love some input! Thanks! ![Bicycle hubs and the street grid in ArcGIS.](http://i.stack.imgur.com/IQ58H.png)


title: What kind of QT widget is the QGIS drawing area?

text: I mean when open QGIS , and add a vector layer or other type of data. The content will show on the "Drawing Area", i want to ask, the drwaing area is some type of widget or just in the QMainWindow? ![enter image description here](http://i.stack.imgur.com/7ICjS.png) I download the qgis source code(1.8 or 1.9). and searche the src/ui/ folder, do not find any named similar to "drawing " *.ui file. And i use qt-creator and open the "qgisapp.ui" and judge by the size of the window, the drawing area is in this window. But i am not sure, because i am not familiar with qgis code. ![enter image description here](http://i.stack.imgur.com/cGlaG.png) What i want to do is want to mime QGIS to code little GIS drawing toy using Qt. So anyone could answer my question?


title: can't Find file chooser widget in QTdesigner for QGIS file name edit widget custom form

text: I designed a custom form (*.ui) for my shape file under QT Designer. I want the user to choose a file through a QGIS "file name" widget but I can't find the control in QT designer. Regards,


title: Are there any alternatives to ArcGIS Locator Report?

text: I'm a pretty new GIS programmer and I've found ArcGIS/Arcpy to be extremely unstable. I'm trying to run a Locator Report in arcpy with no success. I'm wondering if there are any other alternatives? Are there other software packages that can accomplish the same thing but with open source code in case I encounter bugs?


title: Which software to use to create professional fly-over/ fly-through animations?

text: I will need to create a high quality fly-over animation which penetrates the surfaces and reveals subsurface geology. The intended use is for demonstrations, presentations and a company website so it cannot include copyrighted data. _**Available Data:_** I have high resolution satellite & aerial imagery, 20m resolution DEM for the entire area of interest, 10 x 10 km 1m resolution LiDAR data, subsurface geology in 3D DXF and ArcGIS multipatch features; basically all the data I would like to "fly over and through". _**Software:_** I have ArcGIS 10 ArcEdtor with Maplex, 3D and Spatial Analyst, Access to MapInfo Professional (no experience), Google Sketchup Pro, Leapfrog Mining and CorelDraw X5 and Adobe Master Collection CS5 (limited experience) at my disposal. I've seen some ArcGlobe and Google Earth animations and from what I have seen seems too choppy and the camera moves too fast and turns abruptly. Here is an example of something I would like to create except that I need to also explore features beneath the surface. I will need some pop-up balloons to indicate points of interest similar to what's in the linked video as well as a compass or an arrow showing approximate north orientation. I am only looking for suggestions how to create the "guts" of this animation; the actual fly over and into the data with the pop up labels and direction. I would really appreciate any suggestions on:   * Recommended workflow for such animations (preferably using the tools I already have)   * Other software commonly used to create such animations Please share links to any good examples.


title: Failed creating a feature class inside a file geodatabase

text: I'm having trouble creating a feature class within a new file geodatabase. I can create shapefiles and sde databases, but not feature classes. Here's the message I get from arcmap 10.1. Any ideas? ![enter image description here](http://i.stack.imgur.com/XFQWk.png) In Catalog on a network share I create a file GDB by right clicking and go to new > File Geodatabase. That works fine. Then I click on the newly created File GDB and New > Feature Class. This is where it gives me the error message. When I recreate this process on my desktop it works fine. It seems there's something wrong with the network share/privileges.


title: Mapserver PHP mapscript: how to set attribute value of layer?

text: **EDIT** : solved, see working code below How can I set attribute values of a layer object with MapScript? I'm using MapServer version 5.0.3, the php_mapscript.so lib is loaded in my app context I have a simple mapfile:               MAP         SIZE 500 500         EXTENT -180 -90 180 90         UNITS dd         CONFIG "CPL_DEBUG" "ON"         CONFIG      "MS_ERRORFILE" "/var/www/Data/Data/mapserver/mapserver.log"         SYMBOLSET "/var/www/Data/Data/mapserver/symbols.sym"         FONTSET "/var/www/Data/Data/mapserver/fontset.txt"          #--------------------------------------------------------------------------------     # load lines     #--------------------------------------------------------------------------------              LAYER             NAME "loadlines"             TYPE POLYGON             DATA "/var/www/Data/Data/mapserver/data/loadlines/loadlines.shp"             STATUS ON             OPACITY 100             DEBUG 5                  CLASS #default coloring                 STYLE                   COLOR 0 0 0                   WIDTH 6                   #PATTERN 40 10 END                 END                 NAME "gust"                 OUTLINECOLOR 254 254 254                 COLOR 100 100 100                 #OPACITY 50             END                  TEMPLATE "blank.html"             DUMP TRUE # allow GML export             PROJECTION                 "init=epsg:4326"             END             METADATA                  "wms_title"     "ffg_test"                      "wms_srs"       "EPSG:54004 EPSG:4269 EPSG:4326"              END          END      END      Which I want to view by calling this php script:               <?php     $map_path="/var/www/Data/Data/mapserver/routeguardTemp/test/";          $map = ms_newMapObj($map_path."loadlines.map");     $layer = $map->getLayerByName('loadlines');     $layer->set('opacity', 10);     $image=$map->draw();     $image_url=$image->saveWebImage();          ?>           <HTML>       <HEAD>           <TITLE>Map</TITLE>       </HEAD>       <BODY>           <IMG SRC=<?php echo $image_url; ?> >       </BODY>      </HTML>      The script should alter the opacity value of the layer, but there is no effect. How can I do this? Thank you.


title: Figuring out what coordinates supposedly using 'NAD31 projection' refer to

text: I'm working with a GeoJSON file that contains sets of coordinates associated with street corners in the city of Montreal, Canada. The spec sheet of the API delivering those GeoJSON specifies that the coordinate system uses the _NAD31_ projection datum, which in spite of my search I doubt exist. Because of the street names associated to the coordinates in the GeoJSON, I am able to deduce that the erroneous _NAD31_ coordinates `300247.509025, 5040897.253` refer to the _WGS84_ latitude/longitude `45.5079942,-73.5582652`. Does anybody know what might the geodetic datum _NAD31_ refer to, and how I can convert it to _WGS84_ coordinates?


title: Cell Statistics in QGIS

text: Cell Statistics in ArcGIS will calculate the mean (or statistic of choice) among multiple rasters on a cell-by-call basis. Can QGIS do this? If so, how?


title: Looking for an API that will autocomplete adresses, and split them up into parts

text: I'm kinda new to the whole GIS thing, and found that is necessary for onw of my web apps. It's come to my attention that some people don't know what a proper physical address looks like, and so I'd like to make life a little easier. I have an app that requests a person's physical address. Instead of typing in the full thing, I'd like to use an autocomplete service. I'm in South Africa, so AfriGIS seems to be the best bet. However, I can't seem to find a service that will autocomplete an address for me, and split that same address (result) into it's parts. For example, in JSON:               {         'address': {             'streetNumber': 12,             'streetName': 'Main Street',             'suburbTown': 'Kensington',             'city': 'Johannesburg',             'province': 'Gauteng'             'streetCode': '2094',             <etc>         }     }      Are there any services that support this functionality (for South Africa, of course)? ### Update: This is very close to what I'm looking for: https://google- developers.appspot.com/maps/documentation/javascript/examples/full/places- autocomplete-addressform However, it is not very accurate in terms of complete addresses. In the above case, Google will not return the suburb. In South Africa, suburbs are areas, and not city outskirts. Also, the Zip codes it returns are wrong.


title: Proj4s.js Library Transformation Accuracy Documentation

text: I am using proj4s.js for transforming between coordinate systems for a GIS web application. Is there any documentation out there on transformations accuracy provided by the library?


title: Grass GIS and Python Error

text: I am trying to use Python with Grass GIS. To get started I try to read in a raster TIFF. The import of grass.script works fine:               import grass.script as g      But if I am trying to run a command I get the following error               sloperaster = 'slopet.tif'          a = g.run_command('r.in.gdal',                           input = sloperaster,                           output = 'test')          >>>'r.in.gdal' is not recognized as an internal or external command,     operable program or batch file.      Any ideas what is causing that and how to resolve it? I am on Windows 7 (64bit), Python 2.7.5 (32bit), Grass GIS 6.4.3 (32bit)


title: How to convert multiple csv files to shp using Java

text: I have some 100 CSV files. I want to convert them in shp format using Java. The csv file contains lat long information in WGS84. I can understand it is a batch processing and it is difficult to do in ArcGIS. is there any spatial plugin for eclipse that will do the task? I had a look at Geotools but found it bit complicated. Anything straight forward?


title: Unable to to run any geoprocessing tools on ArcGIS10. Log file shows that file is missing. Happened spontaneously

text: I have ArcGIS 10 installed on my computer connected remotely to my companies license server. Everything has worked fine for months, but within the last few days no geoprocessing tool will run. The results window just shows the dreaded 99998:Unexpected Error. Looking further into the log file I get this:               <Msg time='2012-10-31T22:33:50' type='ERROR' code='5000' target='GPJobServer.GPJobServer' thread='3872'>Job Catalog initialization failed. Directory C:\Program Files (x86)\ArcGIS\Desktop10.0\ArcToolbox\server\system\jobreg\ not found.</Msg>     <Msg time='2012-10-31T22:33:50' type='ERROR' code='1016' target='GPJobServer.GPJobServer' machine='KYLE-PC' thread='3872' elapsed='0.33800'>Server Object instance creation failed on machine KYLE-PC.   Job Catalog initialization failed. Directory C:\Program Files (x86)\ArcGIS\Desktop10.0\ArcToolbox\server\system\jobreg\ not found.  </Msg>     <Msg time='2012-10-31T22:33:50' type='ERROR' code='1008' target='GPJobServer.GPJobServer' thread='8072'>Configuration GPJobServer.GPJobServer can not be started.</Msg>     <Msg time='2012-10-31T22:33:50' type='ERROR' code='1006' target='GPJobServer.GPJobServer' thread='8072'>Configuration GPJobServer.GPJobServer can not be stopped.</Msg>     <Msg time='2012-10-31T22:33:50' type='WARNING' code='1034' target='GPJobServer.GPJobServer' thread='8072'>Warning: Changes to server object configuration GPJobServer.GPJobServer could not be saved to disk and will not be available when SOM restarts.  Reason for failure: Operation on file C:\Program Files (x86)\ArcGIS\Desktop10.0\ArcToolbox\server\user\cfg\GPJobServer.GPJobServer.cfg failed. File not found.</Msg>      Does anyone perchance know what might be going on? Where did the jobreg file go? Do I need to reinstall the software? This all happened completely randomly, and I never touched any of the AppData files, nor the Program Files. Thanks.


title: Problem with the where clause in ArcGIS when dealing with variables

text: I am writing a python script for ArcGIS. I have an input feature class which contains multiple pipelines which is saved as a variable `InputPipeline`. I would like to iterate through the pipelines and run some run some processes on them individually. The script is set up to prompt the user to choose the field which contains a unique identifier for each pipeline. I am saving this field as a variable called `PipelineIDField`. I run a search cursor to get the actual pipe name in that field and I'm saving that as `fullRouteName`. I'm using: `arcpy.MakeFeatureLayer_management(InputPipeline, 'SelectedPipe', where_clause)` and the where clause is defined as: `where_clause=""""%s" =""" %PipelineIDField + """ '%s' """ %fullRouteName` I would assume that I could then successfully use 'SelectedPipe' to run my processes on the individual pipe. This is successfully working when I run the where clause in the python window within ArcGIS, however, it doesn't select out the individual pipes in my script. The where clauses are identical. I have tried using escape characters "\"" and "\'" with the same result. Any idea what is going on here? I'm using ArcGIS 10.1 service pack 1 on Windows 7 64-bit. Thanks.


title: How to calculate the intersect two shapefiles

text: I have two shapefiles.   1. the clinic buffers 500m area.   2. the postcode area. Is there a simple way to calculate the intersect two shapefiles ?   I'm using QGIS Dufour. I tried to use fTools to intersect. However, it didn't show me actually how much every piece of the green area. It just combined the fields of two shapefiles.   Please let me know how to calculate the intersect two shapefiles.   Thank you in advance.   ![intersect two shapefiles](http://i.stack.imgur.com/UofDu.png) ![It just combined the fields of two shapefiles](http://i.stack.imgur.com/P6iaH.png) I used "Field Calculator" to calculate the intersection with $area. But those numbers are bigger than original area.![enter image description here](http://i.stack.imgur.com/esIgj.png)


title: Are there any open source software for obtaining and processing data from the Velodyne HDL-32E lidar?

text: Does anybody know if a FOSS package exists? The best I could find are the scripts at binarymillenium. I was just hoping somebody here knows of a complete package for obtaining and processing the data.


title: Port "Producing Building Shadows" Avenue code to ArcGIS 10

text: Whuber provided an answer here Producing building shadows? which required using Avenue code. Any idea how to make it work in Arc 10? Thanks in advance.


title: Rotated Labels not Respecting Feature Weights: Maplex Issue

text: When setting labels to rotate based on a field, annotation will overlap other features even if the weight is set to 1000. Does anyone know if there is a way to force it to follow both the rotation value and the feature weight? Even if the label engine had to increase the offset or shift the label to the other side of the feature or resize the labels within a set range to sqeeze it in? Open to any ideas or tools/scripts that can get this done. Using ArcMap 9.3 & 10 at the ArcINFO level with Maplex currently. Any ideas are appreciated. Thanks.


title: gdal_polygonize error with Arc/Info binary grid raster file

text: I am trying to convert an Arc/Info binary grid raster file to a shapefile using the gdal_polygonize.py utility. But the operation fails with the following error > Warning 1: Unsupported Arc/Info Binary Grid tile of type 0x80 encountered. > This and subsequent unsupported tile types set to no data value. Here is a screenshot of the error ![enter image description here](http://i.stack.imgur.com/UBlEs.png) I am able to open and view the binary grid file just fine in ArcMAP etc. Here are the results when i run gdalinfo on the binary grid file. It seems to be fine               Driver: AIG/Arc/Info Binary Grid     Files: /home/vish/Sample_Data/nh_priority            /home/vish/Sample_Data/nh_priority.ovr            /home/vish/Sample_Data/nh_priority/vat.adf            /home/vish/Sample_Data/nh_priority/w001001x.adf            /home/vish/Sample_Data/nh_priority/w001001.adf            /home/vish/Sample_Data/nh_priority/hdr.adf            /home/vish/Sample_Data/nh_priority/dblbnd.adf            /home/vish/Sample_Data/nh_priority/prj.adf            /home/vish/Sample_Data/nh_priority/sta.adf     Size is 4987, 9645     Coordinate System is:     PROJCS["NAD83 / New Hampshire",         GEOGCS["NAD83",             DATUM["North_American_Datum_1983",                 SPHEROID["GRS 1980",6378137,298.257222101,                     AUTHORITY["EPSG","7019"]],                 AUTHORITY["EPSG","6269"]],             PRIMEM["Greenwich",0,                 AUTHORITY["EPSG","8901"]],             UNIT["degree",0.01745329251994328,                 AUTHORITY["EPSG","9122"]],             AUTHORITY["EPSG","4269"]],         UNIT["METERS",1],         PROJECTION["Transverse_Mercator"],         PARAMETER["latitude_of_origin",42.5],         PARAMETER["central_meridian",-71.66666666666667],         PARAMETER["scale_factor",0.999966667],         PARAMETER["false_easting",300000],         PARAMETER["false_northing",0],         AUTHORITY["EPSG","32110"],         AXIS["X",EAST],         AXIS["Y",NORTH]]     Origin = (228040.838957945175935,311318.187724410032388)     Pixel Size = (30.000000000000000,-30.000000000000000)     Corner Coordinates:     Upper Left  (  228040.839,  311318.188) ( 72d35'2.81"W, 45d17'53.85"N)     Lower Left  (  228040.839,   21968.188) ( 72d32'41.44"W, 42d41'39.84"N)     Upper Right (  377650.839,  311318.188) ( 70d40'35.98"W, 45d17'51.67"N)     Lower Right (  377650.839,   21968.188) ( 70d43'8.53"W, 42d41'37.85"N)     Center      (  302845.839,  166643.188) ( 71d37'52.26"W, 44d 0'0.06"N)     Band 1 Block=256x16 Type=Byte, ColorInterp=Palette       Min=0.000 Max=2.000        NoData Value=255       Overviews: 2494x4823, 1247x2412, 624x1206, 312x603, 156x302       Color Table (RGB with 3 entries)         0: 255,255,255,255         1: 211,255,190,255         2: 0,84,3,255     <GDALRasterAttributeTable>       <FieldDefn index="0">         <Name>VALUE</Name>         <Type>0</Type>         <Usage>5</Usage>       </FieldDefn>       <FieldDefn index="1">         <Name>COUNT</Name>         <Type>0</Type>         <Usage>1</Usage>       </FieldDefn>       <FieldDefn index="2">         <Name>STEW_STATUS</Name>         <Type>2</Type>         <Usage>0</Usage>       </FieldDefn>       <Row index="0">         <F>0</F>         <F>7209043</F>         <F>non-eligible</F>       </Row>       <Row index="1">         <F>1</F>         <F>7672024</F>         <F>eligible, non-priority</F>       </Row>       <Row index="2">         <F>2</F>         <F>11779591</F>         <F>priority</F>       </Row>     </GDALRasterAttributeTable>


title: How to set spatial reference to OSGB_1936_British_National_Grid in ArcGIS

text: I have two datasets in ArcGIS. One is from the Ordnance Survey (UK National Mapping Agency) and, when I view the properties of the dataset it says that the Spatial Reference is               OSGB_1936_British_National_Grid      I want to set another file to the same spatial reference. I know I can use the _Import_ function to import it from the first file, but how would I select this co-ordinate system from the built-in ArcGIS ones? I can find `British National Grid` under Projected/National Grids/Europe and I can find `OSGB 1936` under Geographic/Europe, but I can't seem to find the combination of the two anywhere? I'm not hugely experienced with co-ordinate systems, but I think that BNG is the co-ordinate system, and OSGB 1936 is the datum. Is that correct? I know that when I select spatial references in ENVI I can choose a co-ordinate system of `United Kingdom` and a Datum of `OSGB1936` - is there a way to do this in ArcGIS?


title: How to create a single depth contour line from a bathymetric raster?

text: I have bathymetric data (raster form) that I'm trying to tease out specific depth contours. For example, I would like to "create" the 150-m depth contour from the bathymetric data in the Gulf of Alaska. By using the Raster Calculator, I was able to create a raster of just the 150 m depth cells. However, I would like these to be connected as a "smooth" contour line. I'm having a hard time thinking of how to go about this. I hope this makes sense!


title: How to determine what Level a new Extent would be, without zooming to it

text: In the ArcGIS Server JavaScript API 2.4, there are some related methods for getting and setting the map's extent.   * `map.getLevel()` determine's the map's current level (when using a Tiled layer)   * `map.setExtent(extent)` changes the map's extent to the new extent   * `map.setLevel(level)` changes the map's level to the new level Is it possible to determine what **level** the value of `extent` would equate to, **before** changing to that extent? That is, **before** I pass in `map.setExtent(x1, y1, x2, y2)` I want to know what level this would display, so that I can adjust the value as necessary. I can set a listener on `map.onExtentChange`, verify the new level, and use `map.setLevel(map.getLevel() + 1)` to fix the extent if necessary - but this means the map may draw twice (once after `map.setExtent` then again after `map.setLevel`) As far as I can see, there is no equivalent of `extent.getLevel` - is there a way to calculate this without actually changing the map level first?


title: UnIntersected part of gdal-geometry?

text: i have a rectangle and polygon geometries . polygon.Instersection(rectangle) gives me the polygon as it is inside the rectangle. But , what i want is the unintersected(non-intersecting) part of both geometries, basically (rectangle - polygon ).Is this possible in Gdal ?.


title: QGIS 2.0 Raster Interpolation

text: I am trying to use the Rasterinterpolation plugin to pick up elevations from a Raster DEM for a series of vector points. Once loaded, the Plugin is accessible through Plugins Menu/ Raster Interpolation and seems set up to do exactly what I need. However it will not let me select the correct Band of the Raster Layer. The plugin sees the layer as only having Band 0, whereas the Raster only has Band 1, and as a result does not populate the 'height' field in the Vector layer. SciPy and NumPy are both loaded I have tried the same plugin with a multiband Raster layer and it still only indicates that the Raster has a single band, Band 0. Any help greatly appreciated.


title: Find polygon that contains all linestring records in PostGIS table

text: I'm still pretty new to this stuff, so bear with me. In PostGIS how could I find a polygon that would contain all the linestring records in a given table?


title: How to point out a city from a map image of world with longitude and latitude?

text: have a map image 800px wide. How can I figure out a city. like New York has 40.75 latitude and -73.98 longitude. I have to pin a mark over this point. Before I have done with x_axis and y_axis which gives problem on more cities figure because i manually had to input values. now i have database with all cities with latitude and longitude. or is there any way that i can convert latitude and longitude to x and y axis. Thank you **EDIT**               x = image.width * (longitude + 180) / (2 * 180)     y = image.height * (latitude + 180) / (2 * 180)


title: How could I label features and convert them to annotations with arcpy?

text: how could I do the following steps with arcpy:   1. right-click on the layer with the features   2. Label Features   3. Convert Labels to Annotations ... The annotions should be saved in a new shapefile. How could I do this in python?


title: Display a local file in Google Maps

text: How can I get a KML/KMZ file to display on Google Maps without a public facing web server?


title: Import GPX file to QGIS and correct apparent CRS problems

text: I am using Easy Trails software on an iPhone to output a gpx file of a walk. When I import the file (import vector layer) to QGIS, into a project with OSGB information and the correct CRS for OSGB data, the imported track appears well to the West and South of where it should be. I have tried changing the CRS in both layers, tried changing the project CRS, enabled and disabled on the fly projection but no combination of what I try will get the track into the right place. I feel I must be doing something fundamentally wrong but cannot work out what it might be. Am I missing some vital step out?


title: Building a thin mapping client with OpenLayers

text: I have installed GeoServer, inside of tomcat, imported my postgis tables into geoserver and have created my styles etc. I can preview them using openlayers in geoserver. Now, i want to build a simple mapping client using OpenLayers. So i download the openlayer.x zip file, i unzipped it to my geoserver/data/www folder. My geoserver sits inside the web apps folder of tomcat. I also copied the wms.html file from the openlayers example folder into this www folder as well and changed parameters like the src to point to my geoserver, the layer and the openlayers.js path. This is following the instruction on this site: https://www.e-education.psu.edu/geog585/l8_p5.html However when i try to open mylayers.html, i dont see the layers which i specified to show. I only see the text and no map shown. Below is the code for my client (mylayers.html) which is saved in the www folder: i have removed the style parameter because of space limitation here, but the url for style.css is localhost:8080/geoserver/OpenLayers-2.12/theme/default/style.css               <?xml version="1.0" encoding="UTF-8"?>     <html xmlns="http://www.w3.org/1999/xhtml">              <script src="http://localhost:8080/geoserver/OpenLayers-2.12/OpenLayers.js"   type="text/javascript">         </script>         <script defer="defer" type="text/javascript">             var map;             var Layer1, Layer2, Layer3, Layer4, Layer5, Layer6, Layer7, Layer8, Layer9, Layer10, Layer11, Layer12;             var pureCoverage = false;                  // pink tile avoidance             OpenLayers.IMAGE_RELOAD_ATTEMPTS = 5;             // make OL compute scale according to WMS spec             OpenLayers.DOTS_PER_INCH = 25.4 / 0.28;                  function init(){                 // if this is just a coverage or a group of them, disable a few items,                 // and default to jpeg format                 format = 'image/png';                 if(pureCoverage) {                     document.getElementById('filterType').disabled = true;                     document.getElementById('filter').disabled = true;                     document.getElementById('antialiasSelector').disabled = true;                     document.getElementById('updateFilterButton').disabled = true;                     document.getElementById('resetFilterButton').disabled = true;                     document.getElementById('jpeg').selected = true;                     format = "image/jpeg";                 }                      var bounds = new OpenLayers.Bounds(                     143.971463002, -8.04166701399993,                     151.699865035, -3.99888201799996                 );                 var options = {                     controls: [],                     maxExtent: bounds,                     maxResolution: 0.0301890704414062,                     projection: "EPSG:4326",                     units: 'degrees'                 };                 map = new OpenLayers.Map('map', options);                           Layer1 = new OpenLayers.Layer.WMS("Trial1:province - Untiled",    "http://localhost:8080/geoserver/wms",                     {                         layers: 'Trial1:province',                         styles: '',                         format: 'format'                         width:'800'                         height:'300'                         srs: 'EPSG:4326'                     },                     {                        singleTile: true,                         ratio: 1,                         isBaseLayer: true,                     }                  );                      layer2 = new OpenLayers.Layer.WMS(                     "Trial1:vct_sites - Untiled", "http://localhost:8080/geoserver/wms",                     {                         layers: 'Trial1:vct_sites',                         styles: '',                         format: 'format''                         width:'800'                         height:'300'                         srs: 'EPSG:4326'                         transparent: true                     },                     {                        singleTile: true,                         ratio: 1,                         isBaseLayer: false,                          }                  );                      Layer3 = new OpenLayers.Layer.WMS(                     "Trial1:road_network - Untiled", "http://localhost:8080/geoserver/wms",                     {                         layers: 'Trial1:road_network',                         styles: '',                         format: 'format''                         width:'800'                         height:'300'                         srs: 'EPSG:4326'                         transparent: true                     },                     {                        singleTile: true,                         ratio: 1,                         isBaseLayer: false,                          }                  );                       Layer4 = new OpenLayers.Layer.WMS(                     "Trial1:major_towns - Untiled", "http://localhost:8080/geoserver/wms",                     {                         layers: 'Trial1:major_towns',                         styles: '',                         format: 'format''                         width:'800'                         height:'300'                         srs: 'EPSG:4326'                         transparent: true                     },                     {                        singleTile: true,                         ratio: 1,                         isBaseLayer: false,                          }                  );                       Layer5 = new OpenLayers.Layer.WMS(                     "Trial1:llg - Untiled", "http://localhost:8080/geoserver/wms",                     {                         layers: 'Trial1:llg',                         styles: '',                         format: 'format''                         width:'800'                         height:'300'                         srs: 'EPSG:4326'                         transparent: true                     },                     {                        singleTile: true,                         ratio: 1,                         isBaseLayer: false,                          }                  );                      Layer6 = new OpenLayers.Layer.WMS(                     "Trial1:health_centers - Untiled", "http://localhost:8080/geoserver/wms",                     {                         layers: 'Trial1:health_centers',                         styles: '',                         format: 'format''                         width:'800'                         height:'300'                         srs: 'EPSG:4326'                         transparent: true                     },                     {                        singleTile: true,                         ratio: 1,                         isBaseLayer: true,                          }                  );                      Layer7 = new OpenLayers.Layer.WMS(                     "Trial1:hc_ap_links - Untiled", "http://localhost:8080/geoserver/wms",                     {                         layers: 'Trial1:hc_ap_links',                         styles: '',                         format: 'format''                         width:'800'                         height:'300'                         srs: 'EPSG:4326'                         transparent: true                     },                     {                        singleTile: true,                         ratio: 1,                         isBaseLayer: false,                          }                  );                      Layer8 = new OpenLayers.Layer.WMS(                     "Trial1:district - Untiled", "http://localhost:8080/geoserver/wms",                     {                          layers: 'Trial1:district',                         styles: '',                         format: 'format''                         width:'800'                         height:'300'                         srs: 'EPSG:4326'                         transparent: true                     },                     {                        singleTile: true,                         ratio: 1,                         isBaseLayer: false,                          }                  );                      Layer9 = new OpenLayers.Layer.WMS(                     "Trial1:census_units - Untiled", "http://localhost:8080/geoserver/wms",                     {                         layers: 'Trial1:census_units',                         styles: '',                         format: 'format''                         width:'800'                         height:'300'                         srs: 'EPSG:4326'                         transparent: true                     },                     {                        singleTile: true,                         ratio: 1,                         isBaseLayer: false,                          }                  );                      Layer10 = new OpenLayers.Layer.WMS(                     "Trial1:areastore_hc_links - Untiled", "http://localhost:8080/geoserver/wms",                     {                         layers: 'Trial1:areastore_hc_links',                         styles: '',                         format: 'format''                         width:'800'                         height:'300'                         srs: 'EPSG:4326'                         transparent: true                     },                     {                        singleTile: true,                         ratio: 1,                         isBaseLayer: false,                          }                  );                      Layer11 = new OpenLayers.Layer.WMS(                     "Trial1:airstrips - Untiled", "http://localhost:8080/geoserver/wms",                     {                         layers: 'Trial1:airstrips',                         styles: '',                         format: 'format''                         width:'800'                         height:'300'                         srs: 'EPSG:4326'                         transparent: true                     },                     {                        singleTile: true,                         ratio: 1,                         isBaseLayer: false,                          }                  );                       Layer12 = new OpenLayers.Layer.WMS(                     "Trial1:aidpost - Untiled", "http://localhost:8080/geoserver/wms",                     {                         layers: 'Trial1:aidpost',                         styles: '',                         format: 'format''                         width:'800'                         height:'300'                         srs: 'EPSG:4326'                         transparent: true                     },                     {                        singleTile: true,                         ratio: 1,                         isBaseLayer: false,                          }                  );                           map.addLayers([Layer1, Layer2, Layer3, Layer4, Layer5, Layer6, Layer7, Layer8, Layer9, Layer10, Layer11, Layer12]);                      // build up all controls                 map.addControl(new OpenLayers.Control.PanZoomBar({                     position: new OpenLayers.Pixel(2, 15)                 }));                 map.addControl(new OpenLayers.Control.Navigation());                 map.addControl(new OpenLayers.Control.Scale($('scale')));                 map.addControl(new OpenLayers.Control.MousePosition({element: $('location')}));                 map.zoomToExtent(bounds);                      // wire up the option button                 var options = document.getElementById("options");                 options.onclick = toggleControlPanel;                      // support GetFeatureInfo                 map.events.register('click', map, function (e) {                     document.getElementById('nodelist').innerHTML = "Loading... please wait...";                     var params = {                         REQUEST: "GetFeatureInfo",                         EXCEPTIONS: "application/vnd.ogc.se_xml",                         BBOX: map.getExtent().toBBOX(),                         SERVICE: "WMS",                         INFO_FORMAT: 'text/html',                         QUERY_LAYERS: map.layers[0].params.LAYERS,                         FEATURE_COUNT: 50,                         Layers: 'Trial1:province',                         WIDTH: map.size.w,                         HEIGHT: map.size.h,                         format: format,                         styles: map.layers[0].params.STYLES,                         srs: map.layers[0].params.SRS};                          // handle the wms 1.3 vs wms 1.1 madness                     if(map.layers[0].params.VERSION == "1.3.0") {                         params.version = "1.3.0";                         params.j = parseInt(e.xy.x);                         params.i = parseInt(e.xy.y);                     } else {                         params.version = "1.1.1";                         params.x = parseInt(e.xy.x);                         params.y = parseInt(e.xy.y);                     }                          // merge filters                     if(map.layers[0].params.CQL_FILTER != null) {                         params.cql_filter = map.layers[0].params.CQL_FILTER;                     }                      if(map.layers[0].params.FILTER != null) {                         params.filter = map.layers[0].params.FILTER;                     }                     if(map.layers[0].params.FEATUREID) {                         params.featureid = map.layers[0].params.FEATUREID;                     }                     OpenLayers.loadURL("http://localhost:8080/geoserver/Trial1/wms", params, this, setHTML, setHTML);                     OpenLayers.Event.stop(e);                 });             }                  // sets the HTML provided into the nodelist element             function setHTML(response){                 document.getElementById('nodelist').innerHTML = response.responseText;             };                  // shows/hide the control panel             function toggleControlPanel(event){                 var toolbar = document.getElementById("toolbar");                 if (toolbar.style.display == "none") {                     toolbar.style.display = "block";                 }                 else {                     toolbar.style.display = "none";                 }                 event.stopPropagation();                 map.updateSize()             }                  // Tiling mode, can be 'tiled' or 'untiled'             function setTileMode(tilingMode){                 if (tilingMode == 'tiled') {                     untiled.setVisibility(false);                     tiled.setVisibility(true);                     map.setBaseLayer(tiled);                 }                 else {                     untiled.setVisibility(true);                     tiled.setVisibility(false);                     map.setBaseLayer(untiled);                 }             }                  // Transition effect, can be null or 'resize'             function setTransitionMode(transitionEffect){                 if (transitionEffect === 'resize') {                     tiled.transitionEffect = transitionEffect;                     untiled.transitionEffect = transitionEffect;                 }                 else {                     tiled.transitionEffect = null;                     untiled.transitionEffect = null;                 }             }                  // changes the current tile format             function setImageFormat(mime){                 // we may be switching format on setup                 if(tiled == null)                   return;                      tiled.mergeNewParams({                     format: mime                 });                 untiled.mergeNewParams({                     format: mime                 });                 /*                 var paletteSelector = document.getElementById('paletteSelector')                 if (mime == 'image/jpeg') {                     paletteSelector.selectedIndex = 0;                     setPalette('');                     paletteSelector.disabled = true;                 }                 else {                     paletteSelector.disabled = false;                 }                 */             }                  // sets the chosen style             function setStyle(style){                 // we may be switching style on setup                 if(tiled == null)                   return;                      tiled.mergeNewParams({                     styles: style                 });                 untiled.mergeNewParams({                     styles: style                 });             }                  // sets the chosen WMS version             function setWMSVersion(wmsVersion){                 // we may be switching style on setup                 if(wmsVersion == null)                   return;                      if(wmsVersion == "1.3.0") {                    origin = map.maxExtent.bottom + ',' + map.maxExtent.left;                 } else {                    origin = map.maxExtent.left + ',' + map.maxExtent.bottom;                 }                      tiled.mergeNewParams({                     version: wmsVersion,                     tilesOrigin : origin                 });                 untiled.mergeNewParams({                     version: wmsVersion                 });             }                  function setAntialiasMode(mode){                 tiled.mergeNewParams({                     format_options: 'antialias:' + mode                 });                 untiled.mergeNewParams({                     format_options: 'antialias:' + mode                 });             }                  function setPalette(mode){                 if (mode == '') {                     tiled.mergeNewParams({                         palette: null                     });                     untiled.mergeNewParams({                         palette: null                     });                 }                 else {                     tiled.mergeNewParams({                         palette: mode                     });                     untiled.mergeNewParams({                         palette: mode                     });                 }             }                  function setWidth(size){                 var mapDiv = document.getElementById('map');                 var wrapper = document.getElementById('wrapper');                      if (size == "auto") {                     // reset back to the default value                     mapDiv.style.width = null;                     wrapper.style.width = null;                 }                 else {                     mapDiv.style.width = size + "px";                     wrapper.style.width = size + "px";                 }                 // notify OL that we changed the size of the map div                 map.updateSize();             }                  function setHeight(size){                 var mapDiv = document.getElementById('map');                      if (size == "auto") {                     // reset back to the default value                     mapDiv.style.height = null;                 }                 else {                     mapDiv.style.height = size + "px";                 }                 // notify OL that we changed the size of the map div                 map.updateSize();             }                  function updateFilter(){                 if(pureCoverage)                   return;                      var filterType = document.getElementById('filterType').value;                 var filter = document.getElementById('filter').value;                      // by default, reset all filters                 var filterParams = {                     filter: null,                     cql_filter: null,                     featureId: null                 };                 if (OpenLayers.String.trim(filter) != "") {                     if (filterType == "cql")                          filterParams["cql_filter"] = filter;                     if (filterType == "ogc")                          filterParams["filter"] = filter;                     if (filterType == "fid")                          filterParams["featureId"] = filter;                 }                 // merge the new filter definitions                 mergeNewParams(filterParams);             }                  function resetFilter() {                 if(pureCoverage)                   return;                      document.getElementById('filter').value = "";                 updateFilter();             }                  function mergeNewParams(params){                 tiled.mergeNewParams(params);                 untiled.mergeNewParams(params);             }         </script>     </head>     <body onload="init()">         <div id="toolbar" style="display: none;">             <ul>                 <li>                     <a>WMS version:</a>                     <select id="wmsVersionSelector" onchange="setWMSVersion(value)">                         <option value="1.1.1">1.1.1</option>                         <option value="1.3.0">1.3.0</option>                     </select>                 </li>                 <li>                     <a>Tiling:</a>                     <select id="tilingModeSelector" onchange="setTileMode(value)">                         <option value="untiled">Single tile</option>                         <option value="tiled">Tiled</option>                     </select>                 </li>                 <li>                     <a>Transition effect:</a>                     <select id="transitionEffectSelector" onchange="setTransitionMode(value)">                         <option value="">None</option>                         <option value="resize">Resize</option>                     </select>                 </li>                 <li>                     <a>Antialias:</a>                     <select id="antialiasSelector" onchange="setAntialiasMode(value)">                         <option value="full">Full</option>                         <option value="text">Text only</option>                         <option value="none">Disabled</option>                     </select>                 </li>                 <li>                     <a>Format:</a>                     <select id="imageFormatSelector" onchange="setImageFormat(value)">                         <option value="image/png">PNG 24bit</option>                         <option value="image/png8">PNG 8bit</option>                         <option value="image/gif">GIF</option>                         <option id="jpeg" value="image/jpeg">JPEG</option>                     </select>                 </li>                 <li>                     <a>Styles:</a>                     <select id="imageFormatSelector" onchange="setStyle(value)">                         <option value="">Default</option>                     </select>                 </li>                 <!-- Commented out for the moment, some code needs to be extended in                       order to list the available palettes                 <li>                     <a>Palette:</a>                     <select id="paletteSelector" onchange="setPalette(value)">                         <option value="">None</option>                         <option value="safe">Web safe</option>                     </select>                 </li>                 -->                 <li>                     <a>Width/Height:</a>                     <select id="widthSelector" onchange="setWidth(value)">                         <!--                         These values come from a statistics of the viewable area given a certain screen area                         (but have been adapted a litte, simplified numbers, added some resolutions for wide screen)                         You can find them here: http://www.evolt.org/article/Real_World_Browser_Size_Stats_Part_II/20/2297/                         --><option value="auto">Auto</option>                         <option value="600">600</option>                         <option value="750">750</option>                         <option value="950">950</option>                         <option value="1000">1000</option>                         <option value="1200">1200</option>                         <option value="1400">1400</option>                         <option value="1600">1600</option>                         <option value="1900">1900</option>                     </select>                     <select id="heigthSelector" onchange="setHeight(value)">                         <option value="auto">Auto</option>                         <option value="300">300</option>                         <option value="400">400</option>                         <option value="500">500</option>                         <option value="600">600</option>                         <option value="700">700</option>                         <option value="800">800</option>                         <option value="900">900</option>                         <option value="1000">1000</option>                     </select>                 </li>                 <li>                     <a>Filter:</a>                     <select id="filterType">                         <option value="cql">CQL</option>                         <option value="ogc">OGC</option>                         <option value="fid">FeatureID</option>                     </select>                     <input type="text" size="80" id="filter"/>                     <img id="updateFilterButton" src="http://localhost:8080/geoserver/OpenLayers-2.12/img/east-mini.png" onClick="updateFilter()" title="Apply filter"/>                     <img id="resetFilterButton" src="http://localhost:8080/geoserver/OpenLayers-2.12/img/cancel.png" onClick="resetFilter()" title="Reset filter"/>                 </li>             </ul>         </div>         <div id="map">             <img id="options" title="Toggle options toolbar" src="http://localhost:8080/geoserver/options.png"/>         </div>         <div id="wrapper">             <div id="location">location</div>             <div id="scale">             </div>         </div>         <div id="nodelist">             <em>Click on the map to get feature info</em>         </div>     </body>      Using the chrome developer tool - console, i see the following errors. Failed to load resource: the server responded with a status of 404 (Not Found) Failed to load resource: the server responded with a status of 404 (Not Found) Uncaught SyntaxError: Unexpected identifier mylayers.html:150 Failed to load resource: the server responded with a status of 404 (Not Found) Failed to load resource: the server responded with a status of 404 (Not Found) Uncaught ReferenceError: init is not defined mylayers.html:632 What could be the problem here? as im sure my folders and files are where they should be. Any help would be very much appreciated. Thanks Barbara


title: How to get raster image as array in Python with ArcGIS

text: When starting to work with Python and ArcGIS 9.3, I assumed there would be a simple way to get a raster image into a Python array so that I can manipulate it before storing it back as another raster image. However, I can't seem to find out how to do this. Is it possible? If so, how?


title: How to use map composer in a stand-alone script?

text: I'm trying to follow the map rendering section from the pyqgis cookbook, but I'd like to test this as a standalone application. I can do the first part, using simple rendering, but I'm a bit stuck doing the second example using the map composer as a standalone script. Here is a standalone example for the bit I can do:               from qgis.core import *     from qgis.gui import *     from PyQt4.QtCore import *     from PyQt4.QtGui import *     from PyQt4.QtXml import *          QgsApplication.setPrefixPath("/usr/", True)     QgsApplication.initQgis()          fh = open("eg.csv","w")     fh.write("""     x,y,name     153.0278, -27.4679, Brisbane     144.2500, -23.4500, Longreach     145.7753, -16.9256, Cairns     """)     fh.close()          uri = "eg.csv?delimiter=%s&xField=%s&yField=%s" % (",", "x", "y")     layer = QgsVectorLayer(uri, "eglayer", "delimitedtext")     QgsMapLayerRegistry.instance().addMapLayer(layer)     img = QImage(QSize(800,600), QImage.Format_ARGB32_Premultiplied)     color = QColor(255,255,255)     img.fill(color.rgb())     p = QPainter()     p.begin(img)     render = QgsMapRenderer()     lst = [ layer.getLayerID() ]  # add ID of every layer     render.setLayerSet(lst)     rect = QgsRectangle(render.fullExtent())     rect.scale(1.1)     render.setExtent(rect)     render.setOutputSize(img.size(), img.logicalDpiX())     render.render(p)     p.end()     img.save("render.png","png")      What I'd really like to do is the same, but use `QgsComposition`, and save as for example pdf. The cookbook says: > When using composer in a standalone application, you can create your own map > renderer instance the same way as shown in the section above and pass it to > the composition. This bit I couldn't do, all my attempts either get an empty map, or a segfault. I'm running linux mint 13, using qgis 1.8.0. It would be great if someone could show me how to modify the simple example to one that uses the composer.


title: How to fill in the parameters for kernel density estimation?

text: I am desperately trying to calculate animals home range with kernel. Unfortunately, the result is obviously incorrect and I cannot identify my mistake (I would like to post a screenshot but its not allowed to me yet). is there an instruction how to fill in the parameters for kernel density estimation? Is there a concept how the attribute table has to look like? (i have 3 columns: ID, longitude, latitude in decimal degrees) Many thanks in advance for help! QGIS 1.8


title: GCS_WGS_1984 Map in Esri Javascript

text: I created a dynamic map service in Arcmap and published to my server. The spatial reference used: GCS_WGS_1984. I used the javascript without basemap to display in a html. I cannot do a centerandzoom(). I suspect it's the spatial reference because if I add a basemap from esri the centerandzoom() works. How Can I solve this?


title: Convert Arcpy results to feature class in ArcMap 10.1

text: I am using ArcGIS 10.1 and am a newbie in Python. Currently we are executing stored procedure (from sql server) through ArcPy in ArcMap. I want to create feature class from executing following python script. I researched on Google but unable to find proper way/parameter in ArcPy. Currently going through following parameters GetParameter & GetMessage Any help/steps will be great! My ArcPy Script :               import sys     import arcpy     from arcpy import env          try:         # Make data path relative         arcpy.env.workspace = sys.path[0]          # Two ways to create the object, which also creates the     # connection to ArcSDE.     # Using the first method, pass a set of strings containing     # the connection properties:     #   <serverName>, <portNumber>, <version>, <userName>, <password>     #   arcpy.ArcSDESQLExecute("locolhost","5151","#","Admin","password")     # Using the second method pass the path to a valid ArcSDE connection file     sde_conn = arcpy.ArcSDESQLExecute(r"C:\Users\Administrator\AppData\Roaming\ESRI\Desktop10.1\ArcCatalog\test.sde")          # Get the SQL statements, separated by ; from a text string.     #sql_statement = 'select Name from all_points'     sql_statement = 'Exec [dbo].[sp_test]'     sql_statement_list = sql_statement.split(";")          print("+++++++++++++++++++++++++++++++++++++++++++++\n")          # For each SQL statement passed in, execute it.     for sql in sql_statement_list:         print("Execute SQL Statement: {0}".format(sql))         try:             # Pass the SQL statement to the database.             sde_return = sde_conn.execute(sql)         except Exception as err:             print(err)             sde_return = False              # If the return value is a list (a list of lists), display         # each list as a row from the table being queried.         if isinstance(sde_return, list):             print("Number of rows returned by query: {0} rows".format(                 len(sde_return)))             for row in sde_return:                 print(row)             print("+++++++++++++++++++++++++++++++++++++++++++++\n")         else:             # If the return value was not a list, the statement was             # most likely a DDL statment. Check its status.             if sde_return == True:                 print("SQL statement: {0} ran successfully.".format(sql))             else:                 print("SQL statement: {0} FAILED.".format(sql))             print("+++++++++++++++++++++++++++++++++++++++++++++\n")


title: Merge separate features into single feature in a shapefile programmatically

text: I have a shapefile in which some features shares same values for all the attributes. I can merge them manually by going into edit session and, selecting features to be merged, and then using merge button under Editor. But, since number of such features are very large, a programmatically solution is required.Any help is appreciated. Thanks in advance. looking forward for some genuine help.


title: How to Pass Multiple Dissolve Fields in VBA Code of ibasicgeoprocessor.dissolve Method

text: When we use the ibasicgeoprocessor.dissolve method it supports only one dissolve field but if we pass more than one dissolve field then it gives an exception: `value does not fall with in the range exception` Please give me any VBA, VB.Net, C#.Net code example of the ibasicgeoprocessor.Dissolve method having multiple dissolve fields.


title: Trying to convert coordinates from WGS84 Web Mercator Auxiliary Sphere to WGS84

text: I've been struggling with this for way too long now. I've got some data that was georeferenced in the WGS84 Web Mercator Auxiliary Sphere that should have been done just in WGS84. The Auxiliary Sphere (ESPG: 3857) is the projection commonly used by Google, Bing, and recently ArcGIS Online. The projected coordinates I have are in meters, as per Aux Sphere specs, but I've got to convert them to degrees and bring them into WGS84. I know the Aux Sphere is just a projected version of WGS84, but I can't seem to wrap my mind around converting coordinates from a projected coordinate system to a geographic coordinate system. I have some code that will convert from WGS84 Web Mercator (WKID: 102113) to WGS84, but because that system uses a spherical globe, rather than the ellipsoid used by WGS84 and the Aux Sphere, my results are a bit off. I've tried modifying it, but I can't seem to get it just right. Here's my code if that is helpful:               public static void ToGeographic(ref double[] mercator)         {             if (Math.Abs(mercator[4]) < 180 && Math.Abs(mercator[5]) < 90)                 return;                  if ((Math.Abs(mercator[4]) > 20037508.3427892) || (Math.Abs(mercator[5]) > 20037508.3427892))                 return;                  double x = mercator[4];             double y = mercator[5];             double num3 = x / 6378137.0;             double num4 = num3 * 57.295779513082323;             double num5 = Math.Floor((double)((num4 + 180.0) / 360.0));             double num6 = num4 - (num5 * 360.0);             double num7 = 1.5707963267948966 - (2.0 * Math.Atan(Math.Exp((-1.0 * y) / 6378137.0)));             mercator[4] = num6;             mercator[5] = num7 * 57.295779513082323;                   }


title: SDE Replication works through ArcMap but not through ArcCatalog

text: I'm running 10.1 on a 10.0 SDE database. We're trying to replicate our database (one-way parent to child) to a file geodatabase. Ideally I'd like to have this automated using python or .bat files but before I do that I want to verify that all the feature classes are being transferred over correctly. This is where I run into a problem, with one feature dataset which contains 5 feature classes.. if I run this replica through ArcCatalog only two of these feature classes are copied over... if I run it through ArcMap all of the feature classes are transferred over as well... They all have global IDs and are versioned... What am I missing? What else can I check for to account for this weird behavior! I'm at my wit's end. Thanks.


title: How to update attributes of a shapefile using list values?

text: I am using an update cursor to update the attributes of a shapefile using values from a list. It seems very easy, but I couldn't make it run correctly. When I run the following code, it updates all rows based on the final value in the list. Please suggest how to fix it.               cur1 = arcpy.UpdateCursor(fc)     for v in cur1:         for s in list:             val = s             v.Val = val             cur1.updateRow(v)


title: How to georeference a raster based on identical, already georeferenced raster

text: I have two identical rasters, one has been georeferenced, the other has not. I would like to georeference the second raster automatically using the first one. I do not want to create links by hand and the raster has been warped so simply shifting xy position does not work. Let me know if you have any ideas... most comfortable working in model builder but would love any solution. Thanks!


title: Cannot add ArcGIS service from a Linux based server in ArcMap (admin)

text: I have installed ArcGIS Server 10.2 on Ubuntu 13.04 (with few tricks which I would be happy to share) and got the server up and running to the extent that I can manage the server through the webinterface and even add it in ArcMap for accessing the maps. My problem is that when I add the server for management, it asks for a hostname and when I enter the same name as it appears on the ArcGIS web interface (where the machines managed by the server is listed), ArcMap complains about not being able to find the machine with the given hostname. Any ideas?


title: HTML Popup Database Actions

text: I'm working on a collection of points on the map, each comes with a bunch of properties. I succeeded using a customized XSL template to show the information when I click on any point. And now what I wanna achieve is to add a "next" & "previous" button in the window, which navigates to the information of next adjacent point. I don't know if this idea is possible at all, or do I have to write a add-in for this function? Any help would be appreciated.