title: Why aren't the info windows showing up on my cartodb layers?

text: I made a cartodb map with 2 layers. The layers draw from the same table but I use filters to display some of the data in one layer and some of the data in another. Without the layers my info windows were working and now they're not. I'm not sure how to start troubleshooting this. The SQL expression on layer 1 is `SELECT * FROM table_name WHERE useful_column > 0` and on layer 2 it is `SELECT * FROM table_name WHERE useful_column IS NULL` -- the goal is to use a different info window altogether on the regions where `useful_column` has no data.


title: How to read out the geoPDF data downloaded from USGS ned database?

text: I am new to GIS, so this might be a stupid question. Google foo din't get me anywhere, so here goes.. How do I read the geoPDF data that I download from the National elevation database from USGS? All I am looking for is lat/long/height for each point as much as the resulution would allow for a specific area in the map. Are there ways to get this data in a better/more friendly format than geoPDF? (I am good with python, so any solutions that are py-friendly would be awesome, feel free to toss some search words at me that would help me further my knowledge) thank you for your help


title: Does Collada allow to georeference the 3D coordinates within a DAE file?

text: Can I assign lat/long/height values to each vertice of a Collada 3D model?


title: How to use Time as a parameter in Network Analyst Service Area and Closest Facility

text: I have a road network downloaded from Ordnance Survey (Meridian 2, with A roads, B roads and Minor roads Merged into one shapefile and then built into a Network Dataset for use in ArcGIS' Network Analyst) and would like to use it to determine Service Areas for 4 hospitals (preferably with drive times of 30 minutes and 60 minutes) and Closest Facilities (same 4 hospitals and directions from 539 points). I've used the Network Dataset to find the Service Area and the Closest Facility but the only parameter I can use is distance (metres). Question: **How do I set up the network to use drive time, not distance, and work out Service Area by, say, 30 minutes?** When I try to create the Network Dataset from the downloaded road network, I try to add Drive Time as an Attribute but I can't get the Network Dataset Wizard to accept the evaluators. I've tried to follow the ESRI help pages and I've Googled around but I'm lost. I'm not sure what the time setup should be or if it's something I can just do or something I need to set up from the start. ![enter image description here](http://i.stack.imgur.com/5B2DG.png)


title: How to measure the accuracy of latitude and longitude?

text: I am very new to GPS stuff. Suppose I have latitude and longitude as `19.0649070739746` and `73.1308670043945` respectively. In this case both coordinates are `13` decimal places long, but sometimes I also get coordinates which are `6` decimal places long. Does fewer decimal points affect accuracy? What does every digit after the decimal place signify?


title: How can I tell the accuracy of coordinates conversion by cs2cs?

text: For example, when I try this:               cs2cs +proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs +to +proj=tmerc +lat_     0=22.31213333333334 +lon_0=114.1785555555556 +k=1 +x_0=836694.05 +y_0=819069.8 +     ellps=intl +towgs84=-162.619,-276.959,-161.764,0.067753,-2.24365,-1.15883,-1.094     25 +units=m +no_defs      When I put > 114.1689 22.4518 It outputs > 835447.18 834705.40 3.11 Can I tell from the result, the accuracy and possible numeric error of the result? Some meters or some percent?


title: Precision of x decimal points in meters?

text: Is there a handy reference that corelates the number of decimal points in a coordinate to the precision in meters? I had discussions like "4 decimal places will place you within a 500 meter radius" and "6 decimal places will place you in a house". I want a reference that can clear up such facts and requirements.


title: Can you suggest some well-written introductory texts about coordinate system projections?

text: Even after some time working in this GIS field, I still don't have a solid understanding of spatial references, how they are projected, and most importantly, **if you still need to worry about the fact that coordinates are geographic once they have been projected.** In which cases should one pay attention when handling coordinates in a projected coordinate system? For example, assume I am working in a projected coordinate system. I have a polygon whose vertices are now specified with X / Y coordinates. Let's say I want to manually calculate the polygon's area -- could I now take any normal algorithm (for Cartesian/Euclidean coordinate spaces) for computing the polygon's area, or do I still have to account for the fact that these are _geographic_ coordinates? What would I have to do to get an accurate result for my area computation? It's questions such as the above example that I still don't understand. Any cues to introductory texts explaining such details are welcome. (I was already considering looking at the PROJ.4 catrographic projection library and how it is used to deduce from that how to deal with projected coordinates, but that would've probably been too early.)


title: Publish Vector and Raster (GeoTiff) layer in Geoserver SLD

text: ![enter image description here](http://i.stack.imgur.com/iJdMH.png) ![enter image description here](http://i.stack.imgur.com/b2VzS.png) Base layer and overlay layer are same when opening in ArcMap But When publishing the same layers in geoserver , distraction of 27 meter. I have attached an example image to show the problem.


title: Fixing address on TeleAtlas, how long before Google Maps has this information?

text: I recently submitted information to TeleAtlas and they told me it was accepted and has been added. When I search on Google Maps the address still does not appear. Does anyone know how long Google Maps take to update with new information or it just a matter of TeleAtlas making the data available for Google?


title: Calculate area with FME

text: I have coordinates in decimal degrees and need to calculate area for polygon objects. I tried to use AreaCalculator and set up criteria to find all these bigger than 100 sqm. But result is not as I expected. Should I convert first from decimal degree to metric?


title: QGis Save Raster as Rendered Image

text: In QGIS 1.9.0 Master, when you right click on a raster in the Layers Panel and select "Save As", you can select the output mode to be "Raw Data" or "Renedered Image". When selecting the Rendered Image option is saves the raster with the layers current styling. How would I do this via the python console? I am able to script gdal_translate, however I am not sure how to preserve the current styling for the layer. Any help would be much appreciated. Regards Gavin


title: Get Intersecting roads (ArcGIS Server10 with silverlight API)

text: I am developing a tool where when I click on a road in the map I have to list the roads intersecting the selected road using ArcGIS Server10 with silverlight API to develop this tool I got the roadID of the selected road using IdentifyQuery Then I got the envelope of that road using QueryTask I passed that road envelope to IdentifyTask again and queried the intersecting roads but its lists me all the roads inside the envelope but I need only intersecting roads also its very very slow :( I am stuck here Please guide Please guide


title: How to design a database for continents, countries, regions, cities and POIs?

text: I'm brand new to GIS programming and I am designing a GIS application. Target is to create system with continents, countries, regions (including states, sub-regions, provinces), cities and places in cities. Each of this elements will contain some text information and related stuff. As database we are going to use PostgreSQL with PostGIS. My question is how to design database for this system? I was thinking of 2 tables `polygons` and `points`, but I'm not sure if it's good way of thinking. What we need is some hierarchical base for relationship between countries, cities, regions, etc. Main goal of the application will be collecting tourist data from many sources about specific cities, regions, countries and so on, and store it in database. Let's say we will have city Rome, we collect some info about this city into database from couple sources. And we need to know that Rome is in province Rome, sub-region Lazio in region Lazio, country Italy. So system should be flexible to allow create such relation from real world.


title: Conversion of coordinates of "Bochum coordinate system"

text: I need to convert historic coordinate data into more recent Gauß-Krueger data. Specifically, I need to convert "Bochum coordinates" used till 1960 in the local coal mines into DHDN/Gauß-Krüger zone 2 (EPSG 31466, ID 2646). Is there in QGis a tool (similar to that that seems to exist is ESRI ArcGIS) that I can use for this job?


title: Workspace or data source is read only

text: I am using ModelBuilder in ArcCatalog 10.1 to convert an ASCII files to TIFF and add spatial reference. I have hundreds of files to convert so I'm using the Iterators. After doing 10-20 files successfully I get the following error.               Executing (ASCII to Raster): ASCIIToRaster "D:\Users\Documents\#Data\# New Data\LIDAR\todo\nz4626_dsm_1m.asc" "C:\Users\Data\# New Data\Lidar2.gdb\nz4626_dsm_1m" FLOAT     Start Time: Thu Aug 15 14:59:06 2013     ERROR 999999: Error executing function.     Workspace or data source is read only.     Workspace or data source is read only.     No spatial reference exists.     ERROR 010067: Error in executing grid expression.      Failed to execute (ASCII to Raster).     Failed at Thu Aug 15 14:59:06 2013 (Elapsed Time: 0.00 seconds)      After using command prompt to remove read only tag from the source and destination folders and restarting ArcCatalog the process can be run again for a few more files before failing again. What is causing the folders to become read only and how do I stop it?


title: Display 3d building using OpenLayers

text: Is there any way to display a 3d building model as a layer under OpenLayers ? i have the building models drawn using Google sketchup , is that useful ?


title: Divide map into features with invisible divisor lines

text: I'm working with QGIS 1.8.0. Is there a way to divide a map with a grid but make the divisor lines invisible? I have a shapefile as the one in the left and I want it to see as the one in the right but still having the divisions (with the 82 different features). Please let me know. Regards! ![actual images](http://i.stack.imgur.com/3N3l9.png)


title: Adding a scale using SDK GXP OpenGeo Suite

text: Followed the tutorial using the OpenGeo suite SDK, and looking at adding a scale to my maps however the tutorial doesn't cover adding widgets only plugins. So I am not sure about how I would add http://gxp.opengeo.org/master/doc/lib/widgets/ScaleOverlay.html wherever I place it fr example the tools section I seem to just get a white screen, which tells me something is wrong. I believe I should call it using `* @require widgets/ScaleOverlay.js` and then as with plugins something like (below) in the tools part of app.js                {     ptype: "gxp_ScaleOverlay",      actionTarget: "map.tbar"     },


title: ArcGIS in PhoneGap

text: I am trying to set up a PhoneGap application with an ArcGIS Map using a Windows Phone template in Visual Studio 2012. I have tried using an example from the ArcGIS Javascript API samples to incorporate a Mobile ArcGIS.com map in my application: http://help.arcgis.com/en/webapi/javascript/arcgis/jssamples/mobile_arcgis.html This has not worked. I am getting the following error in Visual Studio when I try to link to the external CSS files: > Unable to edit > 'http://serverapi.arcgisonline.com/jsapi/arcgis/3.4/js/esri/css/esri.css'. > The style sheet must be part of the current project and the project must be > a Web Site or Web Application project. I found a similar example of setting up a PhoneGap project for Android with an ArcGIS Map, but this has also not worked for me. Does anyone know of a working sample PhoneGap project to incorporate a map, using the ArcGIS Javascript API? Alternatively, does anyone know how I can get around the error with linking the stylesheet?


title: Why doesn't dissolve finish? It just stops without error message

text: I have a shapfile of the whole UK that contained Postcode area, district, and full postcode. I was going to try and dissolve the regions into seperate business areas by joining a lookup table that identified, for each postcode which business area it's assinged to. But i couldn't get the business area to appear in the dissolve popup window. So instead i modified the DBF file, adding the new business area field. Now when i go into the dissolve function, the dissolve field does appear there but i can't get it to succesfully finish. It just stops. No error messages. Nothing. I think what i'm trying to do is simple but QGIS just doesn't seem to be happy. This is the very first thing i'm trying to do with QGIS so please be patient with me.


title: How to make line with arrow head in QgsVectorLayer using python?

text: I need to draw line with arrow head or after drawn the line i need to make it as arrow head line by using python. or can we set any arrow head symbol for line features by using python? Kindly help me with some sample codes.


title: Changing field type of Excel attribute table for a Join

text: I have a string field type of shapefile already in ArcMap. I need to import an excel file in and join with the shapefile. The only problem is the field type of the shapefile is a String and the field type for the excel file is a Double. The excel field is all numbers, but the shapefile field is mostly numbers and a few letters mixed in, I think thats why it came out as a string. I know you can't add a field or use the field calculator into a table exported from excel, I'm currently using it as a .csv.


title: Why do I get a NameError: name 'gdbSuccessCount' is not defined?

text: I'm using a 412-line script that gathers and merges online tile data. However, I'm receiving a Traceback error, NameError: name 'gdbSuccessCount' is not defined. I've included the portions of the code that gdbSuccessCount is included in (only 3 times it appears as is). I, and my work, will appreciate the help to find the resolve the problem. I apologize for the code mess, the indented code is below it (first time on StackExchange). The problem is the last occurrence, in line 406 as follows: if merge and gdbSuccessCount > 1: Here are lines 333-412:                   r = 0         if GDB:             gdbSuccessCount = 0             for g in gdbList:                 tileGDBID = g                 q250k = "q" + string.split(tileGDBID,"-")[0]                 srcGDB = dataPath +"/" + q250k + "/geodatabase/"+ tileGDBID +".gdb.zip"                 destGDB = dest+"\\geodatabase\\"+ tileGDBID+".gdb"                 destGDBzip = dest+"\\geodatabase\\"+ tileGDBID+".gdb.zip"                 if not os.path.exists(destGDB):                     arcpy.AddMessage("Copying Geodatabase " + tileGDBID)                     if not os.path.exists(destGDBzip):                         try:                             urllib.urlretrieve(srcGDB, destGDBzip)                             gdbSuccessCount += 1                         except:                             arcpy.AddMessage("\t WARNING: Geodatabase for tile " + tileGDBID + " not found - moving on to next one...")                      if os.path.exists(destGDBzip):                         zf = zipfile.ZipFile(destGDBzip)                         zf.extractall(dest+"\\geodatabase")                         del zf                         os.remove(destGDBzip)                 else:                     arcpy.AddMessage("\tGeodatabase " + tileGDBID + " exists, moving on to next tile...")              if LAS:             # this is a batch file that will be used to unzip the LAZ files....             f1 = open(dest + "\\las\\unziplaz.bat", "w")             # now go through the list and download the LAZ files...             lasSuccessCount = 0             for l in lasList:                 q250k = "q" + string.split(l,"-")[0]                         srcLAZ = dataPath + "/" + q250k + "/laz/"+ l +".laz"                 destLAZ = dest + "\\las\\"+ l +".laz"                 destLAS = dest + "\\las\\"+ l + ".las"                 if not os.path.exists(destLAZ):                     arcpy.AddMessage("Copying LAZ " + destLAZ + "\n")                     try:                         urllib.urlretrieve(srcLAZ, destLAZ)                         lasSuccessCount += 1                     except:                         arcpy.AddMessage("\t WARNING: LAZ file for tile " + srcLAZ + " not found - moving on to next one...")                                 if not os.path.exists(destLAS) and os.path.exists(destLAZ):                     f1.write(dest + "\\las\\laszip.exe -v -i " + destLAZ + " -o " + destLAS + "\n")                  #close the batch file                 f1.close()             if lasSuccessCount > 0:                 # add a message telling the user what we are doing...                 arcpy.AddMessage("Unzipping LAZ Files....please be patient")                 # download LASZIP.EXE from the mngeo site...                 lazsrc = "ftp://ftp.lmic.state.mn.us/pub/data/elevation/lidar/tools/lastools/laszip.exe"                 lazdest = dest + "\\las\\laszip.exe"                 urllib.urlretrieve(lazsrc,lazdest)                      #call the batch file to unzip the LAZ files...                 os.system(dest + "\\las\\unziplaz.bat")                      # if the user has ArcGIS 10.1 then build a LAS dataset.                     if arcpy.GetInstallInfo()['Version'] == '10.2.2':                     lasDName = dest + "\\" + "las_data.lasd"                     if not os.path.exists(lasDName):                         arcpy.AddMessage("Creating LAS Dataset - " + lasDName)                         try:                             arcpy.CreateLasDataset_management(dest+"\\las",lasDName,"NO_RECURSION","#","#","COMPUTE_STATS","RELATIVE_PATHS")                         except:                             arcpy.AddWarning("LAS files downloaded but attempt to bring into LAS dataset failed...")                     else:                         arcpy.AddMessage("LAS Dataset - " + lasDName + " found....")             else:                 arcpy.AddWarning("No LAS files found in requested area....")              # if the user wants to merge the downloaded data then do that now....         if merge and gdbSuccessCount > 1:             arcpy.AddMessage("Merging tiles to "+ dest + "\\ELEVATION_DATA.GDB")             merge_tiles(dest)         elif gdbCount == 1:             arcpy.AddMessage("Only one GDB found...no merging necessary")         elif gdbCount == 0:             arcpy.AddWarning("No Geodatabase Tiles found in area requested....")


title: Joining a table with thousands of transactions to one point (location) in Arc Map

text: I have a large transaction data set in the millions and around a 100 positional points. I need to join the transaction data set to these points. Arc map won't let me join them without errors or if it does allow it null values are applied to the join. Any help would be appreciated. They have a common field, a meter number for the machines. I have tried using a simple join with the attribute tables. My end goal is to be able to upload this data to an enterprise and integrate it into an operations dashboard. This dashboard would be able to tell the user the number of transactions with the machine within a time period. The transaction data set also has the date attached to it. ArcGIs-10.0 _update_ I can't comment so this will have to do. I found a problem the one web service map I am using has the meters as strings so in the transactions attribute table I added a new field and made it a text. Using the field calculator I added the meter number I was then finally able to perform the join, but only 153 values are in this join meaning it didn't join everything I wanted. How can I get all of the values to join?


title: only 1 tile visible in GE after gdal_translate -of KMLSUPEROVERLAY

text: I have a 3-band ASTER scene in any of several formats (PCIDISK, GeoTIFF, Imagine), already reprojected to the geographic coordinate system. I can view and manipulate the scene successfully in PCI Geomatica and ESRI ArcMap. I want to create a KMZ version of it, and following in part this other GIS StackExchange post, I tried using:               gdal_translate -of KMLSUPEROVERLAY myscene.tif myscene.kmz      but the result is invisible in Google Earth (although it zooms to the correct area). It remains invisible when zooming way out and way in, except that the lower-left tile of the 16 tiles it created **does** show up when zooming in. It, however, is completely black as it is in the no-data collar around the image data. Also, if I unzip the KMZ and open any of the individual tile KMLs that were created, **they** open just fine, and the ones that are not in the collar show the correct image data. So I suspect there's something wrong in the way the output KMZ was formed (which uses Regions and NetworkLinks, KML elements I have not worked with before). But I don't know whether this is due to input data or format problems, to incorrect or missing parameters to gdal_translate, or something else. I've tried several input formats with the same results each time. Tips and suggestions (if not outright answers!) much appreciated!


Is there any way to select only the layer of which the features have been selected? That is, if I have selected some roads then I want only the road layer to be selected.


title: How to fill Field Mappings parameter in script tool on updateParameters?

text: My goal is to get a list of fields choosen by an user from existing fields in feature class. I think that I can show all fields in Field Mappings dialog, an user deletes unwanted fields and I get a required list.   The problem is how to populate Field Mappings when adding a feature class as a first parameter. I can't select anything in 'Obtained from' property.   The code `self.params[1].parameterDependencies = [0]` in updateParameters function does't work in this case (it worked when 'Field' parameter data type was chosen). What should I use to automatically populate Field Mappings when filling first parameter? * * * EDIT:   Fields list is applied as follows: ![enter image description here](http://i.stack.imgur.com/HGDsE.png)


title: calculate building density in ArcMap

text: I have two questions related to calculate building density.   1. I need to calculate the building density in a map which Data is in shape format. The buildings and the area are in two different layers.   2. Same question as above but the map/area is raster format. Thank you.


title: Land Administration Domain Model example implementations

text: Does anybody know of a sample implementation? I am currently trying to adapt it to my country but the LADM is rather abstract. Anybody have any idea how to go about this?


title: How to perform volume calculation in QGIS

text: EDIT 07/12/12 Ok, meanwhile I came to another problem I want to discuss. I just tried around with the volume calculation with some layers and came to a point where I receive the following error. Screenshot: ![error](http://i.stack.imgur.com/XTQNg.jpg) I did some research and discovered this explanation: http://blog.webforefront.com/archives/2011/02/python_ascii_co.html The problem is, that I have absolutely no Idea about python programming. So my question is: Does anyone know how to solve this described error problem? Note: I left the older questions below in order to see the process I made. * * * * * * * * * I m currently trying to make a volume calculation in QGIS. I realised, that using GRASS or Sextante is the right choice to do this. But since there are barely useful tutorials, that go into depht nor any other guides or hints that may help in this context, I want to seek help here. So I want to ask these questions: 1)I want to calculate the volume of a specific map area. So I discovered the algorithm "r.volume." 1.1.)How does it work exactly?The Help text in that menue is not very useful for me as a Grass beginner. 1.2) Why do I need raster layers for this and no vector layers? 1.3) What certain criteria do these rasters have to fulfill? 1.4) I added a screenshot. What exactly are these "clumps" and what rasters have to be chosen in the drop down list? 1.5) Can I have an accurate example of how to perform such an r.volume algorithm?I d be thankful for any detailed help/ raster files to test this. 2)I learned, that GRASS (in QGIS) was integrated into the Sextante Plugin for QGIS. 2.1) Does this mean I can only use Sextante and still have all the functions of Grassgis?I d welcome that because I always had trouble with importing data into Grassgis. Sextante doesn t (obiously) require this.(?) 2.2) A strengh of **gvSIG** is Sextante, so I read; Is the Sextante Plugin for QGIS equal to gvSIG Sextante concerning performance/ possibilities? thanks a lot, Shepherdjo * * * EDIT:07/11/12 8.08am: I know this is not easy to answer but it would be great too, if you can help me out with the other questions (2.1/2.2). * * * ![r.volume menue](http://i.stack.imgur.com/XUIH7.png)


title: Error when importing shapefile to Postgis database

text: I was trying to import shp to postgis database and i get the following error:               Importing with configuration: adresy92, public, geom, C:\Users\tomek\Desktop\adresy92, mode=c, dump=1, simple=0, geography=0, index=1, shape=1, srid=0     Shapefile type: PointZ     PostGIS type: POINT[4]     Failed SQL begins: "SET CLIENT_ENCODING TO UTF8;     SET STANDARD_CONFORMING_STRINGS TO ON;     BEGIN;     CREATE TABLE "public"."adresy92" (gid serial,     "obrot" int2,     "dane_kod" int2,     "xy_dane" varchar(20),     "status_kod" int2,     "xy_status" varchar(50),     "ulica_kod" int2,     "miejsc_kod" int2"     Failed in pgui_exec(): ERROR:  AddGeometryColumns() - invalid SRID     CONTEXT:  SQL statement "SELECT AddGeometryColumn('', $1 , $2 , $3 , $4 , $5 , $6 )"     PL/pgSQL function "addgeometrycolumn" line 4 at SQL statement     Shapefile import failed.      The stange thing is that I have already put a couple shapefiles to this database without any problems. I am importing shapefile through SPIT tool from Quantum GIS 1.8.0. I have also tried to do the same via shp2pgsql plugin in pgAdmin with the same result. PS I found a similar question on gis.stackexchange but it didn't solve my problem. I would be appreciate for your help.


title: raster algebra in python with rasters of different extents

text: I am trying to find out how to use GDAL & numpy modules to average a set of rasters....which are Sigma0 values from satellite passes at different times within a year. Each raster has been mosaicked from a number of smaller images using GDAL_merge. Because the orbits differ each time they pass over the area of interest the merged datasets are different extents and not regular squares. So I want to average the values from each pass....baring in mind that when theoretically stacked on top of one another, there is sometimes overlaps of one/two/three images. I imagine the best way of me getting around this is making all rasters the same extent ... and using a 'no data' value for areas of the raster where there is no data...and then ignoring these values during the calc of the average... If this indeed the best way, how do I go about making them all the same extent when they are not regular (rectangle)? ...or is there a better way than this? I am new to using GDAL/numpy ....my impression is that calculations using multiple rasters is best done in numpy arrays? again, is this correct. Thank you in advance Becky


title: How to avoid colliding parallel lines in QGIS Composer

text: I am preparing a map which is concerned with showing two line features and there relative location to each other. The line features are 20 m apart and at a scale of 1:50 000 or bigger on A3 I want to show which the lines a little separated from each other. The features follow a common path, but change sides irregularly, i.e. one feature is on the south for while, then the north, then th common path changes north south and the line features then alternate between sitting on the east or west of the path. I would like to do this in QGIS. Is there some sort of automagical dynamic offset control for lines? I can imagine how to rotate symbols based on an adjacent feature, but this is a bit more complicated. I can imagine a method of breaking one line up into segments and then creating an index of those segments compared to the other line acting as an anchor. The index might then be used to specify the offset, however this is a sort one time run solution. Rather just looking for the elegant cartographic solutions, if they exist. ta


title: Get XY coordinates in UTM from web-mercator where UTM zone is not known

text: How to get XY(northing/easting) coordinates in UTM where data is WGS_1984_Web_Mercator_Auxiliary_Sphere(WKID: 3857) coordinate system and UTM zone is unknown. For e.g. India falls under multiple UTM zones and map is in Web_Mercator_Auxiliary_Sphere. I am trying to replicate Google earth's UTM coordinate display using ArcGIS Javascript API. Thanks. Nagesh


title: How to get a collection of georeferencing control points in an ArcMap document?

text: I have created control points (i.e. links) using ArcMap's Georeferencing toolbar. How can I programmatically retrieve a list of these points? I'm looking for something like an enum of the links you would see in the Link Table (which is visible through the Georeferencing toolbar). Preferably, I would then be able to step through these points and delete any one. I am working with the ArcObjects SDK for C# .NET.


title: choosing the right parameters to interpolate a DEM form contours

text: this is my first question on the forums, so be gentle :) I have a 5 meter interval contour lines feature class which i would like to interpolate to a DEM. I am well aware of at least 2 tools ("tin to raster" and "topo to raster") which can help me with that. In order to produce the proper DEM i need to define first its purpose. Long story short - i need to track the progression\advancement of small streams (AKA: gullies). they are typically 1 meter wide and 2 m long. My algorithm which will be using the DEM as a basis, must have as much as accurate DEM as possible. how accurate? - detailed enough to see an advancement of a 30 cm length for example. To do that i must interpolate (from the mentioned above contours) a most pleasing DEM. The questions is: what course of action would you suggest me to implement? in order to choose the correct parameters? I know trial and error would be the most advisable answer, yet there are so many parameters (smoothing, cell size, method type, enforce drainage, toelrance) to consider it might take me weeks until i get the correct DEM. Cheers, David


title: Cutting image using gdalwarp to extent of data

text: I am currently generating a raster from point data and applying inverse distance algorithm. This image is then cut with :               gdalwarp -cutline <<src shape file>> -crop_to_cutline -srcnodata 0 results_in.tiff results_out.tiff      Is there a way to cut the output to the extremities of the original data source file rather than a shape file? Many thanks J EDIT : gdal_grid command gdal_grid -a invdist:power=2:smoothing=0.001:radius1=0.0005:radius2=0.0005 -of GTiff -zField z_value -ot Float64 -l data VRT.vrt EM_RESULTS.tiff


title: Is it possible to add a gvSIG document on Html?And if yes, how?

text: I have completed a project in gvSIG and now I would like to be able to open it from my webpage, like a hyperlink or path or something.Can I really do that somehow? And if yes, could you please explain how? Thanks a lot. Eva


title: Looking for average temperature data for Afghanistan

text: Hey all I am looking for a average temperature data for Afghanistan monthly over a year, precipitation would work too. In a shapefile, if anyone has any or knows anywhere that I could access one please let me know.


title: What format to use for highly compressed vector data?

text: Does anyone have information about possible formats for highly compressed vector data? Some basic requirements.   * Must be possible to read using the larger commercially available platforms. (ArcGIS, Luciad, MapInfo etc.)   * Must have tools that allows us to create files/format from other data formats (i.e. Shape)   * Speed when reading the format is important as the data is among other things used as a live map in mobile platforms. Compression speed is less important (although nice)   * Must be possible to use a map product in a wide scale range, a 50k map shall be possible to read at a 1 000 k zoom level   * Ability to do selective updates on area, feature-ID or other search criteria would be nice, but not a must.   * Memory footprint must be small, it needs to work with Windows XP and/or Windows 7 allocating less than 1GB RAM so some kind of smart cache management must be implemented To give an idea of the required compression and as input for memory footprint requirements, I need to store approximately 1TB of shape data to a 100GB disk. Does anyone have suggestions for good products? Best regards, Johan


title: How to summarize population by neighborhood?

text: I'm a novice with PostGIS and am currently working on a lab assignment where I have two geometric tables - one containing demographic data by blockgroup (blockgrp2010) while the other identifies various neighborhoods (nbo_hood). I'm trying to write a query that lists the population for each neighborhood. For now, I'm assuming that if the centroid of the blockgroup is inside the neighborhood, then all the population lies within that neighborhood. I drafted the query below:               select nbo_hood.name, blockgrp2010.pop10     from blockgrp2010     inner join nbo_hood on ST_Within ((ST_Centroid (blockgrp2010.geom)), nbo_hood.geom)     group by nbo_hood.name, nbo_hood.geom, blockgrp2010.pop10;      This returns a table listing the neighborhoods and corresponding population for each blockgroup. However, many of the neighborhoods are listed under multiple rows with a different population assigned (obviously because there are multiple blockgroups within the neighborhood). My question: how do I modify the above query (or should I start from scratch) to list each neighborhood only once with the population sum of all blockgroups? I've tried variations of the SUM function with the SELECT statement but nothing as worked. Thank you,


title: Is co kriging appropriate for interpolation of three spatial co variates?

text: I want to estimate an interpolation model(s) for daily minimum temperature, daily maximum temperature and daily rainfall. I have samples of approximately the same size for all three variables at different locations although for certain days and locations only 1 or 2 of the 3 covariates are sampled . All three variables have some correlation, especially minimum temperature and maximum temperature. I read that co-kriging is especially adequate in cases where:   * one wishes to interpolate one variable that is relatively sparsely sampled   * one has available another covariate which is more densely sampled This is not my case. Rather I have available samples for three covariates and want to interpolate all three. I have the following questions:   1. Is it better to develop a separate ordinary/universal kriging model for each variable or to develop a co-kriging model for all three.   2. I am using days and locations with no missing values to estimate the parameters of the interpolation model , but want to use it to interpolate missing values in different days and locations. Is it possible to do this in cases where:     * certain locations have all three variables missing     * certain locations have one o two out of the three variables missing I am working in R in case there are any special considerations to heed.


title: How to reset ObjectID after running Delete_Features on SDE feature table?

text: I have a python script that deletes all features on a SDE feature table and then inserts new features to it. I noted that the ObjectID keeps incrementing from where the last feature was with before deleted. Just wondering if there is any reset function or tool to make ObjectID start with 1 every time the script runs? Thanks.


title: OpenLayers.Filter.Logical not working

text: I have the following code:                  var map = {};        var googleMercator = new OpenLayers.Projection("EPSG:900913");        var wgs84 = new OpenLayers.Projection("EPSG:4326");        var lyr_osm = {};        var texas = {};        var lll = {};        var query = {};        var fin = {};        var fout = {};     function init(){         OpenLayers.ProxyHost = "proxy.cgi?url=";     //--------------------------------------------------map         $(document).ready(function () {         map = new OpenLayers.Map('map',                         {                         numZoomLevels: 18,                         projection: googleMercator                         }            );         });     //--------------------------------------------------layers definition         var lyr_osm =new OpenLayers.Layer.OSM("OSM");     //--------------------------------------------------setting wms layer            texas = new OpenLayers.Layer.WMS("Texas",                "http://54.228.242.237:80/geoserver/us_admin/wms",                {                'VERSION' : '1.1.0',                layers: "us_states_gen",                format: "image/png",                transparent: true,                cql_filter: "state_name LIKE 'Texas'",                },                     {isBaseLayer: false,                     visibility: true,                     sphericalMercator: true,                                             }                );     //---------------------------------------------------filters                                  fin = new OpenLayers.Filter.Comparison({                                          type: OpenLayers.Filter.Comparison.GREATER_THAN,                                          property: "store_id",                                          value: 500                                  });                                  fout = new OpenLayers.Filter.Comparison({                                          type: OpenLayers.Filter.Comparison.LESS_THAN,                                          property: "store_id",                                          value: 100                                  });                                  query = new OpenLayers.Filter.Logical({                                         type: OpenLayers.Filter.Logical.AND,                                         filters: [fin, fout]                                  });              var awfspc = new OpenLayers.Protocol.WFS({             version: "1.1.0",             url: 'http://loclahost/geoserver/stores/ows',             featurePrefix: 'stores',             featureType: 'stores',                 defaultFilter: fin,                 srsName: "EPSG:3857"                 });              lll = new OpenLayers.Layer.Vector("Points", {                  externalGraphic: 'http://localhost/black_marker32.png',                   graphicHeight: 20,                   graphicWidth: 20,               strategies: [new OpenLayers.Strategy.Refresh(), new OpenLayers.Strategy.Fixed()],               protocol : awfspc,         });          //--------------------------------------------------------------------map controls         map.addControl(new OpenLayers.Control.LayerSwitcher());         map.addControl(new OpenLayers.Control.MousePosition());     //--------------------------------------------------------------------add Layers         map.addLayers([lyr_osm, texas, lll]);     //--------------------------------------------------------------------map position         var mapCenter = new OpenLayers.LonLat(-99, 31.3);         mapCenter.transform(wgs84, googleMercator);         map.setCenter(mapCenter,6);     //  map.zoomToExtent();     }      When I use fin filter or fout filter, the filters are applied and they return good. When I use querry filter, nothing is returned, no error and no data. Any clue? Thanks!


title: Maplex: try different fields

text: Let's say, a road has the name "William-Bill-Gates-Road". If the road is too short, the label can only be abbreviated by using a dictionary in Maplex. Is it possible to make Maplex try several label fields, if one is not fitting? Structure: name | name_1 | name_2 | ... William-Bill-Gates-Road | William-Bill-Gates-Rd. | W.-Bill-Gates-Rd. | W.-B.-Gates-Rd.


title: Ecuador Census Data

text: I am looking for Ecuador census data (digital GIS if possible) at the sub- parroquia (parish) level. The geographies and the tabular data would be great but at this point I am mostly in need of the digital boundaries. I am looking for the full demographic suite of data with attributes if possible. 2010 would be great but I am also interested in 2001 and 1990. Any digital GIS format will suffice for the boundaries and any tabular format for the tabular data. Non-free options are fine as well. The geographies are known as Sector and Zona (Zone) in the Ecuadorian Census.


title: Calculate end point Height above the MSL of the target point given a starting point, distance, azimuth and elevation

text: We need to find the end point height (from the earth surface) given the start point lat, lon, azimuth (270 deg), distance (500 km) and elevation above the horizon of start point (10 deg). The question **How To Calculate Height of Point B?** does not seem to consider the earth's curvature:   * Here the angle (elevation) at point a is 10 deg    * so the calculation referred in the above question is only usual cartesian calculation (i.e h).  I would like to know the mathematical formula for the geodetic calculation of the height POINT B above the MSL at POINT B (i.e h +H) in the picture shown below. Also I would like to program the formula using VB.NET on a independent GIS application. Accuracy in the order of meter would be good enough. Please provide the mathematical formula to calculate the height above earth surface at the end point? ![enter image description here](http://i.stack.imgur.com/2X7kO.jpg)


title: How to debug QGIS installation problems on XP?

text: I'm installing Qgis (Standalone 1.8) on all the PCs of the town planning office where I work. All good for 4 computers, on the fifth installation does not work. Computers have all XP with SP3. In particular: the installation seems to proceed normally until the end, then unlike the other pc, ends without the window appears as Microsoft Visual C++ and bypasses all the post-install. After installation, the desktop icons are missing Grass, MSYS and OSGeo4W. There are Qgis Qgis desktop and browser, but will not start. The same problem occurs with the previous version. I can not figure out what is missing on this computer, which is mine. Is there anyone who has had the same problem and fixed it?


title: How To Do OpenLayers.Mobile Popup

text: It seems that the OpenLayers.mobile.js library (v 2.12) doesn't have a popup class. Does anyone know how to do a popup with the mobile version? I'm using some code that i ported over form normal OL web map. When i call new OpenLayers.Popup.Anchored i get an undefined error on OpenLayers.Popup class. It's not in the OpenLayers.mobile.js. I'm lost


title: Creating a custom, local Route Task

text: I'm trying to replicate something similar to the ArcGIS WPF Runtime SDK (C:\ArcGIS\WPF1.0\SDK\Samples\Sample Application\Sample Application.exe) "Routing with Barriers Online" sample with two major modifications ... 1. I have my own network dataset (created from a feature class using ArcObjects classes) 2. I want it to run locally and not on an ArcGIS server instance. Can I exploit LocalServer for this? First, can this be done? If so, what is the best approach? I have both the WPF runtime and ArcObjects SDK to work with.


title: merge shapefile using ogr2ogr and retrieve the attribute of both

text: Ciao, I have two shape files that I want merge (overlay/intersect) having in the destination the attributes of the two sources file.               ogr2ogr file_merged.shp tl_2010_us_uac10_conus_albersUA.shp     ogr2ogr -update -append file_merged.shp tabblock2010_10_pophu.shp -nln file_merged      In this way I obtain the polygons of the two shapes, but I lose the attributes of the `tabblock2010_10_pophu.shp`. My final aim is then make an sql statement where I can retrieve polygon in common in the both shapes. Any suggestions?? Thanks in advance


title: ArcGIS 10: convert labels to points

text: is there a way to convert labels, created by Maplex labeling engine to point features? In MapInfo and CAD software, it's possible to convert labels to points. Those point features include a text and a rotation for this point. When text is horizontal, it includes all the text. When the text is curved, it included one char for each point with the rotation in degrees. That makes it possible to use the labels in other software or store them in a database. ArcGIS 10 only allows to convert labels to annotations, which result in polygons... Thanks for help :)


title: Interacting with the Map after selecting a button/tool

text: I am using the ESRI Javascript API 2.0 and Dojo dijits. Does anyone know how to wire up a dojo button (or toggle button) in a toolbar to a mouse click event on the map? I want the user to have to select a tool/button before clicking on the map and NOT just have the map listen to any old mouse click event as this is too open i.e. I want to constrain the functionality to the selecting a tool from the UI. Many thanks


title: Detect whether there is a spatial polygon in a spatial extent

text: Given a spatial extent (for example `cbind(c(0,0),c(10,10)`),I want to see whether there is a part of a SpatialPolygons falls in this extent? How can I do this? Thank you very much!


title: Python 3 with QGIS on Mac?

text: I need to start learning Python, and I want to learn with version 3.x. I understand that QGIS--at least some plug-ins--relies on Python and Tcl/Tk, and probably on specific versions. Assuming that QGIS is dependent on Python 2, can anyone tell me how to set up parallel installations of Python 3 and Tcl/Tk such that QGIS will have what it needs and I can use Python 3 for other purposes?


title: Converting Coverage files into shapefiles?

text: I am attempting to convert coverage files into shapefiles using arcpy. This is not the e00 files, but the bin files. ArcGIS help suggests that it can be done using **FeatureClassToFeatureClass_conversion** , but I am getting an error saying "parcels does not exist or not supported".               fileName = 'C:\\May_2014_Parcels' + "\\" + filename + "\\" + subfilename + "\\parcels"                        arcpy.FeatureClassToFeatureClass_conversion(fileName,            r'C:\May_2014_Parcels\Output\parcels.gdb', filename + "_" + subfilename)      How do I actually get the featureclass out of the coverage file so I can perform this task?


title: Show M values on mappoint /mouseover on a polyline - ArcGIS Server 10

text: I want to show the M values i.e the chainage of the road when I move my mouse over it. I am using mappoint to show some static attributes but I also want to show the dynamic M values along the road when I mouse my mouse pointer over it.. Note: I am using Visual Studio 2010 with Silverlight API 2.2 and ArcGIS Silverlight 10 Please guide...


title: how to filter point features based on the filterted polygons with OpenLayers

text: I have one WFS polygon layer and one WFS point layer. The polygon layer has an attribute user_id. The web map has two requirements.   1. Filter the polygons based on the attribute user_id. That is, a user can only see polygons which are assigned to him. This can be done like this: filter: new OpenLayers.Filter.Comparison({ type : OpenLayers.Filter.Comparison.EQUAL_TO, property : "user_id", value : user_id })   2. Filter the points based on the filtered polygons. That is, a user can only see points within the filtered polygons from step1. I checked the  OpenLayers.Filter.Spatial.INTERSECTS. it only takes OpenLayers.Bounds || OpenLayers.Geometry. I have no idea how to achieve this. Please help. Thanks.


title: How to return the first row of feature class as a list in ArcPy?

text: I need to return the first record of a feature class as a list. Is this possible? I can't seem to figure it out using `row.getValue()` as it does not allow me to put a wildcard in there.


title: How to get route plan with certain distance with pgRouting?

text: I am working for a project. I have a road network which has been imported in PostgreSQL/PostGIS database. I can query shortest path between two points. I want to get route plan with certain distance. Suppose, I want to go from point A to point B. There may have different way to go from point A to B. But, I want to fix the length, suppose, 10 km. That means, I want to know the route plan to go from point A to point B but in 10 km range. Is it possible to get this kind of query using pgRouting? Or I should use any other tools? If it is possible how should I write the query? Note: in my database table 'route', i have 'gid', 'source', 'target', 'cost' etc columns. Thanks, in advance.


title: Vehicle Routing Problem - Cost Calculation According to the furthest point in route

text: Certain logistics companies do not calculate their costs according to the total distance travelled but according to the furthest point in the route. In our problem, Vehicle Routing Problem do not really help as almost all of the cases, furthest point is always have less travel distance than the whole route. Do you have any idea on how to tweak Network Analyst in ArcGIS, or know any open-source solution that I can solve this problem with? Thanks a lot in advance, Anil. EDIT:![enter image description here](http://i.stack.imgur.com/mf8rc.png) As seen on the graphic above, it really does not matter how many points you stopped between start and end point. However, we still would have other constraints, such as travel time, vehicle capacity, multiple routes and depots.


title: Field Collection Apps for Windows 8

text: I was wondering if any of you out there have experience with any apps that work on windows 8 tablets free or paid?We have a client that wants us to use their tablets for data collection. i was wondering if some of you out there have used apps for the windows store that you can recommend. Thank you


title: QGIS style save issues

text: I have a map that I want to do graduated chloropleth maps for each year over 10 years, and I would like to use the same graduations for each year. I know this was asked earlier and therefore I know that I save it as a style and then load it for each new year, but I'm having problems with the save style/load style function in qgis - every time I load the style the data converts back to the original year. Has anyone else had this problem and figured out a solution? I'm using qgis 2.2.0 valmiera (desktop).


title: Geoserver shutdown after load big imagemosaic

text: catalina.log: 24.05.2013 15:33:38 org.apache.tomcat.util.net.JIoEndpoint processSocket SEVERE: Error allocating socket processor java.lang.OutOfMemoryError: unable to create new native thread at java.lang.Thread.start0(Native Method) at java.lang.Thread.start(Thread.java:657) at org.apache.tomcat.util.net.JIoEndpoint$Worker.start(JIoEndpoint.java:513) at org.apache.tomcat.util.net.JIoEndpoint.newWorkerThread(JIoEndpoint.java:744) at org.apache.tomcat.util.net.JIoEndpoint.createWorkerThread(JIoEndpoint.java:723) at org.apache.tomcat.util.net.JIoEndpoint.getWorkerThread(JIoEndpoint.java:757) at org.apache.tomcat.util.net.JIoEndpoint.processSocket(JIoEndpoint.java:789) at org.apache.tomcat.util.net.JIoEndpoint$Acceptor.run(JIoEndpoint.java:355) at java.lang.Thread.run(Thread.java:679) ...................... and Caused by: java.lang.NoClassDefFoundError: Could not initialize class sun.font.FontManager at sun.java2d.SunGraphicsEnvironment$2.run(SunGraphicsEnvironment.java:179) at java.security.AccessController.doPrivileged(Native Method) at sun.java2d.SunGraphicsEnvironment.(SunGraphicsEnvironment.java:163) at sun.awt.X11GraphicsEnvironment.(X11GraphicsEnvironment.java:253) at sun.reflect.GeneratedConstructorAccessor109.newInstance(Unknown Source) at sun.reflect.DelegatingConstructorAccessorImpl.newInstance(DelegatingConstructorAccessorImpl.java:45) at java.lang.reflect.Constructor.newInstance(Constructor.java:532) at java.lang.Class.newInstance0(Class.java:374) at java.lang.Class.newInstance(Class.java:327) at java.awt.GraphicsEnvironment.getLocalGraphicsEnvironment(GraphicsEnvironment.java:82) at java.awt.image.BufferedImage.createGraphics(BufferedImage.java:1152) at javax.media.jai.TiledImageGraphics.getBogusGraphics2D(TiledImageGraphics.java:266) at javax.media.jai.TiledImageGraphics.(TiledImageGraphics.java:192) at javax.media.jai.TiledImage.createGraphics(TiledImage.java:755) at javax.media.jai.ROIShape.getAsImage(ROIShape.java:1148) at com.sun.media.jai.opimage.MosaicOpImage.(MosaicOpImage.java:374) at com.sun.media.jai.mlib.MlibMosaicOpImage.(MlibMosaicOpImage.java:49) at com.sun.media.jai.mlib.MlibMosaicRIF.create(MlibMosaicRIF.java:74) ................... and more ERROR [geotools.map] - Call MapContent dispose() to prevent memory leaks .............. After run layer preview(big imagemosaic) geoserver(and other opengeo) send error 502.


title: How to scale/reduce my polygon without changing the central lat & long?

text: I am writing a VBA code to resize my polygon, I get a input as vector points(Shape file). The code below resizes it as according to my scale factor but my central latitude and longitude are also distorted. I want my new polygon to have the same central lat and long as of the orignal shape file. As if its overlaid.               lVct = LBound(vct)     uVct = UBound(vct)          ScaleFactor = 0.59       j = 0     i = 1     Set PPolygon = New Polygon      For i = lVct To uVct               vct(i).y = vct(i).y * ScaleFactor               vct(i).x = vct(i).x * ScaleFactor               PPolygon.addPairOfCoordinates vct(i).y, vct(i).x   'add coordiantes into the new Polygons            Next i      Describe if I have to use to some coordinate reference system or projections !!!


title: Which geodatabase type is most appropriate for both open source and ESRI use?

text: Main Question: Is there a geodatabase type that seamlessly works with both ESRI and open source products? Details: My company currently passes around geodata in various vector and raster file formats and isn't particularly organized. Most processing is done using open source tools, but some customers demand that we host our data in ESRI Map service format. We have an ArcGIS server instance already that hosts these services, but we still do our processing upstream using open source tools before dumping the results in an ESRI arcsde geodatabase. From what I can tell, ogr/gdal have somewhat limited abilities to read/write from an arcsde geodatabase, but not enough that I could reasonably store all of our data in an arcsde geodatabase and expect our processing tools to be able to use the arcsde geodatabase as the primary data store for the company. Am I wrong? If I switched to an open source geodatabase such as postgis, our processing tools would work fine but I don't think that I could use a postgis database as the data store for hosting ESRI services. Am I wrong? Is there a geodatabase type out there that can service both open source and ESRI seamlessly? If not is it common practice to maintain parallel geodatabases, one for open source and one for ESRI? Additional wrinkle: The answer I got was to use PostgreSQL/PostGIS . The only potential wrinkle with that is that I'd like to host this database on Amazon Web Services (AWS) and if possible I'd like to use Amazon Relational Database Service (RDS) to minimize our database management overhead (backups, load balancing, etc). However, RDS only currently supports MySQL, Oracle and SQL Server as database engines. I'm thinking that if RDS starts supporting PostgreSQL that this would be the ideal scenario, but just to double check, PostgreSQL is my only option for ArcGIS + gdal/ogr-based open source interoperability correct? If so I'll likely create a PostgreSQL database on AWS, I'll just have to spend more resources managing it.


title: Use FME to delete records from PostGIS before loading new records

text: I have some shapefiles. Each shapefile contains the data on a topic for a county. I have data for 26 counties. In the filename of the shapefile, the county ID is given. I load these shapefiles into PostGIS using FME. The data for all counties is loaded into one table (so 26 shapefiles get loaded into 1 PostGIS table). This data gets updated every now and then, and I can request the data for the county where the data has been changed. I get supplied with a new shapefile with all data for that county (not just the changed records). I would like to delete the data in the PostGIS table for the county that I have the new data for before I load the new data. How do I do this? Is this possible if I get the county ID from the filename of the shapefile and each record in the PostGIS table has a county ID attached to it? My plan is to just add the new shapefile to a specified folder, and run my FME workspace. Edit 31/05/2013 Just to clarify, I have multiple records for each county (i.e. all roads, buildings etc in a county). I basically want to bulk delete a counties data, and replace it with the new data.


title: POSTGIS query of point feature.  How to find nearest point to query location?

text: I am querying a point feature using a POSTGIS request fired in Javascript. The decoded URL looks like this and returns great results for polygon features.               http://test.tablelocation.com/api/v1/sql/?q=select '+tempCoi+' from '+tempTableName+' where ST_Intersects(the_geom, ST_GeomFromText('POINT('+clickLng+' '+clickLat+')',4326))&format=json&diagnostics=true&_maxage=86400&callback=?      As you can see, there are a few dynamic Javascript variables that are injected into each URL request. Two of these variables (clickLng and clickLat) specify the lat/lng of my query. They are generated based on where the user clicks on a google map. I am returning one attribute from the point layer referenced via the tempCoi variable. When I try to query a point feature, no data is returned. I am guessing that this is because the location of my query is too specific (i.e. 6+ decimal points). To make a long question shorter, is there a way to fuzz my data request. For example, if I click at 45.123456N and -113.654987W how could I look for the closest point to that query? Hope this is clear.. thanks in advance.. -j


title: How to determine if a LYR file has more than one Layer using ArcObjects?

text: In using arcpy, I have found that there are just a few lines of code needed to get a list of layers in a LYR file if more than one layer exist. How can I accomplish this same feat using ArcObjects in C#? I know that ILayerFile is the way to start this off and I can open a specific layer file using ILayerFile.Open() but how can I get back a list of layers in this layer file if more than one exist?


title: How can I extract raster values around points?

text: I'm using ArcMap 10. I have a point dataset and a .flt raster. For each point I want the sum of the raster values that lie within a certain radius of each point. i.e. I want to extract values within a given radius around each point and then sum them for each point. I just can't work out how to do it! Any help would be really appreciated... Thanks!


title: View Markers on Both Sides of OpenLayers Dateline -- Tips and Tricks

text: I'm looking for advice on how to display vector markers on both sides of the dateline. I'm using a GeoRSS feed. Currently only one OpenLayer's "globe" worth of data shows at a time. You have to pan east/west of dateline to get that "globe's" data. I've read about some tricky hacks that might work -- though they'd be expensive to do depending on the size of your dataset. So i'm wondering what other hacks people here have come up with. Please share your ideas.


title: Openlayers popup not closing

text:
This is a special case, I have many popups attached to objects and they work fine and I am having a problem in only one case. All layers and popups are handled by select control's `onFeatureSelect`, but this one is handled by layers `featureSelected` event. On opening popup it throws: > Uncaught TypeError: Cannot read property 'lat' of null. Expanding error shows error on function `OpenLayers.Bounds.OpenLayers.Class.determineQuadrant` but the popup opens. On Clicking the close button, the click event is also fired and `onPopupClose` function is also called ( I have checked it with monitor and monitorEvents in console) but the popup does not close. The popup in other layers is still function normally. Here is the code:               select_Control = new OpenLayers.Control.SelectFeature([business, schools, healthfacilities, open_space], {         onSelect: onFeatureSelect,         onUnselect: onFeatureUnselect,         click: true     });          open_space.events.on({         "featureselected": function(e) {             var PopupPos1 = e.feature.geometry.getBounds().getCenterLonLat();             popup_content = '

sifal ground

';             var popup1 = new OpenLayers.Popup.FramedCloud("chicken",                 PopupPos1,                 null,                 popup_content,                 null, true, onPopupClose);                  e.feature.popup = popup1;             map.addPopup(popup1);             popup1.setBackgroundColor('#FFBBBB');             popup1.updateSize();             popup1.draw();         }     })          function onPopupClose(e) {         select_Control.unselectAll();     }          function onFeatureSelect(feature) {         //debugger;         var PopupPos1 = new OpenLayers.LonLat(feature.geometry.x, feature.geometry.y)         popup_content = feature.attributes.name;         var popup1 = new OpenLayers.Popup.FramedCloud("chicken",             PopupPos1,             null,             popup_content,             null, true, onPopupClose);         feature.popup = popup1;         map.addPopup(popup1);         popup1.setBackgroundColor('#FFBBBB');         popup1.draw();     }          function onFeatureUnselect(feature) {         map.removePopup(feature.popup);         select_Control.unselectAll();     }


title: GCM selection for future climates

text: This is not a coding question but rather an inquiry of which Global Climate Model(GCM) to use when projecting species distribution models for future climates. With the release of the IPCC AR5, and the multiple options of GCMs WorldClim GCM options, I dont know which one would be most appropriate for my work in the Australian tropics and in the North American deserts. Does anyone have any literature suggestions for selecting the appropriate GCMs?


title: I can not download gpx from GPS through GPS tool with GPS babel on QGIS

text: I have already install GPS tool plug-in and GPS Babel. Garmin GPS map 62s is connected my PC through USB cable. I changed device setting according to below site. http://nickmcw.wordpress.com/2010/07/30/quantum-gis-using-a-usb-gps-device/ I set below condition on "Download from GPS"tab. GPS device: Garmin USB( ←I create by myself) Port:usb Feature type:way point Layer name:F:\Garmin\GPX output file:test At least, I got new GPX file on layer of QGIS. Altgough the new GPX file is expected include existing gpx information from GPS, there is no row in attribute table. Colum names were put automatically.followings are. Name/elevation/symbol/comment/description/source・・・・・ Why I can not get gpx file from GPS???


title: ArcMap trace the creation/moving/deletion events by windows user

text: I want to check whether there is a way to trace the creation, moving and deletion in the editing session by windows currently logged in user. What I only know now is how to trace the creation/moving/deletion events by database user and displayed in attribute table of the layers under the monitoring. I am using ArcMap 10.2.2.


title: Simplest Way to Create a Custom Choropleth Map for a Website?

text: I am working on a choropleth map for a real estate site. My current workflow is:   1. Edit and save SHP file's DBF in Libre Office Calc.   2. Load SHP file into Quantum GIS. Set colored breakpoints based on data.    3. Load modified SHP file in TileMill (MapBox). Layer it over basic map.   4. Push the map to MapBox host.  This whole process is very tedious considering that the choropleth map requires frequent updating. Now, I am trying to find a more efficient workflow to produce the same result. My intuition is that leveraging Google Maps and Google Spreadsheets may be a better way to go. However, I haven't found a good tutorial that accomplishes my specific needs. Any advice is greatly appreciated.


title: What's the best projection for rasterization of random latitude and longitude data in the northern Atlantic?

text: I have a dataset in the WGS84 geographic coordinate system that I would like to interpolate along a grid using Python. It appears `gdal_grid` is not available through the Python bindings, so I plan to use another 2D interpolation method after rasterizing the data. To rasterize the data, I plan to re-project it using pyproj, then write it to a raster using the gdal python bindings. The UTM projection coordinate system seems like it would be the best projection, but my data spans the North Atlantic (i.e. multiple UTM zones). **Update:**   I'm open to suggestion on the interpolation method, but I had planned to use an inverse distance interpolation (as shown here). These are presence absence sightings data from ship cruise tracks, which do not necessarily follow lines of latitude. They span from latitudes in the mid 70's to mid 50's. Thanks! ## What projection would be best to use?


title: How to get the top left and lower right coordinates of polygon in ArcObjects?

text: First of all, draw the rectangle               ESRI.ArcGIS.Display.IScreenDisplay screenDisplay = Map1.ActiveView.ScreenDisplay;     screenDisplay.StartDrawing(screenDisplay.hDC, (System.Int16)         ESRI.ArcGIS.Display.esriScreenCache.esriNoScreenCache);     ESRI.ArcGIS.Display.IRgbColor rgbColor = new ESRI.ArcGIS.Display.RgbColor();     rgbColor.Red = 255;          ESRI.ArcGIS.Display.IColor color = rgbColor;     ESRI.ArcGIS.Display.ISimpleFillSymbol simpleFillSymbol = new         ESRI.ArcGIS.Display.SimpleFillSymbol();     simpleFillSymbol.Color = color;          ESRI.ArcGIS.Display.ISymbol symbol = simpleFillSymbol as         ESRI.ArcGIS.Display.ISymbol;     ESRI.ArcGIS.Display.IRubberBand rubberBand = new         ESRI.ArcGIS.Display.RubberEnvelope();     ESRI.ArcGIS.Geometry.IGeometry geometry = rubberBand.TrackNew(screenDisplay,         symbol);     screenDisplay.SetSymbol(symbol);     screenDisplay.DrawRectangle(geometry as ESRI.ArcGIS.Geometry.IEnvelope);     screenDisplay.FinishDrawing();      How to get the coordinates (top left, lower right) of this rectangle?


title: Is it possible to display an OpenLayers map upside-down?

text: This is more fun than a serious question, but I'd like to know if it's possible to display an OpenLayers map upside-down as per the uncommon, but interesting upside-down maps. I guess it's mostly a matter of finding a source of the tiles and turning them 180 degrees, but I figure it's something someone has already thought through. I'm not the only one who's interested in the stories that are (or are not) being told through maps.


title: How to get latitude /longitude in Degree Minute Second (DMS) using QGIS field calculator?

text: I have come up with the following expression to convert 5.1234 to 5° 7`24.24" **toint ($x) || '° ' || toint((($x) - toint ($x)) * 60) ||'`' || substr( (tostring(((($x) - toint ($x)) * 60) - toint((($x) - toint ($x)) * 60)) * 60),1,5) || '"'** Here’s the problem : for some points such as 5.1234, it works. but for others, it doesn’t. I suspect the problem is the integer conversion which round up the decimal figures instead of truncating it. Is there any other option? thanks.


title: Openlayers zoomToExtent pause until image is retrieved

text: Problem:   1. layer transition effect: resize   2. Do programmatic map.zoomToExtent() to somewhere OUTSIDE of current wms map bounds   3. OpenLayers setCenters to target bounds and zooms in.   4. While new wms image is being downloaded, OpenLayers displays 'enlarged' map   5. Once wms image is loaded, it's transitioned into map viewport Is there a way to 'prevent' OpenLayers from set-centering at #3 until wms image is downloaded ? What I want:   1. as above   2. as above   3. as above   4. OpenLayers doesn't setCenter nor zoom-in (or just doesn't update map)   5. if new WMS image is loaded, then update map.


title: Can a spatial View in ArcSDE be viewed in ArcCatalog?

text: Can a spatial view in ArcSDE created by "sdetable create view" command be viewed as feature class in ArcCatalog, or is it just another plain table? I have a feature class table and an attribute table (non spatial) stored in ArcSDE database. They can be related by a command column. I plan to create a spatial view in ArcSDE joining those tables and also publish the view as a feature service in ArcGIS Server. Is it possible? And, if the view can be published as a feature service, can feature operations such as Add, Edit and Delete Features be applied to all the view fields? Thanks for any idea! Thanks Alex


title: Holux GPS logger .itm fileformat

text: I'm looking for a spec sheet or a library to open and decompress .itm files from the Holux GPS logger. Has anyone worked with the proprietary .itm format already and can provide me with some usefull information? Thanks! Martin


I'm new in pgrouting. Is it possible to get 10 (some) ordered shortest paths ordered by distance (cost) with PgRouting (ex: SELECT TOP 10 SHORTEST)? Update: why "second shortest" needed? "second shortest" or "10 shortest path", this is the because the shortest path is not the best, for example if all people go to work with the shortest path, it will become a traffic jam, so "the shorted path is not the best". This is can be applied too with "10 cheapest path". Is this impossible? -aris setyawan


title: How to achieve Activity Sizing in QGIS

text: I have plotted 30,000+ points in QGIS by postcode, various data items for each. How do I resize these plots based on the activity at each of these points i.e. small dot for <100, larger dot for 100-200 and larger dot for 200+?


title: Version editing for Postgresql 9.2.4 / postGIS 2

text: **Background:** I am "building" a GIS with PostgreSQL 9.2.4 / PostGIS 2 as the database and QGIS 1.8.0 as the geographic application. There will be multiple users accessing, editing and viewing the data at the same time. **Question:**   1. Is version editing automatically enforced in PostgreSQL 9.2.4 / PostGIS 2?   2. Is it possible to enable/enforce version editing and if yes how? I have had a look at pgVersion for PostgreSQL 9.2.4 / PostGIS 2 and the QGIS pgVersion plugin. However, if I understand correctly (reading the help files and testing it), pgVersion only allows editing of cell values of existing columns and also no new geographic features can be added. Thanks


title: How to display pop up image when click point in open layer web map

text: I have basic WMS web map showing parcels and points layer which representing buildings of parcel.Then I want to add building images for each and every point. Then when a user clicks on a point image should appear as a pop up window. I am using OpenLayers 2.10 Beginner's Guide book and can't find any method to do the above thing. Is this possible with open layers? If so please give me some advise or some links to help me. [All layer are in WMS form and i am using geoserver and tomcat.] Thanks lot for your valuable time.


title: Creating shape file from file with coordinates

text: How can I make a shape file of an csv file with coordinates?


title: GeoServer WPS Aggregation Functions Example

text: The scenario: I want to present a user with names of polygons from a layer in GeoServer. I then display the polygon(s) whose name(s) the user selected. What I've done: I have a search box, which allows a user to enter name of a polygon they'd like to see. This works perfectly fine, even for wildcards. I've implemented it using CQL filters. What I want: Instead of a user entering a name of a polygon, I want to pull all the polygon names and display them to the user for selection. After some research, I found out that my best bet is to use GeoServer's WPS aggregation function `gs:Aggregate`. I've been looking for examples on how to use this, but no success so far. Can anyone who's had this experience point me in the right direction. I'll also be grateful for any alternative approaches. My environment: GeoServer 2.1.3, PostGIS 1.5, OpenLayers 2.11.


title: What are some good introductory books or articles about Open Source GIS for students coming from an ESRI background?

text: I am preparing an intermediate GIS course for graduate and undergraduate students who have likely not been in touch with anything but ESRI software. It is an existing course in the curriculum that I will be teaching for the first time. Currently, ArcGIS is the weapon of choice for the lab sections and practical assignments. I want to tweak the course a little to include an introduction to Open Source GIS alternatives. For now, this part of the course will only be two to four weeks (I'm thinking a kind of extended epilogue) so I won't be able to dig too deep. I hope to branch this off into a full Open Source GIS course next year, but curriculum constraints prohibit me from doing that right away. Here's some reading I have ben considering to support the Open Source GIS part of the course, to give you an idea of what I'm (not) looking for:   * The Geospatial Desktop is a book I'd love to use for a full Open Source GIS course but is too much to cover in a few weeks. The chapter 'Survey Of Desktop Mapping Software' looks like something I could use.   * The Dekstop GIS book is of similar breadth but currently out of print.   * The Grass Book seems too focused on GRASS. I think GRASS will put students with an ArcGIS mindset off. Also, I want them to learn to appreciate the breadth of the OS geospatial software spectrum.   * 'How to go from GIS novice to Pro without spending a Dime' takes a good, practical approach to delving into OS GIS, and has good links for further reading.   * A white paper from OpenGeo talking about markets for geospatial software, and how that landscape is changing. **Can you suggest other articles and / or books that would be useful to ease students coming from an ESRI / ArcGIS background into appreciating and using Open Source alternatives?**


title: Where can I download a free dataset containing major ports across the world?

text: I am working on an analysis and would like to incorporate major maritime ports from across the world. Some initial searching turned up a dataset produced by General Dynamics, however it will be prohibitively expensive. I would like to find a free dataset to use, preferably in shapefile or some other Arc friendly format. The only information desired beyond the location (point or polygon) would be the name, however even that is not required. Any help is much appreciated.


title: ArcObjects: Setting image of a command button to a partially transparent bitmap

text: In an ArcMap 10 add-in, I am attempting to create a toggle button of sorts whose image changes based on the button's state. This answer describes how to set an add-in button's image programmatically and I am using the same technique. The problem is that this method does not appear to support full alpha channel transparency, only single-color transparency, so those fancy PNG icons they added at ArcGIS 10 don't look quite right. The button's image is set initially by the add-in framework (using a 32-bit PNG file with alpha channel) -- the image looks as you would expect with proper partial transparency. To create the image toggling behavior, in the button's `OnClick` method, I check the state of a boolean variable and set the button's image (using the technique in the above-mentioned answer) to one of two bitmaps stored as resources in the add-in assembly. This all works fine except that the button's image is given a gray background instead of a transparent background. I have searched around a bit and come to the conclusion that ESRI's `GetIPictureDispFromBitmap` method converts from a GDI+ bitmap to a GDI bitmap and in doing so loses the alpha channel of the original bitmap. I have found a solution to this issue and will post it as an answer shortly.


title: How do I implement 'weighted' kriging in R?

text: I'm trying to interpolate the values of an attribute (called "disease_rate") that i've sampled at a set of spatial points and i'd like to add 'weights' to the resulting kriging estimates. Basically, I have about 1500 points (a random sample of 50 is shown below). I have an attribute vector (the rate of a particular disease) and a vector of absolute population sizes ("population"). I'd like to interpolate about 10,000 new points for "disease_rate", and at the same time weight the estimates by absolute population size (as areas with higher absolute populations have absolutely more people with the disease, and should probably have greater influence on the kriging estimates). I'm using the packages "automap" and "gstat", with the following code:               library(automap)     vgm <- autofitVariogram(disease_rate ~ 1, , input_data = dat)     plot(vgm)          kr.ok <- autoKrige(disease_rate ~ 1, input_data = dat, model = "Exp",         verbose = TRUE, maxdist = Inf)      plot(kr.ok)      I'm stuck on two points: 1) how to add the weights? Is there an argument for this within the autoKrige() or krige() functions? I don't see a way to do it. 2) when using autokrige() from the package "automap" a set of 5000 point locations for interpolated values is generated by default, by using I believe a convex hull. In "gstat", how do I get a similar sample of point locations? At the moment I'm doing this:               prediction_locations <- spsample(SpatialPolygonsDataFrame, n = 10000,         type = "regular")     prediction_locations@coords <- jitter(prediction_locations@coords)      the object "SpatialPolygonsDataFrame" is too big to include here, but provides a reference area within which to sample the points. However, I always get the following error:               "chfactor.c", line 130: singular matrix in function LDLfactor()     Error in predict.gstat(g, newdata = newdata, block = block, nsim = nsim,  :        LDLfactor      Any help/suggestions for the above two problems is much appreciated.               # sample of data (50 points)     dat <- new("SpatialPointsDataFrame"     , data = structure(list(population = structure(c(13L, 37L, 9L, 25L, 5L,      2L, 40L, 31L, 10L, 43L, 22L, 7L, 32L, 19L, 12L, 28L, 23L, 30L,      44L, 6L, 39L, 1L, 4L, 35L, 15L, 16L, 20L, 21L, 27L, 24L, 18L,      41L, 42L, 26L, 14L, 45L, 47L, 33L, 36L, 3L, 46L, 4L, 4L, 11L,      29L, 38L, 48L, 34L, 8L, 17L), .Label = c("10323", "1070", "1142",      "1218", "1226", "1228", "12858", "1583", "1619", "1645", "1737",      "1832", "1944", "2003", "2125", "2140", "2248", "273", "2816",      "2868", "2986", "3093", "322", "3338", "3707", "371", "3841",      "4151", "4184", "4189", "4198", "434", "4528", "491", "546",      "560", "65", "654", "7138", "735", "776", "82", "845", "879",      "91027", "943", "948", "989"), class = "factor"), disease_rate = c(0,      0, 154.4163064, 140.2751551, 228.3849918, 0, 0, 71.46260124,      0, 0, 226.3174911, 107.5091726, 0, 159.8011364, 0, 120.4529029,      0, 45.35688709, 227.5312856, 0, 77.59122357, 33.5322916, 0, 0,      0, 93.45794393, 239.0914525, 0, 0, 0, 0, 0, 0, 0, 0, 189.7113215,      0, 0, 0, 175.1313485, 0, 369.4581281, 0, 0, 103.5691523, 0, 0,      0, 0, 111.2099644)), .Names = c("population", "disease_rate"), class = "data.frame",     row.names = c(351L,      251L, 477L, 349L, 56L, 434L, 299L, 248L, 301L, 262L, 115L, 325L,      402L, 238L, 73L, 376L, 454L, 470L, 6L, 359L, 473L, 233L, 448L,       121L, 35L, 253L, 324L, 369L, 492L, 463L, 400L, 365L, 245L, 157L,      264L, 101L, 367L, 142L, 109L, 389L, 158L, 254L, 457L, 68L, 219L,      395L, 450L, 380L, 474L, 199L))     , coords.nrs = c(4L, 3L)     , coords = structure(c(634873.420766408, 541790.625187116, 649805.450219927,      536333.978570239, 581957.093009378, 607202.361088042, 617240.615621174,      505627.220834258, 620543.999367496, 592414.016643991, 579995.994396474,      636851.153664358, 593891.033349101, 594650.473384273, 607103.681876253,      659509.454474378, 604222.534339966, 592766.490473358, 625515.430661865,      545043.782766906, 663659.275576225, 626566.467170621, 687406.871345609,      579464.084578232, 594577.343169983, 650081.936891944, 661853.556699763,      620028.294881262, 622881.614147092, 714608.702621971, 620991.249432193,      531188.877201674, 654835.610827768, 527241.805893339, 604428.976031172,      688457.08368477, 613888.674693582, 616783.067305461, 500082.617682936,      609515.575391672, 611609.283708876, 689412.437797476, 675561.054104031,      601548.525491549, 656484.375080421, 677341.365311871, 616862.413097344,      644065.67844776, 475677.631939326, 617193.359923705, 5144526.77927643,      5328308.33322789, 5216468.95632244, 5278378.21530125, 5242407.33311755,      5107204.72128753, 5104424.13808425, 5346313.09064684, 5291125.71228228,      5260872.67506553, 5296414.86191041, 5153372.04174738, 5107141.20462901,      5126046.0142604, 5188835.84829657, 5241495.91305329, 5120799.1233737,      5077612.68066607, 5173630.04242822, 5351027.28834978, 5231260.93265272,      5321245.1139323, 5239845.72123907, 5257533.19157413, 5093198.30043019,      5176787.1924563, 5208692.02318177, 5301924.71404294, 5307513.36239273,      5170246.64046221, 5129845.77701215, 5284862.55527842, 5147875.41139468,      5295710.09388113, 5293345.39302848, 5189299.58825685, 5165031.11113162,      5136185.92439977, 5334176.39641511, 5166895.32812857, 5279273.17144778,      5236974.38111571, 5222457.79639944, 5279521.25697787, 5171514.42644342,      5244570.31251435, 5113000.12179841, 5170833.89650194, 5346485.72226854,      5068270.63886515), .Dim = c(50L, 2L), .Dimnames = list(NULL,      c("long.proj", "lat.proj")))     , bbox = structure(c(475677.631939326, 5068270.63886515, 714608.702621971,      5351027.28834978), .Dim = c(2L, 2L), .Dimnames = list(c("long.proj",      "lat.proj"), c("min", "max")))     , proj4string = new("CRS"     , projargs = NA_character_     )     )


title: is my open source web GIS software selection ok?

text: Which are the best open source tools for building a web gis application having PAN, Zoom, Query/search, Label etc options ?. I have decided to use Postgre sql for database and geoserver as middleware. I am a bit confused in selecting the software to design the GUI for the client side. I am a novice programmer , will geoexplorer be a good option or are there some other tools which can be easily customized for my work. ??


title: Watermarking photos using their geotag and create a shapefile

text: I have geotagged photographs and I would like to place a watermark of the photos location on the front for easy reference as well as creating a shapefile from the coordinates. Does anyone know of any open source software which could do this or a python or arcpy script? This link is an example of what Iam trying to acheive Thanks.


title: How to truncate field names to 10 characters in model builder

text: Apologies that I have not posted an image of my model, I am away from my computer for a while but am trying to sort this out anyways. I have a number of different animal GPS location points, and a bunch of different habitat characteristics. I have built a simple model that uses the Near tool to calculate the distance from each point, create a new field using Add Field, naming it after the habitat characteristic shapefile name (using Parse Path and inline variable substitution), and then populates that added field with the NEAR_DIST values. It then continues on for all of my habitat characteristics. I built the model a while ago and have since changed the names of the habitat characteristics (and there are a bunch of them) and unfortunately in the renaming them they are all now above the 10 character limit for a field name. What I was wondering is if there was any way to truncate the habitat characteristic shapefile title to 10 characters so that the add field tool doesn't error out. Thanks in advance, and if you need any more information please let me know. I tend to leave things out by mistake when I ask questions on here. Edit - Examples of habitat names are - 2001 Upland, 2001 Harvest, 2001 Deciduous Essentially I need to remove the year and the space and they will all be under 10 characters. Also, at home now so here is the tool I created. It runs through but doesn't successfully "Add Field" because its trying to name it after the "Near Feature" (the habitats whose names are too long) and so it just skips it. ![enter image description here](http://i.stack.imgur.com/AHEPh.png)


title: How to load map tiles according to the drag/touch event

text: We are building map apis for javascript integrated with our own map data server(tiles). (Though we can use openlayers or something else, but we can not for some thrid-part reasons). In a word, we are re-writing something like openlayers. The following is what we have done:   1. When a map is created, it will have the container(dom element),center(LatLng), and zoom . Then we will orginaize the tiles inside the container according to the center and zoom.    2. When user drag the container. We will record the increasement of the mouse position. Then caculate the new center for the map, and then re-caculate the required tiles and orginaize them.  Like this:               onMouseMove(dx,dy){       var dinx=caculateDegreeInCurrentResolution(dx);       var diny=caculateDegreeInCurrentResolution(dy);       this.center=new LatLng(this.center.lat+diny,this.center.lnt+dinx);       this.refresh();     }      It works as expected, but we found that it is hard to drag the map. It is not smooth like openlayers or other map library during the drag operation. So I wonder if we are going the wrong way?


title: Analyse spatial data and extract results in tabular form in QGIS

text: I'm sorry if this question is badly framed, I only started with GIS and QGIS 2 days ago so I don't know any terminology yet - I'll try to explain as best I can. I have 2 layers in QGIS (1.8, on Win7), one with geological information, and the other with the state boundaries. I want to get out a report which tells me the composition of each state, i.e. a table of the form:               State name    Rock type 1    Rock type 2     State1            32 %           68 %           State2             5 %           95 %           I guess this is a pretty basic thing to want to do, but I'm having trouble finding how on my own because I don't know what to search for. Sorry if bad question, thanks for any help.


title: How to declare ArcSDE workspace for Geoprocessing tool in ArcGIS 10.1 using ArcObjects and c#?

text: What is the proper way to declare an ArcSDE workspace for Geoprocessing tools in c# .net? I am using ArcGIS 10.1 SDE (MS SQL) on SQL Server 2008. I want to write code so that every night the server will fire a program and compress one of my ArcSDE databases. I thought I could do this in c# using the Compress tool in the data management namespace but have run into a snag. I am told my workspace is invalid. Esri help docs are not helpful, as they just show generic **help** and my internet search has lead nowhere. Can anyone tell me how to define an ArcSDE workspace in 10.1 in net? Right now I am using @"Database Connections\CME.sde". Here is the error: Executing: Compress "Database Connections\CME.sde" Start Time: Fri Oct 12 13:40:47 2012 Failed to execute. Parameters are not valid. ERROR 000837: The workspace is not the correct workspace type. Failed to execute (Compress). What should I type for the in_workspace? * * *               private void button1_Click(object sender, EventArgs e)             {                      Compress Compressor = new Compress();                 object in_workspace = @"Database Connections\CME.sde";                                     Compressor.in_workspace = in_workspace;                      Geoprocessor GP = new Geoprocessor();                 RunTool(GP, Compressor, null);                       }             private static void RunTool(Geoprocessor geoprocessor, IGPProcess process, ITrackCancel TC)             {                      // Set the overwrite output option to true                 geoprocessor.OverwriteOutput = true;                      // Execute the tool                             try                 {                     geoprocessor.Execute(process, null);                     ReturnMessages(geoprocessor);                      }                 catch (Exception err)                 {                     MessageBox.Show(err.Message);                      ReturnMessages(geoprocessor);                 }             }             private static void ReturnMessages(Geoprocessor gp)             {                 if (gp.MessageCount > 0)                 {                     for (int Count = 0; Count <= gp.MessageCount - 1; Count++)                     {                       Console.WriteLine(gp.GetMessage(Count));                          }                 }                  }         }


How does ArcGIS determine which transformation to use when projecting on the fly? Where this would come in would be going from NAD 83 to WGS 84. Multiple transformations could be considered so how is one picked?


title: ArcView 10.1 PDF and Popup

text: When converting your map into a PDF from ArcView 10.1, is there a way to make the PDF Map to have a pop-up when you click on a tract?


title: Make Openlayers request WMS tiles only within a certain bounds

text: Is there a way to actually tell openLayers not to request tiles out of particular bounds? I want the rest of map and other layers still accessible, so maxExtent on map is not an option for me. I just want to restrict one WMS layer to particular bounds. If view port is inside the bounds and user zooms in fine, openlayer requests new tiles! but if viewport is outside, the WMS should only fill the bounded area with tiles from that WMS layer. Does that make sense or is that even possible?


title: Does qgis_mapserv.fcgi segfault because ECW data is used?

text: I've compiled qgis from source, and when I call qgis_mapserv.fcgi in my browser with REQUEST=GetProjectSettings (and my map file set via map=) is returns the XML correctly. But if I then try to actually view a map, it seg faults qgiswebclient.html?map=/var/www/html/testmap.qgs               mod_fcgid: error reading data, FastCGI server closed connection, referer: http://www.ground2control.com/test/site/qgiswebclient.html?map=/var/www/html/GIS/Stonelees.qgs     Premature end of script headers: qgis_mapserv.fcgi, referer: http://www.ground2control.com/test/site/qgiswebclient.html?map=/var/www/html/GIS/Stonelees.qgs     mod_fcgid: process /var/www/cgi-bin/qgis_mapserv.fcgi(27652) exit(communication error), get unexpected signal 11          kernel: qgis_mapserv.fc[27692]: segfault at 0 ip 000000000043dc02 sp 00007fffa9967bf0 error 4 in qgis_mapserv.fcgi[400000+104000]      Actually, it works fine with the helloworld.qgs file supplied with QGIS. The main difference I can see is that my qgs file uses some ECW files. I've made sure that qgis was compiled with ECW support, but I wonder if that's the problem. Any ideas on how I can discover what is causing the seg fault? UPDATE: Running the mapserv CGI through GDB:               Program received signal SIGSEGV, Segmentation fault.     0x00007ffff59445a3 in qHash(QString const&) () from /usr/lib64/libQtCore.so.4          (gdb) bt     #0  0x00007ffff59445a3 in qHash(QString const&) () from /usr/lib64/libQtCore.so.4     #1  0x000000000044c2a3 in QHash<QString, QHashDummyValue>::findNode (this=0x87eac0, akey=..., ahp=0x0) at /usr/include/QtCore/qhash.h:875     #2  0x0000000000477b92 in QHash<QString, QHashDummyValue>::contains (this=0x87eac0, akey=...) at /usr/include/QtCore/qhash.h:867     #3  0x0000000000475797 in QSet<QString>::contains (this=0x87eac0, value=...) at /usr/include/QtCore/qset.h:86     #4  0x000000000045cc7e in QgsProjectParser::mapLayerFromStyle (this=0x87e9e0, lName=..., styleName=..., useCache=true) at /root/qgis/qgis-2.0.1/src/mapserver/qgsprojectparser.cpp:1365     #5  0x00000000004ae4a0 in QgsWMSServer::applyOpacities (this=0x82d4d0, layerList=..., vectorRenderers=..., rasterRenderers=..., labelTransparencies=..., labelBufferTransparencies=...)     at /root/qgis/qgis-2.0.1/src/mapserver/qgswmsserver.cpp:2045     #6  0x00000000004a2363 in QgsWMSServer::getMap (this=0x82d4d0) at /root/qgis/qgis-2.0.1/src/mapserver/qgswmsserver.cpp:674     #7  0x000000000043eb60 in main (argc=1, argv=0x7fffffffe5a8) at /root/qgis/qgis-2.0.1/src/mapserver/qgis_map_serv.cpp:533


title: Trying to compile QGIS but it complains it can't find QGIS library

text: I'm trying to compile QGIS in Windows7 with a MSYS, with the objective of making a QGIS Windows App using the API. I am using this document. I'm using CMake-GUI and it gives the error > Warning Target "qgis_core" requests linking to directory "c:\OSGeo4W\bin". > Targets may link only to libraries. CMake is dropping the item.] Screenshot I think its looking for qgis_core.lib, but I don't think it exists because I want to compile it... so this error makes no sense to me.


title: How do you auto populate a field when a new point feature is digitized?

text: Is there a way to populate a field automatically when a new point is dropped on the map? We have our map in State Plane Feet, and we want to the database to automatically add WGS xy coordinates to fields in the table when a new point is added. Something like a calculated field, but always there? Forgive me if this is basic, as I am an systems admin and google has not helped me.


title: How to move line segment to certain direction using ArcGIS for Desktop?

text: How can I move a Line Segment to a certain direction? In my case i want to move the Line up north, so i need to change the y-coordinate somehow, I'm just not sure how to do this in ArcMap. Im using ArcGIS 10.1


title: Giving properties to Zoomify layers in OpenLayers?

text: I am using `OpenLayers.Layer.Zoomify` to show large image files and this works like a charm. The only thing is that I would like to add some options to that layer ( opacity and transitionEffect) but it doesn't seem to work. Code:               var zoomify = new OpenLayers.Layer.Zoomify( "Zoomify",       zoomify_url,        new OpenLayers.Size( zoomify_width, zoomify_height ),       { opacity: 0.5, transitionEffect: 'resize'}     );      Isn't this the right way of doing it? Or is it not possible to add these properties to Zoomify layers? Nothing happens to my zoomify layer.


title: ArcGIS field calculator: VBScript get data from previous row

text: When using the field calculator, how do I get data from the previous row in VBScript? In Python I'd do something like this:               for index, row in enumerate(rows,start=1):         row[index-1] += 5      I'm using ArcGIS 9.3.1. Let's assume I have a column `[Distance]` and all it's rows are `0`. How do I increment every row with 5 starting at row 1. So the first row remains `0`, second row'll be `5`, third `10` and so on.


title: Looking for good load testing tools to test my web GIS

text: I am looking for good load testing tools to test my web GIS which runs on Apache/ESRI AU's Dekho application (Flash based, so would require AMF support within the load testing tool). I am currently working SoapUI/LoadUI, but I am having a tough time creating steps for it, as there are numerous AMF calls from the client to the server. Some calls even have multiple calls embedded into it. Has anyone used any load testing tool with GIS? Many Thanks   Ujjwal