title: Why does Merge then Dissolve return ERROR 000732: Input Features: Dataset ... does not exist or is not supported?

text: I'm not sure why I'm running into this problem. I thought it might have something to do with the merged file being locked or open, but I don't see a "LOCK" file being created in Windows Explorer. I'm simply trying to Merge a bunch of files and then Dissolve the merged file. Here's my script:               workspaces = arcpy.ListWorkspaces("*", "Folder")     for workspace in workspaces:              arcpy.env.workspace = workspace         shps = arcpy.ListFeatureClasses("*", "ALL")              ## Process: Merge Layers         descWorkspace = arcpy.Describe(workspace)         workspaceName = descWorkspace.baseName              mergeOutput = workspaceName+"_merge"         arcpy.Merge_management(shps, mergeOutput)              ## Process: Dissolving Features              <--HERE IS WHERE I'M HAVING PROBLEMS         disInput = mergeOutput         disOutput = mergeOutput+"_dis"         arcpy.Dissolve_management(disInput, disOutput)      This is the error that I'm getting after the merge runs successfully: `ExecuteError: Failed to execute. Parameters are not valid. ERROR 000732: Input Features: Dataset Bus04_merge does not exist or is not supported` I'm obviously missing a step here... (or using incorrect parameters, but I can't figure out in what way) Thanks for your help!


title: GRASS python script runs twice

text: Is there a reason that every script I have written for GRASS in python runs twice? I start them using the GUI script launcher tool, and they all begin with:               #!/usr/bin/python     import grass.script as grass      So nothing fancy. Still, every script always runs through twice, first one run through and then another, with the first run through locking the GUI until it finishes, and the second unlocking the GUI & displaying status update bars, etc. I've confirmed that they're running twice both by using log files & with ps -e. This happens both in GRASS 7.0 and in GRASS 6.4.2 (Ubuntu 12.04). Am I missing something obvious? If not I'll go ahead & submit a bug ticket.


title: how to search a vector/wfs layer

text: I'm using Geoserver/PostGIS/openLayers I've succesfully set a vector layer. And now I'm trying to create a search functionality. But something goes wrong. First I set up the layer and add some basic style for the features. So far, everything is ok. Then I create a function, using javascript, which comes after the layer. This function gets triggered after the user completes a HTML form and hits the "Search" button. Here's what's inside the function:               //check for empty fields in the form-if not pass the value in a var           if (document.getElementById("typeselect").value!=" ") {var filter_a = document.getElementById("typeselect").value;}          if (document.getElementById("eraselect").value!=" ") {var filter_b = document.getElementById("eraselect").value;}          if (document.getElementById("nm").value!="") { var filter_c = document.getElementById("nm").value;}          if (document.getElementById("descr").value!="") { var filter_d = document.getElementById("descr").value;}          if (document.getElementById("chrono").value!="") { var filter_e = document.getElementById("chrono").value;}          //create a style           var vsd = new OpenLayers.Style();          //put together the final filter, using the vars from above           //filter_a and filter_b are numbers           //filter_c to filter_e are texts           var filter_final = new OpenLayers.Filter.Logical({ type: OpenLayers.Filter.Logical.AND,           filters: [ new OpenLayers.Filter.Comparison({           type: OpenLayers.Filter.Comparison.EQUAL_TO, property: "p_type", value: filter_a }),           new OpenLayers.Filter.Comparison({ type: OpenLayers.Filter.Comparison.EQUAL_TO, property: "p_era", value: filter_b }),           new OpenLayers.Filter.Comparison({ type: OpenLayers.Filter.Comparison.LIKE, property: "p_name", value: filter_c }),           new OpenLayers.Filter.Comparison({ type: OpenLayers.Filter.Comparison.LIKE, property: "p_descr", value: filter_d }),           new OpenLayers.Filter.Comparison({ type: OpenLayers.Filter.Comparison.LIKE, property: "p_chron", value: filter_e }) ] });                    //add the filter as rule to the style           vsd.addRules([filter_final]);                    //apply the style as default           var vst = new OpenLayers.StyleMap({'default': vsd});                    //apply the new style to the layer (named pins)           pins.styleMap = vst;      This is not working at all. I've tried editing my code, I don't know what to do, I've searched for an answer and nothing really helped me. Please advise... Thanks. EDIT I also added this line at the end and still nothing `pins.refresh();` Do I need to use the WFS protocol, combined with the `read()` ? EDIT #2 I've replaced `pins.refresh();` with `pins.refresh({force:true});` and all the features dissappear. I also tried something like this and still nothing. This is a torture...


title: Getting coordinates of pixels

text: I want to compute an area of certain object on a picture using Monte Carlo method. The picture shows some shape on white background. To do so I figured out that I'd get the coordinates of non-white points and then I'd launch the rest of the routine in e.g. Excel. My problem however is that I don't know about any application or software that performs what I want - i.e. getting coordinates of pixels of some specified colour. Can someone suggest some application or software which can do what I want?


title: Get rid of raster grid cells that are less than 1

text: I'm using QGIS to work with raster data. I have a raster dataset of a country where each grid cell is assigned a number 1-100. However, the raster is in the format of a rectangle based on the extent of the file, so all of the cells outside of the borders of the country are assigned a 0. I cannot figure out how to create a GeoTIFF that excludes these cells that have a value less than 1. Any ideas? Thanks!


title: How to format spreadsheet to input point data for later import to PostgreSQL

text: I want to create a spreadsheet to hold data from our cemetery. I will use the spreadsheet to get the data into a PostgreSQL database. What I dont know yet is the format that the point (positional) data will be in. Is there a 'standard' format for it? Do I need one column that will hold the lat and long positional data together or do I need two columns one for the lat dat and one for the long data. If two columns is there a 'standard' for which comes first in the spreadsheet and should the positional data be at the beginning, left side of spreadsheet, or the end, right side of the spreadsheet? Any other advice would be useful and appreciated.


title: Determining if Polygon Contains Negative Space using ArcGIS SOAP Web Service?

text: I am running queries against a polygon layer using the ArcGIS SOAP Web Service. A user selects a point within that layer and the resulting feature set contains PolygonN objects. Occasionally, the polygons that come back contain negative space. These are also known as donut or island polygons. Here is an example: ![enter image description here](http://i.stack.imgur.com/AnvKy.png) In this case, I would only want the geometry contained in the purple area, not the green area. When the user selects a point within the purple area, the resulting PolygonN object contains 2 rings, representing the outer and inner areas, as it should. The inner ring, as far as I can tell, has no way to indicate that it is actually negative space. So the resulting polygon looks like it has 2 solid shapes, when in fact it is only one shape with a hole in the center. Does anyone know how I can determine that a polygon's ring is actually negative space?


title: Having Trouble with ST_Buffer and ST_Union

text: I am trying to build a coverage map, and measure how well a given set of points covers the target area. I started by writing a query like this:               SELECT cartodb_id,        ST_Transform(ST_Buffer(the_geom::geography,400)::geometry, 3857) as the_geom_webmercator     FROM set_of_points      That gave me a nice map of the areas covered by the set_of_points. But now I want to combine those buffers into a set of polygons where the overlapping buffers are merged into a single polygon. What's the best way to do this?


title: Linear Referencing Methods: Station Equation origins

text: DEAR PEOPLE FROM THE FUTURE: Here's what we've figured out so far… I am working in the pipelines field, and I get in touch with the concept of the linear referencing methods, specifically with station equation method. I had to do research about this method, I'm quite confident about my understanding of the topic and the concept associated to it, but a set questions came to my mind: Where did this topic originate? Who proposed it? When was it proposed? After an intensive Google search I have not found a clue about the history of station equation method. I got a lot of technical documents that describe the concepts and how this concept works, but not its origins. Any book or paper sugestion will be appreciated.


title: CenterPoint of Polygon in ArcEngine .Net 10

text: I am currently trying to get the center point of a polygon. In GIS, i can see the "x" mark in the center of the polygon, but I do not know how to get it from the IPolygon object. I tried getting the enevelope and getting the midpoint from that but that didn't work out very well, as getting the envelope does not always ensure to get the midpoitn where marked "x". Thanks and Regards, Kev


title: How do make a CSV file visible in QGIS?

text: I am using QGIS.   I have a .csv file containing 3 pieces of information on each row; name, longitude and latitude.   I would like to turn each row of the .csv file into a vector point. The file will load into QGIS when I open it as a vector layer, but whatever I do, the screen remains blank. Here is the metadata for the file: General:   Storage type of this layer: CSV   Source for this layer: G:/VECTOR FILES for qgis/Maritime Traffic Visualisation/ships2.csv   Geometry type of the features in this layer: No geometry   The number of features in this layer: 20189   Editing capabilities of this layer: Add Features   Extents:   In layer spatial reference system units : unknown extent   Layer Spatial Reference System: +proj=longlat +datum=WGS84 +no_defs I'm fairly new to GIS and this is probably something very simple, but your help would be appreciated. * * * Here are the first 10 rows:               name    lon lat     NAPOLEON I  -79.5231    8.66208     ZZ-9 VERTROUWEN 3.91968 51.64689     MIN LONG YU 66628   118.4183    24.12905     ZZ-5 CONSTANT PETRUS    3.92033 51.64691     MIN JIN YU 8401 118.574 24.33855     MIN LONG YU 60209   118.4849    24.15916     JINTANGYU03817  117.7681    39.07499     GOLDEN SCEPTRE PD50 -1.7716 57.5045     MIN  LONG  YU  62503    118.4402    24.27062     THILANT 4.43238 51.88438     VIKING  4.431762    51.23607     TYCHA   4.12002 51.9541     ILHAN AGA   27.981  43.199     ANKIE   -0.7144283  53.60629


title: qgis2 html annotation docs?

text: The HTML annotations in qgis2 seem to fit the bill of what I'd like to do: have a full-custom data-based drawing in a callout. Problem: I can't find documentation anywhere on how to write the html file to incorporate placeholders to take data values from the annotated object. Web searching shows a few discussions that say something about placeholders, which is the only evidence so far that it could actually do the trick. Can someone please point me to documentation on writing the html file used by html annotations? Maybe I'm just missing something obvious? Thanks


title: How do you install updates/patches for bug fixes in QGIS?

text: I'm trying to fix this bug in Qgis 2.0. I suppose the code, under the "Files Changed' tab needs to be copied/installed somewhere? Could you please advise on how to do so? I cannot find any help page on this - but may have been Google-ing the wrong terms. More broadly (but not my immediate concern) - why are these changes just not corrected in the original source code?


title: how can i use the postgresql's Geometric Types

text: How can i implement Geometric Types in postgresql db.. Postgresql db is hosted in windows os.. and how can i get whether point is in circle or not? getting errors when i am using predefined functions... Is there anything do i need to install further .. i am using postgresql 9.1 version. I have seen PostGIS for Linux... what is that .. how can will it be useful? Thanks in advance Aj..


title: Why is the imported TIFF file distorted with noise in Quantum-GIS?

text: I've a TIFF-file. Loading the raster image into QGIS does not display correctly. (Using the `Add raster layer` button.) It seems to be distorted with noise. See image below: ![qgis \(l\), image viewer \(r\)](http://i.stack.imgur.com/AtE9h.jpg) The source of the TIFF file is a 300 dpi scanned document. Other scanned documents display correctly. Projection doesn't matter as the image is not georeferenced yet. I'm using qgis-1.7.0 in a sabayon-linux environment. What's the source of this problem? Did anyone else encounter this before? Thanks for any hints or experience on this topic.


title: Algorithm to place the minimum of rectangles to cover a polygon

text: I'm trying to cover a polygon with the minimum set of rectangles like explained in the picture below, the used rectangles have the same size a*b (a != b) , the rectangle could be rotated by 90°(otherwise the square size could be a*b or b*a). in the best case if we could have some overlap area, on the edge of adjacent rectangles, it would be great. i would like to know if there's an algorithm could help in this problem. any help is welcomed. ![enter image description here](http://i.stack.imgur.com/5noFh.png) I'm using arcpy in arcgis 10.1


title: Automatically updating a table column with the name of table

text: Does anyone know a straightforward way in MapBasic to automatically update a table column with the name of the table the column exists in? Going from the MapBasic reference guide it seems the TableInfo function should serve this purpose but I can't seem to get it to work. In the MapBasic window whenever you use that function it says not unrecognised command. Any ideas? Cheers, Jon


title: Filter contour lines in QGIS

text: I have a shapefile of contour lines, with an `ALTITUDE` attribute with the height ASL in metres. The source data is in 0.5m intervals, but I'd like to filter it down so I get contour lines at (for eg.) 5m intervals. I tried doing a query on the layer using `"ALTITUDE" % 5 = 0`, which almost worked except I get two contour lines where I would expect one - I'm guessing this is because `%` is an integer operation so `245.5 % 5` and `245.0 % 5` both = 0? Is there a better way to filter for 5m (or some other interval, say 1m or 2.5m) data?


title: Use a Tiledmap service as basemap instead of an Arcgis.com map in AMD style coding

text: How will I define a arcgisTiled map service as the base map (instead of “topo”) as against using arcgis.com map services while defining a new map in AMD style?               map = new Map("map", {               basemap: "topo",               center: [-122.445, 37.752],               zoom: 14             });      I would like to replace "topo" with a map tiledservice from my server:               http://ourserver/arcgis/rest/services/tiledmapservice/MapServer      Any thoughts?


title: How to build a confusion matrix for image classification using SAGA or QGIS?

text: I did a supervised classification of a Landsat image subset using my training data polygons. I used SAGA GIS for that and everything worked out fine. Now I am stuck since I have no clue how to prepare the confusion matrix using my test polygons! I know that I need for each polygon the number of pixels of the different classes, i.e. a cross tabulation. How could one achieve that? And yes, I know that there is this semi-automatic image classification plugin for QGIS and so forth, but I would like to understand the process of generating the confusion matrix itself from scratch. thank you! Jens


title: How to increase maximum memory usage cap for running Arc 10.1?

text: Not sure if this is strictly GIS related or perhaps belongs on StackOverflow...it is regarding ESRI ArcGIS 10.1, so I'm sticking it here. I'm looking for options to increase processing capacity, and I've read unhelpful bits and pieces about increasing the cap on memory usage for ArcGIS. I'm working on ESRI ArcGIS 10.1 (all licenses) on a 32-bit Windows 7, 3.33GHz Intel DuoCore, 4GB RAM computer. I'll eventually be upgrading to a more powerful, 64-bit computer with 8GB Ram, and would like to know what to do in the meantime, or when the new one gets in. Thanks!


title: How to calculate the shared area between a raster and a shape file in ArcGIS 10?

text: I have a raster file with 5 classifications developed from habitat modelling algorithms. I'd like to know how much of 2 classes of the raster is covered by the network of protected areas (Shape file). Could any one advise me how to do so?


title: How to convert Landsat L1G raster to Geotiff?

text: Is it possible to convert a Landsat L1G raster file like L71001070_07020000728_B30.L1G to Geotiff using free software?


title: How to calculate area of lat/lon rectangle?

text: (Strange that this Q seems to be unanswered...) Given a rectangle like this, say (from Java, but that should be immaterial):               double minlat;     double maxlat;     double minlon;     double maxlon;      How to calculate the area of that rectangle? A math formula for square meters would be fine.


title: Select features by attribute if in Python list

text: I believe this will amount to largely a sql syntax question. Essentially, I am trying to complete a select by attribute in Python but based on the query of whether an attribute is present in a list. Such a query at its simplest should be something like this:               qry = " \"OBJECTID\" in oid_list"     arcpy.SelectLayersByAttribute_management(inft, "NEW_SELECTION", qry)      but that approach returns an invalid expression error. In the past, I've had to use more complicated sytax for this type of query, such as:               sqlQuery2 = "nid in (" + ','.join(["'"+x+"'" for x in delta_list]) +")"      but an adaptation of this snippet doesn't seem to work for me either, ie.:                "OBJECTID_1 in (" + ','.join(["'"+str(x)+"'" for x in oid_list]) +")"      So, I refer to you SQL experts: got any suggestions? What am I missing here?


title: Does linux suport "like '%c'" filter(cql_filter)

text: I can execute the following request: http://My.host:8080/geoserver/wfs?request=GetFeature&version=1.1.0&typeName=topp:states&cql_filter=STATE_NAME like 'California' But I can't execute the following one: http://My.host:8080/geoserver/wfs?request=GetFeature&version=1.1.0&typeName=topp:states&cql_filter=STATE_NAME like 'C%' Here goes the error message: XML Parsing Error: no element found Location: .http://My.host:8080/geoserver/wfs?request=GetFeature&version=1.1.0&typeName=topp:states&cql_filter=STATE_NAME%20like%20%27C%%27 Line Number 1, Column 1: I think something is wrong when I use % he encode the request and something happens... GeoServer:V2.1.4 Browser:FireFox SO: Ubunto


title: Append feature layer from one personal geodatabase to another using .Net

text: I currently have a client that still uses personal geodatabases (Access files). I currently have to open up ArcCatalog and manually append any FeatureLayers that are contained in the original personal geodatabase to the newly created one. I now need a non ArcObjects/ArcView way of doing this. The FileGeodatabase API does not support personal geodatabases (as far as I know). You cannot use a common OleDb connection to the personal geodatabase and transfer the feature layer tables correctly (must use the ArcObjects OleDB provider which requires ArcObjects). The above are the only things I have tried. Is there a way to transfer these personal geodatabase feature layers from one geodatabase to another, without a license for ArcView or ArcObjects? I am open for using the filegeodatabase api if there is a hackish way to do it, but I can't find one. I prefer a .Net solution, but would be able to create a wrapper for a C++ one. **Update** I am currently looking at the GDAL OGR drivers to see if it is possible to use them and accomplish what I am trying to do. Has anyone used them to access personal geodatabases?


title: How can I fix Print Composer errors due to huge file size?

text: I have put together a map using the British OS Strategi data set which is pretty densely detailed. It is surveyed at 1:250,000 and so contains roads down to track level, rivers, villages, railway lines - you name it. There's a lot of data. When I try and export an image version of the map (at about 1m x 1m), I always get an 'Exception: bad allocation' error. This happens when I try to export as a raster image or PDF. I have tried printing to a PDF generator (CutePDF) but the same thing happens. The print spooler shows the file size at 1.6Gb. I am running QGIS on a dual core 2ghz Windows 8 system with 4Gb of RAM. Is it just a case of the file being too big? Is there anything I can do? Or is there an alternative, free GIS package that would do a similar job?


title: How to choose raster interpolation method?

text: I work with a dataset of ~1.3 million residential buildings; each of them is described by a normally distributed variable ranging from 0 to 100. The dataset covers whole Switzerland so there are obviously areas with both very dense and very sparse point density. The limitation imposed by our data provider is that we cannot release datasets more detailed than 5 buildings. I was thinking about converting this point data into a raster surface. My main goal in this exercise would be to create the surface that minimizes the error that will arise when somebody tries to get the point value from this raster by extracting data in their's point locations of interest (as compared to extracting it from the original dataset of point locations). I'd like to know what kind of method should be best suited to do that. There are quite few to choose from in the Spatial Analyst toolbox and even more in Geostatistical one so I'd very much appreciate some help in where to start.


title: plpgsql | before update trigger: insert the new pseudo-record in another record and leave the old one as it was

text: I've got the following table:               CREATE TABLE alumbrado_puntos_luz     (       id numeric(8,0),       "GEOMETRY" geometry,       revision_actual numeric(10,0),       revision_expirada numeric(10,0),       CONSTRAINT alumbrado_puntos_luz_id_key UNIQUE (id)     )     WITH (       OIDS=TRUE     );      When updating the table i would like to update just the revision_expirada value with the next serial and insert the NEW pseudo-record's values in a new record on the same table. I mean, every time i have a change in the drawing (GEOMETRY) i will create a new record with the changes and leave the old one as it is. I don't know if a trigger is the best way of doing this or creating a rule would be a better idea. Regards,


title: Query DEM Elevation at Point in QGIS

text: I'd like to be able to pull up a DEM (or create a mosaic of them) and simply query the elevation at points I click. Is this possible using QGIS?


title: When computing flow direction, the D8 algorithm is used - are there major differences if others are used?

text: I am a student working with Hydrological models. Within these models we can route water across a DEM (within a delineated watershed or just within an un- delineated DEM) using many different flow routing algorithms: D8, MFD (sometimes FD8), Rho8, DEMON, D-Inf, and so on, with different variations within even the ones I listed possible (original MFD vs. MFD with a power term). But, when working within a delineated watershed, a flow routing algorithm (primarily the D8) has already been applied to get the flow directions grid - typically the first step in a delineation. Has there been any research looking at the difference between delineations performed on flow directions grids created using the other algorithms? I know since these other algorithms allow multiple flow directions, that it might be difficult to visualize this data (unless there is some compounding into one value) - and that that alone might make it impossible. It doesn't seem like the delineations would be any different - if point y is in the watershed of pour point x, water will flow from point y to pour point x, whichever way it may get there. The flow may be split into an infinite number of paths (infinite dispersion), but they will still all end up at point x, right? So even though D8 might not be perfect for predicting actual flow paths, it should produce the same delineation as any other, right? Thanks, Tom


title: Use GDAL command line to copy projections

text: I have two images, a georeferenced Tiff and a second Tiff without a reference. The second image is a copy of the first, but it lost its reference. Now, I want to copy the reference from the first Tiff to the second Tiff using GDAL using the command line. I had a look at several gdal utilities like gdalwarp, gdal_translate etc. but it seems that they only work on one image at a time and cannot "copy" a projection. Does anybody had a similar problem? Thx! Some background: I am using ImageJ for image analysis (segmentation) on 600 Tiff-images, all formerly geotiffs. However, during import ImageJ sets the xmin,ymin to (0,0) and looses all spatial information. Now, I try to use GDAL within an ImageJ macro to reload the spatial information so that the segmentation output also has a proper spatial reference


title: How to style points in QGIS based on multiple attributes

text: While trying to find a useful way having earthquakes displayed, I noticed that it isn't that easy to manage. I'd like to be able to recognize the "Magnitude", "Date", and "Depth" from the style of the point. Something like:   * size: based on magnitude   * depth: basic color (ex. <15 red, 15-30 yellow, >30 blue)    * date: intensity of color as shading from 100% color to white/... Is there a way to do so without defining dozens of rules?


title: How can I work with very large shapefiles (~1 GB) in QGIS?

text: I'm performing some basic spatial analyses on a large database of traffic injuries in California between 2003 - 2011. Because the dataset is large (nearly 1 GB of points), I'd like to first cut it down to a specific geographic region using a spatial query, but I find that QGIS consistently freezes and hangs if I try to use a spatial query or filter the layer. What are some ways I can work with this dataset? Are there more efficient data formats I can use?


title: Getting a list of zip codes with its adjacent zip codes

text: I need a list in the form of adjacency list that contains all zip codes in the United States that each has a list of adjacent zip codes. I tried searching for it on the internet. Many of them only provide a list, but none of them had a list showing the adjacent list of zip codes. How do I find something like that? The closest thing I found is a list of counties and their adjacent counties. There are several websites that I believe that they have the information but they are using it for applets like this website, http://www.freemaptools.com/find-zip-codes-inside-radius.htm If no data is available, could I use the counties list to generate the said list? EDIT: OP here, I was new to the whole GIS system. I don't know which software to use. I decided to try ArcGIS trial. I downloaded ArcGIS and that .lpk file. But I have no idea how to properly load .lpk in a python script and run PolygonNeighbors_analysis function. ArcPy Documentation and Google are not giving me a clear answer. Can anyone show me how to load .lpk package and use the function? Thanks EDIT2: I think I should have mentioned that I want that list in a simple .txt file. I think I could do this with ArcPY to collect all zip codes and their neighbors, then store it in the .txt file. My current obstacle to this is that I have no idea how to load .lpk file to ArcPY. EDIT3: I got it to work! By following what Barrett described, I got a raw .txt file of zip codes and their neighbors. I cleaned it up with my small program. It looks good. Thanks for help!


title: Labels with white (or other color) background

text: I tried to give my labels somekind of box and a backgroundcolor (or boxcolor) in QGIS, but i couldnt find the right option. My labels are not readable because my layers are too dark. I found offset but would prefere a small box in white. Can u help me?


title: ArcGIS: inside buffer with attributes?

text: How can I create buffers for polygons that a) lie only inside the polygon boundary, and b) know who their parent is? An example of where this is useful is First Nation Traditional Territories. There are overlaps all over the place and it’s nigh on impossible to sort out which is which, unless we leverage colour and/or pattern. The Buffer Wizard is the usual way we’ve been creating buffers. It’s relatively fast and convenient and allows one to save the result in the .mxd as well as in the file system (the former means you don’t have to think about where to save it, and just get on with the job). It’s major deficit though is that the buffers are dumb. They don’t know who they belong to or what their name is. In the example of FNTT, this means to draw Carcross Tagish FN differently from Teslin Tlingit Council and we have to create buffers twice, each time using selected features only. Suddenly no longer convenient. With attributed buffers can we use Definition Queries to display only certain FN’s and/or show each one with a separate symbol all from a single feature class. ![alt text](http://www.environmentyukon.ca/maps/media/uploads/preview/fntt.jpg)


title: Georeferencing with Quantum GIS and reprojecting to EPSG:900913

text: I'm trying to reproject a map to EPSG:900913 The image I'm reprojecting is the following http://i.imgur.com/HmfLW.png I'm using the cities as reference points obtaining the coordinates from Google Maps. I'm entering the coordinates for each of the reference point this way: ![Reference points](http://i.imgur.com/o8bJw.png) So X is longitude in degrees and Y is latitude in degrees. I'm using the following transformation settings and target srs: ![Trasformation settings](http://i.imgur.com/zAFnO.png) ![Target SRS](http://i.imgur.com/vDnDf.png) Given these settings the "Generate GDAL script" action gives me the following script:               gdal_translate -of GTiff -gcp 39.084 83.0534 12.9639 55.6032 -gcp 569.16 509.313 24.0161 49.838 -gcp 109.924 489.771 14.3701 50.1065 -gcp 420.153 334.656 20.9839 52.2413 -gcp 426.26 581.374 21.2256 48.7345 -gcp 549.618 142.901 23.8623 54.9019 miasta.bmp miasta_translated.bmp     gdalwarp -t_srs epsg:900913 -r near -order 1 -co COMPRESS=NONE  miasta_translated.bmp miasta_modified.tif      It is strange it is not adding the -t_srs parameter to the gdalwarp command, so I've added it myself. I expect the miasta_modified.tif to contain the Google Maps projection of the original map but instead it contains the following which doesn't look like the Google Maps projection the most noticeable difference being that it is compressed in the vertical direction: http://i.imgur.com/9DMyT.png (I've converted it to PNG for the image hosting to accept the file) What am I doing wrong? Thanks!


title: How to map our data with Google Map data

text: I want to add all the bus route data and create a route finding app for my country since Google map has no data about bus routes in my country but with JavaScript API i could get the shortest path between two cities. e.g if i am going to D from A. it says from A to B then turn left go to C and then go to D likewise. My idea is to get those points and i have my database with bus stops and route numbers, map the data and create my app. but the problem is the points which gives by Google map and the bus stops are not always the same. so i want to know whether there is a way in Google maps to add my bus stop data and bus route information along a particular road then i think it can easily be searched. Thanks..


title: Why does GPS positioning require four satellites?

text: I have a question on the GPS positioning algorithm. In all books I've read for 3D positioning we need four satellites, and I don't understand why. We need to calculate three variables: x, y, z. We know when satellite send the signal to earth and when we receive it we can measure the time the signal travel to earth by checking the shift in PRN generator. For what purpose do we need four satellite?


title: Help with customizing GeoExplorer

text: I'm just beginning to play with BoundlessGeo's GeoExplorer interface. I'm trying to figure out if I can customize it for use in-house. We basically need an intra-net type web mapping app to display various data layers (client project info, location, etc.). We need to do some custom queries, and there;s interest in a tabbed identify display. Anyone know if it can be customized in such ways? I'm finding the documentation kind of lacking, and this seems like a logical place to turn for answers. Any help appreciated! -monkey5672


title: Performance Tuning on CartoDB

text: What kind of tuning can be done on your own cartodb instance to speed things up (outside of sql query optimization and indices)? Is there any documentation on this? If the records do not change, could that make a difference?


My test scenario is to confirm the image below, I select the points on the line in order of the direction the line was drawn. ![Order points](http://i.stack.imgur.com/OZcnf.png) I used the sql:                   Select       testepoint."id"      FROM         testepoint     inner JOIN         teste     ON       ST_DWithin(teste.geom, testepoint.geom, 1)      WHERE         teste."id"= 1      ORDER BY       ST_Line_Locate_Point(teste.geom, testepoint.geom)       I placed the order by to sort the way I want, as the line was drawn                   ORDER BY       ST_Line_Locate_Point(teste.geom, testepoint.geom)       And the error occurs:                ERROR: line_locate_point: 1st arg is not a line      This was because my geometry is a MultiLineString. How can I do this sort of which is a MultiLineString?


title: QGIS WCS client with clip parameter

text: Can be defined bbox parameter in a WCS service? I require trimming the layer coordinate in a box, I'm trying to prove the following WCS, but shows the entire layer globally: http://geobrain.laits.gmu.edu/cgi-bin/gbwcs- dem?bbox=-76.70429,7.26988,-74.59818,9.56054&crs=epsg:4326 regards, Samuel Mesa.


title: How to use Fields in the Field Calculator?

text: I'm using the Field Calculator in ArcMap and I need to create a unique ID for every storm drain in my county. An ID Should look something like this: 16-I-003 The first number is the municipal number which is in the column/field titled "Munic" The letter is using the letter in the column/field titled "Point" The last number is simply just 1 to however many drains there are in a municipality. So far I have:               rec=0     def autoIncrement()     pStart=1     pInterval=1      if(rec==0):       rec=pStart      else:       rec=rec+pInterval     return "16-I-" '{0:03}'.format(rec)      So you can see that I have manually been typing in the municipal number, the letter, and the hyphens. But I would like to pull the first number and letter from the fields: Munic and Point so I don't have to manually type them in each time it changes. I was hoping I could just put !Munic! and !Point! in my return statement before the .format(rec) or somehow work it in there. I'm a beginner when it comes to python and ArcMap, so please dumb things down a little. The simpler the better!


title: Unique ID in ArcMap 10.2

text: In an attribute table, how should I go about adding a Unique ID into a field? For example, I have a Water Meter points layer with an attribute table but I need a field that list them numerically with WM attached to the numbers. I would like to have a field displaying WM1001, WM1002, etc. Is there any way to accomplish this? Thanks.


title: Retrieving the data from SQL database is taking much more time than expected

text: I am creating a block wise vector layer and related data stored in a sql database while retrieving the data from the SQL database through querying and would represent the result on a map. This is taking a lot of time and there is no problem with the data structure because it performs all statistical operations very quickly. What is taking time is the processing required to represent the result on the map. Database is sql server, software used Map Suite GIS Editor, no wms, Does anybody have suggestions as to what the problem might be and what is the possible solution?


title: How to assign values to a polygon layer based on nearest point in a point layer?

text: I have a point file depicting groundstations for rainfall measurment. I also have a polygon file representing districts of a state. I want to assign rainfall values in a particular district based on the station located there. I understand this can be achieved through a spatial join. However, I also want to assign these values to near by districts based on a value taken from the nearest station in the surrounding districts. Please note that   * a single district can theoretically have multiple stations.    * I want to assign the exact value, and not an interpolated value. I am working on ArcGIS 9.2. So any solution based on Arc9.2 or QGIS 1.6 is highly appreciated.


How to get pixel size of canvas in QGIS? It is possibly by python?![example](http://i.stack.imgur.com/ERjLJ.jpg) I mean this read rectangle area.


title: Isolate And Save Raster Features From Existing Multivalue Raster Map

text: I have Cyprus Corine map with 44 different types of land use. I want to run a raster calculator query or any other way to isolate the areas with codes 10 to 22 and 26 to 34, all other values 1 to 44 will be excluded. When I open the Corine Raster layer and I choose the Identify Features Tool I always see the Band 1 under the Feature Category, the number under the Value Category it change every time I chose a different place, when I try to run the query it choose all the values. I found the map from: https://dl.dropbox.com/u/13604252/dem_gr/data.html at the end of the page on the right choose: CLC.r.data.UTM36. How can isolate the values I want to save as a new Raster map?


title: Should GDAL be set to produce GeoTIFF files with compression? Which algorithm should be used?

text: I have a folder of GIS data that consists mainly of GeoTIFF files. The whole set weighs in at about `1.2 GB`. I noticed that if I pack the contents into a tarball, it smashes down to about `82 MB`. I would like to check the set into a revision control system sot it can be worked on by other people and it looks like there is some space that can be squeezed out. The GDAL GeoTIFF driver page lists plenty of options that may be used to create compressed GeoTIFF files. There are also plenty of options that affect the way each algorithm works. The help page does a good job at describing the options but doesn't elaborate on how to select an algorithm or the tradeoffs that are associated with the varying level of compression. This leads to the following questions:   * The pros of using compression are a dramatic savings in space. What are the cons? Is information lost when the image is compressed?   * How should one go about choosing an algorithm and compression level. Do some types of images lend themselves to a certain algorithm?


title: ogrinfo SQL GROUP BY Clause

text: Does `ogrinfo` support a _GROUP BY_ clauses? I could not find any good documentation of all the SQL syntax features it supports. I tried the following:               ogrinfo . -geom=no -sql 'SELECT lsad FROM tl_2011_us_county GROUP BY lsad'      And got the following error:               ERROR 1: SQL: Failed to parse SELECT statement, extra input at BY token.     INFO: Open of `.'           using driver `ESRI Shapefile' successful.      What I really want to do:               SELECT lsad, count(*)     FROM tl_2011_us_county     GROUP BY lsad      What I can do:               SELECT DISTINCT lsad     FROM tl_2011_us_county      And then:               SELECT count(*)     FROM tl_2011_us_county     WHERE lsad = 'foo'      ... but that is annoying.


title: OGR and DXF files: Is it possible to handle layer names? (only "entities")

text: I'm using Python and OGR to extract the layers of a DXF file and convert them into SHP. I began with:               import ogr, os, sys     driver = ogr.GetDriverByName('DXF')     datasource = driver.Open('test1.dxf', 0)     numberLayers = datasource.GetLayerCount()     for i in range(0, numberLayers):         layer = datasource.GetLayerByIndex(i)         layerName = layer.GetName()         numberFeatures = layer.GetFeatureCount()         print 'Layer=%s|Features=%s' % (layerName, numberFeatures)      but I get one only layer: **entities** (as expected: http://www.gdal.org/ogr/drv_dxf.html) Obviously, it's the same information as if I execute: `ogrinfo -so test1.dxf entities` My DXF file contains several CAD layers: "ROADS", "ELEVATION", ... Is it possible to handle these layer names with OGR? Thank you very much.


title: Is there any site which provide classified image download for soil,vegetation or water resources?

text: (its particularly for indian states,i need.....,but kindly quote if any other info you have) the raster data may also work for me in above aspects.


title: What is the best practice for visualisation of a spatio-temporal data set?

text: I'm investigating visualisation of a spatio-temporal data-set. What examples are available, online or offline, of high-quality spatio-temporal data visualisation?


title: How to access vector coordinates in GRASS GIS from python?

text: Is there a way to iterate over lines and their vertices using python? I'm new to GRASS GIS, and it looks like there are ways to read raster (for numpy), but I don't see a direct way to access coordinates and topology for vectors. I know that I can access binary data but it does not seem right. I can convert vector feature to db for later use with vector_db_select but it is an extra operation for conversion and I have no data what was what (how can I walk a graph?). Also when I try to run `v.to.db` like               v.to.db map=merged_single@PERMANENT layer=1 option=coor units=meters columns=cat      I'm getting               This option requires at least two columns     Finished with error      For instance in ArcGIS I can easily iterate over elements using SearchCursor. I wonder if there is anything similar in GRASS GIS? **Update 2012-06-24** After reading GRASS GIS manual, it looks like one has to use C library and functions prefixed with `Vect_`. Another possible choice probably is to use OGR as it somewhat supports GRASS GIS, but its python bindings are not perfect. It all seems a bit awkward as QGIS so nicely works with all formats from GUI so the inability to access all that from python (which is the official scripting language) seems like a huge shortcoming unless I missed something. Perhaps it will be better in the future.


title: Why i cannot see my globe service?

text: I want to explore it in ArcMap or in ArcExplorer but both cannot "see" it. For the other resources like map services everything is workikng well. I can add that the service is started and in a preview view eveything seems to be ok.


title: Is it possible to have pages in mapbook link to different data frames on 1 layout page in ArcMap?

text: I am trying to create a map series of approximately 80 sheets, but I need to have 2 sheets per layout page. Basically, I need to know if there is a way to have the pages in the map series go to different data frames inside my map? Any help would be greatly appreciated. thank you!


title: Cannot upgrade User Schema Geodatabase when Master Schema Geodatabase is at previous release?

text: I am getting this error when trying to upgrade ArcSDE from 10.0 to 10.1: > Cannot upgrade a user schema geodatabase when the master schema geodatabase > is at a previous release. > > The master schema geodatabase must be upgraded first. > > ERROR: Geodatabase schema object install not completed (-522). However, I am updating the ArcSDE master geodatabase only, so why am I getting this error?


title: UAV Image collection Procedure

text: Do we need to collect ground control points before flying drones or do the drones collect control points during their flight ?


title: What does it mean "GML is an XML grammar"?

text: I'm trying to get my head around the theoretical side of XML and GML. The OGC web page for GML states that: > The Geography Markup Language (GML) is an XML grammar for expressing > geographical features. But what does it mean by a " _grammar_ "? I can't see that phrase (as a noun) used anywhere else. Is it a unique GML meaning? I see the terms "language", "schema", "format" and "standard". Are they the same thing? So can I say that GML is an XML schema? Or would that mean something different? If so, is there a single XSD file that defines the entire GML specification?


title: ArcGIS 10 Addin - Code Doesn't Execute on Another Computer

text: I've created an ArcGIS 10 Addin that consists of a toolbar containing one button and one tool. The button opens a form, and the tool is disabled until the user enters some information on the form. This addin works great on my development computer and on other computers in my office. But when I deliver the addin to the client, toolbar is there with the button and tool, but when he clicks the button, nothing happens. Not knowing if this had something to do with the form, I added a third button to the toolbar which displays a simple message when the button in clicked. For the client, a click on this button also results in nothing happening. The addin was created using vb.net and Visual Studio 2010. I've checked that he has the same ArcGIS 10 service pack installed that I have, and that he has the correct .net framework installed. What else could be missing on his computer? This is driving me a little nuts, so I would appreciate any and all suggestions. Thanks.


title: Is QGIS suitable for noise modelling?

text: Hi I've never used any type of GIS software before so I'm trying to find out if it is suitable for what I'm trying to do: I've written a c++ program which estimates ground level noise from a distant source as a function partially of distance. There are three things I'm looking to be able to do from QGIS,   * select a position in 3d space,    * find the distance to the terrain below this point, and    * plot the resultant noise onto the terrain. Is this sort of functionality achievable without significant work? If so do you have any suggestions for what should be looked up in the docs? Many Thanks


title: Does ArcEditor use ArcEngine?

text: If I use ESRI ArcEditor / ArcView / ArcInfo is it using ESRI ArcEngine under the hood?? I think the answer is no, but please enlighten me.


title: How to combine information from raster and vector layers?

text: in QGIS I have 1 raster-layer (elevation) and 1 vector-layer (administrative units). Now I would like to be able to query one of the administrative units and see the mean and min/max of the elevation within its area. Which QGIS tool can do that? THX, stn


title: Copy Features giving ERROR 000210: Cannot create output in ArcPy?

text: I have an issue where I'm trying to save separate feature classes for each route in an already solved closest facility problem, and one seemingly arbitrary route won't save. I tried to look at the original feature class that had all the routes, and I couldn't find anything particularly notable about this one route. Furthermore, just to get an idea of whether there was more than one feature that had a problem, I tried a try-except in arcpy, but every feature class after the one seemingly arbitrary broken one saved just fine. Moreover, when I used Model Builder to create the feature classes, everything worked just fine, without any hitches. But I have to use Python, and it's just not working. Can anyone give any pointers? Here's the code, although I'm not convinced it's even going to be helpful. (Note: Obviously, the code looks kind of rough, because a good bit of it was exported from Model Builder, and I just reworked what wouldn't export. Also, I'm well aware that turning the features into layers, doing a select, and copying the layers back to features is inelegant in comparison to simply doing a select_analysis, but that's just the last version of the code. Neither way worked.):               # Import arcpy module     import arcpy          # Check out any necessary licenses     arcpy.CheckOutExtension("Network")          # Script arguments     Building_network_1_nodes_1_1 = r"C:\\Users\\Me\\Desktop\\TestFolder\\TestGeo.gdb\\Test01\\Building_network_1_nodes_1_1"           # Local variables:     Test01_ND = r"C:\\Users\\Me\\Desktop\\TestFolder\\TestGeo.gdb\\Test01\\Test01_ND"     Test01 = r"C:\\Users\\Me\\Desktop\\TestFolder\\TestGeo.gdb\\Test01"     Test01__3_ = r"C:\\Users\\Me\\Desktop\\TestFolder\\TestGeo.gdb\\Test01"          # Process: Make Closest Facility Layer     outNALayer = arcpy.MakeClosestFacilityLayer_na(Test01_ND, "Closest_Facility_Model01", "Length", "TRAVEL_TO", "", "2", "", "ALLOW_UTURNS", "", "NO_HIERARCHY", "", "TRUE_LINES_WITH_MEASURES", "", "NOT_USED")          # Process: Select     arcpy.Select_analysis(Building_network_1_nodes_1_1, r"C:\\Users\\Me\\Desktop\\TestFolder\\TestGeo.gdb\\Test01\\Test01_FModel", "node_type = 3")          # Process: Add Locations     arcpy.AddLocations_na("Closest_Facility_Model01", "Facilities", r"C:\\Users\\Me\\Desktop\\TestFolder\\TestGeo.gdb\\Test01\\Test01_FModel", "", "5000 Meters", "", "Building_network_1mz_1_1 SHAPE;Building_network_1_nodes_1_1 SHAPE;Test01_ND_Junctions NONE", "MATCH_TO_CLOSEST", "APPEND", "NO_SNAP", "5 Meters", "INCLUDE", "Building_network_1mz_1_1 #;Building_network_1_nodes_1_1 #;Test01_ND_Junctions #")          # Process: Select (2)     arcpy.Select_analysis(Building_network_1_nodes_1_1, r"C:\\Users\\Me\\Desktop\\TestFolder\\TestGeo.gdb\\Test01\\Test01_IModel", "node_type = 1")          # Process: Add Locations (2)     arcpy.AddLocations_na("Closest_Facility_Model01", "Incidents", r"C:\\Users\\Me\\Desktop\\TestFolder\\TestGeo.gdb\\Test01\\Test01_IModel", "", "5000 Meters", "", "Building_network_1mz_1_1 SHAPE;Building_network_1_nodes_1_1 SHAPE;Test01_ND_Junctions NONE", "MATCH_TO_CLOSEST", "APPEND", "NO_SNAP", "5 Meters", "INCLUDE", "Building_network_1mz_1_1 #;Building_network_1_nodes_1_1 #;Test01_ND_Junctions #")          # Process: Solve     Solved = arcpy.na.Solve(outNALayer, "SKIP", "TERMINATE", "")     arcpy.CopyFeatures_management("Routes", r"C:\\Users\\Me\\Desktop\\TestFolder\\TestGeo.gdb\\Test01\\Test01_Routes")     arcpy.CopyFeatures_management("Facilities", r"C:\\Users\\Me\\Desktop\\TestFolder\\TestGeo.gdb\\Test01\\Test01_Facil")     arcpy.CopyFeatures_management("Incidents", r"C:\\Users\\Me\\Desktop\\TestFolder\\TestGeo.gdb\\Test01\\Test01_Incid")          arcpy.env.overwriteOutput = True          # Script arguments     Routes = r"C:\\Users\\Me\\Desktop\\TestFolder\\TestGeo.gdb\\Test01\\Test01_Routes"     Facilities = r"C:\\Users\\Me\\Desktop\\TestFolder\\TestGeo.gdb\\Test01\\Test01_FModel"     Incidents = r"C:\\Users\\Me\\Desktop\\TestFolder\\TestGeo.gdb\\Test01\\Test01_IModel"          # Local variables:     Test01 = "C:\\Users\\Me\\Desktop\\TestFolder\\TestGeo.gdb\\Test01"     NATest01_071414_gdb = "C:\\Users\\Me\\Desktop\\TestFolder\\TestGeo.gdb"     Test01_A = "C:\\Users\\Me\\Desktop\\TestFolder\\TestGeo.gdb\\Test01_A"          u = [row[0] for row in arcpy.da.SearchCursor(Routes, "IncidentID")][-1]          for v in range(1, u + 1):         arcpy.MakeFeatureLayer_management(Routes, "Route" + str(v))         arcpy.SelectLayerByAttribute_management("Route" + str(v), "NEW_SELECTION", """ "IncidentID" = """ + str(v))         arcpy.CopyFeatures_management(Routes, r"C:\\Users\\Me\\Desktop\\TestFolder\\TestGeo.gdb\\Test01_A\\Route_" + str(v))         print "Route " + str(v) + " exported."      The whole rest of the script would work until the last part. Here's an abridged version of the print log:               Route 1 exported.     Route 2 exported.     ...     ...     ...     Route 37 exported.     Route 38 exported.          Traceback (most recent call last):       File "C:/Users/Me/Desktop/TestFolder/Test_01", line 53, in <module>     arcpy.CopyFeatures_management(Routes, r"C:\\Users\\Me\\Desktop\\TestFolder\\TestGeo.gdb\\Test01_A\\Route_" + str(v))       File "C:\Program Files (x86)\ArcGIS\Desktop10.2\arcpy\arcpy\management.py", line 2429, in CopyFeatures     raise e     ExecuteError: ERROR 000210: Cannot create output C:\\Users\\Me\\Desktop\\TestFolder\\TestGeo.gdb\\Test01_A\Route_39     Failed to execute (CopyFeatures).      I'm at a loss. I have no idea why I'd get that error for that particular feature class. It seems so, like I said, arbitrary. Like I said before, too, it works just fine in Model Builder, but I get this hiccup with arcpy. If anyone has ideas, I'd like to hear them. I feel like I've done all the obvious things, so far.


title: How to check if a point is near or inside polygon?

text: I am trying to find a simple library which can do point in polygon. Earlier i did hit Java Spatial Index library ...could not figure out. Tried openmap..looked simple didn't work. I am totally new to GIS..my application has polygons. i need to search if my point isin a polygon or new polygon(s) in open map i wrote as follows.               package abc.poi;          import com.bbn.openmap.geo.Geo;     import com.bbn.openmap.geo.OMGeo;          public class PoiTest {              public static void main(String[] arg) {             Geo geo1 = new Geo(1,1);             Geo geo2 = new Geo(3,2);             Geo geo3 = new Geo(4,4);             Geo geo4 = new Geo(3,4);             Geo geo5 = new Geo(4,3);             Geo geo6 = new Geo(1,1);                  Geo listGeo[] = {geo1,geo2,geo3,geo4,geo5,geo6};                  OMGeo.Polygon poly = new OMGeo.Polygon(listGeo);                  System.out.println(poly.isPointInside(new Geo(2,2)));                   }     }      It always returns false. What i am missing ?


title: Selecting and extracting data from TIFF files

text: I have some tiff files (2 for example) and a shapefile. The real format of the tiff files is geoTiff. I want to get the points inside the shape boundary. Each point should have 4 parameters: geom(coordinates), red, green, infrared (we have tiff files with such info). At the end of the operation I want to get a vectorized shape with points. Could you please give me some advice on how can I do that programmatically. Are there any libs/tools? It is all in the context of geoserver/postgresql/geomajas framework. Upd: now I have more information. My tiffs will be one of the Landsat combination: http://gif.berkeley.edu/documents/Landsat%20Band%20Information.pdf The main question now is: to get RGB channels info in exact point, should I parse image points on RGB,or I can get this info from geoTiffMetadata (there is http://download.deegree.org/deegree3/nightly/core/javadoc/org/deegree/coverage/raster/io/imageio/GeoTiffIIOMetadataAdapter.html ). I didn`t find really useful info while debugging the code and calling the methods of GeoTiffIIOMetadataAdapter.


title: How to use MasterMap address layer with QGIS?

text: I have a copy of QGIS. I have a MasterMap address layer for my area. How do I get the map onto my computer and open it please?


title: How can I remove old partially installed postgis 2.0 from Mac?

text: I am trying to install postgis in mac using brew. However it shows that postgis-2.0.0 is already installed. Actually I tried to install it from source. Up to configuration it was fine but when I did make, it showed some errors and didn't continue. Further I installed postgresql from brew so I think it would be better to install postgis from brew as well. How can I remove this old partially installed postgis of version 2.0.0


title: Is it possible to manually correct a raster DEM (OpenSource)?

text: Is it possible, to correct a raster DEM by changeing the raster-cell-values manually, like by "painting"? My DEM (ASTER) have one small error where the elevation is about 300m below the actual elevation. It should be a hill, and not a hole... SRTM 4.1 is good, but has other errors (the hill is very steep and small, so SRTM "ignores" it (<90m). ![enter image description here](http://i.stack.imgur.com/qddkR.png)


title: Calculate Moving Time

text: I am wondering how web apps like Garmin connect or Strava calculate a "moving time" from a gps track. For example I might record a GPS track on my Garmin GPS that has an elapsed time of 5:34:26 but they calculate the moving time as 2:33:39. I can't seem to find an discussion or examples on how this is done, but maybe I'm searching for the wrong things.


title: How can I get a popup window into the ArcGIS online widget_Legend sample?

text: I am trying to use the sample code from ArcGIS javascript API called widget- Legend. I want to add the popup widget to this sample so it will include both a legend and a popup. The legend works fine with my hosted feature layer, but the popup window does not. I tried to paste the code from the popup window sample into the code I am using and the screen just turns completely white! It only crashes when I put "infoWindow = popup" inside the map constructor. I am a student and new to code. Does anyone know how to solve this issue?...Here is the code starting from ...               <style>     html, body {       height: 97%;       width: 98%;       margin: 1%;     }               .esriPopup.myTheme .titlePane,     .dj_ie7 .esriPopup.myTheme .titlePane .title {       background-color: #E68A2E;       color: #333333;       font-weight: bold;     }          .esriPopup.myTheme .titlePane {         border-bottom: 1px solid #121310;     }     .esriPopup.myTheme a {         color: #d6e68a;     }     .esriPopup.myTheme .titleButton,     .esriPopup.myTheme .pointer,     .esriPopup.myTheme .outerPointer,     .esriPopup.myTheme .esriViewPopup .gallery .mediaHandle,     .esriPopup.myTheme .esriViewPopup .gallery .mediaIcon {         background-image: url(./images/popup.png);     }          .esriPopup.myTheme .contentPane,     .esriPopup.myTheme .actionsPane {       border-color: 1px solid #121310;       background-color: #424242;       color:#ffffff;      }     #rightPane {       width: 20%;     }          #legendPane {       border: solid #97DCF2 1px;     }          .esriPopup.myTheme .titlePane,     .dj_ie7 .esriPopup.myTheme .titlePane .title {       background-color: #E68A2E;       color: #333333;       font-weight: bold;     }      var map; require([ "esri/map", "esri/layers/FeatureLayer", "esri/dijit/Legend", "dojo/_base/array", "dojo/parser", "esri/dijit/Popup", "esri/dijit/PopupTemplate", "dojo/dom-class", "dojo/dom-construct", "dojo/on", "dojox/charting/Chart", "dojox/charting/themes/Dollar", "dijit/layout/BorderContainer", "dijit/layout/ContentPane", "dijit/layout/AccordionContainer", "dojo/domReady!" ], function( Map, FeatureLayer, Legend, arrayUtils, parser, Popup, PopupTemplate, domClass, domConstruct, on, Chart, theme ) { parser.parse();                 var popup = Popup({           titleInBody: false       },domConstruct.create("div"));            map = new Map("map", {         basemap:"topo",         center: [-86.60, 35.63],         zoom: 8,         infoWindow = popup                 });           //add the legend       map.on("layers-add-result", function (evt) {         var layerInfo = arrayUtils.map(evt.layers, function (layer, index) {           return {layer:layer.layer, title:'Tornado Intensity (F Scale)'};         });         if (layerInfo.length > 0) {           var legendDijit = new Legend({             map: map,             layerInfos: layerInfo           }, "legendDiv");           legendDijit.startup();         }       });            domClass.add(map.infoWindow.domNode, "myTheme");           var template = new PopupTemplate({          title: "Tennessee Tornados",          description: "{FATALITIES} people died in this tornado",          fieldInfos: [{ //define field infos so we can specify an alias            fieldName: "DATE_",            label: "Date"          },{            fieldName: "FSCALE",            label: "Intensity"          },{            fieldName: "FATALITIES",            label: "Fatalities"          }],            mediaInfos:[{ //define the bar chart              caption: "Fatalities to Intensity",              type:"barchart",              value:{                 theme: "Dollar",                 fields:["FATALITIES","FSCALE"]              }            }]         });            var featureLayer = new FeatureLayer("http://msgis-664-14.spatial.redlands.edu:6080/arcgis/rest/services/Micah/664_Assignment_3/MapServer/0", {         mode: FeatureLayer.MODE_ONDEMAND,         outFields:["*"]       });           map.addLayers([featureLayer]);     });      html, body { margin: 0; } </style> This pane could contain tools or additional content


title: Land Parcel and JSON Data with onclick event

text: I currently am working on a school project and have been introduced to javascript with little knowledge. What the project entails is a webapp of land parcels that contain information about each one, I currently have that data in JSON format converted from an MDB. Both parcels and JSON data have a unique key that connects them. Anyways, the end goal is to be able to click on the parcel and query the JSON data to display any information connected by the unique keys (this can be plural). This can be in an info window or panel. I have gotten to a point where when I load the map, it will load the data from the JSON and I can choose which row to show in Firefox console. I haven't been able to figure the best route to include an onclick event and to be honest have not found a good resource for displaying data from the JSON into a display. Below is a snippet of the code that loads and displays data in the Firefox console.               function initialize() {          var latlng = new google.maps.LatLng(43.044131, -79.283196);          var myOptions = {              zoom: 12,              center: latlng,              mapTypeId: google.maps.MapTypeId.HYBRID          };          var map = new google.maps.Map(document.getElementById("map_canvas"), myOptions);         var layer;              var tableid = '11WmrR1ElZMrqGaavOHV9VGMvTcBsW4GZtFycLpmy';              layer = new google.maps.FusionTablesLayer({                    query: {                          select: 'geometry',                          from: tableid                    },              });              layer.setMap(map);                   var xhr = new XMLHttpRequest();         var MVRows;         xhr.open('GET', 'MV.json', true);         xhr.onload = function() {             MVRows = JSON.parse(this.responseText).dataroot;             console.log(MVRows._x0031_993to1996[0]);         };         xhr.send();          }       I hope I was concise enough, thanks all! Also, I don't want anyone to do my project for me, just some helpful hints or advice. -Dave


title: Desire Lines - 080040 : Error while executing Desire Lines, anyone know how to fix this?

text: Using ArcGIS with Business Analyst: I have a customer data file that contains around 723K customer points for 79 stores that I am trying to calculate the drive times for. I have put my customer and store files into the Business Analyst Store and Customer setup. Whenever I try to run desire lines I keep getting this error: "080040 : Error while executing Desire Lines". Anyone know of a solution? I really need to be able to run all the stores and customers at once as doing one store at a time would be cumbersome and too time consuming.


title: OGR Info - Check for spatial index?

text: If we want to get information about a vector layer (shapefile for example) we could use ogrinfo:               ogrinfo -so myshapefile.shp mylayer      This returns a lot of pertinent information about the shapefile. Is it possible to return whether or not a spatial index exists? If so, can we access info about the spatial index? Note: I realize it is trivial to open the folder containing the shapefile and looking for a .sbn or .sbx and so the question is specific to ogrinfo.


title: Generating regular sample points in a polygon

text: Is there a way I can generate regularly distributed sample points within a polygon using ArcGIS 10? I originally used Hawth's Tools to do this in 9.3, however I don't think this is possible in 10. Thanks!


title: Is it possible to render tile-based web mapping projects in any other projection than web mercator?

text: Web-mapping technology has to a large degree adopted the web-mercator projection as the standard projection for displaying geographic objects with interactive map-clients. This is due to many factors including the mathematical simplicity of the cylindrical projection which enables fast tile creation, the relatively good preservation of shapes and angles which attracted seafarers to the projection presented by mercator in 1569, and of course it´s adoption by google, bing, ArcGIS Online which further strengthens it´s position as a standard web-mapping projection as in order to create slippy maps, one has to conform to the base-map projection as given. I am however creating a geoportal for a region of the world which has absolutely no interest in displaying other regions of the globe. The map will also never be used for sea navigation as it is completely landlocked. Having reliable map resolutions which do not change with latitude is also significant due the necessity of presenting formal map objects(e.g Planning documents) with a given scale which is a massive disadvantage of the web mercator projection. My question is: Is it at all possible to render tile-based web mapping projects in any other projection than web mercator and what are the implications regarding performance and interoperability with other WMS services if this possible? This would be advantageous for regional web mapping projects in order to offer added functionality which is often missing e.g reliable distance measurement over relatively large areas and also support of human understandable resolutions and therefore scales instead of zoom-levels giving users the ability to render maps at 1:1000, 1:10,000 etc.


title: Extract Longitude and Latitude from Geotiff (ASTER 30m DEM)

text: my problem is apparently quite simple, but I cannot really find an easy answer: in QGIS I've created a raster layer by merging different GeoTiff images from the ASTER 30m DEM dataset. The merging is ok: by moving the cursor I can see that the coordinates make sense and so the elevation. What I'd like now is to extract the longitude and latitude informations from each pixel (generating a matrix the size of the layer with Lat values and one with Lon values), but I really don't know how to do it. I also tried to import the image in MATLAB using the geotiffread function, but I cannot extract such information... Thanks to all!


title: Raster slightly displaced after exporting from PostGIS database

text: With the following queries I've exported a raster (GTiff) from my postgis database: 1.Creating a large object:               WITH buffer AS (         SELECT             id,             ST_Buffer(geom,50) AS geom         FROM polygon_table),          clip AS (         SELECT              partnumber,             ST_Union(ST_Clip(a.rast, 1,b.geom, TRUE)) AS rast     FROM aspect a             JOIN buffer b             ON ST_Intersects(a.rast, b.geom)                 GROUP BY partnumber)          SELECT          oid,          lowrite(lo_open(oid, 131072), png) As num_bytes         FROM (             VALUES (lo_create(0),                    ST_AsGDALRaster(((SELECT ST_Union(rast) FROM clip)),'GTiff')                )) As v(oid,png)      2.Exporting it to a GTiff file (e.g. with oid=176148)               SELECT lo_export(176148, '/tmp/demo_rast.tif')      (3. Deleting the large object: `SELECT lo_unlink(176148)`) The aspect raster data in the postgis database has the following information:                rid | scalex | scaley |  srid     -----+--------+--------+-------       1  |   25   |  -25   | 31468      The aspect raster resolution is 25m. Tiling is 100x100. When load the raster into QGIS it is slightly displaced. The colored aspect raster is the original. The grayish one is the exported one. ![enter image description here](http://i.stack.imgur.com/M2d0t.png) The amount of the displacement is about 2.65m to the west in x direction and about 1.35m to the north in y direction. Where does this displacement come from? **UPDATE** The problem appears exclusively in QGIS. When I load both rasters into GRASS GIS there is no displacement. I didn't used the option `Override projection (use location's projection)` when I import the rasters into GRASS GIS. As you can see here: ![enter image description here](http://i.stack.imgur.com/C3NFr.png) So is there a general problem (bug) which appears here in QGIS???


title: PostGIS geography query returns a string value

text: I have a strange issue. The lonlat column on my app works well on the development server –– its output is in the form of POINT(X Y). But when I move the data to the production server, the output is strange! ActionView::Template::Error (undefined method `lon' for "0101000020E6100000541B9C887E7A52C02920ED7F80614440":String): The lonlat value, which is encoded with SRID: 4326, is being read as a string. I am almost certain that there was a corruption in the data during migrating it from development to production because this was not a problem before the migration. Does anyone know what about the database schema or column may cause this issue?


title: Convert OSM File to get addresses using Osm2pgsql

text: I am working to covert OSM File to get all building's informations in the city. I am using osm2pgsql package to transform the data. Now I am confused as I have 5 tables :                  - planet_osm_point        - planet_osm_line        - planet_osm_roads        - planet_osm_polygon      https://github.com/openstreetmap/osm2pgsql I am confused now, how can I get buildings' informations from these tables. What do every table represent ?


title: How can I optimize ArcMap Editor's streaming mode for use with a tablet such as the Wacom Intuos?

text: I have the Wacom Intuos 4 Large tablet which I use occasionally outside of ArcGIS. I now need to digitize a large volume of polygons from aerial imagery and am finding the tablet's ArcGIS functionality cumbersome. The pen tip works fine in a regular mode where I tap it to create a vertex (albeit far too easy to mistakenly engage a double click) but **I would like to be able to delineate features as if I were doing it on paper - a fluid hand motion**. I turned on streaming and set it to 5 map units which is OK but when I reach a point where I need to lift and reposition my hand on the tablet or pan the page while still digitizing, streaming does not stop when I lift the pen from the tablet. Ideally, I would like the streaming to behave exactly as it does in Google Earth - streaming pauses when the tip is lifted and resumes when a pressure is again applied to the tip. GE also generalizes sketches on completion not while sketching as ArcMap does - this makes the finished product simple but precise. Alternatively, I used ArcMap's Tablet toolbar which streams beautifully to produce "ink" sketches. However, when I lift the pen from the surface, line is interrupted and while this is fine for polylines, the line interruptions create artifacts when creating polygons. **How to I set up ArcMap Editor streaming so that streaming pauses when I lift the pen tip from the tablet's surface?** OR **How can I draw an ink sketch and pause drawing when the pen is lifted and have the sketch resume in a new location where the pen touches the surface?** OR **How can I draw an ink sketch and snap to another ink sketch start/end node?** I am off to give it this a try in QGIS...


title: ArcGIS Heads up digitizing - Rooms and doors

text: I am digitizing the floorplans of various schools in our county. This is for the Fire Department so they want to be able to see where doors are within the building. The doors do not necessarily have to be indicated; a gap in the wall is sufficient. It would be much more efficient for me to use the 90 degree tool and knock out these classrooms in full. Is there some sort of eraser feature that might split a line segment in two wherever it intersects? That way I could basically create a rectangle for the room and erase where the door might be to leave a gap. Or are there any other tools that could make this process easier? I've got a few dozen schools to digitize, and by far the most time consuming process is creating small segments around the doors just to leave that gap. Cheers!


title: OpenLayers crashes stock browser on Android 4.0.4 devices

text: * * * **Update:** The issue appears to be with OpenLayers Release 2.13 dev. I mistakenly mentioned 2.12 as being problematic but have retested with no problems. Also, no problems with 2.13rc3 ## **Update 2:** This issue appears to exist again in 2.13.1 We're developing a mapping tool using Openlayers ~~2.1.2 (issue also confirmed on~~ 2.1.3 ~~).~~ I've tested on 4 devices running Android 4.0.4. The map panning/zooming works for a few seconds before the browser completely locks up - no longer able to pan/zoom or navigate to another page. The browser then either crashes or I need to force stop in application settings. Devices running 4.1.1 have no problems. This is also an issue with OpenLayers own mobile example page: http://openlayers.org/dev/examples/mobile.html Again, after a few seconds, the browser freezes. I see no reports of this online. The OpenLayers documentation says it is supported on Android 2.1+ Has anyone else run into this issue?


title: What is difference between simultaneous users and simultaneous connections in ArcGIS for Server 10.1?

text: What is the difference between simultaneous users and simultaneous connections. We're looking at Workgroup standard vs Enterprise Standard. Here's what we want to do: Use Arc Gis Server to serve data to app. 30 vehicles using Ipads via Sierra Wireless 4g modems. This data would be maps for Right of Way, outage locations etc. We would also like to be able for users to send data to other users. Thanks, Russell


title: Is it possible to use UTF-8 characters for shapefile attributes?

text: I heard that according to a specification shapefile attributes stored in dbf- files cannot have UTF-8 as encoding. And then I found out that an additional file with the ending CPG can be used to define the code page to identify the character set used in the dbf file. For me this sounds contradictory. So my question is now - is it possible to use UTF-8 characters for shapefile attributes? and - additional question, since encoding and character sets are always an unpopular topic, do You know a page where the topic is described in detail?


title: Best way to get the overlapping segment of two lines

text: Could anyone help me figure out the best way to get the overlapping line segment of 2 line sdo_geometry objects? my query looks like this:               select SDO_GEOM.SDO_INTERSECTION(SHAPE, :SHAPE, 0.5) as geom     from MyTable      where MDSYS.SDO_OVERLAPS(SHAPE, :SHAPE) = 'TRUE'      the geometry objects in the table are all lines as is eth parameter being passed in, but the geometry being returned for all results is a point, not a line.


title: travel survey multi stop routing problem

text: I have a travel diary survey dataset that I am attempting to transform into more useful information. Currently, the data comes in the format of:               Survey# ,  PlaceVisited#,  TravelMode, TripPurpose, LocationType          1,      1,       NA,        NA,         home          1,      2,       auto,  doctor,     office          1,      3,       auto,  shop,       grocery_store          1,      4,       auto,  NA,     home                      where each place the survey respondent visits is a separate record (row) in the survey table, but this is really one multipurpose trip. What I would like to do with the data is separate each actual trip into a route with various points and purpose information included as (i think) a line using shortest path on a road network (I don't have the route actually taken). Currently, if a respondent leaves their house, and travels to several different stores by car, each store is entered as a separate record, as is their home starting location. What I would like to have is one "trip" for each time the respondent leaves their house and comes back, regardless of how many places they visit while preserving the information about each place they visit. I am wondering what the best way to tackle this problem would be. I have so far struck out finding prepackaged tools to perform such a task in QGIS, GRASS, ArcGIS, Maptitude, and PostGIS. Should I approach this as a scripting problem using SQL, scripting a shortest path tool in a GIS, or is there another more efficient way to do it. Many thanks in advance!


title: OSM XML: Should I generate contact:phone or phone?

text: I am generating OSM XML data for all POIs of Wikivoyage, the collaborative travel guide. When I want to express a phone number, should I generate `contact:phone` or `phone`? Both mean the same thing. According to a 2012 note on this page, the only difference is that `contact:phone` is newer but `phone` is more popular. Which one should I use, that so my data is readable by most 2013 apps?


title: how can the ogr2ogr help to merge a shp file?

text: I have over 30 shp files which is named "street.shp" into the sub-folder. There is over 30 sub-folders. In the sub-folder, there is another shp files which may be named a building.shp...etc. I only want to merge the street.shp files. The structure of all street.shp is exactly same! I tried this script -               for /R %f in (*.shp) do ogr2ogr -update -append merge.shp %f -f “esri shapefile” -nln merge       How can the script find the shp file named "street.shp" only? Thanks a lot.


title: Trigger a custom ArcToolbox tool outside Model Builder

text: I have created a model in a custom Arc Toolbox with Model Builder in ArcMap 10.1. Is it possible to trigger the model outside Model Builder and ArcMap? For example, use a Python script to trigger the model to operate?


title: ArcToolbox model run but warning output in Python

text: I have created an ArcToolbox and a Model in the toolbox. It uses a lot of Spatial Analyst tools and deals with raster images. When I imported the toolbox in a Python script and ran the model from Python in Windows command prompt, I got a few warning messages as below but the model got run successfully.               <string>:1: SyntaxWarning: import * only allowed at module level      It was shown several times. I am unable to tell where this error comes from. Any idea? Cheers


title: How to identify valleys on DEM in ArcGIS 10.1

text: First, What is the best and easiest way to identify valleys in a DEM raster? Second, how do you create a least-cost pathway for each valley in a region based on low angle slope and proximity to rivers? Any help would be massively appreciated!!!


title: SAGA raster calculator, ifelse command syntax problem?

text: I am trying to get the saga raster calculator in QGIS to spit out all values of a raster that are inferior to 10, and make all other values zero. I think the syntax is correct but I keep getting the original raster as an output which is definitely not what I want. Here is the code I am using:                   ifelse(lt(a,10),a,0)      I'm thinking there must be a problem with the way I am referencing the raster? I read in the SAGA help pdf that "if you have four grid map layers in your ‘Input’ list, they will be referred to by a, b, c, and d." So I'm not sure what I'm doing wrong.


title: Coordinate conversion (EPSG:3857 to EPSG:4326) using OpenGIS/JTS too slow

text: I am using the Java based libraries of OpenGIS and JTS to convert coordinates from EPSG:3857 to EPSG:4326 to proceed to compute the distance using the Haversine formula. The distance features in a constraint on whether to zoom in further on a map or not. We are using a tangible user interface and depend on the reactivity of the system. Unfortunately, the conversion is quite slow which induces even more delay than loading the map from a server. Would there be a less costly way to proceed with the conversion using local resources?                       CoordinateReferenceSystem targetCRS = CRS.decode("EPSG:4326");             CoordinateReferenceSystem sourceCRS = CRS.decode("EPSG:3857");             MathTransform transform = CRS.findMathTransform(sourceCRS, targetCRS, false);             GeometryFactory geometryFactory = new GeometryFactory(new PrecisionModel(), 3857);             Point minPoint = geometryFactory.createPoint(new Coordinate(xMin, yMin));             Point minTargetPoint = (Point) JTS.transform(minPoint, transform);             Point maxPoint = geometryFactory.createPoint(new Coordinate(xMax, yMax));             PointmaxTargetPoint = (Point) JTS.transform(maxPoint, transform);


title: Problem with precision when convert from WKT to shapefile

text: I have a problem with converting from WKT to shapefile in pygdal. In wkt I have fields with precision (for example real 7.2). When I create a shapefile I have field numeric with no precision (numeric 7.0). I can't find problem. Do you have any suggestions?


title: How to standardise colour-value graduations across multiple maps

text: In QGIS I have a vector dataset (imported csv vector + csvt) spanning multiple years, and a shape file shape which is essentially a grid. By joining the csv file with the shape file I can apply a colour to each grid square and hence a gradient scheme for the whole grid. I have created a series of maps - one map per year, each with its own colour gradient. However the colour gradients are not comparable between maps because the range of values is different among years. How can I make the gradient scales standardised across all maps? I've tried splitting the csv file into years, and with each year individually, it does not seem possible to create a 'custom' range of values for the gradient, say 0-100, even tho the values span say 20-70. I've also tried working with one csv file containing all years, thinking that I could somehow 'trick' QGIS into applying a gradient that spans the range of the whole dataset, while at the same time only map individual years. Categories doesn't seem to be an option as if I have 200 grid squares I get 200 categories each with own unique value. Any help would be much appreciated. It may be that there is a database-based solution which I am open to. I am also open to GRASS solutions, but I am less familiar with GRASS as this stage. Many thanks! Paul


title: Are there GIS tools for simulating landscape fragmentation?

text: I am currently looking for any GIS packages (or add-ons) that allow the simulation of landscape fragmentation on raster data. Maybe also other software than GIS tools? Does anyone of you have had experience with that issue already and can suggest a good tool for that task? Info on tutorials etc. are also welcome! ThanX! Jens


title: Map Refresh in Infrastructure Map Server 2015

text: I have a MapGuide application build in MapGuide 2011 where I refreshed my map using the following script               function MapRefresh() {            debugger;            ExecuteMapAction(20);             parent.parent.mapFrame.Refresh();              }          Now I am upgrading my application to Autodesk Infrastructure Map Server 2015 where the above script does not work, Please guide me as to how I can refresh my main map from a panel?


title: How to get elevation profile from .tif dsm on the web map?

text: I have my dsm .tiff file, pixels are floating point 32bit values giving the terrain elevation in metres. I need to create a web map where user can draw a path and get elevation profile for that path. Raster has attribute 'GRAY_INDEX' where values are stored. What are the options to achieve this and is it even possible? I thought about Geoserver `GetFeatureInfo`, but I'm not sure how that would work as I can get values for only one point at a time. I thought when user draws a path, I make a number of `GetFeatureInfo` requests exactly as points in that path, but I assume that would lead to slow performance. Is there a better option?


title: Custom style a particular road using CartoCSS/TileMill/OSM

text: For a particular city, I want to make the label for its main street labeled as prominently as Motorways. i.e. although this road would normally only be labeled from about zoom level 15 and up, I would like it to be visibly labeled from zoom level 12. Is there a syntax I can use in TileMill's CartoCSS to allow me to custom adjust label properties for this particular road? Ideally it could be based on road name (e.g. "Mary Street") but pretty much any solution could work such as using its ID number, because we only have to do it once.


title: GeoServer + ImagePyramid: Request is empty where data exists

text: I have set up an _image pyramid_ based on 5 levels of tiling. The data is served using the _ImagePyramid plugin_ for _GeoServer_. When I do a _layer preview_ the data is visualized just like I expect them to be. But when I zoom in a couple of levels, I get a _white_ tile in return. The _GeoServer log_ tells me the following: > 2012-05-14 08:07:07,495 DEBUG [gce.imagemosaic] - Reading mosaic from > file:/e:/GeoServer/data_dir/data/pyramid/2/2.shp   >  2012-05-14 08:07:07,495 DEBUG [gce.imagemosaic] - Highest res > 0.6400000000000005 0.6400000000000005   >  2012-05-14 08:07:07,495 DEBUG [gce.imagemosaic] - The computed > CropBoundingBox ReferencedEnvelope[534687.5 : 535156.25, 6225625.0 : > 6226093.75] Doesn't intersect the coverage BoundingBox > > ReferencedEnvelope[-1.0 : 0.0, -1.0 : 0.0] resulting in an empty request   >  2012-05-14 08:07:07,495 DEBUG [gce.imagemosaic] - RequestedBBox empty or > null   >  2012-05-14 08:07:07,495 DEBUG [gce.imagemosaic] - Request is empty: > RasterLayerRequest description:   >      RequestedBBox=ReferencedEnvelope[534687.5 : 535156.25, 6225625.0 : > 6226093.75]   >      RequestedRasterArea=java.awt.Rectangle[x=0,y=0,width=768,height=768]   >      RequestedGridToWorld=PARAM_MT["Affine",   >          PARAMETER["num_row", 3],   >          PARAMETER["num_col", 3],   >          PARAMETER["elt_0_0", 0.6103515625],   >          PARAMETER["elt_0_2", 534687.8051757812],   >          PARAMETER["elt_1_1", -0.6103515625],   >          PARAMETER["elt_1_2", 6226093.444824219]]   >      ReadType=JAI_IMAGEREAD   >  2012-05-14 08:07:07,495 DEBUG [gce.imagemosaic] - The response is empty. > ==> returning a null GridCoverage   > If I zoom in a couple of levels further, I get data again. The extent of the data is [530000, 6220000, 540000, 6230000] and the projection EPSG:25832, which the GeoServer admin interface agree on. **Then what could be the reason for GeoServer to look inside the BBOX [-1, 0, -1, 0], when my data is inside the requested map?**


title: osm2pgsql missing coordinates

text: I am quite new to GIS and OSM and am trying to get my head around exporting custom OSM region and uploading them into PostgreSQL. I have started with a very small map area to export with only two shapes - an island and a meadow - (http://www.openstreetmap.org/#map=19/43.86436/15.34266) which export the following .osm data:               <?xml version="1.0" encoding="UTF-8"?>     <osm version="0.6" generator="CGImap 0.3.3 (23630 thorn-02.openstreetmap.org)" copyright="OpenStreetMap and contributors" attribution="http://www.openstreetmap.org/copyright" license="http://opendatacommons.org/licenses/odbl/1-0/">         <bounds minlat="43.8635800" minlon="15.3411900" maxlat="43.8650000" maxlon="15.3445500"/>         <node id="63910084" visible="true" version="3" changeset="17709161" timestamp="2013-09-06T20:07:53Z" user="idelac" uid="183676" lat="43.8644519" lon="15.3432012"/>         <node id="63930603" visible="true" version="3" changeset="17709161" timestamp="2013-09-06T20:07:54Z" user="idelac" uid="183676" lat="43.8642369" lon="15.3431567"/>         <node id="63968158" visible="true" version="3" changeset="17709161" timestamp="2013-09-06T20:07:57Z" user="idelac" uid="183676" lat="43.8645365" lon="15.3424617"/>         <node id="1122490202" visible="true" version="2" changeset="17709161" timestamp="2013-09-06T20:07:44Z" user="idelac" uid="183676" lat="43.8640377" lon="15.3426171"/>         <node id="63879158" visible="true" version="3" changeset="17709161" timestamp="2013-09-06T20:07:52Z" user="idelac" uid="183676" lat="43.8639894" lon="15.3430113"/>         <node id="2448417703" visible="true" version="1" changeset="17709161" timestamp="2013-09-06T20:06:52Z" user="idelac" uid="183676" lat="43.8641440" lon="15.3427266"/>         <node id="2448417704" visible="true" version="1" changeset="17709161" timestamp="2013-09-06T20:06:52Z" user="idelac" uid="183676" lat="43.8641542" lon="15.3426165"/>         <node id="2448417706" visible="true" version="1" changeset="17709161" timestamp="2013-09-06T20:06:52Z" user="idelac" uid="183676" lat="43.8642112" lon="15.3428451"/>         <node id="2448417707" visible="true" version="1" changeset="17709161" timestamp="2013-09-06T20:06:52Z" user="idelac" uid="183676" lat="43.8642458" lon="15.3424979"/>         <node id="2448417708" visible="true" version="1" changeset="17709161" timestamp="2013-09-06T20:06:52Z" user="idelac" uid="183676" lat="43.8643394" lon="15.3428931"/>         <node id="2448417709" visible="true" version="1" changeset="17709161" timestamp="2013-09-06T20:06:52Z" user="idelac" uid="183676" lat="43.8644045" lon="15.3425008"/>         <node id="2448417711" visible="true" version="1" changeset="17709161" timestamp="2013-09-06T20:06:52Z" user="idelac" uid="183676" lat="43.8645022" lon="15.3428874"/>         <node id="2448417712" visible="true" version="1" changeset="17709161" timestamp="2013-09-06T20:06:52Z" user="idelac" uid="183676" lat="43.8645144" lon="15.3425403"/>         <node id="2448417713" visible="true" version="1" changeset="17709161" timestamp="2013-09-06T20:06:52Z" user="idelac" uid="183676" lat="43.8645937" lon="15.3427745"/>         <node id="2448417714" visible="true" version="1" changeset="17709161" timestamp="2013-09-06T20:06:52Z" user="idelac" uid="183676" lat="43.8645978" lon="15.3426673"/>         <node id="1122490153" visible="true" version="2" changeset="17709161" timestamp="2013-09-06T20:07:43Z" user="idelac" uid="183676" lat="43.8645802" lon="15.3431461"/>         <node id="1122490171" visible="true" version="2" changeset="17709161" timestamp="2013-09-06T20:07:43Z" user="idelac" uid="183676" lat="43.8646997" lon="15.3426923"/>         <node id="1122490322" visible="true" version="2" changeset="17709161" timestamp="2013-09-06T20:07:44Z" user="idelac" uid="183676" lat="43.8647101" lon="15.3429354"/>         <node id="1122490460" visible="true" version="2" changeset="17709161" timestamp="2013-09-06T20:07:45Z" user="idelac" uid="183676" lat="43.8641347" lon="15.3424630"/>         <node id="1122490524" visible="true" version="2" changeset="17709161" timestamp="2013-09-06T20:07:45Z" user="idelac" uid="183676" lat="43.8642836" lon="15.3424149"/>         <way id="8837893" visible="true" version="3" changeset="7098380" timestamp="2011-01-26T21:37:50Z" user="goranT" uid="217969">             <nd ref="63879158"/>             <nd ref="63930603"/>             <nd ref="63910084"/>             <nd ref="1122490153"/>             <nd ref="1122490322"/>             <nd ref="1122490171"/>             <nd ref="63968158"/>             <nd ref="1122490524"/>             <nd ref="1122490460"/>             <nd ref="1122490202"/>             <nd ref="63879158"/>             <tag k="natural" v="coastline"/>             <tag k="place" v="island"/>             <tag k="source" v="PGS"/>         </way>         <way id="236881851" visible="true" version="1" changeset="17709161" timestamp="2013-09-06T20:07:26Z" user="idelac" uid="183676">             <nd ref="2448417709"/>             <nd ref="2448417707"/>             <nd ref="2448417704"/>             <nd ref="2448417703"/>             <nd ref="2448417706"/>             <nd ref="2448417708"/>             <nd ref="2448417711"/>             <nd ref="2448417713"/>             <nd ref="2448417714"/>             <nd ref="2448417712"/>             <nd ref="2448417709"/>             <tag k="landuse" v="meadow"/>         </way>     </osm>      Running `osm2pgsql -S default.style -d gis -U postgres island.osm` uploads the data into my gis database with one row in the planet_osm_line table (for the island) and one row in the planet_osm_polygon table (for the meadow). Now, when issuing sql query `SELECT ST_AsText(way) FROM planet_osm_line`, pqsql returns               "LINESTRING(1707992.39 5444456.51,1707997.34 5444489.7,1707915.02 5444502.76,1707932.32 5444425.75)"      Here I don't understand two things- **Questions**   1. Why does the Linestring only contain four points when the original coastline of the island has 10 distinct nodes? Similarly when issuing the sql query on the planet_osm_polygon table to see the polygon for the meadow shape, the polygon has fewer coordinates than its original shape described by the nodes (7 instead of 11).   2. _(answered just below in update) How are the original lat,lng coordinates projected? I know osm2pgsql defaults to mercator projection. However, I don't understand how it computed the values in the Linestring above. Is the width of a tile (256px?) and some zoom (z=?) taken into account in the projection? Would LatLng[0,0] be projected as a Point[0,0]?_ **Update** The answer to my second question is that osm2pgsql projects from latitude and longitude coordinates to mercator in meters. LatLng[0,0] will indeed be projected as Point[0,0] where the point coordinates are the distance in meters from the origin on a mercator plane. Here is the projection function in case someone stumbles upon this question and looks for a solution               var degrees2meters = function(lon,lat) {             var x = lon * 20037508.34 / 180;             var y = Math.log(Math.tan((90 + lat) * Math.PI / 360)) / (Math.PI / 180);             y = y * 20037508.34 / 180;             return [x, y]     }      I still don't know why some of the nodes are missing in the resulting shapes in the database. Has anyone experienced this?


title: Viewing map tiles generated from Maperitive in leaflet or openlayers

text: for few days i have been working on developing my own osm tile server. I tried it on linux as well as on windows, no complete success till now. I Have to make tile server on both windows and linux. for windows i used Maperitive and succeeded in generating tiles and then i used this guys code               http://gis.stackexchange.com/questions/84697/viewing-map-tiles-from-maperitive-in-leaflet-or-openlayers      now the prob is tiles are not showing up Do i need to instal some Apache module for it ? My tiles are kept here               C:\Profiles\Maperitive\Tiles      Where as i am opening html file from server               http://localhost:8080/mwork/maperitive_test1/index_leaf.html      which is kept in here               C:\Profiles\ApachePHP\apache\www\mwork\maperitive_test1


title: OGR: How to best use GetFeature on a selection?

text: I think the following is very confusing               sel_str = "ND_14 = '108'"     l.SetAttributeFilter(sel_str)     l.GetFeatureCount()     >>> 1       In a for-loop this gives me               for f in l:        g = f.GetGeometryRef()        print g.GetEnvelope()     >>> (8.39746, 8.95788, 46.52757, 46.99331)      Whereas GetFeature() gives me               >>> feat0 = l.GetFeature(0)     >>> feat1 = l.GetFeature(1)      >>> g0 = feat0.GetGeometryRef()     >>> g1 = feat1.GetGeometryRef()     >>> g0.GetEnvelope()     (8.21206, 10.65409, 56.55045, 57.75181)     >>> g1.GetEnvelope()     (10.60014, 10.6172, 57.48933, 57.4979)      Apparently `GetFeature()` is based on the FID of the whole set, not the selection. But using a for-loop on a selection of length one is pretty ridiculous. How do I get the FID of my selected feature then? Or what are other ways of using `GetFeature()` with your current selection?


title: data range changes when warping

text: When the source band having pixel values in the range of 52 and 207 is warped using the following command:               gdal_warp -t_srs wgs84 -r near in_file out_file      The resulted band has the pixel values in the range of 0 and 207. How can I make the data range unchanged during warp?


title: QGIS 2.4: Basemap.at is displayed in poor quality

text: I have a huge problem displaying the austrian Basemap.at WMTS Service in QGIS 2.4. When I add the service, qgis logs multiple errors on failed tile requests. When I zoom into the map, it is displayed in a very poor quality (see screenshot).   I really don't know how to fix this problem - can anyone tell me how to increase the quality of the tiles? WMTS URL: http://maps.wien.gv.at/basemap/1.0.0/WMTSCapabilities.xml ![Basemap.at Scrennshot: Poor Quality](http://i.stack.imgur.com/gDW9M.png) UPDATE 1:   I tried setting the resampling methods in QGIS to Cubic/Average, as you guys advised me to and compared the result with the same zoomlevel in ArcGIS. In this screenshot you can see the result: At the same scale level QGIS just skips an entire tile level, whereas arcgis loads this level correctly. - Any other ideas? ![enter image description here](http://i.stack.imgur.com/YJJMs.jpg)