title: Is it possible to open a shapefile using a ArcGIS Runtime for Windows Mobile application

text: I need to build a tablet application using ArcGIS Runtime for Windows Mobile. Is there a way to open a shapefile (or a pgdb) without having to convert it to a .db file first?


title: Sencha / Extjs 4 with GeoExt 2 - How to use

text: I'm working on Sencha Architect 2 - extjs 4 i found out that GeoExt doesn't work with extjs 4. but i really need to use it so i searched and end up here https://github.com/geoext/geoext2 I'm trying to something like examples so i just add :               Ext.Loader.setConfig({         enabled: true,         disableCaching: false,         paths: {             GeoExt: "../geoext-2/src/GeoExt",             // for dev use             Ext: "http://cdn.sencha.io/ext-4.1.0-gpl/src"             // for build purpose             //Ext: "extjs-4.1.0/src"         }     });      to my app.js file. but i keep getting "GeoExt is not defined" error message. anyone could help me with that ?


title: How to split lines at vertices in ArcGIS for Desktop with Basic level license?

text: Is there a way to work around splitting lines at vertices in ArcGIS for Desktop without an Advanced level license? I don't have the advanced license and I have a model that uses this tool. If need be, I can also use QGIS and then modify my model but I prefer an ArcGIS method.


title: Why outside border of a service area is always not reasonable?

text: In ARCGIS 10.1 network analysis, why results of service area are different if number of break points is changed? For example, the grid is 300m*300m, I try to produce service area around a node. settings are like this... ![enter image description here](http://i.stack.imgur.com/oiqoV.jpg) ![enter image description here](http://i.stack.imgur.com/BpK3y.jpg) Then I got this result ![enter image description here](http://i.stack.imgur.com/JBNlp.jpg) then I **only change breaks to 300,600** , result: ![enter image description here](http://i.stack.imgur.com/Sdlti.jpg) Why the outside border is always in strange shape? Why the **service areas of 600m** are different in these 2 cases? I try to change polygon generation to "detailed", it still generate this kind of results. Really confused, and I have problems understanding any service area result because of this. Highly appreciated to any ideas!


title: How do i install QGIS with FileGDB?

text: Please help me i am very new at this and the last time I asked i was sent to a link about an older version of QGIS. I need to know how to install QGIS 2.0.1 for Windows 8 64 bit so it can recognize and open ESRI .gdb files. I know its not hard for you all, but pretend i know nothing about anything, what are my steps?? What do i download? What settings do i change?? Please do not just send me here: How to get ".GDB" (Esri File Geodatabase) support in Quantum GIS (osgeo4w - qgis) As the information at this link was posted two years ago, I am looking for an update!!


title: Date fails in String -> Date in CartoDB

text: I'm trying to convert a column in CartoDB from String to Date. My dates are the following format:               2014-01-02T00:02:00     2014-01-02T08:05:00     ...      I've tried to append +00:00, so they looks like:               2014-01-02T00:02:00+00:00     2014-01-02T08:05:00+00:00     ...      matching CartoDB's own dates. However, when converting from String to Date the time is lost and they become:               2014-01-02T00:00:00+00:00     2014-01-02T00:00:00+00:00     ...      What to do?


title: Creating Network Dataset programatically

text: I'm trying to code to create a network dataset using a single line feature class with minimal attributes I'm looking for a nice tutorial which can guide me through the process.I have already looked at the examples provided by esri (http://help.arcgis.com/en/sdk/10.0/arcobjects_net/conceptualhelp/index.html#//000100000466000000) but they are really complex and specific to the question, its really hard to get the concept from it.


title: What are the appropriate parameters to convert .asc file to geotiff using gdal_translate?

text: My starting file of interest is an ASC file with the following gdalinfo               Driver: AAIGrid/Arc/Info ASCII Grid      Files: sd_1_3_110805.asc      Size is 11880, 9720      Coordinate System is `'      Origin = (-117.900000000000010,33.399999999999991)      Pixel Size = (0.000092592592593,-0.000092592592593)      Corner Coordinates:      Upper Left (-117.9000000, 33.4000000)      Lower Left (-117.9000000, 32.5000000)      Upper Right (-116.8000000, 33.4000000)      Lower Right (-116.8000000, 32.5000000)      Center (-117.3500000, 32.9500000)      Band 1 Block=11880x1 Type=Float32, ColorInterp=Undefined      NoData Value=-9999       To translate to a tiff I have tried:               gdal_translate sd_1_3_110805.asc lj.tif      Subsequently I would like to change the projection to UTM flat to ultimately create a 3d model of flat vs geocentric terrain. Unless I have specified using a geocentric model to osgdem (VirtualPlanetBuilder) to build the 3d terrain model, my model always comes out mirror image (north is south etc...). So I desire to use a flat terrain model... what parameters should I be using to get to a geotiff with the correct projection to start with?


title: Display XY after Button Click

text:
I am creating a map that will generate an XY coordinate based off the location that a user clicks. I only want this location to be generated after the user clicks a button to enable the draw toolbar and then clicks the map. The coordinate will be passed to another web page and ultimately stored in a database so I need to make sure that a point can't accidentally be created. Code below as well as a link to what I currently have running. I may be wrong but I think it should have something to with a dojo.connect but I may be wrong. But you can see the commented out options that I have tried.               <!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01//EN" "http://www.w3.org/TR/html4/strict.dtd">                     <meta http-equiv="Content-Type" content="text/html; charset=utf-8"/>         <meta http-equiv="X-UA-Compatible" content="IE=7" />         <meta name="viewport" content="width=device-width,user-scalable=no"/>                  <meta name="viewport" content="initial-scale=1, maximum-scale=1,user-scalable=no"/>         <title>Maps Toolbar</title>         <link rel="stylesheet" type="text/css" href="http://serverapi.arcgisonline.com/jsapi/arcgis/2.2/js/dojo/dijit/themes/claro/claro.css">

           <button dojoType="dijit.form.Button" onClick="projectToLatLong(evt); toolbar.activate(esri.toolbars.Draw.POINT);">Show XY</button>

                       http://geoville.org/viewers/XY/ Thanks for your help. -Mike


title: Convert polygon data into extreme lon/lat coordinates using QGIS

text: I'm very new to QGIS. I have a shape file of census blocks and want to use the data to plot points within each boundary using Google Maps api. I figure the best way to do this is to find the extreme north, south, east and west lat/long for each polygon (so later I can plot points within those boundaries). I'm wondering:   * Is there a way to convert the poloygons into lat/lon coordinates and then find the 4 extreme boundary points?   * Is there a better way to go about plotting points in boundaries that I'm missing?


title: Three-dimensional interpolation of temperature profiles using gstat in R?

text: I'm working on a problem in which I'd like to interpolate a set of oceanographic temperature profiles over a three-dimensional matrix of unsampled x, y, & z locations. The data to interpolate is housed in a SpatialPointsDataFrame (WOD) in which I've stored temperature data (WOD$TEMP) in association with coordinate information in the form of latitude (WOD$LAT,in degrees),longitude (WOD$LONG, in degrees), and depth (WOD$DEPTH, in meters). The data consists of multiple temperature measurements collected on different depth levels at each station (i.e. multiple TEMP and DEPTH values for each LAT and LONG pair). Here's an sample of the data which is drawn from the World Ocean Database (I can send along a more complete dataset as a .csv if you're interested)                  head(WOD[WOD$YEAR%in%2013 & WOD$TYPE%in%"PFL",c("LAT","LONG","DEPTH","TEMP")],100)[seq(1,100,2),]                  LAT    LONG   DEPTH   TEMP       4004511 27.644 -78.732    3.68 23.151       4004711 27.644 -78.732   14.20 23.133       4004911 27.644 -78.732   22.75 23.128       4005111 27.644 -78.732   33.77 23.116       4005311 27.644 -78.732   42.41 23.119       4005511 27.644 -78.732   53.04 23.115       4005711 27.644 -78.732   63.96 23.086       4005911 27.644 -78.732   72.50 23.045       4006111 27.644 -78.732   83.82 23.043       4006311 27.644 -78.732   92.36 23.022       4006511 27.644 -78.732  103.29 22.977       4006711 27.644 -78.732  116.59 22.967       4006911 27.644 -78.732  137.54 22.618       4007111 27.644 -78.732  176.15 20.855       4007311 27.644 -78.732  222.69 19.743       4007511 27.644 -78.732  272.30 18.979       4007711 27.644 -78.732  320.31 18.334       4007911 27.644 -78.732  371.48 17.855       4008111 27.644 -78.732  445.34 16.469       4008311 27.644 -78.732  543.75 13.591       4008511 27.644 -78.732  663.89  9.110       4011811 26.952 -77.340   10.33 23.660       4012011 26.952 -77.340   20.07 23.699       4012211 26.952 -77.340   30.00 23.693       4012411 26.952 -77.340   39.83 23.683       4012611 26.952 -77.340   50.16 23.664       4012811 26.952 -77.340   59.50 23.641       4013011 26.952 -77.340   69.53 23.601       4013211 26.952 -77.340   80.05 23.580       4013411 26.952 -77.340   89.79 23.552       401369  26.952 -77.340   99.52 23.625       401382  26.952 -77.340  109.45 23.349       4014011 26.952 -77.340  129.60 22.105       4014211 26.952 -77.340  158.99 20.615       4014411 26.952 -77.340  198.49 19.640       4014611 26.952 -77.340  249.10 18.740       4014811 26.952 -77.340  298.21 18.209       4015011 26.952 -77.340  347.80 17.910       4015211 26.952 -77.340  397.58 17.292       4015411 26.952 -77.340  496.51 15.465       4015611 26.952 -77.340  594.80 13.275       4015811 26.952 -77.340  744.02  9.963       4016011 26.952 -77.340  892.25  6.338       4016211 26.952 -77.340 1090.28  5.297       ...               WOD=WOD[WOD$YEAR%in%2013 & WOD$TYPE%in%"PFL" & WOD$REGION%in%c("TOTO","SEPC","NWPN","NWPS"),]             coordinates(WOD) = ~LONG+LAT+DEPTH                  #Based on the demo(gstat3D) example I have fitted the following variogram models:             library(sp);library(gstat);library(raster);library(rgdal)             WOD_vg=variogram(TEMP~1,WOD)        WOD_vgm=fit.variogram(WOD_vg,mod=vgm(psill=2,"Exp",range=1))        WOD_vgm1=fit.variogram(WOD_vg,mod=vgm(psill=1.8,"Sph",range=1.0,nugget=0.6))        WOD_vgm2=fit.variogram(WOD_vg,mod=vgm(psill=2,"Sph",range=1.5,nugget=0.6))        WOD_vgm3=fit.variogram(WOD_vg,mod=vgm(psill=2,"Mat",range=1.5,nugget=0.6))             plot(WOD_vg$dist,WOD_vg$gamma)        lines(variogramLine(WOD_vgm,1.5)$dist,variogramLine(WOD_vgm,1.5)$gamma)        lines(variogramLine(WOD_vgm1,1.5)$dist,variogramLine(WOD_vgm1,1.5)$gamma)        lines(variogramLine(WOD_vgm2,1.5)$dist,variogramLine(WOD_vgm2,1.5)$gamma)        lines(variogramLine(WOD_vgm3,1.5)$dist,variogramLine(WOD_vgm3,1.5)$gamma)             # and attempted to krige the data and variogram model over three dimensional SpatialPixels defined in terms of ~LONG+LAT+DEPTH             grid3D <- expand.grid(LONG=seq(from = -79.00, to = -76.50, by=0.05),                              LAT=seq(from = 23.25, to = 26.75, by=0.05),                              DEPTH=seq(from = 0, to = 1200, by=50))        gridded(grid3D) = ~LONG+LAT+DEPTH             res3D <- krige(formula = TEMP ~ 1, WOD, grid3D, model = WOD_vgm2)      #This produces the following error:               > res3D <- krige(formula = TEMP ~ 1, WOD, grid3D, model = WOD_vgm2)     [using ordinary kriging]          "chfactor.c", line 131: singular matrix in function LDLfactor()     Error in predict.gstat(g, newdata = newdata, block = block, nsim = nsim,  :        LDLfactor      This error I gather from reading lots of earlier posts on various forums probably relates to having multiple observations (at different depth levels) with same LAT and LONG coordinates (i.e. horizontal distance between these observations = 0). I'm relatively new to the world of geostatistics so I have a number of questions: 1) Is 3D kriging using gstat an appropriate method of interpolation in this context? If not is there an alternative approach that I should pursue? 2) How should I deal with having observations on multiple depth levels at the same LAT and LONG coordinates? 3) Do latitude, longitude and depth need to be in the same units (e.g. meters) for accurate estimation of distances in variogram calculation and kriging interpolation? 4)How should I set initial parameters for the anisotropy term given that temperature varies by almost an order of magnitude (4-30 C) over 1000m in the vertical direction but <=1 degree C over >100 km in the horizontal dimensions (more with latitude than with longitude)?


title: Get IFeatureLayer from ILayer using ArcGIS Engine?

text: I try to select object from layer by name. I think i found sollution here. First problem i got is that i use GroupLayer. So i made research how to get subLayer from GroupLayer. Sollution is here. Now i have a problem how to get IFeatureLayer if i already have ILayer from where i need to select object.               var mapControl = form.axMapControl2.ActiveView.FocusMap;          var basemapCompositeLayer = mapControl.Layer[0] as ICompositeLayer;     var basemapSubLayer = basemapCompositeLayer.Layer[0];     var featureLayer1 = basemapSubLayer as IFeatureLayer; // always NULL      Or this is because that Layer isn't Featured by default?


title: gdal's polygonize function (via qgis) is not turning monochrome rasters into readable shapefiles regardless of inversion of color scale

text: I've tried polygonizing this monochrome bitmap as this: ![enter image description here](http://i.stack.imgur.com/mktPO.png) and this: ![enter image description here](http://i.stack.imgur.com/nthfu.png) and either way it ends up looking like this when (zoomed in) after gdal's polygonize operation: ![enter image description here](http://i.stack.imgur.com/hIdty.png) Is there any way I can tell it which binary component I want to be the shape and which I want to be the void?


title: Per-pixel (statistical) calculations on a raster stack using GDAL

text: In R, it's relatively trivial to perform per pixel calculations based on a raster stack (e.g., get std.dev for each pixel on a 12 layer GeoTIFF). Unfortunately, the speed is less than desirable when working with large rasters. Recently, I've been trying to use GDAL for this purpose, but I can't seem to get the expression down. I've tried the following, but it doesn't seem to work as advertised (or how I think it should work):               gdal_calc.py -A tmean_monthly_1980.tif --calc="std(A)" --outfile=tstdev_1980.tif --allBands=A --calc="std(A)" --NoDataValue=-9999      In the above example `tmean_monthly_1980.tif` is a raster composed of 12 layers (Jan-Dec) of mean monthly temperatures in 1980 (stacked using `gdal_merge`). Although the GDAL command produces an output, the results are definitely not the standard deviation per pixel. Another option might be to use single band files, but `gdal_calc.py` appears limited to "only" 26 inputs (A-Z; i.e., "`-A file -B file -C file` etc."), and the resulting command line/script would be fairly cumbersome when working with large, multi-layered rasters (e.g., time series). **Is there a specific syntax for performing per-pixel raster calculations on a multi-band image using GDAL? If only individual files/bands can be specified using`gdal_calc.py`, please provide a scripted example using Bash, Python, etc.** Specifically, my goal is to calculate the **standard deviation** and **cv** for **each pixel** in a multi-layer raster stack, or using individual files that make up the same raster stack. **Note** : I'm holding out for a GDAL specific answer, but other FOSS solutions are welcomed and will be upvoted, especially if they are relatively compact and adaptable to large number of inputs.


title: Delete duplicate spatial point/points in mssql from a perticular location

text: I am working with a Spatial data which has a street intersection points.... I get the intersection points from spatial query on street/roads datain MMSQL . but in my case some streets/roads are multi lane.. which gives me duplicate intersection points for same streets/roads for each lane about 3-5 meters away from each other.. now i have the point data with Duplicate intersection points with a geometry differing in about 3-5 meters.. Now i want to delete the duplicate intersection points and keep only one intersection point.. i tried taking buffer of a point n den deleting the points which intersects with this buffer.. the the query is very slow and takes ages to complete.. i tried using spatial indexes but it hardly helps... is there a algorithm or any method in QGIS ARCGIS OR POSTGIS Database which i can use to delete theses duplicate points.. or any other method?????


title: WFS, ArcGIS Server 10.2 and WFS maxFeatures?

text: I'm experiencing intermittent 'out of server memory' errors on our public ArcGIS Server 10.2 server. From what I can glean from the OGC documentation the server side maxfeatures parameter was only available in the WFS 1.1.0 spec but not 1.0.0 I suspect that our server is temporarily being overwhelmed by WFS 1.0.0 requests. IE some feature types have in excess of 90,000 records We need WFS capability on our services for integration with another system. The other system is being coded to perform WFS 1.1.0 requests so the maxfeatures will apply. Is there any way I can 'turn off' WFS 1.0.0 support in ArcGIS Server 10.2 without turning 1.1.0 off? The only thing ESRI support and I came up with was to create a custom getcapabilities xml version 1.0.0 file and delete all feature types in that minus the one with least records... Regards, Richard


title: How to add new menu item to ArcCatalog using ArcObjects?

text: I'm trying to a add a new menu drop down to the ArcCatalog menu bar. I tried to use the built in ArcGIS base menu template in Visual Studios 2010 but no luck. I also tried to manually add the dll create with the base menu template in the custom dialog box but my menu does or menu items do not appear in the list. I think my best approach is to use the IExtension but all the example I see online is using ArcMap. The examples retrieve the main menu by getting it from a new or open document using the "IApplication...CommandBar.Find" but I don't think there's a "document" in ArcCatalog? Any ideas would be grateful. Thank you. I tried the following with no luck. I assume the document return from the IApplication is a GxDocument. I might be wrong.                   public void Startup(ref object initializationData)         {             ESRI.ArcGIS.RuntimeManager.Bind(ESRI.ArcGIS.ProductCode.Desktop);             ICommandBar topMenu = null;             UID topMenuUid = new UIDClass();             object index = null;             ICommandBars commandBars = null;             IDocument document = null;                  if (initializationData is IGxApplication)             {                 m_application = (IApplication)initializationData;                 topMenuUid.Value = ARCCATALOG_MAIN_MENU_UID;                 document = m_application.Document;                 ICommandBars commandBars = document.CommandBars;              }         }


title: Find nearest point (in file A) for every point (in file B)

text: I have 2 csv files both with 3 columns: variable, lat, lon. Let's call them file A and file B. File A has 100 values for variable (in this case city), File B has 100,000 values for variable (station numbers in this case). The goal is to assign a city (from file A) to every station number (from file B). So for each city, I will have a station number which is the closest station number to that city (by euclidian distance using latitude longitude). Does anyone know how I can do this? I've looked at using Nearest but that doesn't seem like what I'm looking for. I can use ArcPy also if that helps, but I'm not too experienced in ArcGIS.


title: FME - How to write all attributes when multiple source files have the same name?

text: I am using FME to write all the MapInfo TAB files in a directory and it's subdirectories to a File Geodatabase. I'm using dynamic readers and dynamic writers as I don't know the geometry or schema of the input files beforehand. So far so good. Now, suppose that I have a file called **`"Pipeline.TAB"`** at the root of the directory and I have another file called **`"Pipeline.TAB"`** in a subdirectory. No problem. FME can merge the two files together and I get the geometry from both files in my new feature class called **`"Pipline_line"`**. My issue is, suppose those two original files have different attributes. Say, the file in the root folder has a field called "Name" and the second file has a field called **`"Pipe_Name"`**. My resulting feature class gets the **`"Name"`** field but not the **`"Pipe_Name"`** field. The features that came from the file in the subdirectory have a _NULL_ value in the **`"Name"`** field. Instead, what I would like is to have both fields come through. I'm happy if the **`"Name"`** field is _NULL_ as long as the **`"Pipe_Name"`** field is there and has the correct values. It appears that the schema is determined by the first feature that is being read, whereas, I need it to use the schema from all of the features with the same name. Is this possible? **EDIT** I have posted this question on FMEPedia as well. Here is some additional information which may help: The example I gave was simple in that it only had two files called Pipeline.TAB. In reality, I'm looking at dozens of folders with thousands of files. there might be 20 files called "Pipeline" and another twenty called "Valves", etc. I won't know the schema of all of them at the start. What I would like the final attribute table of my destination feature class to look like would be something like this:               Name       | Pipe_Name     Pipe1      | NULL     Pipe2      | NULL     Pipe3      | NULL     NULL       | Pipe 4     NULL       | Pipe 5      In the final file, there may actually be many fields. But I won't know all of them at the start:               Name      | Pipe_Name  | Diameter | Material | Type     Pipe1     | NULL       | 500      | Steel    | Gas     Pipe2     | NULL       | 1200     | Steel    | Gas     Pipe3     | NULL       | 800      | PVC      | Water     NULL      | Pipe 4     | NULL     |NULL      | Water     NULL      | Pipe 5     | NULL     |NULL      | Water      In response to @BradNesom's answer, I have looked at the SchemaMapper, but I have way too many files to look at for that to be feasible. I also have looked at the other link you mentioned and I tried adding a reader as a resource but that didn't seem to work either.


title: Summarize values from a raster brick by latitude bands in R

text: I would like to sumarize (sum and/or mean) values from layers on a raster brick based on latitude bands (been using raster package). Been looking at the package documentation and didn't find a way to do it. Can anyone give any advice on how to do it? I guess one option would be to use extract by spatial polygons? or is there a way to do it with another raster where the values would be classes of the latitudes of cells?


title: Torque for chloropleth and bubble maps?

text: Torque looks like a brilliant tool. However, can I use the same tool for chloropleth maps for visualizing monthly/annual changes of GDP, health expenditure, population growth etc. over the past 10 years? Will it work with bubble charts (changing the size of bubbles rather than just mapping singular events?).


title: How to replicate trunc function on QGIS?

text: I've search a bit before asking this question, but I haven't found the answer. Is there a Trunc function on QGis? I've tried using floor, but it doesn't work (like I want) if the number is negative, for instance: If `$x=-8.639` Floor will give me `floor($x)=-9` And what I'm looking to achieve is -8 (it's the trunc function). (I'm passing angles to DMS notation) I think this may be resolved with a conditional statement for negative values and positive values I also found this expression that makes the job `regexp_replace($x, '\\\\..*', '')`... but I guess there should be another (cleaner) way to do it. I'm using QGis 2.0 on a 64bit Win7 machine.


title: How to compute the routing through N addresses/coordinates?

text: I have a set of _N_ coordinates (addresses). My goal is to compute the shortest driving route through the set of coordinates, where the algorithm (not the user) decides which of the _N_ coordinates should be the start point and which of the _N - 1_ coordinates (excluding the start point) should be the end point. **Can this be done and which algorithms are used for this kind of routing?** By assuming that the routing can be done using the arial distance between the coordinates instead of roads, **how could I then solve my problem?** _At this stage of my questioning, I am trying to determine how to solve the problem in theory but the tools at hand are a PostGIS database and anything I can implement my self in C# .NET. If any free libraries can help me, then they can be considered too :)_


title: Using field variable to calculate values in ModelBuilder?

text: I am using ModelBuilder and trying to get a user defined field from a user defined polygon. From this, I want to later populate a different field with calculated values in a model. Is there a way to get a list of fields from the user defined polygon and use that value to populate the other field? How would I go about using the user defined field in a VB expression within the calculate value tool? `newfield = [%field%]`?


title: Selection subquery using another table in database ArcGIS

text: I need to 'extract' a list of variables based on the IDs from another table. Using a subquery seems to be the best method. I have a large list of variables and have outside calculated a second list. This second list does not have coordinate data attached to it, so I need to select the same values in the original list to pull the data out. The shapefile and table are both in the catalog. I try to do a simple select by attribute. It gives me the               SELECT * FROM US_Maize_All WHERE: "pointid" = (SELECT "pointid" FROM Maize_Table)      but this only gives me a SQL expression error. The "pointid" parameters are selected from the option box. As far as I can tell this is the exact format as any other info I can find online. Has to be something I'm overlooking why this will not work.


title: Is there a GeoDjango tutorial for GeoRSSFeeds?

text: I'm going through GeoDjango basics right now. Ideally, I'd like to publish a GeoRSS from a set of data I have that is already in Django's database. I know how to get a regular Feed working, but when I try GeoRSSFeed, I'm running into problems. One of my main problems is:               Exception Value: __init__() takes at least 4 arguments (1 given)      In reference to my urls.py file when I try to make my Feed to GeoRSSFeed. Is there no tutorial for getting GeoRSSFeeds to work? I've seen this but it really isn't helping me.


title: Ontario Buildings Data

text: I am looking for a buildings layer for Ontario, particularly the Chatham-Kent area to properly assess setback for wind turbine development in the region. Currently, I am utilizing the buildings data provided on Geogratis; however, I would like to know it's accuracy and if other buildings data is available free or for purchase.


title: digitize esri 9.2 or 9.3

text: Is it a way to digitize indicating a curve (colored on a white or black or different colour background) from point to point ? I think i've been told that, but i don't KNOW to do it ! Please consider a quite low GIS knowleadge....(sorry)


title: Computing an averaged latitude and longitude coordinates

text: How can I compute the average between several latitude and longitude spots ? Should I just compute the Arithmetic mean for both lat and lng ? thanks


title: Error in using/installing v.digit: GRASS GIS

text: When I wanted to edit a shapefile in GRASS I got this message: > Unable to start wxGUI vector digitizer. Do you want to start TCL/TK > digitizer (v.digit) instead? > > Details: No module named grass6_wxvdigit Therefore, as suggested in this link: http://osgeo- org.1803224.n2.nabble.com/wxpython-gui-and-v-digit-td3019439.html and https://svn.osgeo.org/grass/grass/trunk/gui/wxpython/README - 6. VECTOR DIGITIZER I have been trying to compile/install vdigit component. I could install GRASS with C++ (--with-cxx), Python (--with-python) and wxWidgets (--with-wxwidgets) support. But in the next step when I change to the directory "cd gui/wxpython/vdigit" and run make I get the following error:               ..............     cats.cpp:456: error: ‘Vect_rewrite_line’ was not declared in this scope          cats.cpp:461: error: ‘G_debug’ was not declared in this scope     cats.cpp:467: error: ‘Vect_destroy_line_struct’ was not declared in this scope     cats.cpp:468: error: ‘Vect_destroy_cats_struct’ was not declared in this scope          error: command 'gcc' failed with exit status 1     make: *** [OBJ.x86_64-unknown-linux-gnu/_grass6_wxvdigit.so] Error 1      I use Ubuntu 10.04 LTS with GRASS 6.4.1svn. I tried installing `python-dev`, `python2.6-dev` packages as suggested in many forums but the error continued to appear.


title: GIS projects with a raspberry pi

text: Does anyone know of a project using a raspberry pi and gis that would work for a research project for a computer science major? I am looking to do a research project involving a raspberry pi and incorporating gis into it. I am also looking to use python but all ideas are welcome. Very broad but looking for any ideas. Possibly a mapping application.


title: Dialog statement

text: I am trying to write a code which has dialog including such as toggle button, edittext, pushbutton, popupmenu, etc. I have no problem with all, but when the dialog screen which is named "Find Name" is seen on Map Window, I could not do something without closing dialog screen. **I want to work with "Find Name" dialog screen and also use all the tools in MapInfo such as "Radius Select", "Grabber", etc. without closing dialog screen.** How to do either using MapInfo tools and appearing "Find Name" dialog screen on the Map Window at the same time. Thanks in advance.


title: Creating a geoserver layer with curl (but with a bounding box)

text: I'm attemping to load a layer into Geoserver from a PostGIS database using cURL. It works, but I also need to set the bounding box to -180,-90,180,90. Does anyone know how to add the bounding box to the below command? curl -v -u xxx:xxx -X POST -H "Content-type: text/xml" -d `"<featureType><name>Airports</name></featureType>"` http://xxx.xx.com:8080/geoserver/rest/workspaces/PostGIS/datastores/PostGIS/featuretypes Thanks for your help!


title: ArcGIS ArcPy: Iterate over shapefiles and report field "Desc" Length

text: I've just started learning to script using python, i.e. no programming experience. I need to iterate over a folder that contains multiple shapefiles that all start with nlc... I need to iterate over each shapefile and then report the length of the field "Desc". The reason that I need to do this is that the field lengths are not the same and I've generated a model in modelbuilder that iterates over the folder, generates a feature layer (in_memory) filtering the output layer by using a sql expression to return a specified set of records and then append the output records into a single feature class. I only realised half way through running my model that the values in certain shapefiles are larger than the initial field length that I set for the new feature class. An example to get me started, would really be appreciated, thanks in advance. Peter The following is the initial python script that I wrote, not realising that I didn't need to loop through the fields list:                   import arcpy     ... arcpy.env.workspace = r'E:\Projects\H109009\data_received\DOA\NLC2000\tiles'     ... fcs = arcpy.ListFeatureClasses("nlc*","polygon")     ... for fc in fcs:     ...     fields = arcpy.ListFields(fc,"Desc*")     ...     for field in fields:     ...         print("{0} is a type of {1} with a length of {2}"     ...           .format(field.name, field.type, field.length))


title: Arcpy zoom to selected features not working

text: I'm writing a python script to automatically take some snapshots of my ArcMap in 10.2, but I can't get the zoom to selected features tool to work properly. I have a layer called carereport, which I then make a feature layer of by providing a gid value. If I right-click on this new layer and click 'Zoom to layer', the tool works fine. But when I run the following arcpy script, it stays very zoomed out. It's as if some other features are selected, but I can't see where.               mxd = arcpy.mapping.MapDocument(r"CURRENT")     df = arcpy.mapping.ListDataFrames(mxd, "Layers")[0]          arcpy.MakeFeatureLayer_management("carereport", "Siteslyr", """gid = 2739""")          addLayer = arcpy.mapping.Layer(r"Siteslyr")     arcpy.mapping.AddLayer(df, addLayer, "AUTO_ARRANGE")          arcpy.RefreshActiveView()     arcpy.RefreshTOC()          lyr = arcpy.mapping.ListLayers(mxd, "Siteslyr", df)[0]     df.zoomToSelectedFeatures()          df.scale *= 1.1     arcpy.RefreshActiveView()     arcpy.mapping.ExportToJPEG(mxd, r"C:\caredev.jpg", df,                                df_export_width=1600,                                df_export_height=1200,                                world_file=True)     del mxd, df


title: NameError: name 'QImage' is not defined

text: Traceback (most recent call last): File "", line 1, in File "F:/QGIS/bin/test.py", line 6, in img = QImage(QSize(800,600), QImage.Format_ARGB32_Premultiplied) NameError: name 'QImage' is not defined I encountered this error when I was trying for Map Rendering and Printing on python console given in QGIS. Plz help


title: How do I set interactivity on a sublayer that is created in CartoDB?

text: I am trying to set the interactivity on a sublayer that I created on CartoDB web interface so that my events and tooltips receive the correct fields as arguments of the callback functions. Here is how I have tried to do this:                   layers[1].setInteraction(true)         poverty_layer = layers[1].getSubLayer(0)         poverty_layer.set("interactivity", ["cartodb_id, namelsad10, hh_median"])              poverty_layer.on "featureOver", (e,ll,pos,data)->           console.log data      I assumed that by doing this I would see three data fields –– cartodb_id, namelsad10, and hh_median –– but I only see cartodb_id. I already managed to set interactivity successfully on sublayers that are created on the clientside. This is the one use case where it hasn't been successful. Appreciate any help!


title: Run PostGIS command line functions in pgAdmin?

text: Is it possible to utilize the PostGIS command line functions, e.g. shp2pgsql, from within pgAdmin III? I am using pgAdmin to connect to a remote server.


title: Free Geoserver SLD for 2009 Tiger Line Data (preferably google maps color scheme)

text: I need to render a base map using the Open layers. I have downloaded 2009 Tiger line Data. I need to get up and running quickly. I think I can save some time if i use ready made SLDs, instead of designing them my self. I am looking for some decent looking SLDs especially for the road network and its labels. In fact i prefer if the SLDs are google map style color scheme! :) Please point me to links where i can download these SLDs. Thanks Shaunak


title: How can I tell what licensing level my ESRI GIS server is from client tools?

text: From client tools, like ArcCatalog, is there a way to know if the GIS server I am connecting to is running under a workgroup or enterprise license? Also, is it possible to tell if the edition is basic, standard, or advanced? I ask because we are trying to migrate a 9.3.1 database to 10 (migrate but not upgrade yet) using ArcCatalog. When we add a new database server in ArcCatalog the only options are to add a personal or workgroup database server. It seems like it is only looking for SQL Express instances, but it should be also be able to see full SQL Server (Enterprise) servers. We have ESRI enterprise licensing, but I want to know if it is possible to discover that from the client tools. Help, About ArcCatalog shows ESRI ArcCatalog 10, ArcGIS Desktop 10 SP2, license type ArcInfo. This is the ArcCatalog (client) licensing, not the server.


title: How to select multiple Features from multiple Layers in OpenLayers?

text: Is there a way to select multiple features from multiple vector layers by using clickTolerance in click? I checked the documents of getFeature control and selectFeature control, seems not work in this way. I know that it is possible to select multiple features by holding down Shift key. That's OK in a desktop environment, but it's not desirable in mobile/tablets. So I'm looking for a different solution.


title: How can I turn islands into individual features in QGIS?

text: I have what seems like a pretty straightforward task in QGIS: I have a shapefile that, in some cases, groups multiple separate landmasses into a single feature. For example, the Hawaiian islands are all one feature. How can I split those up into multiple features? This seems like it should be easy but it's proving impossible. Using the "Split Feature" tool doesn't do anything when I try to draw a split line between the islands, or surround one of them. The "Singlepart to Multipart" option in the vector tools might do the job, but it's not selective, it just splits up EVERY multipart feature, when I only want to split up a few of them. How can I pick a polygon that's one part of a multipart feature, like a single island in a chain, and split that out as a separate feature?


title: ArcGIS Desktop - Web Services Definitions

text: Aside from `.LYR` files, is there a file format (preferably open) for defining web service endpoints for Arc products? For example, if I publish a WMS feed, can I provide a downloadable resource for Arc users? EDIT -- Note endpoint is not ArcGIS Server - I'm looking for a simple method to allow ArcGIS Desktop users load OGC services for specific layers from a download. EDIT x 2 -- What about `.FDL`? EDIT x 3 -- It seems the `.NMF` (ArcGIS Explorer) formats can do the job. Are `NMF` files supported by ArcGIS Deskop 10?


title: How to get raster layer pixels contained in a shape defined in a vector layer

text: I am new in QGis python plugins dev and I have a "maybe stupid" question (in the title :-) ). I try to retrieve pixels (from a raster layer) which are "contained" in a polygon defined in a vector layer. How can I achieve that via a python plugin ? The user has to select the vector layer and the raster layer using the GUI of the plugin (which is already coded). Thanks in advance !


title: What needs to be done to ArcGIS Server if the server machine/host name is changed?

text: We have a Windows Server 2008 where ArcGIS Server 10 was installed. Not much has been done with the machine yet, but the name of the machine had to be changed after the ArcGIS installation. The old server name was DAYTONA and the new name of the server is PIGIRON (names changed to protect the innocent). In the server activity log there is a message: "Machine daytona is not a valid server container. Machine daytona is not found." I also see an error in the system event log that reads: "DCOM was unable to communicate with the computer daytona using any of the configured protocols." I found a reference to daytona in the file C:\Program Files (x86)\ArcGIS\Server10.0\server\system\Server.dat, but when I remove the machine reference in this file, save it, and restart the SOM, the machine reference gets put back into this file. This leads me to believe that there are other things that need to be updated because of the server name change. I found something on ESRI's website about doing this on a Linux or AIX box and the three step solution was uninstall ArcGIS, rename the server, install ArcGIS. I might end up doing that, but thought there might be a way to do it with out reinstalling everything. The SOC is not listed in ArcCatalog under this GIS Server. We tried adding it anyway, and then removing it, but that didn't work.


title: Asking about OFTInteger data type in OGR/GDAL?

text: I created the shape file by using the feature field type: OFTInteger. But when I open the file in Arcgis (ArcCatalog) the data type is: **Double (not integer)**. Do you know how to fix this (I want it is Int) ? Please see my codes:               OGRFieldDefn oSHToBDField( "SHBANDO", OFTInteger);     poSHPLayer->CreateField (&oSHToBDField);


title: How can I load stricly tiles within the current viewport in Leaflet.js?

text: I have tiles that are not available for the whole map. How can I stop Leaflet.js/Mapbox.js from making requests to these tiles? There are dozens of 404's that I'd really like to get rid of. Thank you.


title: extract country code from lat long submission - offline

text: 1) need to extract the country code 3166 from a lat long submission (example http://ws.geonames.org/countryCode?lat=-16.03&lng=167.2 2) want to localize that data set of country coordinates if I'm offline and have a cache of submissions


title: SLD MaxScaleDenominator question

text: I have quite a problem, hope you'll be able to help. I've got a very large shp layer that I'm loading via wms to my openlayers project. Due the fact that the layer is very big (2 millions polygons) It causing my geoserver to crash. I have 2 solutions to this situation:   1. I wrote a js code to load the layer only in certain zoom levels. This solution works but i don't want it to be executed via the js code, but via the sld.   2. write a MaxScaleDenominator statement. when I try this, I know it executes (I see the polygons only at the selected views), but geoserver still crashes. As I can see it, It still loading all the polygons and just after it filtering it, causing it to crash.  Any idea how can I make the zoom filter without loading all the layer ? Thanks!


title: configuring cartodb for running on local machine

text: How to configure cartodb to run fully on local machine. I have setup a development version of cartodb. In this the tables are limited to 5, the disk space is limited and the view. Is these variables can be configured to be run on local server. I have looked the database.yml and app_config.yml. Is there any docs available for cartodb server configuration on local machine.


title: Is it possible to use GDAL without admin rights?

text: I am trying to do some image analysis in Python but am unable to install GDAL due to limited rights. I have attempted copying some of the files and folders into my python folder but am constantly encountering the same error               ImportError: No module named _gdal      The searching that I have done always seems to point to some sort of installer for GDAL, but I am only able to copy the needed .py folders into my python folder. Is there a way to do this on Windows and if not is there an alternative image reader that can properly open rasters? Thank you. EDIT: sys.path results: `['C:\\\Users\\\ûser\\\Documents\\\Python Image Extractor\\\Image Interpreter\\\src', 'C:\\\Python26\\\ArcGIS10.0\\\lib\\\site- packages\\\xlutils-1.5.2-py2.6.egg', 'C:\\\Python26\\\ArcGIS10.0\\\lib\\\site- packages\\\xlrd-0.8.0-py2.6.egg', 'C:\\\Python26\\\ArcGIS10.0\\\DLLs', 'C:\\\Python26\\\ArcGIS10.0\\\lib', 'C:\\\Python26\\\ArcGIS10.0\\\lib\\\plat- win', 'C:\\\Python26\\\ArcGIS10.0\\\lib\\\lib-tk', 'C:\\\Python26\\\ArcGIS10.0', 'C:\\\Python26\\\ArcGIS10.0\\\lib\\\site- packages', 'C:\\\Python26\\\ArcGIS10.0\\\lib\\\site-packages\\\PIL', 'C:\\\Python26\\\ArcGIS10.0\\\lib\\\site-packages\\\win32', 'C:\\\Python26\\\ArcGIS10.0\\\lib\\\site-packages\\\win32\\\lib', 'C:\\\Python26\\\ArcGIS10.0\\\lib\\\site-packages\\\Pythonwin', 'C:\\\windows\\\system32\\\python26.zip', 'C:\\\Program Files (x86)\\\ArcGIS\\\Desktop10.0\\\bin', 'C:\\\Program Files (x86)\\\ArcGIS\\\Desktop10.0\\\arcpy', 'C:\\\Program Files (x86)\\\ArcGIS\\\Desktop10.0\\\ArcToolbox\\\Scripts']`


title: Slow processing with QGIS Raster Calculator- Size issue?

text: I am first trying to determine what size raster is too large for processing in QGIS. I have 2 raster datasets that are over a GIG a piece (both have ige files) and I need to see how they are different. They are both 30 meter datasets. Raster Calculator (Raster1 - Raster2). After about 10-15 min I am only on 2%. I am guessing this is a size issue, but am not sure? I have used the set extents feature of Raster Calculator and processed tiny ( difference of .3y & .3x dd) and small sections(5y x 4xdd) . Each taking approximately 5 min. However, both output files are always the same 14 gigs. I have been performing calculations with Raster Calculator with 100 meter data and have not had an issue, but those datasets are small in comparison (< 50mb). I am using QGIS version 1.8.0-Lisboa QGIS code revision 6416f38 Compiled against Qt 4.7.1 Running against Qt 4.7.1 Compiled against GDAL/OGR 1.9.1 Running against GDAL/OGR 1.9.1 GEOS Version 3.2.2 PostgreSQL Client Version 8.3.10 SpatiaLite Version 3.0.1 QWT Version 5.2.1 Suggestions? or is this just a size issue, in which case I will not be able to use QGIS to determine differences in my raster files.


title: Creating Cluster "Zones"

text: I've looked at the "Questions that may already have my answer" and it didn't appear that the questions were specific to my questions. Although I am sure I am not uncovering some brand new quetsion which no one has ever asked before. FYI I am running ArcGis 10.1; have access to PostGIS (not familiar), SagaGIS (familiar to a degree), and am using excel for my data storage and manipulation capabilities. ![map of residences in a neighborhood](http://i.stack.imgur.com/q7ANu.jpg) I have a map of a neighborhood roads, in a lion shapefile format put out by my local department of city planning, and point data of residences which has been mapped to those roads. I also have a spreadsheet in .xls format which contains further data for all of the points. **The question is: what is the most efficient way to display these points in a visual context, and one that can be used to quickly look up where certain people live?** A couple of strategies I have thought of: 1) Label the points; this does not work very well, because there are approximately ~1500 points in a small space, the map quickly becomes over-crowded with labels. 2) Create 'Zones'; IE, this zone contains points 1-10, the next one 11-20, and so on, and have another table which has 'point 3445 = john doe", "point 346 = jane doe", etc. The problem is that I would want to look up the information by name as well; to say that "jane doe lives at point 346", and I do not know how that can be done in ARCGIS, it seems to be something for excel instead. -> A further problem is that as I got the data, it was not ordered by FID, but rather the last name of the person - but that seems to be more of an excel question rather than an ARCGIS question. **I guess the question really is, whats the best way to translate Field ID to last name of person, and is ARCGIS the most efficient program to do so, or would this work be better done in excel, or other database software?** ' 'Any thoughts on the matter much appreciated


title: How to install the GDAL plugin for QGIS 2.0 on Linux Mint?

text: What do I need to access the GDAL plugin in QGIS 2.0? I've installed QGIS on my Linux Mint 16 64bit system, but can't access the GDAL Contour tool. The problem seems to lie with my package dependencies:               sudo apt-get install python-gdal     [sudo] password for user:      Reading package lists... Done     Building dependency tree            Reading state information...      Done     The following packages were automatically installed and are no longer required:     libgsl0ldbl libqwt5-qt4 libspatialindex1 libspatialite3 libxerces-c28     qgis-common qgis-providers-common     Use 'apt-get autoremove' to remove them.     The following extra packages will be installed:     libgdal1h libkml0 libspatialite5 liburiparser1     Recommended packages:     proj-bin     The following packages will be REMOVED     libgdal1 libqgis2.0.1 qgis qgis-providers qlandkartegt     The following NEW packages will be installed     libgdal1h libkml0 libspatialite5 liburiparser1 python-gdal      Obviously I don't want QGIS uninstalled, but how can I access the GDAL plugins without python-gdal? I'm sure this is a simple oversight but it's been driving me just a little bit crazy! Thanks in advance for your assistance!


title: Why doesn't CartoDB display my data?

text: I have loaded a table into CartoDB with a column for latitude, longitude, date/time and a number. CartoDB says that the data is georeferenced by the designated columns but when I open the map, I see nothing but the base map. My data is not being displayed. Can someone tell me some basic mistakes I might have made?


title: Import CSV: Is there a cap for number of row

text: I am trying to import several CSV files which contain more than 1,000 rows respectively. However, but only 129 records were only imported each file. Can't see any options that deal with this? I am using Mac(10.7.5), QGIS 1.8.0-Lisboa, GDAL/OGR 1.9.1. Anyone can advise me? Many thanks.


title: Does Geowebcache uses JAI?

text: Do I need to install JAI when I use Geowebcache? Or can I just use one or the other? I've installed JAI and uses Geowebcache, I'm not sure if I take out the JAI step it would have given the same results. Thanks


title: Creating responsive interaction with tiled maps

text: I think I may have two conflicting aims but... I need to create responsive Choropleth map very similar to http://leafletjs.com/examples/choropleth.html or http://www.texastribune.org/library/data/us-senate-runoff-results-map/ I need to take it down to county level of detail - like the texastribune example - but unlike the two examples above I cannot restrict the map to either a state or even a country. It has to have the show the detail globally. To describe the whole world in GeoJson polygons down to the county/ postal area level would involve shipping an insane amount of data to the browser which has taken me down the route of using tiles and leaflet.js to orchestrate. This works well for providing a basic Choropleth map but I can't then see a way to display the additional data when you mouse over an area. Is there any hybrid way of getting the best of both worlds


title: Python: "Import Script" greyed out

text: I have a tool that I created from a Python script which I update from time to time. Usually I will right-click the tool and choose Import Script, but lately this has been unavailable, as shown below: ![enter image description here](http://i.stack.imgur.com/jEZIa.png) This has been happening to almost all the tools I create, usually after a few imports. The problem is that if the script is cached the only workaround is to create a whole new interface, which can be pretty time-consuming. Does anyone know why this happens? I found this post on the Esri forums describing the same problem: http://forums.arcgis.com/threads/70829-Python-Import-Script-grayed-out Thanks!


title: Arcpy: calculating a field based on feature class geometry, eg 'if FC =="POINT" calculate xy coord

text: I have a script that populates new fields based on user input. I would like it to also add the XY coordinates in DD if the Feature Class is a point feature, else, it moves on to the next command. I would have thought that using the Describe function with an IF statement that it would work (see code at line 32). However, although I don't receive an error message, the fields are not populated. The script works fine without the IF statement on a point featureClass but (obviously) fails if it is not a point feature class. I've placed the entire code in just in case. Any thoughts?                   # -*- coding: utf-8 -*-     # ---------------------------------------------------------------------------     # HQ_Std_fields.py     # Created on: 2014-04-14 14:49:12.00000          # ---------------------------------------------------------------------------          # Import arcpy module     import arcpy     from arcpy import env     env.workspace = "CURRENT"     from time import strftime     import numpy     # strftime("%y%m%d")     #--------------DEFINE VARIABLES---------------         in_featureClass = arcpy.GetParameterAsText(0)         de_classe = arcpy.GetParameterAsText(1)         de_objet = arcpy.GetParameterAsText(2)         de_source = arcpy.GetParameterAsText(3)         de_remarque = arcpy.GetParameterAsText(4)     #de_traitment = arcpy.GetParameterAsText(5)     #num_sequence = arcpy.GetParameterAsText(6)          #-----------------DEFINE EXPRESSIONS-----------------     ClaExpression = str(de_classe)     objExpression = str(de_objet)     deSouExpression = str(de_source)     coDate = strftime("%Y-%m-%d")     deRemExpression = str(de_remarque)               #----------------------------------------------------     #add fields xDD and yDD with WGS84 spatial reference     #populate xDD and yDD with decimal degree coordinates.     #source: http://gis.stackexchange.com/questions/29478/how-to-calculate-xy-coordinates-in-decimal-degrees-from-projected-dataset-in-mod     #source: http://gisstudio.wordpress.com/2013/04/29/calc-the-latlong-of-a-point-in-python-using-new-arcpy-tools/          desc = arcpy.Describe(in_featureClass)     type = desc.shapeType     if type == "POINT":         sr = arcpy.SpatialReference(4326)         updCursor = arcpy.UpdateCursor(in_featureClass,"",sr)         for row6 in updCursor:             pnt = row6.Shape.getPart(0)             row6.xDD = pnt.X             row6.yDD = pnt.Y             updCursor.updateRow(row6)         del updCursor, row6          #-------------------------------------------------------     ClasseTemp = de_classe     for ClasseAbrev in ClasseTemp:                    #reads feature class name     #determins geometry and plots xy coordinates for points and centroid for polygons and lines?     #reads class abbreviation and populates class field     #-------------------------------------         if ClasseTemp == "activite autochtone":             ClasseAbrev = "Activite autochtone"         elif ClasseTemp == "agriculture":             ClasseAbrev = "Agriculture"         elif ClasseTemp == "aire designee":             ClasseAbrev = "Aire designee"              if ClasseTemp == "archeologie et patrimoine":             ClasseAbrev = "Archeologie et patrimoine"         elif ClasseTemp == "batiment":             ClasseAbrev = "Batiment"         elif ClasseTemp == "aire designee":             ClasseAbrev = "Aire designee"              if ClasseTemp == "equipement":             ClasseAbrev = "Equipement"         elif ClasseTemp == "faune aquatique":             ClasseAbrev = "Faune aquatique"         elif ClasseTemp == "faune terrestre":             ClasseAbrev = "Faune terrestre"              if ClasseTemp == "faune avienne":             ClasseAbrev = "Faune avienne"         elif ClasseTemp == "mesure d'attenuation":             ClasseAbrev = "Mesure d'attenuation"         elif ClasseTemp == "milieu biologique":             ClasseAbrev = "Milieu biologique"              if ClasseTemp == "milieu humain":             ClasseAbrev = "Milieu humain"         elif ClasseTemp == "modele numerique":             ClasseAbrev = "Modele numerique"         elif ClasseTemp == "milieu physique":             ClasseAbrev = "Milieu physique"                   if ClasseTemp == "modele":             ClasseAbrev = "Modele"         elif ClasseTemp == "vegetation":             ClasseAbrev = "Vegetation"         elif ClasseTemp == "voie de communication":             ClasseAbrev = "Voie de communication"              if ClasseTemp == "donnee auxiliare":             ClasseAbrev = "Donnee auxiliare"         elif ClasseTemp == "hypsometrie":             ClasseAbrev = "Hypsometrie"         elif ClasseTemp == "hydrographie":             ClasseAbrev = "Hydrographie"          #--------------CALCULATE TEXT FIELDS---------------------------     with arcpy.da.UpdateCursor(in_featureClass,["de_cla"]) as cursor:         for row in cursor:             if row[0] == None or row[0] == ' ':                 row[0] = de_classe                 cursor.updateRow(row)          del row     with arcpy.da.UpdateCursor(in_featureClass,["de_obj"]) as cursor:         for row2 in cursor:             if row2[0] == None or row2[0] == ' ':                 row2[0] = de_objet                 cursor.updateRow(row2)          del row2          with arcpy.da.UpdateCursor(in_featureClass,["de_sou"]) as cursor:         for row3 in cursor:             if row3[0] == None or row3[0] == ' ':                 row3[0] = de_source                 cursor.updateRow(row3)     del row3          with arcpy.da.UpdateCursor(in_featureClass,["co_sou"]) as cursor:         for row4 in cursor:             if row4[0] == None or row4[0] == ' ':                 row4[0] = coDate                 cursor.updateRow(row4)     del row4     with arcpy.da.UpdateCursor(in_featureClass,["de_rem"]) as cursor:         for row5 in cursor:             if row5[0] == None or row5[0] == ' ':                 row5[0] = de_remarque                 cursor.updateRow(row5)     del row5                    arcpy.CalculateField_management(in_featureClass,"co_ver","'V'" + strftime("'%Y-%m-%d'"),"PYTHON_9.3","")


title: Can I use python 2.6 with ArcGIS 9.3?

text: Is it possible to use Python 2.6 (or 2.7?) with ArcGIS 9.3? The official stance is no, not supported. There are a number of reports in the wild though that it can be done. How is this accomplished? and is there anything to watch out for?


title: Displaying rectified raster (topos) in Google Earth

text: Looking for tool or an easy way (limited programming skills here) to display geotiffs in Google Earth. Specifically Canadian topos like this This Geogratis site has the same topos in regular tif and pdf formats but can't seem to locate a way to display them in Google Earth. I found this script but don't know how to implement it. Any help would be appreciated!


title: What are the implications of using .NET Generics to store ArcObjects?

text: I'm wondering what are the implications of using Generic types in .NET to store ArcObjects, such as ILayer, IField, etc. The compiler throws me a warning when using these kind of values               public class Foo     {         private List<ILayer> fooLayers;              public List<ILayer> FooLayers         {             get { ... }             set { ... }         }     }      Here is the warning: Warning 15 Type library exporter warning processing 'ArcMemorialCore.Topography.IMemorialDocument.set_ProfessionalsEnvolved(value), ArcMemorialCore'. Warning: Type library exporter encountered a generic type instance in a signature. Generic code may not be exported to COM. ArcMemorialCore My concern is not only the warning existence, but of good design practices, performance, etc, also.


title: How to turn off tool tip for the Identify button?

text: When using the esriControls.ControlsMapIdentifyTool tool, the (i) button, whenever you click it, it says "Click on or drag a box over a feature and place on the map where you want to identify" in a window... Is there a way to remove this auto tool tip? It's totally unnecessary because if you close it and click on a point on the map it will respawn itself with real data. I'm using ArcEngine 10, C#, & VS2010.


title: user select fields to be deleted in a feature class

text: How do I make the script below accept parameter arguments correctly in ArcGIS 10? I added it to a toolbox as a script, created parameter called _Fields to be Deleted:_ and type as Field, but when I double click it to run it, the display dows not allow for adding fields (to be deleted)? PLEASE HELP:(               # Import arcpy module         import arcpy         arcpy.env.workspace = "c:\\test\\Mydata.mdb"         Drop_Field = arcpy.GetParameterAsText(0)          # Process: Delete Field (36)         arcpy.DeleteField_management("ArcAddress", Drop_Field)


title: Call FWTools (ogr2ogr) from ArcGIS python script

text: In ArcGIS, I want to use a Python script (and a standard ArcGIS script tool UI) to call functionality from FWTools. How would I make the call to FWTools? I tried using os.system(C:\WINDOWS\system32\cmd.exe /K "C:\Program Files\FWTools2.4.7\setfw.bat") which is what the windows shortcut calls, but that is not it. It seems I am probably missing something rather simple, but I can't figure out how to make the right call through python to get FWTools & ogr2ogr operating on the string I pass. My UI looks like this: ![enter image description here](http://i.stack.imgur.com/EZb1B.jpg) and the pseudo-code behind it is something like:               'ogr2ogr -f "' + outformat + '" "' + outfile + '" "' + infile + '" ' + ogrvars


title: Replace Polygon geometry in SDE Database using FME

text: I'm still trying to wrap my head around using FME and encountering another issue. This is my scenario: I have a master SDE database that i am trying to update using the data from a GDB file. I was able to successfully update the SDE database on the Points data but when i tried to do the same process on Polygon's, the update breaks my SDE data. FME Translation Setup:   * I have 2 readers (GDB file reader and SDE DB reader)   * I created 3 writers (INSERT, DELETE, UPDATE) *there could be a better implementation   * I am just reusing these readers and writers and created an update process for each feature inside the GDB file. Process: ![enter image description here](http://i.stack.imgur.com/Icsw5.png)   1. Read through the SDE and GDB data.   2. Run CRC on each dataset.   3. Merge resulting datasets using FeatureMerger      * Type: Attributes and Geometry   4. Run GeometryCoercer and 2DForcer.   5. Tester, check if both CRC fields are the same.      * If yes, do nothing     * If no, update the row. The translation succeeds but when i try to view the feature in ArcCatalog, it gives me this error: ![ArcCatalogError](http://i.stack.imgur.com/rVf6M.png) I might just be doing/setting something stupid that is causing this issue. **UPDATE:** Also tested the updated geometries before the insert using the Data Inspector and everything seems to be fine. I think my issue is somewhere in the WRITER.


title: Can't get dataProvider object outside of QGIS python Interpreter

text: Within PyCharm               import qgis.core as qgis     vector_lyr = qgis.QgsVectorLayer("Point", "distance nodes", "memory")     pr = vector_lyr.dataProvider()     print pr     prints None      Within Python Interpreter in QGIS               import qgis.core as qgis     vector_lyr = qgis.QgsVectorLayer("Point", "distance nodes", "memory")     pr = vector_lyr.dataProvider()     print pr     prints <qgis.core.QgsVectorDataProvider object at 0x000000000DAFF0D0>      I've tried this Getting dataProvider from vector layer outside QGIS and still no luck. Thanks. [Edit] I found this: Access Sextante (processing) in standalone QGIS app? is this the case for dataProviders? Is there a workaround? Thanks. [EDIT2] Tried this to no avail. I'm trying to create an in memory vector layer, write values to it then do something else after.               import qgis.core as qgis          qgis.QgsApplication.initQgis()          vector_lyr = qgis.QgsVectorLayer("Point", "distance nodes", "memory")     pr = vector_lyr.dataProvider()     print pr     prints None


title: How to compute percent of a field

text: I have a layer with a column conaining population values. I want to write an expression to compute in a new field the percentage of population each polygon has in relation to the total population. I know I can use basic stats to get the SUM of area column, but I want an automatic process to repeat over several fields).


title: Rotating points at intersections toward branching vector

text: I've successfully rotated points on based on the azimuth of a line, following the lead of this post. However, I now have reached an obstacle for one of my layers that occurs only at line intersections. Would it be possible to get points rotated only towards lines that branch off of principal vectors?


title: What sub-meter geolocation solutions exist and how can I integrate them into a web and mobile mapping workflow?

text: I am the manager of seven water systems. What I am trying to do is map out every meter and valve on to a map and then my workers can have the maps on their iPhones so when they have to go fix a water leak or something they can just look on their map and know where the meter and the valve is. I've plotted a few meters on my phone and then got the longitude/latitude and transferred that into Google Earth, but there are a few problems with that:   1. My phone is usually 10 or so yards off on my "current location". I need to be within a few feet to plot water meters.   2. It takes forever to take the points that I have plotted on my phone and key them into Google Earth.   3. Once they are on Google Earth my workers still do not have them on their phones so it does them no good while they are on the job. So does anyone know of a much easier way to do this? I just need to plot the points very accurately and then be able to access the points from different iPhones and computers.


title: Split Lines with Points Layer Tool

text: I asked a question a while ago to figure out how to split a line layer with a points layer using QGIS, but cannot, however, get the accepted answer to work due to limitations of my knowledge and the fact that i'm working on Windows 7. I noticed this feature request that dates from about a year ago, requesting a tool that accomplishes exactly what I need. I've been looking for it, don't see it, and am wondering if it is currently an available tool and i'm just not looking in the right spot or not.


title: how to extract network of city routes from city blocks shape file?

text: I have a shapefile which contains urban blocks. So, I am going to extract network shapefile of city routes and calculate total length of network arcs. I am using arcgis 10 for this work. Is there any tools for this work? My idea is to use arcmap tools that can process city blocks and draw routes arc in gap between block parcels of cities!! If the question is unclear please tell me to describe more. tnx


title: How to convert raster negetive pixel value to positive

text: I was working on ground water Level. I used raster math minus function and the Output raster grid data contains negetive cell values. ma resulting raster contains values ranges from High : 0,187477 to Low : -1,14927. But I do not want any negetive values.I just need to alter the negetive sign to positive vales.How can I do that using arcgis10.2


title: Segmenting road from Excel spreadsheet based on specified distance?

text: I have created a model that turns an excel spreadsheet containing coordinates into a Kml/Kmz file of a road that is segmented based on the X_from,X_to,Y_from,Y_to Columns, and the image below represents the format of the excel spreadsheet that is the models input. The model works 100 percent all the time, and is not the problem. As you can see in the image the grey and white rows and columns determine the points for the road the "Breadcrumbs", and the rows and columns in blue determine the line itself and how it is segmented. To segment the line previously I just made the To and From coordinates, each of the Nodes, I would like to find out if there is a formula, calculation, anything like that, I can use to determine the To and From coordinates based on a specified distance of my choosing. Example: I want to segment the road every 20 metres. I am using **ArcGIS 10.1 for Desktop**. ![enter image description here](http://i.stack.imgur.com/ID3UW.jpg)


title: Postgis query not using the spatial index

text: I've been using postgres8.3/postgis1.5 for some time, and it works well. I've moved my entire database to a new, faster machine and upgraded to postgres9.3/postgis2.1 using pg_dump/postgis_restore.pl, which also worked well. Everything else seems to be running faster, but a particular query is now taking extremely long to complete - so naturally I checked that it was correctly indexed and the tables were vacuum+analyzed. VACUUM ANALYZE public.places (the_geom) doesn't change the result. Using EXPLAIN, I see that the 9.3/2.1 instance isn't even using a spatial index. The public.places table is identical in both databases; same structure and exactly the same records (236,884 rows). I don't doubt that I can write a more efficient query, but before I launch into that I'd like to at least get it running as fast as it was in 8.3/1.5 (700ms vs 580000ms). Suggestions? Does anyone see my error? What has changed in postgis2.1 that keeps the spatial index from being utilized in this query? The sql:               SELECT max(foo.pname) as pname, split_part(attribute, '|', 2) as zipcode FROM places,      (SELECT features.pname as pname, ST_Buffer(the_geom, -.001) as the_geom FROM places p,          (SELECT pname, pid FROM search WHERE gid = '50873' LIMIT 1) as features      WHERE p.pid = ANY(features.pid)) as foo      WHERE places.ptype = 'Zip Code Areas'      AND ST_Intersects(places.the_geom, foo.the_geom)      GROUP BY zipcode      ORDER BY zipcode;      The table schema:               CREATE TABLE public.places     (       pid integer NOT NULL DEFAULT nextval('places_pid_seq'::regclass),       ptype text,       pname text,       the_geom geometry,       sid integer NOT NULL,       pcountry text,       woeid integer,       attribute text,       CONSTRAINT places_pkey PRIMARY KEY (pid),       CONSTRAINT enforce_dims_the_geom CHECK (st_ndims(the_geom) = 2),       CONSTRAINT enforce_srid_the_geom CHECK (st_srid(the_geom) = 4326)     )     WITH (       OIDS=FALSE     );     ALTER TABLE public.places OWNER TO xxxx;     GRANT ALL ON TABLE public.places TO xxxx;     GRANT SELECT ON TABLE public.places TO xxx;          CREATE INDEX places_pcountry_idx ON public.places       USING btree (pcountry COLLATE pg_catalog."default");          CREATE INDEX places_ptype_idx ON public.places       USING btree (pcountry COLLATE pg_catalog."default");          CREATE INDEX places_ptype_key ON public.places       USING btree (ptype COLLATE pg_catalog."default");          CREATE INDEX places_spatial_index ON public.places       USING gist (the_geom);          CREATE INDEX places_woeid_idx ON public.places       USING btree (woeid);      EXPLAIN 8.3/1.5: ![EXPLAIN 8.3/1.5](http://i.stack.imgur.com/Nzpdu.png) EXPLAIN 9.3/2.1: ![EXPLAIN 9.3/2.1](http://i.stack.imgur.com/r2CaN.png)


title: empty file geodatabase error (nowhere feature class)using arcpy

text: i want to import a python(2.7) script to arcmap 10.1 my script has 2 inputs variables one for workspace and another raster dataset and i want o convert input raster as feature class polygon . i want to create a file geodatabase and import in her a export feature class,run this script and finaly show me **complete message(not error)**... but created the basis is empty,nowhere feature class.my feature class= **Catchment** and i want to import in **gdb** (my file geodatabase) but if i run this script show me complete message but the **gdb** (my file geodatabase) is empty               import arcpy     import sys     import os     from arcpy import env as e     from arcpy.sa import *     arcpy.ImportToolbox("C:\Program Files (x86)\ArcGIS\Desktop10.1\ArcToolbox\Toolboxes\Arc Hydro Tools.tbx", "archydrotools")     try:      e.workspace = arcpy.GetParameterAsText(0) # import workspace      Cat = arcpy.GetParameterAsText(1)         # import raster dataset      ws = e.workspace      gdb = "teleia"                            # import name of geodatabase      # Local variables:      arcpy.CheckOutExtension("spatial")      arcpy.env.overwriteOutput = True      arcpy.CreateFileGDB_management(ws,gdb)      out = ws + os.sep + gdb + os.sep      # Process: Catchment Polygon Processing      arcpy.CatchmentPolyProcessing_archydrotools(Cat, Catchment)      Catchment =out + "Catchment"     except:      arcpy.GetMessages()      any idea ?


title: Reprojecting kml in NAD27, 83 or WGS84 using QGis

text: I have downloaded a kml file from Google containing location of shops. The coordinates in that file are in lat/long. I have other shp files containing roads and territories. I am trying to have all those files in the same projection and since I need to work with buffers and distances, I figured NAD27 (EPSG: 26718) or WGS84 (EPSG:4326) were good choices (units in meters). However, using QGis, I have tried to reproject the kml file with no success. Even though I can see that they have the same CSR in the metadatas, the kml file do not overlap with the other two shp files. For example, with NAD27 Zone 18 N set as the CSR, the kml file coordinates range from xMin,yMin 362234.57,45.4667 : xMax,yMax 362234.77,45.5717 while the shp file range from xMin,yMin 3271320.50,5039271.30 : xMax,yMax 305262.50,5064239.70. Sooo would anyone give me a much needed helping hand?


title: Is it possible to customize openFeatureForm of pyqgis?

text: I am a beginner in PyQgis. I am trying create a plugin, which will be used to AddFeatue,Delete Feature, and Edit the Feature. It is more customized than which is provided by Qgis (2.2). I am capturing the Points using RubberBand API and then Add the attributes for the feature using openFeatureForm. The openFeatureForm takes default values as NULL. What I want is to customize this form. Like populate some values with combobox (to allow user to select the required values). Also for fileds like date I want to put DateTimePicker widget in the form. So Is it possible to customize the form provided by openFeatureForm()? Thanks.


title: How to export raster from ArcGIS for Desktop to MATLAB using Tiff without losing projection?

text: How do I export data from ArcGIS for Desktop and import to MATLAB using Tiff without losing projection (lat lon) information? ASCII is too slow for my application. I manage to export a raster dataset (.img) to Tiff and tfw (worldfile) using ArcGIS for Desktop. However, when I imported the data in MATLAB, I lost my projection. I tried both `imread` and `geotifread`.


title: How to create SQL statement to only return chainage values every 500m?

text: I have a list of chainage values every 100m (as a whole number) and I am trying to create an SQL statement to only return values every 500m. Any ideas? Using ArcGIS 10.1


title: Read .DEM (USGS) FIles with C#

text: I have a huge number of .DEMs and have been investigating how to process them into something useful. Does anyone know how, or even if it is possible to read coordinates and elevation data from .DEMs using .NET? Specifically C#. Thanks in advance. C


title: Creating MBO files with MapBasic and Notepad++

text: [Note: I have been using Notepad++ (with some customisations as per this page), rather than MapBasic itself. This may make things a bit trickier, but so far the time saved in development has been worth the extra hassle during setup.] I am trying to start using project files with MapBasic. As I understand it the process is as follows: compile your modules from .mb files into .mbo files. Then process the .mbp file which links the .mbo files together and, if there are no errors, creates a .mbx. This makes sense to me except for one thing...how do i get my .mb files to compile to .mbo instead of .mbx? What am I missing?!


title: Formula to get polygon centroid

text: I need to compute the centroid position of an irregular polygon with lat and lon coordinates in Java. Which formula should I use to correctly compute a centroid?


title: How to draw nodes of a graph in Python and link them to a table?

text: I would like to draw some points with a Python GUI, like a graph, and link those points to a table.


title: Is there a public shapefile with the London Underground stations and lines?

text: Is there a public shapefile containing the stations and lines of the London Underground? I've been doing some intensive searching but I can only find the coordinates of the stations, or the lines in way too much detail on OpenSteeetMap. Ideally, I want each line to have a single track, not multiple like with OSM.


title: How to create routable network by pgrouting with assign_vertex_id?

text: I'm trying to create the routable network as per http://underdark.wordpress.com/2011/02/07/a-beginners-guide-to-pgrouting/and use assign_vertex_id. Do I still need to create the node table? How would code for creating network look like?


title: Copy Raster Catalog from .mdb to .gdb, but only within a certain extent with arcpy?

text: I am stumped again, this website has really helped me before so here is another shot. I am trying to copy a raster catalog from a personal geodatabase to a newly created file geodatabase but only within a certain extent. Using the script below gives me a new file geodatabase with a new raster catalog but no rasters in it. Any help is greatly appreciated!!!               # Import arcpy module     import arcpy          # Set Geoprocessing environments     arcpy.env.extent = "-119.142433995 54.1505080220001 -118.783607274 54.2814004220001"          arcpy.env.overwriteOutput = True          # Local variables:     NewGDBOutputDirectory = "C:/GDB"     InputRasterCatalog = "C:/GDB/old.mdb/oldrastercatalog"     OutputRasterCatalog= "C:/GDB/new.gdb/newrastercatalog"          # Process: Create File GDB     arcpy.CreateFileGDB_management(NewGDBOutputDirectory, "new", "CURRENT")          # Process: Copy Raster Catalog Items     arcpy.CopyRasterCatalogItems_management(InputRasterCatalog, OutputRasterCatalog, "", "0", "0", "0")


title: How can I show custom titles on my table fields with OpenLayers?

text: I have a server with Geoserver, Apache and OpenLayers and it displays the map okay; I can see my data with an `OpenLayers.Popup.FramedCloud` call and it shows the content of my layer. I want to know if it's possible to show custom titles for my data instead of the ones in the table. Example: I have **COSTRUC** field and I want to show **Construction** So there is any chance to do it with Openlayers?


title: How to find embedded Coordinate Reference system in a Layer

text: I want to find CRS in an added shape file. Ex. In Arc GIS if we go to Layer Properties and Source we can find embedded CRS to a layer. How to find in QGIS, embedded CRS ia a layer/ Shape file??? I think i QGIS, if we move to Layer properties and General, it doesn't show extract CRS of layer. _Because I Added shape file with undefined CRS. In that case it shows EPSG:32244 - WGS 72 / UTM zone 44N. But I added undefined one????_ I'm using qgis 1.7.4, System is Windows 7 **How to find in QGIS, embedded CRS ia a layer/ Shape file???**


title: Algorithm to find the polygon ID where a point belongs to

text: I am wondering what is the simplest way to perform a look up of polygon Id based on the point location. Input: a point with coordinates, a table of polygons Output: the polygon Id where a point is located. I am thinking to preload the polygon table into PostGIS (The size is about 220,000 entries). What is the best way to process all these data? Thanks.


title: Why do I get wrong length using STLength on MultiLineString?

text: I have 2 MultiLineString representing the a route (first is driving there and the second is the return way) .               MULTILINESTRING ((34.787315938035292 31.318235314177119, 34.788431415005249 31.336622682320776, 34.794144933818309 31.355201806839506, 34.794207563998924 31.371163684259038, 34.792967344484573 31.375581921229493, 34.789226146087763 31.390623411898762, 34.788274952156939 31.408575008023988, 34.786368903544009 31.42299930464733, 34.786328720608147 31.42193697540446, 34.786297827387109 31.422105913663092, 34.786280312670428 31.421992838029546, 34.786190597348 31.421582726711659, 34.786161861793758 31.421905927160644, 34.784898002487694 31.424062281290652, 34.773401914392345 31.434990008872095, 34.765590772924938 31.446642271859293, 34.77342394487551 31.466799593779264, 34.77975547041099 31.490976936335215, 34.782630232155611 31.506823501958817, 34.792679583505226 31.531054387260554, 34.800484354694987 31.551364145061932, 34.806546173040267 31.574946250797954, 34.813199698392914 31.597257951527482, 34.814700882755012 31.619577056013927, 34.809541545607935 31.657018778594132, 34.80638283046715 31.674769044083895, 34.815697398246741 31.69878758883678, 34.821170796592632 31.721351512377144, 34.822676820645704 31.740659141022942, 34.826406253975371 31.760366608188313, 34.832777976436191 31.795874550345438, 34.840133770521163 31.803493762964539, 34.857643693856282 31.809682398342275, 34.876558349351328 31.809001817509152, 34.899775788981991 31.810874151671349, 34.915550260013873 31.814641705371283, 34.92078885084598 31.819327945165384, 34.94872743329605 31.826006370096646, 34.968653641494242 31.831077555261356, 34.977436258261335 31.8348796713872, 34.988078880595516 31.834460104995458, 34.99298330536876 31.826840136006965, 35.000267546905079 31.826135241324245, 35.005568854286416 31.825923237904671, 35.017740547302338 31.8215162498958, 35.02155196628523 31.818288829834177, 35.026913523519241 31.813341112120234, 35.03143746024918 31.809795839880959, 35.038159534410212 31.80619724694477, 35.045420139243944 31.804745994112629, 35.053171381816789 31.804966181835955, 35.06486714082564 31.804141661709476, 35.075841806779032 31.800281337346696, 35.088383165072827 31.800341312456926, 35.106401658005012 31.801455597207628, 35.126101243739683 31.800681742893538, 35.146976345262694 31.797221642316991, 35.162196277950279 31.793296817012131, 35.171594439546396 31.797658570113885, 35.187671634314036 31.797928602082614, 35.188362774999014 31.796533416394119, 35.1887256549307 31.7961346714363, 35.190472346375131 31.794962708212374, 35.192562710517549 31.794612374668478, 35.192707852005682 31.793637348918882, 35.193583458825614 31.793129756791721, 35.194747974527196 31.793406487587923, 35.196095096006964 31.79164621288934, 35.196955729849186 31.790773244292048, 35.1978874456438 31.790160915552967, 35.198711620213004 31.789879413000829, 35.199295088482693 31.787229378909604, 35.19850752277577 31.785375424789216, 35.199416396367077 31.784687679081976, 35.199241688426945 31.782313977270359, 35.198391385845284 31.781402127099458, 35.201245278823734 31.776056353428189))          MULTILINESTRING ((35.200768733043091 31.77473722793188, 35.204546738625524 31.774719539653546, 35.208965075961345 31.776185106153822, 35.208640148891355 31.782029315838962, 35.209755124306163 31.782465782224509, 35.211377925506667 31.781691930524566, 35.213549710915274 31.780813498217512, 35.214769003020763 31.78168291327713, 35.214692303698442 31.780946053301573, 35.212572272348339 31.781205360240914, 35.208643904703017 31.783054725633296, 35.208298951212782 31.783020204739561, 35.203248270523765 31.782224483116469, 35.198547661079523 31.786334626536405, 35.199845310730879 31.787516017026629, 35.200144134813939 31.788041155898213, 35.197250692115119 31.790485693913038, 35.191107763219378 31.795067815242238, 35.181610401764608 31.801656830805612, 35.1687011405429 31.795676931971681, 35.159244819963 31.795794581845737, 35.119424750248832 31.798775201079778, 35.080393585902307 31.800856179659302, 35.063209196166355 31.805088407881986, 35.048056788398441 31.804031639631894, 35.032927570632985 31.808478772457441, 35.020924279147472 31.81877639146, 34.990129787588472 31.829259174717787, 34.988047951875991 31.839040689444193, 34.977978798901404 31.835101996624207, 34.9498888706375 31.825967830705292, 34.92736197558343 31.823111392882268, 34.91695636557349 31.815625240675068, 34.903688415937665 31.811906632219383, 34.882969413185968 31.810547012219338, 34.86002652341891 31.810052256918588, 34.8400240232646 31.803567891209664, 34.830159736601907 31.79428193887949, 34.832461959577664 31.781116935605244, 34.827114623267633 31.762582638502867, 34.822782721105142 31.743736272787874, 34.820540303486339 31.720919484741522, 34.814209114195862 31.695453086060326, 34.806275341573269 31.675224452093136, 34.81001464412418 31.653220249705704, 34.814807584458386 31.623882775804965, 34.813291228448193 31.60270186644005, 34.8080051217612 31.582810947810554, 34.803038127301519 31.560726326313713, 34.797198877391246 31.54151166321618, 34.786754234512 31.523229452759789, 34.782952679629339 31.504485052320216, 34.776143426121912 31.483992375769422, 34.77342018413583 31.467251749582854, 34.766455832382029 31.450437751509558, 34.765905289603033 31.442231443933725, 34.785259642730011 31.422175881020067, 34.789073528508013 31.401593669076053, 34.791297259262187 31.381621922859551, 34.795717357038875 31.365271563573277, 34.792079987685526 31.350441990935515, 34.786453724639145 31.329829087974357, 34.788396830251187 31.309436171497104, 34.790810764157918 31.289914211434404, 34.7976451665849 31.287943033661829, 34.80366544178834 31.286754464613146))      The 1st route is around 95km and the 2nd is 124km (approx). When doing STLength() on them I get the following results :               0.909006861532263     0.983569104560351      Any idea why ? And how can I get the correct length of the route ? Interesting case - when I serve the line to OpenLayers and I do feature.geometry.getLength I get the following results :               111km     120km      * * * Here is how I got the length using OpenLayers: I've upload the feature from SQL to OpenLayers using OpenLayers.Format.WKT() The data in sql is in EPSG4326 (and so the OSM map) And then I did the following :               layer.features[0].geometry.getLength()      ![enter image description here](http://i.stack.imgur.com/L0xOt.png) Thing is that the length is still not correct. More over , How can I calculate it correct on SQL ? What do I need to do ? (The geometry is in EPSG4326 as I mentioned before)


title: How to open an ArcGIS Desktop .lpk file using QGIS on a Mac?

text: I downloaded an LPK file from the ESRI site that contains USPS defined zip codes. Because I'm working with QGIS on a Mac, I converted the file to a .zip. There's a .lyr file that opened up in the zip, but I'm not able to open that up on QGIS. Is there a way to convert the contents of an ESRI layer package file (.lpk) into a shapefile with the associated data, and that I can open with QGIS on a Mac? I looked into using utilities like 7-zip, but it doesn't seem to work on a Mac. Apparently one of the .xml files in the layer package should contain a link to the actual data under the packagelocation tag, but I wasn't able to find that tag in the .xml files. Does anyone know if that tag name might have switched? I downloaded the layer package file from **here**.


title: How to use arcpy.setValue for Title casing?

text: I am trying to figure out how to use fieldName.title() method while updating rows using updateCursor for multiple fields. I was trying to improve the code of this post to fit my purpose: arcpy.UpdateCursor for multiple fields               def ChangetoTitleCase(sourceTable, fieldNames):           try:             rows=arcpy.UpdateCursor(sourceTable)             for row in rows:               for fieldName in fieldNames:                 uField='"%s"' % fieldName                 row.setValue(uField, "uField.title()")                 rows.updateRow(row)             del row, rows                except Exception, e:             import traceback, sys             tb = sys.exc_info()[2]             print "Line %i" % tb.tb_lineno             print e.message      I am working with Python 2.6.5, Arcpy 10 and Win 7 64. Please let me know where I am wrong? Thanks in advance for any help!!!


title: Convert lat/long four corners data to raster

text: I have a data set in the following format.each row represents the coordinates of the four corners of a cell and the last column is the value associated to the cell.               lon1 lat1 lon2 lat2 lon3 lat3 lon4 lat4 TotHg     ========================================================     -44.48, 49.88,  -44.69, 50.51,  -43.71, 50.64,  -43.5,  50.02, 8.50          -42.53  50.14   -42.71  50.77   -41.71  50.88   -41.54  50.25   9.02          -40.56  50.35   -40.71  50.98   -39.7   51.08   -39.56  50.44   9.21          -38.56  50.52   -38.68  51.16   -37.66  51.23   -37.56  50.59   8.75          -36.55  50.64   -36.64  51.28   -35.61  51.33   -35.54  50.69   8.53      I would like to convert it to a raster to be used in GIS. At the moment I do not know the coordinate system but I am investigating. All I know is that the resolution is supposed to be 2.5°. How can this be done in either R or in ArcMap? I have read a similar thread but I am not sure if it is the same situation since I have the coordinates for the four corners of each cell.this is the other question: How can I convert data in the form of lat, lon, value into a raster file using R? I have to resuscitate this question because I have received some information about projection. I am not sure if this will change the answer or help in any way. I am still quite confused at the moment. Here is the information I have received: My model is define on a 2-d plane by using a polarstereographic projection with a griddistance of 150 km at 60N. The projection is define by: X_in_grid_units=x_coordinate_of_northpole+6370km/150km*(1+sin(60N))*tan(45-latitude/2)*sin(longitude-(-32)) y_in_grid_units=y_coordinate_of_northpole-6370km/150km*(1+sin(60N))*tan(45-latitude/2)*cos(longitude-(-32)) where 6370km is the used diameter for the projection, 150km as told before is the griddistance at 60N, -32=32W is central longitude for the projection, eg. the longitude parallel with y-axis. The griddistance as a function of latitude is defined as d=150km*(1+sin(latitude)/(1+sin(60N)), eg. 160.77km at northpole, and 80.38km at equator.


title: Issue loading GML with Openlayers no obvious error

text: i'm newbie using openlayers, there are no error in my firebug,   but, my GML can't show in my browser, what should i do? this is my GML from result of converting my shapefile with ogr2gui :                   <?xml version="1.0" encoding="utf-8" ?>     <ogr:FeatureCollection          xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"          xsi:schemaLocation="http://ogr.maptools.org/bataskabupaten_poly.xsd"          xmlns:ogr="http://ogr.maptools.org/"          xmlns:gml="http://www.opengis.net/gml">       <gml:boundedBy>         <gml:Box>           <gml:coord><gml:X>478266.71875</gml:X><gml:Y>9225540</gml:Y></gml:coord>           <gml:coord><gml:X>529920.5</gml:X><gml:Y>9291891</gml:Y></gml:coord>         </gml:Box>       </gml:boundedBy>            <gml:featureMember>         <ogr:batas_kabupaten_poly fid="batas_kabupaten_poly.0">           <ogr:geometryProperty><gml:Polygon><gml:outerBoundaryIs><gml:LinearRing><gml:coordinates>512850.53125,9268024.0 512921.875,9267914.0 512966.40625,9267855.0 513016.0625,9267800.0 513045.09375,9267757.0 513064.40625,9267704.0 513089.5625,9267643.0 513184.46875,9267479.0 513202.0625,9267448.0 513261.875,9267363.0 513375.9375,9267200.0 513443.3125,9267104.0 513483.9375,9267049.0 513549.03125,9266957.0 513584.625,9266902.0 513665.09375,9266795.0 513732.75,9266705.0 513816.84375,9266578.0 513863.21875,9266514.0 513953.5625,9266381.0 513975.65625,9266348.0 514071.03125,9266243.0 514125.59375,9266176.0 514181.0625,9266088.0 514204.59375,9266056.0 514278.8125,9265970.0 514394.25,9265813.0 514424.59375,9265777.0 514456.125,9265746.0 514478.9375,9265712.0 514505.34375,9265662.0 514542.375,9265609.0 514570.5625,9265581.0 514603.5,9265547.0 514683.65625,9265463.0 514750.5625,9265393.0 514913.75,9265245.0 514971.4375,9265201.0 515054.34375,9265119.0 515124.21875,9265049.0 515177.5,9264984.0 515229.71875,9264924.0 515343.90625,9264818.0 515399.78125,9264763.0 515462.40625,9264718.0 515640.28125,9264551.0 515788.46875,9264438.0 515829.5,9264407.0 516030.40625,9264271.0 516166.1875,9264196.0 516248.625,9264151.0 516327.28125,9264103.0 51... </gml:coordinates>     </gml:LinearRing></gml:outerBoundaryIs></gml:Polygon></ogr:geometryProperty>               <ogr:FIRST_PROP>Jawa Tengah</ogr:FIRST_PROP>               <ogr:AREA>1582363704.750</ogr:AREA>               <ogr:PERIMETER>264009.826</ogr:PERIMETER>               <ogr:ACRES>391009.023</ogr:ACRES>               <ogr:HECTARES>158236.370</ogr:HECTARES>             </ogr:batas_kabupaten_poly>           </gml:featureMember>         </ogr:FeatureCollection>


title: Editing 2007 MS Access .accdb using OLE DB connection?

text: Can you edit a 2007 MS Access .accdb file format using the OLE DB connection? I have made the connection using the following: Connecting to a 2007 Microsoft Access database (.accdb) in ArcGIS I can add and view the tables from this file in ArcMap, but when I try to start an edit session on this table, it errors out with the follwoing 2 errors:   * No editable layers   * Check to see that you have the appropriate permissions I created a .accdb file from scratch, so it is not read only. Under the Advanced tab in the Data Link Properties for the OLE DB connection, I have set the ReadWrite option check mark and unchecked the rest and I am still unable to edit this table in ArcMap. ![enter image description here](http://i.stack.imgur.com/kJAYa.jpg) Click image for larger version Name: OLE_DB_Advanced.jpg Views: 14 Size: 44.7 KB ID: 25411 There isn't much documentation on the OLE DB connection in the help file other than how to connect as posted above. Can anyone point out what I may be doing wrong? Thanks, Mike


title: How do I initialize the ArcGIS proxy via server side .NET code (C#)?

text: I know how to do this via JavaScript code:               esri.config.defaults.io.proxyUrl = "../proxy.ashx";      But, how do I do this in .NET? I am very new to the GIS space.


title: Trying to install SpatiaLite 4.3 with SQLite 3.8 on Windows 7

text: I'm new to the open-source GIS database world, and am trying to have SpatiaLite installed on my Windows 7 machine. I downloaded SQLite 3.8.2 and unpacked the 3 files in a newly created folder c:\Program Files\SQLite. I added the folder to the PATH env variable, and managed to get that going. To install SpatiaLite, I went to webpage https://www.gaia- gis.it/fossil/libspatialite/index, downloaded libspatialite-4.1.1.zip , which I believe is the latest SpatiaLite version for Windows. I unpacked this zip- file into the aforementioned SQLite folder, and then tried to run .load 'SpatiaLite.dll' from the SQLite prompt. But this gives me an error. Investigation shows that there is no such thing as 'SpatiaLite.dll', even though I thought I was just following the download and install instructions. Where am I getting wrong?


title: SLD file already exists

text: Using Geoserver 2.1.3 I do delete sld using rest api with manually delete the sld file. But when I issue a POST new sld with the same name I got an error "SLD file test.sld already exists." Event there is no test.sld in the styles folder.


title: How to convert images from the South African co-ordinate system to WGS84 lat/long?

text: I have a number of orthophotos for an area in South Africa, but they are projected in the South African co-ordinate system which measures a location in meters from a central (uneven) meridian. These particular images should use Lo31. However, QGIS does not seem to recognise the South African co-ordinate system and I am battling to display the images correctly. Any ideas how I can solve this dilemma?


title: General GeoDjango Opinion

text: I'm forced to crawl out of my intranet and create a public GIS site with quite an extensive dataset. All of my data is contained within PostGIS. I started dabbling with Geoserver/Postgis/OpenLayers and find it quite easy to work with but perhaps a little more difficult to merge into a full blown website with security, etc... So I've come across GeoDjango and I have enjoyed the python programming I've undertaken thus far. Has anyone used a GeoDjango/Postgis/OpenLayers combination? How fast is GeoDjango in serving data? Do you know of any professional websites I can browse that are built with GeoDjango? I've found it hard to find any demonstrations on the net.


title: Implementing Select Features by drawing on the Map in ArcGIS iOS SDK

text: I am developing an app using ArcGIS iOS SDK 10.1.1 with help of some map services hosted on ArcGIS server. I need to build a following functions on my app. User select a tool and draw a figure like rectangle, polygon etc, then a particular layer features are selected on the map, which are contained / touched the user drawn figure. Can anyone figure out how I can implement this?


title: Which of EGM96 geoid or WGS84 ellipsoid fits the earth better?

text: ![Figure 1](http://i.stack.imgur.com/QHX85.jpg) ImageCredit **Question** If we separately add up |H| and |h| over the earth surface, which of the sums come out smaller? * * * **Background** There is a phrase in the tag description of "ellipsoid" that says "... the geoid, the truer figure of the Earth ..." I am not sure "true" here means less absolute error but I am interested in this question nevertheless. The following passage seems to suggest that the geoid may not approximate the earth surface particularly well in the land areas: "The dominance of the low- harmonic degrees in the geoid power spectrum indicate that the dominant shape of the geoid is controlled by structures deep within the mantle."