title: Conditional reclassification of a raster

text: I am attempting to do a conditional reclassification of an NLCD land-use raster. I specifically want to reclassify all developed areas (classifications 22,23,24) that are near a road (classification 21) into distinct new classes: "developed land near road". My workflow so far using ArcMap has been:   1. Extract the road class into a new raster using _"Extract by Attributes"_   2. Convert that raster to point using _"Raster to Point"_   3. Buffer around those points using _"Buffer"_   4. Reconvert the resulting polygon back to a raster using _"Polygon to Raster."_ This file has a single value for being within distance of a road and NA's otherwise/   5. ...   6. profit! So step 5 is where I have trouble. I need to do some sort of a conditional reclassification on each of the 3 housing classes. In raster calculator I think I worked out that it needs to be of the form: `Con(IsNull("buff1_21_Buf60Rast"),"lu_1.asc",Diff("buff2_21_Buf30Rast","lu_1.asc"))` I'm suspecting that I will need to run it sequentially for each of the 3 housing values. But I'm not really sure how or which of the logical math operators to use. Thoughts?


title: How to calculate volume from a raster file?

text: I have a raster file containing the depth of water. I need to find the volume of water using that raster file. I need to find the area of each cell and multiply with the depth. How can I do it in Python using ArcGIS 10 ? * * *               rasvol = 0     rasterArea=gp.searchCursor(reservoir_area)     rasRow=rasterArea.next()     while rasRow:          rasCount=rasRow.getValue('Count')          rasValue=rasRow.getValue('Value')          rasvol=rasvol + rasCount*rasValue*cellarea          rasRow=rasterArea.next()      This code is correct or not for calculating the volume ?


title: How to have two OpenLayers maps on the same page

text: Specifically, my question is: how to assign the maps to different divs? In my page it kept placing the two maps under each other, although I defined the divs to be horizontally aligned.   So I have been searching for the answer to this question, when I discovered it by mistake, and I think it is not obvious so I thought I might share it.


title: Calculate area of intersecting polygons in QGIS

text: I have two collections of shapefiles: 1 - postal codes, administrative and political areas 2 - various land use categories like parks & allotments, and points like trees I want to produce statistics on the area of those land use categories broken down by postal code, admin & political area, and the number of points by the same means. For example, there are 4.1ha of park in SE15, 879 trees in SE22, and so on. I am using QGIS. I've processed my raw data to get all the shapefiles ready, but after much fruitless searching around on this and other web sites I can't see a way to do this.


title: GeoDjango: transform() not working?

text: I have a point: `pnt = GEOSGeometry('POINT({0} {1})'.format(lon, lat))`, and tried `pnt.transform(900913)`, and it gives error:               GEOSException: Calling transform() with no SRID set is not supported      Then I tried transform() with my points in the db:               qs = T.objects.all().transform(900913)      and it works fine:               >>print qs[0].geom     >>POINT (-10887469.2136471699923277 3564159.8204717375338078)      Here's my model T:               class T(models.Model):         name = models.CharField(max_length=200)         geom = models.PointField(srid=4326, blank=True, null=True, default=None)         objects = models.GeoManager()      Any Idea? Thanks very much!


title: How to convert from .SHP to .MAP for mapsforge?

text: We are using **mapsforge** , a Java component for Android. This component use a proprietary format file (.map) that can be download from OpenStreetMap (OSM). I need to use .shp maps, so I want to know how can I convert a .shp file into a .map file for mapsforge for Android


title: Lock during ApplyEdits

text: Have a JavaScript application trying to push data to a FeatureServer on ArcGIS Server 10.1. I get this error occasionally: {Error code '400': 'Unable to complete operation.' Lock request conflicts with an established lock [dbo.DEFAULT]}


title: GeoDjango: Only numeric values of degree units are allowed on geographic DWithin queries

text: I'm trying to get N points within 5000 meters of a certain point. My database is Postgres(PostGis). I have a model like this:               class Theatre(models.Model):         geom = models.PointField(srid=4326, blank=True, null=True, default=None)         objects = models.GeoManager()      Then I ran a couple lines in django shell: > > > Theatre.objects.filter(geom__dwithin=(GEOSGeometry('POINT(30.111199 > -97.309990)'), D(m=5000))) and I get error:               ValueError: Only numeric values of degree units are allowed on geographic DWithin queries.      I found the document in here and it says calculating the geometry distance.. Please help.. Thanks very much!!


title: GPS-Plugin and Garmin data

text: I am an avalanche specialist and I use a GPS Garmin 62 for fieldwork. I for example use waypoints to measure snow depth in avalanches for volume calculations. When I add a waypoint in the garmin unit it allows me to enter digits in the predefined field called depth. This field depth can be seen in the attribute table for GPX in Garmin Basecamp. However the QGIS does not import this field, for it only imports standard GPS data. Does anybody have a solution on how to import this field into Quantum GIS?


title: adding polygon to database in openlayers

text: I want to add a polygon to the database. I used the following code. The problem is that I am to able to save in database. Whatever polygon I draw on map it gets save temporary and as soon as I reload the application all the polygon vanishes. Please help               wfs_layer = new OpenLayers.Layer.Vector("wfs_layer", {            strategies: [new OpenLayers.Strategy.Fixed(),saveStrategy],            visibility: false,            protocol: new OpenLayers.Protocol.WFS({            version: "1.0.0",            srsName: "EPSG:4326", // this is the default            url:  "http://localhost:8089/geoserver/IIRS1/wms",            featureType: "india_state",            featureNS: "http://www.opengeospatial.net/IIRS1",            extractAttributes: true,        geometryName: "the_geom",     schema: "http://localhost/geoserver/wfs/DescribeFeatureType?version=1.1.0&;typename=IIRS1:india_state"         })                                  }     );          vector = new OpenLayers.Layer.Vector('vector');     map.addLayers([wfs_layer,vector]);      // saving               function saveSuccess(event) {               alert('Changes saved');               }                 function saveFail(event) {                alert('Error! Changes not saved');                   }                  var saveStrategy = new OpenLayers.Strategy.Save();             saveStrategy.events.register('success', null, saveSuccess);             saveStrategy.events.register('fail', null, saveFail);      //adding draw polygon button to tool bar                   action = new GeoExt.Action({         text: "draw_poly",         control: new OpenLayers.Control.DrawFeature(             vector, OpenLayers.Handler.Polygon),         strategies: [saveStrategy],         map: map,         toggleGroup: "draw",         allowDepress: false,         tooltip: "draw_polygon",          });     actions["draw_poly"] = action;               toolbarItems.push(action);      toolbarItems.push("-");      // Saving button               action = new GeoExt.Action({         text: "Save",         control: new OpenLayers.Control.Button({                trigger: function() {           if(feature) {                       feature.selectControl.unselectAll();                       }           saveStrategy.save();       }         }),         map: map,         // button options         toggleGroup: "Save",         allowDepress: false,         tooltip: "Save",          });     actions["Save"] = action;     toolbarItems.push(action);      toolbarItems.push("-")


title: Problem running the pgr_trsp algorithm from pgrouting

text: I was studying the pgRouting Documentation 2.0. I was trying to run the pgr_trsp algorithm as shown in the tutorial: SELECT seq, id1 AS node, id2 AS edge, cost FROM pgr_trsp('SELECT id, source, target, cost FROM edge_table',7, 12, false, false); But i have this result: seq|node|edge|cost 0|-1|0|0 I check the edge_table table a it´s works for other algorithms such as pgr_astar, but for pgr_trsp don´t work. I also tried with other edges table, but I always get the same result. The sample data is the provided by the documentation 2.0: DROP TABLE edge_table CREATE TABLE edge_table ( id serial, dir character varying, source integer, target integer, cost double precision, reverse_cost double precision, x1 double precision, y1 double precision, x2 double precision, y2 double precision, the_geom geometry ); INSERT INTO edge_table (cost,reverse_cost,x1,y1,x2,y2) VALUES ( 1, 1, 2,0, 2,1); INSERT INTO edge_table (cost,reverse_cost,x1,y1,x2,y2) VALUES (-1, 1, 2,1, 3,1); INSERT INTO edge_table (cost,reverse_cost,x1,y1,x2,y2) VALUES (-1, 1, 3,1, 4,1); INSERT INTO edge_table (cost,reverse_cost,x1,y1,x2,y2) VALUES ( 1, 1, 2,1, 2,2); INSERT INTO edge_table (cost,reverse_cost,x1,y1,x2,y2) VALUES ( 1,-1, 3,1, 3,2); INSERT INTO edge_table (cost,reverse_cost,x1,y1,x2,y2) VALUES ( 1, 1, 0,2, 1,2); INSERT INTO edge_table (cost,reverse_cost,x1,y1,x2,y2) VALUES ( 1, 1, 1,2, 2,2); INSERT INTO edge_table (cost,reverse_cost,x1,y1,x2,y2) VALUES ( 1, 1, 2,2, 3,2); INSERT INTO edge_table (cost,reverse_cost,x1,y1,x2,y2) VALUES ( 1, 1, 3,2, 4,2); INSERT INTO edge_table (cost,reverse_cost,x1,y1,x2,y2) VALUES ( 1, 1, 2,2, 2,3); INSERT INTO edge_table (cost,reverse_cost,x1,y1,x2,y2) VALUES ( 1,-1, 3,2, 3,3); INSERT INTO edge_table (cost,reverse_cost,x1,y1,x2,y2) VALUES ( 1,-1, 2,3, 3,3); INSERT INTO edge_table (cost,reverse_cost,x1,y1,x2,y2) VALUES ( 1,-1, 3,3, 4,3); INSERT INTO edge_table (cost,reverse_cost,x1,y1,x2,y2) VALUES ( 1, 1, 2,3, 2,4); INSERT INTO edge_table (cost,reverse_cost,x1,y1,x2,y2) VALUES ( 1, 1, 4,2, 4,3); INSERT INTO edge_table (cost,reverse_cost,x1,y1,x2,y2) VALUES ( 1, 1, 4,1, 4,2); INSERT INTO edge_table (cost,reverse_cost,x1,y1,x2,y2) VALUES ( 1, 1, 0.5,3.5, 1.999999999999,3.5); INSERT INTO edge_table (cost,reverse_cost,x1,y1,x2,y2) VALUES ( 1, 1, 3.5,2.3, 3.5,4); UPDATE edge_table SET the_geom = st_makeline(st_point(x1,y1),st_point(x2,y2)), dir = CASE WHEN (cost>0 and reverse_cost>0) THEN 'B' -- both ways WHEN (cost>0 and reverse_cost<0) THEN 'FT' -- direction of the LINESSTRING WHEN (cost<0 and reverse_cost>0) THEN 'TF' -- reverse direction ELSE '' END; SELECT pgr_createTopology('edge_table',0.001); Other additional info: select * from pgr_version(); "PostgreSQL 9.3.2, compiled by Visual C++ build 1600, 64-bit" select postgis_full_version(); "POSTGIS="2.1.1 r12113" GEOS="3.4.2-CAPI-1.8.2 r3924" PROJ="Rel. 4.8.0, 6 March 2012" GDAL="GDAL 1.10.0, released 2013/04/24" LIBXML="2.7.8" LIBJSON="UNKNOWN" TOPOLOGY RASTER" select version();   "PostgreSQL 9.3.2, compiled by Visual C++ build 1600, 64-bit" Can anybody help with this issue?


title: Exporting lists of Feature classes and fields from SDE

text: I'm using matt wilkies' arcplus module http://gis.stackexchange.com/a/5943/16793 to export a list of feature classes within an SDE database to a text file. I've also written a tool to write field names to a csv file:                  """The source code of the tool."""       import arcpy, os, csv               fclass = "SDEConnection\\"+Feature     fieldnames = [f.name for f in arcpy.ListFields(fclass)]     writepath = r"C:\Users\ME\Desktop\TableFields\\"+Feature+".csv"     with open(writepath, "w") as output:             writer = csv.writer(output, lineterminator='\n')             for val in fieldnames:         writer.writerow([val])      Ideally, I'd like to loop through the feature classes and fields and write a write file containing something like below: Feature Class or Table-----------------------------------Fields WELLS.WELL_DATA-------------------------------------WELL_KEY, SITEID, DEPTH, DIAMETER I'm a little stumped on nesting the loops and formatting the output like this. \-----------EDIT----------------------------------------------------------- Using the Walk method with a script like this               import os, arcpy, csv     def inventory_data(workspace, datatypes):     for path, path_names, data_names in arcpy.da.Walk(             workspace, datatype=datatypes):         for data_name in data_names:             yield os.path.join(path, data_name)      obtains a unicode string like **"C:\Users\user\AppData\Roaming\Esri\Desktop10.2\ArcCatalog\Connection to oracle.sde\BASINS\BASINS.ONE"** while               fn = [f.name for f in arcpy.ListFields(fc)]     print fn       extracts a list **[u'FNODE_', u'TNODE_', u'LPOLY_', u'RPOLY_',]** My aim is to combine Walk with ListFields in a loop and write the FeatureClass with its accompanying fields (delimited) to a file.


title: Create an X,Y latitude and longitude grid from an raster

text: I have a satellite image (LandSat 5), I want create a raster with X,Y coordinates, this is possible in Model Builder of ArcGis 10.0? or other option. I know that the option $$YMAP ($$XMAP) not work more in ArcGis 10.0 I want make the same. But my input is a raster (landsat image, UTM projection), I need obtain (output) a grid of latitude and a grid of longitude (two grids: latitude and longitude). First, this code is valid, and if, is possible work with ModelBuilder. I'm novel in arcgis and python, I working with evapotranspiration and I need this grids for calculate it. This is the code of video (www.youtube.com/watch?v=_-6apNB18gk): ![enter image description here](http://i.stack.imgur.com/mU55D.jpg)


title: No data appears when using gdal_calc.py but source had data values

text: I'm in the process of calculating topographic exposure using a mix of GDAL, grass, R and python in Windows. I need to process a very large tiff file using the following calculation: arctan(topex_e.tif)*57.295779513082323 I attempted to perform this using gdal_calc.py by running the following in from GRASS: `python C:/OSGeo4W/bin/gdal_calc.py --creation-option="COMPRESS=DEFLATE" --creation-option="BIGTIFF=YES" -A C:/GIS/topex_e.tif --outfile=C:/GIS/topex_e_deg.tif --calc="arctan(A)*57.295779513082323"` This runs fine with no errors and mostly produces the desired output. However, there are nodata values scattered all over the map, even though the source tif doesn't have nodata values. My mathematical knowledge tells me that arctan should always return a result between -Pi/2 and Pi/2. See: http://docs.scipy.org/doc/numpy/reference/ufuncs.html#trigonometric-functions Is there any reason why the result would contain Nodata values? Edit: If this helps, the source .tif file is a compressed GeoTIFF with Float32 values. with approximately 95000 by 85000 pixels. I've imported into grass and am in the process of using r.mapcalc so I will hopefully work around the issue, but am concerned about the reliability of gdal_calc.py


title: Create fishnet with row Iterator: fishnet extent does not update to reflect each row feature

text: I'm a writing about a geoprocess that I can do on an individual basis but for some confounding reason am unable to do with the row iterator in Model Builder. Given that I want to repeat the same process 68 times, I hope that I can be spared the pain of doing this the long way (one-by-one). I have tried to be as complete as possible in describing the problem, including pictures and detailed commentary, hence the length of this post. The task in question is exceedingly simple. Let's say there is one square. I want to create a fishnet bounded by the extent of the square. In the create fishnet tool, in the " _template extent_ " drop down box, I select the file with the square. Each square happens to be 2200 meters by 2200 meters. So I create 22 rows and 22 columns of 100 by 100 meter cells. The entire square is covered by the fishnet. So far so good. I can do this part. Now let's say I have a shapefile that contains 68 rows. Each of these rows corresponds to one of 68 squares scattered across a metropolitan area (all the square are the same size). This is what I am talking about: one shapefile in the table of contents and 68 squares: ![enter image description here](http://i.stack.imgur.com/8kK4B.jpg) You will see that the attribute table is quite normal: ![enter image description here](http://i.stack.imgur.com/bEhv1.jpg) I want to do exactly what I have described above with one square, except now the file in question contains 68 rows instead of one and I want to use the row iterator function to execute the task on each row of the file. I set up the process like so in Model Builder: the file with the 68 rows is the "input table" to the "iterate row selection" tool. If I understand correctly, each row is treated as a unique item by the darker green "I_FeatureEnvelope_68HHs." I subsequently want to treat each of these unique features as inputs to the "create fishnet" tool. ![enter image description here](http://i.stack.imgur.com/4jZOy.jpg) When I go to connect the "I_FeatureEnvelope_68HHs" to the "Create Fishnet" box I see a pop up box that says _Template Extent_. I click _Template Extent_ , thinking that it means: "treat each feature represented by each of the 68 rows as unique extents to be used in creating 68 different fishnets." (sort of like if you were to open the attribute table, highlight a row and click 'zoom to selected' each row corresponds to an extent) When I double click on the Create Fishnet box, I see this: ![enter image description here](http://i.stack.imgur.com/zdD7b.jpg) _Template Extent_ has already been chosen for me. The selection appears to make sense. I set the cell width and height and numbers of rows and columns: 100m x 100m and 22 x 22 as I stated before. I select polygon geometry type at the bottom. I run this: ![enter image description here](http://i.stack.imgur.com/cTpw4.jpg) Let me tell you what _I think_ I should see. To recap, there are 68 rows. Each of these corresponds to a square. I am creating 22 row x 22 column fishnets for each square which leads to a 484 cells fishnet per square. That means the output file should have 68 x 484, or 32,912 rows (or 484 cells per each square). Let me tell you what actually happened. My process led to the creation of only one fishnet in the lower left hand corner: ![enter image description here](http://i.stack.imgur.com/XcF8Y.jpg) Close up it looks alright, a 22 x 22 grid with 484 cells: ![enter image description here](http://i.stack.imgur.com/k8dAc.jpg) I have a faint idea of what is going on. Consider the file with the 68 squares. The one fishnet that was created corresponds to the lower left hand corner of the full extent of this file - so the Template Extent is not considering each of the 68 rows as unique items. It simply takes the full extent of the 68 row file and focuses on the lower left hand corner. I have played around changing the settings a bit, but to document all that I'll end up writing a novel. To make a long story short, I have found it impossible to use the row iterate tool to create a fishnet that corresponds to the extent of each of the 68 features in the file in question. Any helps, tips, clues, etc... would be immensely appreciated!


title: Qgis - Is there a way to split a polygon layer with a line layer?

text: I am drawing administrative regions that have boundaries follow streets in another layer. As these streets are pretty long with many vertices, I don't want to re-trace them. I am looking for a way to "split" the polygon like "split features" tool (the one with the scissors icon). Before: one whole polygon with a line running through it ![enter image description here](http://i.stack.imgur.com/5kadB.png) After: the line works as a pair of scissors and split the polygon into 2 parts (in the picture below, the polygons were moved apart for illustrating purpose only) ![enter image description here](http://i.stack.imgur.com/s5hYX.png) Thanks a lot.


title: How to create Feature Class based on DBF schema?

text: I have an attribute table schema (about 30 fields) stored in an Excel document, which I exported as a dbf. I can easily make a geodatabase table from this, using the Excel sheet as a template. How can I create a feature class from this schema (either from the Excel sheet or the dbf or the geodatabase table)? The Excel sheet is blank except for the field names.


title: LoadingPanel - Waiting for Vector & WFS loading

text: I have vector layer with WFS protocol. and because I have a huge amount of data, there is a while before features loaded on the map. So i need to tell user to wait a moment until they load by display a "waiting" animation. I found this Waiting for Vector & WFS loading and tried it with no luck so far. here is what I did:   * I placed the LoadingPanel.js   * I added this stylesheet:              <style type="text/css">          .olControlLoadingPanel {                  background-image:url(loading.gif);                  position: relative;                  width: 195px;                  height: 11px;                  background-position:center;                  background-repeat:no-repeat;                  display: none;              }           </style>`        * and finally added the LoadingPanel control : `map.addControl(new OpenLayers.Control.LoadingPanel());`


title: How to iterate over selected layers?

text: In QGis by pressing shift or ctrl keys a user can select multiple layers in toc. How can a python script iterate over selected layers? I'd expect a iface.selectedLayers() but I can find only iface.activeLayer() which return only one layer.


title: Formatting Google Places API GeoJSON and Import to QGIS?

text: Hi Stackexchange I'm trying to export JSON or XML data to my desktop and import it to QGIS. This is the URL I'm using but it appears to be incorrectly formatted: https://maps.googleapis.com/maps/api/place/nearbysearch/json?location=51.503186,-0.126446&radius=9000&keywords=exchange&types=bank&key=AIzaSyDyGA1g71hnK85zIMZ7CVp0CKoiRfG_Fpk My goal is to query the Google API and receive a JSON or XML file that I can then load into QGIS/ArcGIS ... etc as a feature dataset. The JSON text that is received with the API query has incorrect syntax and this fails to load, please advise?


title: QGIS is locking the clipboard from other application? (SitesPlus)

text: We are starting to use QGIS along with Sales Forecasting software Sites Plus by Gravitec. An ugly problem just popped up: Whenever QGIS is open, any attempts to calculate a data table in SitesPlus results in an error stating "Unable to open Clipboard" and the changes are not committed. This happens even when just opening a blank QGIS project, without adding any data. If I close out QGIS, the problem goes away. Is there anything that could be causing QGIS to lock up the clipboard? Anyone else having a similar issue? I've tried clearing the clipboard before attempting edits to no avail. QGIS seems to lock it up once it's open. Thanks! Running QGIS 2.2 Valmiera & SitesPlus 2010 on a Win 7 machine Have duplicated the error on multiple machines.


title: Add band to GDAL VRT dataset and convert to TIFF

text: I'm trying to add a band to a TIFF dataset by converting it to a VRT format and using the `AddBand()` function, then convert it back to a TIFF to access the band and write the data. Here is an example of my code:               const char *formatVRT = "VRT";     const char *formatTIFF = "GTiff";     char *filenameTIFF;     char *filenameVRT;     GDALDriver pDriverTIFF, pDriverVRT;     GDALDataset pTiffDS, pVrtDS;          GDALALLRegister();          pDriverTIFF = GetGDALDriverManager()->GetDriverByName(formatTIFF);     pDriverVRT = GetGDALDriverManager()->GetDriverByName(formatVRT);          pTiffDS = pDriverTIFF->Create(filenameTIFF, cols, rows, bands, GDT_Float32, NULL);     pVrtDS = pDriverVRt->CreateCopy(filenameVRT, pTiffDS, FALSE, NULL, NULL, NULL);          pVrtDS->AddBand(GDT_Float32, NULL);     pTiffDS = pDriverTIFF->CreateCopy(filenameTIFF, pVrtDS, FALSE, NULL, NULL, NULL);      I get the following errors: > ERROR 1: Deleting (filenameTIFF) failed: > > Permission denied > > ERROR 1: TIFFFillStrip:Read error at scanline 25; got 0 bytes, expected > 24840 Any suggestions are appreciated.


title: How to display a "Save As" dialog for raster files?

text: I'm writing a QGIS plugin in python. I need to show a "Save as" dialog window in order to allow the user to select the path of the output raster file. Do I need to create a QFileDialog and manually add all the raster formats supported by QGIS or does it exist a method to load these formats from somewhere? Or, better, is there a method to show a "Save As" dialog window ready to manage raster files? Thanks!


title: converting oracle geometry to postgis with 3d data (wkt?)

text: I am converting Oracle spatial database to postgis. Most of of the spatial data is 2.5d. Meaning we have a Z-coordinate which indicates the depth/height of the element compared to soil level. Getting the 2D elements is pretty easy: I select the WKT presentation of the element and insert that back into postgis. But, in Oracle the `sdo_util.to_wktgeometry` only works on 2D geometries. I can convert geometries to 2D using `sdo_cs.make_2d` but only for those elements that have no height at all. What is the best way to convert the data?   * is WKB exchangeable between databases? (I guess not)   * are there other formats to extract data from oracle in a format postgis understands? Thanks.


title: Count Points in Polygon with Postgis

text: I've got a simple problem: I want to count the number of points within a set of polygons. I have a SQL already but it only gives back the gid of the polygone that actually contains points. My tables: a polygon layer with 19.000 rows and a point layer with 450 rows. The following SQL               select grid.gid, count(*) AS totale FROM grid, kioskdhd3 WHERE      st_contains(grid.geom,kioskdhd3.geom) GROUP BY grid.gid;      return only some 320 polygons that actually contain points. But I want all polygons returned, even thought the number of points is 0. Of course it has to do with my WHERE-clause. Where do I have to put in my st_contains? Thank you Carsten


title: QGIS 1.8.0 Mac where is measure window or measure area symbol?

text: Just trying to measure the area of a selected polygon (a few km2) but can't find out where the measure symbol is (mentioned in some posts here) or how to get the measure window shown in the manual. How?


title: ArcGIS - Only draw data in current layout view

text: I am using county data to create a map of a very localized area (about a city block). The data includes parcel, sewer, water, electric, stormwater, imagery, centerlines, etc. Is there any way to force ArcMap to only draw the data that is within my current layout view? Each time I move the map an inch or add something it seems like I'm waiting a long time for it to draw everything out in the background. I know I can select features manually and create a new layer or shapefile with them individually, however, there are quite a few dozen feature classes that need to be shown on this map. This is something that would probably be useful in general and I'm sure it's a setting I have overlooked.


title: How to refer to Obj file when storing tab file name in a String variable?

text: I need to loop through tab files and use select, am doing the following:                Dim file_name as string      Dim file_name_obj as string           file_name = area1      file_name_obj = area1.Obj      am getting an "Datatype mismatch in expression" message in mapinfo, is it incorrect to have the above assignment? and what is the alternative method to refer to object files while file name is stored in a string variable? Thanks


title: What does the code RTTYP represent in the USA Tiger road files?

text: I'm trying to figure out what a particular code represents in the US Census Tiger Line files for roads. I downloaded roads for the US from here. I've searched the documentation (link to large pdf) and the web; all that I can find is the following:               All Roads (county-based) Record Layout   [top]     File Name: tl_2010_[state-county FIPS]_roads.shp     Field   Length  Type    Description     STATEFP     2   String  State FIPS code     COUNTYFP    3   String  County FIPS code     LINEARID    22  String  Linear feature identifier     FULLNAME    100 String  Concatenation of expanded text for prefix qualifier, prefix direction, prefix type, base name, suffix type, suffix direction, and suffix qualifier (as available) with a space between each expanded text field     RTTYP       1   String  Route type code     MTFCC       5   String  MAF/TIGER feature class code       I assume that `RTTYP 1 String Route type code` means it is some sort of classifier related to `Primary Road` or `Secondary Road` but I can't find that written down anywhere. If I examine a shapefile with this data I see the following: ![enter image description here](http://i.stack.imgur.com/lyLdY.jpg) I understand that `MTFCC` is a classifier from here :               24  S1100   Primary Road     25  R1011   Railroad Feature (Main, Spur, or Yard)     26  R1051   Carline, Streetcar Track, Monorail, Other Mass Transit Rail)     27  R1052   Cog Rail Line, Incline Rail Line, Tram     28  S1100   Primary Road     29  S1200   Secondary Road            so I can figure out what each road segment represents, but I'm curious to know what the `RTTYP` codes mean. I would love to find information about the surface material (whether it is asphalt or concrete, but I assume that information is not included).


title: Difference between Garmin eTrex Legend HCx and eTrex 30?

text: I have used etrex Legend HCx a lot and been satisfied. Now I've seen etrex 30 just for a while and it seems to me that the Legend is much more user friendly, good old menus with a lot of functionality, the display is great, whereas the display of etrex 30 feels a bit junky and cheap - it's kind of blurred.... I don't really understand why, as 30 is newer model. So why 30 seems to suck? But I've not tested the eTrex 30 in the field. **Can you please tell me what are the differences between these models?** Is eTrex 30 more precise? I found no precision specificatin for it. For HCx it should be 15 meters. I tried to find eTrex legend HCx on Garmin.com but I was not successful.


title: ArcGIS Server Object Extenstion REST Endpoint That Supports GET and POST

text: I have a Server Object Extension for ArcGIS Server 10.0, that implements IRESTRequestHandler. The extension contains many individual RestOperations. By default each RestOperation supports HTTP GETs when using the extension. In the RestOperation constructor, there is a parameter called postOnly, which is a boolean. If I set postOnly = true, then the endpoint will only support HTTP POSTs, and will throw an exception if you attempt to use a GET. I would like to create an endpoint that supports both GET and POST. If you interrogate some of the supported methods on a standard MapService like Find and Identify, you will notice these support both GET and POST:   * http://sampleserver1.arcgisonline.com/ArcGIS/rest/services/Demographics/ESRI_Census_USA/MapServer/identify   * http://sampleserver1.arcgisonline.com/ArcGIS/rest/services/Demographics/ESRI_Census_USA/MapServer/find Does anyone know of a way to implement this on a custom Server Object Extension?


title: map does not show up

text: I spent to much time working on a project, I saved it .. I'm sure now when I open it again .. arcmap does not display my map ..the layers there and everything else but the map's not displaying , It freaks me out .. It hapend to me more than one time I don't know why .. my map is so big.. could it be the reason ? any help any ideas .. It really drives me crazy I don't know what to do please help me here thanks


title: CORINE land cover 2000: coordinate reference system

text: The PDF in this link states that the CORINE land cover 2000 used ETSR89 as its coordinate reference system: http://www.eea.europa.eu/publications/technical_report_2006_7 I wanted to check I am reading this correctly: did CORINE 2000 use ETSR89 as its coordinate reference system? In addition, what would be an appropriate coordinate reference system to re- project CORINE 2000 over the UK?


title: Get centroid of SimpleFeatures in GeoTools8

text: I've got some SimpleFeatureCollections and now I want to get the Centroid of every SimpleFeature in the collection. Unfortunally I've only found a centroid method for Geometry. Is there any way to get the centroid?


title: IFieldChecker problem

text: I'm having some problems using `IfieldChecker`. I have this following code, which I'm trying to test permitted/not permitted fields, according to my workspace               public IEnumerable<IFieldError> GetFieldErrors()             {                  IFieldChecker checker = new FieldCheckerClass();             IEnumFieldError fieldErrors = null;                  checker.ValidateWorkspace = this.WorkspaceHandler.GetCurrentWorkspace;             checker.Validate(this.Fields,out fieldErrors,out this._ValidatedFields);                  if (fieldErrors == null)             {                 yield return null;             }             else             {                 fieldErrors.Reset();                 IFieldError fError = fieldErrors.Next();                      while (fError != null)                 {                     // debug only                       Console.WriteLine(fError.FieldIndex.ToString());                     Console.WriteLine(fError.FieldError.ToString());                          yield return fError;                     fError = fieldErrors.Next();                 }             }                  yield break;         }      Heres the test                   [TestMethod]         public void FieldBuilder_TestGetFieldErrors()         {             fBuilder = new FieldBuilder(this.wHandler);                  fBuilder.AddField(name, alias, esriFieldType.esriFieldTypeBlob, length, doublePrecision, doubleScale, false, false, null,                 false, null);                  Assert.AreEqual(2, fBuilder.GetFieldErrors().Count());         }      In my opinion, this should yield 1. The test completes, and it gives me a green light, but it should fail, since I'm testing with a bunch of values (1, 2, etc.); The workspace I'm creating is one on C:\, therefore, should not allow a Blob field type. Any ideas why this is working like this? **EDIT:** I'm working under the premise that after setting the workspace, ArcObjects would take care of validating it and telling me which types of fields are permitted or not. Example: FileSystemWorkspace - cannot have Blob field types; Or, does `IFieldChecker` does not work the way I'm thinking?


title: Is there a way to render/show ACAD drawing in ArcMap correctly?

text: I would like to extract only certain lines from an ACAD drawing but it seems that a lot of lines are being dropped when the dwg drawng is viewed and rendered in ArcMap. (see below) **Is it possible to view an autocad drawing in ArcMap the same way it is rendered in AutoCAD or TurboCAD.** Intersection as viewed in TurboCAD:![enter image description here](http://i.stack.imgur.com/5O6yw.jpg) Same file and the same intersection viewed in ArcMAP: ![enter image description here](http://i.stack.imgur.com/4siAd.jpg) In ArcMap when converted to FGDB Feature Class: ![enter image description here](http://i.stack.imgur.com/jsTxo.jpg)


title: Using "managed" intermediate data in ModelBuilder - ArcGIS - what are the benefits? vs the costs of not using in-line variable substitution

text: Do you have experience or examples of using "managed" data to ensure custom models are able to be shared? What are the costs / benefits of this versus setting the file path using in line variable substitution Most of my models at present use in line variable substitution e.g. %Scratch%\Var1b or %Scratch%\Variable%Value%_result Within the Model Scratch is then set as a model parameter I have found that in trying to make the models more robust and transferable (where there have been issues in sharing the custom tools) that making more of the intermediate data "managed" has solved some issues. However it causes other issues because if all data is set to be "managed" then the in_memory cannot be used and also in line variable substitution cannot be used to name files. I'd be grateful to hear your experience and solutions with using "managed" intermediate data to ensure models are transferable (and any implications regarding using in-line variables / in_memory work space). There is help note here that mentions that in some (unspecified) situations that not using "managed" intermediate data can cause generic error messages in Model Builder http://help.arcgis.com/en%20/arcgisdesktop/10.0/help/index.html#//002w0000005p000000 I am using ArcGIS 10.1, SP1, Desktop Advanced, on Win 8, OS 64 bit, Thanks,


title: Google Time slider - Recreated in MapInfo

text: Would it be possible to re-create the Google Earth Time Slider function with MapInfo? I assume through the use of MapBasic?


title: Is it possible to set a reference scale in QGIS?

text: It's a little tricky to pick label and feature styles (line thickness, colours, etc) in the map canvas because they don't reflect what I would see in my print composer which is at a different scale because of paper size and such. I haven't found any options in the QGIS options or project properties.


title: how to use the TRIM function in field calculator

text: I need to remove leading spaces in a field in my attribute table, but not how to do it from the field calculator


title: Convert from Geometry to SDE Binary storage type in SDE database

text: We have been battling performance issues between SQL 2012 and ArcGIS 10.1. One of the most common solution is to convert to SDE Binary storage. We know that this storage type is not as dynamic as the Geometry type. But we just opted to try to convert atleast one table in the database (assuming that it is possible). Problem is, we are not sure on the correct procedure to do this. Do we do it on the database side? or in the ArcGIS side? Thanks.


title: How to build attribute table from raster .tiff in ModelBuilder?

text: I have built a model that goes through a folder of raster files (TIFF), calling each time one file and running a routine. The problem I am facing is that imported raster files do not have an attribute table. I know how to create it manually for each one, but I would like to incorporate it to ModelBuilder. When I create manually the attribute table, ArcGIS creates some auxiliary files (XML, OVR and DBF). Does anyone know if this can be automatized in ModelBuilder? and how?


title: Doesn’t cal AGSSOM program to manage AGS mapservices from python script

text: We are using ArcGISServer 10.0 on a windows server. Because of database management we need to stop several ArcGISServer Services. For this purpose I use the AGSSOM program, which is called from the following Python program:               import os, subprocess          def manageMapServices(command, svc, svcType):         PathToAGSSOM = r"c:\Program Files\ArcGIS\Server10.0\AGSSOM_GESTAO\AGSSOM"         server = "AGSSOM"         p = subprocess.Popen(PathToAGSSOM + " " + server + " " + command + " " + svc + " " + svcType)          # Calling AGSSOM - stop servico     manageMapServices("-x", "GeoPlaneamento/AMAE2001", "MapServer")      When calling this program I receive the following error: *Could not find an available ArcGIS Engine or Desktop product to initialize. Could not connect to server. Make sure this user is a member of the AGSADMIN group. Details: The RPC server is unavailable. (Exception from HRESULT: 0x800706BA) AGSSOM [usage]: {server} [-s|start, -x|stop, -r|restart -list|list_services -listtypes|list_server_types, -describe|describe],-publish|publish from MX D path as second parameter [servicename (or _all_ )] {servicetype|default:MapServer, GeocodeServer, GPServer, GlobeServer}* This message points to problems with permissions, however the user has all the correct permissions. And when trying to call the AGSSOM program directly from the command prompt I do not receive any error. Does this anyone know what might be the problem causing this error?


title: Upstream to a river and tributary of a river - What is the underlying difference

text: What is the difference between a stream/ river upstream to a river A and a tributary of River A? I know that all tributaries are upstream to its river (A) and feeds into river A, but can there be a stream / river / lake upstream to river A but does not flow into it. In other words, if a waterbody is not connected to river A , can you still call it upstream to River A. Edit: I am interested in the semantics of these terms. This question is not related to any software, however modelling streams as either Tributaries of a river or just as upstream to a river is a GIS task (outside the scope of this question). If I say a particular stream "stream 1" is a tributary to a river A, then it automatically implies that "stream 1" is upstream to "river A". If so then modelling explicitly that "stream 1 is upstream to river A" is redundant, unless there is indeed a case where "stream 1 is upstream to river A, but its not a tributary of River A"


title: Error running mapserver export in quantum gis

text: I have a project in Quantum GIS that I need to export to a .map file using de plugin mapserver export. I fill all the forms but I get this error when running the app. Any ideas of what is happpening> Here is a copy of the error. Thanks in advance.               Traceback (most recent call last):   File "C:/PROGRA~2/Quantum GIS Lisboa/apps/qgis/./python/plugins\mapserver_export\mapserverexport.py", line 132, in ok_clicked         self.saveMapFile()   File "C:/PROGRA~2/Quantum GIS Lisboa/apps/qgis/./python/plugins\mapserver_export\mapserverexport.py", line 148, in saveMapFile         print "Creating exporter using '%s' and '%s'" % (self.dlg.ui.txtQgisFilePath.text(), self.dlg.ui.txtMapFilePath.text()) UnicodeEncodeError: 'ascii' codec can't encode character u'\xed' in position 49: ordinal not in range(128)          Versión de Python:     2.7.2 (default, Jun 12 2011, 15:08:59) [MSC v.1500 32 bit (Intel)]               Versión de QGIS:     1.8.0-Lisboa Lisboa, 6416f38          Ruta de Python: ['C:/Users/SergioVM/.qgis//python/plugins\\sextante', 'C:\\Users\\SergioVM\\.qgis\\python\\plugins\\elevation', 'C:\\Users\\SergioVM\\.qgis\\python\\plugins\\cswclient', 'C:/PROGRA~2/Quantum GIS Lisboa/apps/qgis/./python', 'C:/Users/SergioVM/.qgis//python', 'C:/Users/SergioVM/.qgis//python/plugins', 'C:/PROGRA~2/Quantum GIS Lisboa/apps/qgis/./python/plugins', 'C:\\Windows\\system32\\python27.zip', 'C:\\PROGRA~2\\Quantum GIS Lisboa\\apps\\Python27\\DLLs', 'C:\\PROGRA~2\\Quantum GIS Lisboa\\apps\\Python27\\lib', 'C:\\PROGRA~2\\Quantum GIS Lisboa\\apps\\Python27\\lib\\plat-win', 'C:\\PROGRA~2\\Quantum GIS Lisboa\\apps\\Python27\\lib\\lib-tk', 'C:\\PROGRA~2\\Quantum GIS Lisboa\\apps\\qgis\\bin', 'C:\\PROGRA~2\\Quantum GIS Lisboa\\apps\\Python27', 'C:\\PROGRA~2\\Quantum GIS Lisboa\\apps\\Python27\\lib\\site-packages', 'C:\\PROGRA~2\\Quantum GIS Lisboa\\apps\\Python27\\lib\\site-packages\\PIL', 'C:\\PROGRA~2\\Quantum GIS Lisboa\\apps\\Python27\\lib\\site-packages\\win32', 'C:\\PROGRA~2\\Quantum GIS Lisboa\\apps\\Python27\\lib\\site-packages\\win32\\lib', 'C:\\PROGRA~2\\Quantum GIS Lisboa\\apps\\Python27\\lib\\site-packages\\Pythonwin', 'C:\\PROGRA~2\\Quantum GIS Lisboa\\apps\\Python27\\lib\\site-packages\\wx-2.8-msw-unicode', '~/.qgis/python', 'C:\\Users\\SergioVM\\.qgis\\python\\plugins\\openlayersov\\logic', 'C:\\Users\\SergioVM\\.qgis\\python\\plugins\\openlayersov\\gui', '/usr/share/qgis/python', 'C:\\Users\\SergioVM\\.qgis\\python\\plugins\\QuickMultiAttributeEdit/forms', 'C:\\PROGRA~2\\Quantum GIS Lisboa\\apps\\qgis\\python\\plugins\\fTools\\tools', 'C:/SIG-SIXAOLA/Cartograf\xc3\xada/base_datos_final/Quantum_GIS_BD']


title: Accounting for Spatial Autocorrelation in ZINB Model?

text: I'm doing an analysis on disease in Chile. I'm using the rate as the `Dependent Variable` and both remote sensing and housing-type data as the `Independent Variables`. I have several questions regarding my methods:   1. The data appears to follow the zero-inflated negative binomial distribution. I checked with the countfit tool in Stata, and the data best fits this distribution. Is it ok to fit a ZINB model to continuous data?   2. I need to account for spatial autocorrelation in both the dependent, and the remote sensing data. I'm not sure exactly how to do this. After some reading, I thought that I needed to generate a spatial weights matrix (using ArcMap). If this is correct, how do I implement it? I generated one for the disease rate, and used it in Stata: `[pweight = weight]`, which greatly improved the model fit. However, I'm not sure if this is the correct usage. Can anyone give me some pointers, or readings, or anything? I really want to complete this project!


title: How to Find if set of Polylines are Connected?

text: I am using ArcGIS 9.3. I have a Polyline Feature class with attribute field as Street_Name. I want to Check if for say Street_Name='My Street' (say around 10 Features from the layer) all the polylines are Connected to each other. Edges in the Layer can be complex edge with multiple vertices. Any idea ?


title: Automatic bulk snapping with FOSS

text: I'm looking for a tool that would be able to automatically snap features closer than a given distance. Is there such a tool/plugin for QGIS/GRASS/OpenJump? Would it be possible to achieve that with PostGIS? I have a river network where some of tributaries don't touch the river they flow in and I need to get them snapped.


title: Perpendicular lines at start and end of a line

text: I have some lines. I want to draw perpendicular lines at the start and end of the lines, on both sides of the lines. Can somebody help me in Python? Also, what are the variable settings of the script?


title: Error on rendering linestring

text: I am using Mapnik 2.0.1rc0 on a Windows 7 with python 2.7 32bits inside a GeoDjango App using xml. The datasource come from spatialite 4.0. My tiling script can render polygon and point features correctly but linestring are just not displayed. The tiles are empty, without error. I can see the background of the tiles, but no trace of any lines when zooming in and out. Does anybody have an idea what could be the problem?


title: ArcGIS Add-In Project to Standard Library Project

text: I have been developing an ArcGIS 10 Add-In tool, but requirements changed halfway through development. Fortunately, this isn't changing any of the code, but I just need to be able to compile the project as a standard COM extension instead of an Add-In. I removed the `esriaddinx` file, and when I go to compile, it states that it is an invalid Add-In project. I haven't seen any properties within the project that declare it as an add-in project. Does anyone know how ESRI flags their projects to compile as add-ins, and is there a way to reverse it without creating a new project and importing files?


title: Can one edit an ArcSDE Oracle spatial view?

text: By 'edit' I mean edit the geometry in the feature class, not any attributes in the joined table. Unversioned. I really had it in my head that you can, as long as the objectid column is included in the view along with the shape column, some half-remembered reference. So I was forging ahead and decided to double-check if there was fine print, and now it seems I was whackadoodle in the first place. I've seen nothing to suggest it is supported, and seen in Administering ArcSDE for Oracle, page 225: "Spatial views cannot be edited through an ArcSDE client." But then I went ahead and tried it today, and I could edit a spatial view in ArcMap, but given the above I question how safe it is. Anybody have any experience to share?


title: mapnik error - fatal flex scanner internal error--end of buffer missed

text: I am using Mapnik 2.0.1rc0 on a Windows 7 with python 2.7 32bits inside a GeoDjango App. I often got this error: `fatal flex scanner internal error--end of buffer missed` that makes python to crash (sometimes python crash without error message) when zooming in and out of mapnik tiles rendered in OpenLayers map. The tiled data are stored in spatialite 4.0. I know such setup could provide a lot of source for this crash. Anybody has encountered this error before or may know where to look to make it works without crashing? **Edit:** GeoDjango is used as a Tile Map Server. I start by creating the xml in python then I use mapnik python bindings to generate the tiles. The data is already stored in spatialite. So each time I zoom in and out in the map, mapnik query spatialite like this:                   query = '(select ' + geometryField + ', "id" from "shapefile_feature" where' + ' shapefile_id=' + str(shapefile.id) + ') as geom'         datasource = mapnik.SQLite(file='C:\mygeosite\sqlite3\sql3.db',                                         table=query,                                         srid=4326,                                         geometry_field=geometryField)


title: Opening netcdf in arcmap as a raster with uneven coordinate spacing?

text: I have SMOS CATDS 2010-2013 daily data that is in EASE-grid netcdf format. I have 100 sites that I need specific daily soil moisture values for from this dataset. I have daily AAFC CATDS soil moisture data in tif format that I already extracted the soil moisture data I needed from, which I did using extract values to points in ArcMap based on the AAFC CATDS tifs and my in situ point shapefile sites. I need the same data from the SMOS CATDS NC files, but I don't know how to import these files so that the output raster data should have 586 rows by 1383 columns with a nominal cell size of 25km. Apparently, ESRI supports this format, so I just need to find a way to read in the data and have it output to a raster in this format without resampling the data (ie my output raster should be 586 by 1383, and I should achieve this without having to upscale or downscale the data). I tried using the multidimension toolbox, specifically netcdf to raster, but it wont let me because it says the spacing isnt consistent. So, the tool error help said to make it a feature layer and interpolate to raster, but I should be able to import it as a raster without doing this, as it changes the values and becomes inaccurate. I opened the nc files fine in Panoply, and it looks like how I want it to open in ArcMap, but I can't extract the points I need in it. I was thinking i might have to convert lon/lat to row/col on ease-grid, like it suggests on the CATDS website, http://www.cesbio.ups- tlse.fr/SMOS_blog/?page_id=749, but i thought it already was in this format. I was also hoping to avoid using MatLab in case there's a simpler way on ArcMap. I am not familiar with scripting aside from basic importing tools to IDLE and iterating them, so ncdump or others would not be really helpful, unless it was very guided/detailed. I also tried looking through http://www.unidata.ucar.edu/software/netcdf/docs/software.html#QGIS, a list of netcdf programs, but I want to keep it simple, and I'm most familiar with ArcMap. This is my first time posting, so sorry if this post is too long or unclear.


title: Multibeam bathymetry xyz data 53000000 features arcgis 9.3

text: Im working with an xyz file that has 53000000 points, i have only been able to see this in Microsoft access, because in txt its says that the amount of data its too large. So what i´ve tried to do:   1. Change .xyz to .txt ---> 3D Analyst tool, conversion, from file, ASCII 3D to Feature Class: So when i try to convert to natural neighbor i cant see the z data so i cannot do the interpolation. I also try using the Tin tool, but it just say that an error happened when it try to draw it.   2. I´ve opened .txt with Microsoft access (its when i saw the amount of data that i have), and i saved this table in .mdb, but i have no idea on how to use this to make an interpolation. Please i need some help to work with this data!


title: What is the Login to the ESRI FTP site for Arc Hydro Tools?

text: ESRI used to have a forum sticky that contained the login instructions to the FTP site for the Hydro tools. They have since deleted all of their forums. Does anyone know the login instructions for the ESRI Hydro FTP site?


title: How to choose a coordinate reference system

text: I have a jpeg image of a reserve's layout and would like to use it to plot animals movements on i.e. GPS points. I have no idea which coordinate reference system to choose during the "georeferencing" process.


title: How to Disable the add-in button when ArcMap is started?

text: Question 1:How to Disable the add-in button when the ArcMap started,then I use the SelectByRectangletool select features,then the button enabled! I am a fresh man about the add-in. Question 2:I made a menu as show ,As the figure shows, how to avoid the part of in red rectangle,because the menu do not have the child menu again,but it is show the arrow ,and part of the child menu. I means the tool’s skin is gray，and it is not activate ![question 1:tool gray， not activate](http://i.stack.imgur.com/sWgil.png) ![question 2 :My  custom toolbar with the menu click here ](http://i.stack.imgur.com/IHyHK.png) I am sorry about my poor engilsh here is the question2's config.esriaddinx file:               <ArcMap>       <Commands>         <Button id="Yangtze_University_GRP_Mult" class="Mult" message="new" caption="test" tip="tip" category="Add-In Controls" image="Images\1.png" />         <MultiItem id="ESRI_SelectionSample_MultiItem" caption="MultiItem AddIn" class="ZoomToLayerMultiItem" >         </MultiItem>       </Commands>       <Menus>         <!--main menu-->         <Menu id="Yangtze_University_GRP_My_Menu" caption="MainMenu" isRootMenu="true">           <Items>             <Menu refID="Yangtze_University_GRP_My_Menu_2" />             <Menu refID="Yangtze_University_GRP_My_Menu_1" />           </Items>         </Menu>         <!--submenu-->         <Menu id="Yangtze_University_GRP_My_Menu_2" caption="New file" isRootMenu="false">           <Items>             <!--<Button refID="Yangtze_University_GRP_Mult" />-->             <Menu refID="Yangtze_University_GRP_My_Menu_3" />             <Menu refID="Yangtze_University_GRP_My_Menu_4" />             <Menu refID="Yangtze_University_GRP_My_Menu_5" />             <MultiItem refID="ESRI_SelectionSample_MultiItem" />           </Items>         </Menu>         <!--third menu-->         <Menu id="Yangtze_University_GRP_My_Menu_3" caption="test1" isRootMenu="false" isShortcutMenu="true" />         <Menu id="Yangtze_University_GRP_My_Menu_4" caption="test2" isRootMenu="false" isShortcutMenu="true" />         <Menu id="Yangtze_University_GRP_My_Menu_5" caption="test3" isRootMenu="false" isShortcutMenu="true" />       </Menus>       <Toolbars>         <Toolbar id="Yangtze_University_GRP_My_Toolbar" caption="toolbar" showInitially="true">           <Items>             <Menu refID="Yangtze_University_GRP_My_Menu" />           </Items>         </Toolbar>       </Toolbars>     </ArcMap>


title: How to set default bad layers handler in the python plugin?

text: I am trying to write a python plugin for QGis. My python plugin will copy a template qgis project into data catalog with ESRI shapefiles. Because data catalogs can vary in its content there will be possibility that in the output catalog may be less layers than in the project template.   Is it possible to set default bad layer handler inside python plugin every time when I load programmatically copied project form output catalog?   I am new in programming qgis plugin and I suppose that I should use something like that:               proj = QgsProject.instance()     handler = proj.setBadLayerHandler() #I don't know what kind of arguments I should put in this method.      Thanks,


title: Is it possible to export an image that can be zoomed-in on?

text: I am using QGIS on a Mac. I would like to export my map so that the viewer can zoom in on it, just as I view the map on my QGIS canvas. I am not sure how to word what I want. If you see the map from the zoomed-out perspective, the symbols all bunch together, you need to be able to zoom in to see the details. Anyone have any idea if this is possible? I tried to export a high res image in the proint composer, but its not nearly enough resolution. Thanks!


title: ZoomTo is not working

text: I have a graphics layer. When the user selects an address from our software I want the map to ZoomTo that address. This works, but the user has to have an interaction for the map first. I thought this might be because the Map Extent was being reset, but when stepping though the code, it seems to Initialize first. Below is the code I am using. I separated the zooming code into another method just to see if it made a difference, but it did not. I also tried taking the extent off of "MyMap"               public void ShowMap(string account_num, string street, string city, string state, string zip)             {                 Locator FindlocatorTask = new Locator("http://geocode.arcgis.com/arcgis/rest/services/World/GeocodeServer");                 FindlocatorTask.AddressToLocationsCompleted += ShowMapLocatorTask_AccountToLocatonCompleted;                 FindlocatorTask.Failed += ShowMapQueryTask_Failed;                 AddressToLocationsParameters addressParams = new AddressToLocationsParameters();                 Dictionary<string, string> address = addressParams.Address;                 string formatAddress = " ";                 formatAddress = street + " " + city + ", " + state + " " + zip;                      OStreet = street;                           if (street != null)                 {                     RelatedFeatureGrid.Visibility = Visibility.Visible;                     OptionsPortal.Visibility = Visibility.Collapsed;                 }                      accountAddress = account_num; //used to find the associated meter                 RelatedAddressTextBox.Text = formatAddress;                      if (!string.IsNullOrEmpty(formatAddress))                     address.Add("Address", formatAddress);                           FindlocatorTask.AddressToLocationsAsync(addressParams);                  }                  private void ShowMapLocatorTask_AccountToLocatonCompleted(object sender, ESRI.ArcGIS.Client.Tasks.AddressToLocationsEventArgs args)             {                      CandidateListBox.Items.Clear();                 _candidateGraphicsLayer.ClearGraphics();                      try                 {                     AddressCandidate candidate = args.Results[0];                     Graphic graphic = new Graphic();                     graphic.Symbol = LayoutRoot.Resources["DefaultMarkerSymbol"] as ESRI.ArcGIS.Client.Symbols.Symbol;                     graphic.Geometry = _mercator.FromGeographic(new MapPoint(candidate.Location.X, candidate.Location.Y));                     graphic.Attributes.Add("Address", candidate.Address);                     string latlon = String.Format("{0}, {1}", candidate.Location.X, candidate.Location.Y);                     graphic.Attributes.Add("LatLon", latlon);                     _candidateGraphicsLayer.Graphics.Add(graphic);                          ZoomToLocation(graphic);                 }                 catch                 {                     MessageBox.Show("Address Cannot be Found", "Error", MessageBoxButton.OK);                 }             }                  public void ZoomToLocation(Graphic graphic)             {                      MapPoint candidatePoint = graphic.Geometry as MapPoint;                 double displaySize = MyMap.MinimumResolution * 30;                      ESRI.ArcGIS.Client.Geometry.Envelope displayExtent = new ESRI.ArcGIS.Client.Geometry.Envelope(                         candidatePoint.X - (displaySize / 2),                         candidatePoint.Y - (displaySize / 2),                         candidatePoint.X + (displaySize / 2),                         candidatePoint.Y + (displaySize / 2));                           MyMap.ZoomTo(displayExtent);             }      No errors are being thrown. Any help would be appreciated. I cannot think of any other problems that would cause this. Important Comment: I actually think I know what it is but I have no idea how I am going to handle it. It seems that the Initial Extent of MyMap is being set AFTER the Extent of the ZoomTo. The ZoomTo basically beats the initialize event, even though it is fired first. I have tried wrapping the code in an "if (MyMap.IsInitialized)", but it still is a no go. HOWEVER, if I put a messagebox in the ZoomTo method that forces the user to click ok, it works every time. I really don't want to have to do that. Thread.Sleep also does not work because it stops the Map from loading completely


title: Looking for a QGIS solution to make aerial imagery ArcSDE compatiable

text: Would like to avoid having to use 3d Analyst and use free raster tools in QGIS or elsewhere to make aerial rasters ArcSDE compatible. This 18 page work flow is what I am trying to use. Any suggestions on improving this workflow are also much appreciated.


title: Creating Irregular ellipses using SQL Spatial 2012 Geometry

text: I have stumbled across a problem in SQL Server 2012 - whereby I need to create different shape ellipses from the following data points: centre_x | centre_y | max_x | max_y | min_x | min_y I have tried the create curvepolygon approach but this doesn't work too well, especially if one side of the ellipse is different to the rest, it just seems to create two circles and merges them. Does anyone have an idea how I can tackle this approach? This is for a project where tree caponies need to be displayed, where the respective shape can be different for each feature. If there is a way whereby I normal ellipses can be made, this may suffice.


title: How can I calculate the bearing between two points in PostGIS?

text: How can I find the bearing between two points in PostGIS? Please specify in your answer whether or not the method producing a bearing on the spheroid, or a planar bearing.


title: How to draw features with openlayers?

text: i have a DrawFeature tool.               create= new OpenLayers.Control.DrawFeature(myVecLayer,     OpenLayers.Handler.Polygon);     app.mapPanel.map.addControl(create);      i have a store.                                           st = new GeoExt.data.FeatureStore({                                     layer: myVecLayer,                                     fields: [                                         {name: 'id', type: 'int'},                                         {name: 'mydata', type: 'String'},                                     ],                                      proxy: new GeoExt.data.ProtocolProxy({                                         protocol: new OpenLayers.Protocol.WFS({                                             version: "1.0.0",                                             srsName:"EPSG:900913",                                             url: "http://localhost:8080/geoserver/wfs",                                             featureType: "filedata_temp",                                             featureNS: "http://www.opengeospatial.net/cite",                                             filter:Ffilter                                         })                                     }),                                     autoLoad: true                                 });      I start drawning               for(var i=0;i<create_cont.length;i++){      create_cont[1].activate();     }      After i put attributes to feature and save result Openlayers.Strategy.Save().save() But its not saved. POST               <wfs:Transaction xmlns:wfs="http://www.opengis.net/wfs" service="WFS" version="1.1.0"  xsi:schemaLocation="http://www.opengis.net/wfs http://schemas.opengis.net/wfs/1.1.0/wfs.xsd" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance">      <wfs:Insert>       <feature:filedata_temp xmlns:feature="http://www.opengeospatial.net/cite">        <feature:the_geom>         <gml:Polygon xmlns:gml="http://www.opengis.net/gml" srsName="EPSG:900913">          <gml:exterior>           <gml:LinearRing>            <gml:posList>6756846.305626001 7731074.3336396 6756578.776027001 7730338.6272425 6757696.667564999 7730118.8707862 6758441.928590998 7730978.7873543 6756846.305626001 7731074.3336396</gml:posList>           </gml:LinearRing>          </gml:exterior>         </gml:Polygon>        </feature:the_geom>        <feature:id>10</feature:id>         <feature:mydata>10</feature:mydata>       </feature:filedata_temp>      </wfs:Insert>     </wfs:Transaction>      ANSWER               <ows:ExceptionReport version="1.0.0" xsi:schemaLocation="http://www.opengis.net/ows http://localhost:8080/geoserver/schemas/ows/1.0.0/owsExceptionReport.xsd" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:ows="http://www.opengis.net/ows">      <ows:Exception exceptionCode="InvalidParameterValue">       <ows:ExceptionText>Error performing insert: Error inserting features</ows:ExceptionText>      </ows:Exception>     </ows:ExceptionReport>      What i doing wrong?


Where can I download openlayers mobile api? Can i do everything with openlayers mobile like openlayers api? Thank you very much


title: Worsened image quality and false colours after georeferencing

text: When I georeference a map of vegetation types (about 20 colours, each representing a certain vegetation) the quality of the georeferenced map is really bad and some colours have changed (e.g. black becomes white), which makes it almost impossible to recognise certain types. I did the "usual" procedure: transformation type: linear, method: nearest neighbor. For source I tried jpeg & tiff, same bad result. Any suggestions how to improve the quality of the map and the originality of the map? Thanks!


title: Spatial Rendering using ArcGIS Restful APIs

text: Always when I read the documentation, I always hear the term Restful APIs and ArcGIS Server using Restful Service to render the mxd hosted on the server. I am a bit unclear about what exactly does Restful API mean? Also how does Servers like Geoservers and ArcGIS Server, render these spatial data? What is it that is different about these servers from normal servers we have? I know they are able to render spatial data, unlike normal servers, but what makes them have this capabilities and is it or not, related to Restful APIs?


title: How to fix undershoots in data converted from Microstation using ArcMap 10.0?

text: I have converted data from Microstation to ArcGIS. The converted data has undershoots at intersections, so I tried to snap lines using Spatial Adjusment. I've added some control points but ArcMap didn't fix them. **How can I fix the undershoots automatically?** I have many lines in my data, and fixing them manually would be too time consuming. ![enter image description here](http://i.stack.imgur.com/kIZmn.jpg)


title: MapInfo - WFS service with multiple SRS codes?

text: I've posted to the mapinfo forum (https://groups.google.com/forum/?fromgroups=#!topic/mapinfo-l/HfA2zna9Rao) but hoping for better response here. Will someone share an example configuration or documentation of a mapinfo WFS service which supports multiple SRS codes? The only way I've been able to do this is by creating a .TAB file for each projection.


title: "Horizontal" and vertical datums terminology

text: I have heard people referring to as a horizontal datum to models such as WGS84 and as a vertical datum to geoids. The former is primarily used for determining latitude and longitude but also the height (and thus X, Y Z can be calculated) above the ellipsoid while the latter one is used for determining height above an equal-potential surface only. I wonder why one would call "horizontal" a full 3D datum. The fact that you have a pure vertical datum (i.e. geoid) does not mean you have to call "horizontal" the remaining one. Here are two articles where you can find the term: http://geography.about.com/od/geographyintern/a/datums.htm http://www.esri.com/news/arcuser/0703/geoid1of3.html


title: How to Build Raster Attribute Table from .tif with pixel type of float or double?

text: I have a raster with the extension .tif. It's a subset of 87x87 pixel with values from some 6,66 to -4,27. If I right click on it there is no Attribute Table available. However, I can change the Symbology to Unique, Streched, Classified or Discrete Color and also classify manually, everything just fine. However, I need this attribute table, as I want to calculate the statistics of the values (not only the simple like mean and min, i got those, but a bit more advanced). The ArcGIS help site ('Build Raster Attribute Table (Data Management)') mentions, that: 'It is not possible to build a raster attribute table for a raster dataset that has a pixel type of floating point or double precision.' My values, however, ARE 32bit floating point. Does any of you have a clue what I can do?


title: Is there a function like XY-to-Line that calculates distance along a network for many XY pairs?

text: So:   I have a set of XY origins and a set of destinations. I know that XY to Line will give me the Euclidian distance between a given set of origins and destinations, and that an ODS cost matrix can return the distance to a given set of destinations for a given set of origins along a network. Unfortunately, an ODS cost matrix seems like a bad decision given the size of the data set I'm working with. So, I'm looking for a function that is exactly like XY to Line, but that will give me distance along a network between many sets of two points, instead of outputting Euclidian distance. Any ideas?    PS   Why does the size of my file factor into this? It seems to me that:   I am calculating the distance between (households of people on food stamps) to (retailers that accept food stamps). If I calculated the distance between every food stamp recipient in our sample (~55,000) to every retailer in our sample (700) for every month that we are concerned with (17), the calculations start to add up (~655 million). But if I only calculate the distance between the recipients (55,000) and the stores they actually shop at (about 5 per household), it gets a lot more managable, processing time-wise (~4.7 million).


title: How to implent a "maketrans" python function on field calculator

text: i'm traying to implent python string translate function at field calculatore http://www.tutorialspoint.com/python/string_translate.htm in order to fix hebrew imported from cad. here is the code -               def Ttxt (RefName):       from string import maketrans          intab = "asdfghjkl"       outtab = "שדגכעיחלך"       return trantab = maketrans(intab, outtab)          Ttxt( !RefName! )      Any sugestion about how to right this function ?


title: Failed to seed layer in GWC of GeoServer when Zoom Start is selected greater than 4 using a user-defined Grid Set

text: Everybody, I'm undergoing development on a application which is constructed by GeoServer and OpenLayers. Considering that there is a large set of layers (8 layer groupS and totally 121 layers) to display, we decided that user access the map in the cache way rather than dynamic generating. Because there are only 2 optional Grid Set in GeoWebCache from, I searched from google to create a new user-defined Grid Set 'EPSG:2326' which is aimed for HongKong, following are steps I did: 1, Newly added a Grid Set in 'geowebcache.xml' which is under path 'data_dir/gwc':                 <gridSets>         <gridSet>           <name>EPSG:2326</name>           <srs>             <number>2326</number>           </srs>           <extent>             <coords>               <double>806677.9759</double>               <double>802420.9858</double>               <double>8876817.5647</double>               <double>853384.7903</double>             </coords>           </extent>           <scaleDenominators>             <double>25000000</double>             <double>1000000</double>             <double>100000</double>             <double>25000</double>           </scaleDenominators>           <tileHeight>256</tileHeight>           <tileWidth>256</tileWidth>         </gridSet>       </gridSets>       <layers>      2, Added new optional Grid Set into layer group's XML config file which is under the path 'data_dir/sample_group' and named for 'sample_group.xml':               <entry key="GWC.gridSets">EPSG:4326,EPSG:900913,EPSG:2326</entry>      Finally, I go to the GWC Seed Form page to generate cache png, and seeding task can done successfully only on the condition that Stop Zoom is lower 4, otherwise will be failed and below is the related error message:               java.lang.IllegalStateException: Found no range bounds for z level 4: {0=[J@16a4cf2, 1=[J@1e1b0e7, 2=[J@2516dc, 3=[J@1000323}             at org.geowebcache.storage.TileRange.rangeBounds(TileRange.java:169)             at org.geowebcache.seed.SeedTask.tileCount(SeedTask.java:217)             at org.geowebcache.seed.SeedTask.doActionInternal(SeedTask.java:107)             at org.geowebcache.seed.GWCTask.doAction(GWCTask.java:76)             at org.geowebcache.seed.MTSeeder.call(MTSeeder.java:36)             at org.geowebcache.seed.MTSeeder.call(MTSeeder.java:25)             at java.util.concurrent.FutureTask$Sync.innerRun(FutureTask.java:303)             at java.util.concurrent.FutureTask.run(FutureTask.java:138)             at java.util.concurrent.ThreadPoolExecutor$Worker.runTask(ThreadPoolExecutor.java:886)             at java.util.concurrent.ThreadPoolExecutor$Worker.run(ThreadPoolExecutor.java:908)             at java.lang.Thread.run(Thread.java:619)      I guess that means no Bound Box for Zoom Level 4, I'm now confused and hope anybody out there give me some ideas? Thanks a lot.


title: How to Display Pie Chart when clicking on point feature in map?

text: I want to display pie charts from MS SQL Server database data. When I click on a point in the map, I wish to see a pie chart over the point. Can you please tell me how to implement such pie charts? Do you know any map site that using charts?


title: Problems configuring geoext tree

text: I am facing some difficulties configuring properly my geoext tree. I'll appreciate if anybody can help me. When I try to open the "images" folder, the node stays rounding, trying to load something forever. Here is my code:               var DNI2008_dec = new OpenLayers.Layer.WMS("2008_DNI_dec",             "http://localhost:8080/geoserver/wms", {                 layers: "personal:2008_DNI_dec",                 transparent: true,                 format: format,                 tiled: true,                 tilesOrigin : map.maxExtent.left + ',' + map.maxExtent.bottom             }, {                 buffer: 0,                 displayOutsideMaxExtent: true,                 ratio: 1,                 opacity: 1,                 visibility: false             }         );          map.addLayers([DNI2008_dec]);              var LayerTree = new Ext.tree.TreePanel({         title: 'Map Layers',         renderTo: 'layers',         loader: new Ext.tree.TreeLoader({                         applyLoader: false                     }),         root: {             nodeType: "async",             text: "Layers",             leaf: false,             expanded: true,             children: [                 new Ext.tree.TreeNode({                     nodeType: "folder",                     text: "images",                     leaf: false,                     children: [                            new GeoExt.tree.LayerNode({                                text: "test",                                nodetype: "gx_layer",                                layerStore: mapPanel.layers,                                layer: "2008_DNI_dec",                                leaf: true                            })                     ]                 })             ]         }     });      (CORRECTING A MIS-EDIT)


title: Openlayers Vector Event Popup with Textbox

text: The Openlayers map here: http://www.dassh.ac.uk/demonstrations/valmer/demo5.html has an event to add a vector layer, plus creates a textbox on the page below. Ideally I'd like the click event to add the marker, plus display a popup with a textbox/textarea where the user can fill in information on the location. This would then be submitted to the textbox below. Can anyone help? Thanks Dan


title: Connecting line features and determining the length of the longest line

text: I have a line feature (see image) representing a river that I created using the Stream_to_Feature tool. The attribute table contains several records representing different lines - the problem is the longest line (easily distinguishable visually) is not represented as a single line in the table, it is actually made up of many smaller lines. The lines appear to be touching, although they do not cross through each other. How can I merge these lines and then determine the length of the longest line using ArcObjects or manual methods that I can convert to ArcObjects? An even better solution would involve getting rid of all of the tributaries and leaving me with just the river channel as a line. ![Line feature - river](http://i.stack.imgur.com/YYZVV.png)


title: Question about quantum gis, and importing csv point data

text: I'm trying to add some point data from a csv file (it's in decimal degree lat lon) and overlay over road data that I already have added. The problem is that I can't get it to line up. I'm sure I'm making some kind of mistake with the coordinate system, but I can't figure out what it is. The data appears to the left and down miles away from my road data. When I import the same data in arcgis, I have to choose GCS_North_American_1983_HARN as my geographic coordinate system to make it line up, but that does not seem to be an option in quantum gis. Any suggestions?


title: How can I access second (or even third) level attributes in a vector layer and add them to the attribute table?

text: I got some .dxf file containing buildings I need for my map. The file loads nicely into QGIS and I converted sucessfully to shape. I still can find the Information I need (the buildings height) in the element identify box. However, it is organized as some kind of second level attribute and I'm running out of ideas how to access the information (screenshot below). ![Object Identifier box showing how the attributes are organized](http://i.stack.imgur.com/FhSnC.png) In this example the only entry in the attribute table is "myid". How can I access second (or even third) level attributes in a vector layer and add them to the attribute table, so I can work with them? Best regards! EDIT: The structure with the sub attributes (a... and (A... already exists in the .dxf. EDIT2: The Attribute table currently only consists of one single field named "myid". ![AttributeTable](http://i.stack.imgur.com/Q6N0d.png) EDIT3: Hm… The categorys (a… and (A… are shown in the Identify box when querying the .dxf. This is what I wanted to say. But of course, you may be right… ![Identify results for the original .dxf](http://i.stack.imgur.com/efSOC.png) If this is rather the perimiter than a height field, it is useless to me anyway. Bad luck. Do you by any means know of a better way to convert .dxf to .shp creating a height field?


title: Relationships in ArcGIS Server 9.3.1

text: I am working on a viewer in which I would like to show two related feature classes. If I select a feature of the first one, the related features of the second one should be selected too. I found a script on the Esri Resource Center: http://help.arcgis.com/en/webapi/javascript/arcgis/help/jssamples/fl_query_related.html My problem is that I can't get the relationship in de Mapserver. We run ArcGIS Server 9.3.1. I created a "Relationship class" in the database. This relation works fine in ArcMap but when I make a service of it, it disappears. Secondly I created a "Relate" in ArcMap, which worked fine there, but disappeared in the service. What is the solution to this problem? Is it possible to create a relationship with ArcGIS Server 9.3.1? The example server from the Resource Center uses ArcGIS SErver 10.


title: How many service instances is optimal?

text: What is the best number of instances? When I increase the number of instances the services will run faster?


title: `gdalbuildvrt` and selecting individual bands

text: I want to produce a set of VRTs (GDAL's internal virtual rasters) that comprise a re-arrangement of bands of a number of individual GeoTIFF files. Say I have my initial set of GeTIFF files, all having 4 bands:               file_01: band_01...band_04     file_02: band_01...band_04     file_03: band_01...band_04     file_04: band_01...band_04     file_05: band_01...band_04      And what I want to have is               vrt_01: file_01:band_01 file_02:band_01 file_03:band_01 file_04:band_01 file_05:band_01     vrt_02: file_01:band_02 file_02:band_02 file_03:band_02 file_04:band_02 file_05:band_02     vrt_03: file_01:band_03 file_03:band_03 file_03:band_04 file_04:band_03 file_05:band_03     vrt_04: file_01:band_04 file_04:band_04 file_03:band_04 file_04:band_04 file_05:band_04      My initial thoughts were:   1. Use `gdalbuildvrt`, creating an "file of files", but I don't think you can specify a particular band, just filenames   2. Use `gdal_translate` to create VRTs of the individual bands for each file, then proceed as in (1) above, but that looks like far too much effort. I'd like to avoid creating the XML file myself because I find that that is quite error prone and requires quite a bit of debugging (they can all have different extents, etc etc etc)


title: Labels on the QGIS WMS layers are not drawn according to the style definition

text: I have an issue regarding the styling of labels of WMS layers using QGIS server. I can't seem to get them to display how I want. What I did:   * create a QGIS project (2.2 on Mac) with a layer from a database   * style the layer, including the labels on QGIS desktop   * save project on the server (Ubuntu) to create the WMS What I get:   * my WMS layer, nicely showing up   * except for the labels that seem to use a default symbolization, whatever the changes I make in the project (understand big black labels) What I wished:   * control the size and maybe also the font of the labels. Is there a nifty parameter that I forgot somewhere or is the labeling engine not working propoperly with WMS? Any hints or similar experiences welcome.


title: gdal_rasterize -at function still won't capture covered pixels properly

text: I'm currently using the gdal_rasterize tool in bash. My code looks like this:               gdal_rasterize -l $(basename $file .shp) -burn 1  -a_nodata 0 -te -180 -90 180 90 -tr 0.1 0.01 -at -ot Byte $file $OUTPUT/rast_2.tif      I was hoping that the `-at` function would give a '1' burn value to any pixel touched by the input shapefile. However, this is what the comparison looks like: ![Areas not rasterized properly](http://i.stack.imgur.com/kEQ0Z.png) Is this the a problem with the -at ALL_TOUCHED function or should I be doing something differently to ensure that all areas are accounted for? Any help much appreciated


title: How to create lines from points pair coordinates with arcpy?

text: I've got some points pair coordinates (start and end points) which I have to transform into lines. Until now, I used an append of both the coordinates in a `pippo.Point()`, a `pippo.CalculateGeometry()` to define the geometry of each piont, and `pippo.append(defined geometry)` to identify the pair of points, and then PointsToLine to obtain my line. Because this is quite time expensive to do for hundreds of lines, I'd really like to know if there's a shorter way to do this, e.g. place starting and ending point of each line in different fields of a single table and import lines directly without passing for points geometry.


title: How to configure Proxy.cgi file for Geoserver On Tomcat?

text: I have a web map application in my local machine(windows xp and dont have IIS) and I have postgis data in Geoserver on Tomcat (Ubuntu). To get the WMSgetfeatureInfo, I have configured the proxy.cgi file in Geoserver (usr/lip/cgi-bin/proxy.cgi) and added my local ip address in allowedHosts = ['xxx.xxx.x.xxx', 'www.openlayers.org', 'openlayers.org', 'labs.metacarta.com', 'world.freemap.in', 'prototype.openmnnd.org', 'geo.openplans.org', 'sigma.openplans.org', 'demo.opengeo.org', 'www.openstreetmap.org', 'sample.azavea.com', 'v2.suite.opengeo.org', 'v-swe.uni-muenster.de:8080', 'vmap0.tiles.osgeo.org','www.openrouteservice.org'] and added the below line in my script.               OpenLayers.ProxyHost = "/cgi-bin/proxy.cgi?url=";      but it throws an error in chrome console:               GET localhost:3890/usr/lib/cgi-bin/….…T%3Dapplication%252Fvnd.ogc.gml%26SRS%3DEPSG%253A27700%26X%3D682%26Y%3D284 404 (Not Found)       I dont know whether its right or wrong. Proxy.cgi should be copied in Geoserver or my machine? in which path? added my IP into allowedhosts is correct? Please anyone suggest me to fix this issue


title: What is proper bearing for plat ties to monuments?

text: I had a 'discussion' with an engineer about the proper way to describe a bearing on a plat (survey of land) that ties the plat to a known survey monument. He stated that it should be from the point of beginning (POB) of the plat to the monument, while I argued that you should start from the known point and go to the POB. His argument was that the POB is a known point. It may be to the surveyor, but the rest of us don't know the location of the nail or rod the surveyor found. I know it's just a matter of reversing directions, not a big deal, but is there a standard for this? Thanks.


title: Intersection of multi-polygons in postgis

text: How can I get the intersection of 3 multi-polygons in postgis? I have three geometry collection and I want a fourth with the result of intersection of the other three.               GEOMETRYCOLLECTION(MULTIPOLYGON(((-49.64435197 -29.23896786187647,-49.68688498 -29.290587315131337, etc etc)))          GEOMETRYCOLLECTION(MULTIPOLYGON(((-38.7333 -15.635996265940811,-38.706199999999995 -15.699596252308577,-38.6708 -15.749496241625828,-38.650999999999996 -15.790396232878372,etc)))          GEOMETRYCOLLECTION(MULTIPOLYGON(((-38.7333 -15.635996265940811,-38.706199999999995 -15.699596252308577,etc)))          GEOMETRYCOLLECTION(MULTIPOLYGON(((result)))      Thanks


title: Converting WKT to WKB

text: How can I convert a WKT string to a WKB byte[] in Java using postgis functions? I need to save that WKT string in a postgis database (geometry column. Thanks


title: export from .shp to .gpx and .kml fails

text: I've just upgraded to QGIS Dufour and I am having some problems with the exportation of a simple track from .shp to .gpx and .kml format. I imported from mi smartphone a .gpx track (SR WGS84 EPSG:4326), i converted it into .shp (EPSG:3857 in order to have it correctly projected over the google maps layer); I modified it and I finally need it to be exported into .gpx and .kml formats. For .gpx: I did the "Save as", chose GPX format, set the RS to EPSG:4326, wrote "GPX_USE_EXTENSIONS=YES" in the Data Source field and "FORCE_GPX_TRACK=YES" in the Layer field. Qgis says it all went ok (and it loads the track in the viewer and it looks fine), but when I load the .gpx track into my gps device I only see a straight line into the sea. For .kml: I did the "Save as", chose KML format, set the RS to EPSG:4326. Qgis says it all went ok, but when I load it in the viewer it looks just empty and in google earth I can see a huge red cross over Africa... and I was looking for a road in Sardinia. Can anyone help me, please? Thanks a lot!


title: How can I cache reverse-geocode results?

text: I'm implementing an Android application that needs to do some reverse- geocoding to resolve a geocoordinate to a city. What is the best way to cache this information on the smartphone so that I don't have to always call an online server to perform this reverse-geocoding?


title: How to copy table in ArcGIS for Desktop and keep same feature ID (objectID)?

text: What geoprocessing tool I can use to copy a table from SDE geodatabase to personal geodatabase, keeping the feature ID (objectID)?


title: How can I make a point buffer map on the web?

text: I am looking for a tool like Google Maps Engine Lite or Mapbox where I can draw a simple (buffer) circle around a point with a known diameter. **Are there any quick and easy web mapping tools to draw a single circle and share?** For example, someone on a forums says, "Show me what a 10 mile ban on kite-flying would look like?"


title: Euclidean Distance Function performed on unprojected data

text: I am looking to perform a simple multi-criteria analysis on a data set comprised of two criterion: occurence of wildfires and road accessibility. First, the project CRS was enabled with "on the fly" projection selected, as well as "meters" selected as the distance measure. Second the data was rasterized, then the euclidean distance was calculated using r.grow.distance in the GRASS functionality. The results show distance in distance degrees as opposed to meters, even though the project properties have been specified in meters. Is there anyway to perform this analysis so that the results are in meters?.


title: NoData Value = 0, problem when creating NDVI, Erdas

text: I try to create an NDVI from Landsat and ASTER satellite images using Erdas Imagine. In Erdas, the satellite images have a nodata value of 0. After creating an NDVI, the nodata value is still 0, but obviously there are also file pixels with the value 0 now (NDVI ranges from -1 to +1). So if I reset the NoData value to 0 (it is not done automatically), some pixels in my raster are missing since they are now classified as "nodata" Any suggestions?


title: List index out of bounds [-1] when making bioclimatic variables in DIVA GIS

text: I am working with climate data and I would like to make 19 bioclimatic variables which I can use for ecolgical niche modeling in Maxent. I am following this tutorial http://www.ccafs- climate.org/downloads/docs/Ramirez_Bueno- Cabrera_2009_tutorial_bcvars_creation.pdf (google: bioclimatic variables + Maxent) After I create the CLM files in DIVA-GIS, I go to data>climate>map and I receive this error: "List index out of bounds [-1]". Has anyone run into such an error and do you know the solution?


title: Save As in QGIS 2.4

text: I am trying to save edits to a raster file in QGIS, but;   1. When you edit a raster file and select "Save As"; no output file is generated.    2. My version is in Spanish; however "Save as raster layer" is displayed. Just as "Output" How do I properly save/export a raster in QGIS?


title: Census data 5-year estimate for the United States

text: I am looking to download demographic data 5-year estimate based on census tract for the whole United States. The census website has the information, however, it does not let me to select the whole country based on census tracts along with the variables. Could you please advise?


title: Error executing tool.: Layer "layer1" Unable to connect to map server

text: I am using the Configurable Javascript Viewer, and have only changed a couple of the configuration settings. One change I made was to use my own Export Web Map Task. In viewer.js, I changed:               printTaskURL: 'http://maps.evansvillegis.com/arcgis_server/rest/services/Utilities/PrintingTools/GPServer/Export%20Web%20Map%20Task'      I then saw this note in the readme.md in the proxy directory: > You will more than likely need a proxy page for printing and other large > cross domain requests. So, I set that up using the provided proxypage_net.zip files, proxy.ashx and proxy.config. In proxy.config, I added:               <serverUrl url="http://maps.evansvillegis.com"                matchAll="true"></serverUrl>      Then I made sure my viewer.js file pointed to the current proxyUrl:               esriConfig.defaults.io.proxyUrl = 'proxypage_net/proxy.ashx';      However, when I try to print, I am getting the following error in the response:               {"error":{"code":400,"message":"Unable to complete operation.","details":["Error executing tool.: Layer \"layer1\": Unable to connect to map server at http://maps.evansvillegis.com/arcgis_server/rest/services/BASEMAPS/BASEMAP_GOOGLE_WGS84_NO_FOOTPRINTS_BW/MapServer.\nFailed to execute (Export Web Map).\nFailed to execute (Export Web Map Task)."]}}      I have seen related posts: Secure Print Service Failed to Execute ArcGIS Online printing Unable to connect to map server? However, the provided information (mind you, limited) was not much help. **The 400 error - is this because my proxy is not setup correctly, or is there something wrong with my PrintTask service?** This is the default PrintTask that is configured with ArcGIS Server 10.2.1. The only other change I have made, is ensuring IIS users and the ArcGIS Server Account can access the ArcGIS Server directories. Any help appreciated! Viewer: http://maps.evansvillegis.com/viewer


title: Converting ASCII grid/NTF to SRTM/DEM

text: I want to use the Ordnance Survey Land-Form Panorama product in some software called Maperitive, which can read custom SRTM DEM files. The data is provided in the following formats: DXF [Contours] ASCII grid [DTM] NTF [DTM] https://www.ordnancesurvey.co.uk/opendatadownload/products.html (search for land-form) Can anyone recommend a free tool or toolchain (runs on Linux or Windows) that I can use?


title: Anyone built a custom GPFunction in ArcGIS10 and having closing issues?

text: I have a custom GP function in ArcGIS 10 written in C#. When the tool is run it takes ArcGIS about 5 minutes to close (as in the application closes but in Task manager ArcMap.exe and ArcSOCP.exe are still running). When closed ArcMap.exe will increase memory too about 170 - 250 MB. To ensure it wasn't my code I created a blank GP Function with no code in the Execute function and return and empty paramarray value. Ran the tool and ArcMap still took ages too close. I then opened ArcMap again and launched the tool but this time pressed cancel instead of go. ArcMap closed fine. So as soon the tool enters the execute method it seems to then take ages to close when the Execute method finishes. If I run one of the out of the box GP tools ArcMap closes fine. Anyone encountered the same problem?


title: Cannot replicate gdaltransform example

text: In the description of gdaltransform this example is given: > gdaltransform -s_srs EPSG:28992 -t_srs EPSG:31370 177502 311865   >  Produces the following output in meters in the "Belge 1972 / Belgian > Lambert 72" projection:   >  244510.77404604 166154.532871342 -1046.79270555763 If I try this I get               gdaltransform -s_srs EPSG:28992 -t_srs EPSG:31370 177502 311865       244422.238366377 165970.245450285 0      If I write the coordinates into a single point-shape, same result when applying ogr2ogr               ogr2ogr -s_srs EPSG:28992 -t_srs EPSG:31370 out31370.shp in28992.shp      Any ideas what I'm doing wrong?


title: What is the best way to mosaic 2 million files?

text: Trying to create 1 mosaic dataset from 2 million files. Would it be better to break the files into regions and then create regional mosaic datasets and then make a mosaic dataset of the mosaic datasets (does that make any sense?)? Or just slug it out and do one dataset of 2 million files? BTW the datatype is CIB


title: What are some examples of model-based upscaling of raster data?

text: Remotely-sensed raster data may be too fine-grained for many purposes. A typical approach to coarsen the grain is to upscale by reducing the resolution. If the data are continuous (e.g elevation, surface temperature) cells on the raster can be aggregated (and therefore reduced in number) by using a mean function. Or if they are discrete (e.g. vegetation) they may be aggregated using the mode of adjacent cells. Moving windows are a second example of upscaling. Here a "window" is slid across the raster and the cells under the window are replaced using some function of their values (e.g. using a mean or mode). I believe that _model-based_ approaches to upscaling will simplify a raster by adopting a more complex model of the raster (i.e. the cells and their adjacencies). Can anyone clarify my definition of _model-based_ approaches, and provide an example? Do the resulting "cells" even need to be regular, could they vary in shape and size? EDIT: To be more specific about what I mean by model-based upscaling. Here is an example of what I think _may_ be considered a model-based approach. Consider using a mathematical graph to represent discrete features on a raster. Then use some method to identify components of these features (i.e. groups of vertices). A Voronoi tessellation, for example, using the components as generators, could then be used to produce an irregular tiling that represents an upscaling of the original raster. ... I could give more detail here, but opt not to, as I am interested in mostly whether or not this would be considered a model-based approach to upscaling ... or alternatively how such operations could be described in generic terms. Is this merely a way of doing interpolation (as defined in an answer below)?


title: How to reproject shapefile using ogr2ogr

text: I have shapefile in EPSG:2834 and I want to reproject it into EPSG:4326. I am trying to use `ogr2ogr` in Windows; however, I get this failure output:               C:\Users\Nik\Desktop\Logan_Lines>ogr2ogr -t_srs epsg:4326 spans_4326 spans     FAILURE:     Unable to open datasource `spans' with the following drivers.       -> ESRI Shapefile       -> MapInfo File       -> UK .NTF       -> SDTS       -> TIGER       -> S57       -> DGN       -> VRT       -> REC       -> Memory       -> BNA       -> CSV       -> NAS       -> GML       -> GPX       -> KML       -> GeoJSON       -> Interlis 1       -> Interlis 2       -> GMT       -> SQLite       -> ODBC       -> PGeo       -> OGDI       -> PostgreSQL       -> MySQL       -> XPlane       -> AVCBin       -> AVCE00       -> DXF       -> Geoconcept       -> GeoRSS       -> GPSTrackMaker       -> VFK          C:\Users\Nik\Desktop\Logan_Lines>      I am hoping someone can tell me what I am doing wrong; I have little experience with this tool.


title: Rounded Grid - Is it possible?

text: I was wondering if it's possible to make a layout, with rounded measured grids, in zoom maps. I'm using ArcMap and rounding the corners of data frame on layout mode to give a magnifying glass effect over the map. Here's a picture example of what I've already got. ![enter image description here](http://i.stack.imgur.com/WKLD5.png) Thanks.! Rounded dataframe ArcMap