slow import via osm2pgsql to postgresql database

I try to import berlin.osm (33Mb) into the postgres database 9.1 via osm2pgsql. My desktop is mem 2Gb,proc 2Ghz,hdd 60Gb, RHEL6u2x64. This process is so slowly in particular when I'm going over pending ways. I see this: _osm2pgsql SVN version 0.80.0 (32bit id space)_ _osm2pgsql -M -d berlin -s -S "/usr/share/osm2pgsql/default.style" -U "postgres" -W -H "localhost" -P 5432 "/opt/maps/berlin/berlin.osm.bz2"_               Processing: Node(1644k 49.8k/s) Way(304k 12.17k/s) Relation(6337 52.37/s)  parse time: 179s          Node stats: total(1644720), max(1968785138) in 33s     Way stats: total(304296), max(186177757) in 25s     Relation stats: total(6337), max(2504060) in 121s     Committing transaction for planet_osm_point     Committing transaction for planet_osm_line     Committing transaction for planet_osm_polygon     Committing transaction for planet_osm_roads          Going over pending ways          Using 1 helper-processes     processing way (10k) at 0.11k/s     processing way (195k) at 0.11k/s     All child processes exited          Pending ways took 1719s at a rate of 114.01/s     Committing transaction for planet_osm_point     Committing transaction for planet_osm_line     Committing transaction for planet_osm_polygon     Committing transaction for planet_osm_roads          Going over pending relations          Using 1 helper-processes          All child processes exited          Committing transaction for planet_osm_point     Committing transaction for planet_osm_line     Committing transaction for planet_osm_polygon     Committing transaction for planet_osm_roads     Sorting data and creating indexes for planet_osm_point     Sorting data and creating indexes for planet_osm_line     Sorting data and creating indexes for planet_osm_polygon     node cache: stored: 1644720(100.00%), storage efficiency: 66.93% (dense blocks: 60, sparse nodes: 1597241), hit rate: 100.00%     Sorting data and creating indexes for planet_osm_roads     Stopping table: planet_osm_ways     Building index on table: planet_osm_ways (fastupdate=off)     Stopping table: planet_osm_nodes     Stopping table: planet_osm_rels     Building index on table: planet_osm_rels (fastupdate=off)     Stopped table: planet_osm_nodes in 0s     Stopped table: planet_osm_rels in 1s     Indexes on  planet_osm_roads created  in 5s     Completed planet_osm_roads     Indexes on  planet_osm_line created  in 21s     Completed planet_osm_line     Indexes on  planet_osm_point created  in 26s     Completed planet_osm_point     Stopped table: planet_osm_ways in 31s     Indexes on  planet_osm_polygon created  in 31s     Completed planet_osm_polygon      _Osm2pgsql took **1931s** overall_ I used the next parameters in my postgresql.conf :               shared_buffers=512Mb     temp_buffers=16Mb     work_mem=64Mb     maintenance_work_mem=128Mb      When I checked it on the Windows desktop it acceptable worked, but if I started it on my linux system I have a problems with a import perfomance. Could someone give me advice what I should do that can fix it?


ENVI to map composition

I'd like to print out 6 images in a single page. How can I do this using map composition tools in ENVI (6 images are year-to-year dynamic images with the same size and the same geo-referencing)? Thanks all


is ellipse inside ellipse in WGS84 ellipsoid

My question is quite clear: Assume that there is an ellipse, namely EL1. Parameters of the ellipsoid are (focus11 [latitude, longitude], focus12[latitude, longitude], radius1). And a second ellipse, namely EL2. Which is defined with the same parameter set (focus21 [latitude, longitude], focus22[latitude, longitude], radius2). The reference ellipsoid is WGS84. What is the method to find if EL1 contains EL2 ? So far the best approach I come up with is to convert the EL1 and EL2 to polygons. But it is quite costly. Is there any idea? Any help is appriciated. Waiting for your replys ...


Reading Points as PolyLines

I'm currently developing an Android app where I need to show some information on a map, which is stored as Shapefiles (.shp + .shx + .dbf). I noticed that when the shapefile is opened in QGis, a nice map appears (apparently all points get merged into lines. This is the expected behavour), but after further inspection of the .shp and .shx header, I discovered that the ShapeType defined is actually Point (0x1) and not PolyLine(0x3). I have no other information about the topology, so I don't know how can I achieve the same result. Can somebody explain to me how does QGis infers the edges that connect the points? I'll glady provide further information, if needed. Edit: ![SS from QGis](http://i.stack.imgur.com/ZAb1Y.png) Edit2: Header from the SHP File: https://gist.github.com/elms1990/586fd308be468cc1df0e


Laying tracking data onto map

I apologies if this is a very open ended question but I need a point in the right direction. I've got GeoServer downloaded and running nicely on my Tomcat server. I want to do something quite simple and make a track from a number of GPS coordinates to show the path of satellite across the earth and save this into a png/gif/jpeg so it can be used by another application. I thought you'd sort of upload your background image and then plot the tracks onto it. But I don't think thats correct. I need to have the map as a shape file to right? So what I'd like to do is:               - Have a map of the world     - Read GPS data from mySQL database.     - Plot cooridnates onto map to generate track     - Output track and map as image file.       Is there away to schedule this so it updates the track every few hours? Thanks in advance for any guidance you can offer.


What does this error mean in Mobile Project Center: 'The underlying connection was closed'

While publishing the ArcGIS Mobile project in Mobile project Center, I am getting this error:! ![enter image description here](http://i.stack.imgur.com/8ZTbA.png) This started coming few days ago. And I can not connect with the published project. Please advise which connection this is referring to and how can I correct it.


Is it possible to compare two shape file if they have overlapped and gap between two data or not in QGIS

I am QGIS 2.0 user. I would like to compare two different shape file in the same region. This two shape file's outside border must be the same position. It means the overlap and gap must be zero. In my understanding, Topological checker can not set this kind of rule. Hopefully, I would like to check on the QGIS itself, not using GRASS. Let me know, if somebody knows how to check. Thank you in advance.


How can I use a custom font in QGIS?

I would like to use a custom font to mark certain features on a map. It needs to be defined as an attribute as the features in question are part of a unified text layer that includes various categories. How can I do this?


Openlayers / QuantumGIS, export or sync features Offline

writing from Spain. I am developing a web application based on OpenLayers , GeoServer and PostgisSQL , packaged under OpengeoSuite . The issue has to do with vector layers and their features. Until now, when you add a new item to the map , it is stored in the database by a php call to the database , all Online. What I want now is to work offline in tablets and then add the database items created offline on the tablet , either passing a file to your pc or creating a script to read the new items . I'm evaluating using Quantum Gis for Android since I cannot mount offline application, but do not know how to pass the pc features once created, and the dbf files associated to every shapefile that creates Quantum do not contain the Geom column necessary to know the coordinates of the point.   * QUESTION : What method or program recommend me to create points with attributes in my maps, offline , in a tablet, and how to recover those points in a pc for saving to a Postgis database ? Either by connecting the tablet to PC, or using online from the tablet when available. Thanks in advance for your time.


Why do points shift when reprojected from WGS84 to a state plane coordinate?

We have a Trimble unit we use to collect equipment. We collect the data in WGS84 and differentially correct it and then reproject it to a state plan coordinate system. The points always shift the same distance in the same direction after the reprojection. What is the reason for this? Thanks.


Uploading raster format(*.tif) to postgis through raster2pgsql

Hi I'm new to postgis and raster file format. i currently working on netcdf files and i want to import it into postgis through curl command line tool but its failed for me.so i jump on to GDAL v1.9.2 first i convert it into gtiff format and then upload it into postgis 2.0. I followed the steps in the below link: http://www.postgis.org/documentation/manualsvn/using_raster.xml.html#RT_Raster_Loader to upload it. But i got some issue in that. I used the below command to convert from single netcdf subdataset to single tif file.               gdal_translate -of Gtiff -a_srs EPSG:4326 NETCDF:"3z.nc":aabbc archv_4d.tif      It successfully created. now the problem starts, I used the below command to upload the .tiff file to postgis :                   \9.2\bin>raster2pgsql -I -C -F -s 4326          -d -t 256x256 project/index_u.tif -b 1 public.td_test >test.sql              \9.2\bin>psql -U postgres -d template_postgis_20 -f test.sql -h localhost -p 5432          the response was:          psql:test.sql:2: NOTICE:  CREATE TABLE will create implicit sequence "td_te     st_rid_seq" for serial column "td_test.rid"     psql:test.sql:2: NOTICE:  CREATE TABLE / PRIMARY KEY will create implicit index     "td_test_pkey" for table "td_test"     CREATE TABLE     INSERT 0 1     INSERT 0 1     INSERT 0 1     INSERT 0 1     INSERT 0 1     INSERT 0 1     CREATE INDEX     ANALYZE     psql:test.sql:11: NOTICE:  Adding SRID constraint     CONTEXT:  PL/pgSQL function addrasterconstraints(name,name,name,boolean,boolean,     boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean)      line 53 at RETURN     psql:test.sql:11: NOTICE:  Adding scale-X constraint     CONTEXT:  PL/pgSQL function addrasterconstraints(name,name,name,boolean,boolean,     boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean)      line 53 at RETURN     psql:test.sql:11: NOTICE:  Adding scale-Y constraint     CONTEXT:  PL/pgSQL function addrasterconstraints(name,name,name,boolean,boolean,     boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean)      line 53 at RETURN     psql:test.sql:11: NOTICE:  Adding blocksize-X constraint     CONTEXT:  PL/pgSQL function addrasterconstraints(name,name,name,boolean,boolean,     boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean)      line 53 at RETURN     psql:test.sql:11: NOTICE:  Adding blocksize-Y constraint     CONTEXT:  PL/pgSQL function addrasterconstraints(name,name,name,boolean,boolean,     boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean)      line 53 at RETURN     psql:test.sql:11: NOTICE:  Adding alignment constraint     CONTEXT:  PL/pgSQL function addrasterconstraints(name,name,name,boolean,boolean,     boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean)      line 53 at RETURN     psql:test.sql:11: NOTICE:  Adding number of bands constraint     CONTEXT:  PL/pgSQL function addrasterconstraints(name,name,name,boolean,boolean,     boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean)      line 53 at RETURN     psql:test.sql:11: NOTICE:  Adding pixel type constraint     CONTEXT:  PL/pgSQL function addrasterconstraints(name,name,name,boolean,boolean,     boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean)      line 53 at RETURN     psql:test.sql:11: NOTICE:  Adding nodata value constraint     CONTEXT:  PL/pgSQL function addrasterconstraints(name,name,name,boolean,boolean,     boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean)      line 53 at RETURN     psql:test.sql:11: NOTICE:  Unable to add constraint: enforce_nodata_values_rast     CONTEXT:  PL/pgSQL function _add_raster_constraint_nodata_values(name,name,name)      line 48 at RETURN     PL/pgSQL function addrasterconstraints(name,name,name,text[]) line 94 at assignm     ent     PL/pgSQL function addrasterconstraints(name,name,name,boolean,boolean,boolean,bo     olean,boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean) line 53 a     t RETURN     psql:test.sql:11: NOTICE:  SQL used for failed constraint: ALTER TABLE public.td_test ADD CONSTRAINT enforce_nodata_values_rast CHECK (_raster_constraint_n     odata_values(rast)::numeric(16,10)[] = '{1.26765060022823e+030}'::numeric(16,10)     [])     CONTEXT:  PL/pgSQL function _add_raster_constraint_nodata_values(name,name,name)      line 48 at RETURN     PL/pgSQL function addrasterconstraints(name,name,name,text[]) line 94 at assignm     ent     PL/pgSQL function addrasterconstraints(name,name,name,boolean,boolean,boolean,bo     olean,boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean) line 53 a     t RETURN     psql:test.sql:11: NOTICE:  Returned error message: numeric field overflow     CONTEXT:  PL/pgSQL function _add_raster_constraint_nodata_values(name,name,name)      line 48 at RETURN     PL/pgSQL function addrasterconstraints(name,name,name,text[]) line 94 at assignm     ent     PL/pgSQL function addrasterconstraints(name,name,name,boolean,boolean,boolean,bo     olean,boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean) line 53 a     t RETURN     psql:test.sql:11: WARNING:  Unable to add constraint: 'nodata_values'.  Skipping          CONTEXT:  PL/pgSQL function addrasterconstraints(name,name,name,boolean,boolean,     boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean)      line 53 at RETURN     psql:test.sql:11: NOTICE:  Adding out-of-database constraint     CONTEXT:  PL/pgSQL function addrasterconstraints(name,name,name,boolean,boolean,     boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean)      line 53 at RETURN     psql:test.sql:11: NOTICE:  Adding maximum extent constraint     CONTEXT:  PL/pgSQL function addrasterconstraints(name,name,name,boolean,boolean,     boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean,boolean)      line 53 at RETURN      addrasterconstraints     ----------------------      t     (1 row)...      the postgis table look with only serial and rid [pk] values **NO** raster data. please help me out to find out the solution. thanks in advance.


How to associate polygons with the nearest point on a point layer? (Or ... how do I assign a postal area to a salesman?)

I have a polygon layer of postal code areas and I want to use them to create sales territories based on where our sales people are based. So I assume that I need to work out for each of the polygons which is the nearest point on the sales person point layer. I suppose that I could do this on the nearest part of the polygon or the centroid ... I don't think it would make much difference. Can someone point me in the direction of what this type of analysis is even called? Thanks


How to fix buffer selection issue in ESRI ArcGIS Javascript API and LGIM applications

I had strange LGIM application behavior where the map spatial search did not highlighted the parcels within the Geometry Buffer. What could be the possible reasons? ![Geometry Buffer Selection Issue](http://i.stack.imgur.com/rWLYj.png) LGIM stands for the Local Government Information Model. It is a geodatabase schema ESRI uses to support the ArcGIS JavaScript applications they have.


How to delete isolated pixels in ArcGIS 10.1?

I want to know that it is possible or not to edit the raster pixel in ArcGIS. I want to delete isolated pixel and I want to delete some pixel by manually. please tell me how can i do it manually or application of tools.


Rendering GML in QGIS 1.9.0

I have a set of GML geometries (see here), which I try to load in QGIS. The set contains both points and polygons, but when loaded into QGIS, only points are shown. I tried to follow this thread, but the piece of code suggested for overriding is not present in .gfs file. Any chance to make it work? I tried using OGR2OGR for converting to .shp, but only the points are converted.


How to install PostGIS on Windows?

I want to create an application platform based on PostGIS data. I read the documentation where I saw the prerequisite to install GEOS, GDAL and Proj4. The main problem is I am using windows platform. Can anybody refer me where I can get a full documentation on it to install PostGIS, GDAL, GEOS, Proj4 on windows with all their library? or is it possible to install all the library on windows?


Permission rights with google maps

I created an image in QGIS using the Google Maps layer and I would now like to publish this as part of an academic journal article. However, the journal editor has asked if I have permission to use the Google Maps base layer. Does anyone know if I need to obtain permission from Google to use this image?


Arcgis: Programmatically export to geotiff in different scales

_I am totally new to ArcGis so please bear with me using the wrong terminology._ I would like to programatically export geotiff files in different scales (e.g. 1:25000, 1:50000 etc) using arcgis10. As input I have some shape files and a DEM file as well as a corresponding style. I tried to understand the different possibilities by reading the online- documentation but I am unsure what Api and Calls to use. C#, C++ and python are all ok but I would like to create a fully automated solution. I think a console application might be the best fit but executing a python script from within arcgis might also be ok. I have the following questions:   1. Should I use the .Net API or is a python script sufficient in this scenario?   2. What are the commands / methods I have to call? A code example would be great but pointers to the right documentation is also nice.


Sort feature class in Modelbuilder and populate new ids

I would like to sort my feature class in a model based on coordinates stored in its attribute table. The sorting must be done on `Y_coord` descending, then `X_coord` ascending. The answer to this question looks like what I need, except I only have an ArcEditor Licence, and free versions of Xtools Pro and ET Geowizards. My model performs several other processes as well, so I need the sorting to be included in it. How can I perform the sort with the above limitations? **edit** Looking at the options in the `Sort` tool, I've realised that what I actually want to do is sort my features spatially, starting from the upper left, but to sort based on the `Shape` field requires an ArcInfo licence. The sort tool from ET Geowizards does not allow sorting on this field either. Are there other tools which could do this?


How to classify agriculture employment data?

I have employment data which are huge - about 500 data for each category, like agriculture field, fishing field etc. Now I want to categorise this data in to 7 ranges. Suppose <50=1, >50-100=2, etc. for producing an agriculture map of a country by showing these categories as a legend. I don't want to change my main GIS data. I only want to show this data into small categories in Map. How can I do it?


how do I get php to connect automatically to GeoServer?? so I do not need to go into GeoServer manually by entering the username and password. And how to queries on PostGIS can directly publish and request layer from GeoServer? -thanks for answering-


Google Drive and Backing-up/Syncing GIS Data

Seeing that many of us have three computers these days (e.g.. work, home and laptop) is anyone using Google Drive as their main file repository for GIS data? I like the idea of backing everything up in the cloud and syncing it to all of my computers. I would like to hear anyone's experiences with backing up and syncing their GIS data in Google Drive or one of the other cloud sync/backup services. Are there any serious limitations in using/storing GIS data in the cloud? Is anyone integrating Google Drive with their work server to access data at home? Thanks for your thoughts. **Update:** I spent the $5/month to upgrade my Google Drive storage capacity to 100GB. My GIS files are now synced to all of my computers and backed up in the cloud. I have been working with Arc for several days off my Google Drive with no noticeable issues--so far so good.


georeferencing aerial photos with only known centroid

I have heaps of aerial photos that need to be georeferenced. I have access to ArcMap and ERDAS. All I am given is an excel sheet with the centroid coordinates of each image. Each photo covers an 18m x 24m region. The photos are in a very remote region with no roads or other structures to georeference them to. Is it possible to do this? If so, please help. Cheers, Kate


Setting geometry types importing to PostGIS with OGR

I'm using `ogr2ogr` to import data into a PostGIS database. Each dataset I need to import has (theoretically) the same exact data structure, but it appears that they are _not_ the same geometry types. Example, I've already imported the first layer, and there's a layer called "streets". Running this command, I import the first feature class from the first database (I have GDAL configured to read `.mdb`):               ogr2ogr -f PostgreSQL PG:"host=localhost user=postgres dbname=database port=5432" city1.mdb streets -nln streets -nlt MULTILINESTRING -lco GEOMETRY_NAME=geometry      When I go to append (using the `-append` flag) the same streets layer from `city2.mdb`, I get this error (presumably due to conflicting geometry types for the two "streets" layers:               Warning 1: Geometry to be inserted is of type 3D Line String, whereas the layer geometry type is 3D Multi Line String.     Insertion is likely to fail     ERROR 1: INSERT command for new feature failed.     ERROR:  Geometry type (LineString) does not match column type (MultiLineString)      Is there a straightforward way with `ogr2ogr` to cast geometries to the same type on import? Or a way to handle this by predefining the schema and geometries for each table beforehand?


Creating a sidewalk map for a walkability survey

If I wanted to do my own walkability survey of my city, what tools could I use? I have map data from OpenStreetMap, but I have to consider that streets have two sides (which aren't always accessible from the other side of the street). I guess I need to make a sidewalk map. I'd like to use path-finding algorithms against the data to provide walking directions that avoid sections of streets with poor walkability. Some blocks have sidewalks, some don't, and some streets have a sidewalk on one side but are too busy to cross. The directions would cross the highway at places where it's easier--very easy at the crosswalk, hard at several other streets, impossible at a few streets. I know that weighted path-finding algorithms exist. I'm just wondering how to create this data. I've toyed with taking a map image, putting it in Inkscape, turning on grid lines, and overlaying paths along the sidewalks. With the grid lines, the points of paths connect for path-finding. Path color indicates difficulty of walking. Aside from the path-finding, it'd be a good map to show walkability in town. Anyone have any ideas?


Deriving separate polylines from coded points using ArcObjects (Figure inside)

I have a single point featureclass that represents the left and right sides of a landform (the points are coded as either 'right' or 'left'). Visually I can pick out and digitize two lines along these points, however, I need to do this automatically using ArcObjects or calls to geoprocessing functions in Arc10. The figure below shows the two separate lines that I want to automatically select (in red). ![Polyline selection from points](http://i.stack.imgur.com/9oQFD.png) It is not as easy as selecting by attribute (e.g. 'right') and using Point to Polyline since there is no guarantee that these points are collected in order - the output polylines often jump all over the place. Using a bounding polygon (see figure again) I have attempted to: 1) Split the polygon into lines at its vertices. 2) Select by attribute on the points (e.g. 'right). 3) Select by location on the polylines where they intersect the selected points. The issue with this is I end up with extra lines selected as the points at the top of each desired line intersect with the "top" and "bottom" lines of my polygon. Am I overthinking this? I just need to select the two longest 'sides' of the irregular polygon. Eventually I need to batch process this task on many point featureclasses, so they may not be as clear-cut/spaced as the above example. EDIT: My desired output is a polyline featureclass containing the two separate polylines coded as right and left, respectively.


What's the single point listed in a GPS track's attribute table?

When I import GPS waypoints to QGIS, it includes coordinates that show up in the attribute table. That's great. What I don't understand is when I import tracks. I still get coordinates, as though it is a single point. What point is it? Is the numerical data available for the whole track? Thanks TVZ


How to resize layer using .prj (WKT) projection? (Geographic Coordinate System)

Briefly: How to parametrize .prj WKT file so that I can perform 7 point tranformation (wiki). I know how `false_easting and false_northing` params work, but how can I adjust scale? I do not mean scale_factor' That's the problem description: I have transportation network (vector layer) saved in non-GIS environment (transport modeling software). Network consists of nodes (points) and polylines (road links). It's done mostly from random backgrounds, regardless any projection, coordinates, etc. I need to set appropriate projection for the network. I have accesss to .prj files (if I'm in an say WGS84 projection I can switch to any other projection) So that's what I'm trying:   I try 7 point Helmert Transformation (http://proj.maptools.org/gen_parms.html). I use towgs84 transformation as a WKT param in .prj file, where I assume that rotation matrix is zero (can I do so?) and I calculate only delta_x, and delta_y, and scale param. However it will not work. This is my .prj , params in TOWGS84 do not affect transformation:               PROJCS["UTM 17 (WGS84) in northern hemisphere.",     GEOGCS["WGS 84",         DATUM["WGS_1984",             SPHEROID["WGS 84",6378137,298.257223563],             TOWGS84[0,0,0,0,0,0,100000000000000000000000]],         PRIMEM["Greenwich",0],         UNIT["DMSH",0.0174532925199433],         AXIS["Lat",NORTH],         AXIS["Long",EAST],          PROJECTION["Transverse_Mercator"],     PARAMETER["latitude_of_origin",0],     PARAMETER["central_meridian",0],     PARAMETER["scale_factor",1],     PARAMETER["false_easting",0],     PARAMETER["false_northing",0]]      So I tried to use `false_norting` and `false_easting` params, and those work good, and transform my network proprely, BUT: It will not chcange scale of my network, only position. So how can I rescale my network using .prj file? Thanks for any hints


How to avoid field map change when altering file locations in a model (e.g. using feature class to feature class tool)?

I have a backed up file GDB which I have to use because the current one had some broken models in it. I got this backup delivered on my local D:\ drive and have to copy models from this GDB to a current one, which sits on a network location. Now the problem is that all the links to data in the models have been reset to this local drive, whereas they should point the network location. When I now edit the model and set back the paths to the correct locations (e.g. in the 'feature class to feature class' tool), all my fieldmaps that were adjusted in the backed up models will change (originally deleted fields appear again, adjusted field names get lost). This is a problem because not only do I have a couple of models where this has to be changed, also the feature classes within them have a lot of fields. This means a lot of tedious work that I want to avoid by all means, especially as the correct field mapping is available in my backed up models. I’m using ArcGIS 9.3, the ‘Store relative path names’ option is checked in all the models. May that have caused the problem? It does however not make a difference if I change to ‘absolute path names’ now, as the field map gets lost as soon as I change the location of my data in the model. This question may be a basic one, but I could not find a solution so far, neither in the ArcGIS help nor on the web. Maybe someone can point me in the right direction?


How to calculate the size a bounding box?

What formula could I use to calculate the size in Kilometres of a bounding box based on a given Southwest latitude/longitude and a Northeast latitude/longitude points? The bounding box format is defined as:               bounds = sw_latitude,sw_longitude,ne_latitude,ne_longitude


Integrate QGIS into Hadoop cloud is it possible?

It might be very basic question. Can I integrate QGIS +Hadoop cloud? Currently I am going through similar question and few research articles & researching on my side. In upcoming ArcGIS 10.2 ESRI announced that there will be cloud extension in ArcGIS. So similarly can I integrate open source desktop gis (like QGIS) with hadoop cloud? **Purpose :** Analyze the social networking trend. Any help/suggestion will be great !! thanks in advance : )


Workflow for determining stream gradient?

As far as data goes, I'm working with NHD .shp files, 10m DEMs, and some LIDAR data. My goal is to determine gradient for 100m segments of a network of streams. I'm already able to do this, but I expect that my workflow is nonideal, especially in that I can't deal with branched networks at all. If you all were going about this, what sort of steps would you use? In addition, I posted about the problem here, where I think I did a much better job describing what my goals are.


Is there a software package that allows users to create map tiles?

We would like to create map tiles to be used via mobile device in the field in an offline environment. The imagery will come from our personal collection of aerial images taken via UAV.


Can I define a transformation between 2D and 3D in Geotools?

I'm using Geotools, and I need to convert coordinates between 2D and 3D coordinate reference systems. I recognize that there's no universally-correct way to do this: 2D -> 3D means "making up" a Z coordinate, and 3D -> 2D means losing information. But, is there some way to teach Geotools that, in my particular context, I want to use a specific transformation. E.g.:               +-----+                 +-----+     |     | --- Z = 15 ---> |     |     | 2D  |                 |  3D |     |     | <-- drop Z ---- |     |     +-----+                 +-----+      I know I can transform any particular coordinate manually this way, but I'd like to integrate it so this Just Works. e.g. I can do `CRS.findMathTranform(a, b)` and if `a` is convertible to my 2D CRS, and my 3D CRS is convertible to `b`, it works.


Access to python addin objects/variables in python window

Following up on this question: Can a Python script be executed from a .NET AddIn and also write its standard output and geoprocessor messages to the ArcMap Python Window? Is is possible to access python addin variables/functions in the python window? Direct access doesn't seem to work (e.g. after running my addin, variables from my addin are not available in the python window). But maybe I can get access a different way?


Analyse slope gradient from LiDAR .asc data

I'm new to using QGIS. I have purchased some LiDAR data which downloads as .asc format, is there a way to analyse the terrain slope angle from this data or does it need to be converted into DEM? Ideally I would like to display the steeper slope angles (>35 degrees) as a different colour to the shallower (less steep) slope angles. I'm using QGIS 2.0 Dufour on windows vista platform.


On install "Couldn't load SIP module." - should I reroute Python config or download sip?

(Please note I'm still relatively new to Linux). When I run QGIS after install I get a message 'Could't load SIP module'. I've looked for help and am conflicting advice on what my next step should be. I have previously installed Anaconda (details below). Some advice appears to suggest to change PYTHONPATH but am unclear why or how. There's no way to change the system settings within the graphical interface and I've seen references to "$ export PYTHONPATH=/Applications/QGIS.app/Contents/Resources/python/" but that would appear that I would no longer be able to use my Anaconda install. Or should install another version of Python that includes SIP (PyGIS) or should I download SIP separately? Any help would be greatly appreciated! Couldn't load SIP module. Python support will be disabled. Traceback (most recent call last): File "", line 1, in ImportError: No module named sip Python version: 2.7.5+ (default, Sep 19 2013, 13:53:30) [GCC 4.8.1] QGIS version: 1.7.5-Wroclaw 'Wroclaw', exported Python path: ['/usr/share/qgis/python', '/home/mita/.qgis/python', '/home/mita/.qgis/python/plugins', '/usr/share/qgis/python/plugins', '/home/mita/anaconda/lib/python2.7', '/home/mita/anaconda/lib/python2.7/plat-i386-linux-gnu', '/home/mita/anaconda/lib/python2.7/lib-tk', '/home/mita/anaconda/lib/python2.7/lib-old', '/home/mita/anaconda/lib/python2.7/lib-dynload', '/home/mita/anaconda/lib/python2.7/site-packages', '/home/mita/anaconda/lib/python2.7/site-packages/PIL', '/home/mita/anaconda/lib/python2.7/site-packages/setuptools-0.6c11-py2.7.egg- info']


How can I add a PostGIS query result to OpenLayers?

I'm doing a query on PostGIS that returns a lot of geoemtries as results (more than 1000 which is why I can't use Geoserver CQL) - the query WHERE CLAUSE is being changed each time by the user (its some kind of filter with UI form) What is the best way to add this dynamic query result to map in OpenLayers?


Robust importing of shapefile into spatialite database?

I am trying to import shapefiles into a spatialite database using spatialite- gui. I am having problems with the types of the variables since spatialite-gui has no obvious way to handle shapes which use `","` as decimal separator, and mangles my date fields which are in `YYYY/MM/DD` format. QGIS 1.7.4 imports the same shapefiles without any issues, handling all the conversions transparently, but will not let me save to a Spatialite database. Is there a more robust way of importing a shapefile into Spatialite?


Legal Issues Concerning Intellectual Property

I'd like your opinion: When the GIS Analyst leaves his job (not that I'm in this position, thank god), is he allowed legally to take with him odds and bits he collected or made during his work (small stuff like e.g. small VBA/Python scripts, examples of complex labels), or Like all his bigger Projects it is copyrighted and belongs to the hiring Company? I mean, it's clear that he can't take programs, layers and big stuff (or secret data for that matter), but what about the small non-specific odds and bits?


WMS Service is not displayed

I try to display a wms-server in qgis. It serves several coordinate systems, but only few are working. Here is the link:   http://www.webatlasde.de/arcgis/services/Maps4BW/MapServer/WMSServer   For example EPSG:3857 is displayed, 31467 (which I need) is not   In ArcGis, the WMS is displayed correctly, so the server should be ok.   Can anyone confirm this problem or help me?


How do I create a leader line for a label in my bar graph in ArcMap?

I am displaying a bar graph of various tornado watches and warnings along with the actual tornado touchdowns themselves. Since the tornado touchdowns are such a small number compared to the large quantity of watches/warnings I would like to label the number of tornado touchdowns with their specific number. The problem is that when I display the label (called a mark within the properties window) it covers my points on the graph. **I would like to create some type of "leader line" to move the labels away from my points.** I have enclosed a picture of something kind of like how i want it - I drew the sloppy little black lines is if they were leader lines and somewhere near the end is where I would like to move the labels so they don't cover my points on the graph. **If there are any suggestions I would appreciate it.** Thanks a bunch. ![enter image description here](http://i.stack.imgur.com/EPXPc.png)


how to open layers with multiple geometry types in qgis?

I recently had some difficulties opening some MapInfo TAB vector layers in QGIS because each MapInfo layer contained two different geometry types: point and line. The MapInfo layers had originally been converted from AutoCAD DWGs. When opening in QGIS, only the point features appeared and not the linear features. Is there a way for QGIS to open both geometries of a mixed-geometry layer? Or to convert or save them into two separate-geometry layers? Failing either of those, is there a way to specify which geometry type should be opened (in my case, I was interested in the lines and not the points)? Thanks, George


How to handle a geotiff map together with points with openlayers (or other libraries)?

I create a lot of printed maps using various applications (ArcGIS, QGIS, Microstation, Illustrator etc). Sometimes I'd like to turn this map somewhat interactive on the web. What I'm looking for is a simple way to turn a huge geotiff into tiles and present them in openlayers (or similar) together with some points or polygons on top. If it wasn't for the points and polygons I guess I could have sticked to OpenLayers using zoomify tiles. I want to stay in the local projection and because of that I guess it has to be OpenLayers since Modest Maps etc only handles the spherical Mercator projection? Is there a tool for slicing up my tiff into tiles and then presented so that OpenLayers can handle it together with other point features in the same projection? Do I need a TMS/WMS? A tileserver? Thanks a lot, Oskar.


Is there any method to create monthly LST images?

Is there a method to create monthly LST images from combining MOD11A2(8days) product? The MOD11C3 product which is delivered monthly has larger spatial resolution than MOD11A2. Im trying to create a script in python to do that but it seems to be above my coding skills


I need help with saving Graphics to the a map

I have everything working just like I need it to, but when I add a new graphic to the map it does not save it. So I will add it, close the application, and when I reopen it the graphic is gone. I have AutoSave set to true on both of my feature maps. If anyone knows what my problem is or possible solutions I would really appreciate the help. I am using WPF and C#. I also have my map published to a map server. I don't quite understand the arcGIS for Server though. I have only been using arcGIS for a couple of weeks so bare with me on this one haha. Below is some of my code. If you need me to clear up anything just leave me a comment. XAML:               <Grid x:Name="AddressLayout" Background="White">             <Grid.Resources>                      <esri:Editor x:Key="MyEditor" LayerIDs="WaterMeterLayer" GeometryServiceUrl="http://tasks.arcgisonline.com/ArcGIS/rest/services/Geometry/GeometryServer"/>                 <esri:SimpleRenderer x:Key="WaterMeterRenderer">                                      <esri:SimpleRenderer.Symbol>                     <esri:PictureMarkerSymbol                                            OffsetX="0" OffsetY="31" Source="C:\Program Files\ESRI\Data Interoperability\icons\V2.png"/>                                      </esri:SimpleRenderer.Symbol>                 </esri:SimpleRenderer>                      <esri:SimpleRenderer x:Key="LotRenderer">                     <esri:SimpleRenderer.Symbol>                         <esri:PictureMarkerSymbol                                             OffsetX="0" OffsetY="31" Source="C:\Program Files\ESRI\Data Interoperability\icons\J3.png" />                     </esri:SimpleRenderer.Symbol>                 </esri:SimpleRenderer>                           <esri:PictureMarkerSymbol x:Key="DefaultMarkerSymbol"                                           OffsetX="0" OffsetY="31" Source="C:\Program Files\ESRI\Data Interoperability\icons\K2.png"/>                 <LinearGradientBrush x:Key="PanelGradient" EndPoint="0.5, 1" StartPoint="0.5, 0">                     <LinearGradientBrush.RelativeTransform>                         <TransformGroup>                             <ScaleTransform CenterX="0.5" CenterY="0.5"/>                             <SkewTransform  CenterX="0.5" CenterY="0.5"/>                             <RotateTransform Angle="175" CenterX="0.5" CenterY="0.5"/>                             <TranslateTransform/>                         </TransformGroup>                     </LinearGradientBrush.RelativeTransform>                          <GradientStop Color="#FF145787" Offset="0.15"/>                     <GradientStop Color="#FF3D7FAC" Offset="0.5"/>                     <GradientStop Color="#FF88C5EF" Offset="1.0"/>                 </LinearGradientBrush>                  </Grid.Resources>                            <esri:Map x:Name="MyMap" Background="#FFE3E3E3" WrapAround="True" MouseClick="MyMap_MouseClick" MouseMove="MyMap_MouseMove" MouseRightButtonDown="MyMap_MouseRightButtonDown">                 <esri:Map.Layers>                 <esri:ArcGISTiledMapServiceLayer ID="StreetMapLayer"                            Url="http://services.arcgisonline.com/ArcGIS/rest/services/World_Street_Map/MapServer" />                 <esri:GraphicsLayer ID="CandidateGraphicsLayer" >                     <esri:GraphicsLayer.MapTip>                         <Grid>                             <Rectangle Stroke="Gray"  RadiusX="10" RadiusY="10" Fill="#77FF0000" Margin="0,0,0,5" >                                 <Rectangle.Effect>                                     <DropShadowEffect/>                                 </Rectangle.Effect>                             </Rectangle>                             <Rectangle Fill="#DDFFFFFF" Stroke="DarkGray" RadiusX="5" RadiusY="5" Margin="10,10,10,15" />                             <StackPanel Orientation="Vertical" HorizontalAlignment="Center" Margin="30,20,30,30">                                 <TextBlock Text="{Binding [Address]}" HorizontalAlignment="Left" Foreground="Black" />                                 <TextBlock Text="{Binding [LatLon]}" HorizontalAlignment="Left" Foreground="Black" />                                  </StackPanel>                         </Grid>                     </esri:GraphicsLayer.MapTip>                 </esri:GraphicsLayer>                               <esri:FeatureLayer ID="WaterMeterLayer"                                    Url="http://localhost:6080/arcgis/rest/services///World_Street_Map/MapServer"                                    AutoSave="True"                                    Mode="OnDemand"                                    Renderer="{StaticResource WaterMeterRenderer}"                                    ValidateEdits="True"                                    DisableClientCaching="True">                              <esri:FeatureLayer.MapTip>                             <Grid>                                 <Rectangle Stroke="Gray"  RadiusX="10" RadiusY="10" Fill="#77FF0000" Margin="0,0,0,5" >                                     <Rectangle.Effect>                                         <DropShadowEffect/>                                     </Rectangle.Effect>                                 </Rectangle>                                 <Rectangle Fill="#DDFFFFFF" Stroke="DarkGray" RadiusX="5" RadiusY="5" Margin="10,10,10,15" />                                 <StackPanel Orientation="Vertical" HorizontalAlignment="Center" Margin="30,20,30,30">                                     <TextBlock Text="{Binding [MeterNum]}" HorizontalAlignment="Left" Foreground="Black" />                                     <TextBlock Text="{Binding [Comment]}" HorizontalAlignment="Left" Foreground="Black" />                                     <TextBlock Text="{Binding [Address1]}" HorizontalAlignment="Left" Foreground="Black" />                                     <TextBlock Text="{Binding [Address2]}" HorizontalAlignment="Left" Foreground="Black" />                                           </StackPanel>                             </Grid>                         </esri:FeatureLayer.MapTip>                     </esri:FeatureLayer>                          <esri:FeatureLayer ID="LotLayer"                                    Url="http://localhost:6080/arcgis/rest/services///World_Street_Map/MapServer"                                    AutoSave="True"                                    Mode="OnDemand"                                    Renderer="{StaticResource LotRenderer}"                                    ValidateEdits="True"                                    DisableClientCaching="True" MouseLeftButtonDown="FeatureLayer_MouseLeftButtonDown">                         <esri:FeatureLayer.MapTip>                             <Grid>                                 <Rectangle Stroke="Gray"  RadiusX="10" RadiusY="10" Fill="#77FF0000" Margin="0,0,0,5" >                                     <Rectangle.Effect>                                         <DropShadowEffect/>                                     </Rectangle.Effect>                                 </Rectangle>                                 <Rectangle Fill="#DDFFFFFF" Stroke="DarkGray" RadiusX="5" RadiusY="5" Margin="10,10,10,15" />                                 <StackPanel Orientation="Vertical" HorizontalAlignment="Center" Margin="30,20,30,30">                                     <TextBlock Text="{Binding [MeterNum]}" HorizontalAlignment="Left" Foreground="Black" />                                     <TextBlock Text="{Binding [Comment]}" HorizontalAlignment="Left" Foreground="Black" />                                     <TextBlock Text="{Binding [Address1]}" HorizontalAlignment="Left" Foreground="Black" />                                     <TextBlock Text="{Binding [Address2]}" HorizontalAlignment="Left" Foreground="Black" />                                      </StackPanel>                             </Grid>                         </esri:FeatureLayer.MapTip>                     </esri:FeatureLayer>                      </esri:Map.Layers>             </esri:Map>      C# Code Example of when I add something:               private void MyMap_MouseClick(object sender, Map.MouseEventArgs e)             {                 if (blnAddInformationClicked)                 {                       locatorTask.LocationToAddressCompleted += LocatorTask_LocationToAddressCompleted;                 locatorTask.Failed += LocatorTask_Failed;                           double tolerance = 30;                 locatorTask.LocationToAddressAsync(e.MapPoint, tolerance, e.MapPoint);                           InstructionGrid.Visibility = Visibility.Collapsed;                 PortalGrid.IsEnabled = true;                      }             }                  private void LocatorTask_LocationToAddressCompleted(object sender, AddressEventArgs args)             {                 locatorTask.LocationToAddressCompleted -= LocatorTask_LocationToAddressCompleted;                 if (Information.ContainsKey(txtNumber.Text))                 {                     MessageBox.Show("This Object has already been added");                     PortalGrid.Visibility = Visibility.Visible;                     InformationGrid.Visibility = Visibility.Visible;                     return;                 }                      Address address = args.Address;                 Dictionary<string, object> attributes = address.Attributes;                           if (radWaterMeter.IsChecked == true)                 {                  Graphic graphic = new Graphic()                 {                     Symbol = AddressLayout.Resources["WaterMeterRenderer"] as ESRI.ArcGIS.Client.Symbols.Symbol,                     Geometry = args.UserState as MapPoint                 };                      string address1 = attributes["Street"].ToString();                 string address2 = String.Format("{0}, {1} {2}", attributes["City"], attributes["State"], attributes["ZIP"]);                 string addressCombined =  String.Format("{0},{1},{2},{3}", attributes["Street"].ToString(), attributes["City"], attributes["State"], attributes["ZIP"]); ;                 //Information.Add(txtObjectNumber.Text, addressCombined);                      ObjectProperties record = new ObjectProperties                 {                     ObjectNumber = txtObjectNumber.Text,                     ObjectComments = txtComments.Text,                     ObjectAddress = addressCombined,                 };                      Properties.Add(record);                                //var recordsAsString = string.Join(" ", Properties.Select(p => p.ToString()));                      //MessageBox.Show(recordsAsString);                      if (blnAddInformationClicked == true)                 {                     graphic.Attributes.Add("MeterNum", record.ObjectNumber);                     graphic.Attributes.Add("Comment", record.ObjectComments);                     graphic.Attributes.Add("Address1", address1);                     graphic.Attributes.Add("Address2", address2);                          FeatureLayer graphicsLayer = MyMap.Layers["WaterMeterLayer"] as FeatureLayer;                     graphicsLayer.Graphics.Add(graphic);                     blnAddInformationClicked = false;                 }                 }                      if (radLot.IsChecked == true)                 {                     Graphic graphic = new Graphic()                     {                         Symbol = AddressLayout.Resources["LotRenderer"] as ESRI.ArcGIS.Client.Symbols.Symbol,                         Geometry = args.UserState as MapPoint                     };                          string address1 = attributes["Street"].ToString();                     string address2 = String.Format("{0}, {1} {2}", attributes["City"], attributes["State"], attributes["ZIP"]);                     string addressCombined = String.Format("{0},{1},{2},{3}", attributes["Street"].ToString(), attributes["City"], attributes["State"], attributes["ZIP"]);                          ObjectProperties record = new ObjectProperties                     {                         ObjectNumber = txtObjectNumber.Text,                         ObjectComments = txtComments.Text,                         ObjectAddress = addressCombined,                     };                          Properties.Add(record);                                   graphic.Attributes.Add("MeterNum", record.ObjectNumber);                         graphic.Attributes.Add("Comment", record.ObjectComments);                         graphic.Attributes.Add("Address1", address1);                         graphic.Attributes.Add("Address2", address2);                              FeatureLayer graphicsLayer = MyMap.Layers["LotLayer"] as FeatureLayer;                         graphicsLayer.Graphics.Add(graphic);                         blnAddInformationClicked = false;                           }                  }             private void LocatorTask_Failed(object sender, TaskFailedEventArgs e)             {                 MessageBox.Show("Unable to determine an address. Try selecting a location closer to a street.");             }      I have two different feature layers. There is a radio button there that determines which feature layer and graphic to use. The graphics/renderers are just .png files that are on my C: drive. I add the attributes to the graphics using a list and the information giving by the user. These attributes are then added to the maptip. The graphics are being added, but it does not save it to the map, so when I reopen the application the graphics that I added are not there anymore. I have looked through all of the resources for arcGIS that I could find and there is nothing there that is helping me much. I feel as though the resources are lacking when it comes to beginners... Thank You for your help! :)


Its possible to catch ServiceException in Openlayers?

Im trying to catch exception in OpenLayers using vector layers in conjunction with http protocol, its possible to intercept an ServiceException from geoserver ... it will be great to display some kinda of error message when it happens. here my snippet:               var layer = new OpenLayers.Layer.Vector("Foo", {       projection: new OpenLayers.Projection("EPSG:4326"),       strategies: [new OpenLayers.Strategy.Fixed()],       reportError: true,       protocol: new OpenLayers.Protocol.HTTP({         url: "/geoserver/ows?service=WFS&version=1.0.0&request=GetFeature&typeName=ns:foo&outputFormat=json&viewparams=foo:var;",         format: new OpenLayers.Format.GeoJSON()       })     });      ideas ?? thks !


Best practices for planning hardware requirements for ArcGIS Server deployment scenario

In this question I would like to ask for sharing of some experience in planning ArcGIS Server hardware part in Web mapping projects (Flex, Silverlight). Let suppose that we are developing RIA application and deploying ArcGIS Server with services. I have noticed that very often on the stage of early planning of GIS projects customers would like to know approximate hardware requirements for ArcGIS Server part. Usually it sounds in these ways:   * We need to serve N editor user and M viewers of spatial data. What server do we need (CPU, RAM)? How many ArcGIS servers?   * We already have dedicated server which we are planning to use for ArcGIS Server. How many users (editors/viewers) will be it capable to serve? It is common that on this stage you don't have access to spatial data so it is not possible to perform any tests. But we can estimate rough amount of Map/Feature and Dynamic/Cached services. I have studied Esri System Design Strategies document but I have not found any numbers or formulas to make such rough estimation. As I understand the best precise way is to obtain at least part of the data and perform real tests. So, what is the best practices in such situation? Is it possible to perform such very rough calculations to answer these questions without tests with real data? May be you already faced such problems and have some experience?


Google earth: adding multiple images to a placemark

I have a couple of images that I want to add to a placemark. The problem is when I add the multiple images, google earth displays them one below the other. Is it possible that the image are displayed side by side whenever I click the placemark icon. And If there a dozen of images for the same placemark can I create a "slideshow" of the imagery that displays only one image at a time and swaps when user clicks the placemark button again? Any help is appreciated. Thank you.


models from QGIS 2.0.2 to Python scripts?

I want to automate some tools for QGIS 2.0.2 using Python. Previously I have used ArcMap for automating tools and now i want to use QGIS. In ArcMap I used the ModelBuilder, and more specifically the command Export to Python Script to generate Python code from my model. Now I want to find something similar that will do the same job for my QGIS model i.e. I want to export my model in QGIS to a Python script. Any ideas?


How set correct options for displaying a grid in QGIS composer?

I am using q. 2.5 4291337 in Windows 7 64 bits. I am working with raster and vector layers in EPSG 26714 (UTM NAD27 14N). In this version of QGIS, you are able to put a grid in a different CRS of the layers in Print composer. I tried to play with this option but I have som odd results. When I set that I want a lat-lon Nad27 grid (EPSG 4267) and set the coordinate format to decimal, everything is OK. Longitudes are shown as negative numbers and latitude are positive. ![enter image description here](http://i.stack.imgur.com/oNvYG.png) When I set to show the labels to deg min sec with suffix the label shows wrong the hemisphere indicator. The sufix for my longitudes are labeled as "S" and for my latitudes are labeled "E". ![enter image description here](http://i.stack.imgur.com/0ovoz.png) I tried to replicate the problem with a project with layers in EPSG 4326 (Lat lon WGS84), and I do not have such a problem. ![enter image description here](http://i.stack.imgur.com/0vCBx.png) My question, then, is how to set correctly the options to use a grid in a different CRS than my layers in print composer? In the project with UTM layers, I chose the CRS for my grid to EPSG 4267. I set the spacing as .01666666 for X and Y (which is equivalent to 1 minute). Also I asked a friend of mine that uses Ubuntu and he does not have such a problem.


Openlayers WFS and Mapfile: Possible to specify multiple online resources and generic URL/variables in OpenLayers file?

I have a working mapfile with an associated working OpenLayers file. The mapfile has a `wfs_onlineresource` like so: `"http://one_of_my_servers:8190/cgi-bin/mapserv.exe?"` Then in the OpenLayers js/html file you have the wfs layer like so:               var wfs = new OpenLayers.Layer.Vector("Airports", {                 styleMap: new OpenLayers.StyleMap({                         externalGraphic: "red_pin.png",                         graphicWidth: 18,                         graphicHeight: 18                 }),                 wrapDateLine: true,                 strategies: [new OpenLayers.Strategy.Fixed()],                 protocol: new OpenLayers.Protocol.WFS({                     srsName: "EPSG:3857",                     version: "1.0.0",                     url: "http://one_of_my_servers:8190/cgi-bin/mapserv.exe?map=/tmp/navhome/navcore/airports_viewer_svn/airports_viewer.map&SERVICE=WFS&VERSION=1.1.0&srsName=EPSG:3857",                     featureType: "Airports"                 })         });      Like I said this works fine. Now the issue: we have quite a few servers here. Say you want to take these two above files and go to server 177. Is there a workaround so you don't have to go into the mapfile and specify/change the `wfs_onlineresource` and go into the OpenLayers file and change the `url:` in the wfs layer? **Basically I'm trying to make a 'universal' edition so matter what the server, it will work.** Not really sure if this is possible after reading the WFS Mapfile Documentation: _See section 12.3.3 of the WFS 1.0.0 specification for the whole story about the online resource URL. Basically, what you need is a complete HTTP URL including the http:// prefix, hostname, script name, potentially a “map=” parameter, and and terminated by ”?” or “ &”._ This may be a longshot, but any thoughts on the matter are appreciated.


GeoServer, OpenLayers and updatable views: unable to edit data of WFS-T feature layer

I am developing an web application, where the user is able to view and update some feature attribute of a WFS layer. My problem: the edited data is not saved at all. Here's the situation:   * GeoServer v2.5   * PostGIS / PostgreSQL database v9.1   * pgAdmin III v1.18.1   * Firefox 31.0 with FireBug extension The web application uses OpenLayers (v2.13.1) to display OpenStreetMap (OSM) data. A WFS layer is added to the OSM base layer:               fwSaveStrategy = new OpenLayers.Strategy.Save();             fwSaveStrategy.events.register("success", '', showSuccessMsg);             fwSaveStrategy.events.register("failure", '', showFailureMsg);                  var layerFootway = new OpenLayers.Layer.Vector("WFS OSM Footway", {                 strategies: [new OpenLayers.Strategy.BBOX(), fwSaveStrategy],                 protocol: new OpenLayers.Protocol.WFS({                     version: '1.1.0',                     url: 'http://<the-server-address>/wfs',                     featureType: 'view_footway',                     srsName: 'EPSG:3857'                 }),                      extractAttributes: true,                 visibility: true,                      styleMap: new OpenLayers.StyleMap({                     "default": new OpenLayers.Style(OpenLayers.Util.applyDefaults({                         fillColor: "#ff00ff",                         fillOpacity: 0.5,                         strokeColor: "#ff00ff",                         strokeOpacity: 0.75,                         strokeWidth: 25                     }, OpenLayers.Feature.Vector.style["default"])),                     "select": new OpenLayers.Style(OpenLayers.Util.applyDefaults({                         fillColor: "#00ffff",                         fillOpacity: 0.75,                         strokeColor: "#00ffff",                         strokeOpacity: 0.75,                         strokeWidth: 25                         }, OpenLayers.Feature.Vector.style["select"]))                 })             });      This layer is a specific view in PostgreSQL: the view combines footways/pedestrian sidewalk vector data from OSM with own data of footway attributes (for example surface, type, inclination etc.). Both data sets are joined via the unique OSM ID that is assigned to every object (street, way etc) inside OSM. ![CREATE view script](http://i.stack.imgur.com/YBzJk.png) For interaction purposes, a user can click/touch on a footway and gets confronted with a popup stating some feature attributes and way attributes. First, the code to achieve this:               var selectControlTouch = new OpenLayers.Control.SelectFeature(layerFootway, {                 hover: false,                 autoActivate: true,                 onSelect: function(feature) {addPopup(feature);},                 onUnselect: function(feature) {saveAttributes(feature); deletePopup(feature);}             });      So, by click/touch on a footway, a popup is initiated and displayed on the screen (see picture below). There the user is able to alter a value of the footway attributes (just a few numbers right now to get it working). If the user clicks elsewhere on the map, the changed attribute value will (hopefully) be saved and the popup will be closed. ![Screenshot 01 - Popup in OSM web application](http://i.stack.imgur.com/P6jjI.jpg) To be able to store data to the view, I added specific rules to it in PostgreSQL: ![Rules for the view](http://i.stack.imgur.com/LUg6k.png) The save function of my application looks like this:               function saveAttributes(feature) {     feature.state = OpenLayers.State.INSERT;     feature.attributes.wn_test = document.getElementById('iWn_test').value;     fwSaveStrategy.save(feature);}      So far, so good. Just a roundup:   * I've created a view in PostgreSQL.   * This view is published via GeoServer.   * In OpenLayers it is initiated as a WFS layer and displays correctly inside my browsers.   * If the user clicks on a footway, a popup comes up with some informations about that way.   * Inside the popup lies an `<input>` field. Goal: write that (edited) data back to the database.   * Once the user clicks elsewhere on the map, the data from the `<input>` field would be saved and the popup closes. BUT:   * If I call `fwSaveStrategy.save(feature);`, the transfer seems to execute correctly, but no data will be transfered though giving me an alert message "Everything went fine".   * If I call `fwSaveStrategy.save([feature]);`, the response from the server is: `view_footway is read-only`. From my web search, the last message seems common among WFS-T users. So I clicked through several pages, tutorials and message boards (including this one) to find a solution. It seemed, that I needed a Primary Key Metadata table to get updatable views. So I added one with the name `gt_pk_metadata_table`to my database: ![Values of metadata table](http://i.stack.imgur.com/PV6CO.png) Then I referenced that table in my GeoServer data store (under "Primary key metadata table") and saved. Unfortunately, that procedure has not solved my problems, it's gotten worse. Before, I could ask for the `osm_id` via `feature.attributes.osm_id` as this field belongs to the view. Now, with the metadata table, the result for requesting the `osm_id` is `undefined`. Therefore, the saving function won't execute correctly:               <ows:Exception exceptionCode=InvalidParameterValue>     <ows:ExceptionText>Error performing insert: Source was null in trying to create a reprojected feature collection!</ows:ExceptionText>      Does anyone has any idea, how to solve this and get the application to save the attributes edited via popup? I am dealing with this issue for 3 days now and am more and more willing to write a .php file with an SQL script instead of using WFS-T.    # UPDATE 26.08.2014 Just played around a little bit without complicated views and specific rules, but all comes down to one error message:   * created a new table containing OSM data (copied several columns from existing table into new one via `SELECT INTO`)   * added a primary key column to that table to get unique ids because osm_id is not unique (via `ALTER TABLE` and `SERIAL`)   * checked GRANTS in pgAdmin and set `public` and the owner to ALL   * tested UPDATE/SELECT in pgAdmin on that new table   * published the table as new layer in GeoServer   * added the new layer in OpenLayers   * all data shows up fine But as soon as the `fwSaveStrategy.save()` function is called, the XML output reads               <ows:Exception exceptionCode="InvalidParameterValue">     <ows:ExceptionText>Update error: Invalid output dimension (must be 2 or 3)</ows:ExceptionText>      I've searched the web for it but found nothing useful. Any ideas?    # UPDATE 27.08.2014 Here's the WFS exception:               <ows:ExceptionReport xmlns:ows="http://www.opengis.net/ows"      xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance"     version="1.0.0" xsi:schemaLocation="http://www.opengis.net/ows http://<serveraddress>/geoserver/schemas/ows/1.0.0/owsExceptionReport.xsd">     <ows:Exception exceptionCode="InvalidParameterValue">     <ows:ExceptionText>Update error: Invalid output dimension (must be 2 or 3)</ows:ExceptionText>     </ows:Exception>     </ows:ExceptionReport>


Exporting GeoTIffs from Data Driven Pages for use on a Tablet

I'm a beginner at ArcGIS 10.2 and have encountered a small problem that I couldn't solve by googling. I'm trying to export a set of raster and vector data for use on a tablet. The aim is to find LiDAR features in the field using the GPS from the tablet and I want to be able to switch between various visualizations and the polygons and lines of the interpretation. For this purpose I try to batch create GeoTiffs with various visualizations of the data and export them using Data Driven Pages, since it lets you very conveniently tile a map and I wanted to create a hard copy mapbook of the data anyway. I found an ArcPy script that batch exports Data Driven Pages to tiffs and this is working fine, with the major problem that it's not possible to include spatial information when exporting from a Layout. Since I'm generating a very lage amount of tiles (over 600) and would like to do so with various visualizations, manually georeferencing them is not an option. Any suggestions, tips or work arounds on how to solve this would be much appreciated.


Is there a PostGIS function to get a tangent and normal line?

I'm looking for a function that will return a tangent (and normal) for a LINESTRING in PostGIS. I think this is the 3rd time I've come across the need for this function and everytime, I end porting whatever I'm working on to c#/ArcObjects in order to use the QueryTangent/QueryNormal method rather than make the effort to research a PostGIS option. Is there an existing functions in PostGIS (I'm not tied to a specific version)? And if it doesn't exist, has anyone already written it and wouldn't mind posting it as an answer?


How to stop ArcGIS Python script printing warnings and other outputs to the console

I've written a Python script using arcgisscripting to do various processing steps on some data. I'm planning to run it from a stand-alone command-line window (ie. not within ArcGIS at all) and, as such, have a number of `print` commands in my code to show its status. However, when I run it from the console I also get a number of other things printed to the console including warnings (such as a warning that my data has no projection) and the full output of the Nearest Neighbour calculation command. These aren't needed, and will clog up the output (I'll be running this as a batch job with the output piped to a file). Is there any way to get rid of these? I was hoping for a gp.quiet = true command or a gp.silent = true, but can't seem to find any of those. Any ideas?


Can an external application start ArcCatalog.exe and immediately navigate to a specific folder?

I'm looking something similar to explorer's command line switches:   `explorer.exe /select,c:\examples\condensed\need\obtains\dispute\under\billion\floor` or `explorer.exe /e,c:\examples\condensed\need\obtains\dispute\under\billion\floor` Can an external application create a folder `c:\examples\condensed\need\obtains\dispute\under\billion\floor` and then start ArcCatalog.exe in that location? ![enter image description here](http://i.stack.imgur.com/5z7vV.png)


How to connect "floating" model variable to the tool that uses the inline variable?

This question likely has a simple answer or the behaviour is by-design and is not worth fixing. The simple model below uses an Expression to 'Select Layer By Attribute." When this model is run, the expression will remain the same but the ID attribute value can changed based on user input. Therefore I have exposed a string variable as a model parameter, and am using inline variable substitution to add the user's value in the query expression (the datatype of the [LeaseID] field happens to be a string and not integer/numeric). I would like to have the "Lease ID" model parameter connected to the "Select Layer By Attribute" tool. This is purely aesthetic and is not required for the model to run successfully. I believe it would look better, be more intuitive, and be easier for new users of Modelbuilder to understand where the model parameter is being inserted into the tool. I could just expose the Expression parameter of the tool as a model parameter...but then anyone running the tool would modify the whole expression and could cause unnecessary error(s). Any thoughts on how to do this...I would prefer that the variable not be floating in free space. Thanks. ![Modelbuilder image](http://i.stack.imgur.com/oc0mS.png) Full size image


Why this two datasets don't match?

I have a postGIS database in SRS 4326 (WGS84) with several Europe layers working well, but when I load a shapefile projected to the same SRS with QGIS, the layers don't match. To explain it better I've made two screenshots: The pink one are the German postcodes from my postGIS database, that are matching with all the rest of my data ![enter image description here](http://i.stack.imgur.com/HeD38.png) The second one, green polygons, are German postcodes too, but from a shapefile. Both of them are in WGS84 (EPSG:4326, I've reprojected it with QGIS, and shows it in layer propperties). I've marked the coordinates below, to see that the it seems that the second one's coordinates are not correct. ![enter image description here](http://i.stack.imgur.com/dVBxF.png) Some Idea that how could I fix the green dataset to overlap with the rest? (pink) The source of the green dataset is this URL: http://arnulf.us/PLZ Thanks!


How to overlay Google Earth to Map Object?

> **Possible Duplicate:**   >  Is it possible to expose Google Maps as background map through ArcGIS > server? i am developing a gis system using mapobjects , i have to prove that i can overlay the google earth in my gis application, how can i do that.? i appreciate if you give a sample program, im using Visual Basic 6. tnks.. Paul


MapQuery (jQuery + Openlayers 2) unbind function

I have small gis app build on MapQuery, but development of this libary was stopped 2 years ago. MapQuery have bind function but do not have unbind one.               bind: function(types, data, fn) {             var self = this;             // A map of event/handle pairs, wrap each of them             if (arguments.length === 1) {                 var wrapped = {};                 $.each(types, function(type, fn) {                     wrapped[type] = function() {                         return fn.apply(self, arguments);                     };                 });                 this.events.bind.apply(this.events, [wrapped]);             } else {                 var args = [types];                 // Only callback given, but no data (types, fn), hence                 // `data` is the function                 if (arguments.length === 2) {                     fn = data;                 } else {                     if (!$.isFunction(fn)) {                         throw ('bind: you might have a typo in the function name');                     }                     // Callback and data given (types, data, fn), hence include                     // the data in the argument list                     args.push(data);                 }                 args.push(function() {                     return fn.apply(self, arguments);                 });                 this.events.bind.apply(this.events, args);             }             //this.events.bind.call(this.events, types, function() {             //  data.apply(self, arguments);             //});             //this.events.bind.call(this.events, types, function() {             //  data.apply(self, arguments);             //});             //this.events.bind.apply(this.events, arguments);             //this.events.bind.call(this.events, types, $.proxy(data, self));             //this.events.bind.apply(this.events, arguments);//.bind(this);             //this.events.bind.apply(this.events, $.proxy(arguments));//.bind(this);             //this.events.bind.apply(this.events, $.proxy(arguments));//.bind(this);             //this.events.bind(types, data, fn);//.bind(this);             //this.events.bind.call(this.events, types, data, fn);//.bind(this);             return this;         },         /**To subscribe to the triggered events, you need to bind to the mapuuu.           map.bind('myEvent', function(evt) {          console.log('the values are: ' + evt.data[0] + ' and ' + evt.data[1])      });      I can't figure out how to write unbind function.


Why can't I perform any operation over a 700MB MIF file containing 5 multipolygons? It takes forever

I'm trying to work with a MIF file containing 5 really big multipolygons, but every operation I try takes hours, or just crash because of memory issues. Mainly using ogr2ogr, to transform MIF file to other format, like SHP or PostGIS. Unsuccessful, so far. The file weights 700MB. I could load it into PostGIS using ogr2ogr. I just left the command running all night and it finished silently. The size of the table in disk is about 160MB. So, I suspect something went wrong (looks too small) **UPDATE** This is the header of MIF file               Version 300     Charset "Neutral"     Delimiter ","     CoordSys Earth Projection 8, 28, "m", -3, 0, 0.9996, 500000, 0     Columns 2       Float Float       String Char(32)     Data      After that, 5 regions with:   * 124222 polygons   * 462270 polygons   * 798815 polygons   * 1124915 polygons   * 1706436 polygons About the PostGIS table, I'm trying to check it. It contains 5 rows, like the original MIF file, but every single operation I try to execute over this table using the geometry column also takes hours. I tried to add a new integer field, and after 4 hours is still running. With ps, the process says "ALTER TABLE waiting". This is the schema of the table                   Column    |       Type       |                         Modifiers     --------------+------------------+----------------------------------------------      ogc_fid      | integer          | not null default nextval('my_table_ogc_fid_seq'::regclass)      wkb_geometry | geometry         |      float        | double precision |      string       | character(32)    |     Indexes:         "my_table_pk" PRIMARY KEY, btree (ogc_fid)         "my_table_geom_idx" gist (wkb_geometry)     Check constraints:         "enforce_dims_wkb_geometry" CHECK (st_ndims(wkb_geometry) = 2)         "enforce_srid_wkb_geometry" CHECK (st_srid(wkb_geometry) = 23030)      So, how could I handle this big MIF file? I can't do anything with it, neither with the result PostGIS table.


How can I add aerial imagery using REST services as a basemap in CartoDB?

I want to add my county's aerial imagery as a basemap in CartoDB, same as I use in OSM's iD editor but it is not loading. I am adding in the XYZ url as this: http://gis.slocounty.ca.gov/arcgis/rest/services/Aerials/2011_Combined_WGSWMAS/MapServer/tile/{z}/{y}/{x} The aerials are hosted here: http://gis.slocounty.ca.gov/arcgis/rest/services/Aerials On a sidenote, is there a way to remove previously added basemaps in CartoDB?


Extract/Split raster dataset by attribute (OBJECTID)

I have a raster dataset with an attribute table generated from a series of image classification (unsupervised) and processing tools (majority filter + boundary clean). The attribute table shows pixels of the classified image, grouped by classification i.e. the "OBJECTID" field. I want to extract the portion of the image by "OBJECTID" to divide the classified image into separate rasters and ultimately, separate feature classes. e.g. take all pixels with OBJECTID = 2 and copy into new raster or convert to feature class This has proved somewhat impossible for me so hopefully there's something simple that I'm missing?


ArcGIS file naming conventions: is there an advantage to using underscores?

In an introductory GIS course, it was insisted that all spaces within file and folder names be replaced with underscores. Why is this taught? Is there a clear benefit of doing so?


ENVISAT ASAR .N1 file is not being opened by ERDAS IMAGINE 2011

I am not able to open ENVISAT ASAR .N1 file in ERDAS IMAGINE 2011. The software has direct read option for files of this format. Whenever I try to open the file, software closes automatically giving an error message. Can anyone point out the problem and suggest solutions.


Different Styles Within Vector Layer

Is there any way to configure OpenLayers styling method? I have a vector layer and I would like to style the features depending on a value of a flag variables, so I can't use rules. Lets say I have two types of features: houses and restaurants. houses flag is "H" and restaurants flag is "R". How can I give them a different style within a vector layer? **EDIT**               function search(term)     {     var flag;           if (regex_3.test(term))   //term is 3 letters              {          flag = "H";     ..     ..     ..         }              else  if (regex_more.test(term))                     {          flag = "R";     ..     ..         }     }      and when I find the feature, I add it to a vector layer.but the thing is I can't find a way to make two styles depending on the flag value.


How to find the overlap of three polygons?

I have three polygon areas showing different land types and I want to create a new layer based on where the three layers overlap each other. I thought this was a subtraction query but havent been able to achieve this in ArcGIS or QuantumGIS. Any help would be greatly appreciated.


Can OpenGeo Suite be deployed in production environment (using WAR package) for free?

Is there a possibility to use OpenGeo Suite in production environment (using the WAR package) for free? Actually we are planning to deploy on a glassfish (alternatively Tomcat) container openGeo Suite. In a first step we think we can manage to set everything ourselves, thus we don´t need any support of OpenGeo. The pricing for the WAR package with glassfish costs 45000$ and what I understood it is only for the support. Can we get the WAR package for free if we don´t need any Support?


What is an extension error?

I am getting an extension error when I try and intersect some points and polygons. ![enter image description here](http://i.stack.imgur.com/OauKN.png) Here is a photo of what I am trying to intersect: ![enter image description here](http://i.stack.imgur.com/C687z.png) Any help is appreciated.


QGIS: how to set data defined text angle without setting x and y coordinate using new style labelling

I have a problem setting angled text labels with the new-style labelling dialog. The situation is this: I have a shape file which only contains only labelling texts. This is a point layer, and it contains all information needed for the label, like text size, font etc. - while the placement of the label is steered by the cordinates of the geometry of the point, so there are no attributes containing the X and y coordinate. This is legacy data, so I have no choice in the matter - the labels aren't linked to objects in other shapefiles, they're just there by themselves, having originally made up a separate toponym layer. If I use the old labelling mechanism, I can enter the text angle field in the appropriate place (in my data it's called TextAngle), and the text is displayed at the desired angle. But if I use the new-style labelling dialog, the text angle field in the 'position' section of the 'data defined settings' tab remains greyed out unless I set the fields 'x-coordinate' and 'y-coordinate'. Yet I don't have fields containing x and y coordinates in my data - they are implicit, being the very coordinates my data points occupy. They appear when I get info for the elements as 'derived', but they aren't in the attribute table. If I try and trick the dialog by entering something else in the 'x-coordinate' and 'y-coordinate' fields, I can enter my text angle field, but it has no effect. Am I missing something? I'd really rather use the new-style labelling, since the buffering is nicer. Help welcome!


QGIS calculate distance of point along a line

I have a polyline (movement path) and points recorded along the line. I would like to calculate for each point the distance along the line in Quantum GIS. I hope for a result which will add a new column to the point attribute table indicating the distance along which it is on the line. I'd appreciate any thoughts. Mike.


R.Mapcalc giving only blank maps as result

No matter what function I put in the r.mapcalc in the GRASS GUI I get a blank map full of null values. I even tried DEM+0 (I can see my DEM) and this also yielded a blank map.


Configuring a WorldImage data store starting from a folder containing tif and tfw files

i have a folder of tif and tfw files and i'd like to show the raster in a map using geoserver. Even reading the geoserver tutorial for this kind of datastore (http://docs.geoserver.org/stable/en/user/data/worldimage.html) i coldn't find the way to show the raster. In the same folder i have also other files (shape, aux, properties, ..) but i don't know if they are useful for publish the map or not. Can you please tell me how to do? Thank you.


Morris Chart X Axis labels disappear

I have created a chart using cartodb and Morris chart library and I have been having problems with the x axis labeling. Originally I used the chart with bootstrap and only one label showed up. So I then implemented the same chart on a blank page with just the morris css and only two labels showed up on the X axis. So I have recreated the code in jfiddle using static data and now I have every other label show up. I am at a loss for why the same code is producing these labeling issues. jfiddle I tried using xLabelAngle with no results unless I set the value to -1 then all of my labels show up however the alignment is off.


Force all labels visible in GeoServer

How can I force all labels to be visible in GeoServer I have tried with both               <VendorOption name="conflictResolution">false</VendorOption>      and               <VendorOption name="spaceAround">-1</VendorOption>      But it doesn't seem to help much. The pointsymbolizer are always visible though. And sometimes it doesn't even feel like its "crouded" when the labels isn't visible. ![enter image description here](http://i.stack.imgur.com/tlGDf.png) The (I guess) interesting parts in my sld file looks like following:               <PointSymbolizer>         <Graphic>                          <Mark>                            <WellKnownName>circle</WellKnownName>                            <Fill>                              <CssParameter name="fill">#FFFFFF</CssParameter>                            </Fill>                           <sld:Stroke>                             <sld:CssParameter name="stroke">#9F0000</sld:CssParameter>                             <sld:CssParameter name="stroke-width">2</sld:CssParameter>                           </sld:Stroke>                                                     </Mark>                          <Size>9</Size>                        </Graphic>                      </PointSymbolizer>                                         <TextSymbolizer>                     <sld:Geometry>                       <ogc:Function name="centroid">                         <ogc:PropertyName>geom</ogc:PropertyName>                       </ogc:Function>                     </sld:Geometry>                                             <Label>                           <ogc:PropertyName>label</ogc:PropertyName>                         </Label>                         <Font>                           <CssParameter name="font-family">Arial</CssParameter>                           <CssParameter name="font-weight">Bold</CssParameter>                           <CssParameter name="font-size">14</CssParameter>                         </Font>                         <LabelPlacement>                           <PointPlacement>                             <AnchorPoint>                               <AnchorPointX>0.5</AnchorPointX>                               <AnchorPointY>0.5</AnchorPointY>                             </AnchorPoint>                             <Displacement>                               <DisplacementX>0</DisplacementX>                               <DisplacementY>-15</DisplacementY>                             </Displacement>                           </PointPlacement>                         </LabelPlacement>                         <Halo>                           <Radius>                             <ogc:Literal>2</ogc:Literal>                           </Radius>                           <Fill>                             <CssParameter name="fill">#9F0000</CssParameter>                           </Fill>                         </Halo>                         <Fill>                           <CssParameter name="fill">#FFFFFF</CssParameter>                         </Fill>                         <VendorOption name="conflictResolution">false</VendorOption>                          </TextSymbolizer>


In OpenLayers can I replace a kml's URL for OpenLayers.Protocol.HTTP with the actual text of the kml?

I had an OpenLayers map that reads in and displayed a kml like this:               var vector_format = new OpenLayers.Format.KML({extractAttributes: true});              var vector_protocol = new OpenLayers.Protocol.HTTP({             url: 'http://www.geos.ed.ac.uk/~s0825955/thomlines.kml',             format: vector_format         });              var vector_strategies = [new OpenLayers.Strategy.Fixed()];             vector_Layer = new OpenLayers.Layer.Vector('More Advanced Vector Layer',{                 protocol: vector_protocol,                 strategies: vector_strategies             });      I am producing the with a perl file and have the kml as a variable. Can I replace the url (`{ url: 'http://www.geos.ed.ac.uk/~s0825955/thomlines.kml',`) with the actual kml text as stored in my scalar variable `$newkml` ? and if so how?


How to convert MODIS (Vegetation Continuous Fields - MOD44B) to (.img and hdr) in ENVI

Please can anyone help me to convert MODIS (Vegetation Continuous Fields) to (.img and hdr) formats? In ENVI, anytime i try using the MODIS conversion toolkit, a dialog box appears telling me to define BAD range. I did input the following (0-100, 200 and 253) but none of them work. I was successful when i opened it using the Generic format tool(HDF) but can not reproject it unto geographic coordinate system. Please I need an assistance. Eric.


How can I show excel file with addresses inside in Google Earth?

I have excel file with 1 column of address , city and country.   I want to "import" it into Google Earth and get placemarks of the addresses (2800 POI).   How can I do it ?


How to cut multiple neighboring raster images

I'm looking for a way to cut a larger number of neighboring aerial raster images (each sized 1x1 km, resolution 10 cm per pxel) into four smaller raster images (size 500 x 500 m). I haven't done this before, is there a (preferably QGIS-based) way to do this (semi-)automatically? I don't need any "slippy map"-like tiles - this is only necessary because my colleagues use a rather exotic and antique CAD software that can't handle larger image sizes... (and our state geodata department doesn't sell smaller custom sized tiles any more) Thanks a lot!


How do I debug a GP Service in ArcServer?

I have a geoprocessing service that was successfully published to our ArcServer instance. I'm using 10.1 which means we could not publish unless it actually worked run as a tool in ArcMap. But now when I try calling form my client I'm not getting anything. And the messages provided in ESRI logging are not very helpful. I tried attaching to the ArcGIS server service but I never had my breakpoints reachable. So what do I attach to in order to listen for exceptions???


Image map plugin not working with point data QGIS

I have a point shapefile and I am trying to create a html image from it using the html image map plugin; however, I can not get it work. I get the error:               Wrong geometry type      Any ideas on this? Also, can you zoom into the layer and render the layer, i.e. have labels turning on at a certain scale. Seems a useful tool.


How do I successfully install the LecoS plug-in for QGIS 1.8

I have been trying to follow the instructions given here, but repeatedly get an error message pop-up titled "Mount" stating Access is Denied. Does anyone have a solution please?


Why do I get 000210 Error trying to wirte output to in_memory workspace in ModelBuilder

My model gives this error message when it tries to do the Clip function ![Error Message](http://i.stack.imgur.com/P7v3C.jpg) Here is the first half of the model showing the parameters used in the clip function ![Model](http://i.stack.imgur.com/jTbmY.jpg) I have successfully run this model when the clip function output was writing to disk but when I changed it to in_memory it stopped working. Any suggestions? I really don't want to write all 3000 files to my hard drive.


How to remove all selection from all registered layers using QGIS plugin?

I want to remove all selection from registered layer in QGIS using Python. I tried to clear the map canvas using:               self.iface.mapCanvas().clear()       self.iface.mapCanvas().refresh()       but the selection persists.


Convert from Gauss-Kruger to Bing coordinates / projection

I am trying to convert a shapefile with `ogr2ogr` from Gauss-Kruger (zone 4) to Bings coordinate system / projection using the advice given in http://alastaira.wordpress.com/2011/01/23/the-google-maps-bing-maps-spherical- mercator-projection/:               ogr2ogr -overwrite -s_srs EPSG:31468 -t_srs Bing-Projection.txt out.shp in.shp      Where `Bing-Projection.txt` contains               PROJCS["Popular Visualisation CRS / Mercator",      GEOGCS["Popular Visualisation CRS",       DATUM["WGS84",         SPHEROID["WGS84", 6378137.0, 298.257223563, AUTHORITY["EPSG","7059"]],       AUTHORITY["EPSG","6055"]],      PRIMEM["Greenwich", 0, AUTHORITY["EPSG", "8901"]],      UNIT["degree", 0.0174532925199433, AUTHORITY["EPSG", "9102"]],      AXIS["E", EAST], AXIS["N", NORTH], AUTHORITY["EPSG","4055"]],     PROJECTION["Mercator"],     PARAMETER["semi_minor",6378137],     PARAMETER["False_Easting", 0],     PARAMETER["False_Northing", 0],     PARAMETER["Central_Meridian", 0],     PARAMETER["Latitude_of_origin", 0],     UNIT["metre", 1, AUTHORITY["EPSG", "9001"]],     AXIS["East", EAST], AXIS["North", NORTH],     AUTHORITY["EPSG","3785"]]      But I get the following output:               ERROR 6: No translation for Mercator to PROJ.4 format is known.     Failed to create coordinate transformation between the     following coordinate systems.  This may be because they     are not transformable, or because projection services     (PROJ.4 DLL/.so) could not be loaded.      See this gist for the full output. Any advice is highly appreciated.


Finding distance (in time) from nearest point out of 4 along a road network

I have four points representing harbors along the coast. I also have a roads layer that has info on each road segment such as length, speed limit, and travel time. I calculated travel time with [(length in miles) / (speed limit in mph)] * 60 for the time it takes to traverse each road segment in minutes. What I want to do now is: Work out from each harbor point along the road networks to get the total travel time to the nearest harbor for each road segment. All the tools I've found that do something similar to this require both a start and end point. Is there a way to do this without turning each road segment into a point itself? Any help is appreciated, thank you. P.S. To better visualize my end goal, I essentially want a sort of "Euclidean Distance" effect, except it needs to stay along the roads layer and will total in time rather than meters/miles/etc.


why is the CDB_QuantileBins function not found in cartoDB?

What is wrong with the following SQL query on my data? I found pretty clear documentation (1,2) that the function is supposed to work.               SELECT CDB_QuantileBins(array_agg(numeric_column_of_mine), 10) from my_table


My fusion table map with search not working

My maps are not working. It shows the search box but not the map:               // ]]></script>

 Postcode        (with space e.g. HP11 2QW) <input id="search-string_0" type="text"> <input        onclick="changeMap_0()" type="button" value="Search">


Openlayers: OGC Filter problem

I have a problem with OGC filter. My layer in OpenLayers has this filter:               filter: new OpenLayers.Filter.Comparison({                  type: OpenLayers.Filter.Comparison.BETWEEN,                  property: "datetime",                  lowerBoundary: "2013-01-10T12:11:50.501Z",                  upperBoundary: "2013-01-10T12:11:55.501Z"                                                          })       and the post to server is:               <ogc:Filter xmlns:ogc="http://www.opengis.net/ogc">     <ogc:PropertyIsBetween>     <ogc:PropertyName>datetime</ogc:PropertyName>     <ogc:LowerBoundary>     <ogc:Literal>2013-01-10T12:11:50.501Z</ogc:Literal>     </ogc:LowerBoundary>     <ogc:UpperBoundary>     <ogc:Literal>2013-01-10T12:11:55.501Z</ogc:Literal>     </ogc:UpperBoundary>     </ogc:PropertyIsBetween>     </ogc:Filter>      but the reply is:               <wfs:FeatureCollection numberOfFeatures="0" timeStamp="2013-01-10T12:11:15.100Z" xsi:schemaLocation="AidAssistant http://127.0.0.1:8080/geoserver/wfs?service=WFS&version=1.1.0&request=DescribeFeatureType&typeName=AidAssistant%3Av_dati_sensori_geom_BNU1 http://www.opengis.net/wfs http://127.0.0.1:8080/geoserver/schemas/wfs/1.1.0/wfs.xsd" xmlns:ogc="http://www.opengis.net/ogc" xmlns:tiger="http://www.census.gov" xmlns:cite="http://www.opengeospatial.net/cite" xmlns:nurc="http://www.nurc.nato.int" xmlns:sde="http://geoserver.sf.net" xmlns:wfs="http://www.opengis.net/wfs" xmlns:topp="http://www.openplans.org/topp" xmlns:it.geosolutions="http://www.geo-solutions.it" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:AidAssistant="AidAssistant" xmlns:sf="http://www.openplans.org/spearfish" xmlns:ows="http://www.opengis.net/ows" xmlns:gml="http://www.opengis.net/gml" xmlns:xlink="http://www.w3.org/1999/xlink">     <gml:featureMembers></gml:featureMembers>     </wfs:FeatureCollection>      I tried to use this filter:               filter: new OpenLayers.Filter.Logical({              type: OpenLayers.Filter.Logical.AND,              filters: [                      new OpenLayers.Filter.Comparison({                          type: OpenLayers.Filter.Comparison.GREATER_THAN_OR_EQUAL_TO,                          property: "datetime",                          value: 2013-01-10T12:25:40.612Z                      }),                      new OpenLayers.Filter.Comparison({                         type: OpenLayers.Filter.Comparison.LESS_THAN_OR_EQUAL_TO,                         property: "datetime",                        value: 2013-01-10T12:25:50.612Z                      })              ]      })       and the post is:               <ogc:Filter xmlns:ogc="http://www.opengis.net/ogc">     <ogc:And>     <ogc:PropertyIsGreaterThanOrEqualTo>     <ogc:PropertyName>datetime</ogc:PropertyName>     <ogc:Literal>2013-01-10T12:25:40.612Z</ogc:Literal>     </ogc:PropertyIsGreaterThanOrEqualTo>     <ogc:PropertyIsLessThanOrEqualTo>     <ogc:PropertyName>datetime</ogc:PropertyName>     <ogc:Literal>2013-01-10T12:25:50.612Z</ogc:Literal>     </ogc:PropertyIsLessThanOrEqualTo>     </ogc:And>     </ogc:Filter>      and the reply is the same, 0 features. There are features in my layer, indeed if I do this request:               http://127.0.0.1:8080/geoserver/wfs?service=wfs&version=1.1.0&request=GetFeature&typeName=AidAssistant:v_dati_sensori_geom_BNU1&srsName=EPSG:404000&sortBy=datetime+D      I receive a lot of features. But if I try this:               http://127.0.0.1:8080/geoserver/wfs?service=wfs&version=1.1.0&request=GetFeature&typeName=AidAssistant:v_dati_sensori_geom_BNU1&srsName=EPSG:404000&filter=%3Cogc:Filter%20xmlns:ogc=%22http://www.opengis.net/ogc%22%3E%3Cogc:PropertyIsGreaterThanOrEqualTo%3E%3Cogc:PropertyName%3Edatetime%3C/ogc:PropertyName%3E%3Cogc:Literal%3E2013-01-10T12:25:40.612Z%3C/ogc:Literal%3E%3C/ogc:PropertyIsGreaterThanOrEqualTo%3E%3C/ogc:Filter%3E      I receive 0 features. Sorry for the long post :( Can someone help me? Thanks. **Update:** Maybe I have found the error. Geoserver to do this query on the DB:               SELECT "node_id","datetime","geom".STAsBinary() as "geom"       FROM "v_dati_sensori_geom_BNU1" WHERE  "datetime" BETWEEN '2013-01-3 16:45:53.000' AND '2013-01-10 16:45:58.000'      but this does not work because the right query on SQL SERVER is:               SELECT "node_id","datetime","geom".STAsBinary() as "geom"       FROM "v_dati_sensori_geom_BNU1" WHERE  "datetime"        BETWEEN convert(datetime, '2013-01-3 16:45:53.000', 120)  AND convert(datetime, '2013-01-10 16:45:58.000', 120)      Workaround for this problem?