title: R generated files getting scrambled on import in ArcGIS

text: I want to put this out there just in case someone somewhere has experienced similar issues. I have a .csv file that I have generated in R 3.0 using this code: `write.table(SubData,paste(RootOutputFileName,"TimeSubsetData_", Date,"_.csv", sep = ""), sep = ",", row.names = F)` When I bring in the resulting csv file into ArcGIS 10.1 it has scrambled the column order, renamed the columns, moved headings, removed data completely, replaced data with NA's, and all sorts of other odd things. Here's what the csv file looks like outside of ArcGIS:               ID             Date_Time    moose   trial   x      y    label   label.1 dog  Reproduction   Position    EvalPt  DistX   DistY   Dist    Heading TimeDif Velocity    Angle   LocType LocNum     3311_1  2010-10-07 00:20    F07011  1   1476024 6620118 before  before  Ludde   No Calf End 1   10  -5  11.18034    2.034444    0.166667    67.082039   2.034444    Obs 1     3311_1  2010-10-07 00:20    F07011  1   1476019 6620136 before  before  Ludde   No Calf End 1   4.528243    12.549702   13.341664   0.346286    0.166667    80.049984   0.346286    rand    1     3311_1  2010-10-07 00:20    F07011  1   1476059 6620089 before  before  Ludde   No Calf End 1   45.067664   -34.189262  56.568542   2.219792    0.166667    339.411255  2.219792    rand    1     3311_1  2010-10-07 00:20    F07011  1   1476018 6620125 before  before  Ludde   No Calf End 1   4.341216    2.480695    5   1.05165 0.166667    30  1.05165 rand    1     3311_1  2010-10-07 00:20    F07011  1   1476014 6620127 before  before  Ludde   No Calf End 1   0   3.605551    3.605551    0   0.166667    21.633308   0   rand    1     3311_1  2010-10-07 00:20    F07011  1   1476020 6620126 before  before  Ludde   No Calf End 1   5.778044    3.407962    6.708204    1.037901    0.166667    40.249224   1.037901    rand    1      Here's what happens after all I do is AddData and then Open Attribute Table:               ID    ID_X   ID_Y              Date_Time    moose   moose_X          moose_Y    trial   x   y   label   label_1 dog Reproduction    Position    EvalPt  DistX   DistY   Dist    Heading TimeDif Velocity    Angle   LocType LocNum     3311_1  1   33.18333333 2010-10-08 09:20    F07011  -147.4600309    9.81849E+15 1   1476079 6620355 before  before  Ludde   No Calf End 1   -13 26  29.06888371 5.819537698 0.166666667 174.4133022 2.677945045 Obs 1     3311_1  1   33.18333333 2010-10-08 09:20    F07011  -147.4600309    9.81849E+15 1   0   0   before  before  Ludde   No Calf End 1   0   0   16.97056275 3.660738768 0.166666667 101.8233765 0.519146114 rand    1     3311_1  1   33.18333333 2010-10-08 09:20    F07011  -147.4600309    9.81849E+15 1   1476092 6620328 before  before  Ludde   No Calf End 1   0   -1  1   3.141592654 0.166666667 6   0   rand    1     3311_1  1   33.18333333 2010-10-08 09:20    F07011  -147.4600309    9.81849E+15 1   0   0   before  before  Ludde   No Calf End 1   0   0   437.7145188 5.561036001 0.166666667 2626.287113 2.419443348 rand    1     3311_1  1   33.18333333 2010-10-08 09:20    F07011  -147.4600309    9.81849E+15 1   0   0   before  before  Ludde   No Calf End 1   0   0   293.4365349 3.457878782 0.166666667 1760.619209 0.316286128 rand    1     3311_1  1   33.18333333 2010-10-08 09:20    F07011  -147.4600309    9.81849E+15 1   0   0   before  before  Ludde   No Calf End 1   0   0   75.92759709 5.98774447  0.166666667 455.5655825 2.846151816 rand    1      I do not know what R is doing because normally I have no issue, but "normally" for me Ive used write.csv and then didn't have to paste all the rest of that information together. I get that it is something with the formatting coding that goes on in the background and I am wondering if it has something to do with the write.table command and then forcing it into csv? But what is ArcGIS doing by creating ID_X ID_Y adding a moose_X and moose_Y and where are those numbers coming from? Why is it removing values (i.e., 2nd row under x and y values now are 0)? Why is it changing the date_time to 9:20? I will also post this type of question in an R forum here - but if someone has experience or tricks to deal with this from the ArcGIS side of things - that would be great. Right now my work around is to take the R output (my csv file) bring it into some "neutral" program like TextPad and then save it as a new csv or txt file, _then_ I can bring it into ArcGIS and seem to have no problems. Im also having better luck just opening the csv file in Excel and saving as xls and then opening that in ArcGIS. Here is the structure of the R input before I write to csv `str(SubData) 'data.frame': 4740 obs. of 21 variables:   $ ID : Factor w/ 29 levels "3311_1","3313_1",..: 1 1 1 1 1 1 1 1 1 1 ...   $ Date_Time : POSIXct, format: "2010-10-08 09:20:00" "2010-10-08 09:20:00" "2010-10-08 09:20:00" "2010-10-08 09:20:00" ...   $ moose : Factor w/ 15 levels "F01001","F070014a",..: 7 7 7 7 7 7 7 7 7 7 ...   $ trial : int 1 1 1 1 1 1 1 1 1 1 ...   $ x : num 1476079 1476090 1476093 1476097 1476129 ...   $ y : num 6620355 6620319 6620332 6620323 6620317 ...   $ label : Factor w/ 11 levels "after","before",..: 2 2 2 2 2 2 2 2 2 2 ...   $ label.1 : Factor w/ 5 levels "after","before",..: 2 2 2 2 2 2 2 2 2 2 ...   $ dog : Factor w/ 3 levels "Ida","Ludde",..: 2 2 2 2 2 2 2 2 2 2 ...   $ Reproduction: Factor w/ 2 levels "Calf","No Calf": 2 2 2 2 2 2 2 2 2 2 ...   $ Position : chr "End" "End" "End" "End" ...   $ EvalPt : num 1 1 1 1 1 1 1 1 1 1 ...   $ DistX : num -13 -2.35 1.19 5.16 36.63 ...   $ DistY : num 26 -9.62 2.93 -6.43 -12.29 ...   $ Dist : num 29.07 9.9 3.16 8.25 38.64 ...   $ Heading : num 5.82 3.381 0.387 2.466 1.895 ...   $ TimeDif : num 0.167 0.167 0.167 0.167 0.167 ...   $ Velocity : num 174.4 59.4 19 49.5 231.8 ...   $ Angle : num 2.68 0.24 3.53 5.61 5.04 ...   $ LocType : chr "Obs" "rand" "rand" "rand" ...   $ LocNum : num 1 1 1 1 1 1 2 2 2 2 ..."`


title: How do I draw on a vector map in GRASS?

text: I have a vector map that has two lines. I want to have the two connected. Is there a way for me to draw a line to connect the two? I am using GRASS gis by the way. The two lines that I want to connect are shown here: ![enter image description here](http://i.stack.imgur.com/pYgh7.jpg)


title: How do I calculate a weighted centroid with respect to a raster, e.g. population weighted centroid

text: In ArcGIS 10.0, how do I calculate the weighted centroid of a collection of polygons, where I want the weight to be a raster layer. In particular I am using a raster of population counts and would like to generate the population weighted centroid as a result.


title: ArcGIS ModelBuilder model works in Edit windows but not from the "open" window - do you have any experience or advice?

text: I am using ArcGIS 10.1, most recent service pack / patch, and on Windows 7 with 8GB RAM I always run Models from ArcCatalog I have a model that validates and runs from the "edit" window. When I run the exact same model with the exact same input data, but from the "open" window it returns error messages. Do you have any experience or advice to share over this issue? Are there certain tools or processes that can cause this issue? The issue that occurs most frequently is with the Focal statistics and the Get Raster Properties Tools - when the model is run from the Open window these tools can fail because the preceding data does not have valid statistics. It is unclear why the model runs OK (with exact same data) from the edit window. Surely this should be picked up by the model validation process? The errors that occur are > ERROR 001100: Failed because no statistics are available. The solution can be to add in "Calculate Statistics" tools into the models, however the behavior appears erratic and this solution has occasionally not worked, even on the same model , same data, same computer. One of the models that runs OK from the edit Window, into which I have added the Calculate Statistics tool, return the following error only when run from the Open Window. > ERROR 999999: Error executing function. The workspace is not connected. > Failed to execute (Calculate Statistics). Any ideas/ thoughts / help would be much appreciated thanks


title: Remove duplicates from Spatialite

text: I have a Spatialite database with points. From time to time now points are added. What would be the easiest way to remove duplicates based on the coordinates? Hints are greatly appreciated.


title: Create HTML report using Python

text: i want to create Analysis report after createing model for generating cut/fill. Once it will done , basedon this cutfill/table, would be display statics in html formate using python. is it possible. pl. tell me how it will done. Thanks


title: How to handle time-based raster data in the ArcGIS Server JS API

text: The TimeSlider functionality in the ArcGIS Server JS API is great for working with time-based vector features - you can scroll forwards and backwards through time to see the vector features which existed at a given time instant. **Has anyone created a similar time-slider for rasters?** I have a series of monthly weather grids which I'd like to display, allowing the user to scroll to any month (eg, show me the average rainfall for June). The time-slider relies on the featureclass being defined in ArcMap with a Time attribute, but this isn't available for raster layers. The best workaround I can think of is:   * create a separate map document for each raster   * build a separate tiled map service from each map document   * add a DOJO slider to the map, with increments to mark each raster   * use JavaScript to switch the appropriate raster on at each mark on the DOJO slider  Is anyone aware of a better/faster approach? This will be quite a laborious process as it will mean creating and maintaining multiple MXDs and MSDs. Thanks


title: Modelbuilder - If precondition is false, continue iterating next input

text: I've built a model in modelbuilder but cannot find an answer to this problem. Please help if you can. Firstly I have an 'iterate feature class' tool drawing inputs out of a folder and then selecting attributes from the input table using SQL. I am then using the "get count" tool (which counts the rows selected) to test if any attributes were selected, this is a precondition to continue the model. However, if this precondition is 'false' then the model stops. I'd like for it to just continue with the next iteration of inputs. Is this possible, and if so, any insight would be great. I was thinking I might need if-then-else logic but not sure how to use it in this setting. Please advise. Here is an image of the model. ![enter image description here](http://i.stack.imgur.com/EWzLi.jpg) Cheers, Ken!


title: How to save a WMS layer at multiple resolutions as .tif using pyqgis

text: How can I export a wms-map-extent as a .tif from QGIS? Is it possible to call the "Save as.." function of a WMS-raster-layer from pyqgis? (E.g. by simulating right-clicking the layer, choosing "Save as" and setting the save-as-options?) Im working on this for days now, but I just couldn't find a suitable solution. I managed to export the extent as an image, but this way I am losing the geo- references. I even tried to call gdal_transform from bash, but QGIS-layer- extent-coordinates (shown in the bottom right corner) do not seem to match the coordinates gdal expects (gdal turns them in some negative src values). What I want to do is:   1. Open a wms-layer in QGIS (in my case it is the city-map of Berlin from http://fbinter.stadt-berlin.de/fb/wms/senstadt/k_baugrund?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&layers=0&styles=&CRS=EPSG:3068&format=image/png&dpiMode=7)   2. Zoom/go to a specific extent of the map   3. Use a one-click plugin to export that extent (in multiple resolutions) I know I can do that by just opening the save-as dialog, but in second step I'd like to be able to open the wms and a shape layer and automatically export the masked/shaped part of the wms. This is the gdal_transform command I tried from bash, which resulted in a blank 1x1 .tif-file.                   gdal_translate -of GTiff -a_srs EPSG:3068 -projwin 17907 16971 16623 -sds WMS:http://fbinter.stadt-berlin.de/fb/wms/senstadt/k_baugrund?SERVICE=WMS&VERSION=1.3.0&REQUEST=GetMap&layers=0&styles=&CRS=EPSG:3068&format=image/png&dpiMode=7 export_gdal_berlin.tif      Any help would be highly appreciated!


title: Specification for US Census GEOID?

text: I'm writing a program that needs to generate GEOID values from 2000 and 2010 Census data, to link the census data to TIGER/Line records. I'm starting with the state, county, region block, etc, values in the census GEO files, and need to create GEOIDs. I need to work with all of the summary levels, and since the structure of the GEOID value seems to vary depending on the type of geography, it would be really nice to have a spec for constructing it. I found this document, but it only specifies the GEOID structure for the redistricting files: http://www.census.gov/rdo/pdf/0GEOID_Construction_for_Matching.pdf I'd like to get the same thing, but for all summary levels. I've considered downloading all of the tiger files and extracting the GEOID10 values and reverse engineering them from the records attached to the shape files, but the programming would cost more than just paying someone to do it manually by inspecting the files or the documentation ( http://www.census.gov/geo/www/tiger/tgrshp2010/TGRSHP10SF1.pdf ) Any other ideas how to most easily generate GEOID from census GEO files, for all records and summary levels? Alternately, can the 0GEOID_Construction_for_Matching.pdf file, which is written for the redistricting files, be used as a spec for the 2010 SF1 summary levels and the GEOIDs of the corresponding TIGER files? Thanks!


title: Export dxf files for use in Sketchup

text: the export to dxf seems like it worked, but Sketchup will not import it. I don't think this is a Sketchup issue, as Sketchup imports other dxf files just fine.


title: how to calculate a field from a shp currently sected in a arcmap session using a Python script in ArcMap 10?

text: Hi I nedd to create a script in Python for use in ArcMap 10. This script must be able to calculate a field from the shp currently selected in ArcMap. Any ideas?, I'm new to Python so any code example will be most welcome. Thank you, I'm confident you will be able to help me


title: PostGIS ST_Distance vs. pgRouting

text: I am looking to create a city directory web app, and the site will need to be able to calculate distances between two geospatial points (like the distance between the user and a restaurant, or the distance between two restaurants). I am a newbie to geospatial function, and recently started using PostGIS. I have been having a heck of a time installing pgRouting into my PostGIS 1.4/PostgreSQL 8.4 server, so am thinking about giving it up for now until I feel like I need it. I have already spent several hours trying to install pgRouting on my OSX 10.6 Snow Leopard machine with little success. With that said, what exactly am I giving up? It seems like ST_Distance will satisfy most of my needs in the short term and things like the Travelling Salesman Problem won't really apply to webapp (at least not for a while, if ever). Is there something else that I am giving up? Another reason to use pgRouting? (As a side note, if you know of a place that makes installing pgRouting easy, please point me in that direction, as I have already tried about 3 pages worth of several google results for tutorials/instructions).


title: How to use the "Iterate Rasters" tool in ModelBuilder to link "Raster to Other Format (multiple)" tool to the "Extract by attributes" tool in ArcGIS?

text: I wish to link the "Raster to Other Format (multiple)" tool to the "Extract by atttributes" tool. The "Raster to Other Format (multiple)" tool only allows you to specify an output workspace, but the "Extract by atttributes" tool requires a specific raster to be named as input. I am attempting to use the "Iterate Raster" tool to provide the specific raster needed as input to the Extract tool but cannot get the 2 tools to link. Any suggestions?


title: Use GetFeatureInfo with OpenLayers and a Bounding Polygon

text: I want to be able to draw a bounding polygon and have my MapServer WMS layer return every feature within that bounding box. I have it working for a single point, but I would like for the query to return every single feature within a polygon. Is that possible with WMS GetFeatureInfo? I am using OpenLayers in Java with GWT. There is a bridge between Java and the OpenLayers JS.


title: Clipping multipart polygons iteratively

text: I have a shapefile of buffers (n=10) and a shapefile of watersheds (n=10). Both shapefile contain overlapping features. I would like to clip each watershed feature with it's corresponding buffer, based on a common ID field (highlighted area below). I'm thinking ModelBuilder should be able to do this, but haven't been able to figure it out. I also have each watershed as an individual file if that helps. ![watershed and buffer clip](http://i.stack.imgur.com/11qKo.png)


title: Style Print dijit menu in esri JavaScript API

text: Trying to add styles to a dijit is difficult enough, as I have to track down the classes it uses after it loads. I've managed to do that, but I'm struggling to find the appropriate classes to alter for the menu that pops out from the esri Print Dijit (`esri.dijit.Print`) that lists the different layout options. This doesn't show up in the structure in the browser's development environment, so I can't track them. Is there a better way to find and alter these style classes?


title: ArcSDE Java API: Find centroid of SeShape object?

text: I've inherited an ArcSDE 9.3 Java app that needs some fixing. The basic task is to find the centroid of a polygon retrieved from an SDE geodatabase, but lately it's been returning the 0th vertex instead. This is how the code looks:               SeShape shape = sde.findShape(aPIN);  // Get the geometry     SDEPoint sdePoint = null;          try {         if (null != shape) sdePoint = shape.generateLabelPoint();  // Find "centroid"     }      ~~I'm not sure why generateLabelPoint was used -- it seems like a hack which worked for a while but has since stopped.~~ According to this post generateLabelPoint() is the correct way to find the centroid of an SeShape, but our client app is still getting a vertex. We use the same shape shape elsewhere in the code and it works fine, so I don't think it's the geometry that's bad. If I could convert it to an SePolygon there's the centroid() method, but those appear to be two different libraries without any way of translating between the two. Does anyone have any insights on this? The test shape I'm using is Los Angeles City Hall. In State Plane CA V Feet, the centroid (found via the `Polygon.centroid` method in ArcPy) is 6488062.24483722, 1841966.80439753. `generateLabelPoint` is returning 6488207.42375, 1842383.41312499, which is a vertex in the corner of the lot. The source of this data is an Oracle 11g database. I have a feeling that the architecture of the database is different than what the 9.3 API is expecting, but still looking into this...


title: Show previously loaded tiles while loading (zoom in or zoom out) map

text: How to remove the white tile part(map loading area). ![enter image description here](http://i.stack.imgur.com/bScHx.png) While zoom in and zoom out, the white tile part is appearing on the map(map loading area). How can we fill that area with the history image? Is there any better way to load maps faster than the GeoWebCache.


title: How to restrict zoom levels based on Bing Map availability?

text: **Background** My site is built on ArcGIS Server 10 with the JavaScript API 2.3, and uses Bing Maps as the background mapping. Although the site is aimed at NSW, Australia we allow the user to zoom anywhere in the world, at any of the supported scales (ie, the tiling scheme as specified in the ArcGIS Help File.) Bing Maps road/satellite datasets are available at varying scale ranges, with populated areas supporting larger scales than empty rural areas. When zooming somewhere with no map/satellite data, a "broken camera" icon is shown. We'd like to avoid showing this ugly icon. **Question** How can we restrict the maximum zoom scale where the map/satellite is not available, but allow full zoom-in where it is available? Google Maps does this quite well - as you move around the world you can see the Zoom Slider expand/contract depending on what datasets are available at that location. They presumably have metadata showing the presence/absence of data, so how can we emulate this when using Bing Maps, with no access to the metadata? Thanks


title: Display Error in Django when querying data from postgis

text:
. i am trying to query my data from postgis through Django .. and i am not sure what is going on with my code. . i just follow processes form online research and docs. . i have my results.html               {% block content %}     {% if results %}       {% for result in results %}         {{ result.id }}         {{ result.barangay }}       {% endfor %}     {% else %}
Please enter a valid UID
         <form method="get" action="/search/">           Search Barangay: <input type="text" name="q" id="id_q" value="{{ query     }}"/>           <input type="submit" value="Search" />         </form>     {% endif %}          {% endblock %}      i have also this in my views.py               from django.template import RequestContext     def search(request):              query = request.GET.get('q')         try:             query = int(query)         except ValueError:             query = None             results = None         if query:                  results = Butuan_Parcel.objects.get(id=query)              context = RequestContext(request)              return render('results.html', {"results": results}, context_instance=context)      i would like to know what am I missing why does the code gave me an error such as TemplateDoesNotExist at /search/ I have this in my urls.py                (r'search', search)      or do I have problems with my tags or anything? Can someone suggest solutions to resolve this?? than you guys for entertaining my question. . BTW, this is the error i've got. ..               TemplateDoesNotExist at /search/          {'results': <Butuan_Parcel: M.J. SANTOS>}          Request Method:     GET     Request URL:    http://localhost:8000/search/?q=3     Django Version:     1.6.1     Exception Type:     TemplateDoesNotExist     Exception Value:              {'results': <Butuan_Parcel: M.J. SANTOS>}          Exception Location:     C:\Python27\lib\site-packages\django\template   \loader.py in find_template, line 131     Python Executable:  C:\Python27\python.exe     Python Version:     2.7.0     Python Path:              ['C:\\CLAIMS_WEB',      'C:\\Python27\\lib\\site-packages\\setuptools-2.1-py2.7.egg',      'C:\\Python27\\lib\\site-packages\\pip-1.5.2-py2.7.egg',      'C:\\Python27\\lib\\site-packages\\django_leaflet-0.10.0-py2.7.egg',      'C:\\Python27\\lib\\site-packages\\django_floppyforms-0.4.7-py2.7.egg',      'C:\\Python27\\lib\\site-packages\\vectorformats-0.1-py2.7.egg',      'C:\\CLAIMS_WEB',      'C:\\Program Files\\ArcGIS\\bin',      'C:\\Windows\\system32\\python27.zip',      'C:\\Python27\\DLLs',      'C:\\Python27\\lib',      'C:\\Python27\\lib\\plat-win',      'C:\\Python27\\lib\\lib-tk',      'C:\\Python27',      'C:\\Python27\\lib\\site-packages']          Server time:    Mon, 10 Mar 2014 17:09:15 +0800


title: earthquake effects on total station surveys

text: What I have been wondering is whether an earthquake can affect a closed system survey. We are using Leica TCR 307s and TCR 407s and I was wondering if an earthquake can throw off your fixed points if staying within a 500m grid? The fault line does not run through the site.


title: Calculating distance

text: I am trying to measure straight-line distance from residence to facility. I have both of these locations in latitude and longitude (degrees) in separate layers in shapefiles. These layers indicate which facility each resident went to. I have converted the projection to be UTM (meters). I am trying to measure distance in kilometers and am not sure how to do this. I cannot do a near feature because we cannot assume that all residents traveled to the nearest facility. I calculated the the hypotenuse of a triangle through a function in excel and have the "distance" in a single number, but this is in degrees, so I am not sure if it can be converted t meters, or if the whole operation should be done differently. Thanks so much.


title: How to reproject raster from 0 360 to -180 180 with cutting 180 meridian

text: I have a geotiff raster image that has a coordinate system with longtitudes from 0 to 360. The horizontal center of the image is 180 longtitude. See image below: ![enter image description here](http://i.stack.imgur.com/uVpfl.png) I want to transform it to EPSG:4326 SRS with -180 180 longtitude range. And I want the center of the image to be at Greenwich meridian (0). I guess this srs is very widely used. I expect the result look like this: ![enter image description here](http://i.stack.imgur.com/WbCRp.png) So I use a gdalwarp command to reproject:               gdalwarp -s_srs '+proj=latlong +datum=WGS84 +pm=180dW' -t_srs EPSG:4326 test_col.tif test_4326.tif      But I only get a tiff with bigger dimensions (more pixels) and EPSG:4326 metadata. The image itself looks the same, as the initial one. But I expect it to swap the hemispheres. The question is - how do I gdalwarp an image to be strictly -180 180 EPSG:4326 with the center in 0 longtitude? This is gdalinfo of my initial file:               Origin = (-0.102272598067084,89.946211604095552)     Pixel Size = (0.204545196134167,-0.204423208191126)     Metadata:       AREA_OR_POINT=Area     Image Structure Metadata:       INTERLEAVE=BAND     Corner Coordinates:     Upper Left  (  -0.1022726,  89.9462116) (  0d 6' 8.18"W, 89d56'46.36"N)     Lower Left  (  -0.1022726, -89.9462116) (  0d 6' 8.18"W, 89d56'46.36"S)     Upper Right (     359.897,      89.946) (359d53'50.18"E, 89d56'46.36"N)     Lower Right (     359.897,     -89.946) (359d53'50.18"E, 89d56'46.36"S)     Center      ( 179.8975000,  -0.0000000) (179d53'51.00"E,  0d 0' 0.00"S)      This is gdalinfo after gdalwarp               Origin = (-180.102727401932952,89.946211604095552)     Pixel Size = (0.091397622896436,-0.091420837939082)     Metadata:       AREA_OR_POINT=Area     Image Structure Metadata:       INTERLEAVE=PIXEL     Corner Coordinates:     Upper Left  (-180.1027274,  89.9462116) (180d 6' 9.82"W, 89d56'46.36"N)     Lower Left  (-180.1027274, -89.9699975) (180d 6' 9.82"W, 89d58'11.99"S)     Upper Right ( 179.8211116,  89.9462116) (179d49'16.00"E, 89d56'46.36"N)     Lower Right ( 179.8211116, -89.9699975) (179d49'16.00"E, 89d58'11.99"S)     Center      (  -0.1408079,  -0.0118929) (  0d 8'26.91"W,  0d 0'42.81"S)


What software or template does Max Fisher use to create his maps for the Washington Post? Do you recognise the software used to make this map? As an example: ![enter image description here](http://i.stack.imgur.com/pnoGJ.jpg) More examples here: http://www.washingtonpost.com/blogs/worldviews/wp/2014/01/13/40-more-maps- that-explain-the-world/


title: Create transparent background for svg marker in QGIS?

text: I have an image that I want to use as a SVG marker in Print Composer. I created it as a .png file in Inkscape and then saved it as .svg. But the background transparency does not come across from the png to the svg. So I guess I have two questions. Can I adjust the transparency of the svg background in QGIS when I set the icon style? (and by default, how?) And, if not, does anyone know how I can preserve the transparency when I create the svg in Inkscape?


title: How to join MySQL database to QGIS polygons

text: I Have an *.sql file that has no spatial information. For every record there is an id. I also have a Qgis project with polygons without any id. How can i connect the polygons with the database? I know that i have to add the same id to the shapes to make the connection. I don't know how to make the connection. Is there any tutorial that can help me?


title: ArcPy and running Python scripts (with parameters) within another Python script

text: A common coding pattern used in AML was to run an AML (with parameters) inside another AML. An application that I'm currently developing would benefit from being able to run a Python script (with parameters) within another Python script. However, this does not seem to be at all straightforward. Using ArcGIS 10, I'm experimenting with wrapping the "inner" Python script into an ArcGIS tool that has the parameters. I thought it would be a simple matter to have the "outer" Python script use arcpy.ImportToolbox to import the toolbox and then run the tool(s) within it. However, in testing so far all my attempts to run the "inner" tool from the "outer" script appear to simply skip the "inner" tool (no error is thrown). Any thoughts on how to do this are welcome. In the meantime I'll try to create some simple test code to illustrate better what I am trying to describe. Test code is now here ... My testinner.py script is:               inputString = arcpy.GetParameterAsText(0)          newFC = "C:\\Temp\\test.gdb\\" + inputString     arcpy.Copy_management("C:\\Temp\\test.gdb\\test",newFC)      My testouter.py script is:               import arcpy          inputString1 = arcpy.GetParameterAsText(0)     inputString2 = arcpy.GetParameterAsText(1)          arcpy.ImportToolbox("C:\\Temp\\test.tbx")          arcpy.testinner_test(inputString1)          arcpy.testinner_test(inputString2)      For testinner.py its tool needs a single String parameter. For testouter.py its tool needs two String parameters The two tools are placed in a test.tbx. The test.gdb just needs a single empty feature class called test. Once you have the above assembled, running the testinner tool with a string like 'abc' passed in as its parameter should result in feature class 'test' being copied to one called 'abc' OK. But when you try running the testouter tool with two strings like 'uvw' and 'xyz' as its parameters, the testinner tool within testouter.py seems to run OK once, but sends ArcMap 10 SP2 on Vista SP2 to a Serious Application Error when trying to use it the second time. The same test using Windows XP SP3 and ArcGIS Desktop 10 SP2 also produces a Serious Application Error at the same point. ================================================================================ The execfile method described by Jason allowed me to rearrange my code to that below and provides a solution to my test problem:               import arcpy          inputString1 = arcpy.GetParameterAsText(0)     inputString2 = arcpy.GetParameterAsText(1)          arcpy.ImportToolbox("H:/Temp/test.tbx")          # Write second Python script to an ASCII file for first parameter & execute it     f = open("H:/Temp/string1.py","w")     f.write('newFC = "H:/Temp/test.gdb/' + inputString1 + '"' + "\n")     f.write('arcpy.Copy_management("H:/Temp/test.gdb/test"' + ',newFC)')     f.close()     execfile("H:/Temp/string1.py")          # Write third Python script to an ASCII file for second parameter & execute it     f = open("H:/Temp/string2.py","w")     f.write('newFC = "H:/Temp/test.gdb/' + inputString2 + '"' + "\n")     f.write('arcpy.Copy_management("H:/Temp/test.gdb/test"' + ',newFC)')     f.close()     execfile("H:/Temp/string2.py")      However, this may prove cumbersome when applied to non-test script(s) which are much longer so I used blah238's generous work that endorsed Dan's approach (Jason also endorsed Dan's approach I may add) and came up with the following final (test) code that does exactly what I need. ================================================================================               # CopyFeaturesTool.py          import CopyFeaturesUtility     import arcpy     outputFCName = arcpy.GetParameterAsText(0)     outputFCName2 = arcpy.GetParameterAsText(1)          CopyFeaturesUtility.copyFeaturesToTempGDB("C:\\Temp\\test.gdb\\test", outputFCName)     CopyFeaturesUtility.copyFeaturesToTempGDB("C:\\Temp\\test.gdb\\test", outputFCName2)      and               # CopyFeaturesUtility.py          import arcpy     import os          def copyFeaturesToTempGDB(inputFeatures, outputName):         """Copies the input features to a temporary file geodatabase.         inputFeatures: The input feature class or layer.         outputName: The name to give the output feature class."""              tempGDB = r"C:\Temp\test.gdb"         newFC = os.path.join(tempGDB, outputName)         arcpy.env.overwriteOutput = True         arcpy.Copy_management(inputFeatures, newFC)          if __name__ == '__main__':         inputFC = r"C:\Temp\test.gdb\test"         outputFCName = arcpy.GetParameterAsText(0)         copyFeaturesToTempGDB(inputFC, outputFCName)      Many thanks to Dan, Jason and blah238 for helping me.


title: pdf creator software to interface with Civil 3D

text: I've experienced a number of problems with various pdf Creator packages throughout my career, especially with different versions of AutoCAD I've used. For example, blobby, illegible text upon export to pdf, conflicts causing the process to crash, converting north arrows to a big black blob, etc, etc. The pdf print driver inherent to CAD is also sub par (in my experience). I found the free pdf Creator at THIS LINK to be the most highly recommended to interface with higher level AutoCAD licences (according to the AutoDesk forum). Feedback would be greatly appreciated regarding this software, or any other reccomendations any of you may have. Purchasing the Adobe Pro is acceptable, however:   1. Is it coded to work well with AutoCAD Civil 3D?   2. Will the free option at the link I posted serve my needs (free is obviously preferred)?   3. Is there a better option than the two I've mentioned? I'm looking for **actual user experience** (with CAD), not random free pdf creator options/suggestions. Thanks in advance **CORRECTION** I noticed at the AutoDESK forum that the answers were not listed by highest rank. There are no vote-ups whatsoever on the suggestions posted there. I assumed the answer at the top of the page was "the best" answer, while in actuality, it was just the "most recent" answer. I should have looked more closely.


title: Do vertex IDs stay the same every time I run osm2po?

text: Does anybody knows if vertex ( _source, target_ columns) in osm2po tables are always the same? I have my **dr_2po_4pgr** table with routing data and **dr_matrix** with routes already made. In this table I have _source,target,cost_ columns. This is for avoid _pgrouting_ functions that penalize the performance (I need to make many calls to these functions). My question is, if when I generate an updated **dr_2po_4pgr** (osm2po routing table) table with a new osm file should also update my table **dr_matrix** or if _source,target_ id's are always the same. Thanks in advance!! Pd: excuse my english


title: How to solve a filename conflict problem?

text: I use a procedure to remove and delete a feature layer with a certain name. Then, using another procedure which is called during the run time of the program (as the remove and delete procedure as well) I try to create a new feature layer which needs to have the same filename as that removed-deleted before. Unfortunatley, the program stops in the code line that tries to create the new feature layer. It is definetely because of the conflict between the previous and the new filename which are the same. How can I solve this problem? Any help please?


title: Resolving Arcpy "GetObject not point" Error

text: I'm new to Python (v 2.6.5) and am trying to write an ArcMap (10) script that reads the vertex coordinates of a polygon shapefile which includes a small hole (doughnut polygon). I've been working from the ESRI Help page on reading geometries. The script I've developed outputs coordinate information until it hits coordinate 1025 of 1088 at which point it fails with the error "Array: GetObject not point". I suspect this is the 'null point object' returned as the separator between my outer and inner rings (described in the Help), but I'm stumped on how to account for this in a way that my script doesn't fail (I thought the "if pnt:" line would handle that...). Here's my script - thanks for any guidance!               import arcpy     infc = arcpy.GetParameterAsText(0)     desc = arcpy.Describe(infc)     shapefieldname = desc.ShapeFieldName     rows = arcpy.SearchCursor(infc)     for row in rows:         feat = row.getValue(shapefieldname)         arcpy.AddMessage("Feature %i:" % row.getValue(desc.OIDFieldName))         partnum = 0         for part in feat:             pointnum = 0             arcpy.AddMessage("Part %i:" % partnum)             for pnt in feat.getPart(partnum):                 if pnt:                     arcpy.AddMessage("Point number " + str(pointnum) + ": " + str(pnt.X)+ ", " + str(pnt.Y))                     pointnum += 1                 else:                     arcpy.AddMessage("Interior Ring:")             partnum += 1


title: Looking for compute intensive GIS vector operations

text: I am looking to make efficient algorithms for those vector GIS operation that take a long time for GIS software like Arc-GIS, Grass etc. I have considered polygon overlay but looking for other operations as well.


title: What methods are there to produce map tiles for google maps?

text: newbie to stack exchange here. I'm currently using the automatic tile cutter found here to make the standard 256x256 tiles required to make a custom map type for my map. The map tile cutter script feels quite clumsy to use, and I quickly use up all of my machine's 4GB of ram to open photoshop and use it. That, and it's very bothersome to pad my images with empty space on the top and bottom so that it's positioned correctly. Does anyone know of other tools which can be used to make tiles for google maps?


title: Creating a multiversioned view with ArcSDE and SQL Server 2008 R2

text: I am using the administration command sdetable -o create_mv_view to create a multiversioned view on a table in SQL Server 2008 R2. Everything appears to be working correctly with one exception: The SDE_get_view_state() scalar function that was included with ArcSDE will only return data for version with name set to 'DEFAULT' and owner set to 'sde'. This is hard-coded in the function. I have some clients who are using the DBO schema for ArcSDE instead of the SDE schema, and therefore they will not have any version data that meets the criteria. Does anyone know why the function is checking for the schema owner to be set to 'sde' and more importantly, what is the best way to get this to work for the DEFAULT layer regardless of the owner?


title: expanding atributte window during typing in ArcPad

text: I want to customize my ArcPad and I need add one ability, which I cannot change: To expand (widen) atributte window during writting in ArcPad. It means, when I'm typing an atributte, when running out of space - windows make wider. Is this function abble or any script to do it?


title: LLA vs ECEF--which to choose?

text: I will be providing points in 3D to an external source via XML. One coworker suggests Lat/Long/Alt and another suggests ECEF. We can do either without much trouble. We have an existing product that uses ECEF, and would like to bolt more things onto it. The new features are natively in Lat/Long/Alt. Should I prefer LLA or ECEF for communicating with the external world? As background, this is an embedded system outputing live updates of items moving in 3D space, at a rate of around 1Hz. The items tracked will be within a 300km radius around our customer's facility. Our customer is writing custom software to parse this on their end and use the data in their existing tool. We can choose any format we wish. I just wondered if there was a preferable format.


title: Arcpy Pythonaddins Equivalent of arcpy.GetParameterAsText

text: I've build a neat little search arcpy script. The user opens any mxd, zooms to a location, chooses a corporate fGDB, and runs the tool. The user gets a list of all the populated feature class that intersects the current extent. Later I aim to write a process that the user can choose which layers they want to add. However I digress. I'm looking at embedding the arcpy into an addin, my first issue is learning what the equivalent of a arcpy.GetParameterAsText is, as I want the user to specify the work-space they need to interrogate.               import arcpy     import pythonaddins     class ButtonClass1(object):     """Implementation for AddIn_Test4_addin.button (Button)"""     def __init__(self):         self.enabled = True         self.checked = False     def onClick(self):         ################## MAP DOCUMENT DATA FRAME EXTENT              StagingGDB = arcpy.GetParameterAsText(0)         #StagingGDB = "C:\\Data\\Projects\\SearchTool\\CorpData1.gdb"         arcpy.env.workspace = StagingGDB         ExtentFC = StagingGDB +"\\ExtentFC"              ############ Clean Up & Pre Check         if arcpy.Exists(ExtentFC):             arcpy.Delete_management(ExtentFC)              ################## MAP DOCUMENT DATA FRAME EXTENT         mxd = arcpy.mapping.MapDocument("current")         df = arcpy.mapping.ListDataFrames(mxd)[0]         currentExtent = df.extent              ################## EXTENT FC         pnt1 = arcpy.Point(currentExtent.XMin , currentExtent.YMin)         pnt2 = arcpy.Point(currentExtent.XMax , currentExtent.YMin)         pnt3 = arcpy.Point(currentExtent.XMax , currentExtent.YMax)         pnt4 = arcpy.Point(currentExtent.XMin , currentExtent.YMax)         array = arcpy.Array()         array.add(pnt1)         array.add(pnt2)         array.add(pnt3)         array.add(pnt4)         array.add(pnt1)         DisplayExtent = arcpy.Polygon(array)         arcpy.FeatureClassToFeatureClass_conversion(DisplayExtent , StagingGDB ,"ExtentFC", "", "", "")         pass      Any ideas?


title: How to Select by Location and assign attributes using Python?

text: I have a polygon layer with a 'lot' field, and a points layer. How can I use python to select the points by location (to find what points are contained in what polygons) and assign the value in the 'lot' field from the polygon layer to the point that is contained in the relevant polygon feature?


title: How does Leaflet layer ordering work with CartoDB layers?

text: All, I'm trying to figure out how Leaflet layer ordering works when you have layers coming from CartoDB. I have a map with 1 polygon layer (cartodb), 3 point layers (cartodb), 1 basemap layer (mapbox) (and in the future hope to add tabletopJS or another google spreadsheet conduit...) The Leaflet Layers example seems straight forward when using cloudmade layer URL's that end in the {z}/{x}/{y}.png format, which will work for Mapbox as well. However, I can't see how to get the {z}/{x}/{y}.png format for a CartoDB layer, and wonder if the same interactivity will be supported (custom infowindows, etc.) So far I have tried this to no avail:               function layersMap() {         var voices = cartodb.createVis('map', 'http://opendetroit.cartodb.com/api/v2/viz/08140c18-8a05-11e3-ae46-0e49973114de/viz.json');             devPoints = cartodb.createLayer('map', 'http://opendetroit.cartodb.com/api/v2/viz/62a64dd6-4d7d-11e3-b48c-db82dabd16be/viz.json');          var points = L.layerGroup([voices, devPoints]);         var map = new L.map('map', {         center : [42.339926, -83.04137],         zoom : 13,         layers: [points]         });          }      The problem with this attempt, however, is the 2nd layer doesn't draw, and there is actually two zoom controls displaying. This is as close a solution as I can find by going through as much info on GIS.SE, Leaflet docs, CartoDB docs, etc. So it seems I might be on the right track by getting 1 CartoDB layer to display...but not sure where to go from here...


title: How can I get the profile graph on ArcMap to display specific points?

text: I currently have 2 layers: 1) DEM of existing terrain 2) Polygon in vector format of future water retention basin I need to draw cross sections across this basin for eventual cut and fill calcs. I already know how to create a profile of the existing terrain based on my DEM. However I don't know how I can have the profile curve display the extents (or at least start and end points) of the basin on the graph. IS there any way to incorporate vector data into a profile drawn of a DEM? Or to have the profile curve display the specific points that are the start and end points of the basin?


title: Invalid geometry message when using Difference tool

text: I am using the Vector \ Geoprocessing \ Difference tool to get the resulting area when taking a) current ground disturbance away from b) historic total ground disturbance to give c) what has been cleaned up over the years. I have a maximum of 41,000 nodes in the larger file (b). I've run the geometry validity checker and noted 4 duplicates. All else seems ok. However when I run the tool I get an invalid geometry error message and only get about 70% of the polygons I expect to get. The OGR log says for these says **"Pointer 'hLayer' or 'hFeat', is NULL in 'OGR_L_CreateFeature'**. I've done this in MapInfo with these polys with no issues (the files have been brought as .TAB's). Are there any other tools / checkers to use? I would assume the Check Geometry Validity would pick up problems or is this limited? I note a **bug #9152** was logged last year for a similiar issue. How do I find out where this is at / status and what I can do about it? Thanks.


title: How to show world map with Lizmap

text: I would like to show world map with Lizmap web client but I don't know how to do. Could you show me how? I follow this tutorial: http://hub.qgis.org/wiki/lizmapwebclient/Installation_en for install lizmap web client. But I just load a part of world. ![enter image description here](http://i.stack.imgur.com/3DVnK.jpg)


title: r.topidx error: required parameter output not set

text: I'm calling GRASS r.topidx from QGIS 2.2.0. I can get other tools to work on the DTM, for example r.watershed. But r.topidx I get the Oops message: Oooops! The following output layers could not be open output: C:/Users/User/Documents/GIS DataBase/Gills/Terrain/TopIdx.tif The above files could not be opened, which probably indicates that they were not correctly produced by the executed algorithm Checking the log information might help you see why those layers were not created as expected This algorithm requires GRASS to be run. A test to check if GRASS is correctly installed and configured in your system has been performed, with the following result: GRASS seems to be correctly installed and configured Looking at the log, I find: C:\PROGRA~1\QGISVA~1\bin>r.topidx input=tmp14037940151265 Output topographic index layer=Output topographic index layer1a269dc6704c40399e5421e3bc50f16e --overwrite Sorry is not a valid option Sorry is not a valid option Sorry is not a valid option Sorry, is not a valid parameter Sorry is not a valid option Sorry is not a valid option Sorry is not a valid option ERROR: Required parameter not set: (Output topographic index map) But, as you can see from the Ooops (very annoying) message, the output file is set to /TopIdx.tif. which should be OK? Any idea what might be going wrong here? Thank you, John


title: Updating Shape File Datastores on Geoserver

text: I am a web developer working with GIS guys who are concurrently working on a vector layer that updates every week or so, however, on our mapping application there are about 30 layers that would need to be deleted and recreated with the Geoserver interface. This process kills about 5 hours every week by hand and I want to automate the process. I am handed a Shapefile and a list of layers which do not change. I imagine I am essentially updating the datastore, but when I try to update just the datastore the layers in my application stop working. How can I update just the datastore without deleting and re-adding every layer?


title: Why is Oracle Spatial Query Layer Not Displaying in ArcMap?

text: I am having an issue with an Oracle query layer using SDO_Geometry not displaying. I have created a polygon layer/table within oracle, registered the table with the SDO_GEOM_METADATA and created a spatial index...twice. Unfortunately, I am still unable to see the polygons display in an ArcMap query layer. I am able to bring the data into ArcMap, open the table, select and identify it and actually see the shape flash, so I know it is there. Also, I can export a shapefile and the shapes render perfectly fine, but for some reason they will not display from the query layer directly. I have checked over the index parameters in the registry and nothing seems off, see:               INSERT INTO all_sdo_geom_metadata                 (owner,                  table_name,                  column_name,                  diminfo,                  srid)          VALUES ('DB',                  'TABLE',                  'GEOM',                  MDSYS.sdo_dim_array (MDSYS.sdo_dim_element ('X',                 -180,180,5.0e-10), MDSYS.sdo_dim_element ('Y',-90,90,5.0e-10)),4326);     index creation:     CREATE  INDEX wm_dim_table_si ON TABLE (geom) INDEXTYPE IS MDSYS.spatial_index      PARAMETERS ('tablespace=DB_INDEX03  layer_gtype=MULTIPOLYGON');      I've done everything, others have prescribed but I still can't view the data directly... Any hints as to what else I should try or should be looking at? Any help is greatly appreciated. My nexts thoughts are that I should check and validate the geometry on the polygons. Also, there are about 30,000 polygons in this table, is there a threshold in ArcMap for query layers? It wouldn't seem so, however I am at a loss right now...


title: Tile script (Mapnik) with OpenLayers - Problem with projection (900913)

text: I use Mapnik python bindings in a geodjango application to render tiles in openlayers with datasource stored in postGIS in 900913 projection. The feature seems to be rendered correctly in the spherical mercator projection, but at a very high zoomLevel in the very left corner of the world map so the features seems really tiny. _EDIT: After several test the problem seems to be in the first lines of the script where minLong, minLat, maxLong, and maxLat are defined. Someone has a similar tiling script in 900913?_ ![enter image description here](http://i.stack.imgur.com/lq2yC.png) This is weird because the bounds of the features are: > -8879435.28405854, 5620161.13498339, -6337967.06663213, 8998594.89464899 but they are rendered around this position in the map: > lat: -19955850.825909 lon: -19952809.28254 There is the script               MAX_ZOOM_LEVEL = 20     TILE_WIDTH = 256     TILE_HEIGHT = 256          dbSettings = settings.DATABASES['default']          def tile(request, version, shapefile_id, zoom, x, y):              try:             if version != "1.0":                 raise Http404             try:                 shapefile = Shapefile.objects.get(id=shapefile_id, created_by=request.user)             except Shapefile.DoesNotExist:                 raise Http404             zoom = int(zoom)             x = int(x)             y = int(y)             if zoom < 0 or zoom > MAX_ZOOM_LEVEL:                 raise Http404             xExtent = _unitsPerPixel(zoom) * TILE_WIDTH             yExtent = _unitsPerPixel(zoom) * TILE_HEIGHT             minLong = x * xExtent - 200037508.34             minLat = y * yExtent - 200037508.34             maxLong = minLong + xExtent             maxLat = minLat + yExtent                  #create de mapnik.map object             map = mapnik.Map(TILE_WIDTH, TILE_HEIGHT, "+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +units=m +k=1.0 +nadgrids=@null +no_defs")             map.background = mapnik.Color("#7391ad")                  #defining the feature layer             geometryField = utils.calcGeometryField(shapefile.geom_type)             query = '(select ' + geometryField + ', id from "shapefile_feature" where' + ' shapefile_id=' + str(shapefile.id) + ') as geom'                  datasource = mapnik.PostGIS(user=dbSettings['USER'],                             password=dbSettings['PASSWORD'],                             dbname=dbSettings['NAME'],                             table=query,                             srid=900913,                             geometry_field=geometryField,                             geometry_table='"shapefile_feature"',                             simplify_geometries=True)                  featureLayer = mapnik.Layer("featureLayer")             featureLayer.srs = "+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0.0 +k=1.0 +units=m +nadgrids=@null +no_defs +over"             featureLayer.datasource = datasource             featureLayer.styles.append("featureLayerStyle")                  #defining the feature layer styles             rule = mapnik.Rule()             if shapefile.geom_type in ["Point", "MultiPoint"]:                 rule.symbols.append(mapnik.PointSymbolizer())             elif shapefile.geom_type in ["LineString", "MultiLineString"]:                 rule.symbols.append(mapnik.LineSymbolizer(mapnik.Color("#000000"), 0.5))             elif shapefile.geom_type in ["Polygon", "MultiPolygon"]:                 rule.symbols.append(mapnik.PolygonSymbolizer(mapnik.Color("#f7edee")))                 rule.symbols.append(mapnik.LineSymbolizer(mapnik.Color("#000000"), 0.5))             style = mapnik.Style()             style.rules.append(rule)                  #add new feature to the map             map.append_style("featureLayerStyle", style)             map.layers.append(featureLayer)             map.srs = "+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +units=m +k=1.0 +nadgrids=@null +no_defs"             #rendering the map tile             map.zoom_to_box(mapnik.Box2d(minLong, minLat, maxLong, maxLat))             image = mapnik.Image(TILE_WIDTH, TILE_HEIGHT)             mapnik.render(map, image)             imageData = image.tostring('png')             return HttpResponse(imageData, image)              except:             traceback.print_exc()             return HttpResponse("")          def _unitsPerPixel(zoomLevel):         return 400075016.68 / math.pow(2, zoomLevel)      In openlayers.html               map = new OpenLayers.Map('map',               {maxResolution: 156543.0399,                 projection: new OpenLayers.Projection("EPSG:900913"),                 units: 'm',                 maxExtent: new OpenLayers.Bounds(-20037508.34,-20037508.34,20037508.34,20037508.34),                 numZoomLevels: 21});     tiledLayer = new OpenLayers.Layer.TMS('TMS',                     "{{ tmsURL }}",                     {serviceVersion: "1.0",                     layername: "{{ shapefile.id }}",                     type: 'png',                     projection: map.projection});     map.addLayer(tiledLayer);      What am I missging? why the features are not rendered at their normal sizes and positions?


title: File size limitation for shape file

text: I am going through File size limitation for shape file esri support article (this article applicable for ArcGIS 9.X version not for 10.X ) I know may be same article will applicable for 10.x version So any update for 10.1 version?? **Purpose:** Above question asked from client, side and I explained and given support article just wanted to know updated information for Windows OS (win 7/8) & Linux.Looking in white paper also (10.1) Thanks in advance : )


title: Creating a systematic rectangular grid

text: How can I produce a point shapefile (or simply point coordinates) that represent a systematic rectangular grid? At least three input variables would be required: distance between points, distance between rows of points and orientation (azimuth of rows). The purpose is creation of sample locations for natural resource surveys. I guess a fourth variable would be required as well: either a starting location or some definition of the area to be encompassed by the grid. Typically an area of interest can be defined by a polygon.


title: Intersection between two LineString gives "Non-noded intersection"

text: I am using NTS in my application and I have two line strings which are as follows: Line String 1: > LINESTRING (51055.227410858584 52102.041513978584, 51073.39084032347 > 52096.46568847251, 51183.32738708461 52062.71727093575, 51280.83632421188 > 52032.78389190314, 51415.62809024075 51991.40539735807, 51575.27507553736 > 51942.39682580468, 51669.91610274913 51913.34384027303, 51745.43773052417 > 51890.16014474777, 51765.32190593836 51884.05608314112, 51855.59683603979 > 51952.36411358454, 51916.20234851714 51998.2222846924, 52004.718294372215 > 52065.199350389405, 52004.718294372215 52065.199350389405, 52033.42616870359 > 52086.921641966816, 52046.98266491563 52097.17939076726, 52057.349397313075 > 52105.02355161466, 52075.69053924701 52118.901682344665, 52119.54979169772 > 52152.08851669904, 52122.739555512315 52154.502104652085, 52133.10628790976 > 52162.34626549949, 52137.89093363166 52165.96664742906, 52137.89093363166 > 52165.96664742906, 52143.4730203072 52170.19042634689, 52227.20432044038 > 52233.54711011433, 52300.56888817612 52289.05963303437, 52392.2745978458 > 52358.45028668442, 52402.64133024324 52366.294447531815, 52415.40038550163 > 52375.948799344, 52428.95688171367 52386.20654814444, 52645.86082110629 > 52550.33052895152, 52739.95885363691 52621.53137356636, 52805.34901183616 > 52671.00992660379) Line String 2: > LINESTRING (51313.564671633 52022.73691510973, 51331.728101097884 > 52017.16108960366, 51441.66464785903 51983.412672066894, 51539.173584986296 > 51953.47929303429, 51673.96535101517 51912.10079848921, 51833.61233631178 > 51863.09222693583, 51928.253363523545 51834.03924140418, 52003.77499129859 > 51810.85554587893, 52131.874967524614 51771.53130283609, 52204.52868538415 > 51749.228000811796, 52310.64135225795 51716.653441276314, 52345.05627124405 > 51706.08871926481, 52361.30776076527 51701.09982275937, 52373.73537039913 > 51697.284784255215, 52391.03842688936 51691.973076799426, 52394.74473141134 > 51690.83531218692, 52397.82849878731 51693.16869616807, 52441.687751238016 > 51726.35553052244, 52444.87751505261 51728.76911847549, 52455.24424745006 > 51736.61327932289, 52460.028893171955 51740.23366125246, 52465.6109798475 > 51744.45744017029, 52521.20645015943 51786.52467937298, 52522.57757869049 > 51819.852573264536, 52526.35931204718 51911.77481482129, 52531.08647874304 > 52026.677616767236, 52531.62085410866 52039.666629161125, 52532.27854686634 > 52055.65310595361, 52532.97734542138 52072.63873754562, 52544.158122302026 > 52344.40884301777, 52549.00860638995 52462.30910936231, 52552.37928177309 > 52544.23980292377) When I try to calculate the intersection between these two line strings NTS gives me an exception that there is a non-noded intersection between them and the exception is as follows: > found non-noded intersection between LINESTRING(51669.9161027491 > 51913.343840273, 51745.4377305242 51890.1601447478) and > LINESTRING(51673.9653510152 51912.1007984892, 51748.1692113208 > 51889.3216323774) [ (51728.806436824, 51895.2656338731, NaN) ] I have tried to replicate the same with JTS TestBuilder but could not do it as it was reporting no errors. Following are the snaps of the result from JTS TestBuilder: ![enter image description here](http://i.stack.imgur.com/WbxnN.jpg) So, in the above snap you can see that JTS is giving the correct intersection points whereas NTS is not. Why it is so? Please help me out. Thanks!


title: How can I create a line with three points with Python in QGIS?

text: I'm new to Python and I'm having some dificulties. I want to create a simple layer with a line in the console of QGIS with Python. How can I do that?


title: How do I change to a Historical version using Python and arcpy

text: I am having an issue trying to change to a historical version on an SDE feature class using Python and arcpy.ChangeVersion_management. I am able to manually change the version to a specific date and time using ArcMap. I am also able to fully automate the process using ModelBuilder. When I export the ModelBuilder model to a python script the script runs without error but the version that is selected is the historical version with todays date and not the date selected as the date parameter. Using Arcmap and the Python window (Geoprocessing -> Python) I can also run the python code without error and the result is the same. The historical layer is selected with the current date and not the date selected as the date parameter. Here is the Python syntax I am executing:               arcpy.ChangeVersion_management("Parcels", "HISTORICAL", "", "7/1/2013 4:30:00 PM")      I also get the same results using:               historyDate = datetime.datetime(2011, 7, 1)     arcpy.ChangeVersion_management("Parcels", "HISTORICAL", "", historyDate)      while the following throws a RuntimeError:               historyDate = datetime.date(2011, 7, 1)     arcpy.ChangeVersion_management("Parcel", "HISTORICAL", "", historyDate)      Is this a known bug or am I just going about this the wrong way? Any suggestions are appreciated.


title: General GIS Questions for Job interviews

text: I am trying to gather general GIS questions that are asked in an interview process, not related to any particular software but questions that can test one's core knowledge about GIS. I will appreciate if you can post any resources/links that are helpful that lists the GIS general questions. Thank you in advance.


title: Convert SqlGeometry data from SQL directly to Tiles / PNG file

text: I would like to know a basic method for 'exporting' SqlGeometry data from SQL Server 2008 to a Tile Server / PNG File generator. I already have a Silverlight Map application which uses a Custom Tile layer to view Tiles. I would like to add another (this generated tiles) transparent layer of PNG files on top of my base layer. The 'Spatial Results' in SQL Server Management Studio 'picture' needs to be converted to 256x256 tiles made up of PNG files, with a quadkey so it knows where to place it on the map. Is there a basic 'quick and easy' way to do this? And, as data gets updated in SQL, the tiles needs to be re-generated. Thanks a lot. Cheerio


title: Mixing geometry types in one PostGIS table

text: I've faced with the following problem. I have to migrate from Oracle database to PostgreSQL+PostGIS. Currently all geometries of all types are stored in one table and each record contains "lid" field which indicates features of the same layer. What are pros and cons of using such way? Should I break up data into multiple tables if I don't need to use database with third-party software? What about performance of spatial query, whether the indexes to help me? Big thanks.


title: Raster topo map "collars" overlap

text: I am using QGIS 2.0. I have downloaded USGS 1:100,000 topographic maps which show adjacent areas. The topographic maps were digitized with the borders (I've seen them called "collars", I think) as part of the rasters. When I load the vertically (north-south) adjacent rasters, the border areas of course overlap the map areas of each other. I have tried to convert one raster to a vector so that I could use the vectorized map as a mask to crop the other raster, but it didn't work. (My computer just kept spinning -- sort of froze up.) I have not found, here or elsewhere, discussions of using one raster to crop or clip another. So, can I use one raster map to crop another, adjacent- but-overlapping raster map? If so, how can I do it? Thank you.


title: How could I insert a 3-D LineString into postgis table?

text: Using:               CREATE TABLE trajectories(     gid SERIAL PRIMARY KEY,     vehID VARCHAR(10),     state INTEGER,     geom geometry(MULTILINESTRING)     )      I get an error, for instance, after inserting:               INSERT INTO trajectories(vehid,state,geom)     VALUES(001140,2,ST_GeomFromText('MULTILINESTRING((1 2 3,4 5 6))'))      it gives: > Geometry has Z dimension but column does not. How could I circumvent this error? thanks!


title: QGIS Export Color mapped GeoTIFF to GeoTiFF

text: I've loaded a GeoTiff in QGIS (Single Band), and applied a Color Mapping. I can view it in the map canvas, but I would like to export the _entire_ mapped GeoTiff as an appropriate 3 Band GeoTiff (i.e to R,G,B values 0-255). I know that you can export the current view -- but this is too low-resolution, and I would like to retain the current resolution. Is there a hook in the Python QGIS module that might enable?


title: GDAL WMS Question

text: I'm new to GDAL and am currently trying to hit a WMS using gdal. I need a way to save the image to my desktop. I currently have an xml file with the WMS description syntax that gdal wants -> http://www.gdal.org/frmt_wms.html. How would I save the returned image to my desktop? Thanks               using System;     using System.Collections.Generic;     using System.ComponentModel;     using System.Data;     using System.Drawing;     using System.Linq;     using System.Text;     using System.Windows.Forms;     using OSGeo.GDAL;     using OSGeo.OGR;     using OSGeo.OSR;     using System.IO;     using System.Drawing.Imaging;     using System.Xml;          namespace GDAL_TEST_WMS     {     public partial class Form1 : Form     {         public Form1()         {             InitializeComponent();         }              private void button1_Click(object sender, EventArgs e)         {             Ogr.RegisterAll();             Gdal.AllRegister();                  XmlDocument doc = new XmlDocument();                  Dataset ds_in = Gdal.Open("<GDAL_WMS><Service name=\"WMS\"><Version>1.3.0</Version><Layers></Layers><CRS>CRS:84</CRS><ServerUrl>http://geoint.lmic.state.mn.us/cgi-bin/mncomp?</ServerUrl></Service><Path>./gdalwmscache</Path><DataWindow><Projection>EPSG:4326</Projection><SizeX>600</SizeX><SizeY>600</SizeY></DataWindow></GDAL_WMS>", Access.GA_ReadOnly);            }      }     }      I'm currently using GDAL version 1.10.1. This is all the code that I have. It's currently filling the Dataset variable because I can see the properties when I debug. If there's a better way to do this please feel free to tell me. Like I said before I am REALLY new to gdal, but I really want to get this done. Thanks in advance! EDIT: Hey I've changed my code a little bit. Now I have my wms request inside of a separate xml file.The request that I'm replicating inside of my xml file is -> http://geoint.lmic.state.mn.us/cgi- bin/mncomp?request=GetMap&service=WMS&bbox=-102,30,-85,51&version=1.1.1&layers=mncomp&width=1000&height=1000&SRS=EPSG:4326. I load that xml file and use gdal.open to access the image that should be returned. using System; using System.Collections.Generic; using System.ComponentModel; using System.Data; using System.Drawing; using System.Linq; using System.Text; using System.Windows.Forms; using OSGeo.GDAL; using OSGeo.OGR; using OSGeo.OSR; using System.IO; using System.Drawing.Imaging; using System.Xml; using GDALRD;               namespace GDAL_TEST_WMS     {                     public partial class Form1 : Form     {              public Form1()         {             InitializeComponent();         }              private void button1_Click(object sender, EventArgs e)         {             Ogr.RegisterAll();             Gdal.AllRegister();                  XmlDocument doc = new XmlDocument();             doc.Load(@"./WMS_Description_File.xml");                  Dataset ds = Gdal.Open(doc.InnerXml, Access.GA_ReadOnly);             string outfile = "./test.jpeg";             OSGeo.GDAL.Driver drive = Gdal.GetDriverByName("JPEG");             Dataset dso = drive.CreateCopy(outfile, ds, 0, null, null, null);         }     }      } However i keep getting this error: ![enter image description here](http://i.stack.imgur.com/wCVna.png) Any suggestions to fix this?


title: How to create a route and add that route as layer to map using ArcPy?

text: I wrote a little script for ArcGis, and when I make it run by GUI it happens that it seems to create routes but the routes do not appear in the TOC. Consequently, when the script tries to execute operation on those routes, the dataset cannot be found. How can I solve this issue? Thanks in advance, Annalisa


title: Best way to create tiles

text: Within a GeoServer framework, does anybody know how I might best create/cache tiles based on layers stored in PostGIS? I am aware of Mapnik and have created a simple TMS using python (+ Django) but this seems 'clunky' and I'd like to know if there are any better ways to go about it. Thanks for any help / advice or experiences shared. John


title: Google Maps basemap not showing up when KML file added as a vector layer in QGIS - only point markers appear

text: I am trying to make a basic XYZ contour plot in QGIS 1.8.0 (longitude, latitude, depth in an excel file) and can't get past the first step of adding a basemap in QGIS. I've tried using the Open Layers plugin and adding a google satellite image, but I can't figure out how to georeference it and define the coordinates. Then I tried exporting a google map of the area as a KML file and adding it as a vector layer, but only the 4 point markers show up. ![This is what the screen looks like](http://i.stack.imgur.com/UbhZa.png) Does anyone know what i'm doing wrong? Thanks so much, leslie


title: How to build GeoExplorer on Windows

text: I am trying to utilize GeoExplorer (made by OpenGeo) to browse data provided by GeoServer. So far, I have downloaded source code from https://github.com/opengeo/GeoExplorer. Some instructions I found mention Ant and Git software to build app., but I can't make it work. Could somebody suggest me how to set up GeoExplorer (non-Suite) on Windows? If you could ouline the required steps, that would be great. Thanks,


title: Clip raster with shapefile

text: I want to clipping raster (Lidar) with shapefile in qgis ver.2.0.1. Which plugin in qgis could I use it ? Thanks for help me.


title: How to register shapefiles to an image?

text: I'm generating images from various sources and the projections are not necessarily known ahead of time. The images are subject to interpolation and other alterations while being rendered, but I can generate 2D arrays of latitude and longitude points that correspond to pixels in the rendered image. I need to be able to draw geo-political boundaries on top of the images using the latitude/longitude arrays as a reference. I don't normally work with GIS type problems, but I've done some research and I think what I want to do is take in a shape file with lat/lon points defining line segments making up geo- political boundaries and interpolate the closest locations for the edges of those line segments on my image. Although I have an idea for a general approach, is there a preferred way to do this in Java, or a library providing the functionality? I found GeoTools which would do what I want if the projection was known ahead of time, but I need to be able to generate overlays on the fly. In addition I need the following features:   * Must work with varying resolution, boundaries should still be legible when zoomed in heavily. That is, when a pixel in the source image is being represented by multiple pixels at the display, boundaries should be drawn gracefully.   * Must be able to render maps in real time, or at least with reasonable latency. Although the projections are not all necessarily known ahead of time, they will not change often. If some amount of automated preprocessing can be done ahead of time to improve performance, that is an option.   * An open source (but commercial friendly) solution is preferred, but I can go with a closed source solution if necessary.


title: How do I only view selected rasters from a geodatabase raster catalog in ArcGIS 10.0

text: I am using a geodatabase raster catalog to hold raster files for multiple countries, mostly for the sake of organization. I don't actually have the need (often) to visualize or draw these all at once. When working on some country specific task, I'd like the current map document to only show the related raster. Also, for spatial analysis, I would want the selected raster to be the only raster operated on. How do I accomplish this?


title: ArcGIS 10/arcpy: Help needed with syntax for CalculateField_management in Python script

text: I need to use CalculateField_management in a Python script, and the calc expression needs to concatenate a runtime variable, a string, and a field (string) value. I have searched all over for examples from which to learn. I have tried many things I have found in my research and some of my own ideas, too, but I just can't get the calc expression to evaluate. I am now completely confused and need some help to dig out of this hole. So, here's the scoop. The field to be calculated, DESCRIP, is a string field in a fGDB feature class with one feature. The runtime variable, _enddesg_ , will be a string. Essentially, what I want to do is concatenate the runtime value with a (longish) string and prepend that to the existing value in the DESCRIP field, storing the result back in the DESCRIP field. I have tried many forms of an "all-in-one" string for the expression, but nothing worked. I came to the conclusion that a pre-logic approach would be better, but I haven't been able to get that to work either. Below is what I have right now.               expr = "getDescription(str(enddesg))"     codeblock = """def getDescription(endstr):         exprstr = endstr + ' Touchdown Zone Elevation;duplicated from existing control point;' + !DESCRIP!         return exprstr"""     arcpy.CalculateField_management(perpptsSX_src, "DESCRIP", expr, 'PYTHON_9.3', codeblock)  #Should be only one feature in the source      (BTW, the _enddesg_ variable is either a string or an integer but is converted to a string as input.) I get the following error:               <class 'arcgisscripting.ExecuteError'>: Failed to execute. Parameters are not valid.     ERROR 000989: Python syntax error: Parsing error <type 'exceptions.SyntaxError'>: invalid syntax (line 2)          Failed to execute (CalculateField).      The following code variation produces no errors, but I also get nothing entered in the DESCRIP field.               expr = "getDescription(str(enddesg), !DESCRIP!)"     codeblock = """def getDescription(endstr, desc):         exprstr = endstr + ' Touchdown Zone Elevation;duplicated from existing control point;' + desc         return exprstr"""      So, I am now thoroughly confused. I work at home and do not have access to another fresh pair of experienced eyes to look at my code. If anyone can help me out, I would much appreciate it. -- Ti


title: Google Places Autocomplete Limit Suggestions to Local Area

text: I'm not sure if this is the right place to ask, but can I limit the suggestions that Google Places Autocomplete find to say 5 km of a certain location? If I'm in Dublin I don't want to see suggestions in frickin Canada, amirite! https://google- developers.appspot.com/maps/documentation/javascript/examples/places- autocomplete


title: ogrinfo sql example

text: I'm having trouble learning how to query my shapefiles for certain features. I have a shapefile that looks something like this when I summarize it using the following command:               ogrinfo -al -geom=OFF tl_2013_us_uac10.shp          OGRFeature(tl_2013_us_uac10):3599       UACE10 (String) = 31450       GEOID10 (String) = 31450       NAME10 (String) = Franklin--Highland Lake, NJ       NAMELSAD10 (String) = Franklin--Highland Lake, NJ Urban Cluster       LSAD10 (String) = 76       MTFCC10 (String) = G3500       UATYP10 (String) = C       FUNCSTAT10 (String) = S       ALAND10 (Real) = 84345960       AWATER10 (Real) = 3320730       INTPTLAT10 (String) = +41.1696387       INTPTLON10 (String) = -074.5456645          OGRFeature(tl_2013_us_uac10):3600       UACE10 (String) = 16171       GEOID10 (String) = 16171       NAME10 (String) = Chester, NY       NAMELSAD10 (String) = Chester, NY Urban Cluster       LSAD10 (String) = 76       MTFCC10 (String) = G3500       UATYP10 (String) = C       FUNCSTAT10 (String) = S       ALAND10 (Real) = 9360429       AWATER10 (Real) = 3420       INTPTLAT10 (String) = +41.3627073       INTPTLON10 (String) = -074.2780785      How can I query for features that have NJ in their Name10 field, for example? I tried something like this:               ogrinfo -al -sql "SELECT * FROM 'tl_2013_us_uac10' WHERE UACE10 = 16171" tl_2013_us_uac10.shp      but I get this error:               INFO: Open of `tl_2013_us_uac10.shp'           using driver `ESRI Shapefile' successful.     ERROR 1: Type mismatch or improper type of arguments to = operator.      Can anyone help me out or point me to a good resource?


title: Trigger a click event on a leaflet map?

text: Is it possible to make a click event occur at a particular lat/lon on a leaflet map? I would like to make a popup appear that is associated with a point at a particular location. I have tried several combinations of `map.fireEvent('click', latLon)` and `map.fire('click', latLon)` but I can't find an examples or references describing this approach. Is this some thing that can be done? One approach that I have tried (and which works) is to identify the object in the leaflet layer and to trigger a click event using `layer[id].fireEvent('click')`. However, this requires that I know the object ID and I would like to be able to just use a lat/lon to trigger this click. (Updating this based on the discussion below). I'd like to open a popup at the clicked location and I don't want to add a marker or add a new popup at this location. What I am trying to do is to fire a click event at a particular lat/lon where a point already exists and to trigger the popup that already exists.


title: How do I convert OS coordinates (X and Y) to longitude and latitude coordinates?

text: How do I convert OS coordinates (X and Y) - Eastings and Northings to longitude and latitude coordinates? For example X and Y below X (Eastings): 347904 Y (Northings): 287484


title: Use existing property pages with new IFeatureRenderer

text: I am trying to override the drawing of a feature within a custom plugin data source. I believe this entails making a custom `IFeatureRenderer`. However, the "renderer" I'd like to design would be 99% the same as a ClassBreaksRenderer, but with a custom Draw override. What would be the simplest manner to be able to reuse the existing functionality of the ClassBreaksRenderer (including and especially the UI property pages) while still being able to provide a custom `Draw` implementation?


title: Create a esri layer file through web map

text: Is there some way to create a layer file through web mapping application. If I expose data through a map service is there a way to create a .lyr file?


title: What causes "RuntimeError: DataFrameObject: Error in setting extent"?

text: I'm having a problem whilst trying to set extent on a python script. mapDoc is created using Web_map_as_Json string. Appreciate if anyone can help.               result = arcpy.mapping.ConvertWebMapToMapDocument(webMap)     mapDoc = result.mapDocument          df = arcpy.mapping.ListDataFrames(mapDoc,'Webmap')[0]        newExtent = df.extent            newExtent.XMin, newExtent.YMin = extentInfo.split(',')[0], extentInfo.split(',')[1]     newExtent.XMax, newExtent.YMax = extentInfo.split(',')[2], extentInfo.split(',')[3]     arcpy.AddMessage("newExtent 1")     df.extent = newExtent         When I try to run this I'm getting following error.               File "c:\program files (x86)\arcgis\desktop10.1\arcpy\arcpy\arcobjects\_base.py", line 87, in _set     return setattr(self._arc_object, attr_name, cval(val))     RuntimeError: DataFrameObject: Error in setting extent      Throws this error from this line. df.extent = newExtent


title: Why is a stored procedures much slower than the single query?

text: I wrote a normal query to snap the next way point from a pair of coordinates in EPSG4326, nothing complex. To optimize the performance a bit I made a boundary condition to minimize the points needed to be checked... Works smooth so far - but then I wrote a stored procedure to make later a neat intergeneration... (quasi wrapper). And here is where my headache starts:   * Single SQL: 650ms   * Stored Procedure: 45xxxms   * Single SQL is 70x faster than my storage procedure! Code below... * * *               CREATE OR REPLACE FUNCTION next_way_point(p geometry) RETURNS next_way_point_result AS $$     SELECT gid, source, target, the_geom, ST_distance(the_geom, ST_SetSRID($1,4326)) as dist      FROM ways      WHERE the_geom && ST_Buffer(ST_setSRID($1,4326), 0.1)     ORDER BY dist LIMIT 1 $$     LANGUAGE SQL;      Called like this:               SELECT * FROM next_way_point(ST_MakePoint(7.26, 47.52));      VS               SELECT gid, source, target, the_geom, ST_distance(the_geom, ST_SetSRID(ST_MakePoint(7.26, 47.52),4326)) as dist      FROM ways      WHERE the_geom && ST_Buffer(ST_setSRID(ST_MakePoint(7.26, 47.52),4326),0.1)      ORDER BY dist LIMIT 1;


title: Extrude between wrong on the edge, ArcScene

text: I'm trying to extrude between two Tins which are equally shaped. However if Extrude Between those tins using input feature created by minimal_bounding from one of them I get a strange error that two sides of my extruded block are inclined and thus the result is wrong (please have a look at the img attached). I have absolutely no clue what could be the reason. I would be very glad for comments. Thank you. Here are some images to illustrate the problem (thanks for the reputation so I can finally post them ;) ): ![a](http://imageshack.us/a/img833/4924/48164856.jpg) ![b](http://imageshack.us/a/img801/706/94039199.jpg) I'm using ArcScene with ArcGIS 10.1. This question can also be found on: http://forums.arcgis.com/threads/77930-Extrude-between-wrong-on-the-edge- ArcScene?p=273925#post273925


title: Display XY Format to Shapefile (Scale)

text: Using ArcView 10.1 and using excel for Lat/Long conversion to a shapefile. The original Lat/Long has 7 numerical characters past the decimal point xx.1143717. When I convert this excel into a shapefile, the numbers are shortened down to 6. The Lat of ex.. xx.1143717 when converted into a shapefile is brought down to 6 numbers past the decimal and rounded up to xx.1143720 Is there a way, or how can I make it where all the Lat/Long numbers are shown, all 7 instead of 6 in the attribute table?


title: Is there a way to turn on labels with ModelBuilder?

text: I have a model that creates a feature class and adds it to the map. The properties of the feature class reference a layer file for the symbology. However, I can't seem to find a good way to pass the labels as well so that the feature class is added with the labels visible. Is there a way to do this?


title: SLD Pointsymbolizer with color ramp based on first random attributes

text: The company I work at does drillings all over the country. These drillings are then presented as a service. So obviously this layer needs an SLD. For the moment every drilling is a green point, but we would like to go a bit more complicated. The latest idea was to the following: From a certain scale we'd like to let predetermined colors be filled for every dossier. One dossier can contain multiple drillings, but every drilling has a field which clearly notes that dossier. I've determined the scale, that's no problem, but I have no clue of how I can add a different color to every dossier grouped point currently in the screen. Anyone have a clue to how this could be worked out?


title: Creating ArcMap context menu's

text: actually i am in need of some requirment regarding context menus. We created our own button in the context menu which will call a exe file whenever that option in selected in the context menu. The problem is that this button is that it only should be enabled when a particular area on our map is selected. Otherwise it should be disabled. As i am new to this arcGIS ..need help on this.


title: Most efficient way to batch clip in FME?

text: I have one large shapefile and 100 smaller shapefiles that are smaller areas within the large shapefile. Using the Clipper transformer, I want to clip the large shapefile using each of the smaller shapefiles (I want to do the clip 100 times), and then write out 100 resulting shapefiles. What is the most efficient way to do this in FME?


title: erase one line layer by another line layer

text: I keep on having this issue and I can't find the "perfect answer": First of all, I work in ArcGIS 9.3.1. I have 2 layers: roads1950, roads2000. I need to make a layer with roads which were built/changed between 1950-2000. The most perfect thing was to have a "erase-like" tool that would cut out the identical to roads1950 features in roads2000. The problem is that the erase tool only works with polygon layers. I could do a merge and sow out the doubles but it's a (very) long process and it isn't accurate. Do you know of a tool/method that could help me ?


title: Creating animated graphics to symbolize selections in Arcgis Explorer

text: I want to create and add some animated graphics to the features based on a query in Arcgis Explorer, like some moving waves at the queried features. Is it possible? If so, please let me know. Thank you.


title: Tiled map service base layer does not display properly in ESRI Javascript map

text:
I am adding a map to a website built in asp.net and Javascript. When I add the esri.Map using a tiled map service layer base map, the map overlaps the map div and what I assume are tiles show up out of order and with spaces between them. Picture a large checkerboard with transparent squares alternating with map squares and Antarctica in the map square to the right of North America and Greenland at the bottom. There also appears to be some metadata text appearing in the div in a transparent square. There is another overlapping div that is turned off and the map div is turned on in response to a button click to display the map. A dynamic map service layer displays properly, but it behaves a bit oddly on zoom, extending to take up the whole page during the zoom, then snapping to the div when the zoom is complete. The map div is being set in css to relative position, and I've tried also giving it a set size in pixels, with no change in behavior. Originally the div was not given any position or size attributes, but the behavior was the same. I get the same result in Opera, Chrome, and Firefox. Nothing shows up in Explorer, but I think that may be a security setting issue. I'll try to attach code and screenshot of the map. Below is the javascript code adding the map. I'll just note that I've tried various extents and spatial references, including copying and pasting from standard "add a map" tutorials, so this is just what's there now. They all look the same. I've also added the basemap with addLayer by url with the same results. Javascript function to add map:               function MapIt(msg) {          //display map box     $("div.selection-criteria").hide();     $("div.map").show();          map = new esri.Map("map", {         extent: new esri.geometry.Extent({ xmin: 38.705772, ymin: -84.820221, xmax: 39.610317, ymax: -83.683907, spatialReference:{wkid:4326} }),         zoom: 3,         basemap: "streets"     });      } ASP code for the content and divs, with the fieldsets for the div that gets hidden when the map is show obfuscated:               <asp:Content ID="Content1" ContentPlaceHolderID="contentMain" runat="server">
          <fieldset class="somefields">     </fieldset>          <fieldset class="someotherfields">     </fieldset>          <input type="button" id="btnMap" class="map-button" value="Map it!" />

      CSS div styles:               div#map{position:relative; width:800px; height:400px;}      ![map]: http://i.imgur.com/65yxznk.jpg


title: Why does merge/mosaic four DEMs into one give ERROR 01042?

text: I'm unable to merge 4 dems into one using ArcGIS Home 10.2. Raster Calculator, adding them, gives me the error: > RuntimeError: ERROR 010422: Error in processing the environment mask. I created a basic polygon and set that as the mask, doesn't help. Using Mosaic to New Raster doesn't error out, but converts all values to 255, so I get a monochromatic output. DEM to Raster won't even allow me to add the DEM, I'm guessing because it doesn't fit the USGS format. These are DEMs downloaded from the MiGDL website (http://www.mcgi.state.mi.us/mgdl/), which I put in a file geodatabase. Years ago I used to be able to use Merge for this, but apparently that functionality is gone. Any ideas? Thanks, appreciate it. Just basically want to create contours and hillshade.


title: Does QGIS Recognize Annotation Feature Classes from File Geodatabase?

text: My work has recently set up QGIS to access an Esri file geodatabase on a shared drive for my co-workers. I have not used QGIS much, but I have encountered an issue where the annotation feature classes are displayed as a polygon instead of text. The geometry type of all the annotations are shown as 'MultiPolygon' when a vector layer from the geodatabase is added to QGIS. When the annotations are added to ArcMap, they appear as they are meant to - as text. Am I doing something wrong? ![annotations_bad](http://i.stack.imgur.com/ic9f4.png)


title: Suggestion about the database to be used with OpenLayers

text: Can anyone suggest me which database to be used with OpenLayers? I have to store the values in a DB table and retrieve them (trivial tasks). which is better MySQL or PostGIS? Can anybody give benefits of these? Thanks in advance


title: How do I represent the whole Earth as a Polygon?

text: This is related to determining the interior of a Polygon. My difficulty is that if I specify a Polygon covering the Earth going from East to West, it will still have a boundary along the antemeridian, which would show up in azimuthal projections. For example, in GeoJSON:               {       "type": "Polygon",       "coordinates": [         [[-180, -90], [-180, 90], [180, 90], [180, -90], [-180, -90]]       ]     }      Is there a standard way to specify a whole-sphere Polygon, with no boundary? I've thought about using an empty exterior ring:               {       "type": "Polygon",       "coordinates": [         []       ]     }      The logic being that I could specify further rings as holes to be punched into the whole-sphere Polygon. Unfortunately, GeoJSON stipulates that LinearRings must have at least 4 points, so this would be straying from the specification. However, I'm still interested in understanding whether this approach is taken anywhere else.


title: Upgrading custom Java script tool (ver. 1.3 to 3.1 ) for ArcGIS Java Script API and getting an error

text: I have found useful tool here (All US States and counties on combo-box) its using Java script API 1.3 version. I have downloaded the tool and trying use but its gives an store error (using version 3.1) but when I change the version to 2.8 its working fine..so how do I implement it for 3.1 version. Is there any document or instructions ?? (Because Esri mentions that only upgrade the version no. ex. 2.8 to 3.1 ).. Any help will be great..and as usual thanks in advance : ) Software that I am using is ArcGIS Java Script API ver 3.0 & ArcGIS Server 10.1,Desktop 10.1


title: How to sync multiple OpenLayers maps on same page?

text: I have 4 maps in one page and want to sync them (zoomend and movened). This is what I ahve done so far: **Approach 1:**                   sync(map[0], map[1], map[2], map[3]);         sync(map[1], map[0], map[2], map[3]);         sync(map[2], map[0], map[1], map[3]);         sync(map[3], map[0], map[1], map[2]);         function sync(m0, m1, m2, m3) {             m0.events.register("moveend", m0, function() {                 m1.setCenter(m0.getCenter());                 m2.setCenter(m0.getCenter());                 m3.setCenter(m0.getCenter());             });             m0.events.register("zoomend", m0, function() {                 m1.zoomTo(m0.getZoom());                 m2.zoomTo(m0.getZoom());                 m3.zoomTo(m0.getZoom());             });         }      **Approach 2:**                  for (var i = 0; i < maps.length; i++) {             map[i].events.register("zoomend", map[i], zoomListener);             map[i].events.register("moveend", map[i], moveListener);         }         function zoomListener(e) {             for (var i = 0; i < maps.length; i++) {                 if (map[i] !== this) {                     map[i].zoomTo(this.getZoom());                 }             }         }         function moveListener(e) {             console.log(this.div.id)             for (var i = 0; i < maps.length; i++) {                 if (map[i] !== this) {                     map[i].setCenter(this.getCenter());                 }             }         }      If I use it for a single map, it is running fine. However it is becoming an infinite loop when used with multiple maps. How can I stop the calling map to stop for further zoom?


title: How to achieve more than basic labelling of Reference Grids (Graticules, Measured Grids) in ArcGIS Desktop?

text: I'm happy to be corrected, but at the moment I do not think it is possible to meet any of the following mapping requirements using either ArcMap or ArcPy:   1. Label measured grids using a format like 6828000mN where 68 is superscript, 28 is normal, 000m is superscript and N is normal   2. Label measured grids using a format that at least enables mN and mE to be included even if all the superscripting above cannot be supported   3. Label graticules using a Degrees, Minutes and Seconds (DMS) format where the minutes and seconds can be written as 00 rather than 0 Does anyone have any thoughts on how to achieve one or more of the above? Note that I have excluded Grid and Graticule layers from consideration here, because even if they can support these grid and graticule labelling scenarios they are not dynamic and I need the labelling above in a Maps On Demand (MOD) / Data Driven Pages (DDP) like application. _**As an update_** it is possible to do some of 1. and 2. using Corner Labels in a static map as long as it is displayed in a Projected Coordinate System. However, this hangs ArcMap (and thus prevents it being saved) when the Coordinate System of the Data Frame is switched to Geographic which is necessary for an application like mine that uses a Geographic Coordinate System to determine the appropriate MGA Zone to use for each map at runtime. Basic labels work for Projected Coordinate Systems just not the fancy ones required in 1.


title: How to subtract Geometry Point Polygons in OpenLayers?

text: I have two json arrays, one for background polygon and one for other polygon. And from them, I got two polygon layers.               backgroundLayer = new OpenLayers.Layer.Vector("Background Layer", {styleMap: backgroundStyleMap} );     ....     totalAreaPolygon = new OpenLayers.Feature.Vector(new OpenLayers.Geometry.Polygon([linearRing]))     backgroundLayer.addFeatures([totalAreaPolygon])      And same for other polygon. How can I subtract the other polygon from background polygon as shown in the picture (right polygon). ![enter image description here](http://i.stack.imgur.com/4Z6vC.png) This is one of the examples from OpenLayers. http://openlayers.org/dev/examples/donut.html


title: Handling geoExt popup when one or more layer fails to load

text: GeoExt info Popup is very easy to use but when one or more layer is missing or fails to load Avoiding Pink Tyles. The popup is not filled with the values. While I can see that request is sent to the server and response is there but popup is blank. Is there any way to handle this condition without building custom component?               OpenLayers.Util.onImageLoadErrorColor = 1;                 OpenLayers.Util.onImageLoadErrorColor = "transparent";                   OpenLayers.Util.onImageLoadError = function(){                      this.src = "/img/blank.gif";                 };          controls.push(new OpenLayers.Control.WMSGetFeatureInfo({                                   autoActivate: true,                     infoFormat: "application/vnd.ogc.gml",                     maxFeatures: 3,                     eventListeners: {                         "getfeatureinfo": function(e) {                                                    var items = [];                             Ext.each(e.features, function(feature) {                                 items.push({                                     xtype: "propertygrid",                                     title: feature.fid,                                     source: feature.attributes                                 });                             });                                  new GeoExt.Popup({                                 title: "Info",                                 width: 200,                                 height: 200,                                 layout: "accordion",                                 map: app.mapPanel,                                 location: e.xy,                                 maximizable: true,                                 collapsible: true,                                 items: items                              }).show();                                                    }                     }                 }));