title: Zoom in to see details

text: I have 2 shape files, one that shows all the counties on an island, the other shows all the villages on the island.   Here is the scenario:   My initial map shows the counties. When you click on a county, the map zooms to that county displaying the villages within. What is the best way to accomplish this using leaflet?


title: R.Neighbors with large neighborhoods?

text: I'm running into problems running the r.neighbors GRASS command in QGIS. Although it will run for small neighborhoods as soon as I increase my neighborhood above a certain number it becomes unbearably slow. I don't understand why Focal Statistics on Arcmap will do this in the blink of an eye while R.Neighbors takes ages?


title: How to generate wind direction raster using its U and V components in arcgis?

text: I have two raster layers representing the U and V components of a vector wind. I could generate the magnitude raster by applying `(u*u+v*v)` in Map Algebra in ArcGIS. Now I need to generate the direction raster with the following function:               if(u > 0 && v > 0,( atan(u/v) + 180),if(u < 0 && v < 0,(atan(u/v) + 0),if(u > 0 && v < 0,( atan(u/v) + 360), if(u < 0 && v > 0,( atan(u/v)+270)))))      but I don't know how to apply it in ArcGIS. Any help please?


title: INSPIRE XSD schema definitions

text: I was looking for `INSPIRE` `XSD` schema definitions in the Web, but nothing found. Only something like THIS. But that resource is "a little" messy... I wanted to use `xsd.exe` to autogenerate `C#` classes that could be easily serialized into `XML` that is `INSPIRE` compliant. The problem is in finding that schema file. Could you point me where to download it?


title: Location Data in TAIP Format to Geoevent Processor

text: We are testing whether we can get the surveyed cordinates in a Trimble GPS transferred online via TAIP format ( Trimble exchange protocol) to ArcGIS Geoevent processor . For this , one main concern is how the device to configure the Device to receive the data in TAIP format and secondly how to push the TAIP format to the Server. The page below provides some insight to the work flow. http://www.arcgis.com/home/item.html?id=6b458b1f4de5436486ffbd7801c6b438 The question is , how can I override the programs in the GPS to configure the TAIP format . I guess this qn should be put to the GPS hardware dealer, but if anyone can offer a suggestion, we make some headway..


title: PostGIS: ST_Clip() throws non-noded intersection between LINESTRING

text: I have a polygon and a set of rasters. I'm trying to clip raster with polygon.   ![enter image description here](http://i.stack.imgur.com/yKhDH.png)   As you see, my rasters is divided into tiles. Any tiles that is contained by the polygon has ST_Overlaps(polygon, rast_geom_4326) = 'f', obviously. And it leads to non-noded intersection error. But, I also have the same problem with this query:               SELECT ST_Clip(rast, ST_Transform(polygon, rast_geom_4326))     FROM landsat8     WHERE ST_Overlaps(polygon, rast_geom_4326)      Is there any function to dissolve the polygon with tiles, to get rid of "non- noded intersection between LINESTRING".   I've read a presentation about this error postgis-power.pdf. But, I can't find the solution for my problem.


title: How can I filling ring in QGIS 2.2?

text: I created a polygon which covers a park. But I want to distinguish some small buildings located inside the park. Consequently, I want to distinguish them and use them as independent patches inside the original polygon. I used do it the "Add ring" function in the QGIS 2.2. But I don't know, how can fill those empty holes, which were resulted. I tried the "Fill rings" options (in QGIS 2.2), but I was not lucky with it. Does anybody have an idea, how can I fix this problem?


title: How to add pages in print composer

text: Is it possible to add more than one page to one QGIS print composer Project? For example one sheet for the map and another sheet for the legend ans some other informations. Thanks in advance, Patrick


title: PostGIS SQL query to find a polygon that intersects with a line

text: I need to build a PostGIS query that would select every polygon inside an area the river is flowing through. I'd better describe with an example below: There is a river ( _blue line_ ) flowing through several violet polygons ( _selected in yellow_ ). Both of the yellow polygons are crossed by the green one (labeled Žulová, which is a boundary of a county), but the river actually flows through the only one of them (labeled 2-04-04-044) and the other one _is only touched_ by the end point of the river. The query I have so far looks like this:               SELECT DISTINCT              poradi              ,ST_Distance(ST_StartPoint(ST_LineMerge(toky.geom))             ,ST_Centroid(pov.geom)) AS dist     FROM              povodi_iv AS pov             ,cr_osy_toku_snap_2m AS toky             ,$table     WHERE              toky.gid = $gid // river id         AND                   $table.uid = $id // county id         AND                   ST_Intersects($table.geom, pov.geom)         AND                   ST_Intersects(toky.geom, pov.geom)     ORDER BY              dist ;      But it gives me the result you see on the picture. Could anyone help me out to build a query that would **return only the violet polygons that are both crossed by county boundary and flown through by a river?** ![Example of what I want to achieve](http://i.stack.imgur.com/txMWP.png)


title: How to extract values from rasters at location of points in R?

text: My question is about extracting values from rasters at the location of points. With the function extract this is very easy, and the function gives me a dataframe with the values of all the variables in the points. I want to have in that dataframe the coordinates of each point. How can I make that happen? Is it possible to say to R that when extracting the values from the raster also add the columns of the location point. This is my procedure:               presencias=read.table("c:/SDM_R/presencias/P_lentiscus_pres.csv",header=TRUE,sep=";")          lista_variables <-list.files(path="Variables_modelizacion/solo_ascii",pattern='*.asc',full.names=TRUE)     variables <- stack(lista_variables)          variables_presencia<-extract(variables, presencias)      the results are something like these:               > bio1  bio12  bio18  bio2  bio4      > 90    875    165    95    4886     > 115   1085   158    83    4075     > 135   1153   153    67    3402     > 85    1026   137    99    5203     > 96    667    128    108   5823     > 98    531    109    113   6305     > 132   450    63     123   6598     > 132   569    104    106   5963     > 95    814    196    98    5571     > 146   474    39     114   6603      But I want two more columns with the coordinate data (but could be others columns from my csv extracting location table). Thank so much.


title: OpenLayers/Leaflet render world borders

text: I am trying to show all world borders. Either using OpenLayers or Leaflet. World borders are taken from here. This code does not work. What I am doing wrong?               var Map = {            map: null,       layer: null,            initializeMap: function() {         map = new OpenLayers.Map('map');         layer = new OpenLayers.Layer.OSM("World Map");         map.addLayer(layer);         map.zoomTo(2);         var geojson_format = new OpenLayers.Format.GeoJSON();         var vector_layer = new OpenLayers.Layer.Vector();          map.addLayer(vector_layer);         vector_layer.addFeatures(geojson_format.read(featurecollection));        }          }


title: GPL-compatible database of country & region names

text: For WordPress's sister project BuddyPress, we are looking for relatively- reputable resources that manage a list of countries and regions. Due to the WordPress OpenSource philosophy, we are only seeking GPL-compatible resources. If no such resources exist, we could possibly consider GPL-compatible services as well. Or there could possibly be a GlotPress-like service (see http://translate.wordpress.org/projects/wp/dev) wherein the community contributes and maintains the database which can be downloaded via a link or API. So far all the databases/services I have found have either been commercial or otherwise incompatible with GPL.


title: How to summarize GIS skills on a resume?

text: I'm a recent graduate and I've used ArcGIS 10.0-10.2 (I just say 10.X) for a variety of tasks:   * Analysis   * Data Conversion   * Georeferencing   * Editing   * Data Management (i.e., creating and organizing geodatabases in ArcCatalog)   * General map creation Does anyone have any suggestions on how best to summarize my GIS skills on my resume?


title: Calling a AddIn from another class or AddIn?

text: I've tried several ways but I can't find anything working. Is there any possebility to call an AddIn from outside? I want to fill a AddIn Combobox when the OnClick-Event of an AddIn Button is fired, but don't know how to adress the Combobox. Have you any Ideas?


title: How many WFS-T clients exist?

text: I would like to know which applications or libraries (proprietary or open source) offer transactional Web Feature Service (WFS-T) client functionality. Those include:   * Desktop or online applications that offer the functionality "out of the box", or   * Libraries that can be used to develop a WFS-T client, like OpenLayers.


title: How to use r.patch with two maps having different colortables?

text: I am using r.patch to combine two raster maps in order to fill in the null values of the other. Is there a way to use r.patch but still maintaining the colortables of the separate raster maps? I do not want to change the colortable of input1 raster map A and input2 raster map B.


title: Background reading on using GPUs for GIS?

text: I'm interested in the emerging intersection between GIS systems and GPUs, which can provide orders of magnitude improvement to certain classes of GIS problems. Do you know of any good resources discussing this area?


title: export shapefile to spatialite hangs qgis

text: I have a shapefile layer containing 9 points. The attribute table has a string column (max 18 chars) and a real column. The .dbf file is ~4KB, all the other files are 1KB. I right-click on the layer, select Save As and set the type to Spatialite, assign a file name, then click OK. The process just runs and runs. No error. If I look at the file, it's growing about 2KB every couple of seconds. I killed it when it passed 2.5MB. I've tried this a couple of times and the same thing happens consistently. Any suggestions about what to try next? [edit]   The problematic shapefile can be found here QGIS version 2.0.1-Dufour   QGIS code revision d94c044   Compiled against Qt 4.8.5   Running against Qt 4.8.5   Compiled against GDAL/OGR 1.10.0   Running against GDAL/OGR 1.10.1   Compiled against GEOS 3.4.2-CAPI-1.8.2   Running against GEOS 3.4.2-CAPI-1.8.2 r3921   PostgreSQL Client Version 9.2.4   SpatiaLite Version 4.1.1   QWT Version 5.2.3   PROJ.4 Version 480   QScintilla2 Version 2.7.2 Windows 7 Pro, 64 bit


title: How to find all features in one table that intersect features in another table in PostGIS?

text: I have two tables, `boundary` polygon table and `location` point table. I need to find all `boundary` polygons that contain at least 1 `location` point. What is the proper way to do this in PostGIS? I tried:               SELECT *     FROM boundary, location     WHERE ST_Intersects(boundary.geom, location.geom)      But this returns each combination of `boundary` and `location` that intersect, when all I want is the `boundary` returned once. For example, I got 76 results, but only have 10 `boundary` entries. Do you know a better way to do this?


title: ArcObjects table views or sqlite3 in Python?

text: Given a larger csv (1000+ points) with attributes, I need to do both attrbute and spatial queries on the data. Would it be faster to read the points into an xy feature class, and do the spatial and data queries as table views and feature layers, or would it be faster to read the csv into a temporary sqlite3 database, do the data queries, then transform the output into an xy class to do the spatial queries?   Also, would it be faster one way or the other if we also needed to include table joins?   Lastly if you have any alternate methods please add them.


title: What's up with OpenJump?

text: OpenJump is a nice GIS, very functional, with many updates and has a beautiful rendering. But it's not very popular, and you don't see it around in discussions, in osgeo comunity, etc. And gis.stackexchange doesn't even have a tag for it.   So, what's up with OpenJump? And what keep it's development? ![enter image description here](http://i.stack.imgur.com/XYQyq.jpg)


title: How to save and reload legend settings using Python?

text: Maybe you can help me with maybe a strange question. I know how to change settings of a legend individually by arcpy (like: Patch, Title, Items, Columns, Frame, Element Size&Position, etc.), but now I would like to use the exact legend settings another time and load them using arcpy. Is there a simple way to save/load these legend settings all together in/from a text file? Thanks in advance for any ideas.


title: Variable XML dataset as reader in FME

text: I have one problem i am thinking how to solve. My workbench is started from python script. Workbench reads xml file, but name of xml file is each time different. So is it possible to solve this with generic reader or with another python script which will read folder where is xml file, find the name of xml file and proceed it to the fme as a reader? Any other idea?


title: OpenLayers for Android: why does a popup panel show up twice?

text: I use Moto Xoom (Android V3.0 touch pad) browse the example "OpenLayers with Sencha Touch", when I click button "LayerList", the associated popup panel doesn't show up smoothly. it will show once, then very quickly show up again. I tested other buttons, for example "Search" button and popup "search" panel, same issue happens. It seems that it's something related to OpenLayers. Because another example "OpenLayers with JQuery Mobile" also has this. I am wondering whether someone could give me some hints about what's happening underlying this issue? appreciate any replies!


title: Access-Control-Allow-Origin problem when loading json file from Amazon S3 to display in OpenLayers

text: This is my code for displaying geojson in OpenLayers, where "geourl" is directed to a json file in Amazon S3.               geojson_layer = new OpenLayers.Layer.Vector("GeoJSON", {         strategies: [new OpenLayers.Strategy.Fixed()],         protocol: new OpenLayers.Protocol.HTTP({         url:geourl,         format: new OpenLayers.Format.GeoJSON()})         });      But I received the error like this:               > Failed to load resource: the server responded with a status of 403     > (Forbidden) "geourl" Failed to load resource: No     > 'Access-Control-Allow-Origin' header is present on the requested     > resource. Origin 'http://localhost:8081' is therefore not allowed     > access. "geourl" XMLHttpRequest cannot load  "geourl" No     > 'Access-Control-Allow-Origin' header is present on the requested     > resource. Origin 'http://localhost:8081' is therefore not allowed     > access.      Nothing changed when I added               self.response.headers.add_header("Access-Control-Allow-Origin", "*")      to my webapp RequestHandler Class.


title: Sort Grid binded to ArcDataBinding.TableWrapper in ArcObjects 10.1

text: I have a Datagridview which is binded to a ArcDataBinding.TableWrapper. So, I want it to be sorted when ColumnHeader is clicked. I tried to do the following on Click event:               ArcDataBinding.TableWrapper twapper = (ArcDataBinding.TableWrapper)((BindingSource)datagridview.DataSource).DataSource;                     ITableSort tsort = new TableSortClass();                     tsort.Table = (ITable)twapper.wrappedTable;     ...                     string column_name = datagridview.Columns[e.ColumnIndex].Name;                     tsort.Fields = column_name;                     tsort.set_Ascending(column_name, false);                     tsort.set_CaseSensitive(column_name, true);                     tsort.Sort(null);     ...                     ArcDataBinding.TableWrapper tableWrapper = new ArcDataBinding.TableWrapper((ITable)tsort);                               tableWrapper.UseCVDomains = true;                     BindingSource bindingSource1 = new BindingSource();                     bindingSource1.DataSource = tableWrapper;                      datagridview.DataSource = bindingSource1;                          datagridview.Refresh();      But it doesn't seem to have any effect. Does anyone knows what is the problem?


title: How to perform a Con( ) evaluation like "If "Inras1" is true then set "inras2" to Value X " correctly in ArcGIS?

text: I want to set infiltration-rates to different Landcover depending on the slope. I have a Landcover Class _Greenalder_ and six slope- classes with different infiltration rates. e.g. greenalder in slope class from 40-50° has a infiltration rate of 0.774mm. In ArcGIS I have a Slope- layer and this Greenalder Layer (shrub_07). Now I want to set the Value of the Greenalder Layer to 0.774 if it is between 40 and 50° slope. I tried it with the following syntax in the RasterCalculator:               Con(("slope" >= 40.0001) & ("slope" <= 50),"shrub_07" == 0.774)      but the result is not 0.774 but always 0. I tried a lot of other possibilities but nothing worked.


title: How to automatically transform between base layer and vector layer in OpenLayers

text: I have a simple example on http://jsfiddle.net/29bsw/ I have a map with an OSM layer as base layer and an additional vector layer. The projection of the map is defined as EPSG:900913 and the projection of the vector layer is defined as EPSG:4326:               //Map     var map = new OpenLayers.Map("map", {         projection: new OpenLayers.Projection("EPSG:900913")     });          //Base Layer     var osm = new OpenLayers.Layer.OSM();        map.addLayer(osm);          //Vector Layer     var vectorLayer = new OpenLayers.Layer.Vector("Simple Geometry", {         projection: new OpenLayers.Projection("EPSG:4326")     });     map.addLayer(vectorLayer);        Now I would like to add some features to the vector layer without transforming any geometry point from EPSG:4326 to EPSG:900913. The point should be lie in Hamburg, Germany, but without transforming (see commented on line 21 in the example) the point lies near by Africa. Because I specified the projection of the map and of the vector layer, is there any way that OpenLayers transform the point from EPSG:4326 to EPSG:900913 automatically?


title: Failed to connect to the specified server, ArcGIS 10.1 connection using Oracle 11 g

text: I am trying to connect an Oracle 11g database to ArcGIS 10, but it's not working and gives the following error: "Failed to connect to the specified server. Operation Failed." While I am sure about the inputs, and it's as the following: Database Platform: Oracle Instance: localhost/dbName Username and password are correct. I tried to connect using sqlplus and it worked properly. I noticed in some posts that if my Oracle is 64-bit, then I should install Oracle client 32-bit and put its directory as the first one in the PATH directory. I've done all the above and nothing worked, so if there is anything that I'm missing please let me know. Thanks in advance.


title: Adding XY coordinates from a point feature class to a one to many related table in ArcGIS 10.0

text: I have a point feature class with a one-to-many relate to a table. I would like to add the XY coordinates of each point in the feature class to the related table where the ID field matches. Of course, it's simple to get the XY coordinates for each point but it's adding them to the table where I'm getting stuck. Thanks in advance


title: Can REST or the ArcGIS API for Flex read VAT tables from a Raster?

text: When publishing a bare earth DEM that has a VAT table, or value attribute table _"A VAT is a table containing attributes for a grid, including user- defined attributes, the values assigned to cells in the grid, and a count of the cells with those values"_ ...The fields are empty in the REST service where there should be values. Anyone know how to consume the VAT table that resides inside the gdb.? I want to be able to use the identify tool to depict elevation values.


title: Architecture issue - ArcGIS Server sourcing data from FME?

text: I'm looking at potential architectures for a new project and the biggest unknown right now is how well ArcGIS Server can communicate with non-ESRI data sources for map services. There are a number of data sources including Oracle Spatial and SDF files. There is no SDE and periodically converting source data into shapefiles or similar is not an option as the data must be live. It has been suggested that FME (server?) may be able to mediate between the two, but this is something I have never tried. Does anyone have experience with this or a similar setup? I'd be interested to know how stable a configuration this is, whether editing the source data would be possible (assuming a capable API was used on the client), or any other alternatives. Any suggestions much appreciated.


title: Monthly Average Solar Maps (R.Sun) in Grass

text: I'm trying to create monthly average maps as described in R.Sun in Grass. I have my DSM. In reading the page on OSGEO it mentions that it is possible without showing how. Ideas? I assume I will need to write a script.


title: ArcIMS - Active layer update

text: I'm trying to change the active layer in ArcIMS from the default top position to index 2 in the TOC. The TOC has 3 grouped layers, but it is the first grouped layer that I want the active layer to be updated to position 2 not the default position 0. I found this code link from esri showing how to update the active layer, however I can not seem to figure it out. Here is the code from the aimsLayers.js file.               //ActiveLayerIndex=<number of the preferred active layer>;         ActiveLayerIndex=1;              if (ActiveLayerIndex>=layerCount) ActiveLayerIndex = 0;         if (!LayerIsFeature[ActiveLayerIndex]) {             var chk = 0;                  for (var i=layerCount-1;i>=0;i--) {                 if (LayerIsFeature[i]) chk = i;             }             ActiveLayerIndex = chk;         }      I added the "ActiveLayerIndex=" and changed it to " **ActiveLayerIndex=1;** " the active layer moved one down to position 1 in the TOC. When I change it to 2 it did not move to position 2 in the TOC and displays no active layer. When I change the index in this line, **if (ActiveLayerIndex>=layerCount) ActiveLayerIndex = 0;** , it does not update the active layer either. How can I get the active layer to be on the third layer (position 2) in the TOC? **UPDATE** I found another code snippet and was wondering by using it to point to the specific layer name in the TOC if it will re-assign the active layer? Looks like it could be called within the TOC.htm               //global reference to mapframe         var t = parent.MapFrame;              function layerName(theFormRef) {             var theLayerNumber = -1;             var theLayer = theFormRef.value;                  for (var i=0;i<t.LayerID.length;i++) {                 if (t.LayerName[i] == theLayer){                     theLayerNumber = i                 }             }             if (theFormRef.checked) {                 t.LayerVisible[theLayerNumber] = 1;                 //activating layer, NB! raster layers needs to be excluded                 setActiveLayer(theLayerNumber);             } else {                 t.LayerVisible[theLayerNumber] = 0;             }             t.sendMapXML();         }      http://forums.esri.com/Thread.asp?c=64&f=786&t=278849 Thanks


title: Changing Names in TOC Doesn't Update Legend Properly

text: This is related to an ArcMap add-in I coded up this past week. I'm trying to rename a layer in ArcMap and have the legend refresh properly. The renaming part is working fine, however I cannot get the legend in the layout to update properly. I call CurrentContents.Refresh() to refresh the name in the table of contents, which works. If I search through the map surrounds and call .Refresh() on the legend, then do an ActiveView.PartialRefresh() on it, the legend gets all out of whack. This is why I am trying to simply trigger the same event as when the user manually changes the name of a layer. The API documentation seems to say that MxDocument.UpdateContents() will trigger the events needed to refresh the legend. However, this does not seem to be the case. In addition to the above methods I have tried ActiveView.ContentsChanged() to no avail. Edit: I am writing this in C# but can read VB if you feel more comfortable writing in that. Any suggestions are welcome. Thanks, Mat


title: PtDistWithin() with Geometry from query

text: I'm trying to run two Spatialite queries in Python — the first one picks a single point from a table, and the second tries to find all points that are within a given distance of it. So my queries look like:               select geometry from points where point_id = 1      and               select count(*) from points where      PtDistWithin(points.geometry, <geometry from first query>, 100)      I'm stuck on how to inject the geometry that I selected in the first query into the second one. If I drop it in as-is it says cannot concatenate string and buffer objects. I tried selecting it `AsText()` and there are no errors, but the count is always 0 (there should definitely be matches within 100 meters). Does anyone know a better way of doing this? Thanks!


title: Better-matching bounding box by using collection

text: I am looking for a way to approximate a geographic region with a few bounding boxes. I know I can get a single box to cover the whole region, but some of my shapes are quite strange (gerrymandered, in some cases). I'd like to know if there's a way to approximate a shape as a (small) collection of axis-aligned rectangles to make a tighter collection of bounding boxes. For example, I know I can draw a rectangular box around a roughly "C" shaped region, but I can make a much better approximation with 3 rectangles like so:                   +--------------+         |              |         +----+---------+         |    |         |    |         | +~~+--------+         +-+~~+        |           +-----------+      ### Constraints:   * The rectangles should contain the entire original shape.   * It's fine if the rectangles overlap. ### Trade-off:   * I'd really like to stay under 10 rectangles. (Maybe an increased penalty for each additional rectangle?)   * I'd also like to have under 25% of the original shape's area as extraneously included by the rectangles, but if it takes too many shapes to do so, I'll go with fewer shapes. Is there a recognized method for doing this to arbitrary shapes? To certain classes of shapes?


title: issues using ogr2ogr to reproject and rename into shp

text: I'm trying to reproject, move from tab to shp and rename all the .tab contained in a folder tree using ogr2ogr. I'm working with Windows 7, QGIS 1.8 installed with advanced OsGeo4w and gdal 1.9.2. I have no experience in batch. Here's the code :               for /R %%f IN ('D:\3.BATCH\FAUNE\*.tab') do C:\Users\Samy\Desktop\SIGMA\15.webmapping\ms4w\tools\gdal-ogr\ogr2ogr.exe -s_srs EPSG:3727 -t_srs EPSG:2975 %%f "ESRI Shapefile" %%f_2975.shp %%f     pause      By calling this .bat into the cmd line, nothing append. When I simplify the code, I have in the best case "Unable to open" as Nik posted in the 3rd ressource below. I wrote this code with the following questions:   1. ogr2ogr convert all shapefiles in a directory   2. Is there any way to reproject multiple layers in QGIS?   3. How to reproject shapefile using ogr2ogr According to resource 1, I translate the -f "MapInfo File" "%~dpnf.tab" "%f" into %%f "ESRI Shapefile" %%f_2975.shp %%f because as it is shape, I can't use dpnf. According to resource 2, I wrote -s_srs EPSG:3727 -t_srs EPSG:2975 to reproject. Finally I checked that the "Hide extensions for known file types" was off, according to ressource 3. I don't know what to do now. There seems to be so many different ways to write this one line of code that I'm kind of disappointed. Please help ! Samy


title: Python script to read coordinate system of an input feature class?

text: How to read the coordinate system of an input feature class in python script.   How we can use "Describe()" function and tool validator for that purpose


title: How to perform a view shed analysis?

text: I want to run a viewshed analysis on a site I am interested in. I have the relevant elevation raster (taken from the OS UK, which is in the ASCII format) and a point feature class from the lat/lon coordinates of the particular site (I have also converted both into the same coordinate system). However when I try to run a viewshed analysis in ArcGIS 10.1, it doesn't work and there is an error message saying that 'there is no valid observation points present'. I was wondering whether anyone could tell me what I am doing wrong and how to correctly create observation points and run the viewshed analysis. I have looked at some books on the subject and have successfully followed their instructions on the data provided by them but I seem to be having trouble when I am using my own data! Any information would be massively appreciated! I have used the Geographic Coordinate system for both: OSGB_1936. I have taken the DTM from OS UK website (which is called Panorama and it is in ASCII format). First I added both to ArcMap and created a Hillshade. Hillshade did work but there was a yellow exclamation mark and obviously I got the error message I mentioned earlier when I tried to do the viewshed. Today I changed the coordinate system of the DTM to the British National Grid (that’s under Projected coordinate systems) but left the XY point the same: Geographic Coordinate system OSGB_1936. When I ran the Hillshade with this combination it worked properly i.e. there was no yellow exclamation mark. Then I ran a viewshed and it seems like this time it worked but the only thing I did different was change the coordinate system for the DTM. I don’t understand why that would make such a difference and I was hoping whether you would be able to let me know if I have done the right thing here (I did try to attach printscreen shots to make the steps I have done clearer but it didn't work - so, sorry about that). I wanted to add 15m and run another viewshed. So I set a OFFSETA (which appeared on the attribute table of my XY point) and another viewshed was generated - indicating that the area of visibility had increased which seems to suggest that it worked. But I was wondering whether you might be able to tell me why changing the coordinate system may have had such an impact and whether I have missed any crucial steps along the way?


title: Qgis 1.9 commitChanges() and Redraw a temporary layer

text: Using QGIS 1.9 and writing a python plugin. I have a function that connects points into polylines and draws them as a temporary, in memory layer. Membership in a polyline is determined by a field in the attribute table (postgis enabled table) of the point. The tool programatically sets up an edit session and alters the correct column in the point attribute table to alter the polylineID. pseudocode:               layer.startEditing()          ...Make edits to the features (not a data provider)...          layer.commitChanges()     remove in-memory layer     redraw in-memory layer by iterating over the point table to building polylines      The redraw fails to see the edits that have occurred. This is after calling `layer.commitChanges()`. What is the proper method to make an edit and propagate it so that an immediately derived layer sees the changes? Tested with data in geojson on disk and via a view on a postgis table.


title: 'DEFeatureClass' vs. 'IFeatureClass'

text: I'm writing an ArcGIS 10.1 .pyt toolbox to accept a feature class as input (not a shapefile). Referencing ESRI's online guide to 'Geoprocessing Data Types of Parameters and Environments' [PDF warning], I can set the datatype to either 'DEFeatureClass' or 'IFeatureClass'. What's the difference between these two datatypes? Which should I use and why? Example code from the getParameterInfo method:               treePolys = arcpy.Parameter(         displayName = "Input tree polygons feature class:",         name = "tpolys",         datatype = "DEFeatureClass",          parameterType = "Required",         direction = "Input")


title: Where can I find municipal-level administrative boundary data for Poland?

text: I am looking for municipal-level administrative boundary data (NUTS 5) for Poland. Since it is public information, I thought the data would be readily available, but I can't seem to find it anywhere. I have checked the usual sources. GADM has data down to the county (powiat) level, and Natural Earth seems to only have data down to the provincial level. I have looked at many of the Polish government's sites, but while I can find plenty of GIS data viewers, I can't find a way to download the administrative boundary data I am interested in. To make matters more complicated, I would really like historic data, but I would settle for current NUTS 5 level data.


title: GeoTIFFs made in GDAL don't have min/max range when loaded in Arc/QGIS

text: I'm creating some floating point rasters in GDAL with the GeoTIFF driver. When I load the resulting images into QGIS or Arc the default symbology puts the min at -3.40282e38 and max to 3.40282e38, so the raster looks exactly gray. Is there a way I can write the actual range into the GeoTIFF directly so when its loaded into a GIS program it automatically scales the histogram nicely? I've tried creating a default histogram like this:               rasterMin, rasterMax = raster.GetRasterBand(1).ComputeRasterMinMax()     raster.GetRasterBand(1).SetDefaultHistogram(rasterMin, rasterMax, 255)      But I get a error that looks like a binding issue:                 File "/usr/lib/python2.6/dist-packages/osgeo/gdal.py", line 846, in SetDefaultHistogram         return _gdal.Band_SetDefaultHistogram(self, *args)     TypeError: not a sequence      Any suggestions on what I could do, or if I've made a mistake?


title: What is Quickest way to Copy Single Selected Polygon from Very Large ArcSDE Feature Class?

text: I have selected a single polygon, which is large but not massively complex, from an ArcSDE 9.3 for Oracle feature class with 3.5 million polygons using ArcGIS Desktop 10.0 SP3. I would simply like to copy it out via arcpy into a new shapefile (or file geodatabase feature class) and to do so in less than 10 seconds. So far I have not been able to get it to happen in less than about 5-6 mins. I think I need to try and take advantage of the processing environment (arcpy.env.extent) during the CopyFeatures but if I set that to the values returned by getSelectedExtent and do the copy the time taken still seems to be 5-6 mins. If instead I could get at a centroid x,y (and one that is guaranteed to fall inside the polygon) then I could use that to buffer a tiny amount and use that to set arcpy.env.extent and I would expect the time to come down to around 10 seconds. However, to get at the centroid I think I still need to use CopyFeatures to get it out into a geometry first which seems Catch 22. (FYI The 10 second target comes from a **_different question_** based on the same dataset)


title: How to count and record the number of entities with the same name in the field calculator?

text: I have a question with ArcGIS 9.2. In an attribute table the column [name] includes the name of every place in this layer. Next to that (in columns [numbers]) there should be recorded how often a place of this name is in the map. How can I calculate this values with VBA in the field calculator? I am really happy about every answer and every help. Thank you.


title: Resolution difficulties in Grass

text: I have gone through similar questions to mine here, but I didn't find a solution. Please help! When cropping a raster image by using `g.region` & `r.mapcalc`, I get a new image but the resulting image doesn't have the same resolution. I have used `r.resample` in the outer raster map but is still different. I have tried to change the resolution using `g.region`, and even though the resolution is the same, the image display is not. Is there a way in grass to cut an image from the main raster image and keep the resolution unchanged? Also, I have tried using zoom to choose the crop area and its works! but as I have to do the same with a lot of images and need to be precise (when zooming), so it is difficult and takes a lot of time. How could I solve this? Many thanks!


title: Are .shp files coordinates in a certain direction?

text: I'm working on a project where I'm importing a .shp file into SQL Server. I'm using the Catfood.dll and I'm wondering if each point in the polygon are sequential (going in the same direction)?


title: Custom tiles in local CRS

text: I want to achieve the following:   * Create an OpenLayers map, that runs totally in a local CRS.   * The basic layer would be a raster tile layer.   * The tiles will be generated by us once for each project (based on TIFFs or other images, provided by our customers) and stored in the file system (the area is very small (~ 10x10 km) and there will be none or not many redundant tiles).   * Make the tile generation process as easy, as possible. We can assume, that the images, provided by our customers have the right orientation and it should be okay, just to say: The upper left pixel of the map is x1/y1 and the lower right pixel is x2/y2 to have our local CRS. Projections are irrelevant for us here. My question is: Is it possible in OpenLayers, to create a base layer, define the URL of the tiles and say: "the upper left corner of the first tile (0/0) is at x: 5271 and y: 12 in my local coordinate system and another point at tile XYZ is at x: 5427 and y: 634". Or do I have to bring the tiles in a predefined format? Or do I have to create a custom projection definition or something like that?


title: What's the projection of this albanian military map?

text: I've received some very useful maps of Albania for my work (scale 1:25000) and I would like to georeference them in order to put them together to one Map for my destined regions. Unfortunately I cannot find the projection of the maps, they don't seem to be written on the map (it is unfortunately in Albanian language). Maybe you guys have some ideas which projection it might be? (I'm using Qgis btw.) Thanks in advance! ![enter image description here](http://i.stack.imgur.com/eNzEZ.jpg) full image


title: Error 000368 selectbylocation in script from model builder

text: Can anyone help? I created a model-has buffer,clip, selectbylocation.... I exported the model as a script but when i try to run it i get this error:               Executing: IGISprogrammingscript C:\Temp\VirtualCampus\IGIS\Thames.shp "C:\Temp\VirtualCampus\IGIS\Super Output Areas.shp" "C:\Temp\VirtualCampus\IGIS\Points of interest.shp" C:\Temp\VirtualCampus\IGIS\Thames_Buffer.shp C:\Temp\VirtualCampus\IGIS\poi_test.shp C:\Temp\VirtualCampus\IGIS\superoutputareas_clip_test.shp     Start Time: Sun Jul 07 20:35:37 2013     Running script IGISprogrammingscript...          Traceback (most recent call last):       File "C:\Temp\VirtualCampus\IGIS\IGISprogrammingscript.py", line 28, in <module>         arcpy.SelectLayerByLocation_management(Points_of_interest, "COMPLETELY_WITHIN", Thames_Buffer_shp, "", "NEW_SELECTION")       File "c:\program files\arcgis\desktop10.1\arcpy\arcpy\management.py", line 6585, in SelectLayerByLocation         raise e     ExecuteError: Failed to execute. Parameters are not valid.     ERROR 000368: Invalid input data.     Failed to execute (SelectLayerByLocation).               Failed to execute (IGISprogrammingscript).      Also i have another question related to this, in my model i have several tools like clip,buffer,selectbylocation....when i run my script with the model parameters it only runs a section of the script, do you know why? Thanks for your help                   # -*- coding: utf-8 -*-         # ---------------------------------------------------------------------------         # IGISprogrammingscript.py         # Created on: 2013-07-07 16:55:49.00000         #   (generated by ArcGIS/ModelBuilder)         # Description:          # ---------------------------------------------------------------------------              # Import arcpy module         import arcpy                   # Local variables:         Thames_shp = "C:\\Temp\\VirtualCampus\\IGIS\\Thames.shp"         Super_Output_Areas_shp = "C:\\Temp\\VirtualCampus\\IGIS\\Super Output Areas.shp"         Points_of_interest = "Points of interest"         Thames_Buffer_shp = "C:\\Temp\\VirtualCampus\\IGIS\\Thames_Buffer.shp"         superoutputareas_clip_shp = "C:\\Temp\\VirtualCampus\\IGIS\\superoutputareas_clip.shp"         Points_of_interest__2_ = "Points of interest"         poi_shp = "C:\\Temp\\VirtualCampus\\IGIS\\POI.shp"         poi_cat = "C:\\Temp\\VirtualCampus\\IGIS\\poi_cat"         pop = "C:\\Temp\\VirtualCampus\\IGIS\\pop"              # Process: Buffer         arcpy.Buffer_analysis(Thames_shp, Thames_Buffer_shp, "300 Meters", "FULL", "ROUND", "NONE", "")              # Process: Select Layer By Location         arcpy.SelectLayerByLocation_management(Points_of_interest, "COMPLETELY_WITHIN", Thames_Buffer_shp, "", "NEW_SELECTION")              # Process: Copy Features         arcpy.CopyFeatures_management(Points_of_interest__2_, poi_shp, "", "0", "0", "0")              # Process: Summary Statistics         arcpy.Statistics_analysis(poi_shp, poi_cat, "CATEGORY COUNT", "CATEGORY")              # Process: Clip         arcpy.Clip_analysis(Super_Output_Areas_shp, Thames_Buffer_shp, superoutputareas_clip_shp, "")              # Process: Summary Statistics (2)         arcpy.Statistics_analysis(superoutputareas_clip_shp, pop, "LSOA_POP SUM", "")      Hi, i've added make feature layer and make table view and i'm still getting an error but it's a different one this time...but still refers to the make feature layer. This is the code:               ~# -*- coding: utf-8 -*-     # ---------------------------------------------------------------------------     # IGISprogrammingscript_1.py     # Created on: 2013-07-09 21:39:46.00000     #   (generated by ArcGIS/ModelBuilder)     # Usage: IGISprogrammingscript_1 <Thames_shp> <Points_of_interest> <Super_Output_Areas_shp> <Thames_Buffer_shp> <superoutputareas_clip_shp>      # Description:      # ---------------------------------------------------------------------------          # Import arcpy module     import arcpy          # Script arguments     Thames_shp = arcpy.GetParameterAsText(0)     if Thames_shp == '#' or not Thames_shp:         Thames_shp = "C:\\Temp\\VirtualCampus\\IGIS\\Thames.shp" # provide a default value if unspecified          Points_of_interest = arcpy.GetParameterAsText(1)     if Points_of_interest == '#' or not Points_of_interest:         Points_of_interest = "Points of interest" # provide a default value if unspecified          Super_Output_Areas_shp = arcpy.GetParameterAsText(2)     if Super_Output_Areas_shp == '#' or not Super_Output_Areas_shp:         Super_Output_Areas_shp = "C:\\Temp\\VirtualCampus\\IGIS\\Super Output Areas.shp" # provide a default value if unspecified          Thames_Buffer_shp = arcpy.GetParameterAsText(3)     if Thames_Buffer_shp == '#' or not Thames_Buffer_shp:         Thames_Buffer_shp = "C:\\Temp\\VirtualCampus\\IGIS\\Thames_Buffer.shp" # provide a default value if unspecified          superoutputareas_clip_shp = arcpy.GetParameterAsText(4)     if superoutputareas_clip_shp == '#' or not superoutputareas_clip_shp:         superoutputareas_clip_shp = "C:\\Temp\\VirtualCampus\\IGIS\\superoutputareas_clip.shp" # provide a default value if unspecified          # Local variables:     Pointsofinterest_Layer__2_ = Thames_Buffer_shp     poi_cat = Pointsofinterest_Layer__2_     poi_cat_View = poi_cat     Pointsofinterest_Layer = Points_of_interest     pop = superoutputareas_clip_shp     pop_View = pop          # Process: Make Feature Layer     arcpy.MakeFeatureLayer_management(Points_of_interest, Pointsofinterest_Layer, "", "", "FID FID VISIBLE NONE;Shape Shape VISIBLE NONE;CATEGORY CATEGORY VISIBLE NONE;NAME NAME VISIBLE NONE")          # Process: Buffer     arcpy.Buffer_analysis(Thames_shp, Thames_Buffer_shp, "300 Meters", "FULL", "ROUND", "NONE", "")          # Process: Select Layer By Location     arcpy.SelectLayerByLocation_management(Pointsofinterest_Layer, "COMPLETELY_WITHIN", Thames_Buffer_shp, "", "NEW_SELECTION")          # Process: Summary Statistics     arcpy.Statistics_analysis(Pointsofinterest_Layer__2_, poi_cat, "CATEGORY COUNT", "CATEGORY")          # Process: Make Table View     arcpy.MakeTableView_management(poi_cat, poi_cat_View, "", "", "CATEGORY CATEGORY VISIBLE NONE;FREQUENCY FREQUENCY VISIBLE NONE;COUNT_CATEGORY COUNT_CATEGORY VISIBLE NONE")          # Process: Clip     arcpy.Clip_analysis(Super_Output_Areas_shp, Thames_Buffer_shp, superoutputareas_clip_shp, "")          # Process: Summary Statistics (2)     arcpy.Statistics_analysis(superoutputareas_clip_shp, pop, "LSOA_POP SUM", "")          # Process: Make Table View (2)     arcpy.MakeTableView_management(pop, pop_View, "", "", "Rowid Rowid VISIBLE NONE;FID FID VISIBLE NONE;FREQUENCY FREQUENCY VISIBLE NONE;SUM_LSOA_POP SUM_LSOA_POP VISIBLE NONE")~      This is the error i get:                 Traceback (most recent call last):       File "C:\Temp\VirtualCampus\IGIS\IGISprogrammingscript_1.py", line 43, in <module>         arcpy.MakeFeatureLayer_management(Points_of_interest, Pointsofinterest_Layer, "", "", "FID FID VISIBLE NONE;Shape Shape VISIBLE NONE;CATEGORY CATEGORY VISIBLE NONE;NAME NAME VISIBLE NONE")       File "c:\program files\arcgis\desktop10.1\arcpy\arcpy\management.py", line 5774, in MakeFeatureLayer         raise e     ExecuteError: Failed to execute. Parameters are not valid.     ERROR 000733: Output Layer: Same as input C:\Temp\VirtualCampus\IGIS\Super Output Areas.shp     Failed to execute (MakeFeatureLayer).


title: Is there any way to get the call stack in Mapbasic?

text: Is there any (supported or unsupported) way of getting the call stack at a given point in a Mapbasic program?


title: How to speed up queries for raster databases?

text: I have a raster database in postgresql/postgis with these columns: **(ID, rast, data_of_data)**. **'rast'** is the column that has raster files in WKT format. An example query to find the DN value of a point in WGS84 system (30.424, -1.66) and for 2002-01-09 is the following:               SELECT           st_value(rast,(st_GeomFromText('POINT(30.424 -1.66)', 4326))) as val     FROM           my_table     WHERE          date_of_data='2002-01-09'      **Is there a method (eg. spatial index) to speed up those kind of queries?**


title: How to add data from a Google map as layer in ArcMap using C#?

text: A want use C# to write a tool: add data google map as a layer in ArcMap. I am using Arcgis 10.0. Can you help me. Thanks and regards.


title: ArcGIS10.0 (32 bit) waits for splwow64.exe indefinitely and stops responding

text: I am running ArcMap 10.0 (32 bit) on Windows 7 (64 bit), while using some extensions the application waits indefinitely for splwow64.exe and then either crashes or stops responding, shall I try running application on a 32 bit system ?


title: How can I get all the latitude, longitude coordinates/points inside a rectangular polygon drawn in google maps?

text: Google maps lets you obtain a valid address by clicking on houses or buildings one by one, but it would be much easier if you could draw a polygon on the map and get a list of all valid addresses within it. Here's some pseudo code:   1. Draw a rectangular polygon on google maps.   2. Generate a fine grid of lat/lon coordinates/points inside of the polygon (or retrieve if a function/algorithm exists).   3. Reverse geocode all the points inside the polygon that were generated to obtain addresses for each and store the data (could be a db, array, xml, text file etc...).   4. Read the data and filter duplicate addresses and anything that's not a valid address.   5. Store the filtered results. I would appreciate any insight. Thanks!


title: ArcGIS XML Import Failed - how to find out why it failed?

text: I am creating an ESRi XML Workspace document in some of my own C# code. The XML doc that is produced validates against the GdbExchange.xsd schema provided by ESRI. However, when I use ArcCatalog to import the XML document into a newly created Personal Geodatabase (.mdb) I get the following error:               Import xml data failed     FDO error: -2147186428 [XML parse error: Unspecified error      [Line: 415, column: 23].]      The line and column indicated above are:                   ...         </DatasetDefinitions>         <Metadata xsi:type="esri:XmlPropertySet">             <XmlDoc />         </Metadata>     </WorkspaceDefinition> <!-- THIS LINE -->     <WorkspaceData xsi:type="esri:WorkspaceData">         <DatasetData xsi:type="esri:TableData">             <DatasetName>PointInspections</DatasetName>             ...      (Column 23 is after the closing angle bracket, ">", an the line marked "< !-- THIS LINE -->") Given that the XML document is valid, how can I get more useful information from ArcGIS? (There is a somewhat similar question asked here: Importing XML file in ArcGIS failed, but that didn't have a resolution that worked for me...)


title: Need advice on how to use Select Junction by Edge Count at pythot script?

text: Is There a way to use the tool Select Junction by Edge Count (selecting all junctions (points) in a geometric network by the number of edges to which they connect) at Python script ? or can i use select by location with some count function ? or i have to use spatial join count results ?


title: Feature.geometry.transform not working

text:
// World Geodetic System 1984 projection         var WGS84 = new OpenLayers.Projection("EPSG:4326");         // WGS84 Google Mercator projection         var WGS84_google_mercator = new OpenLayers.Projection("EPSG:900913");                   var options = {             controls: [                 new OpenLayers.Control.LayerSwitcher(),                      new OpenLayers.Control.Zoom({                     zoomInId: "plus",                     zoomOutId: "minus"                 }),                 new OpenLayers.Control.Navigation(),                  new OpenLayers.Control.MousePosition(             { div: document.getElementById("coordinates") })             ],             projection: new OpenLayers.Projection("EPSG:900913"),             displayProjection: new OpenLayers.Projection("EPSG:4326"),             units: "m",             numZoomLevels: 20,             maxResolution: 156543.0339,             maxExtent: new OpenLayers.Bounds(-20037508, -20037508,                                    20037508, 20037508.34)         };                             map = new OpenLayers.Map('map', options);         // BASELAYER                   var google_maps = new OpenLayers.Layer.Google(             "Google Maps", {                 numZoomLevels: 20             }         );              var google_satellite = new OpenLayers.Layer.Google(             "Google Satellite", {                 type: google.maps.MapTypeId.SATELLITE,                 numZoomLevels: 20                  }         );         // OVERLAYS         tiled = new OpenLayers.Layer.WMS(         "Kaski",         wmsUrl,         {             LAYERS: 'KaskiGrup',             transparent: true,             visibility: true,                  },         {             transitionEffect: 'none',             isBaseLayer: false,                   }     );              // SAVE STRATEGY OLUŞTUR.         saveStrategy = new OpenLayers.Strategy.Save();         saveStrategy.events.register("success", '', showSuccessMsg);         saveStrategy.events.register("failure", '', showFailureMsg);                             // NOKTASAL LAYERS         wfs_layer = new OpenLayers.Layer.Vector("Point",             {                 /*style: {                     externalGraphic: 'http://www.openlayers.org/dev/img/marker.png',                     graphicWidth: 21,                     graphicHeight: 25,                     graphicYOffset: -24                 },*/                 styleMap: new OpenLayers.StyleMap({                     "default": new OpenLayers.Style({                         pointRadius: 10, // sized according to type attribute                         fillColor: "#00FFFF",                         strokeColor: "#33CCCC",                         strokeWidth: 3,                         graphicZIndex: 1                     }),                 }),                      strategies: [new OpenLayers.Strategy.BBOX(), saveStrategy],                      protocol: new OpenLayers.Protocol.WFS({                     version: version,                     // loading data through localhost url path                     url: wfsUrl,                     featureNS: featureNS,                     visibility: true,                          // layer name                     featureType: featureType,                     // geometry column name                     srsName: "EPSG:900913",                     geometryName: geometryName,                     schema: schema,                      })             });      My map,options and layers. Everythings is fine and work. My layers is at correct location. Also when i click on my map, automaticly openning popup. Thats code               function onFeatureSelect(event) {             var feature = event.feature;             selectedfeatureval = feature;             var mapsinirlaridisinda = 0;                            $.mobile.showPageLoadingMsg();                       function mahalleBul(request) {                      var geojson_format = new OpenLayers.Format.GeoJSON({});                 var getMahalle = geojson_format.read(request.responseText);                 if (getMahalle.length > 0) {                     mahalleadi = getMahalle[0].attributes.ADI;                     mahallegisID = getMahalle[0].attributes.GIS_ID;                 }                      function ilceBul(request) {                          var geojson_format = new OpenLayers.Format.GeoJSON({});                     var getilce = geojson_format.read(request.responseText);                     if (getilce.length > 0) {                         ilceadi = getilce[0].attributes.ADI;                         ilcegisID = getilce[0].attributes.GIS_ID;                     } else { mapsinirlaridisinda = 1; }                               function handler(request) {                              var geojson_format = new OpenLayers.Format.GeoJSON({});                         var getbina = geojson_format.read(request.responseText);                         if (getbina.length > 0) {                             binaid = getbina[0].attributes.ID;                             binismi = getbina[0].attributes.ADI;                              }                                   //Binaya Tıklandıysa                         if (binaid != null) {                                  var numfeature;                                  function numataj(request) {                                      var geojson_format = new OpenLayers.Format.GeoJSON({});                                 numfeature = geojson_format.read(request.responseText);                                                if (typeof (numfeature[0]) == "undefined") {                                          $.mobile.hidePageLoadingMsg();                                     fulladres = mahalleadi + " MAHALLESİ " + ilceadi;                                               var content = '
' + binismi + '

'                                         + '
' + mahalleadi +' MAHALLESİ ' +ilceadi+ '

';                                     content = content + '

'                                         +'BildirimVazgeç

';                                          if (content.search("


title: Plugin QGIS error FPDF' object has no attribute 'y'

text: I try to make save as button in plugin QGIS Dufour 2.01 and want to save it as pdf file when i try save my file ( open file destination and then click save ) i got a error message This is the code               def _save(self, simpan):             import fpdf             # Portrait, millimeter units, A4 page size                  pdf=fpdf.FPDF("P", "mm", "A4")             # Set font: Times, normal, size 10             pdf.set_font('Times','', 12)             # Layout cell: 0 x 5 mm, text, no border, Left             pdf.cell(0,5,'Input 1 : ' +  self.ui.lineInput1.text(),border=0,align="L")             pdf.cell(0,5,'Input 2 : ' + self.ui.lineInput2.text(), border=0,align="L")             pdf.cell(0,5,'Recomendation : ' + self.ui.textRec2.toPlainText(), border=0, align="L")             pdf.cell(0,5,'Data 1 :' +  self.ui.lineCond1.text(), border=0, align="L" )             pdf.cell(0,5,'Data 2 :' + self.ui.lineCond2.text(), border=0, align="L" )             pdf.output( simpan+'.pdf','F')      Error message i get               File "C:\PROGRA~1\QGISDU~1\apps\Python27\lib\site-packages\fpdf\fpdf.py", line 615, in cell         if(self.y+h>self.page_break_trigger and not self.in_footer and self.accept_page_break()):     AttributeError: 'FPDF' object has no attribute 'y'     Python version:     2.7.4 (default, Apr  6 2013, 19:54:46) [MSC v.1500 32 bit (Intel)      I have python33 installed but its like QGIS Dufour use python27 so i install it in C:\QGIS Dufour\apps\Python27\lib\site-packages\ with .msi installer for windows I use seven 32bit


title: How to zoom to layer using ArcGIS Javascript API

text: I tried to get the function "Zoom to Layer" done using the following code snippet. But it does not change anything. The layer is a feature service layer, and I can see the full extent object of the layer which is not null. Any idea? Thanks.                                          tree.on("dblClick", function (item, node, event) {                                var layerId=item.id[0];                                  var layer = map.getLayer(layerId);                                        if (!layer)                                   {                                     return;                                   }                                   else                                   {                                     var lyrExt=layer.fullExtent;                                     var gsvc=  esriConfig.defaults.geometryService;                                     var outSR = new esri.SpatialReference({ wkid: 102100 });                                     gsvc.project([lyrExt], outSR, function(newExt) {                                     map.setExtent(newExt[0],true);                                   }                                    });


title: Best way to store polygon geospatial data in postgreSQL?

text: I have to store longitude and latitude coordinates for several places in a PostgreSQL db table. Each row in this table has a field with a set of pairs {latitude1, longitude1, ....latitudeN, longitudeN} wich represent an enclosing polygon. I stumbled upong PostGIS, and I'm wondering if it is necessary to use Polygon datatype for this particular field in my table, like depicted in this example I'm asking because the only thing I'm gonna do with those polygons is drawing them afterwards in a Google Maps, I don't really think I'm gonna do any reasoning on them (yet)...thus, would be PostGIS overkill for this? Besides all the time I have to spend in installing and learning PostGIS... Any toughts? Thanks!


title: What is the most annoying fact about the GIS market/industry?

text: What are the most annoying facts about our industry/market? What makes you angry? Is there anything that you can't stand it?


title: Which ESRI defined geographic transformation do I need to use to convert my " NAD83 Northern Virginia" data into WGS 1984 Web Mercator?

text: My title pretty much tells it all.. I have a number of datasets in which I want to batch project from NAD83 Northern Virginia Datum into the correct Web Mercator projection.. When I do a straight re-projection I get an error message.. this tells me that I need to use a transformation of some kind.. I think.. So, here in lies my problem.. as there is not a web mercator projection transformation listed in the geographic transformation I was thinking I would just use a generic one like... NAD_1983_To_WGS_1984_2 _(there are 8 different ones)_ WGS_1984_Major_Auxillary_Sphere_To_WGS_1984 But I am not sure that is correct.. Which is the correct "geographic transformation"?


title: Switch coordinates in linestring

text: I made a silly mistake: when creating linestrings in my PostGis database i inverted the lon and lat values. Which isn't a problem for postgis itself, but for my visualization and intersections. ;) Of course i only recorgnized this after importing lots of stuff. Instead of simply re-running the import i'm curios if there is a way i could swap the X/Y-Coordinates of every point of a linestring? e.g.:               LINESTRING(lat1 lon1, lat2 lon2, lat3 lon3); // what i got     SomeMagicMethod();     LINESTRING(lon1 lat1, lon2 lat2, lon3 lat3); // what i want


title: call child div function from parent (not using window.opener)

text: This is a a Javascript/JQuery question that is directly related to the ArcGIS Javascript API... I have a map div created using a 'data' tag in HTML, which becomes visible on a button click. (I'll call the map the 'child' and main page the 'parent' although I am not sure that is technically correct in this case?) I would like to force a resize of the map when a 'fullscreen' button is clicked in the parent page. Generally the map resizes on its own but sometimes it does not... How can I do this? I did not create a new window with 'window.open' ... My HTML for the child window is:               <div style='width:100%;  height:100%'><object type="text/html" style='width:100%;  height:100%' data='mapPage.html'>      ..and the JS for the fullscreen button :                   mapFullscreen.click(function () {         window.child.map.refresh(); // this doesnt seem to do anything              $(mapContainer).animate({             "height": $(window).outerHeight(),             "width": $(window).outerWidth(),            "top": 0,            "left": 0             }, 300, function(){         $(".map-fullscreen").fadeOut();       });     });      Thanks!


title: How to digitize from a map with additional detail views?

text: I have a map with a river that includes islands within it. Some islands have a blown up version (detail views) outside of the river showing land ownership details. I would like to make polygons (shapefiles) of the large islands on the outside of the river and then copy these polygons and shrink them to fit within the lines of their matching islands within the river and then safe the result. I suspect that the copy and paste part is possible but can the "shrink or resize" part be done? I want this so that I can do distance queries along the river. Is this possible and how would I do it? Thanks


title: Is there a way to control label density in mapserver?

text: I'm reading a US County Tiger Shapefile and everything is plotting correctly, but I want to slim down the number of labels shown at scales larger than the state level because it begins to overwhelm the map. Is there a way to modify this to do that, but leave the zoomed in look and feel alone?               LAYER     NAME 'us_counties'     TYPE POLYGON     DUMP true     TEMPLATE fooOnlyForWMSGetFeatureInfo       DATA './shapefiles/uscounties_20130327.shp'     METADATA       'ows_title' 'counties_us'     END         PROCESSING "LABEL_NO_CLIP=True"          PROJECTION       'proj=longlat'       'ellps=GRS80'       'datum=NAD83'       'no_defs'     END          LABELITEM 'NAME'     CLASS       NAME 'uscounties_20130327'        STYLE         WIDTH 0.91          OUTLINECOLOR 0 0 0         #COLOR 108 180 98       END       LABEL        FONT verabd       TYPE truetype       SIZE 12       COLOR 0 0 0       OUTLINECOLOR 255 255 255       ANGLE 0       POSITION cc       FORCE true       ANTIALIAS true       PARTIALS false     END      END       END


title: Building Arc add-ins without needing to restart Arc

text: One of the things that annoys me the most with trying to use the Arc add-in wizard is that I need to restart Arc every time I want to test something. This is especially annoying if I've got syntax errors or something little that I need to troubleshoot. What approach would you recommend to negate the need to start Arc repeatedly? It has been noted in comments that this question has already been asked and vaguely answered; however it's totally unclear to me how you integrate this into your script. Can someone give me an example? I've tried to include the following script as noted in the other question linked here.               import os          esri = "C:\Program Files\Common Files\ArcGIS\bin\ESRIRegAddIn.exe"     cwd = os.getcwd()          mapdoc = "C:\Data\Shaun Pilot Area\PilotAreaMXD.mxd"          try: os.system("TASKKILL /F /IM ArcMap.exe")     except: pass          os.system(os.path.join(cwd, "makeaddin.py"))          os.system('"{0}" {1} /s'.format(esri, os.path.split(cwd)[-1] + ".esriaddin")          os.system(mapdoc)      Here's the output I get               ERROR: The process "ArcMap.exe" not found.     'C:\Data\Shaun' is not recognized as an internal or external command, operable program or batch file.     'C:\Program' is not recognized as an internal or external command, operable program or batch file.     'C:\Data\Shaun' is not recognized as an internal or external command, operable program or batch file.


title: How to copy features from a vector layer to another vector layer with only correspondent attributes in QGIS

text: I have two vector layers: A first one with all cadastral parcels of a region. This layer has several attributes, it's only the geometry and 3 of the attributes I'm intrested in. My second layer is a selection of parcels from the same region. This layer has 7 attributes. The first 3 attributes are identical to 3 attributes from the first layer with all cadastral parcels. The 4 other attributes are filled in afterwards. How can I copy a parcel from the first layer to the second one, only with the geometry and matching attributes in the right field of the attribute table?


title: 'D8 Slope Grid' in Sextante-QGIS is not working

text: I can't make 'D8 Slope Grid' work in QGIS-Sextante. 'Pit Filled' and 'D8 Flow Direction Grid' output work okay.   I used the function: Geoalgorithms > Uncategorized > Basic Grid Analysis tools > D8 Flow Directions inside Sextante Box. Then I filled out the required parameters.   What have I missed? Any help is appreciated.   Thanks. \- Al


title: Is it possible to get WMS feature info for the WMS layer in the QGIS web client?

text: When I load a wms layer in qgis, it's possible to get feature information. When I open the web-gis-client on my server, it's possible to get feature information for all layers but not the wms layer. Is it possible to get wms feature info for the wms layer in the client? Maybe I'm doing something wrong or isn't it possible at all. Please help.....


title: Plotting property surveys using AutoCad MAP2013 - Draw arc problem

text: This is likely a simple question that I just can't seem to answer. I am plotting some surveys in AutoCad Map3D but I simply can't figure out how to draw an arc using survey coordinates. The large majority of my surveys are straight lines which I input as [@distance, bearing] which works fine. How do I use arc coordinates (Arc,Radius,Bearing,Distance) supplied by the surveyor in my plotting? * * * I should mention that I am using shapefiles and creating new polygon features to store my boundaries.


title: Assign Z-values of river stream segments

text: I have used Interpolate shape tool in 3D Analyst (ArcGIS 9.2) in order to obtain Z-values of each segment along the stream lines. As input features I use Riverstream.shp file and as raster surface the DEM. With right click on interpolated riverstream.shp after opening the attribute table, I created new column for X and Y coordinates. With right click on each column(for X & Y) using Calculate geometry I obtained X,Y values. I need also to obtain elevation values (Z-values ) for stream segment as well as length of each segment in meters (with Calculate geometry there is not option to obtain Z values). I tried with Field calculator but without success. Can anybody help me : How can I assign Z-values of stream segments. For your help thank you in advance. Subija.


title: Spherical law of cosine what does the 'a' represent?

text: I want to use the spherical law of cosine formula, θ = atan2( sin Δλ ⋅ cos φ2 , cos φ1 ⋅ sin φ2 − sin φ1 ⋅ cos φ2 ⋅ cos Δλ ), to calculate distances between coordinates, but, I cannot for the life of me remember what the 'a' stands for. Can some please help me out?


title: ADC Maps Grid to Lat/Long Coords Conversion

text: While the ADC site apparently has a paid conversion function, it's interactive, and I'm looking for some way to do so on a batch basis, possibly with an API. Any hope?


title: Is there a way to run spatial join within tolerance?

text: I have some data that should be spatially joined, but from small precision differences (less than 1 in) they aren't joined. I've found no way to add tolerance to qgis spatial join operation. Is there a way to config a tolerance value, or a plugin to correct small distances?


title: How to load WMT_MS_Capabilities using Javascript

text: I would like to load and parse an `WMT_MS_Capabilities` using Javascript. I have unsuccessfully tried severaly ways. I am trying with this file but I always get this error:               uncaught exception: [Exception... "Component returned failure code: 0x80004005 (NS_ERROR_FAILURE) [nsIXMLHttpRequest.send]" nsresult: "0x80004005 (NS_ERROR_FAILURE)" location: "JS frame :: file:///home/frowing/Projects/XMLParsing/testXML.js :: <TOP_LEVEL> :: line 10" data: no]      I hope someone can throw some light on this because I have no idea where to go next. Thanks!


title: SpatialRefSys as foreign key

text: Using the working stuff reported here, I've written the following django model:               from django.contrib.gis.models import SpatialRefSys          class LandCover(models.Model):         shp_name  = models.CharField(max_length=100)         geom_srid = models.ForeignKey(SpatialRefSys)      But syncdb gives me this error :               [user@host]$ ./manage.py syncdb     CommandError: One or more models did not validate:     myapp.landcover: 'geom_srid' has a relation with model <class 'django.contrib.gis.db.backends.postgis.models.SpatialRefSys'>, which has either not been installed or is abstract.      Is there a way to solve this?


title: New geographical point won't show on map

text: I did recently ask a question about how to add a new point to a PostGIS table. Now when I've got the new point in the table it should be visible on my map, but it isn't. I'm using Mapnik together with TileStache. The cache for TileStache is temporary turned off, and I've restarted PostgreSQL. It's still not there. The query I used to add the point was the following:               INSERT INTO osm_transport_points(osm_id, name, type, geometry)     VALUES(1, 'name', 'station', ST_SetSRID(ST_MakePoint(lat, lng), 900913));      Is there something I've missed to flush, or am I doing something else wrong?


title: QGIS Field Calculator Results Removed When Saving Edits

text: Mac OSX 10.6.8 QGIS 1.8.0 : Using the field calculator from the attribute table, the $area of each feature is shown nicely in the table until I try to save it, then all the values are replaced by NULL so are no longer available for use. Furthermore, if the column is left there and the field calculator asked to update it, all values remain NULL. Any fix?


title: Automate export of ArcMap legends

text: I've been building legends for my javascript websites that use 9.3.1 services. I know there's the ESRI legend tool for older services but I don't feel that it does a good job of producing legend images for raster layers, especially continuous raster layers. It also drops things like the symbol headings. Ideally I'd like to be able to automatically mimick the look of the ArcMap TOC. I've been trying to build some ArcObjects code to create a Legend to attach to the Page Layout and then export that to an image, but I can't seem to be able to find the code necessary to export the image after I build the legend. the following code is the basic crux of what I've tried so far:                   void OpenMXDViaMapDocument(string path)         {                  IMapDocument pMapDocument = new MapDocumentClass();                              if (pMapDocument.get_IsMapDocument(path))             {                 pMapDocument.Open(path, null);                 IMap pMap;                      IPageLayout layout = pMapDocument.PageLayout;                           for (int i = 0; i <= pMapDocument.MapCount - 1; i++)                 {                     pMap = pMapDocument.get_Map(i);                     UID pId = new UID();                     pId.Value ="esriCarto.Legend";                     IMapSurround surround = pMap.CreateMapSurround(pId, null);                     ILegend iLegend = (ILegend)surround;                     int p = iLegend.ItemCount;                          iLegend.AutoVisibility = false;                          System.Drawing.Bitmap myBitmap = new Bitmap(500,500);                     IntPtr ptr = myBitmap.GetHbitmap();                          ESRI.ArcGIS.Display.SimpleDisplayClass sdc = new ESRI.ArcGIS.Display.SimpleDisplayClass();                          for (int y = 0; y < p; y++)                     {                         ILegendItem iLegendItem = iLegend.get_Item(y);                         LegendFormat legendFormat = new LegendFormat();                         legendFormat.ShowTitle = true;                              iLegendItem.CreateGraphics(sdc,legendFormat);                         IEnumElement enumEles = iLegendItem.Graphics;                         enumEles.Reset();                         IGraphicElement ige = (IGraphicElement)enumEles.Next();                         while (ige != null)                         {                             FileStreamClass fs_class = new FileStreamClass();                             richTextBox1.Text = richTextBox1.Text + "#"+y+" - Graphic Element Type " + identifyGraphicElementType(ige) +"\n";                             if (identifyGraphicElementType(ige).Equals("BmpPictureElement"))                             {                 //All of them are apparently BmpPictureElements                                 BmpPictureElement bmpE;                                 //Java version has bmpE.save(IStream) ?                                  fs_class.SaveToFile("C:\\utils\\newImages\\{0}_{1}.bmp");                             }                             ige = (IGraphicElement)enumEles.Next();                         }                                             }                      }                  }         }      I would also like to generate an image for each map layer showing the title, symbol heading if any and the symbol. How can I do this?


title: Vector Layers are not lining up and are scaled differently

text: Hey I am a first time user of GIS and I am working for a transportation map of North Carolina for a project at my work. So what I have been doing is first placing the Tiger County Boundary by using the add vector layer tool and selecting the directory setting and then select tiger. I then do the same thing for a layer that contains all of the roads. They do not line up and the road layer is significantly larger than the county boundaries. I tried to mess with the CRS, but it did not seem to change anything. I am just taking a guess, but I am assuming that it is a scaling issue? Anyway please help


title: With multi-line labels, how to show the 2nd line only for labels that need it?

text: Mac OSX 10.6.8 QGIS 1.8.0 : I've two columns in the attribute table, name, and name2 for those labels that need to be 2 lines; multiline labelling is checked. But when I use "name"||'\n'||"name2" I only get labels when there is a printable entry for both name and name2. Can you spot what I'm missing to show the 2nd line only for labels that need it?


title: How to manipulate SDF data via java/jsp (MapGuide open source Web API)?

text: I am relatively new to MapGuide open source technology. I would like to parse and modify data that is in SDF format so that my map displays this updated feature data (dynamically). How can i do so? I have MapGuide installed and configured on Linux. I went through the MapGuide Web API documentation but there were no examples on how to manipulate SDF data. Currently i am generating the SDF data in a different environment (in a Windows OS domain) by using scripts to query a spatial database and the query generates SHP files. I then transform those SHP files into SDF files (using scripts) before transferring them to the MapGuide server repository in the Linux OS environment. On the Windows OS end, I am able to edit the SHP files (in particular the dbf file just using a spread sheet) to change the values for the particular feature source. However this process is static and so i was hoping to be able to do so programmatically (using java / jsp) by parsing the SDF files for the particular feature source in the Linux domain. Thanks


title: How to enable label orientation?

text: Mac OSX 10.6.8 QGIS 1.8.0 : I'm mapping an area that's a rectangle slanted at 26 degrees, and neighbourhood maps are always printed so the rectangle is horizontal. I'm trying to make the labels readable when printed that way. In advanced settings the placement 'around centroid' shows a rotation option, but it's greyed out and unresponsive; none of the other options allow specifying a label orientation. The move/rotate tools of ABC are also greyed out and unresponsive. (The 'free' setting puts labels everywhichway and is very awkward to read.) I've added x, y and rotation attribute fields (type real) and enabled editing. Might there be a way of enabling the rotation field that I haven't spotted? Or, something else that will have a similar effect?


title: MapInfo very slow to start-up

text: MapInfo v10.5 is used across my organsiation by dozens of users as a thick client installation. A small group of users have started to have an issue whereby the time taken for a new session of MapInfo to open can last up to 30 minutes. I have logged onto each users’ machine and experienced the same problem as them so this doesn’t seem to be profile related. When checking the Task Manager during MapInfo startup, the MapInfo process doesn’t seem to use more than 1% of CPU and struggles to use above 26mb memory. When I do this on my own machine the CPU usage goes up to 13% and the memory usage quickly rises to 175mb and fully loads in less than a minute. Anyway, just wondering if anyone else has experienced a similar issue with MapInfo. Thanks, Robin


title: Query Feature using OpenLayers + MapServer

text: I'm trying to query a feature in OpenLayers, but it's returning just an image highlighted that is not possible to make zoom, pan and drag. My intention is to query a feature on the map like street name, city and country. Is there any function in OpenLayers? I'm using this code:               var query = new OpenLayers.Layer.WMS( "Utbs",          "http://localhost/cgi-bin/wms.exe",          {         srs: 'EPSG:4326',         width: '800',         styles: '',         height: '600',         layers: 'layer_name',            mode: 'itemquerymap',         qstring: "(value = '03')",         qlayer: 'layer_name',         mapxy: 'shape'           },         {         singleTile: true             }     );          map.addLayers([utbs_query]);     map.zoomToMaxExtent();


title: about updating WFS with popup

text: I want to ask you something about this How to collect data via WFS and Popup I started to build my own WebGis using openlayers and geoserver, but I'm newbie here I want to add form to update the attribute feature data from the client can you explain me more about the php code snipet inside those codes? this is the code below               $res = pg_query ("UPDATE $tabelle SET entryname = '".pg_escape_string ($entryname)."' WHERE gid = '$gid'");      Oh yeah, by the way, I've tried the WMS one with this tutorial : http://anitagraser.com/2011/03/17/picking-a-features-attribute-value-from-a- wms-layer-with-openlayers/ but it's not working


title: How to identify a FRONT LAND?

text: I work at a town hall. I need to identify the front of the land. Here in my country we call it of "TESTADA" or "FRENTE DE LOTE". I did not find the exact term in English. THE FRONT OF A LAND is a line near the entrance to the ground. That is, is the entry of a house. It is where the entrance gate, with the port number. To exemplify following image below. **Map before:** ![before](http://i.stack.imgur.com/qoLCz.jpg) **Map later with the front lines of terrain - As should be:** ![after](http://i.stack.imgur.com/2nOgJ.jpg) A LOT has only one FRONT OF LAND. In corners it's main street. There are old buildings that will stay as they are because it is historic. Currently we do these lines of FRONT LAND manually. We would like to automate part of it. That could be a Model Builder, Python, etc.. **In short:**   * In: LAYER OF LAND, LAYER STREET, BUILDING LAYER.   * Output: LAYER with FRONTS OF LAND LINES. **Another example of front lot:** The term FRONT OF LAND refers to front of property.The input of the house is his front. ![enter image description here](http://i.stack.imgur.com/hTmpC.jpg)


title: How to apply DEM elevation adjustment parameters?

text: I have a nice DEM, created from UAV data, which I want to process further in QGIS, my platform of first choice. In order to represent the evelation of each pixel correctly, I have to adjust the 'elevation adjustment parameters' like Z-scale and Z-offset. In Globalmapper, I was able to do it. Also in the manual of ArcGIS, it is discussed (http://resources.arcgis.com/en/help/main/10.1/index.html#//009t0000003t000000). I do, however, not encounter the way how to do it in QGIS. Does have anyone a clue or advice on how to do that in QGIS? Many thanks, Arnoud Cuppens


title: ArcMap Arcpy scripts: How to handle user password input?

text: I'm currently trying to write a script that has a username & password input which then returns some data if the user has appropriate access. The script works just fine when I hardcode my own credentials in. The problem I have is that I don't know how to mask password input using the "Add Script" wizard. It appears there is an option to accept an Encrypted String as a parameter, which appears to be what I want, but I can't find any documentation on how to actually utilize the GPEncryptedString class within Python. When I try to extract the parameter with arcpy.getParameterAsText(), it just returns a bunch of asterisks as a string value which obviously does me no good. I even tried just popping up my own username/password dialog, but scripts have to be ran with the "Always in Foreground" option for the arcpy.mapping.MapDocument("CURRENT") command to work properly, which constantly steals focus from my dialog. This seems like a really simple problem, but I'm just not finding any documentation anywhere on this issue and I've now spent more time trying to figure out a way to mask password input than I spent on the actual script itself.


title: MODE_SNAPSHOT or MODE_ONDEMAND? Issues with both, which one to focus on?

text: Our web app uses the ESRI JS API. With SNAPSHOT, our feature layers take a long time to show up on the map, especially after using setDefinitionExpression to restrict the feature layer to a date range. If we change it to ONDEMAND, our feature layers show up very quickly, but not reliably. It often takes wiggling the mouse (panning a small amount of pixels) to get the feature layers to show up. Does this sound familiar to anyone? Which of the two modes would be a better starting point to get things working right? SNAPSHOT got much quicker after we turned off TimeInfo in our feature layers, but it's still noticeably slower than ONDEMAND.


title: Attribute Query for First/Last record and every nth inbetween

text: I have a feature class with over 5000 points. There are too many to label so I would like to be able to grab the first, the last and every 30th or 50th in between.


title: How to use QgsOverlayAnalyzer class in pyQGIS?

text: I am very new to Python scripting in QuantumGIS, so bear with me please. I am trying to use the QgsOverlayAnalyzer class in QuantumGIS (http://www.qgis.org/api/classQgsOverlayAnalyzer.html#details) to automate the intersection between different layers (I have plenty of intersections to do). However, I can't get it to work in the python console.. assume I want to intersect layer1 with layer2, and write the results in `C:\shapefile.shp`, I am trying something like this with no success:               QgsOverlayAnalyzer.intesection(layer1, layer2, "C:\shapefile.shp", false, 0)      The error I make must be quite serious, as the message I get is:               NameError: name 'QgsOverlayAnalyzer' is not defined      Any help on this is very much appreciated. Dom


title: Looking for open source Spatial file server that ArcGIS 10.1 will connect and work with?

text: I'm looking for an open source spatial file server to house and serve my GIS data. I'm currently using my Ubuntu Server as a ssh file server and have to download the data to the local computer, manipulate and then re-upload to the server. I would prefer to use ArcCatalog to connect to the server and manipulate the data in-situ. Are there any open source Spatial servers that Arc will communicate with? I'm not very interested in WMS or other web based displays.


title: How can I convert lots of ASCII files (.xyz format) into shapefiles by iterator in the Modelbuilder of Arcgis 10?

text: I have a plethora of xyz format files and I want to convert them into shapefiles separately and it is really time-consuming to do this conversion one by one. I tried to create a model in ModelBuilder to read them from an input folder and does the conversion into shapefiles and then put them into the output folder with the same input name or use a wildcard to choose some of them to convert. * * * ![enter image description here](http://i.stack.imgur.com/j5Zy4.jpg) However, the files that I'm trying to convert are XYZ foramt (e.g. file.xyz) ASCII files. Anyway, I used the following model to convert the files in the folder and what it really does is just converting the first file and then stops. I also ticked on the Recursive but it only converts the first file and stops. Can anyone help with that?


title: Moisture Interpretation from Landsat Images

text: I am working on a **"land suitability for crop"** project. Can anyone let me know how to get any type of information which can be related to moisture index of soil? I have Landsat 5 data with band 3 and band 5 and working on ArcGIS 9.3. Any type of suggestion is most welcomed.


title: How to convert WGS84 to UTM when the study area covers 4 UTM Zones: 43- 46?

text: I'm working in ARCGIS 10.2, my area of interest falls in 4 UTM zones 43,44,45,46 equally distributed so i can't use the option of majority of area's UTM zone. Since i need to carryout varies analysis which involving distance and area i need it to be converted also I need to apply same on my SRTM DEM of the study.


title: Openlayers WMS layer not diplaying in IE windows mobile 6 emulator

text: I am new to openlayers. I tried this example. I have published world shapefile to tiger workspace as 'country' layer to geoserver and want to access with WMS. The world map is getting displayed in IE on Windows XP but not in IE of Windows 6 mobile emulator. I am able to access geoserver's admin page on my PC from the mobile emulator so network connection may not be issue. Kindly suggest anything I am missing.                           <title>OpenLayers Example</title>


title: Developer-friendly Embedded systems GIS library?

text: Having seen the answer to the smallest GIS library for embedded devices, I'm asking a related question. I want to do simple point-in-polygon and other spatial relationship tests against a shapefile. But I also want to invoke REST calls and so I just don't want to be using low-level C/C++ style languages because...well, they're painful. Also, the device can run Windows Embedded or Linux Embedded, but preferably the solution will work on either.


title: How to import shapefiles into MS SQL 2008 and then view that data using QGIS?

text: I'm a QGIS users but not a database person and I need some help on the DB side. I'm currently in the need to import several shape files into MSSQL 2008 and I have no clue how to do that and my IT guy either. Can anyone explain to me the components needed,and the step by step how to proceed, or point me towards an online tutorial for unexperienced DB person like me please. Furthermore, lets assume that all my shape files are imported into MSSQL 2008 that is installed on an internal secured server, Do I need to install another component to have access to my shape files from outside my organization (ie: from a remote location via the web ) ? or will I be able to access my data easily using the QGIS Browser ?. If yes, what do I need to do (or tell my IT guy to do !!) ? As you can tell, I'm a real rookie when it comes to DB servers...I just hope soemone can point me towards the right direction so I can help my IT guy in the process. Thanks for all your help. Yann


title: How to create text label of different scales?

text: I would like to create a different type of map in ArcMAp.I would like to present the frequency of some words. So the dimension of the word (maybe even the color) would symbolize its frequency. I started by created a polygon for each word and in attribute table I created two columns "word" and "frequency". So how can I use the words from the colomn "word" to create labels and give them size according to the colomun "frequency". The result I want to obtain is something like this: http://rising.globalvoicesonline.org/chennai/files/2011/06/figure-1-wordle.png So if anyone has an idea how to create something like this in ArcGIS, I would be very grateful. Thank you very much.


title: Model builder - Nesting models to extract raster values to points

text: ![Random Point model](http://i.stack.imgur.com/sVXuF.jpg) Hello! Hopefully a quick question. I have a simple model below that iterates through a folder full of rasters, taking each in turn to extract data values based on a single shapefile of Random point locations. The rest works over the points table to add a field named after the raster, calculates it from the RASTERVALU field, then deletes the RASTERVALU field. That file is fed back into the model as the next raster comes through to repeat the process. The end result is a series of shapefiles where the final one has as many fields as raster, named after the raster, and with its values. Now I want to do this for many shapes/feature classes of random points. I have tried to use this is a nested iterator model with no luck. When nested, I can pass the name of each random point shp so that it does not over-write the previous run, but in that case the raster iterator only works for the first raster. I realize this might be asking a lot of the model builder, but maybe someone has a better idea fro nesting or a different logic approach. Thanks!


title: QGIS crashes by opening project

text: Suddenly I got a problem opening the project I was working on since a long time. When I open the project, QGIS crashes. I think the problem might be in one layer, but as I cannot open the project I neither can delete this layer in QGIS desktop. Is there any option to do this another way? Thanks!!


title: Do I need to install ArcSDE?

text: I've successfully installed ArcSDE and ArcGIS for Server for 10.1. I'm running SQL Server 2012 on a Windows 2012 Server x64 machine. I'm able to direct connect to the SDE geodatabase and version the data. So I assume this means that ArcSDE was installed and working correctly? But I'm a little confused about the ArcSDE services and what exactly they are used for. I feel like I may not be using ArcSDE for what it was made for without using these services. Is this correct?


title: how to install networkx python in QGIS

text: How to install and how to compile networkx python on QGIS? I always get an error that `no module named networkx` in the console of QGIS I copied the file networkx in `C:\Program Files\Quantum GIS Lisboa\apps\Python27\Lib\site-packages` and there are not change.


title: Is there a public CanTopo Web Map Tile Service (WMS / WMTS) available?

text: Does anyone know of a public CanTopo WMS / WMTS available anywhere? I'm assuming CanTopo is more modern / up-to-date than Toporama? Alternatively, what would be your suggestions for the most pain free and cheapest version to host such a service myself? The data seems to be freely available here: http://www.nrcan.gc.ca/earth- sciences/geography/topographic-information/data/11042


title: Make Route with Python and network analyst in ArcGIS 10

text: Good Day, I'm having trouble debugging my python code to create a network analyst layer in ArcGIS 10.0. Basically when run from a script tool the code completes successfully without actually doing anything or returning any sort of useful information, and when run from the embedded python window in ArcGIS 10.0 i get the following error. "Runtime error : Failed to execute. Parameters are not valid. ERROR 000192: Invalid value for Time Attribute Failed to execute (MakeVehicleRoutingProblemLayer)." The script is supposed to solve the route and create a layer file for the output Here is my code below:               #Import system modules     import arcpy     from arcpy import env          try:         #Check out the Network Analyst extension license         arcpy.CheckOutExtension("Network")              #Set environment settings         env.workspace = r'C:\Layout3.gdb'         env.overwriteOutput = True              #Set local variables         inNetworkDataset = r'C:\Layout3.gdb\NtsFDS\NtsFDS_ND'         outNALayerName = "TestResult"              impedanceAttribute = "CABLE_COST"         distanceAttribute = "TRAVEL_TIME"         timeUntis = "Minutes"         distanceUntis = "Meters"         inOrders = r'\Layout3.gdb\NtsFDS\Orders'         inDepots = r'C:\Layout3.gdb\NtsFDS\Depots'         inRoutes = r'C:\Layout3.gdb\NtsFDS\Route'         outLayerFile = env.workspace + outNALayerName + ".lyr"              outNALayer = arcpy.na.MakeVehicleRoutingProblemLayer(inNetworkDataset, outNALayerName,                                                              impedanceAttribute,                                                              distanceAttribute)              outNALayer = outNALayer.getOutput(0)                   subLayerNames = arcpy.na.GetNAClassNames(outNALayer)              ordersLayerName = subLayerNames["Orders"]         depotsLayerName = subLayerNames["Depots"]         routesLayerName = subLayerNames["Routes"]              candidateFields = arcpy.ListFields(inOrders)         orderFieldMappings = arcpy.na.NAClassFieldMappings(outNALayer, ordersLayerName,                                                            False, candidateFields)         arcpy.na.AddLocations(outNALayer, ordersLayerName, inOrders, orderFieldMappings,"")                   depotFieldMappings = arcpy.na.NAClassFieldMappings(outNALayer, depotsLayerName)         arcpy.na.AddLocations(outNALayer, depotsLayerName, inDepots, depotFieldMappings, "")              arcpy.na.AddLocations(outNALayer, routesLayerName, inRoutes, "", "")              #Solve the VRP layer         arcpy.na.Solve(outNALayer)              #Save the solved VRP layer as a layer file on disk with relative paths         arcpy.management.SaveToLayerFile(outNALayer,outLayerFile,"RELATIVE")              print "Script completed successfully"          except Exception as e:         # If an error occurred, print line number and error message         import traceback, sys         tb = sys.exc_info()[2]         print "An error occured on line %i" % tb.tb_lineno         print str(e)


title: Some terms related to coordinate

text: I was googling coordinate related information all the day. Now I am confused with some concepts and terms. I hope some one can help me to make it clear.   * **Geographic coordinate system**      Geographic coordinate is used to describe the location of the earth in a three-dimension which will base on the ellipsoid and a datum. And I always thought that the geographic coordinate are always the latitude and longitude. Am I wrong?   * **Projected coordinate system**      It seems that the projected coordinate system will contains three points:     1. The projection itself. It means how the earth will be translate to a two-dimension surface. For example the `Cylindrical` and `Pseudocylindrical`.     2. The geographic coordinate system which this pcs working for.     3. The Mathematics manner to translate the geographic coordinate to projection coordinate. BTW, is projection coordinate always a point inside the the Cartesian coordinate system? And is projection coordinate only fit for one kind of geographic coordinate system?    * **Coordinate reference system**      Now I have no idea about what a coordinate reference system is. Is the Coordinate reference system has anything to do with the projection? BTW,is the geodetic latitude is measured on the surface of the referenced ellipsoid?