Way to split and create buffer space between polygons in QGIS or ArcMap

After performing a feature extraction of rooftops, some of my buildings have been joined together as one polygon. I want to split these large polygons so that each new division is only covering one building. I would also like to simultaneously create a small buffer zone or gap between the new polygons. I can do it manually as follows but am looking for a tool/plugin that will do this.   1. Split large polygon with line between the two buildings.    2. Split polygon again parallel to previous split.    3. Delete newly created "sliver" polygon covering the gap between the two buildings. Does any tool or method currently exist in QGIS or ArcMap?


From Raster to Masked and Filtered Polygon Vector

I have a set of rasters (geologic risk) from which I want to extract zones with a specific color, clipping them with a polygon (administrative borders), outputting a shapefile.   I'm using QGis (and GDAL from command line).   I know I can use clipper, polygonize and show only features with the value I need, but:   1. it takes ages to polygonize the raster, and I don't want much of the information in it;   2. clipper doesn't work because the polygon border is larger than each raster, but fits inside the sum of all rasters   3. to make things complicated, the zone I would like to extract don't have the same color value in all rasters (but do have the same rgb color) How can I solve it in a fast way?


Help with gap mask filling and landsat 7 in arcgis 10

i have downloaded images from Landsat 7 ETM+ SLC-off with corresponding gap mask. which is the best method for using it with arcgis 10, i tried the mask and statistical method described in arcgishelp but i have some troubles. (after i aply the mask function, the statistical funcion need the input raster that is somewhere in a temporal file that cant access it.) i will apreciate if someone can help me with step by step guide, or some alternative solution. from already thank you very much.


Choropleth Mapping Using Quantile Method in CartoDB

Does anyone know how CartoDB calculates the buckets for a choropleth map using the quantile method? My understanding of choropleth maps that use the quantile method is that they will assign an equal number of items to each bucket (see http://www.ncgia.ucsb.edu/cctp/units/unit47/html/comp_class.html). A recent map I made with CartoDB generated a map that did not do this. Instead, most items were assigned to one bucket. ![Choropleth Map Example](http://i.stack.imgur.com/LeACr.png) Is CartoDB using a different type of quantile method than what I described? Any help understanding is appreciated.


Why does ArcGIS 10 Buffer not provide accurate geodesic buffer?

I created a shapefile following these criteria as recommended for Performing Geodesic Buffering: The input feature class must be a point or multipoint dataset. (shapefile is a point) The input feature class must have a geographic coordinate system. (GCS: Longitude/Latitude (WGS 84)) The buffer distance must be specified as a linear unit such as kilometers, miles, or feet. (Radius field (String) added with inputs 3 kilometers, 4.5 kilometers...) I am not, however, getting accurate geodesic buffers. If I perform a buffer with the same shapefile and input a Linear Unit (3 Kilometers) instead of specifying by the field, I get accurate geodesic buffers. This was the buffer used with a specified distance which gave me accurate geodesic buffers:               arcpy.Buffer_analysis("test","C:/Temp/testBuffer.shp","3 Kilometers","FULL","ROUND","LIST","Site_Name")      This was the buffer used with a specified field and is NOT giving me accurate geodesic buffers:               arcpy.Buffer_analysis("test","C:/Temp/testBuffer2.shp","Radius","FULL","ROUND","LIST","Site_Name")


Merging styles in ArcGIS, changing default Style location?

Just upgraded to 10.2.1 and I have 'lost' all my personal styles. I have a couple questions.   1. Is there a way to merge my three different style files into one?   2. Is there a way to change the default location of the styles file? (I know that I can change the location each time I create a new style but I want the location (of my choosing) to come up by default so I do not have to keep changing it.)


Placing TIF image on the top of WMS layer

How to place TIF file on the top of WMS layer in QGIS? When I import TIF file, it has different coordinate system, and it is placed far from WMS layer, but I want to put this TIF layer on the top of WMS layer. I know how to make them transparent. Thank You!


Error occured during Geodatabase compressing

When I am doing geodatabase administration in ArcSDE on an Oracle Database, the following error pops up when I run the **compress** tool ORA-00604: error occurred at recursive SQL level 1 ORA-01000: maximum open cursors exceeded. With out compressing, our database performance is coming down.What might be the error and what should be the remedy.


Calculate 2D line of site from line layer and polygon layer using ArcGIS 10

I'm trying to run an analysis where I need to clip a radiating lines layer by land features to get the maximum distance from the central point to the nearest obstacle. What I have is a large line layer containing approximately 40,000 points, each with 32 radiating lines extending out at equal intervals (e.g. 11.25 degrees) and a vector polygon file representing coastlines and reefs. The problem I'm having is that the "Clip" and "Erase" functions remove the sections of line on land but when a line encounters a small obstacle the tool only removes the section over land and leaves the portion of line which continues on afterwards (e.g. adjacent bays, other side of peninsulas, etc). This is especially problematic in regions with complex coastal features. What I need to do is remove the hanging lines after the nearest land. I've looked into line of sight analyses but the only one's I can find rely on 3D surfaces. Does anyone know of a 2D version where all line fragments will be removed after reaching the nearest obstacle? This needs to be as automated as possible as I'm dealing with large datasets (24 to be exact, each with 40,000 points). The end goal is to produce a layer to calculate wave fetch.


Creating an ArcGIS Add-In to Export SLD

If I wanted to create an add-in for ArcGIS that exported the styling of a layer to SLD how much work would it be and what languages could I use? More info: I'm an experienced programmer, and have done a lot of coding for GIS systems but am not a user of ESRI products and have never tried writing an extension for any of their products. I would only need to export basic styling to SLD, not write a comprehensive tool that works for everything.


GeoDjango: read geometry column in geography field

I have a Django application that needs to work with an existent PostGIS database, apart from its own database. I've configured multidatabase support, created a database router, created models for the existent tables in the PostGIS database and told syncdb that, for that models, managed=False (so, it won't create the tables in the database, they already exist) My problem comes because I need to create distance queries using one of those PostGIS tables. The table stores points objects with srid 4326, and I need to know how many of these points are within a given distance of a reference point (like ST_DWithin). I'd like to use geography for the distance queries. So, I could define the geom field in the model that maps to that table with geography=True. But the problem is I just have read access to the PostGIS database. So, I can't change the column type to geography. I don't want to write a raw query, transforming the column to geography. Unless there's no other method. Is there a geodjango- friendly way to deal with this?


How to convert ArcMap color ramps to QGIS?

I have both ArcGIS 10.0 and QGIS 2.0 How can I reproduce a specific ArcGIS color ramp in QGIS? If possible, I would like to simply convert the Esri *.style file but perhaps there is another approach that I am not aware of. I would prefer not to have to manually create a new QGIS color ramp from scratch.


Editing and Deleting Query Tables in ArcMap

How do I edit the query for a Query Table after it has been created? How do I permanently delete a Query Table? Even after I do right-click -> remove, it still appears in dropdown lists.


How to know the required scene required to cover an area?

I want to know how to find the nos. of scene of LISS 4 data required to cover an area(say 84000 sq km)? any suggestion/solution/calculation?


Analyzing spatial differences: conceptual and technical questions

My question is both conceptual and technical, but more conceptual than technical. It's a little long and multi-pronged, but I figured the best answers would consider everything, so thanks for bearing with me. I list specific questions later. Right now I'm research the concept of "academic mismatch" in undergraduate education, or when students attend colleges where their academic background is substantially different than the typical student (e.g., a "little fish in a big pond:" a student with low grades/scores compared to his peers, or a "big fish in a little pond:" a student with high grades/scores compared to her peers). While I'm on solid footing with outcomes analysis (do mismatched students do better/worse in school?), I want to explore the mechanisms of attending a school as a mismatch in the first place. There are a lot of reasons for this, some obvious and well explored (preferential admissions policies), others less obvious and less well explored (where I hope to contribute). One question I want to explore, and many others (including my dissertation advisor) have found interesting, pertains to _where students live_ and their likelihood of mismatching. The basic hypothesis is that when a student has fewer colleges to choose from within his "local college market," the more likely he is to attend a college that is _less selective_ than he could because it's convenient. Likewise, for non-selective schools, students sometimes attend less academically appropriate colleges _because they can get in, even if they're outmatched by their typical peers._ There's some evidence out there in local contexts (e.g., Chicago), but very little in a national context, and almost none in rural contexts (which I/we find most fascinating). This is the basic question I want to answer: **Does the concentration of "good-fit" colleges that are geographically near to a student associate with that students likelihood of attending a college where her academic background matches her peers?** So from this, I need to do the following:   1. Create a measure of the concentration of good fit colleges. I have methods for deciding if a college is a good fit between student and college, but assessing concentration is more challenging.   2. Test whether this measure is associated with students "mismatch status"   3. Create maps demonstrating patterns found Here are my questions:   1. My intuition for concentration is to simply count the colleges that are a good fit within a specific distance (or a market defined by a shape file) for each student and do appropriate analysis. What better ideas do you have?   2. Are there geospatial statistical tests that are better than my plan (logistic regression analysis with outcome being match/not match, independent variable being the measure of college concentration?   3. Software wise, should I bite the bullet and learn ArcGIS? I have experience as a web developer so my intuition is to run to web-based software, but I'm not sure where to start. My problem seems like it could be solved without ArcGIS, and if I can learn software I can use in web development, all the better. As background, I have the following data:   1. A student's home zip code and a number of variables on their academic background (high school gpa, SAT, etc.)   2. A data file of all accredited colleges and their academic profiles (which I can assign geo-coordinates) Everyone is in the USA. I have time to learn, but am in proposal-writing stage, so insights in that regard are especially helpful. As background, I have basically zero experience with GIS, but very good statistical analysis, computer programming skills, and data manipulation skills (for a social science researcher) . Put another way, I'm comfortable with research design and Stata, R, Python, MySQL, but have never done anything with GIS (besides some web analytics using geolocation with IPs). Thanks for the help.


Reprojecting WGS 1984 Web Mercator (EPSG:3857) in Python with GDAL

I am reprojecting rasters in python using GDAL. I need to project several tiffs from geographic WGS 84 coordinates to WGS 1984 Web Mercator (Auxiliary Sphere), in order to use them later in Openlayers together with OpenStreetMap and maybe Google maps. I am using Python 2.7.5 and GDAL 1.10.1 from here, and transforming coordinates using advices from here (my code is below). In short, I imported _osgeo.osr_ and used _ImportFromEPSG(code)_ and _CoordinateTransformation(from,to)_. I tried firstly **EPSG(32629)** which is UTM zone 29 and got this projected raster (more or less fine), so the code seems to be correct: ![utm](http://i.stack.imgur.com/wVeiD.png) Then I used **EPSG(3857)** because I've read this and this questions and found that it is the correct recent **valid code**. But the raster is created with no spatial reference at all. It is far away up in WGS 84 data frame (but will be ok if I switch the data frame to Web Mercator). ![3857](http://i.stack.imgur.com/EQ1yM.png) With **EPSG(900913)** the output is georeferenced, but shifted about 3 raster cells to the north: ![900913](http://i.stack.imgur.com/LUM60.png) When I reproject the raster using **ArcGIS** (export in WGS_1984_Web_Mercator_Auxiliary_Sphere) the result is nearly fine: ![arcgis](http://i.stack.imgur.com/zbiCE.png) And when I use the old code **102113 (41001,54004)** the result is perfect: ![54004](http://i.stack.imgur.com/5Dk7N.png) The summary of my tests using all codes:               3857: far away up (missing georeference)     3785: far away up (like 3857)     3587: far away right     900913: slightly jumped up     102100: python error     102113: perfect     41001: perfect     54004: perfect     ArcGIS (web merc. aux.): good      **So my questions are:**   * Why does the correct EPSG code give me wrong results?    * And why the old codes work fine, aren't they deprecated?   * Maybe my GDAL version is not good or I have errors in my python code? The code:                   yres = round(lons[1]-lons[0], 4)  # pixel size, degrees         xres = round(lats[1]-lats[0], 4)         ysize = len(lats)-1  # number of pixels         xsize = len(lons)-1         ulx = round(lons[0], 4)         uly = round(lats[-1], 4)  # last         driver = gdal.GetDriverByName(fileformat)         ds = driver.Create(filename, xsize, ysize, 2, gdal.GDT_Float32)  # 2 bands         #--- Geographic ---         srs = osr.SpatialReference()         srs.ImportFromEPSG(4326)  # Geographic lat/lon WGS 84         ds.SetProjection(srs.ExportToWkt())         gt = [ulx, xres, 0, uly, 0, -yres]  # the affine transformation coeffs (ulx, pixel, angle(skew), uly, angle, -pixel)         ds.SetGeoTransform(gt)  # coords of top left corner of top left pixel (w-file - center of the pixel!)         outband = ds.GetRasterBand(1)         outband.WriteArray(data)         outband2 = ds.GetRasterBand(2)         outband2.WriteArray(data3)         #--- REPROJECTION ---         utm29 = osr.SpatialReference()     #    utm29.ImportFromEPSG(32629)  # utm 29         utm29.ImportFromEPSG(900913)  # web mercator 3857         wgs84 = osr.SpatialReference()         wgs84.ImportFromEPSG(4326)         tx = osr.CoordinateTransformation(wgs84,utm29)         # Get the Geotransform vector         # Work out the boundaries of the new dataset in the target projection         (ulx29, uly29, ulz29) = tx.TransformPoint(ulx, uly)  # corner coords in utm meters         (lrx29, lry29, lrz29) = tx.TransformPoint(ulx + xres*xsize, uly - yres*ysize )         filenameutm = filename[0:-4] + '_web.tif'         dest = driver.Create(filenameutm, xsize, ysize, 2, gdal.GDT_Float32)         xres29 = round((lrx29 - ulx29)/xsize, 2) # pixel size, utm meters         yres29 = abs(round((lry29 - uly29)/ysize, 2))         new_gt = [ulx29, xres29, 0, uly29, 0, -yres29]         dest.SetGeoTransform(new_gt)         dest.SetProjection(utm29.ExportToWkt())         gdal.ReprojectImage(ds, dest, wgs84.ExportToWkt(), utm29.ExportToWkt(), gdal.GRA_Bilinear)         dest.GetRasterBand(1).SetNoDataValue(0.0)         dest.GetRasterBand(2).SetNoDataValue(0.0)         dest = None  # Flush the dataset to the disk         ds = None  # only after the reprojected!         print 'Image Created'      Thank you for suggestions!


Mapinfo incorrect polyline (boundary) map to region conversion

How to solve issue regarding incorrect boundary to region conversion? I have boundary map plotted on cosmetic layer when I convert it to region it results in line with overlap to other region. The two images are : ![enter image description here](http://i.stack.imgur.com/2WXeI.gif) ![enter image description here](http://i.stack.imgur.com/pFQf3.gif) How do I resolve this issue?


Does Bing Maps have a 3D mode?

The following wikipedia page says that 3d mode in bing maps was discontinued http://en.wikipedia.org/wiki/Comparison_of_web_map_services Whats the reason for that? Is there any official announcement or page regarding this from Microsoft?


Is it possible to reference data using mapped network drives instead of UNC path in arcgis server services?

Has anyone managed to create a service using datasets referenced via a mapped network drive letter on windows? ESRI support claims it is not "supported", but I want to make sure it is actually not possible or just "not supported" as in "not recommended". the permissions are properly set up for my SOC account, since the services created using UNC path to reference datasets are working fine. I guess it all comes down to making the ArcGISSOC account understand those letters. I think I tried logging on as the ArcGISSOC and mapping the drives manually, but it doesn't seem to work. Any ideas? Thanks


Postgis installation error

I am using fedora machine. I want to install postgis on my system. So this is what I have done   1. Installed Postgresql and supporting packages using YUM.    2. Created directory /sources to put all the necessary packages for postgis.    3. successfully installed geos,libxml & proj using ./config,make,make install.    4. successful in make & make install of postgis. But when I give "make check" I am getting the following error(I have already started the postgres server):               [root@localhost sources]# cd postgis/          [root@localhost postgis]# ls           postgis-1.5.0SVN  postgis-1.5.0SVN.tar.gz           [root@localhost postgis]# cd postgis-1.5.0SVN           [root@localhost postgis-1.5.0SVN]# make check           make -C liblwgeom           make[1]: Entering directory `/sources/postgis/postgis-1.5.0SVN/liblwgeom'           make[1]: Nothing to be done for `all'.           make[1]: Leaving directory `/sources/postgis/postgis-1.5.0SVN/liblwgeom'           make -C postgis           make[1]: Entering directory `/sources/postgis/postgis-1.5.0SVN/postgis'           Makefile.pgxs:17: warning: overriding recipe for target `install'      /usr/lib/pgsql/pgxs/src/makefiles/pgxs.mk:92: warning: ignoring old recipe for      target `install'           Makefile.pgxs:63: warning: overriding recipe for target `installdirs'      /usr/lib/pgsql/pgxs/src/makefiles/pgxs.mk:140: warning: ignoring old recipe for      target `installdirs'           Makefile.pgxs:82: warning: overriding recipe for target `uninstall'      /usr/lib/pgsql/pgxs/src/makefiles/pgxs.mk:164: warning: ignoring old recipe for      target `uninstall'           make[1]: Nothing to be done for `all'.           make[1]: Leaving directory `/sources/postgis/postgis-1.5.0SVN/postgis'           make -C liblwgeom/cunit check           make[1]: Entering directory `/sources/postgis/postgis-1.5.0SVN/liblwgeom/cunit'           WARNING:           configure was unable to find CUnit which is required for unit testing.      In order to enable unit testing, you must install CUnit and then re-run configure.           make[1]: Leaving directory `/sources/postgis/postgis-1.5.0SVN/liblwgeom/cunit'           make -C regress check           make[1]: Entering directory `/sources/postgis/postgis-1.5.0SVN/regress'           make -C ../postgis REGRESS=1 DESTDIR=/sources/postgis/postgis-1.5.0SVN/regress/00-regress-install install           make[2]: Entering directory `/sources/postgis/postgis-1.5.0SVN/postgis'           Makefile.pgxs:17: warning: overriding recipe for target `install'      /usr/lib/pgsql/pgxs/src/makefiles/pgxs.mk:92: warning: ignoring old recipe for target `install'           Makefile.pgxs:63: warning: overriding recipe for target `installdirs'      /usr/lib/pgsql/pgxs/src/makefiles/pgxs.mk:140: warning: ignoring old recipe for target `installdirs'           Makefile.pgxs:82: warning: overriding recipe for target `uninstall'      /usr/lib/pgsql/pgxs/src/makefiles/pgxs.mk:164: warning: ignoring old recipe for target `uninstall'      /bin/sh /usr/lib/pgsql/pgxs/src/makefiles/../../config/mkinstalldirs '/sources/postgis/postgis-1.5.0SVN/regress/00-regress-install/lib'      /bin/sh /usr/lib/pgsql/pgxs/src/makefiles/../../config/install-sh -c -m 755  postgis-1.5.so '/sources/postgis/postgis-1.5.0SVN/regress/00-regress-install/lib/postgis-1.5.so'      /bin/sh /usr/lib/pgsql/pgxs/src/makefiles/../../config/install-sh -c -m 644 ./../spatial_ref_sys.sql '/sources/postgis/postgis-1.5.0SVN/regress/00-regress-install/share/contrib/postgis-1.5'      /bin/sh /usr/lib/pgsql/pgxs/src/makefiles/../../config/install-sh -c -m 644 postgis.sql '/sources/postgis/postgis-1.5.0SVN/regress/00-regress-install/share/contrib/postgis-1.5'      /bin/sh /usr/lib/pgsql/pgxs/src/makefiles/../../config/install-sh -c -m 644 uninstall_postgis.sql '/sources/postgis/postgis-1.5.0SVN/regress/00-regress-install/share/contrib/postgis-1.5'      /bin/sh /usr/lib/pgsql/pgxs/src/makefiles/../../config/install-sh -c -m 644 postgis_upgrade_15_minor.sql '/sources/postgis/postgis-1.5.0SVN/regress/00-regress-install/share/contrib/postgis-1.5'      /bin/sh /usr/lib/pgsql/pgxs/src/makefiles/../../config/install-sh -c -m 644 postgis_upgrade_14_to_15.sql '/sources/postgis/postgis-1.5.0SVN/regress/00-regress-install/share/contrib/postgis-1.5'      /bin/sh /usr/lib/pgsql/pgxs/src/makefiles/../../config/install-sh -c -m 644 postgis_upgrade_13_to_15.sql '/sources/postgis/postgis-1.5.0SVN/regress/00-regress-install/share/contrib/postgis-1.5'      /bin/sh /usr/lib/pgsql/pgxs/src/makefiles/../../config/install-sh -c -m 755  postgis-1.5.so '/sources/postgis/postgis-1.5.0SVN/regress/00-regress-install/lib/postgis-1.5.so'           make[2]: Leaving directory `/sources/postgis/postgis-1.5.0SVN/postgis'           make -C ../loader REGRESS=1 DESTDIR=/sources/postgis/postgis-1.5.0SVN/regress/00-regress-install install           make[2]: Entering directory `/sources/postgis/postgis-1.5.0SVN/loader'           Makefile:108: warning: overriding recipe for target `clean'      /usr/lib/pgsql/pgxs/src/makefiles/pgxs.mk:211: warning: ignoring old recipe for target `clean'      /bin/sh /usr/lib/pgsql/pgxs/src/makefiles/../../config/install-sh -c pgsql2shp /sources/postgis/postgis-1.5.0SVN/regress/00-regress-install/bin      /bin/sh /usr/lib/pgsql/pgxs/src/makefiles/../../config/install-sh -c shp2pgsql /sources/postgis/postgis-1.5.0SVN/regress/00-regress-install/bin           make[2]: Leaving directory `/sources/postgis/postgis-1.5.0SVN/loader'      sed 's,$libdir,/sources/postgis/postgis-1.5.0SVN/regress/00-regress-     install/lib,g' ../postgis/postgis.sql > postgis.sql           psql: could not connect to server: No such file or directory              Is the server running locally and accepting              connections on Unix domain socket "/tmp/.s.PGSQL.5432"?           Creating spatial db postgis_reg           createdb: could not connect to database postgres: could not connect to server: No such file or directory              Is the server running locally and accepting              connections on Unix domain socket "/tmp/.s.PGSQL.5432"?           createlang: could not connect to database postgis_reg: could not connect to server: No such file or directory Is the server running locally and accepting              connections on Unix domain socket "/tmp/.s.PGSQL.5432"?           psql: could not connect to server: No such file or directory              Is the server running locally and accepting              connections on Unix domain socket "/tmp/.s.PGSQL.5432"?           dropdb: could not connect to database postgres: could not connect to server: No such file or directory              Is the server running locally and accepting              connections on Unix domain socket "/tmp/.s.PGSQL.5432"?           Something went wrong (no postgis installed in postgis_reg).       For details, check /tmp/pgis_reg_5204/regress_log           make[1]: *** [check] Error 1           make[1]: Leaving directory `/sources/postgis/postgis-1.5.0SVN/regress'           make: *** [check] Error 2           6.Even with the error, i continued to create template_postgis db from postgresql shell:       On giving the command: psql template_postgis -U postgres postgis/postgis.sql I am getting the folowing error:           psql:postgis/postgis.sql:6859: NOTICE:  type "geography" is not yet defined      DETAIL:  Creating a shell type definition.      CREATE FUNCTION      psql:postgis/postgis.sql:6865: NOTICE:  argument type geography is only a shell      CREATE FUNCTION      CREATE FUNCTION      CREATE TYPE      psql:postgis/postgis.sql:6893: NOTICE:  type "gidx" is not yet defined      DETAIL:  Creating a shell type definition.      CREATE FUNCTION      psql:postgis/postgis.sql:6899: NOTICE:  argument type gidx is only a shell      CREATE FUNCTION      psql:postgis/postgis.sql:7744: NOTICE:  aggregate geomunion(geometry) does not exist, skipping      DROP AGGREGATE      psql:postgis/postgis.sql:7745: NOTICE:  aggregate st_geomunion(geometry) does not exist, skipping      DROP AGGREGATE      psql:postgis/postgis.sql:7746: NOTICE:  aggregate accum_old(geometry) does not exist, skipping      DROP AGGREGATE      psql:postgis/postgis.sql:7747: NOTICE:  aggregate st_accum_old(geometry) does not exist, skipping      DROP AGGREGATE      psql:postgis/postgis.sql:7750: NOTICE:  function box2d_overleft(box2d,box2d) does not exist, skipping      DROP FUNCTION      psql:postgis/postgis.sql:7751: NOTICE:  function box2d_overright(box2d,box2d) does not exist, skipping      DROP FUNCTION      psql:postgis/postgis.sql:7752: NOTICE:  function box2d_left(box2d,box2d) does not exist, skipping      DROP FUNCTION      psql:postgis/postgis.sql:7753: NOTICE:  function box2d_right(box2d,box2d) does not exist, skipping      DROP FUNCTION      psql:postgis/postgis.sql:7754: NOTICE:  function box2d_contain(box2d,box2d) does not exist, skipping      DROP FUNCTION      psql:postgis/postgis.sql:7755: NOTICE:  function box2d_contained(box2d,box2d) does not exist, skipping      DROP FUNCTION      psql:postgis/postgis.sql:7756: NOTICE:  function box2d_overlap(box2d,box2d) does not exist, skipping      DROP FUNCTION      psql:postgis/postgis.sql:7757: NOTICE:  function box2d_same(box2d,box2d) does not exist, skipping      DROP FUNCTION      psql:postgis/postgis.sql:7758: NOTICE:  function box2d_intersects(box2d,box2d) does not exist, skipping      DROP FUNCTION      psql:postgis/postgis.sql:7759: NOTICE:  function st_box2d_overleft(box2d,box2d) does not exist, skipping      DROP FUNCTION      psql:postgis/postgis.sql:7760: NOTICE:  function st_box2d_overright(box2d,box2d) does not exist, skipping      DROP FUNCTION      psql:postgis/postgis.sql:7761: NOTICE:  function st_box2d_left(box2d,box2d) does not exist, skipping      DROP FUNCTION      psql:postgis/postgis.sql:7762: NOTICE:  function st_box2d_right(box2d,box2d) does not exist, skipping      DROP FUNCTION      psql:postgis/postgis.sql:7763: NOTICE:  function st_box2d_contain(box2d,box2d) does not exist, skipping      DROP FUNCTION      psql:postgis/postgis.sql:7764: NOTICE:  function st_box2d_contained(box2d,box2d) does not exist, skipping      DROP FUNCTION      psql:postgis/postgis.sql:7765: NOTICE:  function st_box2d_overlap(box2d,box2d) does not exist, skipping      DROP FUNCTION      psql:postgis/postgis.sql:7766: NOTICE:  function st_box2d_same(box2d,box2d) does not exist, skipping      DROP FUNCTION      psql:postgis/postgis.sql:7767: NOTICE:  function st_box2d_intersects(box2d,box2d) does not exist, skipping      DROP FUNCTION      psql:postgis/postgis.sql:7768: NOTICE:  function st_addbbox(geometry) does not exist, skipping      DROP FUNCTION      psql:postgis/postgis.sql:7769: NOTICE:  function st_dropbbox(geometry) does not exist, skipping      DROP FUNCTION      psql:postgis/postgis.sql:7770: NOTICE:  function st_hasbbox(geometry) does not exist, skipping      DROP FUNCTION      psql:postgis/postgis.sql:7771: NOTICE:  function cache_bbox() does not exist, skipping      DROP FUNCTION      psql:postgis/postgis.sql:7772: NOTICE:  function st_cache_bbox() does not exist, skipping      DROP FUNCTION      psql:postgis/postgis.sql:7773: NOTICE:  function transform_geometry(geometry,text,text,pg_catalog.int4) does not exist, skipping      DROP FUNCTION      psql:postgis/postgis.sql:7774: NOTICE:  function collector(geometry,geometry) does not exist, skipping      DROP FUNCTION      psql:postgis/postgis.sql:7775: NOTICE:  function st_collector(geometry,geometry) does not exist, skipping      DROP FUNCTION      psql:postgis/postgis.sql:7776: NOTICE:  function geom_accum(geometry[],geometry) does not exist, skipping      DROP FUNCTION      psql:postgis/postgis.sql:7777: NOTICE:  function st_geom_accum(geometry[],geometry) does not exist, skipping      DROP FUNCTION      psql:postgis/postgis.sql:7778: NOTICE:  function collect_garray(geometry[]) does not exist, skipping      DROP FUNCTION      psql:postgis/postgis.sql:7779: NOTICE:  function st_collect_garray(geometry[]) does not exist, skipping      DROP FUNCTION      psql:postgis/postgis.sql:7780: NOTICE:  function geosnoop(geometry) does not exist, skipping      DROP FUNCTION      psql:postgis/postgis.sql:7781: NOTICE:  function jtsnoop(geometry) does not exist, skipping      DROP FUNCTION      psql:postgis/postgis.sql:7782: NOTICE:  function st_noop(geometry) does not exist, skipping      DROP FUNCTION      psql:postgis/postgis.sql:7783: NOTICE:  function st_max_distance(geometry,geometry) does not exist, skipping      DROP FUNCTION      -bash-4.1$       Any solution for this is of great help. thanx in advance


ArcGIS automatically convert string variables to numeric?

If arcgis is automatically coverting string values to numeric variables in a text document generated by stata, how do you change those variables back into string? I tried opening up the fields under properties in arccatalog and it says I cannot change them because the tables are not empty


How do I set the column type for a CSV file in ArcView?

I have a CSV file that I'm using to create a layer in ArcView 10. The first column contains a mix of numbers in the first few records, with strings further down. On import, ArcView casts the entire column as "long" which makes any record containing a string value display "null" instead. Is it possible to tell ArcView the column type for import? Can I force the column to import as a string rather than a long?


How to automatically email interim results from model builder?

I want to get model builder to email me the results of processes at certain stages so that if there are issues I can remotely login and deal with them as I do with python scripts. I have a python based tool to email the results but don't see a way to get results from individual processes. You can right-click on a process in model builder and view results but how can I get this info to be passed into the email script? Worst case scenario - I guess I can export the model to python and add the emailer script but is there an easier way? ![enter image description here](http://i.stack.imgur.com/P6vAm.jpg) Model Code -this is a simple model to test. The main model to use this on will run over a weekend or more as the raster-polygon take about 6 hours for some tiles.               # ---------------------------------------------------------------------------     # model_emailstatus.py     # Created on: 2011-10-11 16:28:25.00000     #   (generated by ArcGIS/ModelBuilder)     # Description:      # ---------------------------------------------------------------------------          # Import arcpy module     import arcpy          # Load required toolboxes     arcpy.ImportToolbox("D:/SCRATCH/Projects/106/Toolbox.tbx")               # Local variables:     Key_Resource_Area_shp = "O:\\2011\\11\KeyResourceArea.shp"     Key_Resource_Separation_Area_shp = "O:\\2011\\1121\SeparationArea.shp"     Message = ""      This is what I want to fill...in python it would be just Message = arcpy.GetMessages() but what does the same thing in model builder?               Email = "georgec@xxxxxxxx"     Success = ""      How do I fill this with true/false depending on if the merge tool is successful or not?               KeyResourceArea_Merge = "C:\\Users\\georgec\\Documents\\ArcGIS\\Default.gdb\\KeyResourceArea_Merge"          # Process: Merge     arcpy.Merge_management(<snip>)          # Process: EmailStatus     arcpy.gp.toolbox = "D:/SCRATCH/Projects/106/Toolbox.tbx";     # Warning: the toolbox D:/SCRATCH/Projects/106/Toolbox.tbx DOES NOT have an alias.      # Please assign this toolbox an alias to avoid tool name collisions     # And replace arcpy.gp.EmailStatus(...) with arcpy.EmailStatus_ALIAS(...)     arcpy.gp.EmailStatus(Message, Email, Success)      Best,


ArcMap label graphic

I converted an ArcMap VBA code to a Visual Studio dll. In doing so I had to reference a new namespace "stdole" so my graphic label would show in ArcMap. The label works on my computer, but does not work on other computers that do not have VS 2010. I'm guessing it is because the namespace stdole.dll is not on the other computers? Are there any alternative namespaces that I can reference that the older computers might have? Below is the VBA pre and raw VB.NET post code conversion (<<<< this defines the code statement where the stdole namespace was needed). **VBA version** 'Get the point location               Set pPoint = pFeature.Shape     pTextElement.Text = strLabel 'Add the text string to the Text Symbol          Dim myFont As IFontDisp     Set myFont = New StdFont     myFont.Name = "Courier New"     myFont.Size = 1          pTextElement.Symbol.Size = 5#   'Added Text Symbol Size     pTextElement.Symbol.Font = myFont     pElement.Geometry = pPoint     pGraphicsContainer.AddElement pTextElement, 0      **VB.NET version** 'Get the point location               pPoint = pFeature.Shape     pTextElement.Text = strLabel 'Add the text string to the Text Symbol      Dim VB6 As Object VB6 = Nothing Dim myFont As System.Drawing.Font               myFont = System.Windows.Forms.Control.DefaultFont.Clone()     myFont = VB6.FontChangeName(myFont, "Courier New")     myFont = VB6.FontChangeSize(myFont, 1)          pTextElement.Symbol.Size = 5#      pTextElement.Symbol.Font = VB6.FontChangeName(pTextElement.Symbol.Font, myFont.Name) '<<<< stdole namespace needed     pElement.Geometry = pPoint      pGraphicsContainer.AddElement(pTextElement, 0)


How do I keep my callout boxes from looking wrong when the data frame is rotated?

I sometimes like to rotate a data frame in a layout, for example to fit more on the page if my area of interest is rectangular and runs from, say, northeast to southwest. But when I do, something strange happens to any labels I have which are in the callout box style: the carat pointing from the box to the feature gets very long and thin, and it looks terrible. Here's what my labels look like when the data frame is not rotated: ![enter image description here](http://i.stack.imgur.com/X6Yrp.png) I can convert that to annotation and move it around and the carat still looks like it should: ![enter image description here](http://i.stack.imgur.com/YO4rQ.png) But when the data frame is rotated, this is what it looks like, the carat is squashed down to a thin line. The only change I've made here is rotating the data frame: ![enter image description here](http://i.stack.imgur.com/qCG3Y.png) Converting it to annotation and moving it around doesn't improve things much: ![enter image description here](http://i.stack.imgur.com/Zc9Tu.png) Any ideas as to why this happens, and how to prevent it?


How to Receive an Email after a script has run?

I'm working my way through my first attempt at scheduled tasks with modelbuilder and python. So far I have:   1. Created a model that reconciles our SDE versions and then compresses the database.   2. Exported the model to a script so I can set it up in Windows Task Scheduler Its a fairly straight forward model so I expect it to work, however, I'd like to set it up to email me when it is completed and include the geoprocessing messages. I found this support article that includes a script that does just that. My question is: How do I include this or reference this in my reconcile and compress model?


How to get SOAP API Access on an ArcGIS Online Map?

I published a tiled map service to ArcGIS Online. This web site reported that we would be able to access the ArcGIS Online maps via the SOAP API from our applications. I check my first published map and under **Supported Interfaces** , I have REST but not SOAP options. Does anyone know how I can can get SOAP access to my ArcGIS Online Map? Thanks.


ArcSDE 10.1 Spatial View with no records

I can create a spatial view and it renders fine in catalog and ags. But if i modify the sql OR if a situation arises that makes my view devoid of records then all of the sudden sde and ags get really 'stupid' and have no clue what my geometry type or coord sys are for that view. Am i to believe that a spatial view MUST always have at least one record in it? I work with data that expires after a certain time. So as time elapses records become removed from a view such that at some points in time there will be no data that satisfies the view. When that happens everything goes downhill.


How do I hyperlink a to bookmark on pmf or mxd?

Can anyone advise how do I do a hyperlink that directly link to a bookmarked area of a pmf or mxd? Something like what we do to hyperlink to a specific page in pdf. For example, I have a mxd call BD_WORLD with Regional_View as one of the bookmark. SO something like this will be my hyperlink: \\\serverABC\GIS\WORLD_MXD\BD_WORLD.mxd#Regional _View Preferably, I would like to do it in attribute table, if possible. Thanks PS: using 9.3.1, with only Arcview, Spatial Analyst and Publisher.


Proximity Toolset Questions - Description, Examples of Usage, and Limitations

I was wondering if you could help me answer a few questions about the Proximity Toolset. I have been given the task of describing each tool in laymen terms, the typical usage of each tool, and most importantly, **the limitations** of each tool. So far, these are the only limitations that I have found: **1\. Buffer Tool** : Buffering across large regions and datasets, display performance issues when drawing or editing buffers **2\. Thiessen's Polygons** : The polygons are not created based on event importance and there is no facility to weight each point. Polygons can be impacted if the topography of the land changes across the study area. **3\. Generate Near Table Tool and Near Tool** : Point Distance will not work if either coverage contains more than 130,000 points


openlayers get feature info in variable

i m using openlayer and geoserver2.1.3 i have setup proxy issue but what i get               function handler(request) {     alert(request);         alert("respnse : " + request.responseXML);             alert("text : " +request.responseText);             alert("header : " + request.getAllResponseHeaders());          }          var request = OpenLayers.Request.GET({     url: url,     callback: handler     });             from response.text is the php proxy file.. i want to get data of url in variable.Any ideas


Help with Routing..RouteEventArgs cannot be serialized

Hi I created a routing service in ArcGIS Server that allows the Service Area task. I've hit the REST service for this routing service and passed in a single facility XY and it works perfectly fine. The trouble I'm having is in the Silverlight Geocortex application. I'm using the RouteTask activity which simply ask for the Routing service URL (I provided: `http://server/ArcGIS/rest/services/Water/WaterRoadNWwithService/NAServer/ServiceArea/solveServiceArea`), route parameters (I provided multiple variables like facilities, out spatialreference, etc), and lastly a result (which should be of RouteEventArgs). Whenever I run the task, I keep getting the below error and I'm not sure why or what to do next.. Thanks for any advice. ![enter image description here](http://i.stack.imgur.com/78cOo.png) Type ESRI.ArcGIS.Client.Tasks.RouteEventArgs cannot be serialized. Consider making it with the DataContractAttribute attribute and marking all of its members you want serialized with DataMemberAttribute attribute


Case Insensitive WFS Attribute Search with ArcGIS Server?

I'm using ArcGIS 9.3.1 to host WMS and WFS services. I have a map (.mxd file) with a layer whose attribute table includes fields with mixed-case string data. When I open this map in ArcMap and perform a Find (Edit -> Find), selecting the appropriate layer and field, I get all matching results back, no matter what case I use in entering the search criteria. Searching for 'fireside', 'Fireside', and 'FIRESIDE' in STR_NAME all return features whose STR_NAME is 'Fireside Ln'. However, if I take the same map and publish it to WFS in ArcCatalog, when I perform a WFS request against it, I only get matching results when searching for 'Fireside' (not 'fireside' or 'FIRESIDE'). I am performing a WFS request in the form of an XML POST. This is an example request that returns no results:               <wfs:GetFeature service="WFS" version="1.1.0" maxFeatures="1000" resultType='results' xmlns:wfs="http://www.opengis.net/wfs" xmlns:ogc="http://www.opengis.net/ogc" xmlns:gml="http://www.opengis.net/gml">      <wfs:Query srsName="EPSG:4326" typeName="centerlines:street_centerline">       <ogc:PropertyName>STR_NAME</ogc:PropertyName>       <ogc:PropertyName>TYPE</ogc:PropertyName>       <ogc:PropertyName>Shape</ogc:PropertyName>       <ogc:Filter>        <ogc:PropertyIsLike wildCard="*" singleChar="." escape="!">        <ogc:PropertyName>STR_NAME</ogc:PropertyName>        <ogc:Literal>*fireside*</ogc:Literal>       </ogc:PropertyIsLike>       </ogc:Filter>      </wfs:Query>     </wfs:GetFeature>      I have also tried adding a matchCase="false" attribute to the ogc:PropertyIsLike tag (based on a discussion I found on osgeo.org), but this also returned no results. Is there any way to perform a case-insensitive attribute search against a WFS service hosted on ArcGIS Server, either through modifying the request, or by performing some kind of configuration of the layer, map, or service on the server end?


How to get polygon points for all regions inside a city polygon?

Is it possible for me to get all polygon points for all regions inside a city polygon? I am building a website where i need city level polygons as well as polygons for all sectors and estate areas in the city. We tried using the getbyarea api of wikimapia but the box only allows to specify two points (2 pairs of lat long) and not the complete city polygon. This way a lot of regions in the city get missed. Are we doing something wrong or missing a point?


Creating Long & Lat for Polygons

I have a base map image with micro regions, however in the attribute table there is no information regarding the long & lat. Is there a way I can create a long & lat for each micro region? Thanks.


Is there an Offset function in QGIS for Polygons?

IS there an offset function in qgis similar to that in autocad for polygons. I have an irregular shaped pollygon which i want to offset at 100m out from the edges, effectively adding 100m out each time. I can do it buffering but does not show the exact shape.


Composite bands produces a single band image

Running this produces a single band tif and not a three band image. Has anyone seen this before or can suggest ways to trouble shoot it?               Input_Rasters = path + parentdir + "/" + dir + "/" + "Rt_b1dp.TIF"; path + parentdir + "/" + dir + "/" + "Rt_b2dp.TIF"; path + parentdir + "/" + dir + "/" + "Rt_b3dp.TIF"; path + parentdir + "/" + dir + "/" + "Rt_b4dp.TIF" # provide a default value if unspecified          # Process: Composite Bands     print "Process: Composite Bands"     arcpy.CompositeBands_management(Input_Rasters, "Rt_CompBands.tif")


Polygon label repeated for each tile

I'm using the standalone GeoWebCache to serv tiles from a remote GeoServer. My problem is that the polygon label is added to each one of the tiles served, instead of only once in the polygon centroid. I found a post which discusses the issue:   http://osgeo-org.1560.n6.nabble.com/polygon-label-repeated-for-each-tile- td4995203.html The first repley mentuned a possible solution: > "All in all, I suggest to use a tile rendering engine (GeoWebCache, > MapProxy, TileCache) anyway, instead of requesting small image from > GeoServer and have the tile rendering engine do the tile slicing afterwards. > You will have send fewer requests to GeoServer (1 large image instead of > multiple small images), so this speeds up the overall tile cache creation > time." Problem is that I couldn't find how to do that by reffering to the GeoWebCache documentation, and the above mentuned post doesn't explain the way to implement that. I also found a post with an answer that links to the GeoWebCache "Tiled" documentation, but my code allready uses all the necessary attributes and still the label shows up multiple times:               var Layer_1874 = new OpenLayers.Layer.WMS(         'Grundkort',         '/wms10.ashx'         , {             format: 'image/png',             srs: 'EPSG:25832',             layers: 'ballerupkommune_grundkort_bk',             tiled: true,             tilesOrigin: '698804,6173460'         }         , {             displayInLayerSwitcher: true,             isBaseLayer: true,             transitionEffect: 'resize',             displayOutsideMaxExtent: true,             visibility: false         }     );      Anyone has an idea? Kind regards   Amit


Looking for simple software to compare two Landsat Images

At uni I used ERDAS imagine to view two rasters at the same time and it had various viewers for comparing them. I'm looking for something that has a slider so that I can see the differences between two Landsat images (same area different dates) or if not a slider some good functionality for visually comparing them. ArcMap is too slow at switching between the two. Is there anything free that I could use for this?


How to use PostGIS for tile-based aggregation?

How can I create a tile-based aggregation of polygons with PostgreSQL/PostGIS? The input of such an operation could be the database table with the geometry and the data to aggregate, a bounding box, and the tile size and the method which should be used for the aggregation (sum, average, minimum, maximum). The output would be a list of geometries with the outlines of the tiles and the aggregated values. I found this operation in Oracle Spatial but could not find a solution on how to do this with PostGIS.


Aligning two layers captured in different coordinate systems

I have two layers, one of points (locations of fruit trees) and one of Atlanta zip codes, in an ArcMap, but the layers are not aligning. The zip codes are in the top right and the fruit points are tiny in the bottom left corner. ![enter image description here](http://i.stack.imgur.com/bJp7f.jpg) I successfully converted the fruit points to a shp file from a csv file. The fruit points in the csv are decimal coordinates taken from a GPS. The zip code map is in NAD83. **I should have also included that I set the fruit projection to NAD83 (the same as the zip code map). If I understand correctly, this is what is being suggested to me -- but I have already done so. The fruit is a shp layer and the zip codes are as well. How would I go about aligning these layers so the points fall correctly over the zip code map? Using ArcGIS10. Thanks in advance.


Exporting 3D geometry from PostGIS 2.0

I have a table containing 3D geometry in a PostGIS 2.0 database (the geometry is of type `LINESTRING Z`). I need to get the geometry out and into a DXF, DGN or 3D Shapefile. What is the easiest way to do this? (I was sure I had achieved this via QGIS before, but it seems to remove the `Z` value and produce only a 2D file when attempting a 'Save As ...' from a PostGIS layer).


Represent joined data in Info window (ArcGIS JSAPI 3.1)

Currently I am facing following problem I have joined the SQL Server data base to my shape files (In ArcGIS Desktop using “join” option)and published on ArcGIS Server and working on REST service in JavaScript API. After joining the tables to the map… all columns are changed. Let us say info.id, attr_name.ratio etc. Therefore, when I am using info window option I cannot use “info.id” it is considering only “info” it gives an error if I used “info.id”. Any help will be great... Software’s: ArcGIS DT 10.1 and Server 10.1 Code Sample:               'var featureLayer = new esri.layers.FeatureLayer("http://<myserver>:6080/arcgis/rest/services/type/MapServer/1",           {          // mode: esri.layers.FeatureLayer.MODE_SNAPSHOT,           mode: esri.layers.FeatureLayer.MODE_ONDEMAND,           infoTemplate:popupTemplate,           outFields: ['info','attr_name']           }); '


Copying representation symbology from one ArcSDE version to another?

We are having problems with our representations when copying from an ArcGIS 9.3 database to a 10.0. Representation symbologies with dash patterns (dots, dotted lines etc.) are reset to default settings in the version 10.0 database. The problem is the same when you copy from a 10.x to 10.x database. Have anyone experienced the same, and what did you do to make it go away?


How to call GRASS modules in pyQGIS?

I wish to use a Grass module (v.to.db) in QGis, without the grass toolbox. I tried with the Python console, but to no avail :   * Using Sextante plugin but obviously it didn't know all the grass modules. > from sextante.core.Sextante import Sextante > > Sextante.alghelp("grass:v.to.db") > > \--> Algorithm not found   * An other test :  > grass.run_command("v.to.db", map='bl@PERMANENT', layer='2', option='start', > units='meters', columns='X,Y,Z') > > \-->Traceback (most recent call last): > > File "", line 1, in module > > NameError: name 'grass' is not defined Do you know if there is an other way to use the grass module ? I don't want to use it in the grass toolbox because I can select only one colums in "attribute field", and I would like to choose several columns. Thank you for your help Marie


How to custom the LayerSwitcher in OpenLayers to toggle all layers?

I want to have a checkbox in my layer switcher "turn all layers on/off". how can I do that ? I thought about adding an empty layer with that name and somehow bind her an event "On Show" and "On Hide" to toggle the rest of the layers. Any better solution ?


how can I erase a "line feature" (erase feature) from a "polygone feature" (input) in ArcGIS 10? With the tool "erase" I can only use an other "polygon feature" as erase feature for a polygon (input).But I can not use an "line feature" as erase feature for a polygon. I there an other tool in ArcGIS or an other way to solve this problem? Thank you!


Adding Attachments in ArcGIS 10.0 (problem - can’t find the instrument Enable Attachments)

I want to create attachments to shp themes in ArcGIS (doc\pdf files, web- links). In the Help of ArcGIS I found the description of this operation: Adding attachments to features Add Attachments Tools (Data Management) But before attachments can be added using this tool, they must first be enabled using the Enable Attachments tool. So, I have to find Enable Attachments tool in the Data Management Tools. But I don't have this instrument in my ArcGIS Data Management Tools, so can’t add attachments (I tried to do it without use Enable Attachments tool, but in this case I just don't have "Open Attachment Manager button" which is needed for add attachment)). Could you advise me, how I can find the instrument “Enable Attachments tool”? Maybe I have to turn on something in the “properties”? Or this instrument has another name? ![enter image description here](http://i.stack.imgur.com/kGOFa.jpg)


Geowebcache does not displaying map tiles in EPSG:900913 from a remote WMS server

I connect Geowebcache to a remote WMS service through their Getcapability request by adding the url in geowebcache-servlet.xml document. Geowebcache demo page displays the list of wms layers from the remote wms server. When I click on the png link for the EPSG:900913 format , to display the map, only displays the image place holders instead of map tiles. This displays map for the png link for EPSG:4326 formaT. When I issue a Getmap request to the remote wms service directly from my browser, for EPSG:4326 works well but for EPSG:900913 it shows error “SRS Code 'EPSG:900913' is not supported (not found in SRS xml settings)”. When I check their Getcapability document I could not find EPSG:900913 in the SRS list, but EPSG:4326 is available. Since I want display the map in google map api v2 , think i need EPSG900913 format to get this work. I think geowebcache generate an automatic re-projection in EPSG:900913 format for googlemap api. This link is not working for me, I think that could be the reason this is not displying in google map api. could someone tell me a solution for this ? why the reprojection not working ? Here is the link for EPSG:4326 format , this will start disply at a scale 1:272K range 10th zoom This is the link for EPSG:900913 format which is not working for me. This is the link for the Geowebcache demo page I think there is a way to define independent layers and projections in geowebcache.xml document. Not sure how to do this to suit for my situation ?


Custom JSONDecoder for ID Task Result, ESRI Flex API

The problem: When the ESRI API decodes a response from the identify task, there is a few second lock in the Flash application which is due to the decoding of the results. I have used chunk parsing algorithms in action script to prevent the UI from locking, and I would like to implement this for the ID task response decoding. I assume the class I need to extend is the com.esri.serlization.jason.JSONDecoder class. Does anyone know how I can get that code so I can re-write it, and use this extended class to parse the results from the ID task? Or, if I use a third party JSON decoder and re-write it, how can I use this class to parse the results from the ID task? Thanks


Processing for large scale networking animation?

I am looking to create an animation of network connections, something very similar to this: http://vimeo.com/16897839. Let's say I have a huge database (millions of rows), that includes > source location, destination location, timestamp , and I want to create a timelapse animation of all the interactions going on, is it possible to animate this on the fly, or does it have to be saved frame by frame and then rendered into a video? I have been playing around with Processing, and it seems that animating anything > 5000 objects, the animation becomes too laggy to view in real time. I wanted to originally allow the viewer to select certain settings, and have the animation change depending on the filters, but I'm assuming if I have to render the image and record as video than that is not possible.


GRASS GIS: red square around node topology information?

I've been using GRASS since 2005, but I cannot find an explanation to the square boxes that appear around the nodes of a lines layers, when I display topological informations on the map. ![enter image description here](http://i.stack.imgur.com/elyvP.png) I thought they could indicate errors, but my layer appear to be topologically ok. Someone has told that those square lines refer to bboxes, but it doesn't convince me. Nodes' BBOXes? I have many layers and I only see them on this one... Any GRASS user/dev out? Giovanni


How to update the length property of a feature class field?

I have a feature class with a string field that I want to update the length from 10 to 25. Is there a way or a tool to update this property without having to create a new field and using field calculator to move over the records from old to new field? **Update** This script below does what I want, the only thing is it pushes the updated field to the end of the table (not keeping the original field order). http://arcscripts.esri.com/details.asp?dbid=16503 Another option that seemed promissing is to use the Feature Class to Feature Class tool (access this tool by right clicking a layer in ArcCatolog and selecting Export>to Geodatabase single option) . It does create a new layer, however you can update field names and property information. It also appears that you can change the field order using the "Move Input Field Up/Down" arrows, but they seem to not work. There has to be a solution to just edit field properties without having to create a new field? Thanks


Change Extents on Overviewmap

Does anyone know how to change the extents on the overview map in Flex Viewer 2.5? I have a basemap of a single country and my overview map still wants to show me the world. Thanks


Override or add to [GeoExt.ux PrintPreview] rotation

I'm using suite-sdk (GXP) to build my WMS client. I want add to PrintPreview window capability of changing rotation. I find this feature in Heron script but there is other version of [GeoExt.ux PrintPreview] link. So I add               printRotationPage: null,     printRotationExtent: null,     printRotationExtentOptions: null, and               if (this.showRotation) {         this.printRotationPage = new GeoExt.data.PrintPage({             printProvider: this.printProvider         });         this.printRotationExtent = new GeoExt.plugins.PrintExtent(Ext.applyIf({             pages: [this.printRotationPage],             layer: this.initialConfig.layer         }, this.printRotationExtentOptions));     } so file look like this. Now when I click print icon on toolbar I get this in console:               TypeError: this.map is null                                                       Layer.js         resolution = resolution || this.map.getResolution();      from ` getViewPortPxFromLonLat();` function from Layer.js file.   Where I made mistake?


Styling OS Vectormap Local in QGIS?

I am attempting to style Ordnance Survey Vectormap Local on QGIS 1.8.0 by using the styling sheets suplied by Ordnance Survey. The styling sheets are available from the OS website (http://www.ordnancesurvey.co.uk/business-and- government/help-and-support/products/vectormap-local.html). However when I load one of the Vector layers - take 'road centrelines' as an example and then load an SLD (OS VectorMap Local - Road_Centre_Lines.sld) in the style tab, it only displays two rules which don't style the layer. It's as though QGIS doesn't read the SLD correctly? Or am I not loading/applying the SLD files the right way? I was expecting the SLD to configure the symbology so that Motorways would be displayed as thick blue lines for example and minor roads to be displayed as small white lines. I have 16 Vector layers in total and it would be a very long process to manually style all of them so I would rather avoid that approach at all cost. I would be grateful for any suggestions. ![enter image description here](http://i.stack.imgur.com/VMuh6.jpg)


How to stop the opengeo suite services?

I cannot open dashboard. I check ports and found that suite geoserver is working. (I click several times on dashboard.exe) Please, tell me how to open dasboard.exe UPD: when i enter ...localhost:8080/geoserver/web/ in browser i go to the geoserver. It is works.


How do I automate custom edge labels for township range in ArcGIS 10?

I need to label township ranges along the side and top of my map, respectively. I know the grid in the frame properties is a way to label the edges with lat/lon. I Think there is a way to do the same type of thing with dynamic text objects. Anyone have experience with this?


kml does not overlay correctly in QGIS

In QGIS, when I insert a kml contour that was created in Google Earth, it appears misaligned in respect to all the other layers. Note in this image how the red contour is moved upwards with respect to where it should be (use the house and the road going south as reference). I am using EPSG:3003 as project CRS, all the different layers I have in my project from WMS services overlay correctly. The kml layer is set to EPSG:4326 (WGS 84) and I have set "automatic on the fly reprojection" in the options of QGIS. The CRS of the raster is set to EPSG:3003. I have QGIS 2.0.1 but the same problem happens in 2.4. ![example screenshot](http://i.stack.imgur.com/jBuiN.jpg)


insert row into postgis query table in QGIS

I find I am using QGIS (2.4) more and more as a client for PostGIS database tables with geometry, but also for tables without geometry, partly because it's good just to use the one application and partly because it has very good editing support in the form of widgets etc. One thing I don't seem to be able to do is to insert a new row in a particular table. The table is a "query table", but I have set the relevant trigger functions on the query to update the right tables and tested it with normal SQL (it works). When displaying the layer and having editing on, the "Add Feature" button is greyed out. Am I missing a setting or is this a bug or design decision?


GRASS 6.4.0 and GRASS 6.4.SVN

What is the difference between GRASS 6.4.0 and GRASS 6.4.SVN? Thanks


Is it possible to find the Map Services URLs used in ArcGIS Flex or Silverlight web applications?

With the ArcGIS Server JavaScript API, you can find the URLs in the map services in the source code. Do the Flex and Silverlight APIs completely obfuscate the map services URLs? As an example, the County of Ventura has this web application in Flex. However, there doesn't seem to be a way to find the rest services URLs for the map services that one could use in their own applications/mashups. Perhaps obfuscation is the point to using a compiled web application vs the JavaScript API?


How to make the arrows in d.rast.arrow GRASS script visible?

I'm trying to create a Flow or directional raster map  (with arrows) from an SRTM elevation map that I downloaded from CGIAR. I already created an aspectMap to be used for the d.rast.arrow. I used the script: `d.rast.arrow map=aspectMap type=grass` and it only displayed a map that is solid green without any arrows. I also tried the script: `d.rast.arrow map=aspectMap type=grass skip=50` and it displayed dots on the screen for every 50th cell. I already zoomed in to its max, but still I cannot see the arrows. Is there a problem with my script? or is there an alternative way to create the arrows?


ArcGIS Layout View 9.3

I have a question if you are able to, upon rotating the map, is there a way to do the same with the labels? If so, how? Also, besides using the "Rotate Data Frame" to rotate the Layout View, which other method can you do to rotate the layout? I right-click but the rotate options are ghosted, meaning those options are non-clickable. I have ArcView 9.3. Thank You. Z


Is there a performance hit from using a Feature Dataset?

I vaguely remember being taught not to use a Feature Dataset to make things "look neat" in the database, and that the Feature Dataset should only be used when there was some sort of spatial relationship between the component feature classes. The help file lists "organize thematically related feature classes" as a valid reason for using a feature dataset. Is it a bad idea to use a feature dataset to group layers in the geodatabase, when the only purpose is to help navigation in ArcCatalog? Is there a performance hit from storing feature classes in a feature dataset?


Total newbie. Creating map from aerial photo of golf course in QGIS

I have a aerial photo in tif with a tfw file of the golf course that I work at. I would like to create vector data of different features on the course such as greens fairways tees, etc. to calculate areas and distances to start with. So far I have been able to load the tif into QGIS but after that I am completely lost. It seems that every tutorial is using data sets that are complete. I imagine I am creating new data using the raster data? I'm lost. Thanks for any help


Set opacity of cached ArcGISTiledMapServiceLayer

Is it possible to set the opacity of a cached Tiled Map service using ESRI JS API? I am using:               var layer;          layer = new ArcGISDynamicMapServiceLayer("URL",                 {"opacity": 0.5});          map.addLayer(layer);      This should work with 50% transparency being set to the layer, right?


Importing SRTM Data into Grass

I am having issues importing an SRTM HGT file into Grass. I first made a new location as the instructions from this Grass Newsletter Tutorial states. I then went to the download site and downloaded the file to my desktop. Inside GRASS, I went to File > Import raster data > SRTM HGT import [r.in.srtm] I choose the filename of the .hgt.zip file and remove the extension .hgt.zip (as the instructions say). The statement looks like this:               r.in.srtm input=C:\Users\Joe\Desktop\N47E007        The following is the output that I'm getting:               Extracting 'C:\Users\Joe\Desktop\N47E007'...     Converting input file to BIL...     mv: cannot stat `C:\\Users\\Joe\\Desktop\\N47E007.hgt': No     such file or directory     ERROR 4: `C:\Users\Joe\Desktop\N47E007.bil' does not exist     in the file system,     and is not recognised as a supported dataset name.     rmdir: C:\Users\Joe\Documents\grassdata/demolocation/PERMANE     NT/.tmp/8120.0: Directory not empty     ERROR: Unable to import data      At the beginning of the script, the file is definitely there because the program recognizes it and starts converting it. Once it produces an error, I notice the file has been moved to another directory               C:\Users\Joe\Documents\grassdata\demolocation\PERMANENT\.tmp\8120.0      What can I do to fix this problem? Thanks!


Using different colors on Torque CartoDB

We’re creating a new piece and we’re having some limitations using the Torque function. Could you answer these questions, please?   1. Is it possible to insert multiple colors on this function?   2. Is there any way to change the size of the bubbles on Torque (min-max)? I have a programmer guy on my side and he’s trying to put these new features on the css code below               /** torque visualization */          Map {     -torque-frame-count:128;     -torque-animation-duration:30;     -torque-time-attribute:"n_de_participantes";     -torque-aggregation-function:"count(cartodb_id)";     -torque-resolution:2;     -torque-data-aggregation:"cumulative";     }          #testando_o_torque {        marker-width: 12;        marker-opacity: 0.9;        marker-allow-overlap: true;        marker-placement: point;        marker-type: ellipse;        marker-line-width: 2;        marker-line-color: #ffffff;        marker-line-opacity: 1;     }          #testando_o_torque[n_de_participantes=500] {        marker-fill: #A6CEE3;     }     #testando_o_torque[n_de_participantes=150] {        marker-fill: #1F78B4;     }     #testando_o_torque[n_de_participantes=20] {        marker-fill: #B2DF8A;     }     #testando_o_torque[n_de_participantes=600] {        marker-fill: #33A02C;     }     #testando_o_torque[n_de_participantes=70] {        marker-fill: #FB9A99;     }     #testando_o_torque[n_de_participantes=60] {        marker-fill: #E31A1C;     }     #testando_o_torque[n_de_participantes=200] {        marker-fill: #FDBF6F;     }     #testando_o_torque[n_de_participantes=2000] {        marker-fill: #FF7F00;     }     #testando_o_torque[n_de_participantes=300] {        marker-fill: #CAB2D6;     }     #testando_o_torque[n_de_participantes=6] {        marker-fill: #6A3D9A;     }     #testando_o_torque {        marker-fill: #DDDDDD;     }     #testando_o_torque[frame-offset=1] {        marker-width:12.5;        marker-opacity:0.45;      }     #testando_o_torque[frame-offset=2] {        marker-width:14.5;        marker-opacity:0.225;      }


database-backed WMS service runs as transaction

I'm trying to tame a geoserver 2.3.2 installation (backed by Postgres 9.1 and Postgis 1.5) by determining which of the many, many layers are actually being used. To this end, I've set up an audit table that logs SELECT activity on the tables I think might be deprecated. However, when I make a WMS request to a layer that should be audited, I don't get any records in my audit table. Of course, auditing SELECT in PostgreSQL isn't foolproof because INSERTs into the audit table are susceptible to ROLLBACKs. Indeed, when logging all queries from geoserver to postgres, I see that geoserver is querying data in a transaction and then rolling back               bind S_1: BEGIN     execute S_1: BEGIN     parse <unnamed>: SELECT ...     bind <unnamed>: SELECT ...     execute <unnamed>: SELECT ...     parse S_2: ROLLBACK     bind S_2: ROLLBACK     execute S_2: ROLLBACK      There doesn't appear to be any reason to run a SELECT inside a transaction, and even less of a reason to ROLLBACK the transaction as there are no errors. Is there a setting to turn off this transactional behavior when responding to a WMS request? I like the database audit solution because it is easily searchable and reportable. As a work-around I could certainly just turn on verbose logging at either the database or geoserver level and grep the logs for layer names.


The Info Windows content change for Google Maps Engine cannot be adopted

When I edit the Info Window, I modified `<a href="http://cs.sj.gov/pub/prjDetail.cfm?prj_id=+{PRJ_ID}"> View Details </a>` as `<a href="http://cs.sj.gov/pub/prjDetail.cfm?prj_id=+{PRJ_ID}" target="_blank"> View Details </a>` However, in the map preview, the new attribute `target="_blank` is not added into the Info Window. Did I miss anything?


How to split multiple polygon layer into multiple layers by vertices parts?

I have a layer in Arcmap 10 that was brought over from Auto CAD (layer A). This layer has multiple polygon shapes. I need to split these shapes from layer A, and place them into their own layers. When I brought over the layer from Auto CAD, it put all the polygons into one shape. So when I go to select one in editing, it selects them all. When I look at editing the vertices I can see each one is a separate part. Is there a way to split these parts apart to place into their own layers?![Note that there is only one shape in the table, and each part is linked to a specific polygon](http://i.stack.imgur.com/uEIhz.png)


How can you change the bubbles size on Torque in CartoDB

Is there any way to change the size of the bubbles on Torque (min-max)?


How to extract data from a POINT column in PostGIS using JAVA (jdbc)

How do you get a value from a column which has a `point` data type. I've just switched over from MySQL to PostGIS. I have a column which is called point and is of type `Point`. I want to get that value as a result of a SQL query. I've done things like `rs.getString("text")` and `rs.getInt("number")` for other columns and they work fine, but doing               rs.getObject("point")      Causes me to get the error               org.postgresql.util.PSQLException: This connection has been closed.      Here's more code, I thought it might a real-dummies question so I didn't stick it in before. Essentially I have a value in my application which can take on a variety of data types. Not having a truly OO database, I've made a table with lots of columns, depending on which is is populated, I will know what datatype the value is for that record.               public static Set<STriple> readStatements(String subUri, String preUri) {         Set<STriple> statements = new HashSet<STriple>();         String sql = "SELECT sub_id, sub_uri, sub_value, pre_id, pre_uri, pre_value, obj_id, obj_uri, obj_value, number, text, file, url, point FROM triples WHERE sub_uri='"                 + subUri + "' AND pre_uri='" + preUri + "'";         ResultSet rs = Dao.executeQuery(sql);         try {             while (rs.next()) {                 STriple s = new STriple();                 s.setSubject(new SResource(rs.getString("sub_uri")));                 s.setPredicate(new SResource(rs.getString("pre_uri")));                 if (rs.getObject("obj_id") != null) {                     s.setObject(new SResource(rs.getString("obj_uri")));                 } else if (rs.getObject("point").toString() != null) {                     Object pointObject =  rs.getObject("point"));                 } else if (rs.getObject("number") != null) {                     s.setObject(new SNumber(rs.getDouble("number")));                 } else if (rs.getObject("text") != null) {                     s.setObject(new SText(rs.getString("text")));                 } else if (rs.getObject("url") != null) {                     s.setObject(new SURL(rs.getString("url")));                 }                 statements.add(s);             }         } catch (SQLException e) {             e.printStackTrace();         }         return statements;     }


Error message - AdvanceToElement:unable to find xml element

How can I handle the below error when I'm trying to add an XML file into ArcGIS for Desktop using ArcCatalog. > AdvanceToElement:Unable to find XML element (WorkspaceDefinition not found > in XML doc) XML files are provided by the client as a reference files to attach in our ArcMap, this is fist time we are using XML in ArcMap. Our ArcMap version is 10.1. Please, remember that I'm very new to using xml files in arcmap.                   <?xml version="1.0" ?>      - <INVENTORY>     - <ATTACHMENT>       <REF_TYPE>NONE</REF_TYPE>        <ANCILLARYROLE>0</ANCILLARYROLE>        <ATTACHMENT_NAME>SLACK:FIB::961</ATTACHMENT_NAME>        <ATTACHMENT_REF_NAME>20m</ATTACHMENT_REF_NAME>        <CATEGORY_NAME>SLACK-E</CATEGORY_NAME>        <ENABLED>1</ENABLED>        <INVENTORY_STATUS_CODE>IPL</INVENTORY_STATUS_CODE>        <LENGTH_ADJUSTMENT>20</LENGTH_ADJUSTMENT>        <MODEL_UUID>{2937E405-1E38-42F9-991B-4AAA4A6DED77}</MODEL_UUID>        <MODIFICATION_IND>0</MODIFICATION_IND>        <Q_CREATE_WORK_ORDER>TCA_10_4_base_data</Q_CREATE_WORK_ORDER>        <REF_UUID>{5E60CBAE-C234-4A1A-B132-4CAC48AF1FC4}</REF_UUID>        <ROTATION>0</ROTATION>        <GEOMETRY>X="230483.90725" Y="395726.049"</GEOMETRY>        <TRANSMEDIA_NAME>FIB-Q:48F::6406</TRANSMEDIA_NAME>        <TYPE_NAME>FIBER</TYPE_NAME>        <UUID>{ACBBDD7D-1DA1-4F75-84BB-803DB4B04C9B}</UUID>        </ATTACHMENT>     - <ATTACHMENT>       <REF_TYPE>NONE</REF_TYPE>        <ANCILLARYROLE>0</ANCILLARYROLE>        <ATTACHMENT_NAME>SLACK:FIB::962</ATTACHMENT_NAME>        <ATTACHMENT_REF_NAME>20m</ATTACHMENT_REF_NAME>        <CATEGORY_NAME>SLACK-E</CATEGORY_NAME>        <ENABLED>1</ENABLED>        <INVENTORY_STATUS_CODE>IPL</INVENTORY_STATUS_CODE>        <LENGTH_ADJUSTMENT>20</LENGTH_ADJUSTMENT>        <MODEL_UUID>{2937E405-1E38-42F9-991B-4AAA4A6DED77}</MODEL_UUID>        <MODIFICATION_IND>0</MODIFICATION_IND>        <Q_CREATE_WORK_ORDER>TCA_10_4_base_data</Q_CREATE_WORK_ORDER>        <REF_UUID>{5E60CBAE-C234-4A1A-B132-4CAC48AF1FC4}</REF_UUID>        <ROTATION>0</ROTATION>        <GEOMETRY>X="230583.2762" Y="395795.57475"</GEOMETRY>        <TRANSMEDIA_NAME>FIB-Q:48F::6406</TRANSMEDIA_NAME>        <TYPE_NAME>FIBER</TYPE_NAME>        <UUID>{6BB5FB19-B6A5-4DC0-8657-F5B1D76F1757}</UUID>        </ATTACHMENT>       </INVENTORY>      ![enter image description here](http://i.stack.imgur.com/XtJIj.png)


cartodb linked to sharepoint

has anyone tried linking a sharepoint excel spreadsheet with cartodb? my spreadsheet has a url like it would from google drive or dropbox and it is shared with editing and viewing privlidges for everyone (public) however it keeps erroring when I try to import error: Error 1001 Download error What about... The remote URL returned an error. Please verify your file is available at that URL. Any ideas? I would like to use excel as the program for my team to edit my data and I want it synced to cartodb.


How to perform in-line variable substitution for path in ArcGIS ModelBuilder?

I am having difficulties in finding a way to use the input path and set is as the output path for a simple Copy Features tool. I use the substitution by the workspace but it outputs the results into the default ArcGIS folder. I am not sure how to set the workspace to hold the path of the input feature class. As the input feature class is a parameter, I cannot use a hardcoded path. Please see below: ![enter image description here](http://i.stack.imgur.com/SUHaR.jpg)


Convert OpenStreetMap schemas to ArcSDE in PostgreSQL database

I am trying to find a way to store OSM Extracts (or possibly the planetfile) in a PostgreSQL database, keep it up to date with diff files and then process the data further with ArcGIS (e.g. Geocoding, Mapping). The best idea I have come up with so far is to import the osm data into the database with common OSM tools like Osmosis or osm2pgsql. My problem is, that those tools create different schemas than SDE. Has anyone ever tried this before, or a good idea on how to convert between the schemas? I'm open to any comments or suggestions.


copying and pasting files - geodatabase

I have been sent a zipped folder containing files in a geodatabase format. To view them in ArcMap I've created a new file outside ArcCatalog with the extension .gdb and copied the unzipped files into this new file. Outside ArcGIS I can see all the files are present but when I try to preview in ArcCatalog there are no files in the new folder. I've also tried to 'Add data' in ArcMap but still nothing and I've also tried doing the whole copy and paste thing from within ArcCatalog but still no joy. Can someone help please?? What am I doing wrong? Thank you!


Exporting the contents of an attribute table in formats other than an Excel file

I know that it's possible to export the contents of an attribute table to an Excel file, but could they be exported in any other format (from ArcGIS 10.0)? I have a feature class whose attribute table (minus a few fields) needs to be made available to be displayed and read by the human eye as an ordinary table. Could the attribute table be exported into (for example) .pdf format, or do I have to export it to Excel first, and start from there?


Highlighting selected feature in OpenLayers with additional graphic or reticule

I am trying to create an OpenLayers map in which I can highlight a selected point feature by putting another graphic around it, such as this red box: ![Reticule above](http://i.stack.imgur.com/f3Efc.png) What I have done so far is to add a vector layer with a single point feature, and a styleMap that displays the reticule as its externalGraphic. I then move the feature and display it when I want to highlight something, or hide the feature when I don't. This works reasonably well, but I currently have to choose between one of two options, but each has a disadvantage I want to avoid. **Option 1: Correct display, incorrect interaction**. The display is as above, exactly as I want it, but features "covered" by the transparent reticule can no longer be clicked. I achieve this by adding the reticule layer to the SelectFeature control and using setLayerIndex to move it the top. This means that hovering over or clicking on any other features within the red square has no effect. **Option 2: Incorrect display, correct interaction.** If instead I don't add the reticule layer to the SelectFeature control, I can successfully click on any of the selectable features to move the reticule to them. However no matter what I set the z-index of the reticule layer or its feature to, all of the features added to the SelectFeature obscure it, as in this image: ![Reticule below](http://i.stack.imgur.com/64odb.png) * * * So my question is, **how can I achieve the combined effect?** I want the highlighting reticule above all features, but to have no interaction with the mouse events. Is there a better way than defining this separate layer for the highlight image? If this is the best way to go, can I pass any click or hover events through the reticule layer to the feature layers below it?


Geoserver ImageMosaic: error connecting to WorldImage

I would like to use the ImageMosaic plugin to set up a Geoserver layer which refers to a WorldImage, just like in the official tutorial. As a simple test case, I put a custom folder into my Geoserver data directory with one World Image (image + WorldFile + projection). When I specify the folder location to a ImageMosaic Store in Geoserver, it simply displays the error message:               Could not list layers for this store, an error occurred retrieving them: Argument "value" should not be null.      On the other hand, when I put my image into a default folder that comes with the openGeo suite (like `file:data/dem/` with the Medform Digital Elevation Data), the thing works as it should. So I suspect there is something going wrong with user rights, can anybody enlighten me? * * * PS: My Geoserver error log                   INFO [geoserver.web] - Getting list of coverages for saved store file:mosaic         java.lang.RuntimeException: Could not list layers for this store, an error occurred retrieving them: Argument "value" should not be null.         at org.geoserver.web.data.layer.NewLayerPageProvider.getItemsInternal(NewLayerPageProvider.java:123)         ...         Caused by: java.lang.IllegalArgumentException: Argument "value" should not be null.         at org.geotools.util.SoftValueHashMap.ensureNotNull(SoftValueHashMap.java:140)         ...


Programmatical Setting of 'Store Relative Pathways' for multiple MXDs ArcGIS 10.1

I am looking to move large number of mxds to a different location on the same drive. Many mxds haven't got 'Store relative pathways' option checked therefore this operation will result in broken links to data. I've only started using ArcGIS 10.1 and it is now possible to programmatically update relative pathways for all mxds by using arcPy and setting relativePaths to be "True". I had a go at writing Python script which (in theory) loops through a top folder and searches for all mxds and then updates "Store relative pathways" setting to be on. When I run this script in PyScripter, it fails to complete, I get the following error: "exceptions.NameError: name 'mxd' is not defined Can anyone help ?!? Thanks in advance, Magda                   import arcpy, os          folderPath = r"C:\My Folder Path\...\...\..."          for filename in os.listdir(folderPath):         fullpath = os.path.join(folderPath, filename)         if os.path.isfile(fullpath):             basename, extension = os.path.splitext(fullpath)             if extension.lower() == ".mxd":                 mxd = arcpy.mapping.MapDocument(fullpath)                 mxd.relativePaths = 'True'             mxd.save()     del mxd


ArcGIS 10.1 Model cannot be transferred between Win7 and Win8 machines

![screenshot of attempt to "edit" model module](http://i.stack.imgur.com/TrGZz.png) We created a model (created in a file geodatabase) using Model Builder on a trial version of ArcGIS 10.1 (build 3035, License: Advanced). We are trying to now run it on a registered (build 3134, ArcMap 10.1, ArcGIS 10.1 for Desktop SP1, License: Basic, ) version of ArcGIS 10.1 on another computer, but when we try to open it by right-clicking and selecting edit on the model it displays only the outlines of the tool-shapes within the model. We can't open the different modules of the model or anything to be edited. All the model does is calculate fields, so I was not expecting the basic license to be an issue. The same issue occurs when creating a model with the Basic registered version and trying to open it on the trial version, which leads me to think that this is an issue of transferring between computers, rather than a license-related issue. I changed the permissions on the .gdb folder to "full control" for everyone on the network to rule out permissions as an issue. The datasets in the model itself can be opened, but the actual tools (such as "Field calculate") are cannot be accessed. Could it be that the model is trying to reference, for example, the "field calculate" tool on the host computer with the wrong path?? Any ideas what the issue is? Any help would be greatly appreciated!


GeoServer using GeoWebcache and SQL Server whats the process?

I'm using GeoWebCache and GeoServer to tile cache a layer of data from a SQL Server database. If I seed the cache and I send a tiled WMS request what occurs? Does GeoServer make a call to SQL Server or are the tiles taken straight from the GeoWebCache? If the underlying SQL server data changes how does GeoServer know to regenerate the relevant tile? Any help apprciated.


Symbology for Unique Values with Domain

I'm wondering if there is a workaround to this. I want to symbolize a feature class categorically by unique values. The field that I'm doing this on has about 5 unique values. I also have a coded value domain on this field with about 3000 values. When I try and symbolize the field, it creates a symbol for every value in the domain, not just for the values in the feature class. Is there something I'm missing that will limit the symbology to only the unique values that are in the feature class?


QGIS Processing - Additional Creation Parameter for GdalWarp

I am having a hard time getting the "Additional Creation Parameters" to work when using GdalWarp through the processing toolbox. I am trying to set the nodata value to 0 for the output raster.               processing.runalg("gdalogr:warpreproject","InRaster.tif","EPSG:32736","EPSG:32736",10,0,"-dstnodata 0","OutRaster.tif")      Has anyone used this functionality and managed to get it to work?


Possible to add multiple fields in single arcpy statement?

I'm trying to add multiple fields to a table using arcpy and am wondering if it's possible to have multiple fields declared in the AddField management, and their values, separated by commas or something? Here is the code as stands, which is slightly cumbersome:               import arcpy, sys          arcpy.env.workspace = "C:\\gislab2\\Python\\take_home"     shapefile = "USCancer2000.shp"          #Set fieldname variables     field_name1 = "Crdrt"     field_name2 = "Crdrt_min"     field_name3 = "Crdrt_max"     field_name4 = "Ageadj"     field_name5 = "MinAgeAdj"     field_name6 = "MaxAgeAdj"          ###Declare each standard population per age group, and for summed standard population     Pop1Stnd=78784000      Pop2Stnd=150486000     Pop3Stnd=45364000     TotalStnd=274634000                   #Add the fields one at a time          arcpy.AddField_management(shapefile,field_name1,"DOUBLE",     "5","4","","","NULLABLE","NON_REQUIRED","")         arcpy.AddField_management(shapefile,field_name2,"DOUBLE", "5","4","","","NULLABLE","NON_REQUIRED","")         arcpy.AddField_management(shapefile,field_name3,"DOUBLE", "5","4","","","NULLABLE","NON_REQUIRED","")         arcpy.AddField_management(shapefile,field_name4,"DOUBLE", "5","4","","","NULLABLE","NON_REQUIRED","")         arcpy.AddField_management(shapefile,field_name5,"DOUBLE", "5","4","","","NULLABLE","NON_REQUIRED","")         arcpy.AddField_management(shapefile,field_name6,"DOUBLE", "5","4","","","NULLABLE","NON_REQUIRED","")      I know there's a way to create a loop to do it, but was wondering if there was a simpler and more direct way.


Is there a universal map/ layer format?

Just wondering if there is a universal map format that all GIS software can read. ie to quickly share layers with all our partner agencies which have different GIS Software, rather than have to convert it in multiple different formats (from our end or their end). We use ESRI, but have partners using MapInfo, Blue8, CadCorp... I wondering If I was missing a trick.


Techniques for Filtering the Attribute Table in QGIS's Print Composer

For the legend of a trail map I'm working on I'd like to show the trails and their lengths. I have this data in the table. The problem is that not all of the features (lines) have a name or length. Some are ancillary trails that shouldn't be include in the attribute table in the Print Composer. Are there any suggested techniques for filtering the list or creating some kind of "virtual" vector layer that will filter the data and I can point my Attribute Table to that table? Sorting by name descending and limiting the number I'm showing is not really an acceptable solution. ![enter image description here](http://i.stack.imgur.com/XkHZO.png)


Do SLD styles work with KML in Geoserver?

I am trying to use and in the SLD description for a layer on GeoServer. When I download the KML version of the layer everything works just fine in google earth, however when I go to my webpage which uses a KML overlay the `<MinScaleDenominator>` and `<MaxScaleDenominator>` do not work any more, as well as title.ftl. Does anyone have any suggestions or patches that might fix this? Also I would prefer to do a KML overlay and not redo everything as an open layer. Using: Geoserver 2.4 Apache 2.2


how do I create a RESTful PostGIS web Service Framework?

I have PostGIS tables that store vectors that want to show on Google Maps. Read the post h`ttp://gis.stackexchange.com/questions/22354/how-to-display-a- postgis-table-on-google-maps` But one thing that was not clear to me how do I create a RESTful PostGIS web Service Framework for feed data? I have all my data published in GeoServer, I do with it?


How to convert .shp to .dgn (V8i) including attributes AND geometry?

looking for open source solution using GDAL / OGR solution looks more promising , even before someone recommended me but not not sure on how to use then QGIS 2.0.1 Dufour   1. SHP file Left / Right Label Alignment issue in QGIS 2.0.1-Dufour   2. Save label edited shapefile and convert to Microstation DGN


Check if Openlayer Layer is empty, caused by nonexistence of internet connection?

I use Openlayer Google street map as a base-map in my QGIS plugin projects. i want to check programatically if the layer is loaded or is empty and decide if i load my offline base-map instead of Openlayer Google street map. Can someone give me an idea how to check that?


Import GPS data (text documents files) to QGis

I want to import GPS data (text document files) in to Qgis, what file format are you suppose to use ?


Can Google Earth for Android open a local document from a .kml field?

I've been making a .kml file to be used in GE for Android that should open paths from the local device into Polaris Office. I'm using the following attribute tag: `<a href="storage\sdcard0\Download\pcr\1050.xls">Tessier Dr</a>` but not getting anywhere. Has anyone been able to make something likes this work?


importing osgeo module in ubuntu 10.04

I have problem with importing "osgeo" module in ubuntu 10.04; I have installed GDAL 1.7 from source ( 1)./configure 2)make) however, still in python I cant import _osgeo_ module. I appreciate any hints. Basically, I am trying to make the converter script of jvectormap script into work and osgeo module is a prerequisite module. **update 1** _./configure --with-python_ gives me the following               GDAL is now configured for i686-pc-linux-gnu            Installation directory:    /usr/local       C compiler:                gcc -g -O2       C++ compiler:              g++ -g -O2            LIBTOOL support:           yes            LIBZ support:              external       GRASS support:             no       CFITSIO support:           no       PCRaster support:          internal       NetCDF support:            yes       LIBPNG support:            external       LIBTIFF support:           external (BigTIFF=no)       LIBGEOTIFF support:        internal       LIBJPEG support:           external       8/12 bit JPEG TIFF:        no       LIBGIF support:            external       OGDI support:              no       HDF4 support:              yes       HDF5 support:              yes       Kakadu support:            no       JasPer support:            yes (GeoJP2=yes)       ECW support:               no       MrSID support:             no       MSG support:               no       GRIB support:              yes       EPSILON support:           no       cURL support (wms/wcs/...):yes       PostgreSQL support:        yes       MySQL support:             no       Ingres support:            no       Xerces-C support:          yes       NAS support:               no       Expat support:             yes       ODBC support:              yes       PGeo support:              yes       PCIDSK support:            internal       OCI support:               no       GEORASTER support:         no       SDE support:               no       DODS support:              no       SQLite support:            yes       SpatiaLite support:        no       DWGdirect support          no       INFORMIX DataBlade support:no       GEOS support:              yes       VFK support:               yes                 SWIG Bindings:          python             Statically link PROJ.4:    no       enable OGR building:       yes       enable pthread support:    no       hide internal symbols:     no      **update 2** after importing gdal, I am getting the following complains:               >>> import gdal     Traceback (most recent call last):       File "<stdin>", line 1, in <module>       File "/usr/local/lib/python2.6/dist-packages/GDAL-1.7.1-py2.6-linux-i686.egg/gdal.py", line 2, in <module>         from osgeo.gdal import deprecation_warn       File "/usr/local/lib/python2.6/dist-packages/GDAL-1.7.1-py2.6-linux-i686.egg/osgeo/__init__.py", line 21, in <module>         _gdal = swig_import_helper()       File "/usr/local/lib/python2.6/dist-packages/GDAL-1.7.1-py2.6-linux-i686.egg/osgeo/__init__.py", line 17, in swig_import_helper         _mod = imp.load_module('_gdal', fp, pathname, description)     ImportError: libgdal.so.1: cannot open shared object file: No such file or directory      I tried to see where the package is install :               $ /sbin/ldconfig -p | grep gdal         libgdal1.7.0.so.1 (libc6) => /usr/lib/libgdal1.7.0.so.1         libgdal1.6.0.so.1 (libc6) => /usr/lib/libgdal1.6.0.so.1         libgdal1.6.0.so (libc6) => /usr/lib/libgdal1.6.0.so


Change Values of Georeferenced Raster

I am using QGIS to georeference B&W images. I have now done this and would like to set the black surround of the image (as a result of the georeferencing) to transparent but this does not work as bands 1,2 and 3 have the same values. I have noticed there is now a band 4 after georeferencing and would use this to drop out the black surround but I can only do this to bands 1,2 and 3. I am not super experienced with rasters in QGIS so I am having difficulty determining if this question has been addressed. This is a screen grab of the image I have, I am trying to remove the black triangles: ![This is a screen grab of the image I have, I am trying to remove the black triangles](http://i.stack.imgur.com/m6d9P.png) These are the values of the black surround: ![These are the values of the black surround\[\]\[2\]](http://i.stack.imgur.com/f4Gjv.png) These are the values of the black text and lines I would like to keep: ![These are the values of the black text and lines I would like to keep](http://i.stack.imgur.com/6He2F.png) For anyone who reads this who might find this usefull I just worked out how to do it. I can set a band to be transparent. So I just set Band 4 to be transparent and it worked perfect. ![enter image description here](http://i.stack.imgur.com/uUYg6.png)


How to cluster spatial data in R?

_I am very new to R, so please bear with me : >_ I have bunch of data points with latitude and longitude. I want to use R to cluster them based on their distance. I have already taken a look at this page and tried clustTool package. But I am not sure if clust function in clustTool considers data points (lat,lon) as spatial data and uses the appropriate formula to calculate distance between them. I mean I cannot see how they differentiate between spatial data and ordinal data. I believe the distance calculation between two points on map (spatial) and two normal numbers is different. (Is it not?) Also what happens if I want to consider a third parameter in my clustering? Like say if I have (lat,lon) and one other parameter. How is the distance calculated? The other problem I have with clustTool is that it is designed with a GUI in mind. I don't know how I can skip the GUI overhead in the library because I don't need it. **So basically I want to know what options I have in R for cluster analysis of spatial data?** _I really appreciate any hint or suggestion._


ArcGis: Rotated grids and cell size

I am using ArcObjects, namely RasterGeometryProc, to apply rotation on grids. The transformation works fine. However, when looking at the grids properties in ArcGis, the cell sizes has changed. Without rotation, the cell size is 50 X 50 as expected. After applying rotation transformation (30 degrees), I get 72.728 X 64.539. Seems like the new size is the extent of the rotated cell ? Is this a bug, or is it the way it is supposed to be ?


ctypes: bridge in python for a dll

I have a dll and I'd like to call it in python. The dll needs of parameters. Do you know how to pass a radiobutton parameter and a checkbox parameter? thanks