title: Is there a way to dynamically hide annotations that intersect map extent border?

text: I have a feature linked annotation feature class with a number of annotations displayed only partially on the edges of the visible extent of the layout data frame. The link annotation class will be used in a number of other maps of the same scale and paper size but showing different locations. **Is there a way to dynamically hide the partially displayed annotations?** _Out of the box solution preferred, programmatic (ArcObject, VB.Net Addin) acceptable_ ![enter image description here](http://i.stack.imgur.com/rl0iZ.jpg)


title: How to label a lot of nearby points?

text: I am currently working with a lot of points of which every point has to be labeled. I am already using the "Easy Custom Labeling" tool to move the labels and create lines to connect points with labels, as soon as there is not enough space for the labels. Now I have different problems:   1. It takes a lot of work to align the labels manually to reach the goal to label every feature. Is there a faster way to to it? Very often, I let one label disappear, as soon as I move another. I think in ArcGIS there was the option to convert labels to annotations and then align these annotations very close to each other. The labels in QGIS disappear as soon as they are to close together.   2. As soon as I go to the Composer Manager, some labels arrange differently although I choose the same scale as I used when moving the labels manually. This means that I have to arrange the labels manually in my data view, switch to the Composer Manager and refresh the Map to see, whether it finally worked. This process takes even more time. I am happy for every hint to label more effectively! Some maps have a much higher density of points, I just started with Italy. ![Italy](http://i.stack.imgur.com/t9FRc.jpg)


title: Generate tiles for google maps overlays

text: I am trying to show some information on Google Maps which would be in the form of translucent circle overlays. The problem is that some circles overlap, which results in increased opacity and blending of colors if circles are of different color. I searched a lot and came to conclusion that there is no possible way to avoid color blending. The alternate solution I was recommended was to generete png tiles on my server and display them or decide color of each pixel on server side and use HTML5 canvas and render each pixel. I searched a lot but couldn't find any way to do so. It would be very helpful if someone could point me in the right direction. **I am rephrasing my question below to be more clear and concise** I have to display circles on a map. Radius and color of these circles depict their area of influence and type of influence. These circles are translucent but overlapping should result in only the top one being visible without any blending. How do I generate tiles for this information? **Probably I am not clear so I am adding further details.** Check this jsFiddle http://jsfiddle.net/Dv4UT/ .There are two overlapping circles green(zIndex : 10) and red(zIndex : 11) with opacity .5 .The overlapping part is a mixture of red and green while I want the overlapping part to be of top color with same opacity ie:red.If this is possible using map overlays then it would be great.Else I think it is possible using HTML5 canvas but that won't work on older version of IE.So finally I would have to fall back on generating tiles on server side and using them as overlays.I don't have any idea of generating tiles and it would be very helpful, if someone could point me in the right direction.I am attaching a screenshot from trulia commute map, check that circles beneath are being overwritten by top ones.I need similar kind of visualization. ![example image](http://i.stack.imgur.com/fOTUI.png)


title: setLayerDefinitions without defining layers?

text: I am working on a similar project to How to zoom to feature when clicking on drop-down list box item with ArcGIS JS API?, trying to filter the graphics on my map based on the user input into a FilteringSelect. The overall methodology of how my app is set up is that I go from a queryTask into a forEach loop that calls map.graphics.add(graphic) for each encountered feature returned from the queryTask. I was hoping to be able to use the map.getLayer("dynamic").setLayerDefinitions(...) function to filter the map results based on user input, but from what I can tell that's not possible here since I don't define layers. Is there any way to accomplish this filtering with the setup I am currently using?


title: How to speed up the plotting of polygons in R?

text: I want to plot the country borders of North America over a raster image depicting some variable and then overlay contours on top of the plot using R. I have been successful in doing this using base graphics and lattice, but it seems that the plotting process is much too slow! I have not done this in ggplot2 yet, but I doubt that it will fare better in terms of speed. I have the data in a netcdf file created from a grib file. For now, I downloaded the country borders for Canada, USA and Mexico, which were available in RData files from GADM which read into R as SpatialPolygonsDataFrame objects. Here is some code:               # Load packages     library(raster)     #library(ncdf) # If you cannot install ncdf4     library(ncdf4)          # Read in the file, get the 13th layer     # fn <- 'path_to_file'     r <- raster(fn, band=13)          # Set the projection and extent     p4 <- "+proj=lcc +lat_1=50.0 +lat_2=50.0 +units=km +x_0=32.46341 +y_0=32.46341 +lon_0=-107 +lat_0=1.0"     projection(r) <- CRS(p4)     extent(r) <- c(-5648.71, 5680.72, 1481.40, 10430.62)          # Get the country borders     # This will download the RData files to your working directory     can<-getData('GADM', country="CAN", level=1)     usa<-getData('GADM', country="USA", level=1)     mex<-getData('GADM', country="MEX", level=1)          # Project to model grid     can_p <- spTransform(can, CRS(p4))     usa_p <- spTransform(usa, CRS(p4))     mex_p <- spTransform(mex, CRS(p4))          ### USING BASE GRAPHICS     par(mar=c(0,0,0,0))     # Plot the raster     bins <- 100     plot(r, axes=FALSE, box=FALSE, legend=FALSE,          col=rev( rainbow(bins,start=0,end=1) ),          breaks=seq(4500,6000,length.out=bins))     plot(r, legend.only=TRUE, col=rev( rainbow(bins,start=0,end=1)),          legend.width=0.5, legend.shrink=0.75,           breaks=seq(4500,6000,length.out=bins),          axis.args=list(at=seq(4500,6000,length.out=11),                     labels=seq(4500,6000,length.out=11),                     cex.axis=0.5),          legend.args=list(text='Height (m)', side=4, font=2,                            line=2, cex=0.8))     # Plot the borders     # These are so slow!!     plot(can_p, add=TRUE, border='white', lwd=2)     plot(usa_p, add=TRUE, border='white', lwd=2)     plot(mex_p, add=TRUE, border='white', lwd=2)     # Add the contours     contour(r, add=TRUE, nlevel=5)          ### USING LATTICE     library(rasterVis)          # Some settings for our themes     myTheme <- RdBuTheme()     myTheme$axis.line$col<-"transparent"     myTheme$add.line$alpha <- 1     myTheme2 <- myTheme     myTheme2$regions$col <- 'transparent'     myTheme2$add.text$cex <- 0.7     myTheme2$add.line$lwd <- 1     myTheme2$add.line$alpha <- 0.8          # Get JUST the contour lines     contours <- contourplot(r, margin=FALSE, scales=list(draw=FALSE),                             par.settings=myTheme2, pretty=TRUE, key=NULL, cuts=5,                             labels=TRUE)          # Plot the colour     levels <- levelplot(r, contour=FALSE, margin=FALSE, scales=list(draw=FALSE),                         par.settings = myTheme, cuts=100)          # Plot!     levels +         layer(sp.polygons(can_p, col='green', lwd=2)) +       layer(sp.polygons(usa_p, col='green', lwd=2)) +       layer(sp.polygons(mex_p, col='green', lwd=2)) +       contours      Is there a way to speed up the plotting of the polygons? On the system that I am working on, the plotting takes several minutes. I want to eventually make a function that will easily generate a number of these plots for inspection, and I reckon I will be plotting many of these maps, so I want to increase the speed of the plots! Thanks!


title: Basics about Spatial Indexes

text: I was asking on a forum if it could be a good idea to create spatial indexes on subsets of a table when your application display on map only those subsets and never the entire table. I was asking this question because I though that as the subsets do not have the same extents of the whole table, maybe it would be faster to display the subsets with their own spatial indexes. **The answer I received was that spatial indexes do not affect display time but only used for spatial queries such union or intersect. Is it true???** My experience with GIS and database is when a table has no spatial index, the display on map is much slower. I always thought that on display the table was queried to show the features that intersect the current map window extent, so the features outside are not loaded for nothing. Is it really how it works? This is a kind of spatial query no? What is the truth? Is it a good idea to create spatial indexes on subsets?


title: No legend in QGis Web Client

text: I have a QGIS Web Client working, but it lacks a legend and metadata panel as is shown here What do I need to have the referred legend panel?


title: How to smooth a raster in a rectangular fashion using gaussian blur?

text: How do I adjust the code for the Gaussian Blur posted in the thread: What raster smoothing/generalization tools are available? to retain the resolution/gradients in the y direction and blur only in the x direction (rather than over a square area)?


title: georeference IMG file in qgis

text: I have a problem: I'm trying to georeference a satellite image (ASTER, Erdas Imagine file *.img) in QGis. I have georeferenced with QGis before and never had any problems. The file is georeferenced already in UTM 17N WGS 84 but very very badly, so I need to re-do it. I have opened the file in QGis and set the project coordinate system to UTM 17 N, then I have added a basemap with open layers and used the basemap (Bing Maps) to generate GCPs using the Georeferencer tool. I set the output coordinate system to UTM 17 N, too. Everything looks good, but when I finish the georeferencing the file ends up in the wrong place (or in no place at all rather, the basemap disappears when I zoom to the layer). Sometimes it doesn't even display. I have tried different coordinate systems, re-sampling methods, also saved the file as .tiff - no changes. I have tried it with other raster files, they always worked well. Does anybody have any idea what I'm doing wrong? (Just in case: I have tried with ArcGis, that doesn't work because I don't have a BingMap access key and the basemap I can access has a worse resolution than my ASTER image. Same problem in Erdas, theoretically I could use Google Earth for georeferencing there but I'm using a computer at my university and the Google Earth tool isn't working). I'd be very very glad if anybody has a suggestion, this is for my bachelor thesis and quite important, and I've spent two days already trying... :)


title: Looking for a way to change the default shape file label in ArcGIS 10

text: Is there a way to change ArcGIS's default shapefile label `Font: Aerial Size: 8` to something we are using over and over again? As shown in attached screen shot, i'd like to include a label expression in addition to a saved text symbol:Shape file labels example![enter image description here](http://i.stack.imgur.com/JmQAB.jpg)   It would be very helpful to save the couple of steps it takes to change labels as i am working with hundreds of maps / shape files.


title: Is there a comprehensive list of latitude and longitude values for cities around the world?

text: I'm developing an application that calculates the distance between cities using their respective Longitude and Latitude values. But the problem is where to get a comprehensive list of latitude and longitude values of these cities?


title: Where can I download the Google Road network as shapefile?

text: Is there any site from where I can get Google Road network as shapefile? As per my knowledge It is not distributed by google. So anybody can suggest how can I create this road network for a custom use?


title: GMT grid to ESRI readable format

text: I am trying to convert a GMT grid data into a format that is readable in ArcGIS. I have already attempted to load the data into GRASS and use GDAL translate to convert into a tif with unwanted results. Mainly the extent is wrong with a cell size of 1,1 while in reality it should cover the entire world. _**Update_** After installing GMT and utilizing the source code               grd2xyz -Ef input.grd      the resulting ascii grid was converted to a tif format that was read by ArcGIS (ascii to raster). However again the final product does not have a correct extent showing a cellsize of 1,1. Could someone who perhaps is more familiar with GMT check the data above to check if it is viable to use within an ArcMap format?


title: leaflet: Zooming a tile cache

text: I have a NASA BlueMarble tile layer that goes from zoom factors 0-8. I when I go to zoom level 9, I do see a momentary scaling to 9 (stretching pixels). But then, as expected, I get the errorTileUrl. Nothing wrong here. But rather than have to created the tile cache (and massive data store) for zoom factors 9-14 or so, how would one go about having those zoom factors be pixel scaling and stretching? This would be very useful for my offline maps. Ideas? (BTW, any suggestions for how to store offline tile caches? Say 2K to 4K of tiles in PNG format? Would I use an app manifest? (web application)?


title: Filter ArcSDE layers of the TOC into a combobox

text: I've a list of layers of the current workspace in a combobox, but I want to list only the layers of the ArcSDE Connection, not the others I've in the TOC. How can I filter the name of this layers?


title: How to show PostGIS distance in Meter

text: I have calculate distance from one point to a linestring. I have used the following query               WITH data AS (     SELECT geom::geometry AS road_network from linestring)          SELECT  ST_Distance(ST_Line_Interpolate_Point(road_network, ST_Line_Locate_Point(road_network,ST_GeomFromText('POINT(90.402 23.753)'))), ST_GeomFromText('POINT(90.402 23.753)'))  as Dist                FROM data          ORDER BY Dist     Limit 1      Its returning the following result               0.000501079237817185      I want to get the result in Meter. How can I convert it? I am using WGS84 projection system. I know the `ST_Transform` function but I need a perfect query where I use this function. Thanks


title: how can convert geometry type to geometry(Point)?

text: I think my question is stupid but I'm blocking. I just work with postgis ST_Line_Locate_Point function (linestring, point). to find the ends of lines. My problem is that my objects geometry type, not point. So how can I convert the geometry to geometry (point) ..?


title: Google maps and Ordnance Survey derived data licence restrictions

text: **Current Status:** You can't use data derived from restricted licence OS data with Google Maps. There's been an on-going question about whether it is legal to display data derived from Ordnance Survey licenced mapping via google maps. By derived data, I mean geometries which have been copied exactly from, or are composed of parts from, another dataset. The blog post here: http://blog.ordnancesurvey.co.uk/2010/08/what-the-psma- really-means/ outlines the issue. We've spoken to Google and they say that the OS have recently relaxed their licencing, but haven't been able to get any comment from the OS. The last post on the OS blog is from November 2012. Does anyone have any ideas as to how the land lies here? Particularly if anyone has heard anything from the OS on this subject since November 2012. **UPDATE:** 16th April 2013: We've got an account manager at the OS now and a collegue has just spoken to him. Long story short: nothing has changed since November 2012 - any data derived from OS licenced mapping which is not specifically "free" (such as OS OpenData) may not be displayed in Google maps. Our account manager didn't know when the OS was going to come back with an answer. He said he'd chase it up. `*starts holding breath*` **UPDATE:** 2nd July 2013: Chasing our account manager, they reported the following (loosely quoted, do not base decisions solely on this) "The OS's position is unchanged and is not being pursued further. It is our (the data user's) responsibility to check with Google, or any other API provider, to ensure that OS rights are preserved. Rights to derived data cannot be given over to Google, as their T&C's require". Now, talking to Google, they don't think their T&C's violate the OS licence. I'm not sure this is too useful though. Our account manager has promised to send written clarification, which I will post here: a) If it ever arrives and b) if it isn't confidential. The OS blog hasn't been updated yet. **UPDATE:** 2nd July 2013: Just had an email from our account manager which contradicts what we were told earlier: "...have had a chat with our lawyers who confirmed that we are in discussion with Google around the compatibility of their terms with ours. Currently, the position is that Ordnance Survey map data cannot be overlaid on top of Google Maps though we will not have an “official document/wording” available until the current discussions are concluded.". So there! **UPDATE:** 6th September 2013: An email to a colleague from the boss of our account manager includes the line: "...Currently our data cannot be overlayed on top of Google Maps due to their rules around IP ownership". I don't know if this is now their official legal stance or this guy's opinion.


title: How to get the nearest line to a feature?

text: I have features on map, and i want to get the nearest poly-line to each one. How to write spatial query to get the nearest one. I tried to use `ISpatialFilter` but i couldn't reach anything.


title: Vector Analysis: Determining area of county land mass without water features

text: # Here are the layers that I have:   * State / Province polygon file (With Area Field)   * County / Ward polygon (with no area field), but all the polygons are within the state boundaries   * Water ways for the whole province / state. (includes rivers, lakes, creeks, brooks, etc) with no area fields # Here is what I want to do:   * Determine the gross land area of counties without the rivers and lakes. ## Notes   * The area doesn't have to be extremely precise (nearest 5 sq km is fine with me)   * Has to work on ArcMap or QGIS


title: Assign feature actions in QGIS--too many actions run at the same time

text: Using QGIS, I'm assigning actions to a polygon layer whose features are based on extents of hundreds of maps of a variety of files (including DWG, DGN, PDF, jpg, etc). However, since some of the polygons overlap--e.g. those with smaller extents are within those with larger extents or some even have the same extent. Is there a way that I can control which folders or files are opened when clicking a feature that overlaps another? Currently, there are a ton of actions being called at the same time (folders and files being opened). I've looked into ideas from this post: reorder a polygon shape file (not sort), but it seems that even sorting the polygons by area in ascending order doesn't allow a smaller polygon to be selected (using the mouse) individually.


title: How to mask clouds in ER Mapper?

text: I am using ER Mapper. But, I'm still confused how to detect and mask the clouds. I want to reduce or eliminate the cloud to perform supervised and unsupervised classification.


title: error in loading raster tiles on postgreSQL

text: I am trying to import raster data from PostGIS to GeoServer. I follow the instructions given on this link. http://docs.geoserver.org/stable/en/user/tutorials/imagemosaic- jdbc/imagemosaic-jdbc_tutorial.html I used the gdal_retile to make tiles of my GTiff imagery. I configured xml files and placed them at desire place. Then I used the DDL generation utility which create 4 SQL scripts file, then I ran the SQL scripts. But when I start importing the image tiles, it insert few tiles then abort the insertion process. Below is the screenshot. ![Below is the screenshot](http://i.stack.imgur.com/6wEoQ.png) kindly help me.


title: How can I use tiles in an ArcGIS .tpk file with Google Maps?

text: ArcGIS can export a map as server tiles as a compressed folder with the extension .tpk, as described in the ArcGIS documentation. However, it is not clear how I can get Google maps to read these files. I have previously used a ruby script to generate tiles, but ArcGIS can generate them more rapidly. However, my application currently expects to find a set of .png tiles for each zoom level, e.g. for zoom level 3:               [user@server map]$ ls *-3.png     0-0-3.png  1-2-3.png  2-4-3.png  3-6-3.png  5-0-3.png  6-2-3.png  7-4-3.png     0-1-3.png  1-3-3.png  2-5-3.png  3-7-3.png  5-1-3.png  6-3-3.png  7-5-3.png     0-2-3.png  1-4-3.png  2-6-3.png  4-0-3.png  5-2-3.png  6-4-3.png  7-6-3.png     0-3-3.png  1-5-3.png  2-7-3.png  4-1-3.png  5-3-3.png  6-5-3.png  7-7-3.png     0-4-3.png  1-6-3.png  3-0-3.png  4-2-3.png  5-4-3.png  6-6-3.png     0-5-3.png  1-7-3.png  3-1-3.png  4-3-3.png  5-5-3.png  6-7-3.png     0-6-3.png  2-0-3.png  3-2-3.png  4-4-3.png  5-6-3.png  7-0-3.png     0-7-3.png  2-1-3.png  3-3-3.png  4-5-3.png  5-7-3.png  7-1-3.png     1-0-3.png  2-2-3.png  3-4-3.png  4-6-3.png  6-0-3.png  7-2-3.png     1-1-3.png  2-3-3.png  3-5-3.png  4-7-3.png  6-1-3.png  7-3-3.png      but it is not clear to me what the analagous representation of these files is in a .tpk layer. An example of my output is map_tiles.tpk


title: How do I edit info feature attributes in GeoServer with OpenLayers

text: Can you help me? I have a very difficult problem when doing project... I have used OpenLayers and GeoServer. When I click in the map to get feature info, I receive the info for the feature and region clicked. Example: RegionName:Hanoi,Population:500... I want to know how I can change the attribute value "RegionName" from Hanoi to HaThanh and save it in Geoserver. After I load the map I wish to receive the changed information. Example: RegionName:HaThanh,Population:500... Help me! Thanks.


title: I cannot install QGIS 2.4 on my Windows 7, 64 bit machine

text: I am running Windows 7, 64 bit and cannot install QGIS 2.4 after several days of trying. I have also tried downloading 2.2 with no avail. I keep getting this error: ![NSIS error](http://i.stack.imgur.com/61CuA.jpg) The text of the error is: > NSIS Error > > Installer integrity check has failed. Common causes include incomplete > download and damaged media. Contact the installer's author to obtain a new > copy. > > More information at: http://nsis.sf.net/NSIS_Error


title: CGI mapserver with PostGIS layer: pass SQL WHERE clause in URL paramaters

text: I use a third-party software to generate a mapfile and serve them on a OpenLayers based interface from a postgis database. The postgis table, representing the movements of monks in a specific area, has the following structure (geom column trimmed and name starred for data anonymity):                mmid | monk_code |  monk_name   | year |   geom                                                 ------+-----------+--------------+------+--------------         4 | 4910102   | ***          | 2008 | 0102...2E40         1 | 4910102   | ***          | 2002 | 0102...2E40         2 | 4910102   | ***          | 2006 | 0102...2E40         5 | 4910102   | ***          | 2009 | 0102...2F40         3 | 4910102   | ***          | 2007 | 0102...2E40         7 | 4910201   | ***          | 2001 | 0102...2E40         6 | 4910201   | ***          | 2001 | 0102...2E40        25 | 4910302   | ***          | 2005 | 0102...2E40        24 | 4910302   | ***          | 2002 | 0102...3040        26 | 4910303   | ***          | 2006 | 0102...3040      The layer part of the mapfile generated by the middleware looks like this:                 LAYER         CONNECTION "PG: dbname=kc_track user=*** password=*** host=localhost port=5432"         CONNECTIONTYPE OGR         DATA "movements_search"         EXTENT -120 -37 144 51.9         LABELITEM "year"         METADATA           "mm_group"        "Group"           "ows_featureid"   "mmid"           "gml_featureid"   "mmid"           "ows_name"        "movements_search"           "mmDSTN"  "kc_track"           "ows_srs" "EPSG:4326 EPSG:900913"           "ows_title"   "movements_search"           "mmClass" "us"           "gml_include_items"   "all"           "gml_types"   "auto"         END # METADATA         NAME "movements_search"         PROCESSING "CLOSE_CONNECTION=DEFER"      And the PNG is generated through the following URL: http://mapmint.gaia.net.cias.kyoto-u.ac.jp/cgi- bin/mapserv?map=/var/www/mapmint/data/maps/search_testSearch2_movements_search.map&LAYERS=movements_search&FORMAT=image%2Fpng&SERVICE=WMS&VERSION=1.1.1&REQUEST=GetMap&STYLES=&SRS=EPSG%3A900913&BBOX=892784.49024685,-780269.18462653,19399106.279851,5809214.1488647&WIDTH=1891&HEIGHT=673 Which gives the following result: ![request result](http://i.stack.imgur.com/J3isv.png) Here is the output of `mapserv -v`               MapServer version 6.0.3 OUTPUT=GIF OUTPUT=PNG OUTPUT=JPEG OUTPUT=KML SUPPORTS=PROJ SUPPORTS=AGG SUPPORTS=CAIRO SUPPORTS=FREETYPE SUPPORTS=ICONV SUPPORTS=WMS_SERVER SUPPORTS=WMS_CLIENT SUPPORTS=WFS_SERVER SUPPORTS=WFS_CLIENT SUPPORTS=WCS_SERVER SUPPORTS=SOS_SERVER SUPPORTS=FASTCGI SUPPORTS=THREADS SUPPORTS=GEOS INPUT=POSTGIS INPUT=OGR INPUT=GDAL INPUT=SHAPEFILE      I would like to only show the lines respective to a specific monk (e.g. monk_code = '4910102'). Is there a way to include a filter, such a SQL WHERE clause, in the request URL? I have full access to the database but the mapfile is generated automatically by the middleware. If there is no other solution, I can dive into the middleware source myself and change the code but I would prefer a solution that could use only the request part. Edit: The following WFS request returns the correct points: http://mapmint.gaia.net.cias.kyoto-u.ac.jp/cgi- bin/mapserv?map=/var/www/mapmint/data/maps/search_testSearch2_movements_search.map&request=GetFeature&service=WFS&version=1.0.0&typename=movements_search&filter=%3CFilter%3E%3CPropertyIsLike%20wildcard=%27*%27%20singleChar=%27.%27%20escape=%27!%27%3E%3CPropertyName%3Emonk_code%3C/PropertyName%3E%3CLiteral%3E*4910201*%3C/Literal%3E%3C/PropertyIsLike%3E%3C/Filter%3E But the following WMS request with the same filter still returns all the lines: http://mapmint.gaia.net.cias.kyoto-u.ac.jp/cgi- bin/mapserv?map=/var/www/mapmint/data/maps/search_testSearch2_movements_search.map&request=GetMap&LAYERS=movements_search&service=WMS&version=1.1.1&typename=movements_search&filter=%3CFilter%3E%3CPropertyIsLike%20wildcard=%27*%27%20singleChar=%27.%27%20escape=%27!%27%3E%3CPropertyName%3Emonk_code%3C/PropertyName%3E%3CLiteral%3E*4910201*%3C/Literal%3E%3C/PropertyIsLike%3E%3C/Filter%3E&SRS=EPSG%3A900913&WIDTH=1891&HEIGHT=673&BBOX=892784.49024685,-780269.18462653,19399106.279851,5809214.1488647&FORMAT=image%2Fpng


title: Geoprocessing (error 99998)

text: > **Possible Duplicate:**   >  ERROR 999998: Unexpected Error 1.Service Pack - have no Idea. I received the Arc10 Software with my study material.(I just installed it)   1. Shp files (want to clip streets according to a boundary.   2. It does not really matter what tool I use - I get that same error (Merge/ clip/ buffer)


title: Automatically create Text/HTML annotations in QGIS on a Point Layer

text: How would one proceed to generate an HTML annotation or a Text annotation for all points in a shapefile? They are all points. What I want to do is load up a shapefile containing coordinates and values + absolute paths to Photographs and display them in an HTML Field. It would allow all the main information to be available directly at a glance. Since there are over 50 locations, I don't think it is viable to be done by hand.


title: Get km route between two points with lat / lng

text: I'm doing the following query to get the distance between two points with pgrouting:               SELECT             sum(cost) AS cost     FROM pgr_dijkstra(                     'SELECT id, source, target, km as cost FROM at_2po_4pgr',                     50347, 32988, false, false             );      It turns out that the source and destination comes from nominatim geocode coordinates in, and I do not know how to use the coordinates as parameter in sql. Can anyone help me?


title: Adding XY Points to ArcMap 10.2

text: I am trying to import an Excel table with around 10,000 XY points into ArcMap 10.2. I am using State Plane NAD 1983 Alabama West (US Feet) as the Coordinate System in my map data frame. I then add the excel table to my Table of Contents in the map then select Display XY Data. Next, I choose the same coordinate system as my data frame then select OK. Once it creates the Events layer for my points they are not in the correct location. Should I switch everything in my map to WGS84 or is there a solution for my points to be added onto my current map with the state plane coordinate system? Any suggestions would be beneficial! Thanks. ![XY Points in Excel](http://i.stack.imgur.com/GGe8F.jpg) ![enter image description here](http://i.stack.imgur.com/8G7U3.jpg)


title: "TypeError: Unable to draw graphic (geometry:null, symbol:null)" after adding tiled and feature layers

text: I'm seeing the following error when attempting to add a tiled and feature layer to the same map when using the ArcGIS Server JS API 3.1: > TypeError: Unable to draw graphic (geometry:null, symbol:null): _c is > undefined > http://serverapi.arcgisonline.com/jsapi/arcgis/3.1compact/js/dojo/dojo/DeferredList.js > Line 8 After defining the map and adding a basemap, I'm iterating through a list of JSON layer definitions, building each layer based on its definition, then adding it to the map. The pseudo code is:               dojo.connect(map, "onLoad", function {        for layerDef in layerDefs:           if layerDef.type == tiled:              layer = buildTiledLayer           if layerDef.type == feature:              layer = buildFeatureLayer           map.addLayer(layer)     }      If I add just the tiled layer, or just the feature layer, there are no problems - it's only when I add both together that I see the error message. Once the error message has shown, an infoWindow is not shown when clicking on the feature layer. This forum post lead me to suspect the problem may be due to the timing of adding the layers - during debugging I found that if I remove the `map.addLayer` line, and instead push the layers into an array, I can manually enter `map.addLayers(layerArray)` in Firebug, once the map has finished loading. In this case, both layers add correctly. Adding a FeatureLayer ARCGIS JavaScript API isn't relevant as I know the feature layer definition is correct - the feature layer adds without error if it's added without the tiled layer. Similarly, I believe that the layers are correctly defined since each is loaded correctly without the other. Does anyone have any advice as to why this error might be occurring, and how to resolve the error? Thanks


title: Compability problems with OpenGeo Suite and already installed GeoServer and others

text: Sorry I couldn´t find the answer anywhere else. I already have GeoServer 2.3.2, Postgres 9.2, OpenLayers 2.13 installed but in addition I need GeoWebCache and probably GeoExt. Is it wice to download now the whole OpenGeo Suite even if the versions are different or should I try to download GeoWebCache and GeoExt separately? Thanks for help! Laura


title: Openlayers WFS Vector Layer with Popups: how to setCenter on point when you click it?

text:
Working code I have so far is below. Lots of airports on the map, you click an airport icon, the popup appears. No problems there. **Goal: to click a popup, and based on the lat/long data of that particular point, the map will center on that clicked point.**               wfs.events.on({          featureselected: function(event) {                 var feature = event.feature;                 feature.popup = new OpenLayers.Popup.FramedCloud("box",                         feature.geometry.getBounds().getCenterLonLat(),                         null,                         '
Airport: '+feature.attributes.Name +'
 ICAO: '+feature.attributes.ICAO+'
 IATA: '+feature.attributes.IATA+'
',                         null,                         true                 );                 while( map.popups.length ) {                         map.removePopup( map.popups[0] );                 }                 map.addPopup(feature.popup);                 }         });      Not sure how to go about this after browsing through other questions/looking at documentation. Also: `feature.geometry.getBounds().getCenterLonLat()` What exactly is the point of this? Based on the docs it's for `{OpenLayers.LonLat}` but it's not like it is reading the data from there? Any input appreciated.


title: Grass v.clean.snap with Polygons

text: I am doing some topology testing using Grass in QGIS. I am hoping to snap nodes to eliminate gaps between polygons. When I use the v.clean.snap function in Grass it does not work on polygons. When I test it on lines it completes the tasks as I was hoping. Does anyone know of a way to run v.clean.snap or a similar function in Grass on Polygons? Attached are two screen grabs showing the errors as polygons, and the fix as a line.![Grass v.clean.snap with polygons](http://i.stack.imgur.com/a2eYv.jpg)


title: Load information about a feature

text: I must mention first of all that I am very new to this area... I am developing a GIS application. I created a simple map using geoserver with openlayers, it is all working through apache tomcat - and is connected to postgis/postgresql database. When I click on a point (the points represent weather stations - and stored in database) it is set to display the tablewith all the data like: fid, station_name etc. What I want to do is - to store station_name(I clicked on) in a variable and use for my needs. In this case I was trying to break response.responseText in parts, and get it from there but it never worked... response.responseText - returns all the information and creates the tables. I need to only get ONE value from database. Any ideas!? Any help greatly appriciated!!! Here is the code: In this part it is set to read a click and make a response:               map.events.register('click', map, function (e) {                          document.getElementById('nodelist').innerHTML = "Loading... please wait...";                     var params = {                         REQUEST: "GetFeatureInfo",                         EXCEPTIONS: "application/vnd.ogc.se_xml",                         BBOX: map.getExtent().toBBOX(),                         SERVICE: "WMS",                         INFO_FORMAT: 'text/html',                         QUERY_LAYERS: map.layers[0].params.LAYERS,                         FEATURE_COUNT: 50,                         Layers: 'map_with_stations',                         WIDTH: map.size.w,                         HEIGHT: map.size.h,                         format: format,                         styles: map.layers[0].params.STYLES,                         srs: map.layers[0].params.SRS};                          // handle the wms 1.3 vs wms 1.1 madness                     if(map.layers[0].params.VERSION == "1.3.0") {                         params.version = "1.3.0";                         params.j = parseInt(e.xy.x);                         params.i = parseInt(e.xy.y);                     } else {                         params.version = "1.1.1";                         params.x = parseInt(e.xy.x);                         params.y = parseInt(e.xy.y);                     }                          // merge filters                     if(map.layers[0].params.CQL_FILTER != null) {                         params.cql_filter = map.layers[0].params.CQL_FILTER;                     }                      if(map.layers[0].params.FILTER != null) {                         params.filter = map.layers[0].params.FILTER;                     }                     if(map.layers[0].params.FEATUREID) {                         params.featureid = map.layers[0].params.FEATUREID;                     }                     OpenLayers.loadURL("http://localhost:82/geoserver/wms", params, this, setHTML, setHTML);                     OpenLayers.Event.stop(e);                 });             }                  // sets the HTML provided into the nodelist element             function setHTML(response){                      document.getElementById('nodelist').innerHTML = response.responseText;                            };


title: Why do my map labels disappear when publishing as a service on ArcGIS server?

text: I am developing a map service, rather basic, with some polygons and labels. Surprisingly, when I deploy this service on the production server, the labels disappear (see here). I have a second local server - for this one, no problem: the labels are well and truly there. Do you have an idea why my labels vanish that way? Is there some special tuning I should apply on the production server to have them back? Many thanks for your help. **Precision:** It happens with labels with field information. If I use a static text as a labels, the labels are there. The server seems unable to include object fields info within labels. **Additional info:** Apparently, this is a well-known problem in ArcGIS server 10.1. See this thread on ArcGIS forum. No solution is known yet...


title: How to obtain map service url from within Server Object Extension (.NET)

text: I'm having trouble figuring out how to obtain the full url to the map service a Server Object Extension is running on top of. I bascially want:               IMapServer4 mapServer = (IMapServer4)serverObjectHelper.ServerObject;     string fullUrl = mapServer.url       Which would yield: http://sampleserver6.arcgisonline.com/arcgis/rest/services/Census/MapServer Any Ideas? Thanks for any suggestions.


title: Help with Python Update Cursor syntax

text: I'm having trouble running a loop to update the values in a particular field in attribute table by reading from a csv file. Any input is greatly appreciated. The following is the sample code:               import arcpy          featureclass = "C:\\xyz"     table1 = "C:\\pqr.csv"     alpha= list(csv.reader(open(table1,'r')))          crows=arcpy.UpdateCursor(featureclass) # guess here it should be like cursor( , )          for crow in crows:         for row in alpha:             if crow.STNAME==row[3] and crow.ZIP==row[6]:                 if int(crow.STNUM)>=int(row[13]) and int(crow.STNUM)<=int(row[14]):                     if row[4]==crow.STTYPE and row[15]=="B":                        # fill a field in attribute table with a particular row[] value from csv                     else:                         if int(crow.STNUM)%2==0 and row[4]==frow.STTYPE and row[15]=="E":                       # fill the field in attribute table with a particular row[]value from csv


title: How to create a 4096x4096 png from a 1440x720 data set using GDAL?

text: I'm having a bear of a time trying to do something that seems fairly straightforward. I'm trying to take a netcdf file from ftp://eclipse.ncdc.noaa.gov/pub/OI-daily-v2/NetCDF/2011/AVHRR-AMSR/ and turn it into a 4096x4096 png so that I can tile it up for use w/google maps. I've been using the following commands but it always seems to create blank areas on the east and west borders:               /usr/local/bin/gdal_translate -a_srs EPSG:4326 -outsize 4096 4096 -a_ullr 0 90 360 -90 NETCDF:"input.nc":sst -of netCDF output.nc     /usr/local/bin/gdaldem color-relief output.nc -b 1 -alpha colorscale.txt -of PNG output.png     /usr/bin/convert output.png -transparent #FFFFFF output.png      Am I just approaching this the wrong way? Am I using the wrong tools? FWIW - gdal_translate is v1.8.0 Thanks in advance. **UPDATE:** "Answer" below from Dija - colleague of mine - is actually a response to the comments on the question itself. **UPDATE** Accidentally added an answer instead of a comment. Moving it here. I thought I'd weigh in here. I've been working on this issue for a little while now. I've tried many iterations of the process. I scaled these images to the 4096x4096 as it actually creates a better image than if not. Meaning, if I leave the image in it's native size, 1440x720 I get a larger gap at the prime meridian (for this data set) and at the 180th meridian for other data sets. I scaled the output image to be a multiple of 256 (to fit better with google) and get various results. If I scale it to 1536x768 the gap shows up to the right of the 0 meridian. If I scale to 1280x640 it's to the left. If I leave it alone, it's pretty much splits it. I have taken the base image that gdaldem spits out, cut it in half and matched the edges up...they line up perfectly. This leads me to think that there's something weird going on with gdal2tiles.py I wrote my own tiler that I'm getting ready to turn loose on our images to see if the gap remains. I'll put together some screen shots to show what I am talking about with how the gap as well as report on what my tiler does. Thanks for any input you might have. It's appreciated. **UPDATE 2** I ran the maptiler I wrote against the original raster and it works fine with out a gap. gdal2tiles.py seems to be the problem. It introduces a transparent edge to the last tile that touches the 180th meridian.


title: How do use the Shapefile Converter in Oracle 12c to load a Shapefile?

text: **Update:**   I resolved the service name issue. At some point I put in the right combo of parameters (or maybe I omitted some) and this message was displayed:   `The following key/value is missing: -s db_sid OR -sn db_service_name`   The `-sn` option isn't in the docs as far as I can tell. In my case, I replace the `-s GDSRV` with `-sn pdborcl`.               Now I'm getting all sorts of errors:       Connecting to Oracle Database using...       localhost, 1521, null pdborcl, GISD, ********, shapes, tl_2013_us_state, null, 0       Dropping old table...          Table not previously created.       Creating new table...       java.sql.SQLSyntaxErrorException: ORA-00902: invalid datatype            java.sql.SQLSyntaxErrorException: ORA-00942: table or view does not exist          Record #1 not converted.       java.sql.SQLSyntaxErrorException: ORA-00942: table or view does not exist          Record #2 not converted.       java.sql.SQLSyntaxErrorException: ORA-00942: table or view does not exist          ...            Record #56 not converted.     56 record(s) not converted.     0 record(s) converted.     Done.      **Original Question:** I am attempting to use the shapefile converter in Oracle12c and I'm not having much success and I think it is due to the container database (CDB) vs pluggable database (PDB). I need to somehow specify my PDB: `pdborcl` If I provide the SID, I get an invalid username/password....as expected. `C:\temp>java oracle.spatial.util.SampleShapefileToJGeomFeature -h localhost -p 1521 -s GDSRV -u GISD -d gisd -t shapes -f tl_2013_us_state -r 0 -g geom`               host: localhost     port: 1521     sid: GDSRV     db_username: GISD     db_password: ********     db_tablename: shapes     shapefile_name: tl_2013_us_state     SRID: 0     db_geometry_column: geom     Connecting to Oracle Database using...     localhost, 1521, GDSRV null, GISD, ********, shapes, tl_2013_us_state, null, 0     java.sql.SQLException: ORA-01017: invalid username/password; logon denied                  at oracle.jdbc.driver.T4CTTIoer.processError(T4CTTIoer.java:450)             at oracle.jdbc.driver.T4CTTIoer.processError(T4CTTIoer.java:392)             at oracle.jdbc.driver.T4CTTIoer.processError(T4CTTIoer.java:385)             at oracle.jdbc.driver.T4CTTIfun.processError(T4CTTIfun.java:938)             at oracle.jdbc.driver.T4CTTIoauthenticate.processError(T4CTTIoauthenticate.java:480)             at oracle.jdbc.driver.T4CTTIfun.receive(T4CTTIfun.java:655)             at oracle.jdbc.driver.T4CTTIfun.doRPC(T4CTTIfun.java:249)             at oracle.jdbc.driver.T4CTTIoauthenticate.doOAUTH(T4CTTIoauthenticate.java:416)             at oracle.jdbc.driver.T4CTTIoauthenticate.doOAUTH(T4CTTIoauthenticate.java:825)             at oracle.jdbc.driver.T4CConnection.logon(T4CConnection.java:596)             at oracle.jdbc.driver.PhysicalConnection.<init>(PhysicalConnection.java:715)             at oracle.jdbc.driver.T4CConnection.<init>(T4CConnection.java:385)             at oracle.jdbc.driver.T4CDriverExtension.getConnection(T4CDriverExtension.java:30)             at oracle.jdbc.driver.OracleDriver.connect(OracleDriver.java:564)             at java.sql.DriverManager.getConnection(Unknown Source)             at java.sql.DriverManager.getConnection(Unknown Source)             at oracle.spatial.util.SampleShapefileToJGeomFeature.main(SampleShapefileToJGeomFeature.java:341)     Exception in thread "main" java.lang.NullPointerException             at oracle.spatial.util.SampleShapefileToJGeomFeature.main(SampleShapefileToJGeomFeature.java:345)      From SQLPLUS, I would use the PDB name instead to connect and I'm in. How do you connect to the PDB when calling oracle.spatial.util.SampleShapefileToJGeomFeature? `java oracle.spatial.util.SampleShapefileToJGeomFeature -h localhost -p 1521 -s pdborcl -u GISD -d gisd -t shapes -f tl_2013_us_state -r 0 -g geom`               host: localhost     port: 1521     sid: pdborcl     db_username: GISD     db_password: ********     db_tablename: shapes     shapefile_name: tl_2013_us_state     SRID: 0     db_geometry_column: geom     Connecting to Oracle Database using...     localhost, 1521, pdborcl null, GISD, ********, shapes, tl_2013_us_state, null, 0     java.sql.SQLException: Listener refused the connection with the following error:     ORA-12505, TNS:listener does not currently know of SID given in connect descriptor                  at oracle.jdbc.driver.T4CConnection.logon(T4CConnection.java:673)             at oracle.jdbc.driver.PhysicalConnection.<init>(PhysicalConnection.java:715)             at oracle.jdbc.driver.T4CConnection.<init>(T4CConnection.java:385)             at oracle.jdbc.driver.T4CDriverExtension.getConnection(T4CDriverExtension.java:30)             at oracle.jdbc.driver.OracleDriver.connect(OracleDriver.java:564)             at java.sql.DriverManager.getConnection(Unknown Source)             at java.sql.DriverManager.getConnection(Unknown Source)             at oracle.spatial.util.SampleShapefileToJGeomFeature.main(SampleShapefileToJGeomFeature.java:341)     Caused by: oracle.net.ns.NetException: Listener refused the connection with the following error:     ORA-12505, TNS:listener does not currently know of SID given in connect descriptor                  at oracle.net.ns.NSProtocolStream.negotiateConnection(NSProtocolStream.java:272)             at oracle.net.ns.NSProtocol.connect(NSProtocol.java:263)             at oracle.jdbc.driver.T4CConnection.connect(T4CConnection.java:1360)             at oracle.jdbc.driver.T4CConnection.logon(T4CConnection.java:486)      Here's the docs on the 12c converter: http://docs.oracle.com/cd/E16655_01/appdev.121/e17896/sdo_shapefile_converter.htm#SPATL1427


title: pgShapeLoader: how to load a shapefile in a different schema than 'public'

text: I'm trying to load a shapefile in POSTGIS using pgShapeLoader. It works but I can load table only in the "public" schema and i don't know how to configure, the tool to upload my shapefile in another schema. Does it possibile? Any suggestion for workaround? Anyway I can use shp2pgsql and then psql from command line in POSTGIS .... Thank you very much in advance! Cesare


title: model only runs open "modelbuider"

text: I am using ArcGIS 10.2.2 for Desktop. I have a model made in ModelBuilder which currently does:   1. input XY coordinates excel table    2. data display XY event    3. converting a format event shape    4. add field "name1"    5. add field "name2"  So far so good but my problems are:   1. the model only works so edition    2. I cannot add tool to make a selection of my new shapefile ![enter image description here](http://i.stack.imgur.com/42wh5.png)


title: Change in behaviour of layer extents

text: I have a class that inherits a vector memory layer, it all works fine a previous version of QGIS, but on 2.2 it no longer works. It either does not know that the layer is WGS84 or getting the extents is somehow failing.               class DotTrail(QgsVectorLayer):         def __init__(self,name, colour):             QgsVectorLayer.__init__(self, "multipoint?crs=WGS84", name , "memory")             self.pr =self.dataProvider()         self.featId = -1          self.points = []                       # set the colour             rendererV2 = self.rendererV2()             layerSymbol = QgsSymbolV2.defaultSymbol(self.geometryType())             layerSymbol.setColor(QtGui.QColor(colour))             layerSymbol.setAlpha(1)             rendererV2.setSymbol(layerSymbol)             QgsMapLayerRegistry.instance().addMapLayers([self])              def addPointUpdate(self, p):             self.points.append(p)             self.update()              def addPoint(self, p):             self.points.append(p)              def clear(self):             self.points = []             #self.update()              def update(self):             print 'Update'         if self.featId == -1:             feat = QgsFeature()             feat.setGeometry(QgsGeometry.fromMultiPoint(self.points))             res, fid = self.pr.addFeatures( [feat] )             self.featId = fid[0].id()             for feat in self.pr.getFeatures():                 print feat.geometry().asMultiPoint()                 else:                 self.pr.changeGeometryValues({self.featId: QgsGeometry.fromMultiPoint(self.points)})                 feat = QgsFeature()                  self.updateExtents()             self.updateFields()         print 'Dot layer extents'         print self.extent().toRectF().getCoords()      It adds the following points Update [(151.193,-33.8895), (151.193,-33.8895), (151.195,-33.8877), (151.195,-33.8913), (151.191,-33.8913), (151.191,-33.8877), (151.195,-33.8877), (151.195,-33.8877)] But outputs the following extent Dot layer extents (0.0, -33.89132551272559, 151.1951840971511, 0.0) I can't see anything wrong with it.


title: Use of postgis with hibernate

text: We have an existing postgis database. Our application is developed in Java, but queries database through a jni dll, which appears to be not very effective for our purpose. I've managed to migrate to jdbc, and it's working well, but the next step would be to rely on a framework such as hibernate-spatial, which i don't know at all. How mature is it? I've seen some people here using it, what do you think of it? It looks that it supports only cartesian coordinates and it may be a concern for us for distance computation.


title: Modify raster layer display and symbology properties programmatically with arcpy or ArcObjects

text: I have over 1500 NITF files that I need to be able to modify the layer display and symbology of programmatically from within ArcMap. I need to set the properties outlined in red below: ![enter image description here](http://i.stack.imgur.com/4WZFa.png) ![enter image description here](http://i.stack.imgur.com/AMmhj.png) I've looked through the arcpy docs and it doesn't look like these properties are exposed - honestly, I didn't expect them to be. I have a feeling I can get to them through ArcObjects. Anyone had any experience/luck modifying these properties programmatically?


title: A global, grid-like projection for creating heatmaps

text: I'm putting together an application where I need to create a vector grid which will be used to store and display a heatmap. It has the following requirements:   * Can cover the entire planet.    * The vast majority of grid squares won't have values.   * I don't want to have to store the grid itself; I'd like to calculate it on the fly.   * The scale of the data used with the grid could vary massively.   * I anticipate wanting grid squares anything from 1km across to 100km across. (I'm aware how many this will be (~510million for 1km, ~51,000 for 100km)).   * Values will be accumulated/aggregated for each grid square.   * Ideally I'd be able to easily use smaller grid cells to calculate values for larger ones rather than store the larger grid cell values.   * I'll be using OpenLayers to drape it over OpenStreetMap.   * I'll be storing it in SpatiaLite or SQLite, so preferably supported by those natively (i.e. for SpatiaLite = a supported CRS; or for SQLite = a pure numbers based system). So my question is: What projection should I use for this grid? Also - is there a good way to designing this? Does anyone know of a good potential solution to this problem or have solved similar before? Or can point me in a useful direction. **Edit** Use-case - basically I'm aggregating up bounding boxes of various different shapes and sizes. They can be anything from a few hectares to thousands of square km in size. They may also be in different projections. Below is a bespoke version of the sort of thing I'm aiming to accomplish automatically on a larger scale. ![enter image description here](http://i.stack.imgur.com/YEvt9.jpg) Many thanks.


title: QGIS projects - collaboration on different platform

text: on my computer, I have both - Ubuntu and Windows - operating systems. Since I work on both of them and QGIS save path to particular layers using absolute path, it's very inconvenient to store projects on local disk. Moreover, I wish to access QGIS projects and layers from other computers. Is there any possibility to store them on the Internet? What will be the best, free of charge, solution? I thought about Dropbox, but I'm not sure if there's more convenient way.


title: How to change the name of "Base layer" and "Overlays" in LayerSwitcher?

text: I would like to change the name of "Base layer" and "Overlays" in my LayerSwitcher div to something else. I've noticed that the div including the texts only has a class and no ID which means it's little bit tricky to access. I managed to change "Base layer" name with               document.getElementById("layerswitcher").firstChild.firstChild.innerHTML = "new name";      but I don't like this solution. Best would be if I could define the names directly with OpenLayers functions. Is there one? I read something about using language translations but it's not what I'm after.


title: How can I use a neural network in GIS?

text: I'd like to learn how to use a neural network in GIS. I surfed in the internet to find good resources or a tutorial about this but I couldn't find any helpful resources. Are there any good references or tutorials about using neural networks in GIS? I'd like to use an artificial neural network (ANN) in mapping hazard zones, for example flood zones or landslides. I use ArcGIS 10 to process the layers but don't know how can use an ANN.


title: How to calculate the track length of GPX files with PostGIS?

text: I am trying to calculate the length of a GPX track. The original GPS data looks like this:               <?xml version='1.0' encoding='UTF-8' standalone='yes' ?>     <gpx creator="nl.sogeti.android.gpstracker"        version="1.1" xmlns="http://www.topografix.com/GPX/1/1"        xmlns:gpx10="http://www.topografix.com/GPX/1/0"        xmlns:ogt10="http://gpstracker.android.sogeti.nl/GPX/1/0"        xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"        xsi:schemaLocation="http://www.topografix.com/GPX/1/1        http://www.topografix.com/gpx/1/1/gpx.xsd">         <metadata>             <time>2014-04-25T19:20:39Z</time>         </metadata>         <trk>             <name>Track 2014-04-25 19:20</name>             <trkseg/>             <trkseg>                 <trkpt lat="52.51025036" lon="13.3907098">                     <ele>93.30000305175781</ele>                     <time>2014-04-25T19:59:24Z</time>                     <extensions>                         <gpx10:speed>2.75</gpx10:speed>                         <ogt10:accuracy>7.0</ogt10:accuracy>                         <gpx10:course>236.5</gpx10:course>                     </extensions>                 </trkpt>                 <trkpt lat="52.51021522" lon="13.390611">                     <ele>94.0</ele>                     <time>2014-04-25T19:59:26Z</time>                     <extensions>                         <gpx10:speed>2.5</gpx10:speed>                         <ogt10:accuracy>7.0</ogt10:accuracy>                         <gpx10:course>235.1999969482422</gpx10:course>                     </extensions>                 </trkpt>                 <!-- And many more -->             </trkseg>         </trk>     </gpx>      I imported the data into PostgreSQL/PostGIS. The track points have been imported into the table `track_points`:               gpxdata=# \d track_points               Table "public.track_points"            Column       |           Type           |     --------------------+--------------------------+      ogc_fid            | integer                  |       wkb_geometry       | geometry(Point,4326)     |       track_fid          | integer                  |       track_seg_id       | integer                  |       track_seg_point_id | integer                  |       ele                | double precision         |       time               | timestamp with time zone |       magvar             | double precision         |       geoidheight        | double precision         |       name               | character varying        |       cmt                | character varying        |       desc               | character varying        |       src                | character varying        |       link1_href         | character varying        |       link1_text         | character varying        |       link1_type         | character varying        |       link2_href         | character varying        |       link2_text         | character varying        |       link2_type         | character varying        |       sym                | character varying        |       type               | character varying        |       fix                | character varying        |       sat                | integer                  |       hdop               | double precision         |       vdop               | double precision         |       pdop               | double precision         |       ageofdgpsdata      | double precision         |       dgpsid             | integer                  |       gpx10_speed        | double precision         |       ogt10_accuracy     | double precision         |       gpx10_course       | double precision         |       Now I use the following query which I found here.               SELECT gps.ogc_fid, tracks.name,      ST_Length_Spheroid(ST_MakeLine(gps.wkb_geometry),      'SPHEROID["WGS 84",6378137,298.257223563]') AS track_len         FROM (             SELECT ogc_fid, time, wkb_geometry              FROM track_points              ORDER BY ogc_fid, time         )          AS gps, tracks         WHERE gps.ogc_fid = tracks.ogc_fid         GROUP BY gps.ogc_fid, tracks.name         ORDER BY gps.ogc_fid;      The result of the query is:                  | ogc_fid |          name          |    track_len     |        | integer |    character varying   | double precision |     ----------------------------------------------------------      1 |    1    | Track 2014-04-25 19:20 |        0         |      * * * The following query posted by Brad ..               SELECT ogc_fid, time, ST_AsEWKT(wkb_geometry) FROM track_points LIMIT 5;      ... produces this result:               ogc_fid |           time           | st_asewkt     integer | timestamp with time zone | text     ------------------------------------------------------------------------------        1    |  2014-04-25 19:31:13+02  | SRID=4326;POINT(13.38866662 52.54784071)         2    |  2014-04-25 19:31:14+02  | SRID=4326;POINT(13.3886988 52.54778812)        3    |  2014-04-25 19:31:15+02  | SRID=4326;POINT(13.38873776 52.54773285)        4    |  2014-04-25 19:31:16+02  | SRID=4326;POINT(13.38877592 52.5476608)        5    |  2014-04-25 19:31:17+02  | SRID=4326;POINT(13.38881506 52.54760139)


title: Regression analysis - creating dataset

text: ![enter image description here](http://i.stack.imgur.com/895yf.png) I have a vector-based dataset of soils/parent material. The statistics have already been done, with coefficients for each category of parent material. I am not sure how to use this information to create a soils raster layer for analysis within Raster Calculator. For example, to create the model result raster, I'm thinking: = 0.158544748 + (-0.422151664 * LAT) + (01.184392010 * LONG) + (0.001245377 * ELEV) ...and so on for all the variables. BUT, how do I create the raster for the soils/parent material information, and how would that look in the above equation for the raster calculator? (they are labeled PARMAT AL, PARMAT GR, PARTMAT GRHW... in the table) Thanks! EDIT- I'm using ArcGIS 10.1, have Spatial Analyst and many other extensions (not sure what's relevant). Also, I was previously advised to 'create a "ParamEst" column for each of the PARMAT (soils) values, to be populated with the parameter estimates, then convert to raster. Each PARMAT type would then be a parameter estimate and can be added as a standalone layer.' Therefore, would each PARMAT type be exported out as its own raster, with its parameter estimate, and simply added to the model? So instead of + (-1.121263096 * PARMATAL), would it simply be + (PARMATAL) in the equation, because the derived raster is already a parameter estimate? I hope I'm making sense - I only have a basic grasp of statistics, so I'm not sure if I'm asking things properly. ***I'm unable to reach that person again, which is why I haven't had further clarification.


title: postgis -- return records sorted by size of st_intersection

text: I have two sets of spatial records in two different tables. I'm in rails so I'll show my code in that way seeing that its source: one set is CensusPlace records, which are US Census Tigerline shapefiles. Note that Census Places are typically towns and cities (with some sub- municipalities) that can exist in multiple counties. the other set of records are County records. It would be very handy if I could run a postgis command that would return my st_intersects results for a given city like so:               CensusPlace.where(:state => "SC").each do |cp|         puts "\n---------------census place: #{cp.name}"         County.where{st_intersects(:county_shape, cp.census_place_shape)}.each_with_index do |c, i|           puts "#{c.name} county"           puts "\toverlap ranking: #{i+1}"              end      end      such that the actual returning of results would be ordered by the the overlap size. It would be natural for the results to actual return a new "virtual column" or what have you with a 0-1 value of that overlap percentage as well. Thanks


title: How to set symbology based on an attribute of a 1-N table

text: Now, it is possible to edit geographic features and add elements in child tables using 1-N relationships. I want to construct a stacked - symbology which allows me to represent the multiple attributes. For example if I have a line feature "linear park" and my linear park has a linked table (1-N) called "users", i want a stacked or multi-level symbology (with some offsets) that allows me to depict all the "users" of that section of park. Is it possible in QGIS? I know it was possible in ArcGIS using relates. I appreciate your time!!


title: How to indicate a range of zoom on the map?

text: Basically what I need and I want to know is how you can do to tell the map that respects me a zoom range to display the data, ie, I want to see the map from afar without information, and as the zooming you can go see the information. Can anyone help me with this?


title: How to convert .shp into .gpx using QGIS?

text: I'm newbie using QGIS. I would like to convert .shp into .gpx. So, using Qgis, I open my shapefile and I tried to "save as" gpx but it gives an OGR error: `creation of field ObjectId failed (OGR error: Field of name 'ObjectId' is not supported in GPX schema. Use GPX_USE_EXTENSIONS creation option to allow use of the <extensions> element.)` How can I fix this error? Thanks.


title: QGIS - Save as GPX

text: I would like to convert .shp into .gpx. I tried to "save as" gpx but it gives an OGR error: creation of field ObjectId failed (OGR error: Field of name 'Name' is not supported in GPX schema. Use GPX_USE_EXTENSIONS creation option to allow use of the element.) GPX_USE_EXTENSIONS=Y does not work. Checking 'Skip Attribute Creation' does not work.


title: How to open QGIS using VBA?

text: Using VBA, I want to open QGIS (the actual GUI) and execute a simple script through the Python console or the editor.


title: crust index in Erdas Imagine

text: hey people, i have the trouble with CI index in Erdas, when i am useing the image from landsat 5 tm with atmospheric correction it gives me the whole white picture in the results, the formula CI=1-((red-blue)/(red+blue)) also i tryed to change R band with NIR, but......does not work at all when i tryed to insert in the model the bands seperately the reslusts are the same, just white pic, can anyone help with this type of issue?


title: What's a good library for displaying many small maps on a web page

text: I would like to create a visualization on the web similar to the following: http://namemapper.babynamewizard.com/namemapper/ Are there any good libraries for doing this? Is OpenLayers light-weight enough to handle many maps at once?


title: TIGER/Line processing

text: Does anyone know of scripts for merging TIGER/Line shp files? It is for extracting the road features only, but for the whole of USA. They are provided as shp files for one county at a time here: http://www.census.gov/geo/www/tiger/ If anyone sells the data, a link is also appreciated.


title: Comparing two featureclasses using C#

text: I'm sorry, I know this is a pretty simple question but am having trouble getting this to work. Basically, I have 2 linear featureclasses and they both have the same fields and spatial reference. FC1 has 675 records and FC2 only has 505 records. So how can I check what Prop_ID(s) are missing from FC2? I have some code that will copy the missing record from FC1 to FC2 but must first determine the missing Prop_ID. many thanks


title: resize image resolution from 1000 m to 250 m

text: I'm beginner of ARCGIS and RS. I have many problems with the procedure. One of them is   I have NDVI images and Soil map. Both are at the different resolution: NDVI is 250 m. soil map is 1000 m. So I wanna ask how can I change the resolution from 1000 m to 250 m. Thanks a lot


title: Qgis 1.8 fails to install

text: QGIS 1.8 will not start on Windows 7 Basically a similar issue as the above post. Version 1.8 fails to install with the standalone installer; with OSGeo4w, while it doesn't run into trouble during the installation, ultimately fails as well (same error message as below when opening the app). Standalone installer has the following problems: At the end of the installation, it reaches "qgis-common.bat" execution and stalls; meanwhile, Task Manager shows "find string" GREP utility running in the background and repeatedly crashing every minute or so. Ending the process results it immediately starting back up. If the installation is forced to quit or completed, trying to open Qgis results in the error message that **qwt5.dll** is missing. QGREP crash signature:               Problem Event Name:APPCRASH     Application Name:   findstr.exe     Application Version:    6.1.7601.17514     Application Timestamp:  4ce78c0a     Fault Module Name:  findstr.exe     Fault Module Version:   6.1.7601.17514     Fault Module Timestamp: 4ce78c0a     Exception Code: c0000005     Exception Offset:   00003c76     OS Version: 6.1.7601.2.1.0.768.3     Locale ID:  1033     Additional Information 1:   0a9e     Additional Information 2:   0a9e372d3b4ad19135b953a78882e789     Additional Information 3:   0a9e     Additional Information 4:   0a9e372d3b4ad19135b953a78882e789


title: How to capture data for mapping

text: Well it might be an newbie question but I can't find it In my project I need the paths in a campus but these paths are inexistents so I'm going to walk around and mapping (I don't know if it's the correct term) the points to build the correspondent paths, so my question is if exists an easy way to do it or how do you think I can accomplished the work, please write all yours thoughts.. Thanks in advance..


title: QGIS Newbie Rendering Time

text: I am new to QGIS and have just installed 1.8.0 on a Gateway laptop that is about 2 years old. It is running Windows 7, has an Intel Core i3 2.13 GHz processor, 4 GB of RAM. I've just brought in my first layers (235 MB of data in 3 shapefiles) to my first project, and it takes forever to render. By forever I mean more than one minute for each task: zoom, move, rectangle select, changing the order of layers. The data is PLSS data from the state of Oregon as provided through the federal government through BLM. Although I downloaded only the Oregon data, when I right-click on the twnshp layer onto attribute table, I do see records relating to CA and ID as well as Oregon. Will I not be able to run QGIS on this computer? Thanks so much, Yeye


title: Create a dbf table by copying the structure of another dbf table

text: I have a dbf table and I wish to automatically create another dbf table with a different name but with exactly the same structure i.e.fields as the existing table. Is there any ready code in VBA and ArcObjects for that? Thanks


title: Categorized style won't display in QGIS2.2

text: I'm a qgis novice trying to put together a building characteristics map. While 44 out of 46 of my layers have their data categorized and displayed correctly, two layers will not display different colours. Instead I am getting the same colour all over my map, which is the colour assigned to an extra unvalued category that appears at the bottom of my list, as shown in the screenshot. ![Styles](http://i.stack.imgur.com/KcpTH.jpg) Any ideas? Thanks in advance


title: Can you programmatically change the button image for a ESRI.ArcGIS.Desktop.AddIns.Button?

text: I'm new to the AddIn customization model in ArcGIS 10 and I'm trying to figure out how to change a button image for a button. I'm looking for the equivalent to the BaseCommand UpdateBitmap method and haven't had any luck. Is this a limitation of the AddIn model? Or did I miss something obvious?


title: How to list raster bands?

text: I am interested in looping through a directory and listing only the red bands (Band 1) from 4-band NAIP imagery. From there I would like to mosaic all of the band 1 layers into a new single band raster image using Mosaic to New Raster. Is there a way to specify band 1 using `arcpy.ListRasters()` or should I find a new approach?


title: QGIS Spatial Queries

text: I was wondering if I could do the following within QGIS: I want to update a column in my point layer with the attribute information from a specific column in a polygon layer - so I have a layer of postcodes and I want to give it a column of say School Catchment and I want to populate this column with the value of the school name column in the catchment area polygon layer. Any help would be greatly appreciated. Thanks! Steve


title: converting coordinates from unusual DMS format to decimal

text: So I've been shared a dataset with what to me is an unusual format for degree min second data. In their own words:                   The coordinates provided are longitude and latitude in degrees minutes and seconds.         So for example, coordinates recorded in the spreadsheet as          X_Long: -10.0135 Y_Lat:  53.5726, can be taken as Longitude: 10º 01'35" W,          Latitude: 53º 57'26" N using standard notation.      I need to get this into decimal. Any quick way to do this ? Edit: in trying to solve this I made a following code:                   # Fake data         X_Long <- as.numeric(c(-2.4564, -1.0134, -8.01))              long.temp <- abs(X_Long - trunc(X_Long))         long.min <- trunc(long.temp*100)         long.sec <- trunc(as.integer((long.temp*100 - long.min)*100))      Its giving me weird rounding errors:                   > X_Long         [1] -2.4564 -1.0134 -8.0100         > long.temp         [1] 0.4564 0.0134 0.0100         > long.min         [1] 45  1  0         > long.sec         [1] 63 34 99      This seems to be coming from the trunc function - compare results here                   > long.temp*100         [1] 45.64  1.34  1.00         > trunc(long.temp*100)         [1] 45  1  0      Anyone know why this is happening ?


title: GPS tracking and Web mapping?

text: We are developing a tracking software and we are the stage where we test the reliability of our location stigma on the following web-maps   * Google Map   * Nokia HERE    * Open Street Map (OSM) We observe that our location appears accurate on Google maps, however when switched to Nokia or Openstreet there is a discrepancy of around 40-50 metres from our actual location. Most importantly our GPS position accuracy is around 2m. Could it be inaccuracy of map calibration or any similar problems known by anyone?


title: Compare features in PostGIS with updating attributes

text: I have two polygonal layers: "existing" and "modified". ![enter image description here](http://i.stack.imgur.com/GRyFn.jpg) Attribute table of "existing" and "modified" layers have same column "feature_id" - number of feature. Please help me create query which write attributes 't' (if features matches) and 'f' (if features not matches) in 'compare_result' column of "modified" layer using ST_Equals function. Now I receive attributes 't' (if features matches) and 'f' (if features not matches) in pgAdmin result window using query:               SELECT  ST_Equals(modified.the_geom, existing.the_geom) AS compare_result FROM  modified, existing      WHERE modified.feature_id=existing.feature_id      but I don't can create query? which write result in attribute table. Thanks! **EDIT:** Solution without creating additional layer:               UPDATE         modified     SET         compare_result = compare_result_temp          FROM         (             SELECT     r."feature_id" AS idmodified,     ST_Equals(r.the_geom, m.the_geom) AS compare_result_temp             FROM          "modified" AS r,          "existing" AS m     WHERE r."feature_id" = m."feature_id"         ) AS foo      WHERE       foo.idmodified = modified."feature_id";


title: IGeometryBridge.GetPoints throws NotImplementedException. Am I doing something wrong?

text: I am working on some .NET code that processes lots of polygon geometries at the vertex level. In order to reduce COM interop overhead, I want to call `IGeometryBridge2.GetPoints` once instead of calling `IPointCollection.Point` repeatedly:               // using ESRI.ArcGIS.Geometry;     static IPoint[] GetVertices(this IPolygon singlepartPolygon)     {         var pointCollection = (IPointCollection4)singlepartPolygon;         var points = new IPoint[pointCollection.PointCount];         var geometryBridge = (IGeometryBridge2)new GeometryEnvironment();         geometryBridge.GetPoints(pointCollection, 0, ref points);         return points;     }      `geometryBridge.GetPoints(…)` throws a `NotImplementedException`. (I have tried initialising the `points` array with empty `Point` instances before calling `GetPoints`, but this didn't make any difference. I tried `geometryBridge.QueryPoints(…)`, too, and this failed with a different exception.) Am I doing something wrong, or is this more likely to be a bug?


title: Bootstrap3 Resize map height for phone

text: I have a leaflet map in a bootstrap 3 html application. I have content both above and below the the map which has a height of 600px. When I access the application from a phone, the width of the map adjusts to the the width of the screen. However, the height of the map does not update and you end up with a seemingly endless map sandwiched between two chunks of content. The user must then endure endless scrolling. You can see what i mean here: http://punkmap.com/examples/fauxnmap/index.html Anyone know what's the best way to handle this?


title: Comparing field values to update new field using ArcGIS for Desktop?

text: I am currently doing an analysis in ArcMap 10 using python. I have two shape layers, one is a point layer and the other one is a line layer. The point layer have attributes for point's unique ID, as well as X and Y coordinates (expressed as coord_X, coord_Y). The line layer include line's ID, and its has coordinates for starting point and ending point. Here they are expressed as FromID_X, FromID_Y, ToID_X, ToID_Y. There are also two empty column in line layer's attribute table as FromID and ToID for points. All coordinates in both point and line layer are already know by field calculation, but the FromID and ToID of starting point and ending point in line layer is empty. I want to figure out the ID of point from point layer and fill them in the table of line layer, i.e. FromID and ToID. I have an idea to compare the coordinates of FromID_X (line's layer) with Coord_X(point's layer), if they are the same, output the Point's ID to line layer as FromID and repeat this for ToID. But now, I am not familiar with python in ArcMap 10, I am wondering if you can tell me some idea or grammar which I can look into so I could use python in ArcMap to do this work? Or you know some tools for such kind of work?


title: How do I automatically export a release of the EPSG registry?

text: I am trying to use the EPSG Registry API to export a release of the registry for offline use but am receiving a registry error. After following the Developer Guide I am HTTP POSTing the following XML request to the API URL:               <?xml version="1.0" encoding="UTF-8"?>     <GetRecords          xmlns="http://www.opengis.net/cat/csw"          xmlns:ogc="http://www.opengis.net/ogc"          startPosition="1" maxRecords="100"          outputFormat="application/xml; charset=UTF-8"          resultType="results">       <!--           Retrieve the ReleaseObject releated to the current Version History record.         Then use the ReleaseObject identifier to obtain the actual compressed          GML Dictionary file containing the EPSG Dataset via a getRepositoryItem          request.       -->        <Query typeNames="Association=a RegistryObject=release">         <ElementSetName typeNames="release">full</ElementSetName>         <Constraint version="1.1.0">           <ogc:Filter>             <ogc:And>               <!-- find the correct association -->                <ogc:PropertyIsEqualTo>                 <ogc:PropertyName>$a/targetObject</ogc:PropertyName>                 <ogc:Literal>urn:ogc:def:version-history:EPSG::6.15</ogc:Literal>               </ogc:PropertyIsEqualTo>               <ogc:PropertyIsEqualTo>                 <ogc:PropertyName>$a/associationType</ogc:PropertyName>                 <ogc:Literal>                   urn:x-ogp:def:AssociationType:EPSG:ReleaseFor</ogc:Literal>               </ogc:PropertyIsEqualTo>               <!-- now relate it to the ReleaseObject -->                <ogc:PropertyIsEqualTo>                 <ogc:PropertyName>$a/sourceObject</ogc:PropertyName>                 <ogc:PropertyName>$release/id</ogc:PropertyName>               </ogc:PropertyIsEqualTo>            </ogc:And>           </ogc:Filter>         </Constraint>       </Query>     </GetRecords>      The XML is copied from _section 5.3 'How To... Export a Release'_ of the guide, and the URL to send `GetRecords` requests to is suggested by _section 4.1.1 'GetRecordById'_. The XML error response I get is as follows:               <?xml version="1.0"?>     <ows:ExceptionReport xmlns:ows="http://www.opengis.net/ows" version="2.0">       <ows:Exception exceptionCode="InvalidRequest">         <ows:ExceptionText>XML Schema validation error</ows:ExceptionText>         <ows:ExceptionText>     Error:Error      Line:    2      Column:  193      Message: cvc-elt.1: Cannot find the declaration of element 'GetRecords'.         </ows:ExceptionText>                                                      </ows:Exception>     </ows:ExceptionReport>      I am using the `curl` command line program to run the query as follows:               curl -H 'Content-Type: appliation/xml' --data @request.xml 'http://www.epsg-registry.org/indicio/query'      Can you see where I'm going wrong? If not is there any other way to automatically download the registry? I am aware that I can login to the Registry HTML interface and download the registry from there but I would like to automate the process.


title: Australian Topo Map GeoTIFFs

text: Does anyone know if Australian topo map GeoTIFFs are made available anywhere (preferably for free)? LINZ makes them available in New Zealand - http://www.linz.govt.nz/topography/topo-maps/index.aspx ,I was hoping a similar service was available for Australia. Cheers, Gavin


title: Technical comparison of Oracle and PostgreSQL

text: I have searched on google and the archives for a technical comparison of Oracle and PostgreSQL (also MySQL), unsuccessfully. Is there any comparison that covers more than just the basic differences? Maybe you can give me some comparison?


title: geoserver and tiger roads layer overlay problem

text: Been playing with open layers, geoserver and the tiger roads shapefile but cannot get it to line up correctly over new york. I have looked online and followed the instructions regarding setting the projections and setting google maps to sphericalmercator = true however it still doesnt line up correctly. Any ideas?                       var map;             var untiled;             var tiled;             var pureCoverage = false;             // pink tile avoidance             OpenLayers.IMAGE_RELOAD_ATTEMPTS = 5;             // make OL compute scale according to WMS spec             OpenLayers.DOTS_PER_INCH = 25.4 / 0.28;                  function init(){                 format = 'image/png';                      var bounds = new OpenLayers.Bounds(                     -20037508, -20037508,                     20037508, 20037508.34                 );                      var options = {                     projection: new OpenLayers.Projection("EPSG:900913"),                     displayProjection: new OpenLayers.Projection("EPSG:4326"),                     maxExtent: bounds,                     //maxResolution: 0.00075764453125,                     units: 'm'                 };                     map = new OpenLayers.Map('map', options);                     var GoogleMap = new OpenLayers.Layer.Google(                     "Google Satellite", {                     type: G_SATELLITE_MAP,                     maxZoomLevel: 30,                     "sphericalMercator": true                     }                );                        // setup tiled layer                 wmsLayer = new OpenLayers.Layer.WMS(                     "tiger_roads - Tiled", "http://localhost:8080/geoserver/wms",                     {                         layers: 'tiger:tiger_roads',                         styles: '',                         format: format,                         tiled: 'true',                         transparent:'true',                         isBaseLayer: false                     }                 );                      map.addControl(new OpenLayers.Control.MouseDefaults());                  map.addControl(new OpenLayers.Control.Navigation());                 map.addControl(new OpenLayers.Control.ZoomBox({alwaysZoom:true}));                 // Coordinate display at bottom of map                 map.addControl(new OpenLayers.Control.MousePosition());                  map.addLayers([GoogleMap, wmsLayer]);                     var point = new OpenLayers.LonLat(-73.96543, 40.78885);                 // Need to convert zoom point to mercator too                point.transform(new OpenLayers.Projection("EPSG:4326"), map.getProjectionObject());                 map.setCenter(point, 12);                 }


title: Loading graph (.grf fle) into MXD using ArcPy?

text: I have hundreds of pairs of FGDB and associated MXD files. I need to insert a graph in each mxd linked to its corresponding FGDB. No problem creating graph, but stuck on how to get graph into mxd file using arcpy. If the graph creation code is run inside the mxd manually, the graph appears, but I need to do it in the python script loop.               import arcpy          workspace = r'C:\gdb_dir'     workspaces = arcpy.ListFiles('*.gdb')        # list of files to process          for workspace in workspaces:         name_prefix = workspace[0:-4]         current_MXD = name_prefix + '.mxd'       # this file will already exist              dataSrc = r'C:\name_prefix\Sample.gdb\DetailFeatures'         fieldX = 'NEAR_X'         fieldY = 'NEAR_Y'         graph_template = r'C:\temp\ScatterGraphTemplate.tee'         out_graph_name = "ScatterTest"         out_graph = r'c:\graphs\' + name_prefix + '.grf'   # How to insert this into mxd?              graph = arcpy.Graph()              graph.addSeriesScatterPlot(dataSrc, fieldY, fieldX)         graph.graphAxis[0].title = "Near Y"                  graph.graphAxis[2].title = "Near X"         graph.graphPropsGeneral.title = "Test of Flake Property Graph"         arcpy.MakeGraph_management(graph_template, graph, out_graph_name)              arcpy.SaveGraph_management(out_graph_name, out_graph, "MAINTAIN_ASPECT_RATIO", 600, 375)      How do you insert out_graph into current_MXD?


title: Meaning of red cross icon next to Geoprocessing tool/model?

text: I am using ArcGIS 10.1 for Desktop. When I double-click on the toolbox command line, it is not working and appears as a red cross symbol in command tool icon. So please, what can I do?


title: Does EO-01 Hyperion L1GST data show DN values, radiance or reflectance?

text: I am working on Hyperion hyperspectral imagery and I cannot figure out whether L1GST level does do show DN values, radiance or reflectance. * * *


title: Conversion from Namibian Coordinates to Google Maps (WGS84)

text: I have coordinates in Namibia, but don't know how to convert them to use them in Google Maps. The system is Lo22/19°. Then Y: +55 435,89 X: -498 157,28 Can anyone help with a conversion?


title: Apparently MOBAC is not generating MBTiles correctly

text: I am trying to generate MBTiles with MOBAC, to be used in an offline map Android application (trying to use OSMDroid). Now I generate the atlas, but something goes wrong because I get the error > android.database.sqlite.SQLiteException: no such table: images (code 1): , > while compiling: SELECT tile_data FROM images LIMIT 0,1 This is my MOBAC setup: ![enter image description here](http://i.stack.imgur.com/K0i7i.jpg)


title: Erdas python scripting help

text: I want to write scripts for Erdas in Python. I downloaded some Python examples here but I could not make much. Where can I get the full Help material for Erdas python scripting?


title: Creating a map cache without ArcServer?

text: Is there a way to create a tiled map cache without using ArcServer? I have aerial photography in .sid format that takes a while to load, whereas cached jpeg tiles will load much much faster. Any way I can work around this issue?


title: Polygon with the same points over and over

text: I am in the process of reading shape files (provided to us by the vendor dealing with GIS products) using a C# program and loading them into a SQL database. I first read the text for the shape and then use that text to update my geometry and geography types. And I am doing this as part of a SSIS package. My SSIS package failed on one shape file and threw me an exception: > A .NET Framework error occurred during execution of user-defined routine or > aggregate "geometry": System.FormatException: 24305: The Polygon input is > not valid because the ring does not have enough distinct points. Each ring > of a polygon must contain at least three distinct points. After running some manual updates, I found the record that was giving me grief. Then I took out parts of the multi-polygon and ran individual selects on them, creating individual polygons, eg.               SELECT geometry::STGeomFromText('POLYGON(( 118.501323586697 -20.3203577291617, 118.504216161911 -20.3220101539757, 118.502671059623 -20.3212136027048, 118.501323586697 -20.3203577291617 ))', 4326).MakeValid()      Finally, I found the polygon that was the problem shape: `'POLYGON(( 118.860739531873 -20.2274797478397, 118.860739531873 -20.2274797478397, 118.860739531873 -20.2274797478397, 118.860739531873 -20.2274797478397 ))'` I am trying to understand how to fix this so my SSIS pckage runs and translates and loads my SQL tables without me having to fix up anything manually. Any help will be greatly appreciated. Thanks


title: Problem with Coordinate Tolerances when I export my .shp files to .mdb

text: ![enter image description here](http://i.stack.imgur.com/GLRXh.jpg) When I export my shapefile to a personal geodatabase using ArcMap, the result is not the same as the original database. I do not know how to solve this problem.


title: vb.net arcobjects validate topology function only works after I use arcmap topology tool validate topology first

text: I find a strange thing occasionally, I have a validate topology function developed by vb.net arcobjects. It work not properly. The case is:   1. open arcmap   2. start editing   3. use my validate topology function, nothing happen, no error find (there really several errors)   4. set the extent to global and use arcmap topology toolbar 'validate topology in current extent' tool to validate topology   5. open error inspector and get errors list (4 errors)   6. change the current extent and use arcmap topology toolbar 'validate topology in current extent' tool to validate topology   7. open error inspector and get errors list (2 errors)   8. remain the same extent as step 6, use validate topology function developed by vb.net arcobjects, all errors displayed in error inspector. I want to know is it my code bug or operator issue? My code:                   ' Step 1 - access a topology     Dim topoUiD As UID = New UIDClass     topoUiD.Value = "esriEditorExt.TopologyExtension"     Dim m_application As IApplication = My.ArcMap.Application     Dim topologyExt As ITopologyExtension = CType(m_application.FindExtensionByCLSID(topoUiD), ITopologyExtension)  'original code from esri     Dim topology As ITopology          topology = CType(topologyExt.CurrentTopology, ITopology)     MsgBox("successfully access current topology")          ' Step 2 - Finding topology errors     Dim geoDS As IGeoDataset = CType(topology, IGeoDataset)     Dim errorContainer As IErrorFeatureContainer = CType(topology, IErrorFeatureContainer)     Dim eErrorFeat As IEnumTopologyErrorFeature     eErrorFeat = errorContainer.ErrorFeaturesByRuleType(geoDS.SpatialReference, esriTopologyRuleType.esriTRTLineCoveredByLineClass, Nothing, True, False)     Dim topoError As ITopologyErrorFeature          ' check the eErrorFeat length to determine whether has error     topoError = eErrorFeat.Next          topologyExt.ClearActiveErrors(esriTEEventHint.esriTENone)          topologyExt.DelayEvents(True)          Do While Not (topoError Is Nothing)         MsgBox("Topology error exists!")         topologyExt.AddActiveError(topoError, esriTEEventHint.esriTENone)         topoError = eErrorFeat.Next     Loop          topologyExt.DelayEvents(False)      Another thing is if I am not in editing session and use validate topology tool developed by vb.net arcobjects, the map will crash.


title: Select by Attribute --Top 500 values selection in attribute table - ArcGis

text: Is there any way we can select top 500 values in attribute table? I tried by sorting area field ascending/descending... but I can't count big numbers to select. In microsoft excell we can do know, similarly in arcgis, is there any code to keep in attribute selection ??


title: How to call ESRISignAddIn Silently (no UI)?

text: I am creating a build workflow and I ran into some problems. When I call **ESRISignAddIn.exe** from the command prompt it is throwing up a dialog saying "Add-In successfully signed". I was expecting it to report that back within the cmd dialog and not open up a new messagebox. Have anyone else seen this or have worked out how to sign an add-in "silently"? The /s switch does not work. This issue occurs with the 10.2 version of the exe (latest at time) and earlier versions.               "c:\Program Files (x86)\Common Files\ArcGIS\bin\ESRISignAddIn.exe" "C:\temp\my.esriAddIn" /c:"C:\temp\newCert.pfx" /p:myPass


title: Wildlife surveys using tablet and GPS, import/export to/from ArcMap

text: I've found some info about this on a previous post (Advanced android app for mapping and GPS logging), but it didn't provide enough info, I'm hoping some other tech experts can help! I have a tablet (Android) and it has integrated GPS, I would like to run a (preferably free) program into which I can input polygon layers (boundary units) and then track our flight path as we move through those polygon units. We are talking fairly large scale (units are at least 30 km²). There's no wifi, and we require something that keeps up to 100kph. Our current GPSs do this fine, but we can't import polygons into the GPSs. We have used ArcPad in the past, but it is too slow. To summarize - we want a tablet-based program to which we can import GIS polygon layers and then watch our track on the tablet as we fly transects in the polygon units. We would like to export these track to ArcMap when we're done, so pretty much any format will work. I've looked into a couple options (Ozi, Google Earth), but I'd like to know if anyone has any experience with tracking GPS on a tablet at higher speeds and across large areas. Thanks!


title: Trilateration algorithm for n amount of points

text: I need to find algorithm that can calculate centroid **A** (aka gravity center, geometric center, center of mass) from the figure where circles T1,T2,T3,T4,T5,..,Tn intersect AND length of line **R** from centroid to farthest corner of mentioned figure Following information is given:   * T1 Latitude = 56.999883 Longitude = 24.144473 Radius = 943   * T2 Latitude = 57.005352 Longitude = 24.151168 Radius = 857   * T3 Latitude = 57.005352 Longitude = 24.163356 Radius = 714   * T4 Latitude = 56.999042 Longitude = 24.168506 Radius = 714   * T5 Latitude = 56.994226 Longitude = 24.15709 Radius = 771 Result should look like this: A Latitude = XX.XXXXXXX Longitude = XX.XXXXXXX Radius = XX ![enter image description here](http://i.stack.imgur.com/kng2D.gif) As you probably already figured out, I am working on software that can find device location by closest Wifi Access Points or Mobile Base stations, as number of access points or base stations might change, I need an algorithm that can adapt to uncertain amount of points. There are some similar questions here and here, but none of them exactly answers to my question. Thanks!


title: Synching Two ArcMap JavaScript webmaps to respond to panning

text: I am looking for some advice on how to approach a problem I am having. I am looking to build an ArcGIS webmap using JavaScript and HTML 5 that reacts to a user panning on two separate, synchronous maps. Modest Maps provides an excellent example of Loading Two Maps and Synching Their Panning that helps visualize what I need to do. However, I need to accomplish the _Loading Two Maps and Synching Their Panning_ with ArcMap products. The Overview dijit in Arcmap seems to be capturing the essence of what I am trying accomplishing. However, the overview dijit in the widget overviewmap seems to react rather then pan synchronously. The dijit css feeds into this map and can be used as a resource if needed. I couldn't find anything in it that would accomplish synchronous panning. I am open to any ArcGIS suite methodology to accomplish the same as the Modest Maps look and feel. If you have any suggestions of comments I would greatly appreciate it.