title: Geoserver Cascading WMS, Wrong Projection

text: I configured a geoserver store with an external WMS. I used the following WMS url http://www2.dmsolutions.ca/cgi- bin/mswms_gmap?VERSION=1.1.0&REQUEST=GetCapabilities Once I saved the store, Geoserver listed the layers on the screen. I clicked the "Publish" link of the "DEMO" layer and I saved it with the default options. When I display the layer preview, the layer has a wrong projection and is repeated many times. Here you can see the image: http://stackoverflow.com/questions/14336018/geoserver-cascading-wms-wrong- projection The projection of the layer seems to be wrong. I tried to force to another projection or to keep the native one but results are not better. Have you any ideas about the problem?


title: Import list of addresses into MapPoint for route planning

text: I have a list of 40+ points for each work day in a month. I need to somehow calculate the distance between each of those points for each day. I have a list of addresses in a database, house number, street, city, zip and state are all in separate columns so I can generate any output using some simple programming. But I am stuck at importing addresses to MapPoint. I have tried import wizard but all it does is plots points on map and does not add to route planning. Maybe there is an alternative to this? Or other way around it? Thanks much.


title: WMS Problem in ArcGIS Server 10.1

text: I have created a WMS service for big tile map. ( I'm using ArcGIS Server 10.1) Because I don't have enough space in C drive, I create another folder in another drive and copied tile image to that folde and registered that folder. And added the ArcGIS Server users to that folder as well. The published service is working fine. The WMTS, WFS and WCS is working fine. But I had problem to see WMS ! ![enter image description here](http://i.stack.imgur.com/3GwRu.png) Does anyone has this problem before ? I did what ever I can but cannot have any success. There is one more information that there is one service on that service which is belong to some time ago and WMS of that service is fine. But when I tried to create new webservice from the tiles on the new or old address, I don't have much luck!


title: Migrating from ArcGIS for Server on Amazon Web Services (AWS) to local machine?

text: We have a ArcGIS for Server instances hosted by Amazon Web Services (AWS) with many volumes and would like to move all of this from the cloud to a local server. We have a number of instances, and have looked into sending the Amazon Machine Images (AMIs) directly from the AWS instance to the physical machine. We have looked into S3 buckets, but I am not entirely sure how this process works. Any insight on how to move map services, and AWS instances and volumes to a local ArcGIS for Server site is greatly appreciated?


title: How to find the minimum pixel value within a region of interest?

text: I am working on a Python script for dark object subtraction, a radiometric correction method for remotely sensed data.To accomplish this, I need to find the minimum band pixel value of a waterbody within a region of interest (ROI) (screenshot). My initial thought was to create a polygon within a Landsat scene and use Zonal Statistics as Table to find the minimum value. However, I cannot easily call this minimum value later on in the script. Another thought was to convert the band to a numpy array and find the min value (excluding 0), although I am uncertain how to subset the array to my ROI. How can I find the minimum value using a region of interest so that I can easily call this value later on in the script? ![enter image description here](http://i.stack.imgur.com/4Nry8.jpg)


title: how to extract lat long from WKT data?

text: i want to extract lat long form WKT data and stores into arary list i have tried one solution which is proposed by a user but it has some problems                   String str = tv.getText().toString();              str = str.replaceAll("\\(", "");         str = str.replaceAll("\\)", "");                   String[] a = str.split(",");         for (String b : a) {             System.out.println("-" + b + "-");             coordinates.add(b);         }      my data is like `POINT(73.0185097465598 33.6810265037665)` `MULTILINESTRING((73.0111429810373 33.6783008666403,73.0116544995556 33.6785708996273),(73.0119205001311 33.6783781002692,73.0118432996287 33.678392400389,73.0117917999517 33.6784137997571,73.0117403002747 33.6784638002642,73.0116544995556 33.6785708996273),(73.0131224998036 33.679001500419,73.0119635003153 33.678392400389,73.0119205001311 33.6783781002692),(73.0131224998036 33.679001500419,73.0136031002029 33.6783443999742),(73.0136031002029 33.6783443999742,73.0147099372139 33.67685138958),(73.0147099372139 33.67685138958,73.0150124997272 33.6770292997624,73.0154158996241 33.6773507003746,73.0157677998441 33.6776577999676,73.016042399737 33.6779721004322,73.0162998999205 33.6783149004124,73.0167443997343 33.6789772996627),(73.0167443997343 33.6789772996627,73.0171146002587 33.6794933999005,73.0173738999856 33.6798920001158),(73.0173738999856 33.6798920001158,73.0175952995835 33.6801933998056,73.017689300321 33.6802989001741),(73.017689300321 33.6802989001741,73.017964400237 33.680586200293,73.0185097465598 33.6810265037665))` `MULTILINESTRING((73.0185097465598 33.6810265037665,73.0186117996948 33.681108899756),(73.0186117996948 33.681108899756,73.0186413001559 33.681120000088,73.0186946002756 33.681164699991),(73.0186946002756 33.681164699991,73.0187495002893 33.681095800231),(73.0187495002893 33.681095800231,73.0193249000242 33.6803637997496),(73.0193249000242 33.6803637997496,73.0195480002401 33.6800851996723),(73.0195480002401 33.6800851996723,73.0197712002807 33.6797852002269),(73.0197712002807 33.6797852002269,73.0200115000308 33.6794780997346),(73.0200115000308 33.6794780997346,73.020199700256 33.6792233002152),(73.020199700256 33.6792233002152,73.0206966997929 33.679470600288),(73.0209198998336 33.6791562998234,73.0206966997929 33.679470600288),(73.0211343998328 33.6788777995708,73.0209198998336 33.6791562998234),(73.0218726002404 33.6792491997908,73.0211343998328 33.6788777995708),(73.0218726002404 33.6792491997908,73.0218305002774 33.6789161997229),(73.0218305002774 33.6789161997229,73.0218210996641 33.6788205999907),(73.0218210996641 33.6788205999907,73.021797999678 33.6786753001251,73.0217789997012 33.6784876997476),(73.0217789997012 33.6784876997476,73.0217188998074 33.6780734000674),(73.0217188998074 33.6780734000674,73.0217095000934 33.6779850003076),(73.0217095000934 33.6779850003076,73.0217035996415 33.6778538998377,73.0217102995907 33.6776591003872),(73.0217102995907 33.6776591003872,73.0217008998767 33.6775563996082),(73.0217008998767 33.6775563996082,73.0217635997104 33.6773682002823,73.0217961003098 33.6773020002871),(73.0217961003098 33.6773020002871,73.0218382002728 33.6772349998954),(73.0214683032023 33.6770434798131,73.0218382002728 33.6772349998954))` is stores values like `[0] = POINT(73.0166738279393 33.6788721326803) MULTILINESTRING((73.0131224998036 33.679001500419; [1] = 73.0119635003153 33.678392400389;` then again `[2] = POINT(73.0166738279393 33.6788721326803) MULTILINESTRING((73.0131224998036 33.679001500419;` and so on i want to draw this data on google map anyone help me my platform is android


title: How to join different line segments between points into one entity in QGIS

text: I have a .shp file of pipe lines. I want to join all the pipes between two valves to be considered as one segment in QGIS. Help would be appreciated. Thanks


title: Problem analysing movement using cost surfaces (involving ArcGIS's Path Distance-tool)

text: I would like to analyse hypothetical movement (on foot) in a landscape as based on energy expenditure, but I have run into some trouble that I am hoping that you could help me with. I have tried to do this using ArcGIS’s Path Distance-tool in Spatial Analyst using cost surfaces I have created, but they result is not what I would have expected. This is what my elevation surface looks like (downloaded from ASTER GDEM): ![Elevation data with white/purpe areas being the highest and green the lowest.](http://i.stack.imgur.com/ZWeEL.png) Based on the elevation data I created a cost surface that is supposed to contain energy expenditure (metabolic rate in Watts) per map unit (m). For this I used this formula: `M = 1.5W + 2.0 (W + L) (L / W)2 + N (W + L) (1.5V2 + 0.35V * abs(G + 6))` Or put in Raster Calculator terms: `(1.5 * 60) + (2.0 * (60 + 3) * Square((3 / 60))) + (1.2 * (60 + 3) * (Square((1.5 * "movementspeed")) + (0.35 * "movementspeed") * Abs(("slopeinpercent" + 6))))` Where M is the metabolic rate in Watts, W is the weight of the modelled individual, L is the individual’s carried load, N is a factor that describes ease of movement in the terrain (for testing purposes set to 1.2), V is the individual’s movement speed and G is the slope in percent. This created a surface with values ranging between 90 and 25000, with the majority of values between 90 and 1000 (which seems about right, the absurdly high values are most likely a result of flawed slope values, which could easily be fixed). Movement speed was calculated using this formula: `V = 6e^(-3.5 * |s + 0.05|` where s is the slope in degrees. Or put in Raster Calculator terms: `6 * Exp( - 3.5 * Abs(Tan("slopeindegrees") + 0.05))` This created a surface with values between 0 and 5.9 km/h, which seems about right and is consistent with what I expected. Now, these surfaces were used as input in the Path Distance-tool; the DEM as the input surface raster (i.e. in_surface_raster), the surface with energy expenditure as a cost raster and the DEM as the vertical raster to allow the tool to calculate whether or not the modelled individual is moving up or down a slope. For testing purposes, two points in the north-western and south- eastern corners of the DEM were used as source data (i.e. in_source_data). The output was this (red are unintuitively the lowest values and blue the highest): ![Energy expenditure with red being the lowest values](http://i.stack.imgur.com/cjqk3.png) My interpretation of the output is that it pretty much ignores differences in elevation, and the differences in value are simply related to differences in distance. I would have expected the surface to follow the flatter areas in the western part of the region and avoided the mountainous eastern parts, which it clearly doesn’t do. But, I’m still rather new to these types of analyses, and would appreciate the interpretations of others. So, is anyone able to point out any flaws in my methodology/formulas that could cause the strange output? Or, is the output expected and I’m simply misunderstanding what I should expect from a path distance analysis? Any comments, suggestions or ideas are much appreciated :).


title: Move TEXT_ELEMENT using arcpy.mapping and data driven pages

text: I'm customizing an automated mapping solution using arcpy.mapping and data driven pages. Using a looping statement, I'm hardcoding the position for the mapping elements (legend, north arrow, inset map, etc. so that for each map, the mapping elements do not interfere with the area of interest. My code works for all mapping elements except for TEXT_ELEMENTS. Text element stays in the first position assigned for all of the maps. Not quite sure why this is the case, cause I have been able to get the north arrow, a neat line, legend, scale bar, and inset map (dataframe) to move position for 16 different maps. Sample code below showing how I'm moving the inset map and one of the text elements is shown below:               import arcpy, os     # Create an output directory variable     outDir = r"C:\Desktop\Test"            # Create a new, empty pdf document in the specified output directory     finalpdf_filename = outDir + r"\FinalMapBook.pdf"     if os.path.exists(finalpdf_filename): # Check to see if file already exists, delete if it does       os.remove(finalpdf_filename)     finalPdf = arcpy.mapping.PDFDocumentCreate(finalpdf_filename)     # Create a Data Driven Pages object from the mxd      mxdPath = r"C:\Desktop\Test\States.mxd"     tempMap = arcpy.mapping.MapDocument(mxdPath)     tempDDP = tempMap.dataDrivenPages     # Create objects for the layout elements that will be moving, e.g., inset data frame, scale text     dataFrame = arcpy.mapping.ListDataFrames(tempMap, "LocatorDF")[0]            for pgIndex in range(1, tempDDP.pageCount + 1, 1):            # Create a name for the pdf file you will create for each page       tempMap.dataDrivenPages.currentPageID = pgIndex       pageName = tempMap.dataDrivenPages.pageRow.Name       temp_filename = r"C:\Desktop\Test\temp_pdfs\Land_Cover_" + \                                 str(pgIndex) + pageName + ".pdf"       if os.path.exists(temp_filename):         os.remove(temp_filename)            #Code to modify mapping elements for each page            # Code for Map1       if (pgIndex == 1):              #Set inset map position         inset = arcpy.mapping.ListDataFrames(tempMap, "LocatorDF")[0]         inset.ElementPositionX = 3.25         inset.ElementPositionY = 0         # Move text elements         elm = arcpy.mapping.ListLayoutElements(tempMap, "TEXT_ELEMENT", "Westport")[0]         elm.ElementPositionX = 0.3125            elm.ElementPositionY = 3.0227            # Code for Map2       if (pgIndex == 2):              #Set inset map position         inset = arcpy.mapping.ListDataFrames(tempMap, "LocatorDF")[0]         inset.ElementPositionX = 0         inset.ElementPositionY = 3.25         # Move text elements         elm = arcpy.mapping.ListLayoutElements(tempMap, "TEXT_ELEMENT", "Westport")[0]         elm.ElementPositionX = 2.3125            elm.ElementPositionY = 2.0227                 # Code for Map3       if (pgIndex == 3):         #Set inset map position         inset = arcpy.mapping.ListDataFrames(tempMap, "LocatorDF")[0]         inset.ElementPositionX = 3.25         inset.ElementPositionY = 3.25         # Move text elements         elm = arcpy.mapping.ListLayoutElements(tempMap, "TEXT_ELEMENT", "Westport")[0]         elm.ElementPositionX = 1.3125            elm.ElementPositionY = 2.0227            # Code to export current page and add it to mapbook       tempDDP.exportToPDF(temp_filename, "RANGE", pgIndex)       finalPdf.appendPages(temp_filename)            arcpy.RefreshActiveView()          # Clean up     del tempMap          # Update the properties of the final pdf     finalPdf.updateDocProperties(pdf_open_view="USE_THUMBS",                                  pdf_layout="SINGLE_PAGE")          # Save your result     finalPdf.saveAndClose()


title: Postgres function with BOX3D error "doesnt start with BOX3D"

text: everybody. I created function search nearest point (pois) by using pg_routing as folowing               CREATE OR REPLACE function find_pois(inn_point text, distance int) RETURNS float[] AS  --point(x y)     $$     DECLARE         source_id INTEGER;         target_id INTEGER;         row_data RECORD;         i integer;         retVal float[];         temp_length float;     BEGIN       -- search edge of source      SELECT into source_id source ,            ST_Distance(the_geom, ST_GeometryFromText(            'POINT(587314.857692619 2325407.8399756)', -1)) AS dist          FROM ways          WHERE the_geom && setsrid(            'BOX3D(587214.857692519 2325407.8399756, 587414.857692619 2325507.8399756)'::box3d, -1)          ORDER BY dist LIMIT 1;                -- get all point in 500m radius      FOR row_data IN SELECT b.gid, x(the_geom), y(the_geom) FROM station b        WHERE ST_DWithin(the_geom, 'POINT(587314.857692619 2325407.8399756)', 500.0)         LOOP                       -- with each point, search shortest part              -- 1. find taget_id                   SELECT into target_id target ,                ST_Distance(the_geom, ST_GeometryFromText(                '''POINT(' || row_data.x || ' ' || row_data.y || ')''', -1)) AS dist              FROM ways              WHERE the_geom && setsrid(                '''BOX3D(' || row_data.x-100 || ' ' || row_data.y-100 || ', ' || row_data.x+100 || ' ' || row_data.y+100 || ')'''::box3d, -1)  --cast((cast(row_data.x as float))-100 as text)             ORDER BY dist LIMIT 1;         -- 2. tim duong di ngan nhat tu diem nguon toi diem dich         SELECT --rt.gid, ST_AsGeoJSON(rt.the_geom) AS geojson,              into temp_length sum(length(rt.the_geom))              --ways.gid              FROM ways,                 (SELECT gid, the_geom FROM dijkstra_sp_delta( 'ways', source_id, target_id, 100) ) as rt                  WHERE ways.gid=rt.gid;            retVal := temp_length;         END LOOP;           return retVal;               END;     $$ LANGUAGE plpgsql;      Everytime i try excute this function, i have error:               ERROR:  BOX3D parser - doesnt start with BOX3D(     LINE 1: ... || ' ' ||  $4  || ', ' ||  $5  || ' ' ||  $6  || ')'''::box...                                                                  ^     QUERY:  SELECT target , ST_Distance(the_geom, ST_GeometryFromText( '''POINT(' ||  $1  || ' ' ||  $2  || ')''', -1)) AS dist FROM ways WHERE the_geom && setsrid( '''BOX3D(' ||  $3  || ' ' ||  $4  || ', ' ||  $5  || ' ' ||  $6  || ')'''::box3d, -1) ORDER BY dist LIMIT 1     CONTEXT:  PL/pgSQL function "find_pois" line 27 at SQL statement      I guest there are some error on string in BOX3d but cannot fix Please help me. Thank a lot.


title: What are the differences between Anchored, Framed and FramedCloud?

text: From the API document for OpenLayers, I can tell that a `FramedCloud` has more functionalities than `Framed`, which in turn has more functionalities than `Anchored`. I am hoping for someone to describe these differences in a way which a lay person like me can understand. Thanks.


title: Invalid Username When Developing a QGIS Plugin on Mac

text: On a fresh Mac OSX 10.8 install, I've installed Quantum GIS from Kingchaos binaries and it's running fine. I'm trying to develop a plugin for QGIS using Plugin Builder. The skel that Plugin Builder produces should work off-the- shelf, but it does not compile using the makefile. Somewhere python can't find (`pyuic.py`), and it searches for the `/Users/kyngchaos/...` path. This user does not exist on my system:( Here are my ENVS on .profile:               export PATH=/Library/Frameworks/GDAL.framework/Programs:$PATH     export PYTHONPATH=/Applications/QGIS.app/Contents/Resources/python     export DYLD_LIBRARY_PATH="/Applications/QGIS.app/Contents/MacOS/lib:/Applications/Qgis.app/Contents/Frameworks"     export PATH="/Applications/QGIS.app/Contents/MacOS/bin:$PATH"      Here is the ouptput of make:               Arthur@teste:make     pyuic4 -o ui_teste.py ui_teste.ui     /System/Library/Frameworks/Python.framework/Versions/2.7/Resources/Python.app/Contents/MacOS/Python: can't open file '/Users/kyngchaos/Applications/qgis-python27/Applications/QGIS.app/Contents/MacOS/../Resources/python/PyQt4/uic/pyuic.py': [Errno 2] No such file or directory     make: *** [ui_teste.py] Error 2      Thanks for reading.


title: GRASS v.distance - Getting meters instead of degrees

text: I'm trying to measure the distance between points to a coastline. I have a shapefile of the "High Seas" which I got from NOAA's website here. http://www.nws.noaa.gov/geodata/catalog/wsom/html/marinezones.htm I also have a list of points which represent buildings. I've played with different solutions, but I think that GRASS's v.distance is the best solution. I've been able to get my data into GRASS and can plot the map. The problem is that when I run the v.distance tool, I'm getting decimal degrees, but I really want either meters, miles, or kilometers. So the v.distance function seems to be the best option. However, I've been unable to get meters out of this data. In the help guide, for v.distance, I did find this sentence.               In lat-long locations v.distance gives distances (dist and to_along) in meters not in degrees calculated as geodesic distances on a sphere.      Its not exactly clear to me what this means. It sounds like I should set up the location with a EPSG code that uses lat long. I've used a handful to no avail. For example, I tried EPSG 3857, 3786, 4326, all to no avail. I'd would greatly appreciate any help. Additionally, if you think I'm going about this entirely the wrong way, I'm not married to GRASS. I've learned a lot about GIS concepts, so I'm sure I could pick something else up. (As a side note, I did try QGIS's distance to hub that's part of the MMQGIS plugin. It appears that maps to the nearest centroid of a polygon. I'm trying to get the distance to the nearest pixel/border of a polygon. Or the nearest pixel/edge of a line.) Thanks!


title: Why does st_split not work?

text: it seems I am having a problem with ST_Split function. I am trying the following but it does not work. Can anyone spot any problems with the syntax?               select ST_Split(        ST_GeomFromText('LINESTRING(-3.164154576738156 55.926335330437809,-3.16323392764116 55.926283382919387)'),        ST_GeomFromText('LINESTRING(-3.16348554978961 55.92654456162947,-3.163507914241407 55.926112217812459)'));      I have tried some more sql (incl the examples in postgis docs) with st_split and it does not seem to work... Thanks in advance!


title: OpenLayers BBOX sending double requests

text: I've implemented a openlayers bbox strategy on a vector layer. I have my server return a nice jsonp object wrapped in the callback given by the request. The problem is that during normal use of the map I get al lott of errors like:               TypeError: OpenLayers.Protocol.Script.registry.c18 is not a function      where c18 is the number of the request. The map keeps working, but these show up as actual errors in Firebug. This happens very often just playing around for a minute results in a dozen or so of these errors. In the 'NET' tab I can see that it seems to happen every time the map interface issues two requests for data.               GET poi?bbox=-0.597037...cript.registry.c18, 200 OK, localhost, 101.8KB, 127.0.0.1:80, 250ms          GET poi?bbox=-6.831778...cript.registry.c19, 200 OK, localhost, 102.3KB, 127.0.0.1:80, 421ms      I guess the callback corresponding to the first one of the set of two gets deleted by the second request, so when the first bit of data returns, you get the error. So, why does it fire two requests? and how to fix this. My guess is that it has something to do with resizing the map and zooming at the same time and both events result in a request?? Here is how I configured the layer. It's the only layer with a BBOX strategy:               var in_options = {         'internalProjection': config.projection,         'externalProjection': config.displayProjection     };          // add a layer to the map     this.layer = new OpenLayers.Layer.Vector("POI",{         strategies: [             new OpenLayers.Strategy.BBOX({                 resFactor: 2,                 active: false,                 autoActivate: false,                 ratio: 1.3             }),              new OpenLayers.Strategy.Cluster({                 distance: this.clusterDistance,                 threshold: this.clusterThreshold             })         ],         protocol: new OpenLayers.Protocol.Script({             url: 'http://localhost/tests/poi',             callbackKey: 'callback',             format: new OpenLayers.Format.GeoJSON(in_options),         }),         styleMap: new OpenLayers.StyleMap({             "default": style,         })     });     this.layer.strategies[0].calculateBounds = this.myCalculateBounds;     this.layer.strategies[0].activate();      **Edit** I've looked at the layer properties, there are two listeners registered. One for the BBOX strategy and one for the cluster. None of these invoke each other. So I looked at the source of the BBOX strategy. The TriggerRead Method of the OpenLayers.Strategy.BBOX class has one argument: options. If you console log this one, it's just the event object send by the layers eventHandler. It's not an opion object associated with the BBOX instance. That also explains why the setting `noAbort` wasn't working. I changed this line:               if (this.response && !(options && options.noAbort === true)) {      into this:               if (this.response && !(this.noAbort && this.noAbort === true)) {      and ofcourse set the noAbort: true in the layer strategy                new OpenLayers.Strategy.BBOX({             resFactor: 2,             active: false,             autoActivate: false,             ratio: 1.3,             noAbort: true,         })      Did I discover a bug in OpenLayers?


Pretty straightforward question here. At an increasing rate I see crowd sourcing data being mapped to great effect. **Where are the raw data sources for all these maps?** I know they vary but I don't think I have even found a single source. For instance, twitter maps, where is the data?


title: ESRI - Polyline insertPoint issue (JS 2.5 API)

text: I am trying to basically draw a path on an esri map. with circles like so, but when I call Polyline.insertPoint it doesn't show the new segment:   ![path image](http://i.stack.imgur.com/9xbVF.png) I keep adding points as time progresses. My implementation is to have a Graphic object with a Polyline geometry and a SimpleLineSymbol. I call the Polyline's insertPoint method when I need to add a segment. The problem I'm coming across is that when I call insertPoint, the point point is added to the array, **BUT doesn't color in the line**. The interesting thing is that if I remove the graphic from the graphics layer and then add it back on, then the new line segment shows up. This would be sufficient, except I need to have circles drawn on top of the line, so when I remove+add the Polyline it is drawn on top of the circles. So in my current implementation every time I add a line point, I have to:   1. Call Polyline.insertPoint   2. Remove then add the graphic from the graphics layer   3. Remove + add all the circles on the path (because they just got hidden from the above action) The problem with the above steps is that the line and circles flash on and off during this process which is very undesirable. Sorry for the long question, but my question boils down to this. I feel like I am using the graphic objects provided by esri incorrectly because it's so hard to accomplish a seemingly easy task. Do you have any suggestions on an alternative implementation that has a good probability of not having these same graphical artifacts that I'm experiencing?


title: Geoserver - Key Authentication against a MySql Database

text: I am trying to implement a token based security scheme with Geoserver. What we need to do is send the client a license key when they log into our primary application. This license key is generated by another application and is stored in a MySql Database along with the User name and password, Ip Address form login, etc. We would like to be able to pass this License key as part of the WMS Request to ensure that only logged in and validated users obtain access to the data. Has anyone implemented anything similar before, any help, suggestions or code would be appreciated. Edit .. I Think I may have answered my own question, we use Classic ASP Still (Dinasour) and something like the following from here might do the trick.               <%      set objHttp = Server.CreateObject("Msxml2.ServerXMLHTTP")     strURL = Request("url")"      objHttp.open "GET", strURL, False      objHttp.Send     If objHttp.status = 200 Then      Response.Expires = 90      Response.ContentType = Request("mimeType")      Response.BinaryWrite objHttp.responseBody      set objHttp = Nothing      End If %>


title: Watershed extraction using different DEM tiles

text: I would like to extract the watersheds for a particular area using QGIS. However, I only have Aster 30m DEMs to work with, each tile covering one degree square. The problem is that I need to use four different tiles and the results do not all correlate where they meet. How can I resolve this, or am I doing something wrong?


title: Problem Overlaying layers in OpenLayers

text: I have a html page with 2 WMS layers. The first is one Base Layer , & the second is my personal WMS feature. I load it from Geoserver. The problem when I add the 2 layers I do not have the superposition The problem is maybe on the projection. please can you help me this is my code               <script src='http://maps.google.com/maps?file=api&v=2&key=ABQIAAAAjpkAC9ePGem0lIq5XcMiuhR_wWLPFku8Ix9i2SXYRVK3e45q1BQUd_beF8dtzKET_EteAjPdGDwqpQ'></script>       var map;             OpenLayers.ProxyHost = "/cgi-bin/proxy.cgi?url=";             var host = "http://'my_host':8080/geoserver/wms";                  function init() {                      map = new OpenLayers.Map('map', {                         controls: [                             new OpenLayers.Control.Navigation(),                             new OpenLayers.Control.PanZoomBar(),                             new OpenLayers.Control.LayerSwitcher({'ascending':false}),                             new OpenLayers.Control.ScaleLine(),                                                 new OpenLayers.Control.MousePosition()]});                      var reg = new OpenLayers.Layer.WMS(                     "Région",                     host,                     {layers: 'cite:matregion_view', format: 'image/gif',transparent: true}                 );                              var gmap = new OpenLayers.Layer.Google(                     "Google Streets", // the default                     {numZoomLevels: 20}                 );              map.addLayers([gmap , reg]);     //map.setCenter(new OpenLayers.LonLat(146.7, -41.8), 6);             map.zoomToExtent(new OpenLayers.Bounds(-16.884,20.748,-1.026,35.924));                  }      ![Image](http://i.stack.imgur.com/LaAxL.jpg)


title: Is there a better method to fix empty pixels after classification?

text: empty pixels occurs after classification - in arcmap you can use Generalization tools such as majority filter. However is there a script you could write to better fix it? for example - some of the tools are bad at fixing empty pixels on borders. Ok to be more specific. Is it possible to use a python script to select empty pixels and reclassify them based on distance to neighbours? The nibble tool selects No Data pixels and reclassifies them according to distance to neighbours. But I have a question asking if there is a way to better do it - as with current methods often areas such as borders are missed out (due to lack of neighbours) This is for Arcmap10


title: Are there any intelligent travelling salesmen?

text: Jokes aside, I have had a routing problem that is _almost_ a travelling salesman problem (TSP):   * the starting point is defined   * the end point coincides with the starting point   * each node has to be visited   * the total cost should be minimised Two years ago I thought TSP would be a perfect match, so I ran some sample data through `tsp_solve` and Concorde. Luckily, it was quickly obvious that the TSP shortest path is not the real shortest path, since the problem is made easier by unrealistically requiring the nodes be visited **exactly once**. This picture is just a one-step manual attempt at optimisation of the computed solution and already it saves about the distance of the longest used edge. The problem resurfaced, as I'm trying to find optimal routes to subsets of mapping/monitoring sites. Location and road network data is both pretty accurate and precise, so an exercise like this makes sense. I've looked at generalisations of the TSP, but didn't find an appropriate algorithm. Minimum spanning trees don't account for returning from branches (the 1st solution here costs 3 more). From what I understand, the shortest path problem eventually only cares about two nodes and those out of the optimal path would be left out. A special case of the vehicle routing problem seems to fit best, though I don't know if it considers non-direct paths. My question: **is there any settled name, definition for this kind of problem (family)? What algorithm and tool would you use to solve it?** I'm sure it would be computationally heavy, but I'm interested in both general (infinite resources) and practical answers.


title: Keep leading zeros when importing .csv into QGIS?

text: I'm having issues getting QGIS to retain leading zeros when I import a .csv file. I have zip codes I am trying to read in and I need the leading zeros. I am not super familiar with GIS, so please dumb your answer down a bit for me. Any ideas?


title: Performing image rectification to a jpeg using ERDAS Imagine

text: I am trying to rectify some images on ERDAS Imagine. These images are jpegs. They have no projection data. How do I do this?


title: Arcpy: Calculate new field from user input and existing data

text: Can't seem to figure out what is wrong with this create fields from a calculation. The code creates the new field in the shapefile as selected by the user but I keep getting a 9999 error on every attempt I have to go further. I am rather new to Python and ArcPy, if anyone has any suggestions they would be greatly appreciated. The code is to work as a tool from an exported route from a network analysis and provide a fuel cost which changes slightly in relation to the speed limit of the road. Thanks.                       import arcpy             from arcpy import env             import os                  inputFC = arcpy.GetParameterAsText(0)      # Input feature class             distanceField = arcpy.GetParameterAsText(1)      # Name of distance field in                           input fc             speedField = arcpy.GetParameterAsText(2)      # Name of speed field in input feature class             price = float(arcpy.GetParameterAsText(3))      # Input the current price of gas in $/Liter             mpg = float(arcpy.GetParameterAsText(4))      # Input vehicle mpg in MPG                  arcpy.AddField_management(inputFC, 'Cost', 'DOUBLE') #create cost field             rows = arcpy.SearchCursor(inputFC)             #row = rows.NextRow()             #arcpy.AddMessage(row.GetValue(distanceField))             for row in rows:                 arcpy.AddMessage(row.length_12)                 distance = row.getValue(distanceField)                 arcpy.AddMessage(distance)                 lperkm = row.getValue(mpg)                 milageupdate = lperkm * 2.35214  #convert mpg to liters per kilometer                 if speedField <= 30:                     cost =(price * 0.07 * distance)/ milageupdate -(milageupdate * (0.07))                 elif speedField >30 and speedField <50:                     cost =(price  * distance)/ milageupdate -(milageupdate * (0.03))                 elif speedField >50 and speedField <90:                     cost =(price  * distance)/ milageupdate                 row.setValue('Cost', cost)                 rows.updateRow(row)             del row, rows


title: Variogram model fit compatability among geoR, gstat, and nlme packages in R

text: I'm trying to specify the covariance structure parameters in a linear mixed model (using the correlation structure facilities in `nlme`). The parameters are estimated from a `gstat` or `geoR` variogram model fit to the empirical semivariogram. My justification for specifying the `gstat`-derived covariance model in `nlme` is because `nlme` uses only the classical semivariance estimator fit to the full data set. `gstat` and `geoR` permit the fit of variogram models to the robust (Cressie) semivariogram and the setting of a maximum distance of influence. For example, see my question seeking guidance on the selection between the classical vs. robust semivariogram. I can produce nearly identical empirical semivariograms to the residuals from an intercept-only model (first figure below). Likewise, the exponential variogram model fits from `gstat` and `geoR` are essentially identical. However, fitting the exponential variogram model in `nlme` produces bizarre parameter estimates. `nlme` uses a slightly different specification of the variogram models and correlation structures compared to `gstat` and `geoR`: namely, `nlme` standardizes the within-group errors to unit variance and uses a multiplicative rather than additive nugget effect (see, e.g., p. 230-232 in Pinheiro and Bates, 2000, Mixed-Effects Models in S and S-Plus). Nonetheless, in the absence of a nugget effect (the example below) the exponential fits should be easily translatable among the programs. Here is a simple illustration of what is perplexing me. This example uses the data set of heavy metal concentrations along the river Meuse in The Netherlands (`?meuse` in `gstat` package); the autocorrelation is spatial (heavy metal concentrations closer to each other in space are more similar on average than concentrations farther apart). The empirical variograms of the model residuals using `nlme`, `gstat` and `geoR` are essentially identical.               # Load necessary packages     library(gstat) # variogram model fitting     library(sp) # classes and methods for spatial data     library(geoR) # variogram model fitting     library(nlme) # linear mixed models with autocorrelation     library(ggplot2) #graphics          # Use heavy metal data from `gstat`     # Focus on zinc concentrations     data(meuse)     meuse$dGroup <- 1 # dummy grouping variable so we can use `lme` function     m1 <- lme(log(zinc) ~ 1 , random = ~ 1|dGroup, data=meuse)          # Variogram in `nlme`     lVar <- Variogram(m1, form = ~ x + y, maxDist=1700, breaks=seq(100, 1700, 100))     #lVar; plot(lVar)          # Variogram in `gstat`     meuse2 <- data.frame(x = meuse$x, y = meuse$y, resid = resid(m1, type='n'))     gmeuse <- meuse2     coordinates(gmeuse) = ~x+y     v <- variogram(resid~1, gmeuse, cutoff=1700, width=100)     # Convert to `variogram` class for simultaneous plotting with `nlme` variogram     vdf <- with(v, data.frame(variog = gamma, dist = dist, n.pairs = np))      class(vdf) <- c("Variogram", "data.frame")     #vdf; plot(vdf)          # Variogram in `geoR`     geodat <- as.geodata(meuse2)     geodat.v1 <- variog(geodat, breaks = seq(100, 1700, 100), max.dist=1700, option='bin')     # Convert to `variogram` class for simultaneous plotting with `nlme` variogram     gVar <- with(geodat.v1, data.frame(variog = v, dist = u, n.pairs = n))     class(gVar) <- c("Variogram", "data.frame")     #gVar; plot(gVar)          # Consolidate empirical variograms and plot     # Essentially identical empirical variograms     # Fit is a loess smooth, not a variogram model     empirVar <- data.frame(rbind(gVar, vdf, lVar),                             model = rep(c('geoR', 'gstat', 'nlme'), each=17))     theme_update(legend.position = 'right')     empirVarPlot <- ggplot(empirVar, aes(x=dist, y=variog, colour = model)) +          geom_point() + geom_smooth(fill='white', alpha=0) +          xlab('Distance') + ylab('Semivariogram') +          scale_colour_manual(values=c('red', 'blue', 'green')) + theme_bw()     empirVarPlot      ![empirical variogram comparison](http://i.imgur.com/wdCoO8W.png) Estimating the model fit (e.g., an exponential covariance structure) is also trivial in all of these packages. `geoR` and `gstat` produce nearly identical exponential model fits, but the variogram fit from `nlme` is completely different in both fit and estimated semivariogram values.               # Estimate variogram model in `geoR`     geoExp <- variofit(geodat.v1, ini=c(1, 300), nugget=0.1, cov.model='exponential')     geoExp          > variofit: model parameters estimated by WLS (weighted least squares):     > covariance model is: exponential     > parameter estimates:     >    tausq  sigmasq      phi      >   0.0000   1.2457 338.3446           # Estimate variogram model in `gstat`     gstExp <- fit.variogram(v, vgm(1, "Exp", 300, 1), fit.method=1)     gstExp          >   model    psill    range     > 1   Nug 0.000000   0.0000     > 2   Exp 1.246994 340.6168          # Estimate variogram model in `nlme`     # Note nugget value (second value in vector) is specified differently in `nlme`     # In particular it is 'one minus the correlation between two observations taken      # arbitrarily close together' (see ?corExp and Pinheiro and Bates 2000)     m2 <- update( m1, corr = corExp(c(300, 0.7), form = ~ x + y, nugget = T) )     m2$modelStruct$corStruct          > Correlation structure of class corExp representing     >        range       nugget      > 1.051978e+08 4.192519e-07           # Plot them all     plot(geodat.v1, main = 'geoR'); lines(geoExp, col='red')     plot(v, model = gstExp, main = 'gstat')     plot(Variogram(m2, maxDist = 1700, breaks=seq(100, 1700, 100)), main = 'nlme')      ![geoR exponential](http://i.imgur.com/TZEEUIz.png) ![gstat exponential](http://i.imgur.com/ogbx3Be.png) ![nlme exponential](http://i.imgur.com/iWzPR7U.png) Fixing the values of the `nlme` correlation structure to match the `gstat` model seems like an improvement, but the fit doesn't seem to coincide with the semivariogram estimates. The exponential fit is asymptoting to 1, which makes sense given the specification of the correlation structures in `nlme`, but I'm not sure how to modify this behavior or 'translate' the `gstat` or `geoR` parameters into `nlme` terms.               # Fix the correlation structure to match `gstat` variogram model fit     m3 <- update(m1, corr = corExp(gstExp[2, 3], form = ~ x + y,                                     nugget=FALSE, fixed=TRUE))     m3$modelStruct$corStruct          > Correlation structure of class corExp representing     >    range      > 340.6168           plot(Variogram(m3, maxDist = 1700, breaks=seq(100, 1700, 100)), main = 'nlme')      ![nlme with exponential range fixed based on gstat](http://i.imgur.com/2t1qnGy.png) My questions, to summarize:   1. Why the large discrepancy between the exponential variogram model fits between `gstat`/`geoR` and `nlme`?   2. How do I translate a `gstat`/`geoR` variogram fit into `nlme` given the different specifications of the correlation structure and nugget effects?


title: How to dissolve based on the overlap of polygons in the input file?

text: I am using ArcGIS 9.3/10 to dissolve few hundreds of overlapping polygon features. Now, I am seeing that Arcgis is not able to do dissolve for 5000 polygons. Generally, I used to dissolve (unchecked create multipart features) in arcgis based on geometry not with any dissolve fields. I am trying with QGIS version -1.8.0-Lisboa dissolve tool for that 5000 polygons, here it is mandatory to select any one of the dissolve field. note: there is no common field to group the polygons. I selected dissolve, but it generated ONE polygon (with multiparts).. I need multiple dissolved polygons based on the overlap of the polygons in the input file. How to do this in QGIS or any other open source available, (since arcgis is not able to do)


title: How to store raster data outside postgis database

text: I am able to import raster data easily into the database and view it through QGIS. However, I am unable to view the raster layer in Arcgis. Now I want to know how to create a table with a reference to raster files outside the datbaase. However, there isn't much discuss or example with regards to how it can be done. I want to see its benefits given that raster layers can be shared across various applications. Please share your knowledge on this? Thanks in advance.


title: Customizing Identify Dialog in ArcMap

text: Any idea for customizing Identify Dialog in ArcMap.I want to use my control (for example a list view or a custom control) to show attributes of features in same native ArcMap Identify Dialog,mean replacing some controls of that dialog box.If ArcObject doesn't expose any interface for this purpose,is there any way to replace whole dialog with another dialog but still using same tool,I mean using ArcMap native identify tool,we can show a complete different dialog box which simulates same functionalities. Any idea,help or guide line would be highly appreciated.


title: What is a good resource for learning to hand draw a map?

text: I've been wanting to go old school & draw some maps by hand... Something like you might see on the first page of an adventure novel. I've found next to nothing on the web and am hoping that one of you has (or has seen) an ancient cartography text that describes the basics. If yes, please reply with author/title so that I can narrow my search.


title: How to create a function by passing polygon datatype as a parameter

text: I am getting stuck while executing the function,Function is getting created but while i am passing the values getting error.               CREATE OR REPLACE FUNCTION INS_CCX_ply(P_ID INTEGER,P_POLY_ID INTEGER,P_POLY POLYGON)     RETURNS voiD AS      $BODY$     BEGIN     INSERT INTO CCX_poly(id,poly_id,poly) values(P_ID,P_POLY_ID,P_POLY);     END;        $BODY$     LANGUAGE plpgsql;            SELECT * FROM INS_CCX_PLY(1,11,'point((1,1),(8,6),(6,8),(10,9),(1,1))');


title: An example Tissot ellipse for an Equirectangular Projection?

text: I am trying to calculate distortion of an equirectangular projection via Tissot indicatrices. I have tried following the directions on this post (among other things) but it is beyond my understanding as an amateur. So, I am wondering if someone would be so kind as to calculate a single Tissot ellipse for an example equirectangular lat/long (whichever is your favorite and is distorted on an equirectangular projection). I don't understand exactly what the variables are and where they come from, so seeing the equations in action would be very useful. I am trying to understand these equations so I can plug them into a mapping program I am coding. I asked a bunch of general questions in this thread, but I think a specific example will help me figure the rest out. Many thanks, as always. NCashew


title: Symbolizing Layers in PostGIS

text: I have a number of vectors (points, lines and polygons) loaded into PostGIS (not QGIS Server). I would like to provide a colour schema based on attribute values which is saved as a default. What I am hoping to be able to do is have my layers symbolized (ie: blue for lakes, black for roads and a star for capital cities) so that when users from other computers connect to the PostGIS database the layers are loaded with the symbology. Is there a way to save all the symbols within PostGIS similar to a style file? Thanks....


title: From classified image to topology in ArcInfo

text: I have some coastal 4 band imagery that I have classified to designate land, water, and undetermined. I want to convert this to polygons with topology, and significantly simplify the polygons. If I convert the raster to vector, then try to apply a topology, then simplify the polygons, Arc Crashes. (Raster to Polygon>Select By Attribute x3>Export to FeatureDataset as 3 separate classes>Create Topology in ArcCatalog. Started with only Can't-overlap-self rules, but have also tried Can't overlap each other rules, as well as trying one rule at a time. I have tried validating in Arc catalog, with a validation failed error followed by an ArcCatalog Crash. I have tried not validating until in an edit session in ArcMap, which results in a freeze then a crash. I deleted the topology after each error and started a new one.) If I try to split these polygons, I get the empty output bug in ArcInfo. (Raster to Polygon>Select By Attribute x3>Export to FeatureDataset as 3 separate classes> Extract/Split tool using a fishnet as the split features> Processes Polygons, completes without error, copies everything from dataset of input feature to Target Workspace, doesn't split polygons. Have seen threads from older versions having this bug.) If I convert raster to vector, then simplify, I end up with tons of gaps and overlaps left to fix. (Which defeats the purpose of going to a topology first.) I've got to think this is a fairly common conversion. Is there a specific trick or workflow that I am missing here to accomplish this? Or am I at the mercy of ESRI bugs? ArcMap / ArcInfo with spatial analysis version 10.0 SP5. EDIT:Added specific workflows above. Although I've tried a few other variations using mostly the same tools. The usual save and close between steps. Export as a whole before exporting by attribute.


title: Algorithm to fill the gaps adjacent to parcels

text: I would like to find a way to fill the gaps between the road and its surrounding parcels, like described in picture below: I have different parcels like A and B , i want to fill the gaps by adding the gap area to the adjacent parcel, the gap should be determined by the lines of the parcels which intersects the ROAD. The final parcel B = (old) B + D , A = (old) A + C I would like to find an algorithm to do this process. I thought to make this algo :   * Make a convex hull around ROADS and Parcels.   * Erase Roads and parcels.   * Get lines from the result layer and apply a tool like Extend line, but the constraint is that the lines would not be extended because they are already limited by parcel lines. I'm not keen to follow this algorithm, i'm looking for any suggestions. I'm using Arcpy in Arcgis 10.1, Arcobjects. ![enter image description here](http://i.stack.imgur.com/2IInK.png)


title: Open Layers Bounds with Lat / Long

text: I am using a custom map with the TMS in OpenLayers. I know my high/low latitude and longitude points. In other words I know the exact bounds of my map image in Lat/Long. Is there a way to use lat/long in the bounds? Or do I have to convert to mercator? My bounds are pretty small - it only covers maybe 10 square miles and then I want to be able to draw circles on the map at specific lat long locations - like draw a circle to highlight a gas station. I am just a bit confused about the Bounds portion. here is my code this far - use they mercator bounds just to get my image to display.                   function loadmap() {              var get_my_url = function (bounds) {             var res = this.map.getResolution();             var x = Math.round((bounds.left - this.maxExtent.left) / (res * this.tileSize.w));             var y = Math.round((this.maxExtent.top - bounds.top) / (res * this.tileSize.h));             var z = this.map.getZoom();                  var path = z + "_" + x + "_" + y + "." + this.type;             var url = this.url;             if (url instanceof Array) {                 url = this.selectUrl(path, url);             }             return url + path;         }                   var options = {             maxExtent: new OpenLayers.Bounds(-20037508.3427892, -20037508.3427892, 20037508.3427892, 20037508.3427892),             numZoomLevels: 5,             units: 'm',             projection: "EPSG:900913",             displayProjection: new OpenLayers.Projection("EPSG:4326")         };              var map = new OpenLayers.Map(         'wissmap',         options     );              var layer = new OpenLayers.Layer.TMS(         'Aerial',         '/custom/maps/',         {             type: 'jpg',             getURL: get_my_url         }     );                        map.addLayer(layer);         map.setCenter([2048, 2048], 2);     }


title: Selecting only parallel lines along a center line

text: Read through the previous post how to find close and parallel lines but still looking for a simpler solution. Found several ArcGIS forums post on the topic but no answers there either. It seems likely a tool exists in ET GeoWizards, ArcToolbox, or QGIS that will select only parallel lines (with a given tolerance). Any help on location of this tool would be most appreciated. ![enter image description here](http://i.stack.imgur.com/covPo.png)


title: Polygons Shrinking and Moving as if CRS Was Suddenly Out of Sync with Project

text: I long ago set the CRS for a linguistic geography project I'm working on involving Oklahoma data (EPSG:2836 - NAD83(HARN) / Oklahoma North), and I have been closely monitoring that the polygon shapefiles I have been creating are set to the same. I've been working on the project for several days now, and all seemed fine today when I started working. But after hours of work this evening, I saved my project again only to discover that all of the polygons that I made before today have now moved hundreds of miles away to the southwest (while still maintaining their relative locations with respect to one another) and shrunk to ridiculously small size. For instance, a nearly 220 km wide polygon originally in northeastern Oklahoma is NOW ONLY 3 METERS WIDE and located in Southeastern New Mexico! All of the shifted polygons still have the same CRS settings as the others, have same relative position with respect to one another, and their attribute tables are intact. I reiterate that every polygon I've made so far in this project is its own separate shapefile, most of which overlap with one or more others. Oh, yes, I freely admit that I'm very new to QGIS, and GIS in general, only having started to teach myself the basics this summer. But I would **_GREATLY_** appreciate figuring out what I did wrong and how can I fix it (beyond just scrapping those polygons and starting over again). I don't know of anyone else that uses the program, so there's no one I can turn to.


title: Contour command in Spatial Ecology's GME

text: Has anyone been successful using the contour command from GME (Geospatial Modelling Environment)? I'm trying to create contours out of a DEM. I know you can do it in ArcGIS with 3D or Spatial Analyst extension but I don't have them hence giving GME a try. I've tried it and got the following errors: > Error: The specified workspace does not exist. Cannot open raster dataset: > C:\Users\as\Documents\Columbia Test.gdb!Columbia10M_Elevation > > Error: The command text could not be interpreted. Please check the syntax of > the command. > > Error: An important error has occurred. Please include the information below > if you submit a query about this error. The XY domain on the spatial > reference is not set or invalid. I read the help manual and it says that the algorithm hasn't been widely tested. So I'm thinking that it might not be the best tool for getting contours from DEMs. Any suggestions on a fix or a different work around in lieu of Spatial Analyst? **Update 9/8/14:** GME's developer sent me the following: > "Unfortunately GME is not designed to work with GDB-based rasters. You would > have to export to one of the file-based formats (img, tif, or grid) if you > want to use GME. I will have a look at adding that functionality into a > future release."


title: Combining two overlapping rasters (combine cells, retain both values) using ArcGIS for Desktop?

text: Using ArcGIS 10.2, I am trying to combine two rasters which are identically overlapping in shape and extent, but with different values and counts. As an end result I would like each cell to retain values from both rasters, so I can compare the two values within each cell. Does anyone have an idea of how to accomplish this? Update: Thank you to everyone who provided help. I still haven't been able to accomplish my goal. let me clarify. My end goal is to have a table that does a cell by cell comparison, such as: Value 1: 2, 5, 7 Value 2: 3, 3, 9 Difference: -1, 2, -2 I have spatial analyst, and I have been able to find the difference, but am unable to compare the differences to the original values because the object ID is based on values not individual cells.


title: Problem compiling liblas on OSX 10.8

text: I'm trying to compile liblas since I can't find any Mac executables out there(?). I have OSX 10.8 and XCode 4.6.3. Since I'm not a programmer I get really confused when I follow the instructions on liplas.org . In those instructions older version of XCode is being used so with the new XCode I can't find anything.. I get stuck here And if I just skip all the steps and try compile in XCode I just get this message: > make -f > /Users/oskarlin/Downloads/libLAS-1.7.0b2/xcode/CMakeScripts/ReRunCMake.make > make1: *** No rule to make > target`/Users/oskarlin/Downloads/libLAS-1.7.0b2/xcode/CMakeFiles/2.8.10/CMakeCCompiler.cmake', > needed by`CMakeFiles/cmake.check_cache'. Stop. make: *** > [/Users/oskarlin/Downloads/libLAS-1.7.0b2/xcode/CMakeFiles/ZERO_CHECK] Error > 2 Command /bin/sh failed with exit code 2 Has anyone here complied it for Mac and could perhaps give me an updated instruction?


title: NASA SRTM: precipitous contour lines in TileMill

text: I was trying to make a nice looking print map of Bali (Indonesia) in TileMill. After creating a terrain layer (using this walkthrough) I noticed that the island's coastline was precipitous. ![precipitous coastline in TileMill, Kuta Beach](http://i.stack.imgur.com/yTQeY.jpg) Is this caused by GDAL or wrong input data from CIGAR? I believe it can't be fixed within TileMill as the hillshading layer seems also affected. **Edit** : With a bit CartoCSS the contour lines can easily be hidden. But not the slope scale image... ![faulty slope scale, Southern Bali](http://i.stack.imgur.com/xeKVn.jpg) * * * TileMill 0.10.1 (v0.10.1-134-g31027ed, Windows) GDAL: OSGeo4W on Win 32bit


title: Quering a ArcSDE featureclass attribute table after joining

text: I am trying to query a ArcSDE Feature Class after joining using Iqueryfilter I have two SDE feature classes. 1) Subdivisions : which has details of subdivision and districts (parent) 2) Districts: Which has details of districts and States (parent) I want to find out what are the subdivisions in a state. For this I joined Districts table to Subdivision table and now I can do attribute query in ArcMap and get the divisions in a state using select by attribute tool (example : Districts.statename = 'A' selects all the divisions in the state ‘A’) Now I want to do the same using VBA - I am using the following method but it gives me ORA-00904 error : “Districts”.”statename” : invalid identifier. **Anybody have faced same type of issue erliar? What is best way to get this done. Thank you in advance.**               Dim pQFilter As IQueryFilter          Set pQFilter = New QueryFilter          Dim pFeature As IFeature          pQFilter.WhereClause =  " Districts.statename = 'A'"             Dim myFeatureLayer As IFeatureLayer          Set myFeatureLayer =   “Districts”          Dim pFCursor As IFeatureCursor          Set pFCursor = myFeatureLayer.Search(pQFilter, False)          Set pFeature = pFCursor.NextFeature()


How to achieve the statusbar as the same effect with the arctoolbox‘s tool called in ArcMap add-in? i do not know how to use the progressbar as a result of my project involves in a large mount of data calculate and it is very complicated.So I choose the startbar the same effect with the arctoolbox‘s tool(when the tool was called ,the stausbar show ,while the work is finish ,the statusbar is hidden). I want to achieve the same effect as below (picture link)： ![while Arctoolbox’s tool is called](http://i.stack.imgur.com/xlFSU.png) ![while the tool finish the work](http://i.stack.imgur.com/fnBVs.png)


title: How to host a Geoserver based web page on LAN?

text: I am having an issue hosting a GeoServer based webpage over LAN. I am using Geoserver 2.2 (Binary), OS-Windows. I configured my machine as local server using Tomcat as well as WAMP (Apache) and able to see the webpage with published vector/raster data on port:8090 (Tomcat) and on port:80 (WAMP). But when I tried to access the webpage from another machine connected to LAN, it displays a blank frame with no data. I tried my hands on MapGuide Maestro alternatively but got stuck with same problem. Just to test whether the content I am trying to access over LAN is ok I added a simple picture (.jpeg format) on the hosted webpage. The page displayed the jpeg but the frame again is blank. Please guide and let me know if any futher details are required in this regard. Regards James B.


title: Inpho Match-AT (Application Master - Multi photo measurement)

text: I am new to Inpho Match-AT software and trying to do ground control measurement in Multi photo measurement. I am unable to display the outline in the overview and also the image list is not shown after ticking it to display. Has anyone encounter such situation?


title: Coordinate systems conversion

text: I have the following task. This is my sample data.               WGS84 datum (longitude/latitude):     -123.75 36.59788913307022     -118.125 40.97989806962013          Spherical Mercator (meters):     -13775786.985667605 4383204.9499851465     -13149614.849955441 5009377.085697312          Pixels     2560 6144 2816 6400          Google     x:10, y:24, z:6          TMS     x:10, y:39, z:6          QuadTree     023010      I'm using the code from this source http://www.maptiler.org/google-maps-coordinates-tile-bounds-projection/. How should I chain the methods in order to get from google tiles `(x:10, y:24, z:6)` the Sperical Mercator meters               -13775786.985667605 4383204.9499851465     -13149614.849955441 5009377.085697312 ? I have used lots of combinations but I'm getting wrong numbers. Thank you in advance.


title: Messages in JSON output for a Geoprocessing Service from a Model

text: I have a Model created from Model Builder. It was published as a Geoprocessing Service on ArcGIS Server 10.1 When it was triggered from the REST service, the JSON output printed out the "results" and "messages". The "results" can be configured as derived parameters in Model Builder. However, I couldn't find anywhere that "messages" can be configured. Is it possible? Thanks Cheers, Alex


title: How to transform KML from lon-lat to UTM (E-N), using Python?

text: I have a lot of KML files set as lon-lat (geographic coordinates) and want to set them as UTM (E-N) and WGS84 with Python (using arcpy, ogr, or gdal) Does anyone knows how to do this?


title: CartoDB multiple sublayers toggle?

text:
CartoDB LayerActions Toggle Is there a way to make one button that shows a sublayer on one click, then hides on the next click. So one button per sublayer which turns the layer on and turns the layer off. Similar to the functionality built into the CartoDB UI. For example how could I reprogram a toggle button by using some other piece of code, or could I expand the function in the HTML file by using a flag value (0 for an ON state which could be show(), and another value such as 1 for an OFF state which would be hide()) This is my code below:                              <link rel="stylesheet" href="http://libs.cartocdn.com/cartodb.js/v3/themes   /css/cartodb.css" />

     Novant Health     All


title: How to sum raster image values with vector overlay

text: I'm new to this and haven't been able to figure out how to do this on my own. Here's what I want to do ... I have two files: a raster image of land use data and a vector-based shapefile of my territories. The raster image has values like so: 11 = water, 21 = developed - open, 22 = developed - low, 23 = developed - medium, 41 = forest, 71 = grassland, 82 = cultivated crop, etc. Overall, I've got about 30 distinct values. The resolution of this image is pretty big (1m squares, I think) and the filesize is about 16 gigs. I was to add a set of columns to my shapefile, one new column for each of the possible land use values. The column names would be like so: value_11, value_21, vaule_22, value_23, etc. Then I'd like to sum up the number of pixels in each shape that correspond to each value. So, if there are 1200 pixels that have a value of 11, then column value_11 would have 1200 in it. I've used qgis to get the vector and raster files to load at the same time and they're georeffed correctly. It displays beautifully. I'm not sure if I'm making this clear, but I'd like a count for every combination of each vector-based polygon and the unique values from the raster image within that polygon. If I need to do this in multiple passes, that's fine. Using qgis, I was able to get average/median/mode/min/max of raster values within each polygon, but that's useless to me. I need the count of each distinct land use value. Eventually, I want to load the final shapefile into SQL Server 2008 and be able to write queries against it. I'd like to take a series of points, find out which polygon they are in, and determine how many water, developed, forest, etc. pixels are in that polygon. Regarding software, I've played with several free/open source programs (such as QGIS, SAGA, etc), but I just don't know enough about what I'm trying to do to get anything to work. I can download whatever I might need, but I don't have a budget, so if ARCGIS or something expensive is the only way to do this, then I'm out of luck. Once I get the shapefile created, I can get it back into SQLServer and I can write the queries there to do what I want. I just need to get build the shapefile with the land use data added to it. Again, I hope I made sense. Any help is greatly appreciated. Thanks.


title: How to go about creating a custom tool/adding special functionality for a web app

text: In my soon to be developed web app, I will have an image service and also a feature service that is essentially an index to the image service, with all kinds of attributes. I would like for the user to be able to select a feature in the index layer, and then get a prompt to view the corresponding image in the image service. I only want the user to see the images that they have 'selected' by selecting features in the index layer. I know python pretty well but am still new with javascript. I was wondering about creating a custom tool in python and then publishing that as a GP service, but I'm unsure of how to manipulate feature/map services in python. I would assume this could be accomplished in javascript but I dont know enough to know how to do that. In terms of the web app, I would try to build it either in arcgis online or leaflet (which I am just learning). Wondering if anyone has some advice on a good approach to this task? I am going to start by investigating how to use the ArcGIS REST API in python.


title: Manually classifying graduated data outside of range in ArcGIS

text: I have two shapefiles representing grids with values that are symbolized with graduated colors, shown below: ![Original](http://i.stack.imgur.com/BiBRo.png) You can see that the two shapefiles encompass similar ranges of data, but assign slightly different colors to ranges in these data. I would like to symbolize the ranges in each of the two shapefiles using the same colors for the same ranges, enabling easier comparison between the files and allowing a single legend to be used. However, when I attempt to classify the ranges with a manual interval, ArcGIS forces what should be the uppermost range to also include the minimum value from the dataset. You can see this in the "-81.64 - 10.00" range of the left data frame, below. This forces the entire layer to be symbolized with this color. **Essentially, ArcGIS is requiring that the minimum value in the dataset be used as a value in one of the ranges.** ![New Intervals](http://i.stack.imgur.com/ngcZt.png) Is there a way to circumvent this feature?


title: Getting Zonal Statistics of Overlapping Circles?

text: I have a Buildings (polygon) shp file and a DSM. what I wanted to do was to get the pixel values for all the vertices of those buildings, and find another pixel value just around those vertices (say 2m, but exclude values coinciding with the building) that will give the maximum difference when subtracted to the value of the vertex. Finally, retain only the vertex that gave the maximum difference for each building. What I did is I used the 'Feature to Point Vertices' in ArcMap to convert my building vertices into a point shp file. Then I made a 2m buffer around these points and I used the Erase tool using the Buildings as the erase feature. I thought I'd get the value for the vertices using 'Extract Values to Point' tool and then use Zonal Stat as Table to get the minimum values on the buffer circles. This is where my problem starts, ArcMap won't perform the Zonal Stat as Table because I have overlapping circles and performing the process on each circle is time consuming. Does anyone have a suggestion to do it in a more efficient manner?


title: What determines direction in "Hatches"?

text: I created a route layer that uses hatches to help determine addresses along a street centerline. About 90% of the time, the hatches go up from zero in the direction that the line was digitized but about 10% of the time (+/-) the hatches go in (what I would say is) the "wrong direction". I can symbolize this layer with arrows showing the direction of the line which shows me the instances where the hatches are wrong. What determines the direction of the hatches and can I do a better job of controlling this? Wow hatches or hatching or centerlines are new tags ? Wish me luck with this one!


title: Main Menu and Catalog Tree View in ArcCatalog are not visible

text: In ArcCatalog (the program, not the window in ArcMap), the Main Menu (`File, Edit, View... Windows, Help`) at the top and the Catalog Tree View have disappeared and I can't find a way to reconnect them. How do I make them visible again?


title: VB Script can't run in Field Calculator

text: I wrote a simple script in VB. But it failed to run in field calculator. I can't figure out what is wrong with it. Below is a screenshot. Thanks for your help. ![enter image description here](http://i.stack.imgur.com/fHNfd.jpg)


title: Convert Shapefile to CSV using FME

text: I have a simple shapefile (point) that has only ID in the attribute table. I can transform shp to csv but this only gets the ID into the CSV file. What transformer/settings do I need to read the x,y coordinates of each point in the shapefile and pass them to the CSV writer?


title: Layer definition Limit Sliverlight API

text: I am using ESRI Sliverlight API 3.1 using layer definition in dynamic service layer, when I create a definition ID in ('1','2','3'......'1000001') . This definition is not getting applied, when I pass big string in REST unable to perform the query. But I need to set a filter based on ID, but I cannot number of ids. How to resolve the issue?                               LayerDefinition def = new LayerDefinition();                     def.LayerID = layerId;                     def.Definition = "ID" + " in " + idString;                     dynamicServiceLayer.LayerDefinitions.Add(def);


title: OpenLayers - Setting the Min Zoom Level of a Google Layer based on a Bounding Box

text: I have an OpenLayers project that is using Google Maps for the base layer. Maybe I'm going about this the wrong way, but I think what I'm trying to do is pretty simple. I just need to set the minimum zoom level of a Google map dynamically based on a given bounding box. This document describes how to set zoom levels, but it seems that for 3rd party layers like Gmaps, you can't use the minExtent/maxExtent features (this would be ideal). Also, is there a way to figure out which Google zoom level corresponds to a given bounding box **before** the map is initialized? I tried:               myGoogleLayer.getZoomForExtent(bbox);      ... but this throws an exception if called before the map is initialized. I can call it after the map is initialized, but then there seems to be no way to set the minZoomLevel on the Gmap without completely borking the zoom consistency between layers and ending up with something like this (where there is clearly no zoom consistency between layers): See a jsfiddle example here ![enter image description here](http://i.stack.imgur.com/GT3W2.png)


title: Adding new vector layer. Should .shp file only be added or .shp plus corresponding .shx, .dbf, etc.?

text: I am a very definite QGIS newbie. I have been creating a map using files that had been used in ArcGis. Each time I wanted to add a new vector layer I opened an .shp file and specified CRS as WGS 84. The resulting layers look fine but when I try to use the "measuring tool" to measure line length I get an error message as follows: _This map is defined with a geographic coordinate system (latitude/longitude) but the map extents suggests that it is actually a projected coordinate system (e.g., Mercator). If so, the results from line or area measurements will be incorrect. To fix this, explicitly set an appropriate map coordinate system using the Settings:Project Properties menu._   * any advice on how I might determine what might be an "appropriate map coordinate system" ?   * I am also wondering if I should have also opened the corresponding .shx, .sbx, .dbf, .sbn and .qix files also ? Thanks


title: Migrating several shapefiles from one server to another?

text: my organization is a large ESRI user and we are under going some hardware changes. My question is: Is there a way to move several shapefiles scattered across many folders (in one server) to another, completely new server, in an efficient manner without compromising the integrity of associated MXD files? I'm aware that I can cut and paste the folders over but that would take ages as I'd have to reset each MXD file's layers data source to the correct path. This is the main concern. Perhaps there might be python solution however my scripting ability is weak. Any ideas?


title: Why does ArcMap crash if I try this Join Field (which ArcPy does not like either)?

text: * I'm a little uncertain whether this question belongs here or on http://gis.stackexchange.com Have an "AddJoin" that results in no shared records, which is not what I'm expecting. Is it because I'm using a text field for the join? Something else? I'm trying to do this on point features from featureclass layers in a file gdb and a shapefile (both expressed as TableViews) on which I've performed a arcpy.SelectLayerByAttribute_management(). No cursors are being used. Thanks for any ideas. ArcGIS 10.0. Also, ArcMap crashes if I try a JoinField on the same fields.               >>> arcpy.GetCount_management("GagesII")     <Result '1'>     >>> for row in arcpy.SearchCursor("GagesII"):     ...     gageid=row.getValue("GAGE_ID")     ...          >>> gageid     u'01170000'     >>> arcpy.GetCount_management("StreamGageEvent")     <Result '1'>     >>> for row in arcpy.SearchCursor("StreamGageEvent"):     ...     sourcefea=row.getValue("SOURCE_FEA")     ...          >>> sourcefea     u'01170000'     >>> if gageid==sourcefea:     ...     print 'yes'     ...          yes     >>>      >>> arcpy.AddJoin_management("GagesII", 'GAGE_ID', "StreamGageEvent", "SOURCE_FEA", "KEEP_COMMON")     <Result 'GagesII'>     >>>      >>> arcpy.GetCount_management("GagesII")     <Result '-1'>


title: Do map services go stale?

text: This has happened several times over: I publish a dynamic map service, and then use the JavaScript API to allow the user to select features to obtain attribute information as well as to zoom to features using zoom buttons located on a left pane, one per feature listed. I test the map, and both functions work properly. I come back the following Monday, and while the attribute query function still works properly, when I click any one of the zoom buttons on the left pane, it fails to zoom to the feature that is listed. Why the difference in behavior? Does the map service need to be periodically refreshed? Seems kind of silly to have to maintain a map service in such a way. The following is the full code:               <link rel="stylesheet" href="http://js.arcgis.com/3.8/js/dojo/dijit/themes/claro/claro.css">     <link rel="stylesheet" href="http://js.arcgis.com/3.8/js/dojo/dojox/grid/resources/Grid.css">     <link rel="stylesheet" href="http://js.arcgis.com/3.8/js/esri/css/esri.css">                               County


title: Any experiences or "gotchas" integrating Whitebox GAT (free and open) with the ArcMap UI?

text: ![Whitebox Geospatial Analysis Tools](http://i.stack.imgur.com/hvq4g.png) I wrote a number of ArcGIS VBA automations in grad school; however, they're fully dependent on the ESRI SpatialAnalyst extension, which is not only closed-source but expensive to the point of deterrence. Since VBA is deprecated and since some researchers at the U still use my VBA tools, I thought it would be fun to rewrite them in .Net. But now, with more experience, I further realize it would be more appropriate for academic use if those utilities consumed open algorithms. With this in mind, I'm considering Whitebox GAT as a potential stand-in for the SpatialAnalyst hydrology tools, and I'm curious if there are any success stories or time-saving "gotchas" related to ArcGIS/Whitebox integration. I anticipate several folks will want to counter-suggest implementing Saga, GRASS, R, et cetera. If this is your position, please describe why pursuing a Whitebox integration would be unwise. For example, does it only support a few input formats, has it poor handling of large (1-2 GB+) files, is it slower than my hockey skating.. You get the idea. :) **[Update]** I did some playing around with the Whitebox UI, and with the help of their tutorials, it wasn't difficult to pre-process a 30-meter DEM I had laying around. Next, after lining-up the hydro rasters, I created a pour point and rendered its watershed. This was enough to get a feel for the Whitebox user experience. Whitebox is extendible and/or consumable using .Net or Python. Having accomplished some basics in the Whitebox UI, I thought I'd chain together the typical DEM pre-processing tasks with a simple .Net automation (no ArcMap just yet). DEM pre-processing usually means the following:   1. set no data value (Whitebox needs this, but Arc never did)    2. fill sinks   3. create a flow direction raster   4. create a flow accumulation raster I put together the following Windows Form "application" (aka `WhiteboxDaisyChain`). It takes a system directory containing an ArcGIS Grid (.FLT) and performs the tasks noted above. If you want to try this, you'll need to download the compiled binaries, unzip, then copy all the `.dll` files from `..\WhiteboxGAT_1_0_7\Plugins` into your project -- I put everything in `..\WhiteboxDaisyChain\Whitebox`. However, this example only needs the four `DLLs` mentioned at the top of the code sample. * * *               using System;     using System.Collections.Generic;     using System.ComponentModel;     using System.Data;     using System.Drawing;     using System.Linq;     using System.Text;     using System.Windows.Forms;          // 1) Create a new Windows Form     // 2) Put all these in a Whitebox folder in the C# project root.     // 3) Add project references to the following and create using statements:          using Interfaces;       // requires Add Reference: Interfaces.dll     using ImportExport;     // requires Add Reference: ImportExport.dll     using ConversionTools;  // requires Add Reference: ConversionTools.dll     using flow;             // requires Add Reference: flow.dll               namespace WhiteboxDaisyChain     {         // 4) Prepare to implement the IHost interface.         // 5) Right-click IHost, select "Implement interface.."         public partial class UI : Form, IHost         {                  // 6) Add a BackgroundWorker object.             private BackgroundWorker worker;                  public UI()             {                 InitializeComponent();                      // 7) Instantiate the worker and set "WorkerReportsProgress".                 worker = new BackgroundWorker();                 worker.WorkerReportsProgress = true;             }                       // 8) Use some event to set things in motion.. i.e. Button click.             private void button1_Click(object sender, EventArgs e)             {                 progressLabel.Text = "Running..";                      // This is the path containing my ArcGrid .FLT.                 // All processing will unfold to this directory.                 string path = "C:\\xData\\TutorialData\\DemWhitebox\\";                      string[] fltArgs = new string[1];                 fltArgs[0] = path + "greene30.flt";                 // in: Arc floating point grid                      // creates a raster in Whitebox data model                 ImportArcGrid importAG = new ImportArcGrid();                 importAG.Initialize(this as IHost);                 importAG.Execute(fltArgs, worker);                  // out:  path + "greene30.dep"                      // set the nodata value on the DEM                 string[] noDataArgs = new string[2];                 noDataArgs[0] = path + "greene30.dep";              // in: my raw DEM                 noDataArgs[1] = "-9999";                            // mine used -9999 as nodata value                      SetNoData setNoData = new SetNoData();                 setNoData.Initialize(this as IHost);                 setNoData.Execute(noDataArgs, worker);              // out:  path + "greene30.dep"                      // fill sinks in the DEM                 string[] fillSinksArgs = new string[4];                 fillSinksArgs[0] = path + "greene30.dep";           // in: my DEM with NoData Fixed                 fillSinksArgs[1] = path + "greene30_fill.dep";      // out: my DEM filled                 fillSinksArgs[2] = "50";                            // the dialog default                 fillSinksArgs[3] = "0.01";                          // the dialog default                      FillDepsBySize fillSinks = new FillDepsBySize();                 fillSinks.Initialize(this as IHost);                 fillSinks.Execute(fillSinksArgs, worker);                      // create a flow direction raster                 string[] flowDirArgs = new string[2];                 flowDirArgs[0] = path + "greene30_fill.dep";        // in: my Filled DEM                 flowDirArgs[1] = path + "greene30_dir.dep";         // out: flow direction raster                      FlowPointerD8 flowDirD8 = new FlowPointerD8();                 flowDirD8.Initialize(this as IHost);                 flowDirD8.Execute(flowDirArgs, worker);                      // create a flow accumulation raster                 string[] flowAccArgs = new string[4];                 flowAccArgs[0] = path + "greene30_dir.dep";         // in: my Flow Direction raster                 flowAccArgs[1] = path + "greene30_acc.dep";         // out: flow accumulation raster                 flowAccArgs[2] = "Specific catchment area (SCA)";   // a Whitebox dialog input                 flowAccArgs[3] = "false";                           // a Whitebox dialog input                      FlowAccumD8 flowAccD8 = new FlowAccumD8();                 flowAccD8.Initialize(this as IHost);                 flowAccD8.Execute(flowAccArgs, worker);                      progressLabel.Text = "";                 progressLabel.Text = "OLLEY-OLLEY-OXEN-FREE!";             }                            /* IHost Implementation Methods Below Here */                  public string ApplicationDirectory             {                 get { throw new NotImplementedException(); }             }                  public void ProgressBarLabel(string label)             {                 this.progressLabel.Text = "";                 this.progressLabel.Text = label;                    // This is the only one I used.             }                  public string RecentDirectory             {                 get                 {                     throw new NotImplementedException();                 }                 set                 {                     throw new NotImplementedException();                 }             }                  public bool RunInSynchronousMode             {                 get                 {                     throw new NotImplementedException();                 }                 set                 {                     throw new NotImplementedException();                 }             }                  public void RunPlugin(string PluginClassName)             {                 throw new NotImplementedException();             }                  public void SetParameters(string[] ParameterArray)             {                 throw new NotImplementedException();             }                  public void ShowFeedback(string strFeedback, string Caption = "GAT Message")             {                 throw new NotImplementedException();             }         }     }      So far I'm digging this, but I don't yet have a real success story or any show-stoppers to describe.. My next goal will be interactively submitting pour points from ArcMap. Basically, I want to click the map ..get the watershed.


title: Disable leaflet interaction temporary

text: How can I temporary disable zoming/draging the Mapview in Leaflet.js Tried so many ways but without any luck. It's important to make it temporary and I also need the option to enable again.


title: QGIS Processing error with SAGA, "output layer could not be open"

text: QGIS 2.4.0   SAGA GIS 2.0.8 > DESCRIPTION ERROR when processing with SAGA. I used the "Shapes to Grid" function. The error given: "Oooops! The following output layers could not be open". This error seems to be very common, but I can't find any solution that works. I have read various posts, but there is no answer to fix this error that seems to work. I am currently using the latest version of both QGIS and SAGA GIS. > ERROR: Oooops! The following output layers could not be open Grid: C:\Users\Burger\AppData\Local\Temp\processing\d6e5c059ffc247328315ea134947d3da\USERGRID.tif   The above files could not be opened, which probably indicates that they were not correctly produced by the executed algorithm Checking the log information might help you see why those layers were not created as expected   This algorithm requires SAGA to be run. A test to check if SAGA is correctly installed and configured in your system has been performed, with the following result:   SAGA seems to be correctly installed and configured > LOG: Algorithm Shapes to grid starting... grid_gridding "Shapes to Grid" -INPUT "C:/Pieter & Wiesie Burger/BackUP Werk/Overberg Agri/Geografies/Data/Roodebloem/Stroper GPS/2012 Stroop/RB7 - BAR 2012 - 12110600 - filter points.shp" -FIELD "ERTRAG" -MULTIPLE 0 -LINE_TYPE 0 -GRID_TYPE 3 -USER_XMIN 24.53824 -USER_XMAX 14.54312 -USER_YMIN -29.24768 -USER_YMAX -39.24041 -USER_SIZE 10 -USER_GRID "C:\Users\Burger\AppData\Local\Temp\processing\d6e5c059ffc247328315ea134947d3da\USERGRID.tif.sgrd" io_gdal 1 -GRIDS "C:\Users\Burger\AppData\Local\Temp\processing\d6e5c059ffc247328315ea134947d3da\USERGRID.tif.sgrd" -FORMAT 1 -TYPE 0 -FILE "C:\Users\Burger\AppData\Local\Temp\processing\d6e5c059ffc247328315ea134947d3da\USERGRID.tif" C:\PROGRA~1\QGISCH~1\bin>set SAGA=C:/PROGRA~1/QGISCH~1/apps\saga C:\PROGRA~1\QGISCH~1\bin>set SAGA_MLB=C:/PROGRA~1/QGISCH~1/apps\saga\modules C:\PROGRA~1\QGISCH~1\bin>PATH=PATH;C:/PROGRA~1/QGISCH~1/apps\saga;C:/PROGRA~1/QGISCH~1/apps\saga\modules C:\PROGRA~1\QGISCH~1\bin>saga_cmd grid_gridding "Shapes to Grid" -INPUT "C:/Pieter & Wiesie Burger/BackUP Werk/Overberg Agri/Geografies/Data/Roodebloem/Stroper GPS/2012 Stroop/RB7 - BAR 2012 - 12110600 - filter points.shp" -FIELD "ERTRAG" -MULTIPLE 0 -LINE_TYPE 0 -GRID_TYPE 3 -USER_XMIN 24.53824 -USER_XMAX 14.54312 -USER_YMIN -29.24768 -USER_YMAX -39.24041 -USER_SIZE 10 -USER_GRID "C:\Users\Burger\AppData\Local\Temp\processing\d6e5c059ffc247328315ea134947d3da\USERGRID.tif.sgrd" * * * ## ##### ### ### ## ### # ## ## #### # ### ##### ## # # ## ##### # * * * library path: C:/PROGRA~1/QGISCH~1/apps\saga\modules library name: grid_gridding module name : Shapes to Grid author : O.Conrad (c) 2003 * * * go... Load shapes: C:/Pieter & Wiesie Burger/BackUP Werk/Overberg Agri/Geografies/Data/Roodebloem/Stroper GPS/2012 Stroop/RB7 - BAR 2012 - 12110600 - filter points.shp... ready okay Parameters Shapes: RB7 - BAR 2012 - 12110600 - filter points Attribute: ERTRAG Method for Multiple Values: first Method for Lines: thin Preferred Target Grid Type: Floating Point (4 byte) Target Grid: user defined **error:** executing module [Shapes to Grid] C:\PROGRA~1\QGISCH~1\bin>saga_cmd io_gdal 1 -GRIDS "C:\Users\Burger\AppData\Local\Temp\processing\d6e5c059ffc247328315ea134947d3da\USERGRID.tif.sgrd" -FORMAT 1 -TYPE 0 -FILE "C:\Users\Burger\AppData\Local\Temp\processing\d6e5c059ffc247328315ea134947d3da\USERGRID.tif" * * * ## ##### ### ### ## ### # ## ## #### # ### ##### ## # # ## ##### # * * * library path: C:/PROGRA~1/QGISCH~1/apps\saga\modules library name: io_gdal module name : GDAL: Export Raster author : O.Conrad (c) 2007 * * * go... Load grid: C:\Users\Burger\AppData\Local\Temp\processing\d6e5c059ffc247328315ea134947d3da\USERGRID.tif.sgrd... **failed** **error** : Grid file could not be opened. **error** : input file [C:\Users\Burger\AppData\Local\Temp\processing\d6e5c059ffc247328315ea134947d3da\USERGRID.tif.sgrd] **error** : empty input list [GRIDS] Usage: 1 -GRIDS [-FILE ] [-FORMAT ] [-TYPE ] -GRIDS: Grid(s) Grid list (input) -FILE: File File path -FORMAT: Format Choice Available Choices: [0] Virtual Raster [1] GeoTIFF [2] National Imagery Transmission Format [3] Erdas Imagine Images (.img) [4] ELAS [5] In Memory Raster [6] MS Windows Device Independent Bitmap [7] PCIDSK Database File [8] ILWIS Raster Map [9] SGI Image File Format 1.0 [10] Leveller heightfield [11] Terragen heightfield [12] ERMapper .ers Labelled [13] Raster Matrix Format [14] Idrisi Raster A.1 [15] Intergraph Raster [16] Golden Software Binary Grid (.grd) [17] Portable Pixmap Format (netpbm) [18] ENVI .hdr Labelled [19] ESRI .hdr Labelled [20] PCI .aux Labelled [21] Vexcel MFF Raster [22] Vexcel MFF2 (HKV) Raster [23] VTP .bt (Binary Terrain) 1.3 Format [24] Image Data and Analysis [25] NOAA Vertical Datum .GTX [26] NTv2 Datum Grid Shift [27] ARC Digitized Raster Graphics [28] SAGA GIS Binary Grid (.sdat) -TYPE: Data Type Choice Available Choices: [0] match input data [1] 8 bit unsigned integer [2] 16 bit unsigned integer [3] 16 bit signed integer [4] 32 bit unsigned integer [5] 32 bit signed integer [6] 32 bit floating point [7] 64 bit floating point **error** : executing module [GDAL: Export Raster] C:\PROGRA~1\QGISCH~1\bin>exit Converting outputs Loading resulting layers Algorithm Shapes to grid finished


title: Using Spatial References with ArcPy after *.prj files no longer installed by ArcGIS 10.1 for Desktop?

text: I'm writing a python script that uses the CreateFeatureclass_management tool, without a template. One of the required parameters is a spatial reference. I used to be able to point to the system *.prj file in 9.3, but that's no longer an option, and I don't want to point to an existing shapefile's *.prj because that will cause issues when I share the script. Now that Esri doesn't install all the *.prj files with ArcGIS, how can I set the projection in a script without referencing an existing shapefile?


title: How to select WFS features by polygon in OpenLayers?

text: I need to be able to select features on a wfs layer using a polygon instead of a box. I was hoping this would be a standard option in the SelectFeature class. Too bad this isn't the case. Any implementation suggestions?


title: tif tfw to GeoTiff

text: How do I combine a .tif file with a .tfw file to create a GeoTiff? There are many answers telling me to use gdal, but I have no clue. So can someone provide a step by step example of how to do this please?


title: GeoTIFF created using ArcGIS for Desktop not displaying in Google Earth or Photoshop?

text: I originally had an IMAGINE Land Cover raster file that I opened in ArcMap. I used the Extract Feature from Raster to extract only the Developed class and saved it (still IMAGINE). Then I used the Copy Raster to convert the IMAGINE to GeoTIFF, but when I try to view the GeoTiff in Google Earth or Photoshop, it doesn't display anything. However, in ENVI, it displayed successfully. What am I doing wrong?


title: Raster attribute field error

text: I am trying to add jpegs to a raster data type attribute field. I want to this field to be managed by the file geodatabase. When I load the photos and try to save my edits it get various errors but usually I get the error posted below. If i set the field to non-managed it works just fine. I also have an attachment table that I can load photos too and that works just fine as well. This is happening on evey geodatabase that I try. Anybody have any ideas at all? I have post a similar question before but still haven't![enter image description here](http://i.stack.imgur.com/JKnpy.png) found an answer.


title: How to get vector data from OSM to use in QGIS?

text: How do I get data from OpenStreetMap which I can then edit in QGIS layer by layer? What is the best OSM output and what / how do I get to a suitable input for QGIS?   OR...is there another appraach to this? Thanks


title: Where can I find administrative boundary data for Mexico?

text: The USA has states, counties, and cities. How are Mexico's territories split up, hierarchically? It also has states and cities, but does it have an intermediate-level political unit? Where might I find sets of boundary data for such units?


title: How to open QGIS attribute table in form view using Python?

text: I've got a plugin where I load some layers. On the last layer loaded I show the attribute table. I would like to select the view different from the table view (see image). How can I do this in Python? ![Attribute table layout options](http://i.stack.imgur.com/EZZNC.png)


title: How to extract data from Geodatabase?

text: We have been happily using the 2010 Census Demographic Profile data in Shapefile format for some time. We now want to start using the ACS 5-year Estimates Data Profiles which, sadly, are only available in Geodatabase formats. What's the best solution? Do I bite the bullet and buy some version of ArcGIS, are there any services that can convert the files for me or is there some conversion program I can purchase? Note that I am not a GIS specialist. I've learned enough to do what I want with the 2010 Census demographic shapefiles, but this is not the main focus of my website. Recommendations and advice would be greatly appreciated.