title: How much CAD experience for the GIS Professional?

text: As a GIS Professional, one is often expected to have experience with CAD (i.e. AutoCAD, Microstation) to the extent that they can integrate CAD data into GIS workflows and datasets. However, I often see job postings in which employers are looking for candidates that have experience actually using CAD software such as AutoCAD. I would expect that drafting would be left to people with architectural and engineering backgrounds. So my question is, in people's experience, what are employers actually expecting from GIS professionals with regard to CAD experience using CAD-specific software? Secondly, I would appreciate any insight into how I could go about trying to develop those skills in my spare time. Thanks


title: How can I draw custom shapes on Graphics layer?

text: Hello everyone I need help. Any clue, guidance or help that would be appreciated, please bear with me as I am still struggling with Flex 4.   1. What is the simplest way to draw custom shapes using < s : path> on graphic layer. I had data for the already created.   2. How can I pass my own styles in SimpleMarkerSymbol(SimpleMarkerSymbol.STYLE_CIRCLE, 22, 0x009933) in place of STYLE_CIRCLE? I want to create my own styles so that I can explicitly use them there like STYLE_RECT. Can I pass my path data in any symbol like fillsymbol, linesymbol or markersymbol so that they draw the shape?   3. How can I extend these symbol classes so that I can use member functions and properties of these symbols + my MXML code? I am using SFV 2.4. Thanks


title: Are Zonal Statistics possible in PostGIS2?

text: Now with postgis2 which supports rasters by default is it possible to do a zonal statistics analysis? I've google it but I haven't found anything solid? Is there a tutorial outthere to get me started? Can anyone give me an sql example on how to do it? **EDIT :** **An updated (simplified) query according theblog Aragon mention :**               CREATE TABLE sum_pop3 AS       SELECT gid, SUM((ST_SummaryStats(ST_Clip(rast,1,geom))).sum)      FROM perez_grid, ls_den      WHERE ST_Intersects(geom,rast)      GROUP BY gid;


title: Draw a circle in ArcGIS map

text: I have Center point(Latitude,longitude) for a place and radius. How to draw a circle/polygon on an ArcGIS map using Latlong and radius(also need WKT string). I am new to ArcGIS & do not have much idea about ArcGIS. Thanks in Advance.


title: How to reverse-georeference with QGIS?

text: I have a map with a "strange" projection (it's basically hand-painted, not an exact one). I can easily georeference it in QGIS using the Georeferencer plugin and the Thin Plate Spline transformation with a set of some 2000-3000 control points into simple WGS 84 (EPSG:4326) "latlong" projection and process/enhance it further with other data, like from Natural Earth or OpenStreetMap. Now I'd like to perform the _inverse_ operation: Transform one of those raster layers (or rasterised vector layers) _back_ into the hand-painted quasi- projection, ideally using the same set of control points, so that the images match when overlaid in a graphic editor. How do I do it?


title: How does latitude of origin work in Mercator

text: I'm trying to make a coordinate transformation using GDAL/osr python bindings. I have a longitude an a latitude and I want to transform them into the mercator coordinates. So I define the mercator projection as:               srs = osr.SpatialReference()     srs.SetMercator(41.602, 2.403, 1, 0, 0);      and the input point projection as               srs2 = osr.SpatialReference()     srs2.SetWellKnownGeogCS("EPSG:4326")      When I transform the point supposed to be 0,0 in the mercator projection, using               transf = CoordinateTransformation(srs2,srs)     transf.TransformPoint((2.403,41.602))      The x result is 0, as expected, but the y is a high number: 3799198 What am I doing wrong?


title: MSSQL Connection settings for QGIS

text: OK first lets get the blond elements out of the way- I'm not asking for how I connect to **_my_** actual database.   What I'm asking for is an example or screen shot of the entries required **_when setting up a new connection_**. Nothing I try to do to recreate either my 32 or 64 bit ODBC's is seeming to work and test connection just returns > [Microsoft][ODBC Driver Manager] Data source name not found and no default > driver specified QODBC3: Unable to connect Something along the lines of![Create a New MSSQL Connection](http://i.stack.imgur.com/XpKGt.png) _Given that the help files at this point are less than helpful!_![help me](http://i.stack.imgur.com/lL3vb.png)   And I've scoured the documentation for instructions in Lisboa, Dufour and Valmeira and the web all to no avail.   It should be so simple but I can't find anything- presumably because everyone thinks it should be so simple! Pretty please- yes call me a muppet if you want- but this is driving me stark staring bonkers- Can any one help? Alternatively how the &$^% do I find out more about the QODBC3 Drivers? because it is rather irrelevant in this case as I've tried it from multiple configurations but nonetheless QGIS Affected 2.0, 2.2   SQL server: 2008 R2   User environment: Win7 64 and 32 bit tested


title: How to design class breaks for a time series of choropleth maps?

text: I'm working with a time series of choropleth maps showing the same indicator for the same administrative boundaries, but spanning period of nine years. **What would be the most appropriate strategy of designing class breaks for such a time series?** At the moment, I think the quintiles division would serve the purpose the best, since it allows to scan quickly for 20% of areas with highest/lowest values (for instance). That works for one map, however I''m not sure how to porperly choose a scheme for time series. Options that I came up with so far would be:   1. use class intervals created separately for each map - advantage being they reflect the data structure for each year; disatvantage - harder to compare the maps between years.   2. use class intervals of quintiles from the summary map for the whole time period (without yearly split of the underlying data) - advantage being easier comparison of the maps for each year; disadvantage - breaks might not be ideal in certain years.   3. use class intervals of quintiles taken from the whole distribution of the data - all areas+all years - advantages disadvantages will most likely be simmilar to the above scenario. What are your thoughts?


title: How to reduce the opacity of a Google Maps layer in CartoDB.js

text: I have a CartoDB.js-App where I load my vis using `cartodb.createVis` (I do not manually define Leaflet layers). How can I reduce the opacity of the Google Maps base layer? I tried the following, as Google Maps has the layer index 0.               var initializeMap = function() {         // map initialization         cartodb.createVis('map', VIZAPI, {                 shareable: false,                 title: false,                 description: false,                 search: true,                 tiles_loader: false,                 center_lat: 47,                 center_lon: 8,                 zoom: 9,                 cartodb_logo: false             })             .done(function(vis, layers) {                 // layer 0 is the base layer, layer 1 is cartodb layer                 // setInteraction is disabled by default                 window.layers = layers;                 layers[1].setInteraction(true);                 // layers[1] is "all"                 layers[1].on('featureOver', function(e, pos, latlng, data) {                     updateClubInfo(data.cartodb_id);                 });                 layers[1].on('featureClick', function(e, pos, latlng, data) {                     updateClubInfo(data.cartodb_id);                 });                 layers[0].setOpacity(0.3);                 // layers[1].getSubLayer(0).setCartoCSS()                 window.overlays = vis.getOverlays();             })              .error(function(err) {             console.log(err);         });     };      This doesn't work as `setOpacity` is not defined on the Google Maps layer. I also came up with the idea of creating a white, "empty" sublayer between my data layer and the Google Maps layer and set its opacity. But I don't know how to do that with the CartoDB API.


title: Is it possilbe to change SLD from client side?

text: How to change predefile SLD in client side? The reason I want this is I want to turn label on and off on client side.


title: How to accurately join data based on 1x1 degree grid-cell points to a shapefile

text: In the data file each row represents a country and a grid-cell at 1x1 degree resolution, and the given coordinates, in longitude and latitude, label the south-west corner of the grid-cell. The duplicates in terms of coordinate exist if two or more countries' boundaries pass through the same cell. I have tried to join points in the file to a world map shapefile through the Spatial Join geoprocessing tool, but this way creates some problems. For example, many points are lost during spacial join because they locate out of a country's boundary or out of land. Ideally, I want to create a new shapefile in which the point data are saved in a grid; only in land part of a grid if the grid covers oceans; and in part of a grid corresponding to countries if the grid covers several countries. Another reason I want to this is that I want to summarize the data at subnational level. If two regions or more share such a grid, I can split data saved in the grid by whatever algorithm. Thanks!!


title: Is it possible to publish a map service from an SD File purely over REST (no arcpy)?

text: I want to automate the publishing of map services, but without using the arcpy library. Is it possible to publish a map service from an .SD File purely over REST (no arcpy)?


title: Field calculator Python else statement question

text: The codeblock below is run on Field1 and calculates Field1 to equal Field2 if they are not equal. Else it returns Field1. My question is, does the else statement actually recalculate Field1 to equal itself or does it leave Field1 as is. Basically if Field1 is equal to Field2, I don't want Field1 to recalculate to itself.               def update(Field1,Field2):       if Field1 != Field2:         return Field2       else:         return Field1


title: Producing gradients between shapefiles and using original boundarys as isoclines (lines)

text: to produce a bathmertical map of a river bed I draw shapefiles showing the different depths based on satellite images taken during different tidal levels (deepest, 2nd deepest, 3rd deepest, shore etc). to produce a final map would like 2 things: \- draw isoline (lines along the boundaries of different shapefiles \- shade or make gradients from on shapefile into another that also covers little gaps between the individual shapefiles. How can I do this using QGIS/ grass?


title: How does one use leaflet layergroups with CartoDB layers?

text: I'm trying to use leaflet layergroups with a cartodb layer, ala: http://leafletjs.com/reference.html#layergroup but haven't been able to successfully do it. This is what I have so far:                   var map = new L.Map('map', {           zoomControl: true,           center: [41.460, -81.67],           zoom: 11         });              map.locate({setView: true, maxZoom: 18});              var background = L.tileLayer('http://{s}.tile.stamen.com/toner-lite/{z}/{x}/{y}.png', {           attribution: 'Stamen',           maxZoom: 15,           errorTileUrl: 'http://fashiontech.files.wordpress.com/2008/08/plaid_b_single.jpg?w=256'         })          .setZIndex(0),              cdb = cartodb.createLayer(map, 'http://cleveland-metroparks.cartodb.com/api/v2/viz/05af5a7c-668d-11e3-8914-cb7f69afe66f/viz.json')               .on('done', function(layer) {                layer.setInteraction(true);                layer.on('featureOver', function(e, pos, latlng, data) {             cartodb.log.log(e, pos, latlng, data);           });                layer.on('error', function(err) {             cartodb.log.log('error: ' + err);           });         }).on('error', function() {           cartodb.log.log("some error occurred");         }),              labels = L.tileLayer('http://{s}.tile.stamen.com/toner-labels/{z}/{x}/{y}.png', {           attribution: 'Stamen',           maxZoom: 18         })          .setZIndex(600);              var existing = L.layerGroup([background, cdb]);         var baseLayers = {             "Existing": existing         };                      L.control.layers(baseLayers).addTo(map);      This answer hints that it's possible, but my js-fu is really lacking: How does Leaflet layer ordering work with CartoDB layers?


title: ArcGlobe not sitting on Globe Surface

text: I have gone through the process of bringing a SketchUp model into ArcScene utilizing Edit Placement on the 3D Editor window. The model appears to have come in great. However, it is floating slightly above my base parcel layer. Is there a way to get the multipatch 3D layer not to float slightly above another layer? I have even selected Floating on Custom Surface Elevation 30 m and yet it remains hovering above the ground. Any help would be appreciated. ![enter image description here](http://i.stack.imgur.com/lHOmd.png)


title: How do I save as .osm in QGIS?

text: The directions in http://docs.qgis.org/2.2/en/docs/user_manual/osm/openstreetmap.html#saving- osm-data seem to have been written before the OSM plugin was merged into QGIS proper. I can locate no OSM plugin or toolbar, and the vector menu only has options for opening a .osm file in QGIS.


title: Set a RANGE of values as nodata in a .vrt?

text: Is it possible to set a **range** of values to nodata in a vrt file? How? I have some elevation models with incorrect values, down in the -15,000s. I could just clip the rasters along that edge, or apply some raster calculation to transform them to nodata, but that's a fair bit of monkeying around. It'd be nicer to just say "everything less than -100 is nodata" and not duplicate the 99.9% of the cells that are just fine to fix the bad ones. In reading the GDAL Virtual Format Tutorial I see a **ComplexSource** option to add _"a custom lookup table to transform the source values to the destination. ... following form:`<LUT>[src value 1]:[dest value 1],[src value 2]:[dest value 2],...</LUT>`_ but when I try that it transforms all values, not just the ones within the range. Min/max from source: `Computed Min/Max=-16111.000,2618.000`   VRT Lookup transform: `<LUT>-32767:-10,-1000:-1</LUT>`   Out min/max: `Computed Min/Max=-5.000,-1.000` Am I doing it wrong, or is what I'm attempting simply not supported? (by any VRT construct, not necessarily just LUT). Is there some other method of setting a range to nodata that doesn't require duplicating files to get there? Source data:http://www.viewfinderpanoramas.org/dem3/P07.zip (34MB, home page), any tile ending in `W142`.


title: Where can I get data for traffic signals, stop signs, and other sign assets/inventory?

text: I'm working on a project involving roads. The quality of the service would benefit from knowing if a given road is continuous or if it's broken up by stop signs, traffic signals, and other things that effect a motorist. I haven't found anyhing googling around online. I also contacted my local city government, but they only grant access to the data for city projects. I'm guessing most areas would have similar rules. Is there another source, public or commericial, that could help me out? Or am I out of luck? EDIT - I should note that I'm looking for data in the USA.


title: Authoring map with advanced gradient symbol for ArcGIS web service?

text: I need to prepare a web service which will display the gradient of polygon points values. It has to look just like in this picture: link to the picture (please ignore the arrows). Each polygon is constructed always by 9 points. The logical order of points on the polyline forming the polygon is constant. Values of points will differ. The angles and lengths of lines may differ. Basically, it will always form a trapeze. First I played with splitting the 'large' polygon into 4 small ones using its central point. In this way, I could adjust the angles of gradients defining small polygons so that they formed this large polygonal gradient. There were however two issues - connections of the small polygons were looking rather ugly and most important - as the values of points may differ, it only suited the default scenario presented in the picture.


title: How do we add a custom basemap in ArcGIS.Com that can be used in Javascript API?

text: It's possible to add a layer as a basemap in ArcGIS.com by checking the check box use as baselayer: ![Use as base Layer](http://i.stack.imgur.com/aSZCe.png) But when I consume the map using:               var mapDeferred = esri.arcgis.utils.createMap("<guid>", "map", {                     mapOptions: {                         slider: true,                         nav: false                     }                 });      the basemap don't show up, is there any rules that define what content can be used as a basemap?


title: Check Personal Geodatabase attribute against external MySQL database?

text: I have a Personal Geodatabase to which editors enter values from some field inputs. I would like some of the entries, eg Project No ( whose name is unique and has some some related attributes like contractor etc ) to be checked against an MySQL database and pop up properties of that Project number from External DB to be available for comparison during data entry.. This would help to avoid wrong Project name entry , double check whether its the Same project which is currently being entered, based on related properties. I had a conversation with the Developer of the External Web application from where we want to fetch the data - Its a PHP application having a MySQL database at backend. So basically I would need to have a secure Login from the ArcGIS editing Window/Session. Does this warrant a different approach?


title: How to read raster datasets from ESRI geodatabase with FME?

text: I have a file geodatabase (10.2) with gdb feature classes as well as geodatabase raster datasets. I am using the Esri Geodb API reader in FME Desktop 2014 to read feature types from this geodatabase. However none of the raster datasets stored in the geodatabase are being acknowledged by FME. **How can I read raster datasets from an Esri file geodatabase with FME?**


title: Problem inserting new features to file geodatabase feature class using ArcObjects 10 SDK

text: I'm trying to add new feature to file GDB feature class and keep getting this error: Objects in this class cannot be updated outside an edit session. Feature class has attachments enabled (exists relationships). Googled alot about this problem and still can't figure it out. FYI I'm developing a standalone application. My code:               Public Shared Sub CreateFeature(ByVal workspace As IWorkspace, ByVal featureClass As IFeatureClass, ByVal point As ArcGIS.Geometry.IPoint, ByVal data As String)          Dim workspaceEdit As IWorkspaceEdit = CType(workspace, IWorkspaceEdit)          workspaceEdit.StartEditing(True)          workspaceEdit.StartEditOperation()          Dim featureBuffer As IFeatureBuffer = featureClass.CreateFeatureBuffer()          Dim featureCursor As IFeatureCursor = featureClass.Insert(True) ' <- at this point error occurs         Dim featureOID As Object          featureBuffer.Value(featureBuffer.Fields.FindField("Data")) = data          featureBuffer.Shape = point          featureOID = featureCursor.InsertFeature(featureBuffer)          featureCursor.Flush()          workspaceEdit.StopEditOperation()          workspaceEdit.StopEditing(True)      End Sub


title: Proper updating of new memory layer geometry?

text: Current workflow for a toolbar button: -create memory layer -save memory layer to shapefile Then the user uses the built-in QGIS tools to add geometry and finally saves changes via the auto prompt (when done drawing shapes). When this shapefile is later opened in another QGIS session, no geometry features are present. I'm inclined to feel the feature addition must occur prior to the memory-to- shapefile save by way of signal slots? If this is the case, could someone share a simple link or how-to for such?


title: How can I "select by attributes" based on membership in py dictionary?

text: I am trying to write a script in arcpy that checks to see if elevation values in a point file are found in a dictionary, and if they are to produce a new point file with just those entries. Sounds simple enough, but for some reason SelectByAttributes keeps returning the full dataset and not selecting properly. Further, this runs so slowly that I'm convinced there's a quicker way to do this, perhaps by making an array and writing a new point file with InsertCursor. Here's a snippet of my dysfunctional code:               import arcpy, collections     #transects_pts is a point shpfile with fields ID and elev     transect_pts = r'C:\transect_pts.shp'     inputRows = arcpy.SearchCursor(transect_pts,'','','ID; elev','ID')     layer = arcpy.MakeLayer_management(transects_pts,"tlayer")     #tdict is a default dict containing {ID: [elev1,elev2,elev3]}     for row in inputRows:         for m in tdict:             if row.elev in tdict[m]:                 arcpy.SelectLayerByAttribute_management(layer,"ADD_TO_SELECTION")     arcpy.CopyFeatures_management(layer,"new_transect_pts")      If anyone has some insight, or a better way to do this, please help me!


title: What are my hosting options for an app that needs tilestashe along side a django app?

text: I'm building a GIS webapp that uses tilestache to serve MBTiles and a django app to serve geoJson. I'm having trouble understanding exactly what kind of setup I'd need to host this. I'm new to server admin and setup, so something heroku-ish would be nice, but I'd be willing to learn if the solution involved more work.


title: How can I update an MXD file programatically?

text: I'm using ArcGIS 9.3.1. product line and MXD files. At our customer we have multiple environment (DEV, QAS, PROD) which should be more or less identical except they use different datbases and a few other minor tweaks. They alo share the same MXD file published as Map Services. These MXD should be the same, save for the fact their layer point to different data source (database) but with the same name & fields. Now, to update map services we basically manually edit each MXD on all platform manually. This manual process is error prone and has caused trouble in the past. I'd like to edit the MXD document once and simply copy them on other platform. But then I'll need to pass them through some home made made custom tool to edit the data source. I did write such a tool, it read the MXD document, find the used SDF file (connection file) and update the MXD's connetcion properties with the value from the local version of the SDE files. It works well **but for 1 problem**. When I change a layer's data source, its fields loose their aliases and visibility flags, which are reset to the default. I'd like to read the original value before update, update and then set the new values (to the old value)(for field's alias and visibility) I try to cast my layers to ITableField, but I got FieldCount = 0!!! So I can't update anything!!! Below is my current code ( **UpdateSDEWorkspaces** is the method of interest) How could I set up the layer alias and visibility?               public static class MxdUtils     {         public static IEnumerable<IMap> GetMaps(this IMapDocument pMapDocument)         {             int N = pMapDocument.MapCount;             for (int i = 0; i < N; i++)                 yield return pMapDocument.get_Map(i);         }         public static IEnumerable<ILayer> GetLayers(this IMap pMap, UID uid, bool recursive)         {             IEnumLayer pEnumLayer = pMap.get_Layers(uid, recursive);             pEnumLayer.Reset();             for (ILayer pLayer = pEnumLayer.Next(); pLayer != null; pLayer = pEnumLayer.Next())                 yield return pLayer;         }         public static IEnumerable<ILayer> GetLayers(this IMapDocument pMapDocument, UID uid, bool recursive)         {             return                 from m in pMapDocument.GetMaps()                 from l in m.GetLayers(uid, recursive)                 select l;         }              /// <summary>         /// Will update each SDE feature layer to the local sde connection.         /// It will will loop through all layer, find a layer's sde connection, if the local find can be found,          /// connection properties will be updated from it.         /// </summary>         public static void UpdateSDEWorkspaces(string mxdPath) { UpdateSDEWorkspaces(mxdPath, null); }         public static void UpdateSDEWorkspaces(string mxdPath, Dictionary<string, IPropertySet> memory)         {             if (memory == null)                 memory = new Dictionary<string, IPropertySet>();                  IMapDocument pMapDocument = new MapDocumentClass();             pMapDocument.Open(mxdPath, "");                  foreach (var pMap in pMapDocument.GetMaps())             {                 var sdelayers =                     from pLayer in pMap.GetLayers(new UID { Value = "{40A9E885-5533-11d0-98BE-00805F7CED21}" }, true)                     let pFeatureLayer = pLayer as IFeatureLayer                     where pFeatureLayer != null                     let pDataLayer = pFeatureLayer as IDataLayer2                     let pDatasetName = pDataLayer.DataSourceName as IDatasetName                     where pDatasetName != null                     let pWorkspaceName = pDatasetName.WorkspaceName                     let pNewProp = GetConnectionProperties(memory, pWorkspaceName.PathName, mxdPath)                     where pNewProp != null                     select new                     {                         Layer = pLayer,                         DataLayer = pDataLayer,                         DSName = pDatasetName,                         WorkspaceName = pWorkspaceName,                         UpdatedProperties = pNewProp,                     };                      foreach (var item in sdelayers)                 {                     var tf = (ITableFields)item.Layer;                     int N = tf.FieldCount;                     // damn, why is N always 0?                          item.WorkspaceName.ConnectionProperties = item.UpdatedProperties;                 }                      pMapDocument.ReplaceContents(pMap as IMxdContents);             }             pMapDocument.Save(true, false);             pMapDocument.Close();         }              public static IPropertySet GetConnectionProperties(Dictionary<string, IPropertySet> memory, string sdeConnPath, string mxdPath)         {             var path = sdeConnPath;             if (!File.Exists(path))                 path = Path.Combine(Path.GetDirectoryName(mxdPath), Path.GetFileName(sdeConnPath));             if (!File.Exists(path))                 path = Path.Combine(Environment.CurrentDirectory, Path.GetFileName(sdeConnPath));             if (!File.Exists(path))                 return null;                  IPropertySet set;             if (memory.TryGetValue(sdeConnPath, out set))                 return set;                  Console.WriteLine("Reading " + path);             set = GetSDEFactory().ReadConnectionPropertiesFromFile(path);             if (set == null)                 return null;                  memory[sdeConnPath] = set;             memory[path] = set;             return set;         }              static SdeWorkspaceFactoryClass GetSDEFactory()         {             if(SDEDataSourceFactory == null)                 SDEDataSourceFactory = new SdeWorkspaceFactoryClass();             return SDEDataSourceFactory;         }         static SdeWorkspaceFactoryClass SDEDataSourceFactory;     }


I want to convert TIGER data               GEOGCS["GCS_North_American_1983",         DATUM["North_American_Datum_1983",             SPHEROID["GRS_1980",6378137,298.257222101]],         PRIMEM["Greenwich",0],         UNIT["Degree",0.017453292519943295]]      to WGS84. How does ogr2ogr determine which transformation to use and how big would the introduced error be? Study area is the county of Kauai in Hawaii, USA.


title: Adding feature to GeoServer vector layer using WFS-T in absence of fid index?

text: I want to add a feature to GeoServer vector layer with WFS-T. When I use saveStrategy.save(), an error returned that said: > java.io.IOException: current fid index is null, next must be called before > write(). How can I add a geometric feature to GeoServer vector layer without fill feature fid?


title: The "Add delimited text layer" button does not show (QGIS 1.8.0)

text: I want to create a point map layer on the basis of an CSV file with XY data, with the "Add Delimited Text Layer" button. However, this button does not show. Also, it is not possible to add an Excel (.xls) file via "Add vector layer". I'm working in a secondary education school. The computer expert of our school has installed QGIS 1.8.0. On my own PC and laptop (same version of QGIS: 1.8.0), the tool does show and works fine. But on the computers at my school, the tool is simply not available. Why is that? Is it caused by the security settings our our school network system? Can it be that the "Add Delimited Text Layer" tool automatically creates a temporary map layer, which is stored somewhere, and that the security settings of our school network does not allow this? Regards, Tim Favier


title: Why unable to edit ArcSDE from multiple computers on LAN using custom ArcGIS Engine application?

text: I cannot edit an ArcSDE database from multiple computers on a LAN network using a custom ArcGIS Engine application, I have tried both versioned and non- versioned feature classes, but they cannot be edited from more than one computer. Inside the custom ArcGIS Engine application, for every single operation (Create, Remove), I open and close an edit session. Any advice or experiences are welcome as why I may be unable to edit an ArcSDE database from multiple computers on a LAN network using a custom ArcGIS Engine application? I tried to edit the same feature classe from 2 computers in the same time , but when i closed the edit session in the second computer, it gimme this error `Your edits could not be saved because the target version has been updated since you started editing. These edits have been merged into your edit session. Review the changes and save again.` i'm editing feature classes like in this example :                   IWorkspaceEdit workspaceEdit = (IWorkspaceEdit)DB_Access.get_obj().ws;         workspaceEdit.StartEditing(true);         workspaceEdit.StartEditOperation();                  IFeatureClass fc = DB_Access.get_obj().ws.OpenFeatureClass("points");                  IFeature f = fc.CreateFeature();                  f.set_Value(fc.Fields.FindField("name"), "name_of_point");                  f.Shape = geometry;             f.Store();                  workspaceEdit.StopEditOperation();                  workspaceEdit.StopEditing(true);


title: How to get public access to edit a layer in geoserver

text: I want to develop a GIS system with OpenLayers and GeoServer which is installed with OpenGeo Suite 4.1 and I need to edit a layer with no GeoServer authentication. This means that guest members should have permission to edit a layer. How can I do this?


title: problem with the installation of web ADF runtime arcgis server 10

text: I'm trying to install ArcGis server 10 (windows 7 32-bit)and I found problem in the installation of ArcGIS Server for the Microsoft .net Framework - Web ADF Runtime. The error is:                   ArcGis Server 10 - Web ADF Runtime cannot be installed on your machine. web adf runtime requires microsoft internet information services which has not been found on your system.      Thank you .


title: Interpolate shape (3d analyst) problem

text: I am making a code on C# to interpolate z values of a line on a DEM. Below is the code that I've written. I have also tried step by step debugging, the program works fine till the Geoprocessor execution. But on Geoprocessor execution, it crashes the ArcGIS. What could be the problem::: IMxDocument CurrentDocument = (ESRI.ArcGIS.ArcMapUI.IMxDocument)ArcMap.Document; IMap MapSelected = CurrentDocument.FocusMap;                       //The code below is to fetch the Raster Layer having name equal to what selected in combo box             int LayerIndex = 0;             string LayerName;             string comboText = (string)comboDEM.SelectedItem;                  for (int i = 0; i < MapSelected.LayerCount; i++)             {                 LayerName = MapSelected.get_Layer(i).Name;                      if (LayerName == comboText)                 {                     LayerIndex = i;                     break;                 }             }                  Geoprocessor GP = new Geoprocessor();             ILayer RasterLayer = MapSelected.get_Layer(LayerIndex);                  InterpolateShape Interpolator = new InterpolateShape();             Interpolator.in_surface = RasterLayer;                  //The code below is to fetch the Feature Layer having name equal to what is given in textbox             string XS_Line = txtInputCrossSectionLine.Text;             for (int i = 0; i < MapSelected.LayerCount; i++)             {                 LayerName = MapSelected.get_Layer(i).Name;                      if (LayerName == XS_Line)                 {                     LayerIndex = i;                     break;                 }             }                  ILayer FeatureLayer = MapSelected.get_Layer(LayerIndex);             Interpolator.in_feature_class = FeatureLayer;             Interpolator.out_feature_class = @"E:\For Arc GIS\Data\ResultFromTIN";                  GP.Execute(Interpolator, null);


title: How can I merge/collapse nearby and parallel road lines (eg. a dual-carriageway) into one single line?

text: I have a shapefile from OSM which contains all the roads within a city. Larger roads (like dual-carriageways) seem to have 2 parallel lines. I'd like to collapse / combine those into 1 line so I can export and render them nicely in an external 3D application (2 overlapping roads will be complicated and will look weird in 3D). How can I achieve this with either QGIS or PostGIS? I am not worried about a slight loss of accuracy (within a few metres) and would like the resulting, single line to be between (ideally the centre of) the current, parallel lines. Thank you. (here is an example of the dual road lines that I want to combine) ![enter image description here](http://i.stack.imgur.com/Ki460.png)


title: How to get WMS layer title instead of the layer name in QGIS?

text: When I connect to a WMS service in QGIS, the name displayed in the _layers panel_ is the layer `name`. Can I configure QGIS to always display the layer `title` for WMS layers and fallback on the `name` if the `title` is unavailable?


title: How can I clip this shapefile in Arcmap?

text: I'm trying to clip one shapefile with another in ArcMap 10.1, but every time I try I get an error saying "overlay not polygon". I'm pretty new to GIS - is clip the right tool to be using here? I've also tried intersection (as suggested in other threads) but that doesn't work either. ![Here's a picture](http://i.stack.imgur.com/MxUjP.png)


title: Maximum length of WHERE clause with IN statement

text: I need to get all features from feature class where their field value matches to input array. To achive this one may set IQueryFilter.WhereClause property similar to 'FIELD_NAME IN (array[0], array[1], array[2],...., array[n-1])'. But the size of array might be quite big. The questions are how long IN statement can be? And can I define it at runtime? Thanks in advance!


title: How to Clip a Raster dataset so that it is smooth along a boundary

text: I know how to clip a raster, however I need it to be smooth along the edge of a boundary and not jagged as can be seen in the image. ![enter image description here](http://i.stack.imgur.com/mYI7v.png)


title: Encoding in ESRI, unable to encode with UTF-8 when publishing .pmf files

text: I have a map currently made from ArcMap on a Windows OS. Now I would like to publish that map and then allowing it to open/read at a Linux OS (I'm using ArcReader to read on Linux). Everything is fine, the shapes and buildings and roads are all displaying on the Linux OS. But I have some Chinese characters used as layer names that have turned into strange symbols. I would like to keep the chinese character and I suspect it is a encoding problem when publishing the map. The Linux OS is able to display chinese in a .txt (encoding: UTF-8) and also display file names in chinese. I have also tried to used this and publish the map again. But the chinese failed to display also. Are there any ways to solve the chinese display problem in the .pmf file? Thanks!


title: How import a shapefile on-the-fly and have GeoServer render it instantly?

text: I am tasked with adding one feature to an existing app. It was built with OpenLayers and GeoServer. The task is to import a shapefile on-the-fly and have GeoServer render it instantly. However, I can't find any doc or tutorial. Since a user can import any shapefile, the geometry type or the data structure can't be determined ahead of time. Thanks, Fred


title: Modifying Bivariate Moran's I /LISA to include self?

text: I am using GeoDa's Bivariate Moran's I function to examine patterns of collocation between two variables at the county level (call them A and B for now). I know that, having run the analysis I can say things like "Counties with a high value for A tend to have neighbors with a high value of B" What I would rather say is that "Counties with a high value of A tend to locate in regions with a high value of B" The difference being that the latter would seem to include the county where A is high as well as its neighbors, while the former includes only the neighbors. I can see my way towards changing this--simply making every county its own neighbor in the weights matrix, but I was wondering what kind of havoc this would play on my values for I and the rest of the diagnostics. PS. If you happen to stumble on this identical question posted by me here: https://groups.google.com/forum/?fromgroups#!topic/openspace-list/WUL1kQkenWo please note that the answer I received was incorrect as was carefully pointed out by a rather irate reviewer.


title: How to find the attributes of a shapefile from the command line?

text: I'm following Mike Bostock's tutorial "Let's Make A Map." He explains how to filter a shapefile using `ogr2ogr` based on the attributes of the features. I can open a shapefile in QGIS and look at the attribute table to find out what those features are. Is there a way to find that information from the command line?


title: Update Arcpy version?

text: This is probably the answer to my previous question. The version of arcpy that I am using does not include GetNAClassNames > naClasses = arcpy.na.GetNAClassNames(serviceAreaLayer) **AttributeError** : `'module' object has no attribute 'GetClassNames'` According to this help topic, GetNAClassNames is in ArcGIS 10.1. Why is it not in my version? Thank you.


title: Help with Point on Polygon overlay in ArcPy

text: I am trying to do this simple task but not sure how to do this in Python. I have a point shapefile and I want to overlay with polygon layer to get all the polygons for each point. After this would like to loop through all the rectrods to get the attributes. Appreciate any feedback.               wellRows = arcpy.UpdateCursor("H:\points.shp")           for wellRow in wellRows:         wellPoint = wellRow.SHAPE #Not sure how to have shape as input...         landRows = arcpy.SearchCursor("H:\Poly.shp", "", "", "","")              for landRow in landRows:             print landRow.NAME


title: Bilateration using longitude, latitude

text: I have 2 coordinates (Longitude and Latitude) and a distance from them to a third coordinate. How is it possible to calculate the longitude and latitude (probably 2 possibilities) of the third coordinate? I have found these 2 answers: How to calculate the intersection of 2 Circles? How do I find the intersections of 2 circles on earths surface? But did not understand them as I have finished my geometric lessons long ago :) I would be glad if someone could simplify it for me so I can find the third coordinate in an easy way. Thanks!


title: combine raster and polygon values in R

text: My question is similar to this one Get Raster Values from a Polygon Overlay in Opensource GIS Solutions but I think I need another step. I have a polygon layer of ecoregions which I brought into R with "readOGR". It has many attributes such as ecoregion and biome (biomes encompass multiple ecoregion polygons). I have a raster layer (I used "raster" on a tif) that is 0-2. I also have continuous raster layers, that I'll use later so I'm looking for a generalized solution. I would like to be able to do a variety of evaluations like total area and proportion of each biome >0 or when I use the continuous raster, things like mean or range. Ideally I'd like to end up with some sort of attribute table that I can work with that combines the polygon attributes with the raster values. I'm not sure if it's most accurate & efficient to convert both to polygons, both to rasters, or to do summaries as is. I am trying to work in R so that I can script the process; sample code is particularly helpful as I'm a novice. I appreciate any thoughts.   Thanks a lot


title: How to copy a fgdb using Python?

text: I am trying to copy fgdb from one path to another using Python.               copy(base.gdb,dest)      And got the message:               Error Info:      [Errno 13] Permission denied:      How to fix it?


title: Cartographic Representation: how to use different rule for different polygon?

text: I create a file geodatabase to hold a polygon of the land boundary. What I would like to set up the RULE to separate colors. I had to split the polygons into 2 separate polygons and change the ID to 2. How do I add that to the Resprentation Symbology properties ? I already have set up for the north part but need to change the color the south part. Thanks in advance ! here is the picture I want to change the colors. ![enter image description here](http://i.stack.imgur.com/tRdwY.jpg)


title: Is there a tool that allows the same functionality of extend/trim feature that was available in 9.3?

text: In ArcGIS 9.3 there was a tool called Extend/Trim Features that allowed you to select multiple segments and extend or trim them to a feature. In 10.1 that seems to have been moved to the Production Mapping Extension at a hefty price. It was replaced by two tools that don't seem to have the same functionality as the previous versions tool. My question: Is anybody aware of an add in or way that I could create this functionality with python or a custom tool? **EDIT** The old tool can be found referenced HERE, there is also a screenshot of the tool location from 9.3: ![enter image description here](http://i.stack.imgur.com/W6M6N.jpg) Essentially, it did both of the functions that are shown in @Jakub's answer. It looks like the same tool from 9.3 was moved to the Production Mapping toolbar in 10.1. There seems to be an ESRI Idea out there suggesting they bring it back, and that I'm not, in fact, imagining that it used to be this way.


title: Why am I unable to save shapefiles layers? (other layers work fine)

text: I'm a very new user on QGIS. I've managed to find mostly all I needed to use the software, however today I came across some problems similar to those encounted by [someone else][1]. However, when I tried the suggestions in the answers, nothing worked. I added a new shapefile layer to my map, I set the CRS to WGS 84 / Pseudo Mercador EPSG 3857, just like every other layers in the map I am trying to create. I used the Openlayers plugin to load the Google physical basemap. I tried adding it before the other layers or after the other layers. No matter what I do, the basemap always shows up without issues. In fact, my issue is with my shapefile. I use the 'New Shapefile layer' icon, then I choose 'line' and set the CRS right away, and start drawing my line. I pick the color (red), the type of line (dashed) and the size, and everything looks fine. The problem arises when I save the project (or the layers independantly) and I close it. When I open it to work on it again, the basemap is there as well as all the points I created using the 'Add delimited text layer' icon with a .csv file. In the right hand side menu, the shapefile layer for my lines appears on top of the list, the box is checked, but the line color is set to a default purple, but it doesn't show up on the map. It's like my shapefile has not been saved, or that it has been saved before I started working on it. I tried changing the CRS from WGS 84 / Pseudo Mercador to other settings, but nothing works. I can still work on the layer, add to it, but if I try and save it again, I still get the issue when I open it to keep working on it. In the project properties windows, I have tried enabeling and disabeling the 'on the fly' CRS transformations, but it didn't work. I am working on OS X Mountain Lion with QGIS 1.8.0-Lisboa. Thank you in advance for your help!


title: Tilemill reversing land and water at high zoom levels

text: I tried asking this over at MapBox's support forum and got no help.   I am trying to generate a map of Northern france using Tilemill and this set of land polygon data from openstreetmapdata.com (the first download link).   However, at zoom levels 17 the water and land invert colors and styling (see attached output from Tilemill). Any idea of what's going on? Is the data bad? Can someone help me find an alternative? Thanks!! Zoom level 16 Zoom level 17


title: How to calculate a new raster based on nodata values from original raster, in QGIS?

text: I'm using the raster calculator from core in QGIS and I need to generate a new tif raster file (mask) from the nodata values of another tif raster file. I tried isnull(), but no luck. Any one can help? Thanks


title: Estimating cline width at multiple points along a trait gradient

text: In the map below I have 66 (green) points along the midpoint (0.5) for an allele frequency cline across the region (see map). I created the points by placing a point every 25 km along the 0.5 allele frequency contour. The thick red lines represent the 0.22 and 0.78 allele frequency contours. I would like to estimate the width of the cline at each of the 66 points. The "estimated" width should be based on the shortest distance between the two contours when the line runs through the point. This would allow me to compare the cline width at 25 km increments. ![enter image description here](http://i.stack.imgur.com/IyKN5.png) The result would look something like the map below. Note that I hand-drew only of 3 of the 66 desired lines and my lines are probably **not** the shortest distance between the two contours for each of the 3 points. Each of the line must be a straight and unbroken (i.e. no bend at the point). I would also get the total distance of the line, plus the distance for each segment (0.22 to 0.5 and 0.5 to 0.78). ![enter image description here](http://i.stack.imgur.com/2VXMp.png) This is a map created using the "Near" function with my points as the "input feature" and contours as the "near feature". I created the map by running the near function for each contour separately. As you can see the map is not quite right. Any other suggestions? ![enter image description here](http://i.stack.imgur.com/x2WM5.png)


title: raster data expressed as percent

text: I downloaded NDVI greeness data from the USGS Integrated Vegetation Mapping Viewer. In the metadata it says that visual greenness values range from 0 to 100. My data is from the period ending 3/5/12 in the NYC area my data table values only range from 0 to 68 which makes sense since it's winter. I need to display the data in terms of percent, which seems easy enough, but has me hung up. When I load the raster into the dataframe you can see that the values in the table of contents range from 105 to 255 (the BinValues). I think I am having trouble understanding the actual data within the raster and how to manipulate it in order to obtain %. The data is also in GCS (NAD83) and I will have to project it into Albers Equal Area to match my other datasets assigning the cell size to that of another dataset as well (30m) so I have to use a bilinear sampling. I am a newbie at this so I hope this makes sense. Thanks. ![raster](http://i.stack.imgur.com/DIlrO.png)


title: ArcGIS JSAPI selectFeatures on feature layer with "stacked" features only returns one feature

text: I am using the selectFeatures function to launch a popup when clicking on a feature. Because of the nature of my data, many of the features (street segments) have coincident geometry (they are stacked on top of each other). When I click on one of these segments that I know has multiple stacked features, the selectFeatures function only returns one feature (assuming the "top" one). If the features are adjacent (within my 10 pixel click buffer) but not coincident, it does return multiple features. How can I have the selectFeatures return all features at that location, not just the one?               dojo.connect(map, "onClick", function (evt) {         var query = new esri.tasks.Query();         query.geometry = pointToExtent(map, evt.mapPoint, 10); //buffers click point by number of pixels(10)              var deferred = featureLayerSeg.selectFeatures(query, esri.layers.FeatureLayer.SELECTION_NEW);         deferred.addCallback(function (features) {              console.log(features.length); //this returns 1 if the features are stacked on top of each other         });     }      PS- I don't want to use IdentifyTask to a map service because these feature layers have already been filtered by the user based on search parameters.


title: Finding the Nearest Landmass to A Point In The Ocean via API

text: I'm building an iPhone app that often ends up in a point that is in the middle of the ocean. I want to be able to find the nearest landmasses, which in some casses are thousands of miles away. So far I've thought of calling the maps API on the regions NW, NE, SW, SE - i.e. the ocean coordinate +10 lat, +10 long etc. until I come up with a response that isn't ZERO_RESULTS (I could obviously keep expanding the bounding box until I found something). I'm currently trying to find info on the antipode to San Francisco, which is approximately -37, 56. The clear limitation on this is that it is in no way scalable, since Google has a request limit of 2,500 per day and I would probably be making 4+ of them per entry by a user (and I expect a user to make at least a few of these requests per use). Any thoughts? Thanks, Alex


title: Query from plugin that shows result on Map and in Attribute Table

text: I'd like to create a plugin that performs a static (no user input required) query that, based on a condition, selects all the features on the map and displays these in the attribute table. I have the creation of the plugin itself under control. However. 1) Querying the PostgreSQL DB with PyQGIS from the plugin and, 2) based on the result, select the features on the map (I'm not out to create a new layer), 3) and showing the selected features in the Attribute Table, is where I could use some hints. Anyone that has experience with this or knows of one or more plugins where I can get some really good inspiration?


title: Shapes of this entity type not allowed in this layer when CreateFeature() is called

text: I have the following code that calls on a remote geodatabase server. Everything is connected and working, but when I try to add a point and use featureclass.CreateFeature();, it shuts down. At first I realized it was esriMultiPoint, but I fixed this, I remade the command lines, and created the shp2sde again. And this time I can confirm the server says it is esriGeometryPoint.               IWorkspaceEdit wse = ws as IWorkspaceEdit;     wse.StartEditing(true);     wse.StartEditOperation();     Log.Write("Feature class=" + fClass.FeatureType.ToString() + " Fieldtype=" +fClass.ShapeType.ToString());     //Get a reference to the new point object     ESRI.ArcGIS.Geometry.IPoint com_point = (ESRI.ArcGIS.Geometry.IPoint)(serverContext.CreateObject("esriGeometry.Point"));     Log.Write("Create esriGeometryPoint");     ((ESRI.ArcGIS.Geometry.IGeometry)com_point).SpatialReference = ((IGeoDataset)fClass).SpatialReference;     IGeometry geometry_point = com_point as IGeometry;     Log.Write("Check geometry point type=" + geometry_point.GeometryType);     //input coords      object tempObj;     paraMap.TryGetValue("Longi", out tempObj);     com_point.X = double.Parse(tempObj.ToString());     paraMap.TryGetValue("Lati", out tempObj);     com_point.Y = double.Parse(tempObj.ToString());     Log.Write("about to create feature");     IFeature f = fClass.CreateFeature();     Log.Write("Create feature"); //CODE NEVER GETS HERE     f.Shape = com_point;      Here's my **commandline** to create the table in the database... (I put dots for excluding some info).               shp2sde -o create -l mynewdata3,shape -f C:\..bla.shp -a all -i sde:sqlserver:..... -C FID,SDE,1 -s ...... -D ...... -e p -u ..... -p ......      **EDIT, should have used:** shp2sde -o create -l mynewdata3,shape -f C:..bla.shp -a all -i sde:sqlserver:..... -C FID,SDE,1 -s ...... -D ...... -e pn -u ..... -p ...... **The log output and ERROR:**               Feature class=esriFTSimple Fieldtype=esriGeometryPoint     Create esriGeometryPoint     Check geometry point type=esriGeometryPoint     about to create feature     at ESRI.ArcGIS.Geodatabase.IFeatureClass.CreateFeature()\r\n Shapes of this entity type not allowed in this layer      **EDIT:** I checked with logs and saw that the feature class being used, is correct feature class too. **EDIT 2:** Solved Details about the error: http://support.esri.com/en/knowledgebase/techarticles/detail/16135 Solution: http://forums.esri.com/Thread.asp?c=159&f=1707&t=232881


title: How to find locations on a conditional search

text: Hope you can help I have a set of points evenly distributed to 500m intervals, and stored elevation. I want to find from a given location any other point with the lowest elevation within a 1,000m search from the originating location. Do this again to the following 500m away point strictly searching for points with lower elevation to avoid including the previous searching location.


title: An algorithm for locating candidate labelling boxes within a polygon

text: Has anyone implemented this algorithm In PostGIS? Quoting from the paper: > _Vector-based geographic information systems usually require annotation, > such as a polygon number or attribute data, in a suitable location within a > polygon. Traditional methods usually compute the polygon centroid, test the > centroid for inclusion or exclusion, and select some alternative point when > the centroid falls outside the polygon. Two problems are associated with > this approach:(1) the text can be centered on the point, but may be placed > in a visually awkward place, and (2) part of the text may fall outside the > polygon and may overlap other polygon boundaries or other text labels. An > algorithm is presented that circumvents both of these problems, by computing > a number of horizontal candidate labelling rectangles (boxes) within a > polygon from which a suitable selection can be made or from which one may > conclude that the text label does not fit the polygon._ Perhaps in a custom build, or as a combination of several `ST_` functions? Even a command line utility?


title: looking for datum to georeference maps of 19th century Maine

text: I am trying to georeference 3 (1840, 1844, and 1877) northern Maine maps. I really am not able to find which datum to use. Can anyone point me to a list of historic and regional datums? thanks


title: multiple users on one wms layer - style and filter

text: I am using geoserver to show Points data from MySql in OL. The data saved in one table for all the users - 'Points' with a geometry field and Id of the user. I want to be able to let the users to see only their data, and also to style this features, but i didnt find any solution to filter or styling this data **not** from the geoserver panel. I need to filter the data by the Id number of the user, and then he see only his data, it is possible to do it like this? from the js file? if not, there is a other way to do it?


title: What programming languages are used by different GIS platforms?

text: Suppose a programmer has no knowledge of GIS, but wants to leverage his/her programming skills to get into GIS. Please list different GIS platforms, packages, etc. that use a particular language. (I was thinking, as suggested, about doing a simple test of the code formatting capability of StackExchange. I figure this might be done in the context of a useful answer.) When you answer, if the language is one of those listed below, please post a code snippet and test the formatting to make sure it works. Here's an excerpt from Jeff Atwood's answer on meta-stackoverflow. The spec is:               <!-- language: lang-or-tag-here -->              code goes here          <!-- language: lang-or-tag-here -->              code goes here      You may use either a **tag** or a **prettify language code** to specify, though prettify language codes are always guaranteed to work regardless of what language the tag happens to be set to. Use `<!-- language-all: lang-or-tag -->` to use the given highlighting for all the following code blocks. The complete list of supported prettify code languages is:   * lang-none   * lang-default   * lang-bsh   * lang-c   * lang-cpp   * lang-cs   * lang-csh   * lang-css   * lang-hs   * lang-html   * lang-java   * lang-js   * lang-lisp   * lang-lua   * lang-ml   * lang-perl   * lang-php   * lang-py   * lang-proto   * lang-rb   * lang-scala   * lang-sh   * lang-sql   * lang-vb   * lang-xml


title: GPS device taking pictures and saving heading + HDOP?

text: I'm looking for a device that can do all of the following:   * Take picture   * Save lat-long GPS position + _heading_ (angle to North) of picture [so need tilt-compensated 3D compass], in exif or other format (save tilt as well if possible)   * Save _precision-related GPS (NMEA) fields_ , as many as possible (at least HDOP, other if possible)   * If possible: display OSM maps (better if I can put my own tiles) Ideally HDOP and the like should go into both exif of the pictures and gpx track, but post-processing (gpscorrelate or the like) is OK. I would need to know this ASAP to buy the device and learn how to use it: it's for an 80-day alpine OSM mapping party, starting early June 2012! See here: http://www.outdoormaps.org http://wiki.openstreetmap.org/wiki/ViaAlpina There are already some data on this long high-distance trail, but they lack accuracy estimates. Heading will make the pictures much more valuable. Thanks! Mayeul Kauffmann (PS: I've seen this question: Logging GPS and Accelerometer Data - Device Suggestion? but I need heading and HDOP)


title: How to retrieve the year from a date format?

text: I'm attempting to retrieve the year from the following date format: 1/1/1988 I've tried the following but it does not work str(1/1/1988)[4:] and I've also tried str(1//1//1988)[4:]


title: How to set output resolution for gdalwarp to that of source file?

text: I have a large raster, which I want to re-project and crop for a new location that I'm creating in GRASS. I seem to have gotten the re-projection and cropping right with this call to gdalwarp:               gdalwarp -t_srs EPSG:32142 \       NE1_HR_LC_SR_W.tif \       -te -1246270.9107927915 -186894.1315662132 2391103.024689787 2372883.4615327567 \       NE1_HR_LC_SR_W_32142.tif      Unfortunately, the resolution of the target file is much lower than that of the source file. I see in the gdalwarp documetation that the argument `-tr` specifies the target file's resolution, but it is not clear to me how to specify that the target file's resolution should be that of the source file.


title: Configuring output table popup for geoprocessing widget

text: I have built a simple geoprocessing widget that produces a table (record set) as its final product. I am trying to simply change how and which fields are displayed. I am using Flex Viewer 3.5 and Arc Server 10.1. Here is my output: ![enter image description here](http://i.stack.imgur.com/RIJSW.jpg) Basically I want to only show the Value and Acres field. I tried doing this through the flex viewer in while configuring the GP Widget, but it didn't have any popup configuration options. ![enter image description here](http://i.stack.imgur.com/G2iQC.jpg) Any ideas of where I can set these configurations? GP.xml perhaps? Thanks in advanced. Adam


title: set the z-index of an openlayers map

text: I have an openlayers map and I try to implement a "lightbox" image gallery , in the same page. By ""lightbox" image gallery" I mean that the images/arrows/exit button cover the whole screen. The problem is that openlayers map covers the lightbox gallery (see image). I have set `z-index` to the map and to the lightbox divs, like so               <div id='map_element' style='width: 900px; height: 400px; z-index:1;'> </div>          //and then for lightbox     #imagegal2 img {      height:85%;          z-index:90;         }      And this does not work How do I fix this? Thanks in advance ![enter image description here](http://i.stack.imgur.com/H5KFU.jpg)


title: How to use coding for auto-comparing two different fields in ArcGIS for Desktop?

text: I want to compare two group fields in attribute table, which are related to the name of villages. The name of these villages (in both field) are the same but they have in some cases differentiation (because of misspelling). In order to update my databases I need to find different names and edit them. They are many names and it doesn't make sense to check them manually. I think Python or VBA code in ArcGIS is practical. How can this be done?


title: How to resample GeoTIFF images to the same resolution?

text: I have two geotiff images that I would like to resample the same resolution. If I use ArcGIS this would be easy to do. But using QGIS I don't know how to do this. Some help will be very appreciate. Thanks.


title: ArcGIS 10.0 get sde information out of .mxd

text: The overall goal here is to iterate through several directories of .mxd files and document which SQL Server instance of SDE they are pointing to, with the database name and feature class names. All of our connections are direct connect, no SDE services running. I'm trying to get the server/database info out of an mxd using ESRI's example here http://support.esri.com/cn/knowledgebase/techarticles/detail/41112, but it keeps trying to look at the .sde file that was used to create the mxd instead of getting the info straight from the .mxd. specifically, this line               desc = arcpy.Describe(sub_str)      throws an error that it can't find the .sde file. It's not going to since we have many users all with their own .sde connection files. The sub_str ends up being C:\Users\username\etc... depending on who created the .mxd. Is there another way to get the datasource info out of an .mxd that doesn't rely on the .sde file? (ArcGIS 10.0) EDIT: SOLUTION USED On the off chance my horrible python skills can be of use to someone, this is the solution I ended up with. I know it can be cleaned up, but I just needed quick and dirty feel free to critique:               def mxdInfo(inMxd, inCsv):          mxd = arcpy.mapping.MapDocument(inMxd)          for lyr in arcpy.mapping.ListLayers(mxd):         if lyr.supports("SERVICEPROPERTIES"):             servProp = lyr.serviceProperties             lName = lyr.name             dbName = servProp.get('Database', 'N/A')             serviceName = servProp.get('Service', 'N/A')             serverName = serviceName.split(':')             userName = servProp.get('UserName', 'N/A')         else:             print "SERVICEPROPERTIES not supported on " + lyr.name             continue         if lyr.supports("DATASOURCE"):             dSource = str(lyr.dataSource)             fcName = dSource.split('.sde\\')             print "Layer Name: " + lyr.name             print "Feature Class: " + fcName[1]         else:             print "DATASOURCE not supported on " + lyr.name             continue         #write out to csv         if not os.path.isfile(inCsv):             csvFile = open(inCsv, 'wb')             try:                 writer = csv.writer(csvFile)                 writer.writerow(('MXD', 'Layer Name', 'Server', 'Database', 'Feature Class', 'User Name'))                 writer.writerow((inMxd, lName, serverName[2], dbName, fcName[1], userName))             except:                 print "error writing first row of csv"         else:             csvFile = open(inCsv, 'ab')             try:                 writer = csv.writer(csvFile)                 writer.writerow((inMxd, lName, serverName[2], dbName, fcName[1], userName))             except Exception as e:                 print "error writing to csv"                 print e     del mxd


title: Trimble RTK Data to ArcGIS

text: I wondering if there is a way I can take my field data collected on a Trimble R10 RTK GPS and import it into ArcMap for further processing, I have done some research and I know of Trimble GPS Pathfinder Software, is this the only way? If there is a way to bring it into Quantum, that would also be alternatively helpful. Thanks Always


title: Create a new Dataframe using Arcpy in 10.1

text: Does anyone know if there is a way to create a new empty dataframe in ArcMap 10.1 using Arcpy? Thanks


title: How to manage null values in QGIS

text: I made a merge of several AsterGdem blocks near Egypt and Israel. I have portions of the Dead and Mediterranean Seas. Since near Galilea elevations are below sea level this makes the following problem When I define the null value as 0, the resulting image shows the seas as null values, which is correct, nevertheless, I also have some pixels inland forming a kind of a contour with the null value as 0. What I want to accomplish is to have the seas as null values and inland elevation information correct also Thanks before hand Gerardo


title: Looking for a tool for Orthorectification of Archived Air Photos

text: I was wondering if anyone knows of an open source, perferabily something for QGIS and not GRASS, that can be used for orthorectifying scanned archived air photos. I have looked at OSSIM and a few other packages, but they do not allow (from what I can tell) for a single image to be processed based on user input from a geoimage or vectors. Packages like ERDAS and PCI allow for this, but I am looking for an open source solution. The Georefrencer plugin to QGIS is good, but I would like to be able to use a DEM (digital elevation model) as a source for the Z information for terrain relief issues. Any help would be great, or clarification on OSSIM or other products that I may be misunderstanding. Thanks.


title: Using ArcPy update cursor to replace <Null> value?

text: I have been combing the help associated with (None) values in arc, and have not been able to solve the following problem: I have a a field with mixed Null and integer values: 0 100 I would like to transform the null values to 0's in python:                  cursor = arcpy.UpdateCursor(fishJoin, ['myField']):              for row in cursor:             if row is None:                 row = 0                 cursor.updateRow(row)      This does not alter any of the null values. I'm using ArcGIS Desktop 10.0.


title: Is it possible to take automatically a "screenshot" of my selected feature in QGIS using Python?

text: I am implementing a tool which allows the user to create a profile of a feature (in my case urban parcels). Is it possible to automatically take a screenshot of this feature (displayed in a QGIS window) and to display it in the PDF profile? The picture should only show the extent of the feature and maybe a WMS-layer that lies beyond the feature. Are there already some tools which you could suggest?


title: How to relate Slope map and Curvature map so as to have in the end PERCENTAGE PROBABILITY MAP?

text: i am trying to deal with road detection this last month. So i am working with particle filtering that produces particles, that means possible targets(to my case cars)/points, with hidden information of speed). Further on, i implement spatial+statistical analysis through ArcGIS, using the density of these particles/points as an indicator of where could a road possibly be (through a methodology i follow, is possible to extract a road using the centroids od the polygons - small areas of bigger density of points. I have the slope map and the curvature map of the area. I have also converted everything to rasters and what i would like to do is a percentage probability map that will inform the user that "in this area is impossible what you see to be a target - so it's a false alarm". This will happen to my case in an area that the slope is smaller than 30 degrees of ground declination. Are there any ideas of how could i RELATE SLOPE AND CURVATURE and later to proceed to PERCENTAGE PROBABILITY MAP? I am in a good path to achieve it , but i would appreciate any ideas!Thanks a lot!


title: How to calculate raster distance with z value using ArcGIS for Desktop?

text: I'm trying to produce a raster distance map by using euclidean distance function having input feature as flight path (with altitude as Z value). However, my euclidean distance calculation did not consider the altitude (z value of the input feature file). The distance was computed only from xy values. Does anyone know how to calculate distance using all xyz values?


title: Crossover Analysis (Intersection)

text: I have a dense set of flight lines - list of (lat, long) coordinates, and I am interested in determining the intersection (crossover) of the flight lines. I want to correlate the image with the intersection. I know this is a GIS problem. Do you have any suggestions on how to solve this problem? Thanks in advance


title: Cannot Add Address Locator in ArcMap 10.2

text: I have an address locator. I wish to load this locator into ArcMap 10.2 in order to geocode some old addresses. When I go to the geocoding toolbar, drop down the left-most menu, select "Manage Address Locators," select "Add," and then find my address locater and add it, I get the following message: > The address locator is damaged. Would you like to repair it? When I select "Yes," I get the following error message: > Error opening dialog: this address locator does not support the necessary > interfaces needed by this dialog. I am then unable to add the locator this way. I also tried to go to ArcToolbox, Geocoding Tools, Geocode Addresses, input my data, and then when I input the locator I get the following error message: > Error 00732 Input Address Locator: Dataset > C:\Users\person21\Folder\address_locator1 does not exist or is not > supported. I am then unable to add the locator this way. It seems I cannot add the locator at all. Why is this? The locator may have been "constructed" in an older version of ArcGIS. Is this the problem? Thank you for any help.


title: How to limit number of decimal digits in ArcGIS to have exported file in CAD with same format?

text: I am supposed to have some numeric values in ArcGIS 10.2 for Desktop with 3 decimals Digits and then convert them to Microstation or AutoCAD. I did this Setting using field calculator and 'round' function in ArcGIS but when I convert it to CAD again it would come with more than 10 decimals which is not desired. I would appreciate it if you kindly could assist me with this issue.


title: How to replace data source of annotation feature classes using ArcPy?

text: I want to batch replace paths in a series of `.mxd` files using Python. Many of the layers have moved and in some cases the name has changed or its dataset is different. I've written this code which takes an Excel file with two columns (one column for the old path and one column for the new path) and converts to a dictionary. Then I iterate through the layers in my `.mxd` files, and iterate through the dictionary, using `replaceDataSource` to change the data sources. The trouble is this won't work with annotation feature classes. I've added in a try/except with `findAndReplaceWorkspacePath` to see if that may work (I don't believe the annotation layers/names have changed), but the code still doesn't address those layers. One thing to note is that the annotation feature classes are stored in a personal geodatabase, but I don't know if this could be causing any issues. I am hoping someone has ideas on code that will work to change the path of the annotation feature classses! Here is my code:               import xlrd, arcpy, os          mxd_path = r'C:\temp\Chapters'          d = {}     wb = xlrd.open_workbook('U:\GIS\Chevron_Lokern\Projects\Figures_Update_2014\Lokern_Source_Paths_KL.xls')     sh = wb.sheet_by_index(0)     for i in range(293):       cell_value_class = sh.cell(i,0).value       cell_value_id = sh.cell(i,1).value       d[cell_value_class] = cell_value_id            for dirpath, dirnames, filenames in os.walk(mxd_path):         for filename in filenames:           fullPath = os.path.join(dirpath, filename)           if filename[-4:] == ".mxd":             mxd = arcpy.mapping.MapDocument(fullPath)             print filename                  for lyr in arcpy.mapping.ListLayers(mxd):                 if lyr.supports("DATASOURCE"):                     #print lyr.datasetName                     for k,v in d.iteritems():                         #print k,v                         k_path, k_name = os.path.split(k)                         v_path, v_name = os.path.split(v)                         if lyr.dataSource == k and k != v:                             #print 'KEY:',k,'VAL:',v, 'NAME:',v_name                             try:                                 print 'trying to change data source for',lyr.dataSource                                 lyr.replaceDataSource(v_path,'FILEGDB_WORKSPACE',v_name,True)                             except Exception as e:                                 print 'can\'t replace',lyr.dataSource, e                                 try:                                     print 'trying to change data source for',lyr.dataSource                                     lyr.replaceDataSource(v_path,'RASTER_WORKSPACE',v_name,True)                                 except Exception as e:                                     print 'still can\'t replace',lyr.dataSource, e                                     try:                                         print 'trying to change data source for',lyr.dataSource                                         lyr.replaceDataSource(v_path,'SHAPEFILE_WORKSPACE',v_name,True)                                     except Exception as e:                                         print 'still can\'t replace',lyr.dataSource, e                                         try:                                             print 'trying to change data source for',lyr.dataSource                                             lyr.findAndReplaceWorkspacePath(k_path,v_path)                                         except Exception as e:                                             print 'can\'t replace workspace path for',lyr.datasetName, e             mxd.relativePaths = 'True'             mxd.save()     del mxd


title: create polygon adjacency matrix using python

text: I am trying to create an adjacency matrix from a set of polygons. That is, I have a bunch of polygons and I want to identify which polygons have a common edge or "touch" each other. Once I find this information, I want to create an n x n matrix that indicates whether those each polygon either touches or does not touch the other polygon. Something kinda like this below.                     Poly1    Poly2   Poly3          Poly1   0        1       0     Poly2   1        0       1     Poly3   0        1       0      That is the goal anyway. I was wondering if anyone knew of a python package or such that could help me do this. I need to dynamically generate the polygons and then get their adjacency as I go, so using arcgis or qgis not the best solution. Thanks.


title: Temporal GRASS GIS framework

text: I am trying to follow the tutorial posted here. When I get to this part "Extraction of seasonal mean temperature" and run this code:               t.rast.extract \       input=temperature_mean_1950_2011_seasonal \       output=temperature_mean_1950_2011_spring \       where="start_time = date_trunc('year', start_time \              + interval '2 months')" \     expression="temperature_mean_1950_2011_seasonal/100.0" \       base=temp_mean_spring      I get the following error:               ERROR: Unable to get map ids from register table        <temperature_mean_1950_2011_seasonal_PERMANENT_raster_register>     Traceback (most recent call last):     File "/usr/local/grass-7.0.svn/scripts/t.rast.extract", line 91, in <module>     main()     File "/usr/local/grass-7.0.svn/scripts/t.rast.extract", line 85, in main     base, nprocs, register_null)     File "/usr/local/grass-7.0.svn/etc/python/grass/temporal/extract.py", line 95, in extract_dataset     rows = sp.get_registered_maps("id", where, "start_time", dbif)     File "/usr/local/grass-7.0.svn/etc/python/grass/temporal/abstract_space_time_dataset.py", line 905, in get_registered_maps     dbif.cursor.execute(sql)     sqlite3.OperationalError: near "'2 months'": syntax error      I am using the data specified in the presentation (mean temperatures dataset downloaded from http://eca.knmi.nl/download/ensembles/data/Grid_0.25deg_reg/). Please help. Thank you, Bogdan Rosca


title: Searching closest region from Google Maps coordinates

text: Here is my problem. I'm mapping a region in google maps after that I'm saving the coordinates in PostGIS databse as geographic Polygon type.               ST_GeographyFromText('SRID=4326;POLYGON((..))')      The goal of my script is to find the closest region to a given point selected from google maps. Here is my solution which is not working very well. So what i'm doing is to break the region polygon and to find the closest point to my desire point. But the accuracy is not very precise. Some times when my point is closer to a given region it selects a region which is pretty far from it.               SELECT             ST_Distance (                 ST_GeographyFromText('SRID=4326;POINT(" . $point->longitude . " " . $point->latitude . ")'),                 geom::geography             ) AS distance         FROM             ST_DumpPoints (gis_coords::geometry) --gis_coords is the region column.         ORDER BY             distance ASC         LIMIT 1     ) AS distance      When I try to render my region column with QGis and Open Layers Plugin the result is pretty disturbing. **QGis Settings** I'm setting the region layer CSR to WGS 84 (EPSG:4326) and for Open Layers WGS 84/Pseudo Mercantor (EPSG: 3857) ![enter image description here](http://i.stack.imgur.com/IDkoy.png)


title: Reference control on Editform from tableform - ArcPad VBscripting

text: I'm developing an application to record water chemistry data along a profile from the surface of a lake, every 0.5, 1, 1.5 or 2m, down to the bottom of the lake. On the Edit form, the user will record data associated with the sample location (date, time, lake name, UTM co-ordinates and a sampling interval). The sampling interval will vary depending on a few things. Basically I'm trying to execute a small piece of code to transfer the value of the sampling interval during an onload event for the related table, so that I can use that value to calculate a depth for each consecutive sample along the water profile. I've had trouble working with related tables in the past, so maybe I'm running into another limitation of ArcPad. I've tried this piece of code on the related tables onload event for the form and the page with little luck. Essentially I'm trying to grab the sample depth (txtSampleDepth) entered on the editform. Using an SQL statement, insert that value into the related table and update the AXF timestamp and status. Any suggestions to my code or a different approach would be greatly appreciated!                 Dim pDS       Dim pControls       Dim pSamplDpth       Dim sValList            Set pDS = Map.SelectionLayer.DataSource       Set pControls = ThisEvent.Object.Parent.Controls       Set pSamplDpth = pControls("txtSampleDepth")       SampDpth = pSamplDpth.Value            sValList = SampDpth & "','" & "', '" & getDate() & "','1"            pDS.Execute("INSERT INTO DEEPPOINT_PROFILE (DEPTH, AXF_TIMESTAMP, AXF_STATUS) VALUES ('" & sValList & "')")


title: 3D Barriers, ArcGIS Network Analyst

text: I am trying to implement 3D barriers from a feature class into a network using a model, but I can't seem to get the resulting barrier to appear on the proper three dimensional plane. The feature class that the barrier is generated from a 3D feature, and displays properly in ArcScene. Yet when I run my model, the barrier that appears in the route layer doesn't retain the 3D information. I need to have the 3D information, as I don't want to apply the barrier to all features just to the features on that specific plane. I have the barriers set up via a second "Add Location" tool, as a polygon barrier after the Add Location tool that generates the stops. Is there something that I need to set in the tool itself to retain the Z information from the feature class? I do not have "Use Network Location Fields instead of Geometry" selected. Below a picture of the problem. The green is the feature class, red is the barrier that is created after running the model. I should also add that I do not need the barrier to be extruded, or represent a three dimensional box or object, just that it needs to have the correct elevation. I am using ArcGIS 10.1, with the Network Analyst extension, model creating the 3D routes via Model Builder. Thanks, eric ![Barrier Problem](http://i.stack.imgur.com/tN2Ml.jpg)


title: Why does geoserver layer disappear when baselayer is false?

text:
I am trying to overlay a layer from Geoserver on an OpenStreetMap base layer. If I load the Geoserver layer as a base layer, I can see it - but if I set baselayer to false it disappears. My code is:               <!DOCTYPE HTML>                     <title>OpenLayers Demo</title>

             Shows the basic use of openlayers using a WMS layer
                  I've used firebug to investigate the calls being made to geoserver. In the case where I add the WMS as a base layer I see the following:                      whereas when baselayer is false it's                     The first call returns a PNG showing my data, the second one is empty. A final point is that my data in Geoserver has native SRS EPSG:4258 and declared SRS EPSG:4326 - I have set the "Reproject native to declared" option.


title: How to integrate every single street into a network analysis?

text: This is my first question, but I solved problems before thanks to this great community! It is also my first time working with a network analysis (ArcMap 10) and there are some problems which haven't been solved until now. Maybe you can help! I have bus stops and Streets (with length, slope, time from-to for a pedestrian (2.5 km/h)). I want to figure out from which point a person in a village can reach the bus stop within 5, 10 and 15 min. ( I did an ESRI tutorial about this topic, and actually it shouldn't be so difficult, but the tutorial data was well prepared) At first I don't want to add slope restrictions. I just want to generate polygons around the bus stops with 5,10 and 15 min. The result is quite weird. It seems that some streets have been analysed but others haven't been. (see screen) All roads are connected and they have no further restrictions so i wonder why they are not displayed in analysis. Thanks for your help and ideas! ![enter image description here](http://i.stack.imgur.com/MT643.jpg)


title: How to implement mouse down event when the user clicks on the map?

text: Newbie programming in C#.net. I have a very simple desktop application in in arcmap. I'm trying to get the coordinates when a user clicks a button on a toolbar and when click the map. The toolbar is implemented using the BaseToolBar. I use the base tool to grab the mouse down event ( Not sure whether this is the correct way even) Here is the code which it generates.               using System;     using System.Drawing;     using System.Runtime.InteropServices;     using ESRI.ArcGIS.ADF.BaseClasses;     using ESRI.ArcGIS.ADF.CATIDs;     using ESRI.ArcGIS.Framework;     using ESRI.ArcGIS.ArcMapUI;     using System.Windows.Forms;          namespace GIS_FORMS     {     /// <summary>     /// Summary description for GetXY.     /// </summary>     [Guid("b1e4d6f0-53de-44d3-b679-f29432d1fe35")]     [ClassInterface(ClassInterfaceType.None)]     [ProgId("GIS_FORMS.GetXY")]     public sealed class GetXY : BaseTool     {     #region COM Registration Function(s)     [ComRegisterFunction()]     [ComVisible(false)]     static void RegisterFunction(Type registerType)     {     // Required for ArcGIS Component Category Registrar support     ArcGISCategoryRegistration(registerType);          //     // TODO: Add any COM registration code here     //     }          [ComUnregisterFunction()]     [ComVisible(false)]     static void UnregisterFunction(Type registerType)     {     // Required for ArcGIS Component Category Registrar support     ArcGISCategoryUnregistration(registerType);          //     // TODO: Add any COM unregistration code here     //     }          #region ArcGIS Component Category Registrar generated code     /// <summary>     /// Required method for ArcGIS Component Category registration -     /// Do not modify the contents of this method with the code editor.     /// </summary>     private static void ArcGISCategoryRegistration(Type registerType)     {     string regKey = string.Format("HKEY_CLASSES_ROOT\\CLSID\\{{{0}}}", registerType.GUID);     MxCommands.Register(regKey);          }     /// <summary>     /// Required method for ArcGIS Component Category unregistration -     /// Do not modify the contents of this method with the code editor.     /// </summary>     private static void ArcGISCategoryUnregistration(Type registerType)     {     string regKey = string.Format("HKEY_CLASSES_ROOT\\CLSID\\{{{0}}}", registerType.GUID);     MxCommands.Unregister(regKey);          }          #endregion     #endregion          private IApplication m_application;     public GetXY()     {     //     // TODO: Define values for the public properties     //     base.m_category = "GIS_FORMS"; //localizable text     base.m_caption = "get xy"; //localizable text     base.m_message = "get xy"; //localizable text     base.m_toolTip = "get xy"; //localizable text     base.m_name = "get_xy"; //unique id, non-localizable (e.g. "MyCategory_ArcMapTool")     try     {     //     // TODO: change resource name if necessary     //     string bitmapResourceName = GetType().Name + ".bmp";     base.m_bitmap = new Bitmap(GetType(), bitmapResourceName);     base.m_cursor = new System.Windows.Forms.Cursor(GetType(), GetType().Name + ".cur");     }     catch (Exception ex)     {     System.Diagnostics.Trace.WriteLine(ex.Message, "Invalid Bitmap");     }     }          #region Overriden Class Methods          /// <summary>     /// Occurs when this tool is created     /// </summary>     /// <param name="hook">Instance of the application</param>     public override void OnCreate(object hook)     {     m_application = hook as IApplication;          //Disable if it is not ArcMap     if (hook is IMxApplication)     base.m_enabled = true;     else     base.m_enabled = false;          // TODO: Add other initialization code     }          /// <summary>     /// Occurs when this tool is clicked     /// </summary>     public override void OnClick()     {     // TODO: Add GetXY.OnClick implementation     }          public override void OnMouseDown(int Button, int Shift, int X, int Y)     {     // TODO: Add GetXY.OnMouseDown implementation          MessageBox.Show("X position is " + X.ToString() +"Y position is :" + Y.ToString());     }          public override void OnMouseMove(int Button, int Shift, int X, int Y)     {     // TODO: Add GetXY.OnMouseMove implementation     }          public override void OnMouseUp(int Button, int Shift, int X, int Y)     {     // TODO: Add GetXY.OnMouseUp implementation     }     #endregion     }     }      Now I need to get the coordinates out. How can I do this? where should I call the event mouse down event. The message box does not display at all. Pleas help.


title: Are layer files compatible between arcgis 10 and 9.3?

text: If I save layer files in ArcGIS 10, can they be opened and modified in ArcGIS 9.3? Many thanks for your help.


title: Georeference in ArcGIS to a known scale?

text: If I'm georefencing an engineer's drawing and the scale is 1:250 is there anyway I can input this into ArcGIS for accuracy once I've selected a few points to georeference?


title: Best way to find the polygons crossed by a line

text: I'm trying to find all the polygons crossed by a single line (a GPS track). I'm using the OGR library (from python) for computing this, but it's currently a bit 'brute-force' (and slow). For every point of my track, I call the intersect method with all the polygons. The obvious optimization is to check only with adjacent polygons. But I guess this is a classical problem, with an already known solution (which I can't find...). I would like to avoid using a dedicated database as I'm trying to write a standalone software (spatialite is an option if the DB is the way to go). (FYI, the current source code is available here: https://github.com/dkm/airspace-checker )


title: Changing the default symbology field when adding DWG to ArcMap 9.3.1

text: When I add a dwg to ArcMap 9.3.1, the symbology is set, by default, to something the program calls "CAD unique entity values"(CAD Renderer?) - a combination of "Linetype", "Color" and "LineWt" fields. Is there any way to change the default so that any random DWG I load will be symbolized according to the "Layer" field (for example) and not those three field mentioned above?


title: ArcMap vertex problem for Z-value

text: I have an Oracle Spatial dataset with z-coords. I can edit the dataset in ArcMap however the z-values are always rounded to the nearest integer. Why does ArcMap round the z-values? Regards, Pieter


title: How to embed map into access database?

text: I am developing a database for my organisation that could be spatially enabled. I have to use MS Access 2007 on Windows and am therefore looking for any possibilities to encorporate my access DB into a mapping tool to embed inside my access application. Has anyone successfully created an MS Access project which embeds a mapping window? I would ideally like to simply geocode my datasets - digitize a point after inputting my data into the db. If anyone can shed any light on whether this is possible I would be grateful.


title: Changing default precision and scale of number fields in ArcSDE geodatabase table?

text: After I change the precision and scale for my short, long, float, and double field types when creating a table in ArcMap SDE geodatabase, ArcMap will reset my precision and scale to default values. For instance, I wanted my short integer to have a precision of 3 in one table and 2 in the other so I could test some python code. ArcMap changed them both to precision of 5 automatically and the python tool isn't given a proper test. Please note that this is a different but related issue that was addressed here: Arcpy - How do I specify the precision and scale of a new field? I know how to specify the precision, but ArcMap isn't keeping my settings. If this makes sense stop reading here. If you want more detail, please continue reading: I am creating an enterprise geodatabase in ArcMap 10.2.1 for the purpose of testing two tables: a source table and a destination table. I am building a tool that will describe the tables, put the field names as keys in a dictionary and the field objects as values, populate a list of all matching fields, then a series of if statements remove the fields that may have the same name, but incompatible data types, precisions, scales, and field lengths. Incompatible is different than non-matching. I created if statements that allow a short integer in a source table to be compatible with a long integer(or a float, or a double) in a destination table. But this cannot be reversed. Only a short can go into a short. But if they have different precisions, a smaller or equal precision can go into a larger precision, but a larger precision cannot go into a smaller one. This is why I need ArcMap to keep my varying precision values across different fields. I need to be able to test every possibility.               #           Pass SmallInteger types if going into a SmallInteger.                      if source_type == 'SmallInteger' and destination_type == 'SmallInteger':                     pass          #           Filter out incompatible precision.                     if source_precision <= destination_precision:                         pass                     if source_precision > destination_precision:                         matched_column_list.remove(name_1.lower())                         continue      That is my issue. If anyone can get my enterprise geodatabase to behave, please let me know! I'm using ArcGIS 10.2.1.


title: Find minimum elevation along a polyline and return xy or point?

text: Using ArcGIS 10.1 and Python IDLE 2.7.2... I'm trying to find the minimum value of a polyline with a raster layer and turn that LOCATION into a point. I put everything together in model builder, but the raster calculator doesn't like 32-bit floating points, which is what my DEM is. I tried to get around the raster calculator with the Con tool in python, but I believe an attribute table is required for this (which I don't have due to floating point). I could convert the DEM to integer and build one, but this is already a lengthy process that involves many scripts and tools and wanted to see if anyone else out there had any thoughts? I'm open to any anything, really, conversion of polyline to polygon, etc...attaching a screenshot of modelbuilder process that works, but I need a script as I am doing this 1000+ times.![enter image description here](http://i.stack.imgur.com/BF1qm.png) I can turn the polyline into a 3D feature and find the minimum value, but I don't know how to find WHERE that minimum value is. If I could locate it with XY coordinates, I can then generate the point. Thoughts? If I'm wrong about the Con tool and my syntax is wrong, that would be awesome, so please feel free to suggest something. It finds where the min is along the polyline that has been converted to raster and where that min is within the DEM, thus returning a min pixel: tempRas1 = Con(Raster(LG_DEM_0_tif == LG_0_Zonalmin_img), LG_DEM_0_tif)


title: How do I set the xydomain extent and merge point feature classes in a geodatabase

text: I am trying to select point feature classes in a feature dataset and merge them while setting the XYDomain using Python. If I do the merge manually and change the environment settings to use a feature class that encompasses the point data, all works well. However, I have many more geodatabases within this project that I'd like to do the exact same thing to. Using the code below, I get error 000594. I have also tried changing the xyDomain to point to a shapefile but that did not work either.               import arcpy     from arcpy import env     arcpy.env.XYDomain ="-5120900 -9998100 900714804574 900709927374"     env.workspace = "Q:\Project_Data\WSWS 2012 #660\WSWS_4.gdb\ALABAMA_GPS"     fcList = arcpy.ListFeatureClasses('','POINT','')     arcpy.Merge_management(fcList, "MergePNT6")      Code that worked:               import arcpy     from arcpy import env     arcpy.env.XYDomain ="-5120900 -9998100 900714804574 900709927374"     arcpy.env.outputMFlag = "Disabled"     env.workspace = "Q:\Project_Data\WSWS 2012 #660\WSWS_4.gdb\ALABAMA_GPS"     fcList = arcpy.ListFeatureClasses('','POINT','')     arcpy.Merge_management(fcList, "MergeOutputName")