Pan/Scale/Navigation Buttons Disappearing

I am working on an OpenLayers 2.10 map using Internet Explorer 8 (used by work), and am encountering some issues with the pan/navigation buttons disappearing after sepcific events:   1. If I change the layer to another background, the pan/navigation panel buttons lose their picture   2. If I resize the browser window, the pan/navigation panel buttons lose their picture When one of these events occurs, the buttons show the broken image link icon. If I refresh the page, the buttons come back. After perusing some answers, the only similar solutions I came up with were at OpenLayers - redrawing map after container resize using the map.updatesize command. These commands did not seem to work, though I am not entirely sure where to place them. Also, I just started using OpenLayers a couple months ago, and a solution may exist using terminology I am not aware of. I would like a solution that keeps the pan/navigation/zoom panel buttons constantly displayed, even if layers are changed or if the page is resized. It is also acceptable for the buttons to refresh after an event that will cause the button to no longer load. Thanks in advance for any guidance.


Select parcels from SDE layer using table list and exporting to shapefiles

Using a table of county fips codes and state name, I need to select from a large parcel layer stored in SDE and export to individual shapefiles based on the fips code and state name from the table. Can this be done with a Python script or ModelBuilder?


Is there demand for a service to publish web maps with lots of shapes?

Is anyone in the GIS community trying to publish choropleth web maps with a large (10s to 100s of thousands) number of polygons? and if so what's your current approach and what are the shortcomings of that approach? My company is thinking about reviving an old project, but would like to get a sense of the demand before we invest too much time in something that may not be needed. Features of the service would include:   * Support for millions of polygons.   * Ultra fast updates (recolor the maps in < 1 second, our goal is 60ms).   * Realtime query support (mouseover or click a polygon to display it's info).   * Bring your own data (upload your own shapes, classifications, etc)   * Built in Census Data (all census layers and data preloaded)   * Web access (no coding required)   * API access (REST API with clients in Python and Javascript to start)   * Overlay (display your published maps as overlays in Google Maps, Bing Maps, Open Layers, etc) The core of this system is already working, but it would be a pretty major push wrap it in something useable. We're unsure at this point if there is enough demand to justify the effort. My apologies in advance if this isn't the appropriate forum for this sort of question.


Table merging: multiple corresponding objects

I'm working with a shapefile which has all the municipalities of France as areas, identificated with their unique code (`INSEE` code). I want to join a table, where every historical monument is included, it has the corresponding `INSEE` column that I may use to merge. But sometimes there can be more than one monument in a municipality (city) identified as its unique `INSEE` code. And when I merge tables together, I would like it to add automatically objects so that every monument is added. If the city of Paris counts 50 monuments for example, I would like to see a pile of 50 "Paris" on my map, each one having a monument. Because when I try to join tables on MaPublisher or Global Mapper, it is not adding objects, it is just merging tables with a unique `INSEE` code. That way, there is only one monument included, because there is no room for others. More understandably, and with simplifications : 1) I have a `CSV` or `TXT` or `XLS`               INSEE;MONUMENT NAME;     10223;Eiffel Tower;     25685;Castle;     25685;City hall;      2) and a Shapefile with objects as areas               INSEE;CITY NAME     10223;PARIS     25685;LYONS      3) I'm getting a shapefile result :               INSEE;CITY NAME;MONUMENT NAME     10223;PARIS;Eiffel Tower     25685;LYONS;Castle      4 ) But what I want is :               INSEE;CITY NAME;MONUMENT NAME     10223;PARIS;Eiffel Tower     25685;LYONS;Castle     25685;LYONS;City Hall      Can you help me and tell me how may I do ? And you are willing to do that for me, I give you the files for you to download : `PROJECTION : LAMBERT 93 (FRENCH)` -it has extra miscellaneous attributes columns- LINK OF THE FOLDER Thank you :)


Problems with french coordinate systems (i.e. Lambert Zone I, EPSG 27571)

I'm working on a coordinate library implementation and I seem to be stuck with the French Lambert coordinate systems. I'm using this page as a reference (can anybody confirm the results there are accurate?) The individual steps in the conversion seem to fit, but my overall result differs up to 50m from the one computed on that web site. Now step by step: I'm converting as an example the point 2.1°E, 49.8°N from WGS84 to EPSG27571 (NTF Paris / Lambert Zone I)   1. WGS84 -> WGS84 Geocentric: 4122109.79, 151150.90, 4848460.24 (meters)   2. WGS84 Geocentric -> NTF (PM Paris) Datum conversion using Bursa-Wolfe: 4122273.66, 151210.75, 4848135.39 (meters)   3. Geocentric -> Geographic, including meridian shift to Paris: -0.004127387, 0.86917528, -49.58 (radians) I'm quite sure the result is right so far, because these are the same values than said web page returns when using EPSG 4807 as target CS (and converting from grad(!) to radians).   1. Geographic -> Projected coordinates using lambert conformal 1sp, with latitude of origin = 0.86393797974 rad, longitude of origin = 0, False easting 60000, False northing 1200000, Ellipsoid Clarke 1880 IGN: 582974.00, 1233402.04, while expected is 582976.09, 1233397.93. Offset is about 3-5m here.  I have verified the lambert conformal projection using an example for EPSG 24200 JAD96 Jamaica, therefore I _also_ think that the projection itself is right. The error increases the farther away from Paris I get. Not including the false easting or false northing values, the error (quotient between right and wrong values) is in the order of 0.997-0.9998. This looks like some scaling problem, and the WKT description for EPSG 27571 mentions some scale, but I have no clue where to apply that. I do have a scale in the datum, though. Now: Where is the error?


I would like to know how to clean topological errors

I am a new QGIS 2.0 user. I was so impressed the new default tool, Topology Checker, in this new version. But, hundreds errors are bothering me for cleaning up features. Is it possible to clean errors automatically, not manually? Thank you in advance for your concern. * * * Thank you for your quick response, Mr. Solin and Mr. radouxju. If you mind, let me answer additional 2 questions. 1. Is it possible to clean topological error in QGIS not GRASS. 2. According to "v.clean tool", is it possible to clean overlapped areas. In my understanding, there is no specific tool for overlap error  Thank you so much.


How to Map Latitude/Longitude to FAA Chart?

How can I map latitude/longitude to an FAA digital chart? Attached is a sample GEOTIFF. More specifically, how can I get a reference points from the digital data? Once I can get reference points, I have no problem translating Pixel-Lat/Long. In my testing, I have been manually identifying the reference points to the the translation. This is a rather tedious process. It looks like I just need a single lat/long-pixel-x/y mapping to derive them all from this data. The GEOTIFF spec suggests the Model Tie Point is the key but the example in the spec does not match this data. -822739.829757951 -64249.12941973 does not give units not does it give offset from what? In this chart the false origin is located at the Pixel coordinates (6521 3865) from the upper left corner. Thanks               ExifTool Version Number         : 9.53     File Name                       : CF-19 WAC 43 South.tif     Directory                       : .     File Size                       : 47 MB     File Modification Date/Time     : 2013:03:15 10:45:54-04:00     File Access Date/Time           : 2014:02:28 13:45:59-05:00     File Inode Change Date/Time     : 2014:02:28 13:06:22-05:00     File Permissions                : rwxr-xr-x     File Type                       : TIFF     MIME Type                       : image/tiff     Exif Byte Order                 : Little-endian (Intel, II)     Image Width                     : 11888     Image Height                    : 4135     Bits Per Sample                 : 8     Compression                     : Uncompressed     Photometric Interpretation      : RGB Palette     Strip Offsets                   : (Binary data 286 bytes, use -b option to extract)     Samples Per Pixel               : 1     Rows Per Strip                  : 128     Strip Byte Counts               : (Binary data 262 bytes, use -b option to extract)     Planar Configuration            : Chunky     Color Map                       : (Binary data 1536 bytes, use -b option to extract)     Sample Format                   : Unsigned     Pixel Scale                     : 127.09668 127.09668 0     Model Tie Point                 : 0 0 0 -822739.829757951 -64249.12941973 0     GDAL Metadata                   : <GDALMetadata>.  <Item name="PyramidResamplingType" domain="ESRI">NEAREST</Item>.  <Item name="RepresentationType" sample="0">THEMATIC</Item>.  <Item name="STATISTICS_MINIMUM" sample="0">0</Item>.  <Item name="STATISTICS_MAXIMUM" sample="0">255</Item>.  <Item name="STATISTICS_MEAN" sample="0">108.47642098115</Item>.  <Item name="STATISTICS_STDDEV" sample="0">83.473159026412</Item>.</GDALMetadata>.     GDAL No Data                    : 256     Geo Tiff Version                : 1.1.0     GT Model Type                   : Projected     GT Raster Type                  : Pixel Is Area     GT Citation                     : PCS Name = CF-19 WAC     Geographic Type                 : NAD83     Geog Citation                   : GCS_North_American_1983     Geog Angular Units              : Angular Degree     Geog Semi Major Axis            : 6378137     Geog Inv Flattening             : 298.257222101     Projected CS Type               : User Defined     Projection                      : User Defined     Proj Coord Trans                : Lambert Conf Conic 2SP     Proj Linear Units               : Linear Meter     Proj Std Parallel 1             : 46.6666666666667     Proj Std Parallel 2             : 41.3333333333333     Proj False Origin Long          : -69     Proj False Origin Lat           : 44.41667     Proj False Origin Easting       : 0     Proj False Origin Northing      : 0     Image Size                      : 11888x4135


Geocoding house numbers

I have a point feature class with residential information, but all that it has in terms of geographic information is the residence's house number. There is no other information (street, zip, etc). JUST house numbers. The points fall more or less where they are located, but they are somewhat off. I'm working off of a ArcMap basemap and I sometimes I can just move the points to areas that are obvious to me, for example if there's one point in a general area close to a house, and that's the only house, then its obvious to me that that point is that house. However, i need a faster way to do this, and I need help with points that aren't so obvious, like for example when i have to assign a point to a residence and there are several residences to choose from. I thought about geocoding, but I'm not sure how this would work when all I have to work with is a house number. So I thought I could download a shapefile (maybe tiger data) with street names as well as zip codes, and do some kind of spatial join to where the points (house numbers) would be assigned the street name and zip code nearest the point. and then create an address locator to match the points to the actual locations based on those attributes (street name, house number, and zip code). Would this be a good way to tackle this? Does anyone have any other suggestions that might be better? Could you maybe explain to me how I might do this? I'm pretty new to GIS and not very familiar with geocoding.   Also, can someone explain to me more or less how I should perform this spatial join? Thanks!!


Calculate lat lon bounds for individual tile generated from gdal2tiles

I have multiple tile sources generated using gdal2tiles I'd like to present on the same map. So when presenting a tile I need to determine which source to serve it up from, checking against the bounds. Anyone know how to calculate the lat lon bounds for a single tile based purely on zoom, x, and y (from the file structure generated by gdal2tiles)? BTW: I'm using Google Maps API v3 in case there's a need to call some functionality from the API to help with calculations. The reason I'm restricted to purely zoom, x, and y is because the tiles aren't just getting called by an overlay on the map but also some custom print functionality that allows printing outside of the map. Thanks in advance for any help or suggestions anyone has. Cheers, Gavin


How to load Spatialite extension via SQLite3 ODBC Driver?

I am having trouble getting the Spatialite extension to load as part of the SQLite3 ODBC driver...[aka the joys of making software somewhat spatial when it tries so hard not to be]. **What works [loading extension into SQLite Driver]:** I have software that is connected using the SQLite3 ODBC driver (http://www.ch-werner.de/sqliteodbc/). I have previously loaded extensions into this driver without any problems. A previous working example was a trigger that updated columns [using triggers] using some trig functions that are not part of base SQLite functionality - I built the *.so extension and pointed SQLite to it as part of the ODBC connection and everything worked fine - when I connect/export data from the software I got no issues. **What works [loading Spatialite extension into SQLite]:** If I use the SQLite shell I can load the Spatialite (or any other extension) using: select load_extension('libspatialite-2.dll');. That works without fail. **What doesn't work:** However I cannot seem to get the Spatialite extension to load properly as part of the ODBC connection - and I am using a trigger that builds the geometry column of a Spatialite database based on the software providing only x and y values. I get an error that says "Extension C:\windows\system32\libspatialite-2.dll did not load: the specified module could not be found". I also tried the spatialite-4.dll and it doesn't work either. The module and its requirements are located there and the environmental path has (always) been there. If I do not load the extension then I will get an error when the trigger calls a function (makegeometry, makepoint or whatever you want) that is not part of the SQLite driver - it doesn't understand the function. My hack workaround at this point are batch files that run what the trigger should do after the data has transferred from the software to the Spatialite database, but I'd prefer it to be a one-stop shop. Any thoughts would be great. Thank you. I haven't seen this much on the web, other than a Google Groups question (https://groups.google.com/forum/#!topic/spatialite-users/3ImfXm7cdpU) that didn't resolve the problem [I haven't tried to build the CppDB yet - I want to see if there is a more ~~lazy~~ straightforward solution].


How to perform Intersection using ModelBuilder iteration?

![Pic of Model](http://i.stack.imgur.com/UHZrG.jpg)I would like to take the intersection of two areas (one randomly generated), and do this n times. I can do it and get one result, but I can't save all the intersection areas to a file at the end. The problem is that when I run my model more than once it says that my buffer output file already exists, which I can imagine would happen again for the intersection. What I'd like to happen is to run my model where I produce new random points everytime, then buffer them and intersect the dissolved buffer with another area and then get out the area. Then I want to do this several times. If this isn't worded very well then just let me know and I can explain what I mean in another way. Thanks a lot


SQL: Select aggregated Bounding Box of multiple Features

I would like to calculate the aggregated bounding box of multiple features using OGC Simple Feature SQL or SQL API of ArcSDE. The following statement will only query the bounding box of each feature:               select ST_Boundary (shape) from polygon_table where country like '%usa%'      What I needed is one bounding box aggregated over all features, resulting to the condition in the where clause.


Why am I unable to download Landsat 8 geotiff data from the USGS site?

I seem to have an issue with not being able to download geotiff Landsat 8 data from te USGS site. After i find the image i would like to download and proceed to download the Level 1 geotiff option is disabled and I am only allowed to download the LandsatLook jpeg images. I downloaded data from this site before but perhaps I am overlooking something. I tried both http://glovis.usgs.gov/ and http://earthexplorer.usgs.gov/ Screen capture showing the disabled "dimmed" option: ![enter image description here](http://i.stack.imgur.com/dtLjo.jpg)


ArcSDE: connecting via LDAP Authentication

I was wondering if anyone has a sample connection string for SQLPlus using LDAP Authentication. I am also trying to specify in the connection string the database I wish to use. Thanks, CDB


Using pyqgis with existing ArcGIS files

I'm pretty new to GIS. I have some files that seem to have been created using ArcGIS (I got them from here http://www.fao.org/ag/againfo/resources/en/glw/GLW_dens.html). I would like to do some analysis on these maps, preferably using python. I'm wondering if it's possible to use pyqgis to do this kind of thing. I'm really just looking for some info to get me started. e.g. - Is it possible to import these files into a pyqgis application? - Are there any tutorials that might help me get started? Thanks. Edit: The file types I see are .lyr, .aux, .avl, .dat, .nit, .adf Here's a specific example of what I'd like to do. Say I have two maps (map A and map B). I think these are raster images, and let's say that the "colour" in a given cell in map A measures quantity qA. And likewise for map B. Something simple I'd like to do is to produce a third map (map C) whose cells have quantities given by qC = qA*qB.


Control Map Pan per checkbox

i want to add the functionality of switching the option to pan the map and put it in a checkbox. But im not quite sure how i can do that. The reason I want to do this is: I have a layer which displays features and has a select and hover control. So theres no way to pan the map, because the clickEvent triggers direct action from the Handlers (click,hover). I want to add a checkbox where i can activate/deactivate the possibility to pan the map, so i can pan the map when there are features present, which normally would trigger a reaction to the mouseover the feature. The way ive tried to do it is this one but it didnt work. I attached the function to the checkbox but it never worked :/ ` function allowPan() {               if(panAct == false){             highlightCtrl.activate();             selectedControl.activate();             alert("controls activated!!")             panAct = true;         }          if(panAct == true){          highlightCtrl.deactivate();         selectedControl.deactivate();         alert("controls deactivated!!");         panAct = false;         }      } `


sa.fill creates NoData raster in standalone script ArcMap 10.0

I am trying to put together a geoprocessing workflow in Python, and I've hit a stumbling block with sa.fill. When run as a standalone script, the tool produces a raster composed entirely of NoData values. If I run the same code in the Python window in ArcMap, a correctly filled raster is produced. My code snippet is:               outfill = sa.Fill(raster,"")     outfill.save(fill)     print "Sink fill complete."      I've tried copying my entire script into the Python window and adding it to a toolbox, but I get the following error in both cases:               Runtime error <type 'exceptions.EOFError'>: EOF when reading a line      The standalone script does not produce this error. It actually runs through all the steps without incident, but the output is useless as everything after the Fill step is based on NoData values. Does anyone have any thoughts on this?


Creating new shapefile based on its attributes with python

I have a shapefile which has an attribute year. I used to generally copy and paste these features in the dataframe and then using the definition query just show polygons which were established before that year. I want to know how can I do it using python. Like if I can give the input feature class, add all the years I want the different shapefiles for and the workspace where it can create shapefile for the mentioned years with the right naming convention. Thank You.


How to calculate polygon centroids in R (for non-contiguous shapes)

I've spent a little while figuring out the answer to this question. It's not immediately obvious from a Google search, so thought it may useful to post the answer on here. There is also an additional question about non-contiguous polygons. Instant easy answer: **use the command:**               centroids <- getSpPPolygonsLabptSlots(polys)      (This was found in the class description of the SpatialPolygonsDataFrame R data class for the overarching spatial package in R, sp) This seems to do exactly the same thing as               cents <- SpatialPointsDataFrame(coords=cents, data=sids@data, proj4string=CRS("+proj=longlat +ellps=clrk66"))      in the following code, which should be replicable on any R installation (try it!)               #Rcentroids     install.packages("GISTools")     library(GISTools)     sids <- readShapePoly(system.file("shapes/sids.shp", package="maptools")[1],                            proj4string=CRS("+proj=longlat +ellps=clrk66"))     class(sids)     plot(sids)     writeSpatialShape(sids, "sids")     cents <- coordinates(sids)     cents <- SpatialPointsDataFrame(coords=cents, data=sids@data,                        proj4string=CRS("+proj=longlat +ellps=clrk66"))     points(cents, col = "Blue")     writeSpatialShape(cents, "cents")          centroids <- getSpPPolygonsLabptSlots(sids)     points(centroids, pch = 3, col = "Red")      Where cents (blue) and centroids (red) are identical centroids (this should plot should appear after you've run the code): ![centroids calculated by R](http://i.stack.imgur.com/qpCm7.png) So far so good. But when you calculate polygon centroids in QGIS (menu: Vector | Geometry | Polygon Centroids ), there are slightly different results for non-contiguous polygons: ![QGIS generated polygons](http://i.stack.imgur.com/zY5qq.png) So this question is 3-things:   1. A quick and easy answer   2. A warning for people using R to calculate centroids for non-contiguous polygons   3. A question about how it should be done in R to properly account for multi-part (non-contiguous) polygons


building opentripplaner gtfs.zip

I am building a journey planner app using opentripplanner. i have followed their two minutes tutorial in which we just run webapp which uses pre build graph.obj file, it works fine. But in next five minutes tutorial they shows how can i build .obj from a zip gtfs, while i try to do so i got these errors, i try had to figure them out ! ![enter image description here](http://i.stack.imgur.com/FGDrF.png) ![enter image description here](http://i.stack.imgur.com/qkG0T.png) ![enter image description here](http://i.stack.imgur.com/GkjdW.png)


Why doesn't Polygonizer copy the line ids to my polygonized shapefile?

I am trying to convert a line shapefile into a polygon shapefile using Polygonizer plugin, and when my new polygons are created all the cells in my id field are "Null", although I have clicked "copy table field from line layer". I have tried both Polygonizer 2 and Polygonizer 2.1. Any idea why the ids from my line layer are not copied to my polygonized shapefile? Thank you!


Problem looking at tables within ESRI geodatbase

I installed the most recent OSGeo4W with QGIS 1.8.0 and GDAL 1.9.2 on my Windows 7 PC. I can load an ESRI geodatabase (.gdb file) and I can get the window "Select vector layers to add". When I click on any of the layers, however, I get the message "qgis.exe has stopped working" Any idea where the problem might be?


How to get the nonspatial tables in the TOC into a combobox in ArcMap?

I'm new in using Arcobject . i'd like to set dbf of ArcMap TOC to the Combobox Items. How can i set dbf files to items of combobox using C# programming and Arcobject ?


How to add a DEM Band to an existing 3-band TIFF file using QGIS?

I'm quite new in the GIS business and I do not figure out how to do the following thing: I have two geoTIFF Files, one is a 3+1 bands file (RGB+alpha), the second one is a DEM geoTIFF file containing one band with elevation data. For a certain application, what I need is a 1+3+1 bands geoTIFF, containing the elevation Data on the first band, and the optical data plus the alpha info in the 2nd to 4th bands. How can I join or add bands from different geoTIFF's using Quantum GIS, for instance? I have found a similar question here, but there does not seem to exist a definitive answer there. So, I would be very happy to get help from the experts on this platform, Thanks in advance.


Geoserver - problem with MySQL support (impossible to create MySQL store at all)

I need to access spatial data in MySQL database with Geoserver. Official user manual states: just download package, extract it into WEB-INF/lib directory and all will be OK. Instead, just nothing happens - MySQL option doesn't appear when creating new data source, I can't find any relevant messages in logs and Google doesn't provide a hint. What do I do:   1. Download MySQL extension, following recommended link. I have Geoserver version 2.4.3 (uses jre7, downloaded recently). Downloaded archive declares to be of version 2.4.3 (and manual specifically states that versions should match).   2. Unzip archive and copy it contents to `[Geoserver installation directory]\webapps\geoserver\WEB-INF\lib` (also tried other paths mentioned in the manual and extension readme, though they don't make any sense - result the same).   3. Restart Geoserver service (also tried rebooting system), than go to _Data -> Stores -> Add new store_ in Geoserver web interface (logged in as admin). Manual states that MySQL option should be there, but it isn't. As mentioned, I use Geoserver 2.4.3 (latest at current time), OS is Windows 7 x32. Geoserver is installed as a service and starts automatically. Logging profile of server is set to VERBOSE_LOGGING, no traces of "mysql" are found. System logs don't contain anything but messages from Service Control Manager about service being started/stopped. Trying to mess with other JARs in the same directory causes Geoserver to fail. There was installed Geoserver 2.1-beta3 some time ago, but current installation is made into different directory. I'm interested in any advice how to get Geoserver-MySQL combination working, or how to diagnose why it is problematic/impossible.


gtfs-editor compile time error

I'm trying to build the gtfs_editor project by following the instructions on the given page. But, while executing the command **play run** the log at `localhost:9000` says `Compilation error The file /app/controllers/Check.java could not be compiled. Error raised is : The type Check is already defined` Has anyone gone through this problem before? If yes then your help on this issue is much appreciated. Thanks in advance. p.s.- I'm working on `play framework -1.2.5.3`,


Direct Georeferencing with X,Y,Z and Omega, Phi, Kappa

I have images taken from plane, it's longitude, latitude, elevation, omega, phi, kappa. I suppose that earth is flat (I don't need any DEM on that stage) and I want to calculate how image is nailed to the ground. I want to find center and 4 corners of nailed image to the ground. Could you advice some simple formula for calculation? Or free tool that has such calculations just to check my own implementation results. I tried to use GRASS but without any luck, it is too complicated :( Thanks in advance! **UPDATE:** Images were taken from around 3.2k meters without physical GCPs. The precision that we got is around 20 meters with camera angles (Omega, Phi) around 2 grad. Do you think it is possible to get better precision? We are using C# and here is our implementation:                   /// <summary>         ///          /// </summary>         /// <param name="omega">camera angle, in grads</param>         /// <param name="phi">camera angle, in grads</param>         /// <param name="kappa">camera angle, in grads</param>         /// <param name="longitude">camera position</param>         /// <param name="latitude">camera position</param>         /// <param name="altitude">camera position</param>         /// <param name="altOverGround">average flight height, it is defined on flight planning stage</param>         /// <param name="imageLength">it is defined on flight planning stage, calculated from altOverGround and focal lenght, it is in meters</param>         /// <param name="imageWidth">it is defined on flight planning stage, calculated from altOverGround and focal lenght, it is in meters</param>         /// <returns></returns>         public static DbGeography CalculateProjection(double omega, double phi, double kappa, double longitude, double latitude, double altitude,             double altOverGround, double imageLength, double imageWidth)         {             // coordinates of the camera             DbGeography cameraCoordinates = CreatePoint(longitude, latitude);                  // converting angles from grad to radian             omega = ConvertFromGradToRadian(omega);             phi = ConvertFromGradToRadian(phi);             kappa = ConvertFromGradToRadian(kappa);                  // get image center (projected to groud) as offset to camera coordinates             // result of imageCenter is in image coordinate system             Point imageCenter = GetImageCenter(omega, phi, kappa, altitude);                  // calculating real image length for this height             var realImageLength = imageCenter.Z * imageLength / altOverGround;             var realImageWidth = imageCenter.Z * imageWidth / altOverGround;                  // corners offset in image coordinates system             var xOfsfet = realImageLength / 2;             var yOffset = realImageWidth / 2;                  // corners list             List<Point> imageCorners = new List<Point>() {                  new Point(xOfsfet, yOffset, imageCenter.Z),                 new Point(-xOfsfet, yOffset, imageCenter.Z),                 new Point(-xOfsfet, -yOffset, imageCenter.Z),                 new Point(xOfsfet, -yOffset, imageCenter.Z)             };                  List<DbGeography> projectionCorners = new List<DbGeography>();             foreach (var point in imageCorners)             {                 // get image corner projected to the ground                 Point projectedPoint = GetProjectionCoordinates(omega, phi, kappa, altitude, point);                      // add projected point to list of projected corners                 projectionCorners.Add(CreatePointWithOffset(cameraCoordinates, projectedPoint.Y, projectedPoint.X));             }                  // add first point to complete the polygon             projectionCorners.Add(projectionCorners[0]);                  return CreatePolygonFromListOfPoints(projectionCorners);                     }              private static Point GetProjectionCoordinates(double omega, double phi, double kappa, double altitude, Point imageCorner)         {             // converting point coordinates from image's coordinate system to normal coordinate system             Point imagePoint = ConvertImageCoordToNormal(omega, phi, kappa, imageCorner);                  // get coordinates of the projected point             return GetProjectionCoordinates(altitude, imagePoint);         }              private static Point GetProjectionCoordinates(double altitude, Point imagePoint)         {             var x = altitude * imagePoint.X / imagePoint.Z;             var y = altitude * imagePoint.Y / imagePoint.Z;             return new Point(x, y, altitude);         }              private static Point ConvertImageCoordToNormal(double omega, double phi, double kappa, Point imageCorner)         {             double x1, x2, y1, y2, z1, z2, x, y, z;             kappa = Math.PI / 2 - kappa;             x1 = imageCorner.X * Math.Cos(kappa) + imageCorner.Y * Math.Sin(kappa);             y1 = -imageCorner.X * Math.Sin(kappa) + imageCorner.Y * Math.Cos(kappa);             z1 = imageCorner.Z;                  x2 = x1 * Math.Cos(phi) - z1 * Math.Sin(phi);             y2 = y1;             z2 = x1 * Math.Sin(phi) + z1 * Math.Cos(phi);                  x = x2;             y = y2 * Math.Cos(omega) + z2 * Math.Sin(omega);             z = -y2 * Math.Sin(omega) + z2 * Math.Cos(omega);             return new Point(x, y, z);         }              /// <summary>         /// Calculates altitude of the image center         /// </summary>         /// <param name="omega"></param>         /// <param name="phi"></param>         /// <param name="kappa"></param>         /// <param name="altitude"></param>         /// <returns>Altitude of image center</returns>         private static Point GetImageCenter(double omega, double phi, double kappa, double altitude)         {             // calculating coordinates of image center             var xCenter = altitude * Math.Tan(phi);             var yCenter = altitude * Math.Tan(omega);             var imageAltitude = altitude / (Math.Cos(omega) * Math.Cos(phi));                  // recalculating center with kappa value             var alpha = Math.Atan(yCenter / xCenter) + kappa;             var aa = Math.Sqrt(xCenter * xCenter + yCenter * yCenter);                  double imageCenterX = Math.Cos(alpha) * aa;             double imageCenterY = Math.Sin(alpha) * aa;             return new Point(imageCenterX, imageCenterY, imageAltitude);         }              private static double ConvertFromGradToRadian(double omega)         {             omega = Math.PI * omega / 200;             return omega;         }      and Point class:               class Point     {         public double X { get; set; }         public double Y { get; set; }         public double Z { get; set; }         public Point(double x, double y, double z)         {             X = x;             Y = y;             Z = z;         }     }


Print FeatureLayer labels in Esri JavaScript API 3.5 and earlier

I have a need to implement labels on features in a FeatureLayer in version 3.5 if Esri's Javascript API. The labels are from a field returned by a REST feature service. I can't move to 3.7 for various reasons at this time. I have tried using a TextSymbol but my map features just turn to the color of the TextSymbol and no text appears. I may be approaching this in he wrong manner, though. Below is the code I'm attempting to use for labeling with the featureLayer object being my instance of the FeatureLayer I'm adding to the map. Is there a different or proper way to accomplish this task?               featureLayer.on("graphic-add", function (evt) {      var labelColor = new Color([255, 0, 0, 0.25]);     var myLabel = new TextSymbol(evt.graphic.attributes["My Field Name"]);          myLabel.setColor(labelColor);     myLabel.font.setSize("14pt");          evt.graphic.setSymbol(myLabel);          //console.log(evt);               });      Thanks for any help that can be provided!


How can I make use of the $scale function in QGIS data defined properties?

In QGIS 2.0 I can use with success the $scale function in the labeling expressions. When I tried to use it for symbolizing in the data defined properties it has a value of 0. Is there a way to make it perform like it does in labeling (getting the scale)? As I saw in the expression string builder help that               $scale function     Returns the current scale of the map canvas.      Note: This function is only available in some contexts and will be 0 otherwise.      is there a source where I can read more about this e.g. in what contexts it is available? Thanks, Tudor


Are there sources for obtaining pre-made KML data from Google?

I have been asked to investigate methods for obtaining KML files. All I know so far is that it is a format based on XML that Google uses to format their maps (or the metadata). However, are there places or methods that I may avail myself of to download KML data without going through Google Earth? For example, is it possible to get the KML for the vector data of a state in the US (roads, boundaries, etc.) ? To restate it simply: I want Google's vector data so I can plop it into ArcMap. Are there repositories for such things? I doubt it, since Google never struck me as being all handy dandy open source randy, but one can hope.


Torque cumulative creates grid

How can I get Torque's cumulative setting to draw the data points in their correct locations versus artificial grids? http://cdb.io/OAgzRD My dataset of global locations looks fine under non-cumulative setting, but cumulative wants to add the locations in a grid that looks very artificial. C


Toggle PopUp Functionality

Firstly: I've look through every example on Stack GIS.Stackexchange which might deal with my problem. Of course there was good input, but unfortunately I'm not able to fix it. Admittedly my Java Script skills are not the best, so "easy" thing are also very hard to overcome. Until now I've stuck at the following level: http://jsfiddle.net/84vvR/ I'm using nliu's TOC http://gmaps-utility- gis.googlecode.com/svn/tags/agsjs/latest/examples/toc.html and the editor and popup functionalities overed by esri's API. As you can see, theoretically all work. The mapapp contains an editor an the option for popup. The problem is: If the user wants to add new features by clicking on the map, the popup pops out and gives information about the underlying area, which has the defined popup. The same problem, if the user want to use the measurement tool. On the first click, the popup pops out and so it's not possible to add a starting point for measurement. I've tried to implement a workaround with dojo.connect and a clickHandler or onClick like (button dojoType=”dijit.form.Button” onClick=”....(),...();”>Disable Popup Function I assume that it's not too hard for an advanced developer to solve this problem to bind the Popup Functionality to a checkbox or Button. Would be awesome to get some help. Thanks a lot!


Why is the Edit Vertices button (Draw Toolbar) enabled for some annotations but disabled for others?

I keep running into this:   * Export Labels in Layer A & B to annotations (In Map, In current extent only)   * Enable Focus Data frame   * Click on annotation exported from Layer A    * **Click the Edit Vertices button - No problem**   * Click on annotation exported from Layer B   * **Edit Vertices button is DISABLED** What's the problem? When I export the same labels to geodatabase annotations, I can easily edit the vertices on any annotation. Figure 1: Click on Layer A annotation - Edit Vertices button enabled ![enter image description here](http://i.stack.imgur.com/Gcibc.png) Figure 2: Click on Layer B annotation - Edit Vertices button DISABLED ![enter image description here](http://i.stack.imgur.com/qlQb3.png)


Rotate polygon preserving topology in QGIS?

Is it possible in QGIS? What I need is: 1) Rotate a polygon inside another one (bigger), obviously the polygons are not overlapped. 2) Rotate a given polygon in a topologically correct layer (like a cadastral parcels layer) What I would like to avoid is make holes, merge polygon pieces, etc. (the "old way"). Thanks in advance.


How to find the locations of the highest values in a raster?

Using ArcGIS 10, I have a raster where I would like to find the pixel with the maximum value in the raster and return its location (center of pixel) in decimal degrees. I would like to iterate through this process returning the location of the second highest value of the raster, then third, and so on so that in the end I have a list of N locations that have the highest values in the raster in order. I imagine that this might be most easily done using a Python script but I'm open to other ideas if there is a better way. Thanks!


Set Numeric Format for Class Breaks Labels?

Why doesn't this code change the numeric format of the labels generated when I build this renderer? When I execute the set_Break() this nnumeric format is correctly set in the renderer but I still get numbers that look like "110.400000" in the generated label.               INumericFormat numericFormat = new NumericFormatClass();     numericFormat.AlignmentOption = esriNumericAlignmentEnum.esriAlignLeft;     numericFormat.AlignmentWidth = 0;     numericFormat.RoundingOption = esriRoundingOptionEnum.esriRoundNumberOfDecimals;     numericFormat.RoundingValue = 2;     numericFormat.ShowPlusSign = true;     numericFormat.UseSeparator = true;     numericFormat.ZeroPad = false;     IClassBreaksUIProperties classBreaksUIProperties =     classBreaksRenderer as IClassBreaksUIProperties;     classBreaksUIProperties.NumberFormat = numericFormat as INumberFormat;     classBreaksRenderer = classBreaksUIProperties as IClassBreaksRenderer;


Packaging failed in ArcGIS Server 10.2.1

I use ArcGIS Server 10.2.1. on Windows 8 64 bit. ![Packaging failed](http://i.stack.imgur.com/P5bve.png) I'm getting the simple error: "Packaging failed" for each new service I try to publish. Any thoughts what might be causing this? Thanks


How to find network distances between two sets of origins and destinations?

I'm learning how to use Network Analyst, and I'm not entirely sure what to do (or search for). My objective is to measure the nearest distance from the red circles to the blue circles along a road network: ![enter image description here](http://i.stack.imgur.com/9nNjB.png) I have made a network in ArcCatalog using the network wizard and I'm just interested in the simplest network distance. My network looks like this: ![enter image description here](http://i.stack.imgur.com/h70Bt.png) Is this calculation possible to do? My red points are generated as centerpoints of a grid as I am just interested in the average distance for each grid. I realize that these red points are not on the network, so I'm not sure if that is a problem for the calculation. What should be my next steps? (after reading more of the help documentation!)


What is the best way to do reverse geocoding with PostGIS?

Is there a way to do reverse geocoding with PostGIS? There's a lot of information on the Web but I found it very confusing... I have imported osm data into my own Postgres with PostGIS database. Now I want to make geocoding and reverse geocoding with my database. I don't know if there are some Stored Procedures, or a project somewhere that takes care of my issue. The next step that I want to do is make a webservice that can handle this kind of requests, but that is another question. Thanks in advance.


Is possible to Unit Testing using ArcGISLocalFeatureLayer ?

DEAR PEOPLE FROM THE FUTURE: Here's what we've figured out so far… I am developing an application using ArcGIS Runtime for WPF, I want to implement unit test for the functionality in several classes, but some of the workflow includes the use of ArcGISLocalFeatureLayer Class. Trying to analyze the process I came to the conclusion that use this class is going to be an integration test rather than a unit test, due the interaction with the local ArcGIS server that is going to be created. Someone has been able to mock this kind of objects? (please...tell me that you use Rhino Mocks).


Finding parcel polygon features with frontage on multiple streets

The Department of Public Works would like to know how many properties in the city have a sidewalk at both the front and back. To get at this, I'm looking for a method to find parcels with frontage on more than one road. I'm working in ArcInfo 9.3.1 (file geodatabase). I'm looking for a solution that doesn't require scripting or downloading other tools. Street centerlines and parcel polygons are available; there is no sidewalk layer. Using some of the ideas from this forum, I have developed and tested a process. Some of the results are puzzling (corner lots may or may not be included in the final results-I think DPW doesn't want corner lots included) and I'm concerned about scaling it up to the full 42000 parcels (I tested 1223 parcels). Are there steps to add to this method to improve the results? Is there another way to accomplish the same thing that could be more automated through model builder? • Run Feature to Line tool on parcel polygons feature class (two lines are created for boundaries shared by polygons and lines retain parcel id) • Build a geodatabase topology with this newly created parcel line feature class using the rule “must not have overlaps” then validate topology (interior lines will be in error) • In ArcMap, add topology and feature class and open error inspector then search for errors > select all errors > right click > select features > delete (to delete interior lines) • Run Frequency on remaining features in parcel line feature class with parcel id as frequency field – select records in table where Frequency > 1 and export to new table • Join Frequency>1 table to parcel polygons based on parcel id, keeping only matching records – export join results


How to close and fill Gulfs and inlets in QGIS or Grass?

I have a raster layer from which I extracted only zones with a certain color, then polygonized it with QGis. I just generalized the new vector layer in GRASS to smooth things a bit. There are many holes, "gulfs", "inlets" and cuts because the color zone in the raster overlays a black and white cartographic map. I would like to fill this spots, but I can't find a easy way. I already closed valid holes, but with "gulfs" and "inlets" I cannot use the same command. ![expample of cuts and inlets](http://i.stack.imgur.com/HuCvz.png) In that image, only the hole on the right is suposed to exist. Is there a easy way to solve this?


Why my annotation markers are printed smaller?

I annotated one of my maps, for every annotation I set an svg marker whose dimension is either 6 or 8 mm. Everything looks fine in the main window. The problem arises when I try to print my map using the composer, it somehow scales down my markers to a miniscule dimension, not respecting the dimension in mm. This process is glitchy as you can see in the following screenshot: ![composer wrong dimension screenshot](http://i.stack.imgur.com/DZcXN.png) You can notice this particular svg symbol is not drawn correctly: a small portion of it (top left) seems scaled up; this single issue doesn't appear with the other svg markers in my map. EDIT: I just verified that this problem is specific to svg markers, other kinds of indicators are printed just fine


GeoServer layer preview shows only missing images for tiles

Just getting started with GeoServer and already running into problems. When I click on Layer Preview and view any of the preloaded data in OpenLayers (or any other format) it is just blank and has the missing image icon where the data should be. One interesting thing is that I could click randomly in the blank map and see attributes of selected polygons appear so I think the data is partially loaded. Example screenshot. I've installed it on multiple machines and I am getting the same exact problem. I have connected my own data through PostGIS but I got the same result. That is why I am assuming it is user error. Anyone have any ideas?


Multiple field vbscript query

I'm fairly new to arc gis and I'm trying to string along several queries. individually they work in same query they don't. background I want to query all permit types rc, re, re/ap, and db; less than 30 days old, and only in certain counties. I'll be making 3 layers with ranging dates <30, 30-60, 60-90, >90 whats working Layer 1: "PERMITTYPE" ='RC' OR "PERMITTYPE" = 'RE' OR "PERMITTYPE" = 'RE/AP' OR "PERMITTYPE" = 'DB' AND "DATEFILED" >= CURRENT_DATE - 30 ** This returns about 60 results in the below counties Layer 2: "COUNTY" IN ( 'BORDEN (TX)', 'GARZA (TX)', 'LYNN (TX)', 'KENT (TX)', 'SCURRY (TX)' ) Layer 3 (looking for different permit type but works) "PERMITTYPE" = 'OD' AND DATEFILED >= CURRENT_DATE - 30 AND "COUNTY" IN ( 'BORDEN (TX)', 'GARZA (TX)', 'LYNN (TX)', 'KENT (TX)', 'SCURRY (TX)' ) New layer (not working ) "PERMITTYPE" ='RC' OR "PERMITTYPE" = 'RE' OR "PERMITTYPE" = 'RE/AP' OR "PERMITTYPE" = 'DB' AND "DATEFILED" >= CURRENT_DATE - 30 AND "COUNTY" IN ( 'BORDEN (TX)', 'GARZA (TX)', 'LYNN (TX)', 'KENT (TX)', 'SCURRY (TX)' ) any help would be awesome thanks in advance


Workflow with raster in Qgis 2.2

I'm trying to build a workflow with the Qgis 2.2 processing toolbox. I want to clip a raster layer with a vector layer using the GDAL tool "clip with a layer mask". But it doesn't work. The error message in the historic-info is : "GDAL execution console output ERROR 4: `0' does not exist in the file system". The error message mentionned in the historic-algorithm is : "processing.runalg("modeler:model_clipp_raster","D:\Mes Documents\Documents\Qgis processing\Test_LaFleche\la_fleche.shp","D:\Mes Documents\Documents\Qgis processing\Test_LaFleche\72_SCAN_25\SCAN_25_Topo\Extrait_72_Topo.tif",None)". I've done the treatement in the Qgis "classical" interface without problem. It seems to me that problem stand in the fact that the raster can't be loaded. I hope i'm clear enough. If anyone has an idea about the source of my problem and the way to fix it, i would be really grateful. Best regards, Emily.


Problems running model in QGIS

I am trying to make a batch model to convert band 2 and band 1 MODIS 8 day composites into NDVI. The input files are in tif format, but I want the output to be saved in Idrisi .rst format. I am using QGIS 2.2 64 bit on a Windows 8.1 machine. I installed via the OSgeo4W, and installed everything, so I am using the Orfeo-toolbox for the raster calculator. I am new to QGIS and not very experienced in python either. The simple model is presented below. Ultimately, I want to batch this, but for the time being I’m just testing to see if I can get a simple model to work. I have tested the data, using the raster calculator directly in QGIS to make an NDVI and this worked without any problems. I just have a lot of images I need do this for.... ![screen capture of mode in QGIS][1] Below I have copied the OTB execution command and the execution console output from QGIS. My questions are: What does it mean that I am missing mandatory parameter “il 0”? Is there a better way to make NDVI images than using the model I have made? Any help or advice on this matter is greatly appreciated! Thanks, Brian ++++++++++++++++++++ OTB execution command C:/OSGEO4~1\bin\otbcli_BandMath -ram 128 -exp (NIR - RED) / (NIR + RED) -out "C:/8day_NDVI/NDVI_1.rst" +++++++++++++++++++++++++ OTB execution console output ERROR: Missing mandatory parameter: il 0 ERROR: Troubles loading parameter, please check your line argument... This is the BandMath application, version 4.0.0 Perform a mathematical operation on monoband images Complete documentation: http://www.orfeo- toolbox.org/Applications/BandMath.html Parameters: -progress Report progress  MISSING -il Input image list (mandatory) -out [pixel] Output Image [pixel=uint8/uint16/int16/uint32/int32/float/double](default value is float) (mandatory) -ram Available RAM (Mb) (optional, off by default, default value is 128) -exp Expression (mandatory) -inxml Load otb application from xml file (optional, off by default)


arcpy open up a script tool window once user save edits

On our map, after user draw a polyline on map and click "Save Edits", we need user to input the attributes related to the polyline drawn. In case user forget or skip this step (input attribute related), I intent to develop a function to open the "Create New Object" window (shows below) automatically, which means once user click "Save Edits" the "Create New Object" window will be open. The "Create New Object" window is script tool in toolbox. It should be auto triggered in the ArcGIS editor, my question is must I write my own editor, or is there a way to use the existing editor provided by ArcGIS? I am using ArcGIS 10.2. ![enter image description here](http://i.stack.imgur.com/oeebo.jpg)


python script in toolbox add a button on the window

I am getting stuck when I want to create a new python script in toolbox. The function I want to develop is to let user to input the user ID, then retrieve from the database the other values to allow user to do the updating. In this case, I want add a "Retrieve" button in the script to allow user to get their details. Please see the following picture of what I want to give. Is it possible? ![enter image description here](http://i.stack.imgur.com/d03kR.png)


FME - Raster Algebra and Multicriteria Analysis

I'm trying to use FME to do some raster algebra for some Multi-Criteria Analysis. At the moment, I have three rasters, all with the same extents and resolution: Slope, Hydro, and Infrastructure Mask. In the future, I will be adding additional layers (Vegetation, for example). The slope raster is a classified raster representing low to high slope. The categories are:   * 5: Less than 5 degrees (Preferrable)    * 2: 5-10 degrees    * 1: 10-15 degrees   * 0: Greater than 15 degrees (undesirable - No Go) The Hydro raster is a classified raster of buffer distances from waterways. The categories are:   * 5: Within 50m of a waterway   * 4: Between 50-100m of a waterway   * 3: Between 100-500m of a waterway   * 2: Between 500-1000m of a waterway   * 1: Greater than 1000m of a waterway The infrastructure mask has a value of 0 wherever there is some kind of infrastructure (Road, Pipeline, Powerline, Well, etc) and a value of 1 everywhere else. I would like to combine the Slope and Hydro (and eventually other layers) using FME. Ideally, I would like to add the values together, except, if the slope has a value of 0, the resulting raster should be zero. For the moment, I am using a RasterExpressionEvaluator Transformer to combine the layers, but I'm not sure of the right syntax to get what I want. ![enter image description here](http://i.stack.imgur.com/z84kg.png) This is the expression I'm using:               if(A[0]==0,0,($(Hydro_Weight)*B[0]+$(Slope_Weight)*A[0]))      A is the Slope Raster and B is the Hydro Raster. Hydro_Weight and Slope_Weight are parameters that the user enters before running the workspace. Essentially, what I'm doing is: if the slope is 0, then the result should be 0, otherwise, add the weighted values of the slope and hydro rasters. This seems like the correct way to go but all the cells in my output have a value of -128. What am I doing wrong? Is this the correct way to use the RasterExpressionEvaluator?


Attribute table in map composer is limited to 99 rows. How to change it?

I am using Qgis 2.2. I am about to create a cadastral index map, which are showing the parcels in specific gridcells, also there is an attribute table connected controlled by an atlas function. So by switching from gridcell to gridcell the attribute table automatically changes. Every gridcell has more than 99 parcels (Sometimes up to a few 100) BUT: The attribute table in the map composer is limited to 99 rows, so not all parcels can be shown in the attribute table. How can I increase the limitation of 99 rows? Thank you very much for your help.


Split line at a point with ArcGIS 10.1 Basic level license?

I have a point and polyline shapefile and want to split lines using points -- want to emulate the precise functionality of Split Line at Point tool (Data Management). I have ArcGIS 10.1 Basic level license and don't want to use third party tools -- want to achieve this by using Basic level license geoprocessing tools available within ArcGIS 10.1. Any suggestions would be very much appreciative.


Using Python in ArcGIS 10, create a line between the middle of the polylines which are intersected

I have got a river network as a feature class in a gdb. I need to calculate the topologic attributes via the connectivity graph of the river network. A graph-theoretic representation of river networks uses rivers as vertices and river junctions as links of a connectivity graph to understand the topology of river networks (takes rivers as nodes and river intersections as links of a connectivity graph). It will give information about centrality of lines and also simplify the network. In order to do this, I want to draw lines from the middle of each river branch to the middle of intersected higher order river. By this way, I can calculate the topological centrality of a river branch by using from-to node. When doing this manually, I draw a polyline from the middle of a river line to the middle of intersected higher order river line. Please help me to do this automatically by using Python in ArcGIS 10. If I do it manually it will take too long because of the drainage density.


How to resolve an Invalid Topology [Incomplete void poly] error from a geoprocessing service?

I am consuming a geoprocessing service from my Flex application. My geoprocessing service takes the input geometry, and intersects another layer. It then runs a Summery stats on the intersection output and returns the output. When I do a query or identify, and send the resulting geometry to the GP service, it runs fine. But when I Keep it in memory and do other things, which take a lot of time, and then send the geometry to the GP service, it errors out, with this error. The only difference I find in the later geometry, is that the rings are not in the same order as in the geometry I get from the server. What can I do to preserve the order of the rings?


Creating a Google Maps GroundOverlay with UTM NAD83 data

**Edit4:** The solution I posted from Jenner at VBforums was, obviously, in VB, and I used an online converter to port it to C#. Something was lost in translations, and that's why it was 10 miles off. I'm guessing I am just misunderstanding what the algorithms at the Proj.Net discussion boards are for and that's why they weren't doing what I wanted them to. I will attempt to close this question out in two days when I can mark my answer. * * * let me preface this by saying I am by no means a GIS expert. I need to overlay a png image over a Google map. Looking at the v3 api reference, the best way to do this is to use a GroundOverlay, which takes a lat,lng bounds. I am given one UTM NAD83 coordinate - the NW corner. I've found UTM to lat,long conversion algorithms online, but they don't even appear to be very accurate. The lat,lng bounds object in the Google API takes two lat,lng points - NE and SW. (Edit: I found a decent way to get the NE and SW points based on the NW point and the size of the image and the provided meters/pixel, so it appears my only problem is the accuracy of the conversion, if you would like to see the solution I used to get the NE,SW points, look here) (Edit2: I'm using .NET and just found Proj.Net; currently looking for good examples for using it). **Edit3, the current problem** : I am having issues converting the UTM to lat and long. A solution I found here a while ago gives me a point about 10 miles south of where I would expect it. A simple mathematical solution I found here, is giving me unexpected results.               double[] inverseMercator (double x, double y) {          double lon = (x / 20037508.34) * 180;          double lat = (y / 20037508.34) * 180;               lat = 180/Math.PI * (2 * Math.Atan(Math.Exp(lat * Math.PI / 180)) - Math.PI / 2);          return new double[] {lon, lat};     }     double[] toPoint = inverseMercator (686029.702258, 3581213.621173);      I get the following result: long: 6.1627096689832594 lat: 30.602349476368449 I get similar results using the Proj.Net solution provided by D_Guidi in the same thread. Using an online converter, I was able to get something closer to what I am expecting: lat: 32.35238307052292 long: -91.0230710652583 Can anyone shed any light on what I am doing wrong? I'd prefer to stay away from third party software and just find an algorithm to help me accurately overlay the image with the data I'm given, an example of which I've included below:               2.627469     0.000000     0.000000     -2.627469     686029.702258     3581213.621173     UTM     NAD83     15      First row is meters/pixel Any help you may be able to provide is certainly appreciated.


How to add feature class from ESRI personal geodatabase (.mdb) in QGIS

How do you add a feature class in a ESRI personal gdb in qgis. I understand this can be done through the option add vector layer and select geodatabase. But what after that ? It doesn't allow me to browse through .


Creating nanoscale DEM with GDAL

A bit of a strange question perhaps, but let me give you a breif explanation of the background before my actual questions: Atomic force microscopy (AFM) is a method that, in short (and to my limited knowledge) lets researchers scan areas at the micro- and nanoscale. It works by "scanning" an area using a probe of sorts. More that is difficult for me to explain, as I have no real understanding of it. What I do know, and what triggered my curiosity was that the result is in fact a "grid" of "height" values (a matrix of say 512x512 values describing the height of the probe at that point). I then thought: Well, apart from the scale, this is in fact a digital elevation model! And, this means that if I can manage to create a DEM file as understood by GIS tools I could apply GIS analysis to it! As it is, my significat other works at a lab that has a AFM-machine, and is using it in one of her projects. I got some scan files from her, and have managed, using Python (struct and numpy), to parse these binary files and what I have now is a numpy array of size 512x512 filled with int16-values. What I'm planning on next, and what I need some help with is the "mapping to a proper DEM"-part. I have some knowledge about DEMS, but when it comes to actual generation of them I am quite new. What I'm thinking is that I have to georeference my data in some way, and for this I need a custom (planar) coordinate system. I envision my coordinate system would use micro- or nano-meters as units. Then it's just a matter of finding the size of the area scanned with the AFM (this I believe is somewhere in the binary file, assume this known). **update** : I also have several scans at different resolutions, but of the same area. For example i have this information about two scans: larger image:               Scan Size: 51443.5 nm     X Offset: 0 nm     Y Offset: 0 nm      smaller (detail) image:               Scan Size: 5907.44 nm     X Offset: 8776.47 nm     Y Offset: 1486.78 nm      What I', thinking is that my custom coordinate system should have an origin in 0,0 and for the larger image I cam assign pixel 0,0 the coordinate value of (0,0) and pixel 512,512 the coordinate value (51443.5, 51443.5) (Guess you get the picture for the other points needed). Then, the larger image would map pixel (0,0) to (8776.47, 1486.78) and (512,512) to (8776.47 + 5907.44, 1486.78 + 5907.44) **1st Question is then** : How do I create a proj4 def for such a coordinate system? I.e: how do I assign these "real world coordinates" to my custom coordinate system (or, if I follow whubers suggestion and using a local coordinate system and lying about units (i.e. treating my nanometers as kilometers) Then I have to transfer my numpy 2-dimensional Array to a Georeferenced DEM file format. I was thinking of using GDAL (or, rather the Python bindings). **2nd Question is then** : How do I create a georeferenced DEM from "arbitrary" data such as mine? Preferably in Python and using open source libraries. The rest should then be fairly easy going, just a matter of using the right analysis tools. The problem is, this task is driven by my own curiosity, so I'm not quite sure what I actually should do with a nanoscale DEM. This begs the **3rd question** : What to do with a nanoscale DEM? What kind of analysis can be done, what are appropriate tools for DEM analysis and finally: is it feasible to make a map with hillshade and contour lines from this data?:) I welcome all suggestions and pointers, but keep in mind that I am looking for free alternativies, as this is a strictly hobby-based project with no budget or funding (and I do not have access to any lisenced GIS-applications). In addition, I do know that Bruker, the company that sells these AFM machines, ships some software, but using that would not be any fun.


PostgreSQL, PostGIS and GeoServer publishing 1.7 million records

I have a polygon geodata in postgresql and I am publishing out to geoserver the 1.7 million records. How long is too long when processing? EDIT: Forgot to mention that I successfully published a database spatial view subset of this same data. Side question, What is postgresql really doing during this process??? @Devdatta Tengshe, I have a polygon dataset with about 1.7 million records and I am just trying to publish the data to geoserver. I published a spatial view from a small table and a join to the large polygon dataset. This layer seems fine in openlayers. I just published another spatial view from a different small table joined to the large polygon dataset and it published but will not finish loading in openlayers.


unable to install qgis 1.9 on Lubuntu 13.04

I tried to install qgis 1.9 in my laptop. I use Lubuntu 13.04 (64bit). I used the repository               deb     http://qgis.org/debian-nightly raring main     deb-src http://qgis.org/debian-nightly raring main      and I didn't find any error during the installation. When I try to start qgis I obtain this message: > /usr/bin/qgis.bin: error while loading shared libraries: libgdal1.7.0.so.1: > cannot open shared object file: No such file or directory If I try               sudo apt-get update     Sudo apt-get upgrade       I obtain those messages:               Configurazione di qgis-providers-common (1.9.0+git20130529+e0406e7~raring1)...          /usr/lib/qgis/crssync: error while loading shared libraries: libgdal1.7.0.so.1: cannot open shared object file: No such file or directory     dpkg: errore nell'elaborare qgis-providers-common (--configure):      il sottoprocesso installato script di post-installation ha restituito lo stato di errore 127     dpkg: problemi con le dipendenze impediscono la configurazione di qgis-providers:      qgis-providers dipende da qgis-providers-common (= 1.9.0+git20130529+e0406e7~raring1); comunque:       Il pacchetto qgis-providers-common non è ancora configurato.          dpkg: errore nell'elaborare qgis-providers (--configure):      problemi con le dipendenze - lasciato non configurato     dpkg: problemi con le dipendenze impediscono la configurazione di qgis:      qgis dipende da qgis-providers (= 1.9.0+git20130529+e0406e7~raring1); comunque:       Il pacchetto qgis-providers non è ancora configurato.          dpkg: errore nell'elaborare qgis (--configure):      problemi con le dipendenze - lasciato non configurato     dpkg: problemi con le dipendenze impediscono la configurazione di qgis-plugin-globe:     qgis-plugin-globe dipende da qgis (= 1.9.0+git20130529+e0406e7~raring1); comunque:     Segnalazione apport non scritta poiché il messaggio di errore indica la presenza di un fallimento precedente.          Segnalazione apport non scritta poiché è stato raggiunto il valore massimo di MaxReports            Il pacchetto qgis non è ancora configurato.          dpkg: errore nell'elaborare qgis-plugin-globe (--configure):      problemi con le dipendenze - lasciato non configurato     dpkg: problemi con le dipendenze impediscono la configurazione di qgis-plugin-grass:      qgis-plugin-grass dipende da qgis (= 1.9.0+git20130529+e0406e7~raring1); comunque:       Il pacchetto qgis non è ancora configurato.          dpkg: errore nell'elaborare qgis-plugin-grass (--configure):      problemi con le dipendenze - lasciato non configurato     Si sono verificati degli errori nell'elaborazione:      qgis-providers-common      qgis-providers      qgis      qgis-plugin-globe      qgis-plugin-grass     E: Sub-process /usr/bin/dpkg returned an error code (1)      I apologize for the Italian messages. Could you help me? Thank you


How to retrieve PROJECTION value for a map?

I am trying to read the Projection from Geotiff/shape file, using gdal 1.9.0 by the following command in the command line... gdalinfo.exe -proj4 map-path but its retuning empty string for "PROJECTION" parameter.. when I check the same map in QGIS getting the correct projection for this map. Actually I want to retrieve "PROJECTION" value for any map from the C/C++ code , could you please let me know the way to retrieve this ? Thanks, Prady.


How can I choose values from stacked raster images according to location of ROI

Before I choosed ROI on land surface temperature (LST) images and obtained pixel values of LST. Now I need pixel value of another images which are stacked multi-band image. My main purpose is to obtain pixel values of both images "at the same location". I think that I can use ROI file from LST to choose the pixel values of multiband image. Anyone who guide me to accomplish this purpose using ARCGIS, ENVI, ERDAS? I choosed the ROI in rectangular mood in ENVI. When I opened it, it looks like small polygons (not points) spreading on my study area. Thanks a lot


Google Map geocoder returns multiple results against a complete town name (with county, state)

I have the following address: "Chelsea, New York, NY, USA" In the above address: City = 'Chelsea' County = "New York" State = "NY" Country = "USA" When i pass this address to the Google's Geocoder API it returns multiple results, while i am expecting a single result because every part (name, county, state, country) of the town name is specified in the address. Here is the code snippet: (Assuming required files are already included)               var geocoder = new google.maps.Geocoder();          geocoder.geocode({ "address": "Chelsea, New York, NY, USA" }, function (results, status) {                              if (status == google.maps.GeocoderStatus.OK) {                       if (results.length == 1) {                                                  // do something                             }                           else{                              // do something                             }     }      I will appreciate your help. Thank you.


geoserver heatmap tiles doesn't display well (with leaflet)

I am trying to display a heatmap generated with geoserver on my map using the L.tileLayer.WMS class from Leaflet, but the tiles do not fit well together. ![enter image description here](http://i.stack.imgur.com/A4wOE.jpg) Do you have any idea how to fix it?


How to convert arcpy.InsertCursor to arcpy.da.InsertCursor in Python code?

I have some working code using `arcpy.InsertCursor` which I want to convert to `arcpy.da.InsertCursor`. This I want to do primarily as ESRI claims that `da` methods are more optimized. The code I have is follows:               array = arcpy.Array()     point = arcpy.Point()     tie_shp_file = arcpy.env.workspace + "/" + out_name     insert_cursor = arcpy.InsertCursor(tie_shp_file)     feat = insert_cursor.newRow()          for line in all_lines:   #all_lines is of the form array([[0, 3],[0, 7],[1, 3]) indicating line [from, to]         #from         point.X=coords_pts[str(pt[line[0]])][0]         point.Y=coords_pts[str(pt[line[0]])][1]         array.add(point)         #to         point.X=coords_pts[str(pt[line[1]])][0]         point.Y=coords_pts[str(pt[line[1]])][1]         array.add(point)         polyline=arcpy.Polyline(array)# Create a Polyline object based on the array of points         array.removeAll() #Clear the array         feat.shape = polyline         insert_cursor.insertRow(feat)      The code inserts lines to a line shapefile. The lines are created by specifying its end points. How do I convert this? I am specifically stuck on how to pass the to and from points in the new version.


Coordinates conversion doubt

In QGIS do I have a EPSG:4326 shapefile that shows the coordinates "-7148300 -3688500" and I need them in the format -64º -31º for further processing in a geodjango app. However using ogrinfo I get this information from the shapefile: > >     Feature Count: 22 >     Extent: (-64.279736, -31.509785) - (-64.092643, -31.332019) >     Layer SRS WKT: >     GEOGCS["GCS_WGS_1984", >         DATUM["WGS_1984", >             SPHEROID["WGS_84",6378137,298.257223563]], >         PRIMEM["Greenwich",0], >         UNIT["Degree",0.017453292519943295]] >     id: Integer (10.0) >   How can I convert the values? I found useful both solutions: 1) Getting the coordinates of a given feature id, within Geodjango shell, so I have a reference and work from there using a custom conversion function. 2) Know the mathematical equivalence between the two pairs of coordinates. A snippet or idea in geodjango would be nice, but any hint will be appreciated! Thank you in advance.


How to produce a hazard map using slope and geology information

for my project I need to display a landscape's susceptibility to erosion. I created a slope-map out of my DEM and added the bedrock geology as a polygon- feature shapefile. Now I want to tell ArcGIS that areas with a certain slope (lets say >30°) and a soft rock underlying (e.g. schists) are particularly affected by erosion, whereas plain areas with harder lithologies are not. I used the Raster Calculator tool before to produce hazard maps, but I never used a shapefile as input. How can I transform the information stored in my shapefile into a raster? Thanks for helping!


Counting points in a raster

Is there a way to select points in a raster dataset by raster classification? I have a floating point raster created from the points, classified into 4 ranges by standard deviation. I want to find out how many points are in each range. I thought of converting the raster to polygons, but that requires an integer raster. Then I thought ok, convert the floating point raster to integer, but that requires Spatial Analyst, which someone else has checked out. Maybe there's something in Spatial Statistics, but not familiar with those tools. If the raster was in polygon format, a simple Select by Location would work. Thanks. ArcGIS 10.0 sp5. Windows 7 Enterprise


Leaflet and EPSG 4326 Projections

**UPDATE: LEAFLET fixed support for EPSG 4326 (LL84) in the latest betas:https://github.com/Leaflet/Leaflet/issues/1207** Leaflet says it supports EPSG 4326 (Plate Carree, aka ll84) projections. I think it does, but take a look at these two jsfiddles: http://jsfiddle.net/drysg/cBzVc/ - This is the original 3857 web Mercator. As you can see the overlay of both the tiles and the red outline are properly placed near Kraknow. (click the mouse on the map and you get the lat/long as a popup). Ok, now look at: http://jsfiddle.net/drysg/XPP54/ - This is in EPSG 4326 (ll84). The background is a MapCarta ll84 slippy tiles. It seems that leaflet does not know about the 2x1 aspect ratio that OSG wants. That is why the background wraps (http://wiki.osgeo.org/wiki/Tile_Map_Service_Specification#global- geodetic) Well, if they want to support that fine. I am actually happier with with a 1x1 map (I can generate that with SAFE FME). So ignore that for a while. Switch the background to the CloudMade OSM map (ingore that the origin is at 85N, so that the map is 5 degrees shifted to the south). Now everything looks correct, to me. The only thing I need now is to figure out the formula for the SLIPPY tile numbering. The web page at http://wiki.openstreetmap.org/wiki/Slippy_map_tilenames Shows how to do this for Web Mercator. What would the formula be for LL84, where the map bounds are not (-85, -180; 85, 180) but rather (-90, -180; 90, 180)


How to do change detection between two Landsat scenes?

I want to do change detection and calculate the changes step by step using two images: the first image Landsat 1986, and the second Landsat 2000. Actually I have used image differences in ArcMap and I get raster values for RGB ranging from 213 to -69. What I should do to make this result more meaningful?


precision of positionalongline and/or SpatialJoin

I've been using positionAlongLine to create points that are supposed to line up with the lines from which they are derived. When I tried to use `SpatialJoin_analysis` with the `match_argument` set to "INTERSECT" , I'm not getting all of my points represented in the output point set. After zooming in ridiculously far on an example of one of my errant positionAlongLine-generated points and the original linework, I can see that the point is 0.000025 meters away from the line that it is based on. I really just want my `SpatialJoin_analysis` to work with "INTERSECT". I've tried specifying a `search_radius` of 1 with no improvement in results. Short of continuing to use "CLOSEST" as a workaround, which is a bit cheesy, I'm wondering if there is a way to improve the precision of the coordinates I get with `positionAlongLine` or the `NumPyArrayToFeatureClass` I use to bring the coords back into ArcGIS to create the points featureclass. I'm also happy to do something to make the `SpatialJoin_analysis` handle the gap, too. ---Update 1--- Tried @radouxju's suggestion. Can't seem to get `NumpyArrayToFeatureClass` to create anything w/in a FD. Results just appear at top of fgdb (w/no warnings, which seems odd). Copying the source FC to the top of the fgdb and then doing the analysis doesn't remove the shift. Going to isolate the point geometry in my numpy array to see if that's the source of the error. Guessing the solution will not be to fix this, but to improve the SpatialJoin.... ---Update 2--- Checked the coordinates in my numpy array and they do fall on top of the line (at least as close as I can visually verify), but `NumpyArrayToFeatureClass` is not honoring the full precision (dropping decimal places) in the points of the featureclass it creates. ![enter image description here](http://i.stack.imgur.com/eBdbm.png)


what are the all steps to find catchment area for weir by arcmap

what are the all steps to find catchment area for weir on intermittent stream by arc map. I have coordinates and elevation of that point. I am in search of that but I did not find it clearly. i have to find I have not total command on arc map so give me detailed answer. thank you.


how to set workspace of Geoprocessor in .net?

I am migrating some code from VBA to C# for developing an addin. In VBA, the geoprocessor's workspace could easily be accessed by workspace property of Geoprocessor. But how to set the geoprocessor's Workspace property in .NET environment


How do I access QGIS program settings programmatically?

I am trying to temporarily enable the option in QGIS using PyQt4 to suppress the attribute dialog popup after a feature is added. I found an example in the cookbook but it doesn't indicate how to get to existing settings.


Spatial Data from Postgis to Java/REST with RestEasy + Jackson

I'm developing an application where I take a route from google maps, save it in postgresql/postgis. I'd like to take this route (LineString geometry) and send it back to Google Maps. I already can save the map in the database, now I want to show the routes on a map. On the client side I'm using JQuery (AJAX) + Google Maps Api. On the server side I'm using REST services with RestEasy. I tried to work with Jackson to transform the Entity to JSON, but, the Jackson fails to transform geometric data in json. Does anyone know how to transform this spatial data into JSON? **My Entity:**               public class Rota {     ...     private LineString caminho;     public Rota() {     }     public Rota(Integer id, Usuario usuario, String descricao, LineString     lineString) {         this.id = id;         this.usuario = usuario;         this.descricao = descricao;         this.caminho = lineString;     }     // getters and setters     @Type(type = "org.hibernate.spatial.GeometryType")     public LineString getCaminho() {         return caminho;     }     public void setCaminho(LineString caminho) {         this.caminho = caminho;     }     // hashcode and equals method     }      **My REST service:**               @Inject     private RotaBC rotaBC;          @GET     @Path("/rotas/{idrota}")     @Produces("application/json")     public Rota show(@PathParam("id") String id) {         return rotaBC.show(Integer.parseInt(id));     }       **The error:** > 15:05:16,447 ERROR > [org.apache.catalina.core.ContainerBase.[jboss.web].[default- > host].[/sharecar].[default]] (http--0.0.0.0-8080-1) Servlet.service() for > servlet default threw exception: java.lang.ClassNotFoundException: > org.codehaus.jackson.map.JsonMappingException$Reference from [Module > "org.codehaus.jackson.jackson-mapper-asl:main" from local module loader > @1d332b (roots: /opt/demoiselle/server/jboss-7.1/modules)] ... at > org.jboss.modules.ConcurrentClassLoader.loadClass(ConcurrentClassLoader.java:120) > at > org.codehaus.jackson.map.ser.std.BeanSerializerBase.serializeFields(BeanSerializerBase.java:166) > ... at > org.codehaus.jackson.map.ser.BeanPropertyWriter.serializeAsField(BeanPropertyWriter.java:446) > at > org.codehaus.jackson.map.ser.std.BeanSerializerBase.serializeFields(BeanSerializerBase.java:150) > at > org.codehaus.jackson.map.ser.BeanSerializer.serialize(BeanSerializer.java:112)


When a raster is described in general as either a double or an integer how do conversions work for that type

I know a raster format is essentially a grid of cell values that may have values associated with color, but I'm unsure of how the primitive types factor into the format. In addition if I need to write a conversion function for a 'double' raster to an integer raster what is that referring to as in how does the logic work? Is it iterating each cell and simply casting a value? I am using the sextante library to write this java function. I assume the conversion will essentially be one raster converted and then used with a grid reader of some kind.


Is it possible to load large non georeferenced pictures in Geoserver?

I use Openlayers as a viewer of very large pictures, they are tiled and the navigation is exactly what I need. I also can add a vector layer and draw on it. But when it comes to save those drawings, saving them online... I can not find anything. Then I learned about WFS-T and if a map (my picture) and its vector layer are on Geoserver it is possible for the user to draw on it and save the modifications to Geoserver. Now I can not see how to upload a non georeferenced picture on Geoserver ? I got it up and running but that's it. At this point I'm such a newbie I don't even know if my combo is the good solution to my problem : drawing with a vector layer on very large pictures and save online the modifications.


Create Custom Popup for Locate Widget

Without touching the FlexBuilder app (I'm still going through the tutorials for Flex) is there a way to customize the pop-up that appears for the Locate Widget? I am assuming that I could go into the pop-ups folder and create a popup_Locate config file of some sort? Is this correct? Or is there a better way to customize the popups that appear when selecting from the list that populates from a Locate Widget search?


Need to convert from ISpatialReference to SpatialReference

I am attempting to write a method in C# that will accept as input, an ESRI.ArcGIS.Geometry.ISpatialReference object and convert it to an ESRI.ArcGIS.Client.Geometry.SpatialReference object. I need this for converting an IPoint to a MapPoint. Are there any online resources for this? Thanks!


Get Index for a group QGIS Python

I am trying to add new layers into a specific group which I can get to work if I know the index but how can I get python to fid the index of a group so it can load the new layer into the group wherever it is in the TOC ?               qgis.utils.iface.legendInterface().moveLayer( layer, index )


Are multiple separate feature classes more efficient than repeatedly deleting all features from the same feature class?

_Quick Version: Is it more efficient, in terms of database performance, to have multiple separate temporary feature classes or one feature class that you repeatedly empty out and fill?_ I am working with ArcGIS 10.1, Windows 7, using personal geodatabases. I have a python toolbox tool that selects survey cells based on their proximity to a seed-cell. Once a boundary is determined, based on distance from the last selected set of cells, the final selected set is buffered to form the boundary. I have been using a single feature class ( _BoundaryHopper_ )to hold the output of the buffer and then appending that output to a second feature class that holds all the boundaries for a project ( _DerivedBoundaries_ ). _BoundaryHopper_ would be emptied out at the beginning of each loop iteration to avoid repeatedly appending the same boundaries to _DerivedBoundaries_ ] Recently I switched the tool over to using separate outputs for each buffer output ( _TempBoundary1,2,3..._ ) and appending each individual feature class to _DerivedBoundaries_. That solved some issues I was having from trying to overwrite the same feature class all the time. So, my question is: **_Are there any significant performance benefits of one technique over the other?_** Does one technique create less 'bloat' in the database or result in faster processing? Thanks!


Moving a Google Map Marker on a Mac

I have exact lat & long coordinates that I will enter and hit return in Google maps. When I do this it displays this point with a green marker on the map but I can not move this marker. I want to be able to move this green marker nearby and get the new coordinates for exactly where I moved the original marker to. I only need to move it probably less than 1-2 inches on the map. I need someone to walk me through this. I have seen some answers showing java and html but I do not know how to use these to solve my problem. Help


Calculate maximum speed from GPS data

I have raw GPS locations for a track, and need to get maximum speed. Simple methods (calculate speed between 2 points, take max) will give meaningless large numbers due to GPS inaccuracy, the point is jumping around. Can you suggest a good ready-made algorithm to solve this?


Leaflet / Google Map baselayer / markers not visible

New to Leaflet, currently trying to migrate a project from native Google Maps to Leaflet; having been drawn to Leaflet by the layer management and options of using better 'roadmap' style basemaps. Have downloaded the plugin to use Google Map tiles as a base map for Leaflet,and its working great as a base layer. However I have found a problem, and I'm not sufficiently experienced with Leaflet to understand why. I have markers which are visible with any other base layer - but not with Google. Have tried forcing Z order high (1000 etc.) but it makes no difference. Opacity of the markers is full 1.0 Markers are visible on all other layers (e.g. OSM), when I switch base to Google (satellite or map) they're no longer shown, but switch back to a non- Google base and they're there again. The reason I'm pursuing the Google base layer is that Google is (as far as I'm aware) the only provider of detailed consistent aerial photographic tiles for my region of interest. Wish to use other maps but have option to switch to a Google 'satellite' view via the layer manager. Currently using Leaflet v0.4.5 ( but have also tried v0.4 ), same symptoms occur using different browsers (IE, FF etc.). Would be grateful for any assistance.


How do I repair NVIZ error: 'Can't find a usable init.tcl'?

How do I repair the following problem: ![NVIZ error](http://i.stack.imgur.com/k2Ylr.jpg) The error occurred while I was trying to run NVIZ in the GRASS environment in QGIS. NVIZ however runs perfectly in GRASS. On checking the directories listed, init.tcl is located there.


Can one hide the Lat/Lng information sent to Google maps?

most services that use google maps have the Longitude / Latitude stored in the source to be sent to google maps to render their markers. However on this site (sorry Chinese) http://www.dianping.com/shop/500451 I can't find this information for Latitude Longitude of their items. SO how does google know where they are supposed to render the marker?


Backwards compatible python scripts in ArcGIS

As I move through my small collection of python scripts written for arcgis 9.3 I'm discovering that many things work with little or no changes in arcgis 10. For example this pattern "just works":               import arcgisscripting     gp = arcgisscripting.create(9.3)     gp.workspace = 'd:/scratch'      even though the new way is:               import arcpy     from arcpy import env     env.workspace = 'd:/scratch'      So my question is, what it is the best structure to follow so that scripts are only written in the v10 way when they _need_ functionality not available in earlier versions? Or perhaps better phrased as, what is the best structure for backwards compatiblity? I don't want to make my stuff automatically unavailable to people on earlier versions, just when it wouldn't work anyway.


python: Replace Variable in VB-Expression-File (.cal) [ArcGis 9.3]

I'm trying to modify a VB-Expression saved in a cal-file:               dim out     if [D19830701]=-9999 then         out= [D19830501]     else         out= [D19830701]      end if     __esri_field_calculator_splitter__     out      with this simple script:                 import os, sys     o = open("D:\\Expressionnew.cal","a") #open for append     for line in open("D:\\Expression.cal"):         line = line.replace("D19830701","D19830801")         o.write(line + "\n")     o.close()      but all i get is this:                 d i m   o u t                 i f   [ D 1 9 8 3 0 7 0 1 ] = - 9 9 9 9   t h e n                     o u t =   [ D 1 9 8 3 0 5 0 1 ]                 e l s e                     o u t =   [ D 1 9 8 3 0 7 0 1 ]                   e n d   i f                 _ _ e s r i _ f i e l d _ c a l c u l a t o r _ s p l i t t e r _ _                 o u t       No Replacement, the script just added "spaces" between every letter.   I want to use this Replacement-Script to fit the cal-file for several calculations i need to make.      Is there maybe a problem with the encoding (unicode?) when the script opens the file?   Maybe somebody has had the same problem ;)   or a nice workaround.   Thank you for reading.... Markus


Automatically changing field entry in related table when entering it in feature class?

I am working to try and develop a code or process that I can run that will allow me to automatically change an XY field in a related table when I change an associated XY field. Let me explain our situation a little first. We have a system where we are trying to maintain utility data in an SDE, and while doing this data is being entered by technicians in the field and brought back in. When data is entered in the field, there is a FACILITY_ID associated with each item, such as hydrants, values etc. This data is pushed to each cooresponding our feature class. However, when multiple actions are taken on the same feature, we are storing each of those changes in a different table while only reflecting the last action in the feature class. These two have a relationship class. So, there is only one entry for each FACILITY_ID in a feature class, but there may be multiple entries in the 'related' table for each FACILITY_ID. Our problem is what to do if we move a hydrant/value in ArcMap and have to change the XY coordinate. We don't mind recalculating the X and Y values in the feature class table, but we would like a process that could recognize that change in the feature class XY and automatically change all of the X and Y values for the associated related table FACILITY_ID without having to go in and edit them by hand. Is there a simple way to do this? We've been looking at ArcObjects or ArcPy, but does a command like this already exist in ArcMap? I should clarify, the 'related table' is more just a table that has a simple relationship class connected to the feature class. The second table in a non- spatial one.


How do I create a functional spatial index of an ST_Geometry Feature Class in Oracle/SDE?

Copied an SDE feature class (with blob(sdebinary(?) geom) into the same SDE GDB (Oracle 11g/SDE 9.3.1) but with the ST_geometry configuration keyword. I can see that the feature class was created, all the fields come across fine, I can query it with st_geom functions etc., BUT in ArcMap, ArcCatalog, the features don't draw. (*note that when _selected_ I can see the outline of the shape just fine) When I remove the spatial index the features draw, add the spatial index back, _poof_ no draw. It would seem that ArcGIS cannot create a functional spatial index of my ST_Geometry features.


How to create polygon with DBGeography Class in Csharp

I want to create a polygon with csharp 4.5 via DBGeometry class. I found example on csharp to create a point via System.Data.Spatial namespace. (Requires Reference System.Data.Entity) This code is working :               var text="POINT(29.11, 40.11)";     DbGeography.PointFromText(text, 4326);  //I can generate a point with this code      But i want to create Polygon object and these codes dont work . Please help me about which prototype - notation should i use to create a polygon or multipolygon object ? And i can not found examples on Microsoft web site. It contains only prototypes, not examples :(               //Not Work     var polyb = DbGeography.MultiPolygonFromText("POLYGON((29.11 40.11, 29.11 40.0 , 29.30 40.0, 29.30 40.11, 29.11 40.11)", 4326); --          //Not Work     //var polyc = DbGeography.MultiPolygonFromText("POLYGON((29.11, 40.11), (29.11, 40.0), (29.30, 40.0), (29.30, 40.11), (29.11, 40.11))", 4326);--          //Not work     //var reVal3 = DbSpatialServices.Default.GeographyPolygonFromText("POLYGON(127.652000002563     -26.244400002062321, 127.652000002563 -26.194399997591972,127.652000002563 -26.244400002062321)", 4326);      Summary How can i create polygon with csharp DBGeometry class from text method ?


r.neighbors produces no output

I am running               r.neighbors --overwrite --verbose input=bf03pits@PERMANENT output=bf03pitsdenoise method=mode size=5      in the GRASS GUI using the `r.neighbors` dialog. bf03pits has 3955848 all with values either 0 or 1. When I run `r.neighbors` the result has 1 cell with a range of `-nan` to `-nan`. Yet the output shows no errors. Why might `r.neighbors` be failing? How can I go about doing what I want?


How can I fix badly specified geographic coordinate system with units in arc-seconds?

Check out this depth-to-rock soil data. It is linked from here. The spatial reference system is a latlon, but the units are arcseconds. I could not reason out how to override the SRS even after staring at this page in the Proj4 wiki for a while and trying to convert the origin coordinates using `cs2cs` with a `+towgs84` argument. Here's the `gdalinfo` output for that file:               Driver: E00GRID/Arc/Info Export E00 GRID     Files: rockdepm.e00     Size is 6936, 2984     Coordinate System is:     GEOGCS["NAD27",         DATUM["North_American_Datum_1927",             SPHEROID["Clarke 1866",6378206.4,294.978698213898,                 AUTHORITY["EPSG","7008"]],             AUTHORITY["EPSG","6267"]],         PRIMEM["Greenwich",0,             AUTHORITY["EPSG","8901"]],         UNIT["degree",0.0174532925199433,             AUTHORITY["EPSG","9108"]],         AUTHORITY["EPSG","4267"]]     Origin = (-449100.000000000000000,177840.000000000000000)     Pixel Size = (30.000000000000000,-30.000000000000000)     Corner Coordinates:     Upper Left  ( -449100.000,  177840.000) (Invalid angle,Invalid angle)     Lower Left  ( -449100.000,   88320.000) (Invalid angle,Invalid angle)     Upper Right ( -241020.000,  177840.000) (Invalid angle,Invalid angle)     Lower Right ( -241020.000,   88320.000) (Invalid angle,Invalid angle)     Center      ( -345060.000,  133080.000) (Invalid angle,Invalid angle)     Band 1 Block=6936x1 Type=Int32, ColorInterp=Undefined       Min=0.000 Max=191.000       Minimum=0.000, Maximum=191.000, Mean=123.286, StdDev=38.233       NoData Value=-2147483647      Notice that the pixel size is 30 arcseconds. If you divide a given set of coordinates by 3600 you get something that makes sense in degrees. In the course of composing this question I have discovered that the data is available with other SRSs, but in case those are poorly defined as well I thought I would post anyhow.


Southern Newfoundland and New Brunswick buildings data

I am looking for GIS data for southern Newfoundland and New Brunswick. For NFLD: - buildings, property line/polygon, basic contours, topographic items For NB: - buildings Let me know of any data sources or websites that may be of interest.


Where to find books or course materials to help prepare for Esri Technical Certification in ArcGIS Desktop and Developer?

Where I can find books or course materials to help prepare for Esri Technical Certification in ArcGIS Desktop and Developer?


Is it possible to create new polygon with the streaming mode on QGIS1.8?

I would like to create a new polygon with the mode of stream. In Arc GIS, it is possible to change the mode from stream to normal that vertec is placed based on click. Strem mode means vertex is placed automatically at an interval. I'm QGIS 1.8.0 user. Thank you in advance.


Best way to create/draw a map (especially indoor)

I'm currently developing an indoor map/room plan for an university project. I've already searched through stackexchange to figure it out but I couldn't find anything similar. So what I currently do to create a indoor map/room plan is to use QGIS and draw points, lines and rectangles. It has a nice snapping feature and other stuff that makes the map creation slightly convenient. But I have problems when it comes to alignment, squareness etc. That's why I'm searching for a more convenient solution to draw maps, especially indoor maps. What I'm **not** looking for is an app or business solution thats tracks my movement and creates a map with it. I really would like to draw the map on my own. Thank you! =)