title: How do I make a public download page for my CartoDB table?

text: I have some CartoDB data that I'd like to make available to people as shapefile, geojson and csv. I just came across http://opensandiego.cartodb.com/tables/5037/public/map which has a great map preview, with downloads in different format and a cool table view. I have a number of public CartoDB tables, but I can't find how to share my data in that nice page. Is there a button to create the url? Or some way I can figure out what that URL should be?


title: Add a polygon to PostGis from Google Maps

text: i have some experience with Google maps api, I have an app and currently calculate nearest points to locations via distance functions (great circle) and having a lot of waypoints on MySQL. I want to start with PostGis in order to do this a better way or have been able to give a user the ability to build a polygon in GMAPS (point and click on map, saving each lat,lon, connecting the last with the first. I've just downloaded postgresplus to start with this but don't have reference to this: 1\. create a table of polygons 2\. insert via php each polygon to PostGis 3\. be able to query the polygon to see if a point is inside or outside a polygon I know all this is very efficent with PostGis but have 0 experience with it. does someone knows a link or site to learn all this basics of postgis thanks NA


title: Can someone clarify some details about freehand line tracing in QGIS?

text: I have read several tutorials and questions here on this site, but I have some lingering questions/doubts that I would like clarified. My project is tracing roadways from a 1926 map so that I can compare infrastructure then and now. I have successfully georeferenced my 1926 raster map onto a contemporary shapefile, reduced transparency, etc. I have enabled freehand line editing, which is the best way to accomplish this task, it seems. My map has two kinds of roads shown that I want to trace: the paved roads, and the unpaved roads (marked by full, bold lines and staggered, broken lines). So, then, here is what I would like clarified. 1) In the New Vector Layer box, an option is given for Name, Type, Width... I have named the feature line 'roads,' which is relevant for question #2. What exactly is the Type option (text data, whole number, decimal number, etc.) for, what is it not for, and could someone give a quick example of when this option might be useful? 2) When I trace with the freehand edit tool, I realize that when I unclick, several points/nodes/features (whatever they are called!) are placed along the traced line _as if_ I clicked them each individually. OK... so there appears a popup box asking me for 'id' and 'roads'--'id' is blank, and 'roads' is 'NULL.' What exactly is going on here with this? So far, I have simply left the box as-is and clicked 'OK.' Would I regret this later if I wanted to edit the lines or something? What is the point of this box, the information in it, and how is it useful? 3) The paved roads I want to trace are not connected to each other--some paved roads in the north, south, etc. So far, what I have been doing is tracing a line of paved roads as described above, and then moving my location in the map to do the same in the north or other areas. Following this procedure as described, will all the lines traced be counted/recognized as the same feature/network of feature? Is it alright to trace a road coming into a city from one direction, unclick, and then trace the road coming into the city from the other direction, and will it still be part of the same road? 4) How do I trace the unpaved roads once I get the paved roads all traced? (I suspect that I would need to allow for this when creating the new shapefile layer from the start, by adding a 'roads' AND another thing like 'unpaved roads,' is this correct? But if I have not done this from the beginning, am I screwed and should start over?)


title: Geometry column naming convention - 'geom' or 'the_geom'?

text: I'm starting my first PostGIS project, and in various books and tutorials I've seen the geometry column labelled either 'geom' or 'the_geom'. Is one more conventional than the other? Furthermore, is there a good reason to use 'geom'/'the_geom' instead of a more descriptive name for the geometry column? (E.g. 'centre_point' for the central point of a feature.) We don't label conventional db columns 'the_int' or 'the_string', so why label geometry columns this way?


title: ESRI Geoprocessing service and other servers

text: I am new to ESRI GeoProcessing (GP) services so please excuse the naive question. I have an FLEX application running off of ArcGIS Server. I need to push data from the ArcGIS Server to an analytic server that calculates all sorts of interesting metrics using R that I cannot do within ArcGIS. So my question is this: Can a GP accept user inputs (e.g. run on these 5 points) and then push a foreign key to _another_ server for further processing ? The other server would do its magic and then send the results (and foreign key) back to ArcGIS server which displays to the user. Can someone point me to some resources on the web I can read up on this ? Thanks ! -Don


title: Why won't OGR connect to my PostGIS database, but QGIS will?

text: Trying to connect to my PostGIS database with some success... One hiccup is that OGR will not connect to my server from a remote machine, though QGIS can connect to it just fine from that same machine. Here is the connection string I'm using:               ogrinfo -ro PG:"host='machineName' port='5432' user='userName' dbname='postgisDB'"      This connection works in the following scenarios:   * Connect using 'localhost' from PostGIS server machine   * Connect using 'machineName' from PostGIS Server machine   * (as stated) Connect using 'machineName' through QGIS from remote machine I've tried adding a fully qualified domain name, as well as a password, to no avail. Thanks!!! -m


title: Is there a way to convert GeoTIFF to DEM?

text: Simple enough question. The NASA Aster data comes in GeoTIFF format, but the tool I'm using (Sketchup) only accepts DEM. Is there a way to convert? I did a lot of searching, and only found the program GEOTIFF4, but it's obsolete, it's Windows-only, and source code is not available.


title: Batch clip DWG polyline

text: I have about 170 DWG files that need to be clipped to their corresponding tile. Each tile is 14,000ft wide, 8,000ft height. The grid: ![enter image description here](http://i.stack.imgur.com/KYaZj.jpg) Each DWG file overlaps the tile border by 800ft: ![enter image description here](http://i.stack.imgur.com/5TIBh.jpg) Unfortunately, each DWG file is in its own directory, so I'm not sure if that makes coding a batch clip more difficult. The file we're working with will always be named like "###_2ft.dwg" ![enter image description here](http://i.stack.imgur.com/oKHcM.jpg) So my question is: How can I efficiently clip these DWG's and output them as feature class files? I have thought about creating a fGeoDB to copy each DWG polyline file into, however, for whatever reason when expanding the DWG layer it takes ~30 seconds to show the annotation, point, polyline, etc. I've checked out some guides using iteration in ModelBuilder, but I can't get the same results. Specifically posts like this. In that post, a user mentions constructing the geoprocessing command in notepad and using excel to help with filenames, but I'm not sure what the commands would be.


title: Importing UK Land Registry data into Quantum GIS

text: We have purchased data from the UK Land Registry to import into Quantum GIS however when I do so using the Coordinate referencing systems of OSGB 36/ British National Grid (EPSG:27700) the data appears to be stretched compared to our existing paper OS maps. I assume that I have a setting somewhere which is incorrect but being new to GIS and Quantum GIS I do not know where to look. Can anyone help?


title: problem with WFS Data, GDAL, Python, ArcGIS, & Defining Projection

text: I'm pulling in data in gml2 format with a EPSG#4267 from a wfs. I create a datastream in python and then run the following from the python window in ArcMap:               data = datastream.read()          #write the gml data to a gml file     rasterdata = open("E:\\Melo\\USGS" + "\\Roads.gml", 'w')     rasterdata.write(data)     rasterdata.close()          #Run the ogr2ogr to convert to shp     cmd = ["ogr2ogr", "-f", "ESRI Shapefile", "E:\\Roads.shp", "E:\\Roads.gml"]     subprocess.call(cmd)          #Define the projection     infc ="E:\\Roads.shp"     prjfile = "E:\\4267.prj"     arcpy.DefineProjection_management(infc, prjfile)      Unfortunately, when i try to add the layer in ArcMap 10.1 it "appears" blank. Typically, I've only seen this happen when their was something wrong with the projection. I know there is data in the file because I've converted it successfully using alternate methods. I installd GDAL according to these instructions exactly (1500/32 bit/2.7). I did NOT download and install the ArcGIS bindings. Do I need them? Any ideas on what could be causing the layer to appear blank in ArcMap?


title: How can I figure out whether a sphere is intersecting the Earth between a given lon/lat range?

text: Given a lon/lat range that represents a tile on the Earth's surface, I need a computationally efficient way to figure out if a spherical volume defined by a point (assume ECEF coordinates) and a radius intersect with the tile. It can be assumed that the maximum tile side length won't exceed 22.5 degrees and that all lon/lat area coordinates are aligned along traditional map tiling schemes (tiles won't cross the antemeridian, etc). In the following picture, the blue ellipsoid represents Earth and the green sphere represents a volume of interest. ![enter image description here](http://i.stack.imgur.com/o2Efi.png) Given a lon/lat area, you can generate ECEF coordinates for the 4 extreme points and check whether any of these points lie within the sphere through a simple distance calculation. However, this method doesn't always work since you're ignoring curvature. The following picture shows this, with the tile shown in black and a possible sphere of interest in green. ![enter image description here](http://i.stack.imgur.com/6KZEJ.png) You can add more points to test with the lon/lat tile (ie. check the center of the tile along with the 4 corners) but you're just trying to improve an approximation, and this will fail based on the position and size of the sphere, so I want a more analytical solution. I'm looking for a solution I can implement myself in code and not for an application or library. I'd like to use an ellipsoid representation of Earth if possible, but if that ends up being too complex, representing the Earth as a sphere is okay too.


title: Problem loading PyQGIS

text: There is a problem on my system when starting QGIS. I get this error. I am running openSuse Tumbleweed 64bit and openSuse 12.3 64bit. Same problem on both systems. Couldn't load PyQGIS. Python support will be disabled. > Traceback (most recent call last): File "", line 1, in RuntimeError: the sip > module implements API v10.0 but the qgis.core module requires API v9.2 > > Python version: 2.7.3 (default, Apr 14 2012, 08:58:41) [GCC] > > QGIS version: 1.8.0-Lisboa 'Lisboa', exported > > Python path: ['/usr/share/qgis/python', '/home/evert/.qgis//python', > '/home/evert/.qgis//python/plugins', '/usr/share/qgis/python/plugins', > '/usr/lib/python27.zip', '/usr/lib64/python2.7', '/usr/lib64/python2.7/plat- > linux2', '/usr/lib64/python2.7/lib-tk', '/usr/lib64/python2.7/lib-old', > '/usr/lib64/python2.7/lib-dynload', '/usr/lib64/python2.7/site-packages', > '/usr/lib64/python2.7/site-packages/FontTools', '/usr/lib64/python2.7/site- > packages/PIL', '/usr/local/lib64/python2.7/site-packages', > '/usr/local/lib/python2.7/site-packages', '/usr/lib64/python2.7/site- > packages/gtk-2.0', '/usr/lib/python2.7/site-packages', > '/usr/lib/python2.7/site-packages/setuptools-0.6c11-py2.7.egg-info', > '/usr/lib64/python2.7/site-packages/wx-2.8-gtk2-unicode'] The problem would seem to be the API version implemented in the sip module versus the API version implemented in the qgis.core. How do I go about fixing the problem once and for all, as it has been a recurring problem on various occasions over the years. My apologies if this question has been asked and answered here already. Did do a search but could not find this specific problem.                                      ------Edit------------      Used Yast to install Qgis, Grass ect, after adding the required repo's.               <URL: http://download.opensuse.org/repositories/devel:/languages:/python/openSUSE_12.3/     Category: YUM>          <URL: http://download.opensuse.org/repositories/Application:/Geo/openSUSE_12.3/     Category: YUM>      There is listed in the Geo repo develepment versions of both Qgis 2 and Grass 7. I have not installed them at any time on either system.


title: How to set Definition Query on ArcPy layer from shapefile?

text: I'm trying to determine options to set a layer definition query using ArcPy. I know that it is possible when using arcpy.mapping.ListLayers(). However, in this scenario I am loading in shapefiles, and using arcpy.MakeFeatureLayer_management() in ArcPy. I want to set a definition query on the layer that I create from using this geoprocessor. Is this possible?


title: How can I import geometries from an invalid GML quickly?

text: I have the following gml geometry :               <gml:Polygon xmlns:gml="http://www.opengis.net/gml" srsName="EPSG:2154">       <gml:exterior>         <gml:LinearRing>           <gml:posList srsDimension="2">     657983.0000000068 6860322.000000936 657655.0000000069 6860452.000000936 657172.000000007 6859562.000000938 656962.000000007     6859846.000000937 657140.0000000069 6860890.000000934 657169.0000000069 6861180.000000931 657362.0000000069 6861189.000000934     657395.0000000069 6860859.000000935 658016.0000000068 6860653.000000936 657983.0000000068 6860322.000000936 657983.0000000068     6860322.000000936 658016.0000000068 6860653.000000936 657395.0000000069 6860859.000000935 657362.0000000069 6861189.000000934     658110.0000000068 6861132.000000933 658706.0000000066 6861564.000000933 659428.0000000065 6861377.000000932 660042.0000000064     6861175.000000932 660094.0000000064 6860854.000000935 659478.0000000065 6860671.000000935 657983.0000000068 6860322.000000936     659428.0000000065 6861377.000000932 658706.0000000066 6861564.000000933 658110.0000000068 6861132.000000933 657362.0000000069     6861189.000000934 657169.0000000069 6861180.000000931 657099.0000000069 6861841.000000932 657684.0000000068 6862313.000000931     658585.0000000066 6863988.000000925 659117.0000000065 6863683.000000926 659784.0000000064 6864414.000000923 660744.0000000063     6864361.000000924 661242.0000000062 6863433.000000928 661899.000000006 6863173.000000927 661962.000000006 6862500.000000929     659959.0000000065 6861848.000000931 659428.0000000065 6861377.000000932 657983.0000000068 6860322.000000936     </gml:posList>         </gml:LinearRing>       </gml:exterior>     </gml:Polygon>      I'd like to import it as **easily** and **fast** as possible into QGIS.   * I've tried direct GML import => not a valid GML (and i quite agree with that ;))   * I can use QuickWKT => But then I need to put commas every two coordinates   * I've tried "Add delimited Text Layer" => But it does seem to work mainly for CSV file with coordinates on each line. Any ideas ? Thanx in advance


title: Whether to use ArcGIS Platform (from Esri) or Google Maps API for large number of polygons?

text: I am in the planning phases of building a web based application which will be served on IOS devices through a custom application with HTML container. The functionality will be to display a number of icon markers as well as polygons. The core functionality is fairly simple which is to render markers and polygons. The issue I am having is with a feature which I have agreed to which is dynamically adding and removing polygons. I can certainly do this with client side on google maps. However performance of polygons on client side is very poor. There is also a limit of ~10 KML files at any one time to be added to a google map. The polygons need to be in separate KML files so that I can add and remove them. My question is does ESRI GIS maps handle large amounts of polygons being added/removed better than google maps apiv3?


title: Cannot read property 'wrapDateLine' of null - Open layers

text: My code is giving the error: Uncaught TypeError: Cannot read property 'wrapDateLine' of null OpenLayers.js:192 OpenLayers.Map.OpenLayers.Class.isValidLonLat OpenLayers.js:192 OpenLayers.Map.OpenLayers.Class.moveTo OpenLayers.js:187 OpenLayers.Map.OpenLayers.Class.setCenter OpenLayers.js:184 init chapter_8_ex9_using_multiple_maps_custom_event_app.html:38 onload chapter_8_ex9_using_multiple_maps_custom_event_app.html:63                   <!DOCTYPE html>                <meta charset='utf-8' />      <title>My OpenLayers Map</title>


title: Best Practices : ArcGIS for Javascript

text: I have just started developing with **ArcGIS for Javascript** and as I'm learning Dojo I'm seeing a reason to modularize development for a more plug-n- play concept. I can't seem to find any good articles on Best Practices for Web Map development and am wondering if there are anything resembling industry standards for setting up new web-map applications in this environment.


title: Integrating a Canadian WMS with ArcServer

text: We are trying to load a Canadian WMS to our ArcServer and I am told the url should have "arcgis" or "rest" in it for it to work. This WMS-Toporama link works find in ArcMap with "Add WMS Server" tool but our web site person cannot seem to load it with our ArcServer. I am not familiar enough with the web site side of things to give a good answer to this problem. Any ideas here would be most appreciated.


title: Having trouble with clipping tool giving incorrect shape

text: I am trying to clip a shapefile of the Midwest using a Great Lakes shapefile so that it matches my political areas boundary. Whenever I use the Clip tool, it outputs the lakes clipped to states, whereas I need the states clipped to lakes. No matter which input and output I use, it gives the same result.


title: How I can get coordinates by clicking on the QGis Canvas with PyQGIS?

text: Is there are any procedure on QgsMapCanvas (pyqgis, qgis) in order to obtain a pair of coordinates by clicking on the map canvas? Are there any python code I can see?


title: Calculate distance using 2 sets of lat/long values in the same record

text: I have a point layer with over 90K records in it. For each record I have lat/long field for incident, and lat/long for engine location. Is there a way to find the distance between Incident point and Engine Location point? So far, I created points using "Make XY Event layer" using lat/long for incidents. And again I run the same tool for Engine locations. When I run the Near tool in ArcGIS it will find the nearest engine point around the incident and will give me the distance to the nearest engine but not necessarily the distance from the engine location on the same record. How can I calculate distances using the lat/long values for incident and engine location in the same record?


title: How to create randomly points within polygons for each row of a dataframe matching a Polygon ID with R

text: My question is similar to this post, but because there are no "perfect" answer and above all no reproducible example I open this question. I have a "SpatialPolygonsDataFrame" and a data frame. Both have a variable of common factors. I would like to generate a random point for each observation of the data frame within the according polygon. I am starting to learn GIS with R. I read the first chapters of Bivand et al. 2013 but i am still strugling with the basics. I tried with the function spsample from the _library(sp)_ but I am puzzling with the variable n of this function: it's an **(approximate) sample size** In my example there are 3 polygons and 30 observations. The common variable is **MatchID**. I need a SpatialPolygonsDataFrame with 15 points in the first polygons ("abcd"), 10 in the second ("efgh") in 5 in the third ("ijkl"). See               table(df$MatchID)      Example:                library(sp)     # Definition of the CRS     poly.crs <- CRS("+proj=utm +zone=36 +datum=WGS84 +units=m +no_defs +ellps=WGS84 +towgs84=0,0,0")          # Definition of 3 Polygons     poly.a <- Polygon((matrix(c(-4.653724,1.2210259,-3.803160, 1.3799638, -3.245480, 0.1665987,-4.666098, 0.1097523, -4.653724, 1.2210259), nrow=5, ncol=2, byrow=T)))     poly.b <- Polygon((matrix(c(-5.820343, 2.675320,-5.427519, 3.062975,-4.701119, 2.819967,-4.555540, 1.855489, -5.050758, 1.388308, -5.783673, 1.572882,-5.820343, 2.675320), nrow=7, ncol=2, byrow=T)))     poly.c <- Polygon((matrix(c(-3.758639, 2.873654, -3.273958, 3.417311, -2.213099, 2.935320, -1.972031, 1.945992, -3.033510, 1.279312, -3.709359, 1.844633, -3.758639, 2.873654), nrow=7, ncol=2,     byrow=T)))           # Making a SpatialPolygons     polys.a = Polygons(list(poly.a), "pa")     polys.b = Polygons(list(poly.b), "pb")     polys.c = Polygons(list(poly.c), "pc")     Spolys = SpatialPolygons(list(polys.a,polys.b,polys.c), 1:3,  proj4string=poly.crs)          # Making a SpatialPolygonsDataFrame     data.Spolys<- (data.frame(MatchID=c("abcd","efgh","ijkl"), row.names=row.names(Spolys)))     Poly <- SpatialPolygonsDataFrame(Spolys, data.Spolys, match.ID = TRUE)               ## Data frame that should be converted in a SpatialPointsDataFrame thanks to the variable "MatchID"     df <- data.frame(               MatchID=c("abcd","abcd","abcd","efgh","efgh","ijkl"),               V1 = 1:30,               V2 = "a"     )      I need a function to generalize the idea below (looking for the matching with the factor) and giving exactly the number of observation I need (and not an approximate random size)               spdf.a <- SpatialPointsDataFrame(spsample(Poly[Poly$MatchID=="abcd",], n = nrow(df[df$MatchID=="abcd",]), "stratified"), df[df$MatchID=="abcd",])


title: How likely is it that Smartphone Camera GPS coordinates are wrong?

text: According to the coordinates printed directly on this photo of my late son Kelvin here:               https://www.dropbox.com/sc/w6tvwuvbosuh9rk/AADzT-3vDXcFQ9LnKK1BYMZUa      ...and the Exif data from the .jpg:               GPSLatitudeRef: N     GPSLatitude: 24, 58, 34.6     GPSLongitudeRef: E     GPSLongitude: 121, 19, 37.768      ...he was in Taiwan on Jan. 5, 2012 (two and a half months before he died). However, I don't think he was ever in Taiwan, or anywhere near there. Is there any explanation for why a camera would generate the wrong coordinates? Maybe the whole thing is "messed up" because the supposed date of the image, as indicated by all the Exif data extractable:               GPSLatitudeRef: N     GPSLatitude: 24, 58, 34.6     GPSLongitudeRef: E     GPSLongitude: 121, 19, 37.768     Software: HTC Laputa Browser     DateTime: 2012:05:10 18:28:30      ...indicates the picture was taken 5/10/2012, but he died 3/17/2012. I can imagine the date of the camera being off as more likely than the calibration of the coordinates, though. The Bing map for " _No. 23, Xinghua Rd, Taoyuan City Taoyuan County, Taiwan 330_ " (which is what I see when I enter " _24° 58' 34.60" N, 121° 19' 37.77" E_ " into my browser). The "software" entry from the Exif data says "HTC Laputa Browser", and the map does show a couple of "HTC"s on the map, whatever that is. I'm intrigued and confused. Can anybody shed any light on this?


title: (MapServer Osm) Map Accessibility

text: I have a tile server (MapCache) with a personalized map. For client side's rendering I use OpenLayers. I'd like my personalized cache to be accessible only on my web site, i.e. I don't want another website to use the mapcache on my MapCache server. Is there any solution ?


title: Does QGIS Server support QGIS project with Oracle layers?

text: QGIS recently got native Oracle support. I was wondering if QGIS Server supports generating WMS and WFS from a QGIS project with Oracle Spatial layers. I just tested this but I did not succeed...


title: Why do MGI and WGS 84 data not line up in ArcGIS for Desktop?

text: When i reproject raster from MGI 6 zone to the WGS84 it makes a mistake (it is moved for 380 m )See example below:![enter image description here](http://i.stack.imgur.com/N6pne.jpg) I know that the problem is in MGI 6 projection (it is not accurate),but I do not know how to fix it.


title: QGIS translate BAG surface to XYZ

text: I have a BAG (Bathymetry Attributed Grid) from NOAA which I wish to export to XYZ. The grid imports fine as a raster and plays nicely with most of the raster tools that I have tried. It does, however, have a second band which is causing me difficulties: ![BAG_metadata](http://i.stack.imgur.com/jMl1D.jpg) I have been trying to figure out the -b 1 function in gdal_translate, but am not having luck. The function either fails or continues to export both bands, which results in a file too large to edit in excel. A sample of my gdal_translate settings with -b 1 manually set (though I have tried using the -b 1 at various locations throughout the code block with no luck) ![gdal_translate_settings](http://i.stack.imgur.com/fAZJL.jpg) Have tried:   1. gdal_translate -of XYZ "location\H12224_MB_16m_MLLW_combined.bag" C:/Users/bmack/Documents/test7.xyz -b 1 _(exports both bands)_   2. gdal_translate -of XYZ "location\H12224_MB_16m_MLLW_combined.bag" -b 1 C:/Users/bmack/Documents/test7.xyz _(exports both bands)_   3. gdal_translate -b 1 -of XYZ "location\H12224_MB_16m_MLLW_combined.bag" C:/Users/bmack/Documents/test9.xyz _(exports both bands)_ Any suggestions would be greatly appreciated! ETA: I downloaded version 2.2.0 for this.


title: Updating records based on value in another field using arcpy.da.UpdateCursor?

text: I am trying to write file names and paths to two separate fields in a feature class, and the value for each record is dependent on finding a match in a list. When I run my code it writes the same values for each record, so something is going wrong. I think its finidng the first match and writing that value over and over, but I want it to look at the value in the field my search cursor is set to, and then find the match in the list and write that out to the two fields. What am I missing?               import arcpy, os     photos_path = r'P:\Records\GIS\Projects\D04_OHS\OverheadSignStructurePics'     OHS_FC = r'D:\Records\GIS\Projects\D04_OHS\D04_OHS.gdb\D04_OHS'     ID_fld = 'OSMI_ID'     file_nm_fld = 'LINK_FILENAME'     path_fld ='LINK_PATH'     photos_list = []          for dirpath, dirnames, filenames in os.walk(photos_path):         for filename in filenames:             fullPath = dirpath+'\\'+filename             if fullPath not in photos_list:                 photos_list.append(fullPath)          arcpy.AddField_management(OHS_FC,file_nm_fld,'TEXT',100)     arcpy.AddField_management(OHS_FC,path_fld,'TEXT',300)     srch_cursor = arcpy.da.SearchCursor(OHS_FC,ID_fld)     updt_cursor_fn = arcpy.da.UpdateCursor(OHS_FC,file_nm_fld)     updt_cursor_fp = arcpy.da.UpdateCursor(OHS_FC,path_fld)     for row in srch_cursor:         val = str(row[0])         for p in photos_list:             f_name, f_path = os.path.split(p)             if val in p:                 #print val                  for u_row1 in updt_cursor_fn:                     u_row1[0] = f_name                     updt_cursor_fn.updateRow(u_row1)                  for u_row2 in updt_cursor_fp:                     u_row2[0] = f_path                     updt_cursor_fp.updateRow(u_row2)


title: Problems merging OS Terrain50 tiles i.e. gaps at small scales after using "Mosaic to New Raster" Any ideas?

text: I downloaded the OS Terrain 50 dataset in ascii format, then unzipped all the folders and searched for all *.asc files. I then copied and pasted these into a new folder, so that they could be easily located when geoprocessing. I then used the "Mosaic to new raster tool" in ArcMap 10.2 and selected all the asciis as the input rasters, saved to a new geodatabase, defined projection (BNG), pixel type (32 bit float) and number of bands (1) to match the original files. I left everything else as default options. The resulting image has large areas missing in the south east of England (see below). When you zoom in beyond 1:750,000 the gaps disappear. Another strange thing is that when I removed the layer from the table of contents, then re- added it the whole dataset now only shows when zoomed in beyond 1:1,500,000 (and gaps still until 1:750,000). I don't desperately need to show the data at national scale, however just curious as to what might be wrong and whether anyone can help. Thanks, Sarah ![enter image description here](http://i.stack.imgur.com/ho1Vc.jpg)


title: Creating Time-series (Temporal) Map of NDVI?

text: Is there a best practice for a cartographer to create a Time-series or Temporal Map ? I have 83 of the NDVI in an imagine format and should I converted them to polygons and populated them and display them in Monthly series as a field ? This is the one I will be doing a similar to this but it will be mainly the NDVI and how do I put the legend ? NDVI Map using Time-series


title: How to change GDAL version in QGIS to get ECW import to work

text: I’ve previously been successful in importing ecw into QGIS 2.0 but the same methods (for example How to add *.ecw files to QGIS 2.0 on Windows?) have failed with a fresh install of QGIS 2.2 with OSGeo4W on W7 64. Seems to be a GDAL version issue. The OSGeo4W shell tells me I have GDAL version 1.11 but this new installed QGIS is “compiled against GDAL/OGR" version 1.10 according to 'about'. I’ve successfully got the ecw format listed by ‘gdalinfo –formats’ but that’s version 1.11 which I assume QGIS 2.2 can’t see. I note that OSGeo4W64/bin has both gdal110.dll and gdal111.dll. How do I force QGIS to use gdal 1.11 and not 1.10 and so be able to read ecw rasters (which I assume is the problem)? Any ideas gratefully received. * * * Thanks for the response nmtoken. I find that I have a different folder structure from my fresh QGIS2.2 install from OSGeo4W which has no less than 50 gdal files (from “gdal_auth.bat” to “gdalwarp.exe”) in C:/OSGeoW64/bin and one file “gdal_ECW_JP2ECW.dll” which is in C:/OSGeoW64/bin/gdalplugins. No “gisinternals” folder exists for me. I have followed the instructions on the GDAL website and also other answers on gis.stackexchange by manually copying the various ERDAS_ecw_SDKfiles to the appropriate folders. This worked with QGIS2.0 but fails with QGIS2.2. I've put C:/OSGeo4W/bin as the path in the Gdal Tools executables dialog box but QGIS2.2 still fails to recognise ecw rasters. Has anyone successfully imported ecw rasters in QGIS2.2? If so, could they let me (and others) know how they did it! * * * Success! QGIS2.2 64bit now importing all my ecws perfectly. Andre Joost's advice to manually roll back the default Gdal1.11 to Gdal1.10 and reinstall using the OSGeo4W advanced options works. How many users would have worked this one out for themselves I wonder? Either OSGeo4W should default to Gdal1.10 on express install or QGIS should compile against Gdal1.11.


title: Is the label engine or data server responsible for reducing the number of features

text: We are working on an android application which load vector features from the server and render(with labeling) them in the client. Since there maybe a lot of features in a single viewport, so the map would be un-readable if we render and label all the features because there would be so many overlapped features and icons. So reducing the number of the features is necessary. But I am not exactly sure where is the endpoint. The process can be done in the server side or at run-time during rendering. I prefer to do that in the server side because the size the data to be severed would be reduced in this case. But IMO, it is the responsible of labeling engine for conflict detecting and overlap handing, so I am a little confused. What's the acceptable practice?


title: How can I always run an algorithm with standard input shapefiles from mapcanvas in QGIS 2.0.1?

text: I created a model in the graphical modeler and I want to run that model always with two specific polygon shapefiles from my project, e.g. `pst` and `astenot`. Whenever I run my model, it asks me to set the inputs from the list of shapefiles that are in my MapCanvas. To be more specific, I always want `vectorlayer_1='pst'` and `vectorlayer_2='astenot'`. This is the script from my model. Can anyone tell me how I can do that? I'm wondering if there is a way to do that directly in graphical modeler.               ##difference_between_pst/astenot=name     ##vectorlayer_1=vector     ##vectorlayer_2=vector     ##output_alg0=output vector     ##output_alg1=output vector     ##output_layer_alg2=output vector     ##output_layer_alg3=output vector     outputs_0=Processing.runalg("qgis:difference", vectorlayer_1, vectorlayer_2, output_alg0)     outputs_1=Processing.runalg("qgis:difference", vectorlayer_2, vectorlayer_1, output_alg1)     outputs_2=Processing.runalg("qgis:saveselectedfeatures", outputs_0['OUTPUT'], output_layer_alg2)     outputs_3=Processing.runalg("qgis:saveselectedfeatures", outputs_1['OUTPUT'], output_layer_alg3)      Τhanks in advance!


title: kml to shp including a field attribute as name/description (using 2.1.0 Durfour)

text: I have the Durfour version at my disposal. I captured some data in Google Earth, converted it to .shp, made some revision (e.g. added consecutive numbering) and would like to convert it back into .kml. this time however I would like to use the consecutive number field (e.g. objectID) to distinguish the features. I have tried several options: 1: file > "save as kml" option 2-5: file > save as kml > detailing Data Source as follows with ; indicating line division: Name=PVNNo ; Description=PVNType Name=No ; Description=Type Name="No" ; Description="Type" ; (the above with the v-shaped pac man type brackets on either end) 6: Layer > Properties > Fields (selecting the fields I want to incorporate and providing them with a 'dummy' alias (e.g. name, description) All the test files have produced the usual - untitled polygon/point/line when converted back to kml.


title: R for loops in Geospatial modelling environment

text: I am using the isectpolypoly command in GME and would like to iterate over several polygons ("proj_0_ect0207", "proj_1_ect0207"). My code is as follows:               for (i in 0:30)      {          isectpolypoly(in="F:\Users\p3\ResidentialShapefile\Britain_BNG.shp",poly=     ("F:\Users\p3\new.gdb!"Pro_",[i],"_ject0207"), field="GRIDCODE", prefix=[i],         thematic=FALSE, proportion=FALSE, awm=TRUE, min=FALSE, max=FALSE, aws=FALSE);          };      I keep getting the error: > Error: For loop commands must be contained within a pair of braces: { }. I'm following guidance from http://www.spatialecology.com/gme/for.htm Any help appreciated... Cheers


title: Wmslayerinfo null pointer android arcgis

text: yazi = (TextView) findViewById(R.id.txt);                   wmsLayer = new WMSLayer(wmsURL);                 try {           WMSLayerInfo[] wmsLayerInfo = wmsLayer.getLayers();           try {             for (int i = 0; i < wmsLayerInfo.length; i++) {                          yazi.setText(wmsLayerInfo[i].getName().toString());                       }             } catch (Exception e) {                 // TODO: handle exception                 yazi.setText(e.toString());             }     } catch (Exception e) {         // TODO: handle exception         yazi.setText("1"+e.toString());     }      that'is my code and i got this error 1 java.lang.nullpointerexception this error come from wmslayerinfo.lenght. Where is my mistake ? thanks..


title: Field Calculator - replace string function - data value recognised as a column and is "not found"

text: I have the same challenge as GrahamD - Field Calculator - replace string function I have exactly the same challenge. Exactly the same as GrahamD's challenge, p abbreviated for pan. Type_2 as a test column. duplicated. Using underdark and eseglem's formulars/advice here is the copy and paste directly from the expression box:               replace("Type_2", "p", "pan")      All text and quotations in red. ok button greyed out. > Output preview: expression invalid Parser Error: (nothing detailed) Eval > Error: Column 'p' not found I get exactly the same results should I swop Type_2 and p around.


title: Why does the SAGA raster calculator fail when used in a Processing model?

text: I have a problem with the modeler and the saga raster calculator. When I'm using the raster caculator without the modeler it works fine. When I'm using the saga raster calculator in my model I always get the messege: "Oooops! The following output layers could not be open test: C:\Users\miller\AppData\Local\Temp\processing\4a78580e485a40d199a2fa138960042c\RESULTALG2.tif" – Is this a known problem? Some weeks ago everything worked fine.


title: assigning census block demographics to postal data (zip/zip+4)

text: Using Python and R, what would be the best strategy to tie Census Block-level data to zip/zip + 4? I understand we're not matching polygon boundaries (zip codes are not geographies), but somehow service bureaus/lettershops tie addresses to Census demographic data. Do they work off zip+4 centroids? How can this lookup be updated with new zip/zip+4 definitions (i.e. not be bound by the 2010 ZTCA). Thanks for any pointers/links/book references.


title: How to enter quadrant bearings into Copan map check

text: I am attempting to convert, to a 360 degree basis, the normal bearings in a survey. The following are the original survey bearings (2nd column) and on the right (4th column) is the 360 basis amount I obtain, which will be entered in Copan map checking. The 3rd column is my interpretation of what I think the survey bearing should be. ` 1 2 S 05°52'13" E S 05°52'13" E 174.0147 3 N 84°41'49" E S 84°41'49" W 264.4149 4 N 07°36'50" W N 07°36'50" W 352.2310 5 S 50°17'08" E N 50°17'08" W 309.4242 6 N 37°13'00" E N 37°13'00" E 37.13 7 N 52°47'00" W S 52°47'00" E 127.1300 1 S 83°07'18" W N 83°07'18" E 83.0718` The survey plat, with original bearings, is shown under my previous question regarding Map Check file format.


title: TileMill: Can't use both background-color and background-image?

text: Very new to TileMill and am working on my first project. I'm using a pattern image to enhance some areas on the map, water being an example. I've added the pattern to the water areas on my layers, and am now trying to match the background to those areas. Having trouble - can someone please tell me what I'm doing wrong? This renders only the pattern in the map preview, but not the background color as well:               Map {        background-color:@Water;       background-image:url(images/water.png);     }      ![screenshot](http://i.stack.imgur.com/M7WgY.png)


title: Lake Maps with changing water levels

text: Lots of people crash boats and get injured or die every year because of underwater hazards. I'm trying to design lake maps where the topography and data points change with the water level. Any idea on where to begin with this?


Is there a way to convert a shapefile into graph with nodes and edges(except networkx). I'm using qgis and I want to convert a road vector layer into a graph, so that I can use it to apply shortest path algorithms (astar). Please help me to find a solution.


title: Selected output format is not supported with this function? Wikimapia API

text: I have been trying with limited luck to instigate the &format=kml with the place.getbyarea function. I have been trying to isolate any additional paramater that might be conflicting with the exporting process. This has produced the same error every time.               http://api.wikimapia.org/?key=PUT_YOUR_KEY&function;=place.getbyarea&coordsby;=latlon&lon;_min=35.37735⪫_min=33.055148&lon;_max=36.333847⪫_max=34.694027&format;=&pack;=&language;=en&page;=1&count;=50      Is there an issue with any of this input?


title: How to send WPS request using OpenLayers

text: I have generated the WPS request below using GeoServer's `WPS Request Builder`. I wish to append more parameters to this request (geometry selection), and then send it either using GET or POST from within my OpenLayers JavaScript code. I can't seem to get information on how this can be done. How can I append geometry selection parameter to the request, and how can I send the request using OpenLayers? Any assistance will highly appreciated.               <?xml version="1.0" encoding="UTF-8"?>     <wps:Execute version="1.0.0" service="WPS" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns="http://www.opengis.net/wps/1.0.0" xmlns:wfs="http://www.opengis.net/wfs" xmlns:wps="http://www.opengis.net/wps/1.0.0" xmlns:ows="http://www.opengis.net/ows/1.1" xmlns:gml="http://www.opengis.net/gml" xmlns:ogc="http://www.opengis.net/ogc" xmlns:wcs="http://www.opengis.net/wcs/1.1.1" xmlns:xlink="http://www.w3.org/1999/xlink" xsi:schemaLocation="http://www.opengis.net/wps/1.0.0 http://schemas.opengis.net/wps/1.0.0/wpsAll.xsd">       <ows:Identifier>gs:Aggregate</ows:Identifier>       <wps:DataInputs>         <wps:Input>           <ows:Identifier>features</ows:Identifier>           <wps:Reference mimeType="text/xml; subtype=wfs-collection/1.0" xlink:href="http://geoserver/wfs" method="POST">             <wps:Body>               <wfs:GetFeature service="WFS" version="1.0.0" outputFormat="GML2">                 <wfs:Query typeName="learning_workspace:pop_layer"/>               </wfs:GetFeature>             </wps:Body>           </wps:Reference>         </wps:Input>         <wps:Input>           <ows:Identifier>aggregationAttribute</ows:Identifier>           <wps:Data>             <wps:LiteralData>pop</wps:LiteralData>           </wps:Data>         </wps:Input>         <wps:Input>           <ows:Identifier>function</ows:Identifier>           <wps:Data>             <wps:LiteralData>Sum</wps:LiteralData>           </wps:Data>         </wps:Input>         <wps:Input>           <ows:Identifier>singlePass</ows:Identifier>           <wps:Data>             <wps:LiteralData>Yes</wps:LiteralData>           </wps:Data>         </wps:Input>       </wps:DataInputs>       <wps:ResponseForm>         <wps:RawDataOutput mimeType="text/xml">           <ows:Identifier>result</ows:Identifier>         </wps:RawDataOutput>       </wps:ResponseForm>     </wps:Execute>      My environment: GeoServer 2.1.3, OpenLayers 2.11, PostGIS 1.5.


title: Apportioning Points at a street intersection (arcgis or mapinfo)

text: I can't seem to find any evidence of this question, but I'm surprised by that, so pardon if this has been asked. If one is mapping crime location reports, by XY coordinates to streets, to create "hot street" segments, how does ArcGIS handle intersections? Which street (the East/west or north/south) and/or which side of the street, gets counts of events at intersections? Is there a default logic to this? Looking at my map some segments get both intersections, some get one, and they all go to east/west streets, not the perpendicular intersecting streets. For the left/right (east/west) issue, where some segments get both intersections and some get one, my guess was it's based on a true north/south/east/west grid, and whether the intersection is to the east or west, but I can't figure out why the perpendicular streets don't seem to get intersection incidents. Maybe a work around is having those counts be specific to the intersection and not to the line segment? But, I'm not sure how to do that. Thanks for any help in understanding how the system apportions these counts.


title: Trouble when georeferencing CAD files in ArcMap

text: I've definitely got at least a couple of grey hairs while trying to convert some dwg. files to shapefiles. I've seen that this subject has been disscused a lot, but hope someone will have the patience to help me. The procedure I'm following is this: 1) I add the dwg files to ArcMap. 2)TOC- right clik-Convert CAD Feature Layer 3)Data Management Tools-Projections- Define Projection. The conversion happened but then I realized things are not that simple, looked over the files again and apparently I have some "lost points"....![enter image description here](http://i.stack.imgur.com/kxqMx.jpg) My data is acctually in that middle point and it looks like this: ![enter image description here](http://i.stack.imgur.com/KyIR7.jpg) So, I definitely have to georeference the data before converting it to shapefile. But how do I do that? I tried to create some control points, but I'm probably not doing it the right way cause nothing happens. Plus I have no idea where those lost points should actually be. ![enter image description here](http://i.stack.imgur.com/gXxxl.jpg) Forgot to mention that the dwg files have no coordinate system. Thanks for help


title: Where to look for oceanographical/environmental data for the Arctic (Baffin Bay/Davis Strait area)

text: I am looking for oceanographical/environmental variables that are available for Baffin Bay/Davis Strait area over the time period of 1995-2011, however, I am not familiar with oceanographical/environmental datasets. Perhaps some of you may be able to point me in the right direction. Which oceanographical/environmental variables/layers are available from this part of the world and for which time period? Where can I access these? Thank you very much for your help.


title: Unions of polygons (Postgis ST_UNION for Geography type)

text: I am looking for a function that will return the intersection of 2 or more polygons (geography type). I am aware of ST_UNION, ST_COLLECT but it works only for geometry type. Any tip will be really appreciated


title: Identify Error with ArcGIS Server 9.3

text: I am receiving an error ("Error executing identify") when I try to perform an identify operation on map running on ArcGIS Server 9.3. All of the inputs appear to be correct when I check them in the ArcGIS Server Directory. The map I am using has over 90 layers. Could that possibly be causing an issue with the identify operation? Thanks, Seth


title: Google maps, download entire province/region from google maps

text: Is there a way to download a region from Google maps to a PDF or JPEG file? I would like to download as many zoom levels as possible so I can zoom in and out without being online. I need this for practical reasons when I will be traveling to remote China and I need detailed maps of certain provinces/counties.I will be browsing these maps on an iPod touch. I only need the street information, so only the "map" part of gmaps.


title: Is there a way in the Google Maps API to perform a union?

text: I'm looking for a way to take the overlap area of multiple polygons and work with that result as a single entity, ideally all in the API.


title: struggling with GeoDjango install: ERROR: could not access file "$libdir/postgis-1.5": No such file

text: I'm installing GeoDjango on Mac OSX. I've followed the official Mac install instructions, which all went smoothly, and am now creating a spatial database template for PostGIS. However, when I try to load the PostGIS SQL routines, I get `ERROR: could not access file "$libdir/postgis-1.5": No such file or directory`:               postgres$ psql -d template_postgis -f $POSTGIS_SQL_PATH/postgis.sql     psql:/usr/local/pgsql/share/contrib/postgis-1.5/postgis.sql:59: ERROR:  could not access file "$libdir/postgis-1.5": No such file or directory     <snip>     psql:/usr/local/pgsql/share/contrib/postgis-1.5/postgis.sql:7785: ERROR:  type "geometry" does not exist      What's going wrong, and how can I fix it? I think I'm running the right version of pg_config, etc:               postgres$ which pg_config     /usr/local/pgsql/bin/pg_config     postgres$ which psql     /usr/local/pgsql/bin/psql     postgres$ pg_config --pkglibdir     /usr/local/pgsql/lib      If I look in `/usr/local/pgsql/lib`, there is a file called `postgis-1.5.so` there, so I don't understand why it's not found. \-----UPDATE------- Does this tell anyone anything?               delirium:~ ap257$ otool -L /usr/local/pgsql/lib/postgis-1.5.so     /usr/local/pgsql/lib/postgis-1.5.so:             /Library/Frameworks/GEOS.framework/Versions/3/GEOS (compatibility version 4.0.0, current version 4.2.0)             /Library/Frameworks/PROJ.framework/Versions/4/PROJ (compatibility version 7.0.0, current version 7.6.0)             /usr/lib/libxml2.2.dylib (compatibility version 10.0.0, current version 10.3.0)             /usr/lib/libz.1.dylib (compatibility version 1.0.0, current version 1.2.3)             /usr/lib/libSystem.B.dylib (compatibility version 1.0.0, current version 125.0.0)             /usr/lib/libicucore.A.dylib (compatibility version 1.0.0, current version 40.0.0)


title: Spatial query without spatial database module

text: I have the following scenario:   * Oracle database without the space module.   * Application in JEE2 with the following architecture: \-- Entity \-- Persistence (DAO / JPA + JDBC) \-- Business (BC) \-- Vision (MB + JSF) I have an entity called EQUIPMENT must have the attribute COORD. Persist this entity in the database is easy, I'll use JPA. When I enter the system, it should use geolocation to determine the position of the equipment and check if its position is a certain radius from what is recorded in the database. Which APIs should I use to do this type of query considering I do not have the space module of the Oracle? Thanks, Luciano


title: What are the rules that the Gdal rasterize function follows to allocate a value to a cell?

text: How does the rasterize tool decide which value to represent when a cell covers an area which inclueds different values? I need it in this case for polygons but would like to know for lines and points as well. What are the rules that it follows.


title: Create a jgw file for jpg image.

text: I converted ecw file into jpg file with IrfanView, to display the jpg as base layer in openlayers. I understood that I need to create a jgw file to display the jpg file. I found this example:                4.23214625853148 - A      0.00000000000000 - D      0.00000000000000 - B     -4.23214625853148 - E      3404018.70881921 - C      5819863.55539414 - F      for F - i know the left up corner coordinates(in Israel TM projection), i need other projection? I did not understand what is the other parameter:   A is the x scale length in meters? and E its the negative?   And about D and B I have no idea. The file is little bit big: about 300 MB, its will be a problem in openlayers?


title: Where would a perfectly vertical vertical tunnel through the earth come out?

text: I am a programmer. I have a point P1 (lat,lon) given by GPS (WGS84) and I want to know the location of a point P2 if I make a straight line from P1 to the other side of the earth. I have read some topics in this forum but Idon't really know what to do. How can I achieve this if it is possible ?


title: Leaftlet and ArcGIS Server REST

text: I am fairly new to all of this and am looking for some guidance. Any information is appreciated! I am trying to use leaftlet to display shapefiles from an ArcGIS Server using the REST URL. I have seen several examples of this (http://geojason.info/leaflet-vector-layers/demos/arcgis-server/) (http://blog.davebouwman.com/2011/08/04/leaflet-lean-mean-javascript-maps/) and have tried to apply the same formatting to my own code but have not been successful in actually displaying the data on the map. Because I'm not extremely familiar with this stuff I'm not sure what my next step should be to solve the problem.


title: Ordinary kriging with three dimensions in R

text: I want to interpolate temperature for a surface defined by three coordinates: longitude, latitude and altitude. I know one option is to use universal kriging on two dimensions and use altitude as a regressor. ¿is there a way to perform ordinary kriging on three dimensions with the gstat package in R? ¿what option (3d ordinary kriging or 2d universal kriging) is likely to produce more accurate interpolation results?


title: Why is the quantile reclassification method yielding classes of vastly different cell counts?

text: I am trying to reclassify many euclidean distance layers (which are based on lines) each into 10 classes of equal area. When I attempt to do this, using the quantile reclassification method, some classes contain twice as many cells as others. I thought that using the quantile method, the cell counts of each class were supposed to be roughly equal. The "Slice" tool does generate classes of equal cell counts when using the "equal area" method, but I need to know the break values, which the slice tool does not give (But reclassify does). How can I reclassify these layers into ten classes of equal area and know the break values of each class? EDIT: I have found that this can be done by first converting the floating point raster to an integer using the floating point tool.


title: Multiple Centroids within a dataset based on distance

text: Is there a way that we can find out the centroids of clusters within a dataset while differentiating one cluster from other based on distance not on any attribute value. The dataset is present in PostGIS database. Take a look below![enter image description here](http://i.stack.imgur.com/k3ZG7.png)


title: Grass extension from within OSGeo QGIS

text: It is not clear to me how the connection between GRASS and QGIS works. I am aware about GRASS plugin, GRASS tools in processing toolbox and GRASS GUI which comes with OSGeo installation. By default I wanted some good smoothing grid/raster function but here is the chain of problems I am going through:   1. GRASS tools or any other in processing toolbox can't be used when using Windows network login - fortunately I could create local account at work.   2. With local login I successfully tried SAGA tool "Gaussian filter" and some others but was not very happy with the result   3. After that I found here question How to smooth a DEM? which sugest r.denoise and looks pretty cool so I wanted to try it but...   4. I was able to install needed addons extension r.denoise but my GRASS GUI which comes with OSGEO for some reason can't load any raster ( _Unable to fetch interface description for command 'd.rast'._ ) or vector ( _Unable to fetch interface description for command 'd.vect'._ )   5. When I wanted to find out if the problem/bug reading data is reported in OSGeo GRASS I found reporting site down with some python errors... And that brings me to question, is it possible to add/use GRASS addons extension somehow from within QGIS (GRASS plugin or processing toolbox)?


title: Labels disappear zooming IN in Qgis 2.2 Composer?

text: I have a simple map with a polygon Layer of 2 Polygons. Labeling is turned on, so is the option of "Show all labels for this layer, even if they collide". In the print composer the labels keep disappearing if i zoom in! Also when I export as image (I tried with .bmp, jpeg and pdf) only one of the two polygon label shows up!?! What the hell is going on? I get that sometimes labels dissapear when u reduce the scale and the labels run out of space, but by zooming in is new to me! Is there an option in Qgis or the Qgis composer that I am overseeing? Please help, this is driving me nuts! Regards, Martin


title: I need some more information on ASTER for time series analysis

text: See my previous post on ASTER I have a question about that I am not sure about. So, if I wanted to analyze 10 years or so, do I need to download each ASTER file for each year (e.g. 2003 ASTER, 2004 ASTER,...) and then merge all of them into one image? Thank you! Revised on September 29. 2013 So an ASTER revisit a specific area twice a month, so do I need to download two ASTER each month for 2003, 2004, 2005, etc. until 2013. Then run ndvi for each ASTER per month ? I do need to perfrom them crop analysis/ or soil crop health for this from 2003 to 2013. Thank you


title: How do I iterate through every cell in a continuous raster?

text: See this link for more details. ## **The Problem:** I want to loop through a continuous raster (one that has no attribute table), cell by cell, and get the value of the cell. I want to take those values and run conditionals on them, emulating the map algebra steps detailed below without actually using the raster calculator. Per request of comments below, I have added details providing background to the problem and justifying the need to implement a method as such in the section below called "The analysis needed:". **The analysis proposed below, while being relevant to my problem by providing background,does not need to be implemented in an answer.** The scope of the question only pertains to iterating through a continuous raster to get/set the cell values. ## **The analysis needed:** If ANY of the following conditions are fulfilled, give output cell a value of 1. Only give output cell a value of 0 if none of the conditions are fulfilled. **Condition 1:** If cell value is greater than top and bottom cells, give value of 1:               Con("raster" > FocalStatistics("raster", NbrIrregular("C:\filepath\kernel_file.txt"), "MAXIMUM"), 1, 0)      Where kernel file looks like this:               3 3      0 1 0     0 0 0     0 1 0      **Condition 2:** If cell value is greater than left and right cells, give value of 1:               Con("raster" > FocalStatistics("raster", NbrIrregular("C:\filepath\kernel_file.txt"), "MAXIMUM"), 1, 0)      Where kernel file looks like this:               3 3      0 0 0     1 0 1     0 0 0        **Condition 3:** If cell value is greater than topleft and bottomright cells, give value of 1:               Con("raster" > FocalStatistics("raster", NbrIrregular("C:\filepath\kernel_file.txt"), "MAXIMUM"), 1, 0)      Where kernel file looks like this:               3 3      1 0 0     0 0 0     0 0 1       **Condition 4:** If cell value is greater than bottomleft and topright cells, give value of 1:               Con("raster" > FocalStatistics("raster", NbrIrregular("C:\filepath\kernel_file.txt"), "MAXIMUM"), 1, 0)      Where kernel file looks like this:               3 3      0 0 1     0 0 0     1 0 0       **Condition 5:** If **any** one of the adjacent cells has a value EQUAL to the center cell, give the output raster a value of 1 (using focal variety with two nearest neighborhood calculations) ## **Why not use map algebra?** It has been noted below that my problem could be solved using map algebra but as seen above this is a grand total of six raster calculations, plus one to combine all of the rasters created together. It seems to me that it is **much** more efficient to go cell-by-cell and do all of the comparisons at once in each cell instead of looping through each individually seven times and utilizing quite a bit of memory to create seven rasters. ## **How should the problem be attacked?** The link above advises to use IPixelBlock interface, however it is unclear from ESRI documentation whether you are actually accessing a single cell value itself through IPixelBlock, or if you are accessing multiple cell values from the size of the IPixelBlock you set. A good answer should suggest a method for accessing the cell values of a continuous raster and provide an explanation of the methodology behind the code, if not apparently obvious. ## **In summary:** What is the best method to loop through every cell in a CONTINUOUS raster (which has **no attribute table** ) to access its cell values? A good answer need not implement the analysis steps described above, it needs only to provide a methodology to access cell values of a raster.


title: Adding Properties to GeoJSON File

text: I have a GeoJSON file of the country of Haiti with the administrative boundaries for communes (i.e. counties). Here is the file: https://rawgit.com/mattbowlby/bf5243ab097e50eb8864/raw/46818c96bc49629681238d5302d597456227c040/haiti_communes.geo.json That file has a few properties, for example:               ( "properties": {             "id_com": 112,             "commune": "Delmas",             "id_dept": 1,             "dept": "West Department",             "shape_leng": 0.242799419022,             "latitude": 18.5543057016,             "longitude": -72.297644435,             "shape_le_1": 0.242799419022,             "shape_area": 0.00222841371629      Now, I want to be able to add a number of new properties, such as area, population and other things. I have 141 different communes and do not want to enter this information manually. Is there a way to work with a GeoJSON file and add the properties, such as with a spreadsheet or other format? Ideally I would have it in a spreadsheet and add in the new properties as columns.


title: How to automate drawing a polyline from 2 intersecting datasets based on proximity?

text: I've been stuck on a situation recently where my company's data has to be integrated with another.   Basically there are two point shapefiles comprising of >50000 points, and as shown in the image below, one data set lies between the other spatially.   The output I need would be around <100 polylines based on a unique identifier within the two datasets.   The problem occurs because once they are merged, the FID number (which I would usually use on the _Line type_ within the points to line tool) no longer has any spatial use as my own dataset is at the top of the attribute table and the other at the bottom (no longer in spatial order). Can anyone think of a solution to this?   A polyline tool based on proximity would be perfect for this situation but I am unaware of one existing... (I've already thought about sorting by northing/easting, but that creates problems when the data does not go in a continuous orientation) ![enter image description here](http://i.stack.imgur.com/hCg6o.png)


title: How to replace a cyrillic symbol in an attribute table using python?

text: I've tried searching for an answer and I've already seen some similar themes here, but I still haven't been successful in replacing Cyrillic symbols in an attribute table, which is why I'm asking for help. The aim is to find a Cyrillic symbol (if it is in a field) and then change it, otherwise leave other non-Cyrillic characters as they were.


title: Error 1120 and 1172 flex viewer 3.4 and flash builder 4.6

text: I'm adding a thematic widget from GitHub to my project flex viewer 3.4 in flash builder 4.6, but when I add the source path of this widget I have this error: > -Type 1120: Accès à la propriété non définie GRLegendSkin. > > -Type 1120: Accès à la propriété non définie UIComponentUtil. > > -Type 1172: La définition widgets.Thematic.skins:GRLegendSkin est > introuvable. > > -Type 1172: La définition widgets.Thematic.skins:GRLegendSkin est > introuvable. > > -Type 1172: La définition widgets.Thematic.utils:UIComponentUtil est > introuvable. Any help?


title: Couldn't load PyQGIS, QGIS 2.0.1 on Ubuntu 12.04 Precise

text: After upgrading to QGIS 2.0.1 I still get the error: ` Couldn't load PyQGIS. Python support will be disabled. Traceback (most recent call last): File "", line 1, in File "/usr/lib/python2.7/dist-packages/qgis/ **init**.py", line 63, in from qgis import core RuntimeError: the sip module implements API v10.0 to v10.1 but the qgis.core module requires API v8.1 ` I have installed the latest python-qgis package (2.0.1-2+precise1). Am I missing something the right (older) SIP package?


title: Remove slivers from merged polygons

text: I am creating a single polygon by merging different polygons in ArcMap. When I merge them, slivers are created. So I have to delete those slivers manually. I don't have a license for ArcInfo so tools like Eliminate and others are not working. Are there any other options for removing slivers all at once in ArcMap? I am working on ArcGIS 10. ![Example of sliver](http://i.stack.imgur.com/wShVo.png)


title: WMS Request not served to OpenLayers by GeoWebCache off GeoServer

text: I'm developing a GIS app that makes use of GeoServer and GeoWebCache for serving requests, and OpenLayers for frontend. WMS requests to the GeoServer are served with no issues. But when I make a request through GeoWebCache, nothing is returned to the client. I've enabled Direct WMS Integration in GeoWebCache configuration. I'm using GeoWebCache integrated in GeoServer. My request is as shown below:               var wmslayer = new OpenLayers.Layer.WMS(        "Kenya Administrative Boundaries",        "http://localhost:8081/geoserver/gwc/service/wms?",         {layers: 'kenya_admin', format:'image/png', transparent: true},        {opacity: 1.0, isBaseLayer: false, visibility: true}     );     map.addLayer(wmsLayer);      But the layer is not loaded. Any advice on what I may be doing wrong. Thanks in advance.


title: Geodetic Coordinates And Latitude, Longitude

text: Would it right to say that geodetic coordinates (phi & lambda) are the same as latitude & longitude?


title: Google earth deviation from arcgis 10 based maps

text: I really need your help.. I am using arcgis 10 and i mapped some points. with the based maps provided, i tried digitizing the road network. When i exported to kmz file and opened it in google eath. there is a shift. is it that i am not doing something right


title: pgsql2shp encoding problem

text: I have a database and exporting with pgsql2shp various shapefiles. But when I open the dbf the characters are displayed wrong. Example: BroeksterwÔld and this should be Broeksterwäld. This probably has to do with encoding. The database is UTF-8. I tried to export the shapefiles as LATIN1, but I can't figure out how to do so. Hope someone can help me. * * * **UPDATE** I solved this problem by adding 'set pgclientencoding=latin1' (without quotes) to the batch script.


title: C# Library for consuming Geoserver REST?

text: I am currently experimenting with Geoserver and it's capabilities within C#. One thing I would like to do is query the server for basic information (existing workspaces, layers, layer attributes, etc...) I noticed that Geoserver does have a REST implementation, but is there a C# implementation for consuming this service? I am sure a library could easily be written, but I would hate to re-invent the wheel if one already exists. Thank you! **Edit** This can be done using C# WebRequests. However, what I would have created is a C# library for consuming the Geoserver REST implementation. It seems like something many people would use. Another way to ask the question would be, "Does a library full of C# WebRequests exist to consume Geoserver's REST implementation?" It is much easier if a library already exists with classes such as these to consume the webservice, than it is to create an entire library. For Example:               var geoserver = new Geoserver("http://localhost:8080/geoserver/web");     var workSpaces = geoserver.Workspaces;     var layer = geoserver.Workspaces[0].Layers["Foo"];      The above code is the type of implementation I would really enjoy. If I was to create a library that consumed the Geoserver REST service. The above is what I would create for myself to use. If it wasn't created already, I would likely make an opensource project that others could use as well. So that in the future, if someone is wanting a library to easily consume the Geoserver service, they can contribute and use the one that already exists! :)


title: Creating a Geodesic Polyline Buffer in OpenLayers

text: I have a feeling that the answer is going to be 'You can't do that without some very complex coding' but I thought I'd go for it anyway. I recently asked a question on how to create a Polyline Buffer in OpenLayers, and it worked! Well, almost. In our mapping application Dorset Explorer we have a line measurement tool which works using Geodesic measurements. It had to be like this as non-geodesic resulted in massive errors (I measured it against the 100m running track at my old school and it returned 125m, which I know was wrong!). But simply setting geodesic to True in my OpenLayers measure code fixed this. I now want to create polyline buffers of a user specified size in my OL Map. I have managed to get it to make a buffer using the answer to my first question, but the size is wrong. In the attached image I requested a 100m buffer around a polyline going from left to right (which is not visible). I then measured the size of the returned buffer (in red) and came up with the result 125 metres. ![Requested 100m buffer, got 125m back :\(](http://i.stack.imgur.com/sOL8e.jpg) Not good. I wonder if anyone else out there has attempted to make Geodesic polyline buffers. Even a calculation that would allow the user to put in a 'Geodesic' length and return a 'non-geodesic' length for the buffering tool to use, but I can imagine that this is near impossible. It's driving me slightly nuts, it's so close to working but this rather major issue has stopped it. If you want any more info to help answer or point me in the right direction please let me know. Thanks!


title: Spidertracks to GIS

text: I have some data that will be comming to me from helicopter flights and was recorded with Spidertracks software. Does anyone know how to convert flight tracks from 'Spidertracks' gps flight tracks into shapefile or kml? What is the raw spidertracks gps format? I am trying to learn before I recieve the data and I'm not 100% sure what format I will be recieving.


title: OCM Landscape data

text: What is the source of the OCM Landscape data and where can one obtain it? Let me know if I should be speaking with someone else. Thank you,


title: Tabulate Area Errors

text: Having some issues with the tabulate area spatial analyst tool. Not sure why it's not working, but the output table is missing zones from the input feature class that defines them. I have a simple raster with 3 values. I'd like to simply find out the amount of cells that lie in each zone of the polygon. The tabulate tool should do this no problem. But, after the tool runs, the output table only has about 1800 records. The original input zone polygon has around 5000. Any ideas on why this tool is failing so badly??


title: Best definition for GIS today?

text: What is the best definition of Geographical Information Systems? Has the definition changed at all since I was taught it 10yrs back? The technology has changed a lot and it is more powerful than it used to be, but has that changed the overall definition of what it is? e.g. Is stuff like Augmented Reality a separate discipline altogether, or does it also merge into the world of GIS


title: QGIS 2.0 Plugin error: Unable to load GdalTools plugin?

text: I just used the osgeo4w exe to install QGIS and all the other software included in the basic setup. When I open QGIS, I get the following error:               Plugin error     Unable to load GdalTools plugin.     The required "osgeo [python-gdal]" module is missing.     Install it and try again.      I opened the Python.exe in the osgeo4w bin folder and am able to import gdal just fine. Any idea what's wrong here?


title: My Coordinate Reference System listing is in the wrong language

text: ![enter image description here](http://i.stack.imgur.com/upUfe.jpg) I just installed QGIS Desktop 2.0.1 on a desktop (Windows 8-64). All the menus are in English-US (as they should be), but the CRS listings -- when I try to change them for a project are a layer -- are in French. I fiddled with changing the language parameters in Settings | Options | Locale, but no matter what I do, the CRS listings still come up in French. You can see in the attached image that the menus are in English, but the CRS names are in French (i.e., "Projection conique conforme" rather than "Conformal conic projection"). Any idea what I need to do to fix this? Dave


title: How to change polygon display/symbol properties in QGIS while digitizing?

text: I have two questions. It would be great if anyone help me to sort this out. Thanks in advance 1) While capturing a polygon in QGIS, how to remove the red - shade - sort of thing..... ? 2)Though the symbol of the polygon has been set as **No brush - Solid line** ..... While using split tool and when we select the polygon..it's get selected as a whole yellow fill. So if there is background layer which i'll be using as the reference layer to split the polygon...i can't see that one. Split tool works even without selecting the particular polygon. But i wanna know..is it possible to set the select properties like..a simple No fill - Red Solid line ? ![enter image description here](http://i.stack.imgur.com/dpJ0A.png)


title: Add feature to replica

text: I create one-way replica from one geodatabase to another. After that, i add new feature class in parent geodatabase. Is it possible to include new feature class in replica, but i don't want to recreate replica?


title: Tutorials for Geomedia

text: Does anyone know good resources to learn GeoMedia Professional? I am looking for resources which can serve from a future teaching perspective, but other kinds of resources are also welcome.


title: Teaching GIS, what to do and what not to do?

text: Because I will (in the near future) be a teacher of GIS. I need some advice form people who have previously done this. Any clues - what to do and what not to do ? I think it will be very helpful for my preparation.


title: Export list of values into csv or txt file

text: I would like to know how to export in Arcgis a list of values calculated in python script into one of the following data formats: csv, txt, xls, dbase or other. I would also like to know how to create such file in case that it doesnt exist. The list of values looks like res=(1,2,3,...,x). Each value must be written into a new row.


title: Feature Layer from JSON file

text: I would like to project my shapefile information on the maps. I converted my shapefile into a JSON file. How can use the Json file to create a featurecollection object in the ARCGIS Javascript API. The Json file contains 50 thousand features. Have given a sample set of my feature collection below. I do not want to embed the entire JSON file within the javascript file.               { "type": "FeatureCollection", "features": [ { "type": "Feature", "properties": { "ROADNAME": "AAA", "LINK_INDEX": 1}, "geometry": { "type": "LineString", "coordinates": [ [ 163.554342, -77.702882 ], [     163.553160, -77.702327 ] ] } } , { "type": "Feature", "properties": { "ROADNAME": "KIAN TECK WAY","LINK_INDEX": 2}, "geometry": { "type": "LineString", "coordinates": [ [ 162.527145, -77.686591 ], [     162.522273, -77.686598 ] ] } } }


title: How do I fix an error starting QGIS 1.8.0 on Windows XP?

text: I note that the installation of qgis completed successfully. However, the log procedure returns this message: > "...... Estrazione: wvhpgn.gsb Estrazione: wyhpgn.gsb Cartella di > destinazione: C:\Programmi\Quantum GIS Lisboa Creato il programma di > disinstallazione: C:\Programmi\Quantum GIS Lisboa\Uninstall-QGIS.exe Crea > cartella: C:\Documents and Settings\All Users\Menu Avvio\Programmi\Quantum > GIS Lisboa **Impossibile accedere al file. Il file Š utilizzato da un altro > processo.** Running postinstall grass64.bat... Running postinstall > msvcrt.bat... ......" > > System parameters: Quantum GIS Lisboa 1.8.0 Microsoft Windows Professional, > Versione 2002, Service Pack 3 The message: > Impossibile accedere al file. Il file Š utilizzato da un altro processo translates as > Unable to access the file. The file used by another process I then get the following error: > "Impossibile trovare il punto di ingresso locale_charset della procedura > nella libreria di collegamento dinamico iconv.dll" which translates as: > Can not find the entry point of locale_charset procedure in the dynamic link > library iconv.dll


title: Assign unique ID to each polygon

text: I feel like this should be really simple and I am missing something. All I wish to do is assign a unique ID to each polygon since I have multiple identical FID values. In MS excel I would simply type 1 then run this value down the rows of the column and it would automatically assign 1,2,3.....n. Can I do something in Arc10? Thanks


title: How do I randomly alter a raster with som correlation

text: I need to randomly alter an existing raster. The original raster is essentially a topography map that measures its values around a normal distribution ranging between 0.5 and 1.5 The new raster needs to be randomly generated but have strong correlation with the max and min values of the original. correlation reduces as you move further from the max and min values and randomness consequently increases. A similar question is posted below, but I don't think it strickly applies How do I randomly alter a raster map of habitat types? currently using arcgis 10, can use R but would rather keep it in house. any help would be greatly appreciated i have a feeling the raster calculator may be the best avenue, but i have had no real success as yet.


title: Connexion ODBC FilemakerPro10/QGIS 2.2.0

text: I try to establish an ODBC connection between a database FilemakerPro 10 and QGIS 2.2.0 software but I have some problems with installation. At first , I activated ODBC sharing my file in my Filemaker database . I work on a 64bits Operating System which requires me to go through the System32 to recover administrator ODBC data sources. I tried to link with DataDirect SequeLink 6.0 but also with ODBC Filemaker, I loaded with win64 driver, the connection test succeeded in both cases. I encounter problems when I pass on QGIS and I make a connection from the addition of a vector or from the connection evis to a database, in the first case the following error message appears: "Connection to [localhost.name of my database] failed: The specified driver was not loaded due to system error 182 (DataDirect SequeLink 6.0, C: \ Windows \ system32 \ DataDirect \ slodbc60 \ ivslk22.dll) QODBC3. Unable to establish a connection. "And in the second case, QGIS stop funtionning. Would you have a solution for me? I just met another problem with QGIS, I georeferenced cadastral maps on ArcGIS and I would like to import them into QGIS so that my colleagues could view them, unfortunately the plans are not displayed georeferenced. Thank you for your help and I apologize for my English. Gaelle


title: Leaflet: Centering map when zoom in

text: For a project I have to use four sets of JSON datas (and one of the files is bigger than 10M). I am loading these files asynchronously using queue.js. Once it's fully loaded and displayed I am meeting two problems. Firstly, there is no animation when I zoom in or out the window. It's jerky and laggy. Secondly when I tried to zoom with the mouse the the zoom is systematically centered in the top-left corner. I self-hosted my code here: http://adrienrahier.com/test/SFmap/ you can found the sources here: https://github.com/FracArt/ChallengeSF


title: What minimum scale(out beyond) for layer 200k do I set?

text: I have layer Zone9GK200k. So, I set: "In beyond 1:200 000". What minimum scale(out beyond) for this layer do I set (It is not 1:500 000 000)? I do not have top-level layers, except the one layer with the earth.


title: How to install ArcGIS 10.2 for Desktop on Windows 7 (64bit)?

text: I am going to install ArcGIS for Desktop 10.2 but I cannot find a tutorial about how to do that! Does anyone have any advice on how to install this software?


title: CAD annotation guide lines do not show in Arcmap

text: When I open a CAD file in Arcmap the annotation's guid lines and arrows linking to feature do not appear. see snapshots below. how can I show thses guides in armap? ![dwg in autocad](http://i.stack.imgur.com/XsobA.jpg) ![dwg in arcmap](http://i.stack.imgur.com/YS9zl.jpg)


title: downloading data from UNdata onto ArcGIS 10.1

text: http://data.un.org/Search.aspx?q=population I am trying to download data from UNdata and successfully opened one link of raster data, however the data that I really desire won't open on Arc. I hit the download button and it gives me the choices of XML, Comma, Semicolon, and Pipe. I can download the XML to a zipfile to the correct folder and drag it onto the Arc page but then it simply says "could not open specified file". It doesn't like the XML file. This data would be perfect for my project, is there any way I could get it to show up? or is it simply not compatible?


title: Decompress a LZW-compressed GeoTIFF

text: I have a LZW-compressed GeoTIFF file and would like to split it into several smaller regional files and change the format to png. I tried to do the spliting with gdal_translate -srcwin. I only get a black image. I guess that I first have to deal with the LZW-compression, but I cannot find a way to decompress files. To compress the file I have the option -co COMPRESS=LZW. But how do I go the other way around. I am very much looking forward to your answers! :) Here is the command I used:               $ gdal_translate -ot UInt16 -srcwin 6500 3300 1000 1000 input.tif output.tif          $ gdal_translate -of PNG -scale output.tif output.png      This is what gdalinfo says:               $ gdalinfo input.tif     Driver: GTiff/GeoTIFF     Files: EUD_CP-DEMS_2500035000-AA.tif     Size is 40000, 40000     Coordinate System is:     PROJCS["ETRS89 / ETRS-LAEA",         GEOGCS["ETRS89",             DATUM["European_Terrestrial_Reference_System_1989",                 SPHEROID["GRS 1980",6378137,298.2572221000027,                     AUTHORITY["EPSG","7019"]],                 AUTHORITY["EPSG","6258"]],             PRIMEM["Greenwich",0],             UNIT["degree",0.0174532925199433],             AUTHORITY["EPSG","4258"]],         PROJECTION["Lambert_Azimuthal_Equal_Area"],         PARAMETER["latitude_of_center",52],         PARAMETER["longitude_of_center",10],          PARAMETER["false_easting",4321000],         PARAMETER["false_northing",3210000],         UNIT["metre",1,             AUTHORITY["EPSG","9001"]],         AUTHORITY["EPSG","3035"]]     Origin = (2000000.000000000000000,4000000.000000000000000)     Pixel Size = (25.000000000000000,-25.000000000000000)     Metadata:       TIFFTAG_DOCUMENTNAME=EUD_CP-DEMS_2500035000-AA.tif       TIFFTAG_IMAGEDESCRIPTION=File written by egcs_wrgtif 2.1       TIFFTAG_SOFTWARE=IDL 8.2, Exelis Visual Information Solutions, Inc.       TIFFTAG_DATETIME=2013:09:17 18:46:37       TIFFTAG_XRESOLUTION=100       TIFFTAG_YRESOLUTION=100       TIFFTAG_RESOLUTIONUNIT=2 (pixels/inch)       AREA_OR_POINT=Area     Image Structure Metadata:       COMPRESSION=LZW       INTERLEAVE=BAND     Corner Coordinates:     Upper Left  ( 2000000.000, 4000000.000) ( 26d54'5.12"W, 53d35'57.62"N)     Lower Left  ( 2000000.000, 3000000.000) ( 20d45'24.21"W, 45d41'42.74"N)     Upper Right ( 3000000.000, 4000000.000) ( 12d15'58.58"W, 57d13'49.19"N)     Lower Right ( 3000000.000, 3000000.000) (  8d 7'39.52"W, 48d38'23.47"N)     Center      ( 2500000.000, 3500000.000) ( 16d56'50.88"W, 51d31'45.31"N)     Band 1 Block=40000x1 Type=Float32, ColorInterp=Gray


title: How to save a QGIS 2.0 project template?

text: When I attended the QGIS in Action workshop at FOSS4G13 Tim Sutton mentioned there was a new option to save a project as a template which when you opened up a new QGIS project you could pick to choose the saved template. In my QGIS 2.0.1 Windows from the Project menu option there is a "New from Template" but how can I save a project as a template? When I choose Save As... it is just a standard project file. This may be linked with Print Composer but I have not found any guidance on using this option. Thanks


title: How to calculate the distance between two point on the map using ITM coords?

text: good morning, would appriciate any help on how to calculate calculate the distance between two point on the map using ITM coords? i have a map with items on it all with ITM coords and i need to return all items in X meters of radius from my location. the problem is how to convert the distance between the two points on the map to meters or km. thanks in advance.


title: Is there an ArcGIS 10 Amazon Machine Image (AMI) patched with SP4?

text: I know you can patch the SP3 AMI with SP4, but is there a SP4 AMI for ArcGIS server 10 and if it exists, what is the AMI Name (or AMI ID). I know Esri has to specifically share an AMI with my Amazon account before it is discoverable in an AMI search. But you wouldn't know to ask for it.


title: How to use -lco options in gdal_contour?

text: I'm trying to create contours from geotiff raster with gdal_contour. My output format is Postgresql and I would like to ovewrite the existing table. When I do:   **gdal_contour -a 1 -f PostgreSQL -i 50 -nln new_isolines dem.tif "PG:host=localhost port=5432 dbname=mybase user=postgres password=postgres"**   I get:   **ERROR 1: Layer new_isolines already exists, CreateLayer failed. Use the layer creation option OVERWRITE=YES to replace it.** But when I do:   **gdal_contour -a 1 -f PostgreSQL -lco OVERWRITE=YES -i 50 -nln new_isolines dem.tif "PG:host=localhost port=5432 dbname=mybase user=postgres password=postgres"**   I get a wrong usage message. I know I can do it with ALTER/DROP table in my db prior to contouring with gdal, but that's unconvenient, because I want to set some more layer creation options:   **OVERWRITE=YES**   **GEOMETRY_NAME=wgsgeom**   **FID=id**   So how do I set theese options in gdal_contour utility?


title: How to stop features to redraw unwillingly in OpenLayers

text: I have a problem where when editing features from a wfs vector layer. After some changes are made to the geometry, the feature redraws itself and the changes are lost. Usually the problem occurs more often if a feature is huge and has many vertices. Also if the user pans and zooms in/out a lot. My guess is that the browser (FireFox 18) "looses the cache" or something and therefore does a new request to the server and the feature is rendered and the changes are lost. Is there a way to solve this or is the solution to "save often"? Br Mike


title: Legend title in print composer with different font styles

text: I want to design a layout with print composer where the legend Title has two different font styles in one e.g. " I had it once but can´t _remember_ ". I would be grateful for help.


title: How to create list from table

text: I have performed a spatial join on some data and now I would like to take that data and create a list. I'm trying to do this dynamically so that it can be easily repeatable via a script or model. I have two fields in my data that I care about: Station Name and Page Number. The overall purpose of my project is to take a 632 page mapbook and divide it into smaller books based on the Station Name Boundary. Exporting to pdf only the pages that are relevant to each station. My spatial join table has given me One_To_Many list of each page and what station boundary it falls within. Now I'd like to create a list that is something like:               Station1     Pages: 1,2,3,4,5,6...          Station2     Pages: 8,9,10,11,12...      and so on. For what its worth I'm using 9.3.1 and the maplogic extension for the mapbooks.


title: How to convert point data into a gridded dataframe for histogram analysis?

text: I am very new at using GIS data and only modestly experienced with R. I've been reading about how to analyze spatial data using the spatial-analyst.net PDF book, so I'm not completely lost, but I thought I could describe my problem and people might suggest ideas. I have a dataset with about 2000 measurements at different lat/long coordinates, although I will probably subdivide this dataset as the data was collected over 3 years and conditions changed over time. Let's call the variable being measured "IP." I want to create a map of IP in the full area under question using Kriging or some other interpolation method on the sample data. Then I want to create a histogram measuring the amount of land in various IP buckets. I will also need to create a histogram which shows the number of samples in each bucket (note a sample could have a higher or lower actual IP than what kriging predicts for its land). I follow how to load the data into a SpatialPointsDataFrame and run a kriging analysis, **where I'm having trouble is how to convert that data into a gridded dataframe so I can do the histogram analysis.** Any suggestions for converting points into grids? I figure this is something that is easy if you know what you're doing, but unfortunately I don't! Any other suggestions would be much appreciated. Thanks, Zach


title: Is there any extension to "optimize" the creation of maps in QGIS?

text: is there any extension to "optimize" the creation of maps in QGIS? i searched the repository, but found nothing useful, maybe i searched in a wrong way.. there are some things i am really missing:   * rulers to place the map-elements,   * the possibility to use thousand-seperator,   * the possibility of changing the elements propertiers by right-click-menu or at least in an easier way as at the moment would be pretty nice. thx for any answer kreike


title: Efficient algorithms to split and join lines?

text: Referring to the Figure below:   1. How to efficiently join segments into lines?    2. How to efficiently break lines into segments?  The above will be applied million times. ![enter image description here](http://i.stack.imgur.com/Vv6Fs.png)


title: Migration Issue With ArcGIS 10 in ExportOperationClass

text: I am migrating my Application built in Version 9.3.1 to version 10 In my code I am getting execption:               {"Retrieving the COM class factory for component with CLSID {4147A362-EFC4-11D3-      A0A4-00C04F6BC626} failed due to the following error: 80040111."}      at this place               IExportOperation pExportOp = new ExportOperationClass();      I have alreday bind my application with bindmanager code               ESRI.ArcGIS.RuntimeManager.Bind(ESRI.ArcGIS.ProductCode.Engine);      Every thing works fine like creating IWorkspaceFactory InitiallizeLicense, exception arise as this statment is executed.               IExportOperation pExportOp = new ExportOperationClass();      I found at google that beside binding , some times we have to use               ESRI.ArcGIS.VersionManager.LoadVersion(ESRI.ArcGIS.ProductCode.Engine, "10.0");      but I am not able to access ESRI.ArcGIS.VersionManager Is there any solution to avoid exception. Please help. Thanks


title: Best source for high-resolution satellite images free/low-cost?

text: I am doing a household survey in Kenya, and I need to ensure that my survey team visits every house within a given area. I've been going to google maps and right-clicking on roofs to get coordinates via the "what's here" option. This is really an inefficient use of time. I'd prefer to have a high-resolution image that I could upload to google earth engine, and then perform a machine-learning algorithm to identify roofs. I'd want to get a text file of the lat/longs of all the houses, which I could then cluster by distance and give to my surveyors, both in map and in checklist format. But the images in google maps satellite view are copyrighted and not available in earth engine. So I need images at high-enough resolution to identify corrugated iron and straw thatch roofs, with roughly ~3km^2 spatial extent. What are some good sources for such data? Ideally that will give it freely or cheaply to graduate students and/or nonprofits?


title: Pdf layers to Geodatabase

text: We had a drawing as a PDF , with Geolocation and layer information available. We tried to export that as an Xml ( current tree) option in Adode reader and tried to import it in a new File geodatabase. The error we encountered was there is no workspace associated with the Xml. What would be the ideal procedure to convert a Pdf to a Geodatabase ?


title: How to show Label of a new feature Layer in Openlayers...?

text:
I have created a dynamic new feature layer along with a marker using openlayers with popups, but, I want to show labels for each new feature created... Below is the code I used...               //To add new Vector Layer     var style = new OpenLayers.Style({         fillColor: "#ffcc66",         strokeColor: "#ff9933",         strokeWidth: 2,         label: "${type,area}",         fontColor: "#333333",         fontFamily: "sans-serif",         fontWeight: "bold"           });                    var type='MAHAMAYATALA';         var vectorLayer = new OpenLayers.Layer.Vector("SupplyOff",{styleMap: new OpenLayers.StyleMap(style)});          //  For Consumer Query         $qry_xy="select * from ci_master@ltdvlp  where  x is not null  AND GIS_SOURCE_CODE LIKE '01%' AND ROWNUM < 5";          $parsed =oci_parse($db_conn_prod,$qry_xy);     oci_execute($parsed);     $nrows3=oci_fetch_all($parsed,$results3);     echo $nrows3;     for($i=0;$i<$nrows3;$i++){     $x=$results3["X"][$i];     $y=$results3["Y"][$i];     $consno=$results3["CONS_NUM"][$i];     $source=$results3["GIS_SOURCE_CODE"][$i];       ?>                            var x = '<?php echo $x;  ?>';             var y ='<?php echo $y;  ?>';             var cons ='<?php echo $consno; ?>';             var type = '<?php echo $source; ?>';                  //alert(cons);                        // Define markers as "features" of the vector layer:     var feature = new OpenLayers.Feature.Vector(             new OpenLayers.Geometry.Point( x, y ),             {description:'cons: '+cons} ,             {externalGraphic: 'img/marker-gold.png', graphicHeight: 25, graphicWidth: 21, graphicXOffset:-12, graphicYOffset:-25  }         );                    vectorLayer.addFeatures(feature);               map.addLayer(vectorLayer);          // for Pop ups Window          var controls = {           selector: new OpenLayers.Control.SelectFeature(vectorLayer, { onSelect: createPopup, onUnselect: destroyPopup })         };               function createPopup(feature) {       feature.popup = new OpenLayers.Popup.FramedCloud("pop",           feature.geometry.getBounds().getCenterLonLat(),           null,           '
'+feature.attributes.description+'
',           null,           true,           function() { controls['selector'].unselectAll(); }       );       //feature.popup.closeOnMove = true;       map.addPopup(feature.popup);     }          function destroyPopup(feature) {       feature.popup.destroy();       feature.popup = null;     }          map.addControl(controls['selector']);     controls['selector'].activate();      Thanks in advance...


title: How do you ignore one-ways osm2po?

text: I've been using OSM2PO and it's a great system so far. However, in version 4.5.2, I can't seem to ignore the one-way restrictions for calculating walk distances. In the browser version, I check on "Ignore Restrictions" and "Ignore OneWays", but it makes no difference and still obeys the one ways, creating a longer/further walk time than is necessary. I have tried various settings in the config file, such as these, but they don't seem to help:               graph.build.excludeWrongWays = true     graph.support.segments = true     graph.support.edgeflags = true     graph.support.reverse = false      Am I missing some kind of setting that disables oneways? Could I have done something in the build phase which would have prevented this? Thanks for your help!


title: Can I use an Update cursor with a versioned feature class which has attachments in ArcObjects .NET

text: I have a problem with my code (C# ArcObjects 10.1). We are trying to update and insert 10 - 15 rows at a time into a feature class hosted in ArcSDE 10.1 on SQL Server 2008 R2. We have a feature class which is versioned with move edits to base and also has attachments which are photos. I did some proof of concept code a few months back and thought I had everything working - I had tested inserting and updating a versioned feature class and inserting attachments associated with a particular feature. Looking back, I don't think I tried these at the same time - the feature class I tested inserting attachments may not have been versioned. The issue I am facing is that my code was working fine and inserting and updating features, as soon as I added the attachment table - it stopped working and I get this error :               Objects in this class cannot be updated outside an edit session     [SDE_CSO.DBO.DefinitiveList_COPY]      at this line :               IFeatureCursor cursor = featureClass.Update(queryFilter, true);      I am starting an edit session and edit operation before attempting this :               muWorkspaceEdit = (IMultiuserWorkspaceEdit)featureWorkspace;     workspaceEdit = (IWorkspaceEdit)featureWorkspace;               muWorkspaceEdit.StartMultiuserEditing(esriMultiuserEditSessionMode.esriMESMVersioned);     workspaceEdit.StartEditing(true);     workspaceEdit.StartEditOperation();      I did some googleing around and found a few old posts on the ESRI forums about a similar using using ArcPy and ArcObjects which seemed to suggest that using an Update cursor will not work with a table which is involved in a relationship. So my question is.... is this possible with an update cursor - if so how? If not can someone offer an alternative solution?


title: How to get a feature that lies closest to a point in openlayers?

text: Im trying to perform routing with openlayers, geoserver and mysql. In order to get the start and end points, I need to be able to get a feature of the node layer (which has point features) that lies closest to an arbitrary point I click on the map. Does openlayers provide any functionality for this. If not, kindly suggest any alternative methods for achieving this. Thanks in advance. Regards