title: I want to show the information of several layers of a request made GeoJSON geoserver

text: Always enters the error function and shows me the alert 'what happen?'. This is the piece of code that causes problems:               map.addEventListener('click',onMapClick);     function onMapClick(e)     {         var latlngStr = '('+ e.latlng.lat.toFixed(3) + ", " + e.latlng.lng.toFixed(3)+')';         var BBOX = map.getBounds()._southWest.lng + ", " + map.getBounds()._southWest.lat + ", " + map.getBounds()._northEast.lng + ", " + map.getBounds()._northEast.lat;         var WIDTH = map.getSize().x;         var HEIGHT = map.getSize().y;            var X = map.layerPointToContainerPoint(e.layerPoint).x.toFixed(0);         var Y = map.layerPointToContainerPoint(e.layerPoint).y.toFixed(0);         URL = 'http://192.168.60.55:9090/geoserver/ems/wms?request=GetFeatureInfo&service=WMS&version=1.1.1&info_format=application/json&exceptions=application/vnd.ogc.se_xml&layers=ems:tblsevilla,ems:tblmanz,ems:tblcarre&styles=&srs=EPSG%3A4326&format=image%2Fpng&bbox='+ BBOX + '&width='+ WIDTH +'&height='+ HEIGHT + '&query_layers=ems:tblsevilla,ems:tblmanz,ems:tblcarre&feature_count=20&x='+ X +'&y='+ Y;                   function onEachFeature(feature, layer) {             if (feature.properties && feature.properties.PERIMETRO)             layer.bindPopup("Perimetro: "+feature.properties.PERIMETRO);             }         }              $.ajax({             type: 'POST',             url: URL,             dataType : 'json',             success: function (response) {                 alert("2");                 L.geoJson(response, {                     onEachFeature: onEachFeature                 }).addTo(map);             },             error: function (req, status, error) {                 alert('what happen?');             }         });     }      'm Leaflet and geoserver beginner, any help would be welcome. Thank you very much in advance to all.


title: Is there a simple interface from which I can extract average summer temperature and humidity for a list of latitude and longitudes?

text: I am interested in finding the summer average (either climatalogical mean or for a specific year) relative humidity and temperature for a set of globally distributed latitude and longitude coordinates. What I need are mean values for June, July, and August. I would prefer to use actual observed values, but a reanalysis (such as NARR or CRU-NCEP would suffice). Is there a web-based interface from which I can obtain this information? The closest that I have come to this   * Wolfram Alpha, but they only provide rainfall and temperature, not relative humidity   * worldclim.org looks promising, but requires me to download all of the data before using another program to extract information; it also does not have relative humidity.


title: How to batch a model created with "model builder" in arcmap 10?

text: I created a simple model that selects iteratively elements from a point dataset and calculates a cost distance raster for each of these points. It works fine when I run it for just one dataset, but when I try to batch it get the following error. ERROR 000865: Input raster or feature source data: I_atlanta_OBJECTID does not exist. ERROR 000581: Invalid parameters. Failed to execute (Cost Distance). Please any help will be greatly appreciate. Here a link to a snapshot of the model http://img88.imageshack.us/img88/7432/modelmultidist.png


title: Why are steps out of order when running model in Batch?

text: I have an issue using an ArcMap model in batch mode. The purpose of my model is to obtain the result of an OD cost matrix from a selected network and a set of origins and destinations. My goal is to be able to get these results without having to add the data to the display in ArcMap (time and resource consuming), but most of all, I want to be able to automatically process several OD cost matrix one after another. For these reasons, I rebuilt the OD cost matrix routine using ArcMap's Model Builder. ![enter image description here](http://i.stack.imgur.com/6gAj9.png) ![enter image description here](http://i.stack.imgur.com/MJ9KL.png) If I test the model on only only network, it works great. However, when I try to batch the process where each time I have a different network with different locations, I have a problem since ArcCatalog seems to mix the different steps. For example, if I batch with only two different networks... ![enter image description here](http://i.stack.imgur.com/riG94.png) Instead of having:               1. Create OD cost layer (1)     2. Load Origins 1 in network 1     3. Load Destinations 1 in network 1     4. Solve OD cost matrix (1)     5. Create OD cost layer (2)     6. Load Origins 2 in network 2     7. Load Destinations 2 in network 2     8. Solve OD cost matrix (2)      ArcCatalog does this:               1. Create OD cost layer (1)     2. Create OD cost layer (2)     3. Load origins 1 in network 2     4. Load origins 2 in network 2     5. Load destinations 1 in network 2     6. Load destinations 2 in network 2     7. Solve OD cost matrix (1)     8. Solve OD cost matrix (2)      I thought that when batching a model ArcCatalog would do the first task completely before even looking at the second. Why are the steps of my model mixed up? Thanks for your help!


title: Why does ModelBuilder fail to clear batch parameter inputs and only use first batch line's input paramters?

text: I have built a model to convert coverage files to shapefiles and also feature classes in a geodatabase. The model uses a sub model to run an iterator to loop through an initial input folder (Which is set to level 2-102A, etc) which has multiple sub levels of folders:   1. TRIM   2. 102A, 102B,102C, etc   3. 102001P, 102002P, etc.   4. 102001P_hydrology, 102001P_contours, etc.   5. With files of various extensions ( point, tic, arc). Then I switch to the main model and use the select data tool to group the various file extensions and collect values tools. At this point I merge the values and create two paths, one which uses the feature class to shp. and the other which uses the feature class to GDB. I have set the initial folder, wildcard and output Shp folders and GDBs as parameters and when I run the model 1 input at a time it works fine however there are 23 intial input folders each with 7 subfolders so the use of the batch option would really save time however I am running in to problems.   1. The model appears to be retaining some memory of past operations or trails, meaning I tried to batch the entire process using value lists for the parameters, which failed to produce the proper outputs so I returned them all to single value inputs, yet the model still opens with the value lists present.   2. When I ran the entire batch it completed the process however when I checked my output values they were all copies of the original line, just sent to the various output location I had specified. I think this has to do with the fact the initial merge output is set as an intermediate file in a scratch.gdb so it is getting overwritten each time the model is run. I believe inline variable substitution can fix this yet I am unsure of how to do this as the model fails: invalid output name, when I tried the parse tool for the initial input file to use the name as an inline variable in the merge output name.  I have attached the model for clarity and my attempt at the inline variable substitution for the merge output. ![main model](http://i.stack.imgur.com/Fdlaa.jpg) ![output name using inline varaible substitution](http://i.stack.imgur.com/99uKO.jpg) Initially, I also had another pathway which converted the individual files from the output of the collect values tool into feature classes in my geodatabase, but after a while of trying to batch the model our link submodels each selecting for different wildcard the model would no longer create individual outputs it just produced a merge, even though all preconditions were set, and they were separate paths. Eventually I had to pull this path out of the model and run in a separate model. Is there a way to wipe all the past memory without deleting the entire model and beginning anew?


title: How do I pass each parsed file name to the next tool in Model Builder?

text: I cannot get the Parse Path tool to pass each file name to the next tool. In the example shown, the Raster to Ascii tool will take the first entry in the 'filename' variable list as part of the output name. But when it runs the next iteration, the Raster to Ascii tool uses some generic file naming convention. 'filename' is set as a precondition for the R2A tool. Example: raster1.tif raster2.tif Run the Model in Batch mode, using tifs as parameters. 1st output = raster1_asc.txt 2nd output = RastTo_rast.txt <\--- this is wrong. ![Model Builder running as Batch](http://i.stack.imgur.com/ruQc7.jpg) Summary: Raster to Ascii tool sucessfully uses the filename variable as inline variable as part of the output file name. However, the the Raster to Ascii tool uses a generic naming convention for the following files (running Model as Batch). Does anyone have an insight as to why this is happening?


title: is there any free software to create a mapserver that can be installed directly in a subdomain without a server install?

text: i want to create a webmap page using my own shapefiles but all the mapserver softwares that i've seen need to be installed directly in the server. so, is there any software that allows me to create the mapserver just by putting the files in some folder and "execute" them (like php or js, for example)? thanks in advance =^)


title: Is there a way to decrease redraws in ARCGIS?

text: I'm using a huge database on arcgis desktop 9.3.1 and every time I do a pan/zoom I have to wait 0.5-1 minute for the whole database to redraw itself. Is there any way to decrease the amount of redraws ( like making the program redraw only once every five pan actions)?


title: ArcGIS 10 SP1 Custom Overlay Grid: need to add a labeling expression

text: I have a custom overlay grid whose labels are driven by an ugly attribute field which contains unnecessary leading zeroes and other labeling no-nos. How can I introduce a labeling expression to this grid? Thank you!


title: Throwing E_FAIL errors so edit operation is aborted

text: I need to abort an edit operation somewhere, determined by a certain condition. Heres the code I'm using:                   private const int E_FAIL = -2147467259;              private void _Events_OnDeleteFeature(IObject obj)         {             _logger.Info("Iniciando handler de OnDeleteFeature...");             try             {                 // do we need to abort it?                 // this method returns a COMException with code = E_FAIL;                  if (DoINeedToAbort(obj, geaEditOperation.DELETE))                     return;             }             catch (COMException comEx)             {                 throw comEx;             }         }      This is not showing me the appropriate message and it's not aborting the operation. BTW, this event is a IEditEvent, and it's being fired ok. The problem is to abort the operation, without calling _editor.AbortEditOperation(), which gives me an ugly message. This is where I looked it in ESRIs manual (I know that the following applies to class extensions) http://resources.esri.com/help/9.3/arcgisengine/ArcObjects/ao_start.htm#esriGeodatabase/IObjectClassEvents.htm Any tips?


title: sum of raster counts

text: I'm trying to get the sum of "counts" in a raster. I have the code below which can successfully read the counts of each row in raster attribute and stores the counts in a list provided below. > [763461] [407207] [232269] [137428] [83909] [55074] [33340] [24168] [16747] > [11488] [7010] [4323] [4021] [2552] [1497] [775] [72] However, I want to get the sum of the counts. I am providing the code and error below. I would really appreciate your help/ suggestion! code:               import arcpy     from arcpy import env     from arcpy.sa import *          #To overwrite output     arcpy.env.overwriteOutput = True          #Set environment settings     env.workspace = "C:/Subhasis/Test/7-17-14/test"          #set local variables-STI extraction     inRaster ="sti"     inSQLClause = "VALUE >= 9"          #checkout ArcGIS spatial analyst extension license     arcpy.CheckOutExtension("Spatial")          # Execute ExtractByAttributes     attExtract = ExtractByAttributes(inRaster, inSQLClause)           # Save the output      attExtract.save("C:/Subhasis/Test/7-17-14/test/STI-9")          inraster ="sti-9"     rows = arcpy.SearchCursor(inraster,"","","Count","")          for row in rows:         count = row.getValue("Count")         s1=[]         s1.append(count)         print sum(range(s1))      Error:               Traceback (most recent call last):       File "C:/Subhasis/Test/7-17-14/test/test.py", line 31, in <module>         print sum(range(s1))     TypeError: range() integer end argument expected, got list.


title: Where can I download ArcGIS Service Pack 2?

text: I understand SP2 is available for download. http://blogs.esri.com/support/blogs/supportcenter/archive/2011/04/20/arcpad-10-0-service- pack-2-released.aspx > We are pleased to announce that Service Pack 2 for ArcPad 10.0 is now > available to download at the Esri Customer Care website. Maybe I am missing something but when I log into the ESRI customer care website the Service Pack is not listed. Anyone?


title: Creating data driven pages via arcpy

text: We have a data driven pages enabled MXD. Refreshing this mxd in arcmap works and gives me 9000+ pages. Via python however, _mxd.dataDrivenPages.pageCount_ returns 0. Even after _mxd.dataDrivenPages.refresh()_. In ArcGIS resources, they say: > Four elements are required for the creation of Data Driven Pages. They are: > > Detail Data Frame, Index Layer, Name Field, Sort Field In ArcMap I use the FILENAME option in the Name Field. In python, I find that _mxd.dataDrivenPages.pageNameField_ is _None_. But I'm unable to set this property. Any ideas on what to look for to get my pages via python? Thanks! Code is:               import arcpy, oss          logger = oss.log()          mxd = arcpy.mapping.MapDocument(MxdFile)     ddp = mxd.dataDrivenPages     logger.info('Checking required params')     logger.info('dataFrame is : "%s"' % mxd.dataDrivenPages.dataFrame)     logger.info('currentpageid is : "%s"' % mxd.dataDrivenPages.currentPageID)     logger.info('index layer is "%s"' % dpp.indexlayer)     logger.info('pagenamefield : "%s"' % mxd.dataDrivenPages.pageNameField)          logger.info('Refreshing index layer')     ddp.refresh()     logger.info('ddp number after refresh : "%s"' % mxd.dataDrivenPages.pageCount)      This results in:               INFO - Checking required params     INFO - dataFrame is : "<geoprocessing Data Frame object object at 0x0000000010CC     6970>"     INFO - currentpageid is : "0"     INFO - index layer is "VISION.SWG_WPA_MAP_ATLAS_TB"     INFO - pagenamefield : "None"     INFO - Refreshing index layer     INFO - ddp number after refresh : "0"


title: Spatial interpolation with breaklines in R

text: I have a dataset consisting of irregularly-spaced points (x,y,value) that represent an area of ocean that contains some land. I need to interpolate these to a regular grid, ensuring that the interpolation routine does not consider data points that are on the opposite site of a land barrier (e.g. the other side of a peninsular). I understand that this means that I need to insert "breaklines" in the data, and I have polygons of the coastline that I can use for this. Do any R packages allow for spatial interpolation using breaklines? Ideally I think that natural neighbour interpolation would be the most appropriate for my dataset, but I suspect that the breaklines are more important than the interpolation method.


title: Contouring surface with barriers?

text: I'm using ArcGIS Desktop 9.3, and I want to contour a surface with barriers. Every time the result is a polyline contour shape file which covers only a part of the surface not the whole surface. However, when I use "Contour" without barriers it works fine and the resulting contour shapefile covers the whole surface. Can anybody help with this?


title: QGIS documentation in French?

text: I am working on a plugin in qgis and i have to make a document from A to Z about that but on French . and I am looking now for some documentation of python,QT , pyqt4, and QGIS in French.


title: Repaint the map after updating a feature in QGIS

text: I update a feature in a vector layer using the following Python code:               layers = QgsMapLayerRegistry.instance().mapLayers()          layer = layers.get('my_vector20140614125509945')          fid = 1     attrs = { 0 : "hello", 1 : 100 }     layer.dataProvider().changeAttributeValues({ fid : attrs })      The feature is updated but the map is not redrawn. I am using catagorized style. The value of 100 will be shown with a different color.


title: Writing values to feature class in Python script results in strange characters"?

text: The script below does a series of things, and towards the end it is supposed to write values pulled from a list to a new field, if the feature class name is found anywhere in a list element. It does write values but there are strange characters that appear, or it doesn't write the entire list element value. I can't tell what the problem is?               import os, shutil,arcpy,arcgisscripting     gp = arcgisscripting.create()     path = r'C:\Temp\FME_share'     copy_dir = r'C:\test'     gdb_dir = copy_dir+'\\'+'CAD_Export_to_GDB'     BB_gdb_path = gdb_dir+'\\'+'BoundBoxes.gdb'     file_ext_lst = ['.2dm','.2de','.3dm','.3de','.dgn']     ext_to_chng = ['.2dm','.2de','.3dm','.3de']     file_lst = []     TxtVal_lst = []     repr_lst = []     found = dict()     arcpy.env.overwriteOutput = 'TRUE'               for dirpath, dirnames, filenames in os.walk(path):       for filename in filenames:           fullPath = os.path.join(dirpath, filename)           name,ext = os.path.splitext(filename)           if ext not in file_ext_lst:               continue           if name not in found:               found[name] = fullPath          print 'creating unique list of file names...'     for name,path in found.iteritems():       file_lst.append(path)          print 'copying files to local drive...'     for f in file_lst:       file_nm = f.split('\\')[-1]       shutil.copyfile(f, copy_dir+'\\'+file_nm)          print 'changing file extensions of copied files to .dgn...'     for dirpath, dirnames, copiedfiles in os.walk(copy_dir):       for copyfile in copiedfiles:           c_name,c_ext = os.path.splitext(copyfile)           filename_new = copy_dir+'\\'+c_name+'.dgn'           if c_ext in ext_to_chng:               os.rename(copy_dir+'\\'+copyfile, filename_new)          print 'exporting dgn files to geodatabase and exploding complex features...'     for dirpath, dirnames, renamed_files in os.walk(copy_dir):       for rn_file in renamed_files:           rn_name, rn_ext = os.path.splitext(rn_file)           in_dgn = dirpath+'\\'+rn_file           out_gdb = gdb_dir+'\\'+rn_name+'.gdb'           if rn_ext == '.dgn':               gp.ImportCAD_conversion(in_dgn,out_gdb,'','Explode_Complex')          print 'renaming geodatabase output to include original dgn file names...'     for dirpath, dirnames, filenames in arcpy.da.Walk(gdb_dir):       for filename in filenames:           gdb_name = dirpath.split('\\')[3]           orig_dgn_nm = gdb_name.split('.')[0]           if not dirpath.endswith('.gdb'):               arcpy.Rename_management(dirpath+'\\'+filename,dirpath+'\\'+filename+'_'+orig_dgn_nm)          print 'joining dgn text attributes to point feature classes, then adding text values to a list for searching out coordinate system information...'     for dirpath, dirnames, filenames in arcpy.da.Walk(gdb_dir):       for filename in filenames:           fullPath = dirpath+'\\'+filename           if 'Point' in filename:               Point_FC = dirpath+'\\'+filename             arcpy.JoinField_management(Point_FC,'EntID',dirpath.replace('CADStaging','')+'TxtProp','EntID','TxtValue')             values = [row[0] for row in arcpy.da.SearchCursor(Point_FC,'TxtValue')]             unique_vals = set(values)             unique_vals_str = ' '.join(str(uv) for uv in unique_vals)                  if ('CCS 83' or 'CCS83') and 'ZONE 3' and 'METERS' in unique_vals_str:                 print 'Trying to assign coordinate system to',dirpath,'....'                 arcpy.DefineProjection_management(dirpath,26943)             elif ('CCS 83' or 'CCS83') and 'ZONE 3' and 'FEET' in unique_vals_str:                 print 'Trying to assign coordinate system to',dirpath,'....'                 arcpy.DefineProjection_management(dirpath,2227)             elif ('CCS27' or'CCS 27') and 'ZONE 3' in unique_vals_str:                 print 'Trying to assign coordinate system to',dirpath,'....'                 arcpy.DefineProjection_management(dirpath,26743)             elif ('CCS 83' or 'CCS83') and 'ZONE 2' and 'METERS' in unique_vals_str:                 print 'Trying to assign coordinate system to',dirpath,'....'                 arcpy.DefineProjection_management(dirpath,26942)             elif ('CCS 83' or 'CCS83') and 'ZONE 2' and 'FEET' in unique_vals_str:                 print 'Trying to assign coordinate system to',dirpath,'....'                 arcpy.DefineProjection_management(dirpath,2226)             elif ('CCS27' or 'CCS 27') and 'ZONE 2' in unique_vals_str:                 print 'Trying to assign coordinate system to',dirpath,'....'                 arcpy.DefineProjection_management(dirpath,26742)             else:                 print 'Couldn\'t find coordinate system'          print 'creating minimum bounding box for line feature classes...'     for dirpath, dirnames, filenames in arcpy.da.Walk(gdb_dir):       for filename in filenames:           if 'Line' in filename:               BB_in_fc = dirpath+'\\'+filename               BB_out_fc = BB_gdb_path+'\\'+'BoundBox'+filename.replace('Line','')               arcpy.MinimumBoundingGeometry_management(BB_in_fc,BB_out_fc,'RECTANGLE_BY_AREA','ALL')          print 'adding and calculating SHEET_NAME and SHEET_PATH fields to bounding box feature classes, then reprojecting to WGS84...'     for dirpath, dirnames, filenames in arcpy.da.Walk(BB_gdb_path):       for filename in filenames:           sheet = filename.replace('BoundBox_','')           BB_fullpath = dirpath+'\\'+filename           arcpy.AddField_management(BB_fullpath,'SHEET_NAME','STRING')         arcpy.CalculateField_management(BB_fullpath,'SHEET_NAME',"'"+sheet+"'","PYTHON_9.3","#")         arcpy.AddField_management(BB_fullpath,'SHEET_PATH','STRING',200)         for f in file_lst:             if sheet in f:                 f_str = str(f)                 print f_str                 arcpy.CalculateField_management(BB_fullpath,'SHEET_PATH','"'+f_str+'"',"PYTHON_9.3","#")         BB_desc = arcpy.Describe(BB_fullpath)         BB_CS = BB_desc.spatialReference         BB_repr_file = BB_gdb_path+'\\'+filename+'_repr'         BB_repr_CS = arcpy.CreateSpatialReference_management('WGS 1984')         if '1927' in BB_CS.name:             arcpy.Project_management(BB_fullpath,BB_repr_file,BB_repr_CS,'NAD_1927_To_WGS_1984_79_CONUS')         else:             arcpy.Project_management(BB_fullpath,BB_repr_file,BB_repr_CS,'NAD_1983_To_WGS_1984_1')         repr_lst.append(BB_repr_file)     merge_fc = BB_gdb_path+'\\'+'BoundBoxes_All'     print 'finally, merging all bounding boxes into one feature class!'     arcpy.Merge_management(repr_lst,merge_fc)


title: Where can i get Shape Files for AE - U.S. Armed Forces – Europe AP - U.S. Armed Forces – Pacific PR - Puerto Rico?

text: I need shape files for these : AE - U.S. Armed Forces – Europe AP - U.S. Armed Forces – Pacific PR - Puerto Rico Where can i find them ?


title: Need help with GeoExpress 9 by LizardTech/ArcGIS Server 10.1

text: This is a multi-part question and after much reading on different forums I still could use some help/advice. My setup: ArcGIS Desktop 10.1, Dedicated Server for GIS data with SQL Server 2008 R2 and ArcServer 10.1. My problem: I have switched from ArcServer 9.3 to 10.1 which no longer supports ECW format for free. We can't afford the ESRI Image Extension or license for ECW. I need a way to convert my 5 years of Aerial Imagery in ECW mosaics to another format that is compatible with ArcServer, not too large (less than 30GB for each year's imagery). These are County wide datasets so pretty large when not compressed (~250GB for each year of Imagery). I also have individual tiles for each in MrSid, JPEG and Tiff format. I have not found a solution using ArcGIS tools so that the datasets are small, fast over the network/ArcServer, and "seamless". I've recently tried a free trial of LizardTech's GeoExpress 9 to create one MrSid mosaic of all the smaller MrSid files. This would have worked great and is only 24GB with minimal loss in quality, the problem is that the entire mosaic is covered in small semi-transparent hash marks (X X X X) that look like this... Individual tiles: ![Individual Tiles](http://i.stack.imgur.com/y5Q0N.jpg) The GeoExpress 9 (Trial Version) Mosaic: ![GeoExpress 9 version](http://i.stack.imgur.com/UaEgs.jpg) Has anyone seen this before? This is MrSid Gen4 could Gen3 help with ESRI compatibility? Could this be due to only having a trial version of LizardTechs software? I have not received an answer from them yet. Or if someone has used this software maybe I just did something foolish and used a wrong setting? The best solution for me, since I also use the ESRI JavaScript API for our WebApps, would be to use ArcGIS software to create some sort of image or mosaic of each year's imagery. One more option could be an OpenSource solution that supports ECW format? I could bring that service into my Web Apps and use ArcServer for all other layers, but I have not found one yet. I have created raster catalogs, merged the tiles into one raster, etc. with ArcGIS but the original tiles boundaries show through due to slight color mismatching. I'd like my end result to be one "seamless" image. I'm new to imagery processing so any help is much appreciated. Thanks in advance!


title: How to create Series of Rasters that displays Standard Deviation from Mean?

text: Cell statistics question for you. I have 100 rasters of raw rainfall data in which the pixels each represent rainfall for that particular pixel on the landscape for a single year during the 100 year period. I have also created a raster that calculated the mean rainfall received in each pixel over those 100 years. From these raw data rasters, I want to create a new set of 100 rasters (one for each year) that displays the standard deviation from the mean of each pixel (of the entire 100 years) for each pixel. I have a feeling this will be likely solved with python coding, which I have a little bit of experience in, but I am unsure of what tool to use? It seems like Cell Statistics is on the right track, since it focuses on individual cells, but doesn't achieve what I am looking for. Any other ideas on how to tackle this?


title: I am trying to find PDF documentation for Maplex for ArcGIS

text: Does anyone know where to find the PDF documentation for Maplex for ArcGIS 9.3 ? There was a website from ESRI but it was broken link ? I need to review it because I am working on impaired streams and lakes and needs their labels to set it up properly. Trying to do it in Label Manager several times. Thanks,


title: How to choose workspace when uploading SLD to Geoserver with the REST API

text: I have 200+ SLD's I would like to upload to Geoserver using the REST API. I successfully uploaded an SLD using the gsconfig Python module. However, I would like to be able to set the workspace when I upload. Is this possible? Below is my script so far.               from geoserver.catalog import Catalog          cat = Catalog("http://chchchch:8080/geoserver/rest", "username", "pw")          with open("C:\Users\Documents\SLD\my_sld.sld") as f:         cat.create_style("my_sld", f.read())


title: Should I use a CMS for my web mapping application?

text: For my web mapping app, I am using GeoServer as server, PostGIS as spatial DB and OpenLayers on client side (as well as other JS libraries). The map is composed of base maps and vector layers. One of these vector layers is the "Public institutions" layer. For each category of institutions (Ministries, Tribunals...) there are persons _(I called them "contributors")_ responsible for editing data (Tel, fax, location on map...) about these institutions. The contributors are managed by an admin who specify the privileges of each one of them. Some modifications of the contributors need to be validated by the admin (for example deleting an institution)... So I thought maybe using a CMS would help. What do you think? There are some CMS like Drupal that integrates OpenLayers, but it seems the project is still in development link to project page. Any alternative or experience to share?


title: How to get estimates of areas of overlapping polygons in one shp to another?

text: Sorry for my ineptitude, but I have two vector shape files of Africa. One has polygons of countries, and one has polygons of tribes a long time ago. What I want to do is to find the intersection of these two polygons, and get some idea of how many tribes were contained in a country today, and what share of land in the country today that they took up. Hopefully that makes sense. So, I'd like to end up with some data with a row for a country, say, Nigeria, that has: NIGERIA | Tribe 1: 20% of land | Tribe 2: 40% of land | Tribe 3: 40% of land | Something along those lines. Is this possible in QGIS or ArcGIS? My preferred software is QGIS.


title: GDAL.RasterizeLayer not ported to the C# bindings?

text: I am trying to use the `OSGeo.GDAL.Gdal.RasterizeLayer(…)` from _GDAL 1.9.2_ and it requires some `SWIGTYPE` parameters:               OSGeo.GDAL.Gdal.RasterizeLayer(         OSGeo.GDAL.Dataset dataset,         int bands,         OSGeo.GDAL.SWIGTYPE_p_int band_list,         OSGeo.OGR.Layer layer,         OSGeo.GDAL.SWIGTYPE_p_void pfnTransformer,         OSGeo.GDAL.SWIGTYPE_p_void pTransformArg,         int burn_values,         OSGeo.GDAL.SWIGTYPE_p_double burn_values_list,         string[] options,         OSGeo.GDAL.Gdal.GDALProgressFuncDelegate callback,         string callback_data     );      Reading up on some old mails somewhere ( _I've lost the link_ ) regarding ReadRaster and WriteRaster in the C# bindings, **Tamas Szekeres** states that methods requesting `SWIGTYPE` parameters, cannot be used (yet). **Question:**   Does this mean that the RasterizeLayer method hasn’t been ported to the C# bindings or are there some clever way to utilise this method?


title: GPS problem in QGIS

text: I have a very strange problem that I cannot understand. I am trying to put together a demo-laptop with QGIS and gps. On other computers I have it working as expected, but on this particular box something is strange. What I get is a gps position about 680 meters south east of my position. The position is fixed in QGIS, it is not "walking around" as expected. I see the altitude is changing though, but not lon lat. I have tried both connecting directly through /dev/rfcomm0 and through gpsd, same result. In other clients everything looks ok, like Foxtrot GPS and gpsmon. In gpsmon, and Foxtrot gps I get another position than in QGIS and the position in the other clients is moving around. I have moved the laptop and tried at another location and I get the same behavior but with another position some hundred meters away. I did get desperate and re installed the laptop, but get the same behavior. I am totally confused. It cannot be a hardware problem I think since it is working on other clients in the same laptop. It shouldn't be a software problem either since I have re installed everything including the OS. As mentioned I get it working on other computers with the same software. My setup is: OS Ubuntu 11.04 QGIS compiled from trunk a few days ago I can think of two reasons, but don't know how to search: 1) Something in my hardware gives another configuration when I am compiling QGIS than on the other computors 2) I have some software-package present (that I apparently also re installed in that case) that is interfering. Any ideas is highly appreciated Thanks   Nicklas


title: Issues with arcpy.sa.Reclassify

text: I have a raster file with classifcation codes, and want to reclassify valid codes to 1, and all others to 'NoData'.               RMAP2 = arcpy.sa.RemapValue([[3.2, 1],                                  [4.1, 1],                                  [4.2, 1]])          LCA = arcpy.sa.Reclassify(LCA, "LCCODE_USE", RMAP2, "NODATA")      The output however only contains the values 0 and 127, with the 0 value coinciding with just one of the missing values. Any idea what the hell is going on?


title: Help! - ArcSDE is locking our Oracle schemas

text: We changed a password on an Oracle schema, but apparently we have some rogue SDE clients - probably old ArcGIS map services on a dev machine somewhere - attempting to connect with the old password. SDE is causing Oracle to lock out the schema because of the rapidfire attempts to connect with the bad password. Is there any way to detect where the bad login attempts are coming from? I have turned on the ArcSDE intercept log, but I don't see anything related to connect attempts there. We're using SDE 9.3.1.


title: List all active extensions with ArcObjects

text: I am building an ArcGIS 10 Addin that dumps out as much debug information about a system that I can get. I am currently trying to query which extensions are active in the ArcMap session. The following code is to query on the active ESRI extensions, but does not show the 3rd party extensions. Any ideas of what I need to do to get the information on 3rd party extensions?               public string Extensions     {        get         {           var aoi = new ESRI.ArcGIS.esriSystem.AoInitializeClass();           var extensionsEnum = aoi.GetProductExtensions(esriLicenseProductCode.esriLicenseProductCodeArcInfo);           var licenseCode = extensionsEnum.Next();           var ret = "";           while(Convert.ToInt32(licenseCode) != -1)           {              if(aoi.IsExtensionCheckedOut(licenseCode))                 ret += aoi.GetLicenseExtensionName(licenseCode) + ",";              licenseCode = extensionsEnum.Next();           }           return ret+Environment.NewLine;       }     }      Here is the working code:               public string Extensions     {        get        {           var extensionManager = (IExtensionManager) ArcMap.Application;                var ret = "";           for (int i = 0; i < extensionManager.ExtensionCount; i++)           {               var extension = extensionManager.Extension[i];               var config = extension as IExtensionConfig;               if(config == null || config.State == esriExtensionState.esriESEnabled)                   ret += extension.Name + ",";           }                var jitExtensionManager = (IJITExtensionManager) ArcMap.Application;           for (int i = 0; i < jitExtensionManager.JITExtensionCount; i++ )           {              var extensionUid = jitExtensionManager.JITExtensionCLSID[i];              if(jitExtensionManager.IsLoaded(extensionUid))              {                 var extension = ArcMap.Application.FindExtensionByCLSID(extensionUid);                     var config = extension as IExtensionConfig;                 if (config == null || config.State == esriExtensionState.esriESEnabled)                    ret += extension.Name + ",";               }            }                 ret = ret.Substring(0, ret.Length - 1);            return ret+Environment.NewLine;         }      }


title: How to store Large Raster Images of Equipotential Surface Generation Using Monte Carlo Simulation?

text: I have a large (1GB+) raster that I am using as the base for the generation of an equipotential surface. This is to assess the accuracy of DEMs derived from another source. _In situ_ observation is impossible and the equipotential surface allows me to model in error. I have an ArcGIS model that generates a single realization of the surface. The model is outputting beautifully and generates a 1GB raster realization. I am going to generate 100 of these surfaces to start to build a statistically significant pool of realizations (Monte Carlo simulation). How would you go about storing the rasters? How would you visualize them to allow for the selection of the most likely realization (the one at the center of the monte carlo significance envelope)? Here are the options I have considered and I would appreciate any input, links to papers, etc.   1. Convert to point and store each realization as a column in the attribute table. Here I am concerned about the shapefile size limit (2GB).   2. Convert to CSV using gdal2xyz. I tried the conversion and the output of a single raster was 12GB! This is way too inefficient.   3. I considered using R but am not sure which package would facilitate this type of research (spatstat perhaps).   4. Graph either the histograms or cumulative distribution functions of each raster realization for visual identification.


title: JSAPI: How to prevent GetImageUrl from being called twice?

text: I posted this on the ESRI forum a while back but got no response. PROBLEM: GetImageUrl is called twice resulting in two trips to the server if caching is off. QUESTION: How can this be prevented (other than turn browser caching on)? HOW TO REPRODUCE:   1. Go to any sample with a dynamic layer (e.g. https://developers.arcgis.com/en/javascript/jssamples/layout_mapcontainersplitinthree.html)   2. Run a tool like fiddler2 to monitor IP traffic   3. Resize the browser   4. Notice that two identical requests are sent:     * "http://server.arcgisonline.com/ArcGI...%2C506&f=image"     * "http://server.arcgisonline.com/ArcGI...%2C506&f=image" A single call to map._resize results in getImageUrl being called twice (once for the resize and once for the change of extent). I am looking for a way to eliminate the redundant call with a minimal amount of coding.


title: C# version for ArcObjects 9.3

text: Can I use C# 4.0 with the target framework set to .NET 3.5 to develop an extension for ArcMap 9.3? Or does it have to be C# 3.0 or earlier?


title: Cant edit this featureclass. Why?

text: I can successfully create this featureclass in the ArcGIS10 default geodatabase and add it to the project, but it can't be edited. I have tried to add a feature using `FeatureClass.CreateFeatureBuffer()` but get the old cryptic: `{"Error HRESULT E_FAIL has been returned from a call to a COM component."}` If I try to add a feature with the editing toolbar in ArcMap: > The create feature task could not be completed. The featureclass DOES have a spatial reference that matches the data frame and the other layers in the project.   Am I missing something here? Should be easy, any ideas?               private IFeatureClass MakePointFC(IFeatureLayer pLineLayer)     {              try         {                  IGxCatalogDefaultDatabase Defaultgdb = ArcMap.Application as IGxCatalogDefaultDatabase;             Type factoryType = Type.GetTypeFromProgID("esriDataSourcesGDB.FileGDBWorkspaceFactory");             IWorkspaceFactory workspaceFactory = (IWorkspaceFactory)Activator.CreateInstance(factoryType);             IWorkspace pWorkspace = workspaceFactory.OpenFromFile(Defaultgdb.DefaultDatabaseName.PathName, 0);                            IFeatureWorkspace workspace = pWorkspace as IFeatureWorkspace;             UID CLSID = new UID();             CLSID.Value = "esriGeoDatabase.Object";                  IFields pFields = new FieldsClass();             IFieldsEdit pFieldsEdit = pFields as IFieldsEdit;             pFieldsEdit.FieldCount_2 = 2;                       IGeoDataset geoDataset = pLineLayer as IGeoDataset;                       IGeometryDef pGeomDef = new GeometryDef();             IGeometryDefEdit pGeomDefEdit = pGeomDef as IGeometryDefEdit;             pGeomDefEdit.GeometryType_2 = esriGeometryType.esriGeometryPoint;             pGeomDefEdit.SpatialReference_2 = geoDataset.SpatialReference;                            IField pField;             IFieldEdit pFieldEdit;                  //pField = new FieldClass();             //pFieldEdit = pField as IFieldEdit;             //pFieldEdit.AliasName_2 = "ObjectID";             //pFieldEdit.Name_2 = "ObjectID";             //pFieldEdit.Type_2 = esriFieldType.esriFieldTypeOID;             //pFieldsEdit.set_Field(0, pFieldEdit);                  pField = new FieldClass();             pFieldEdit = pField as IFieldEdit;             pFieldEdit.AliasName_2 = "SHAPE";             pFieldEdit.Name_2 = "SHAPE";             pFieldEdit.Type_2 = esriFieldType.esriFieldTypeGeometry;             pFieldEdit.GeometryDef_2 = pGeomDef;             pFieldsEdit.set_Field(0, pFieldEdit);                  pField = new FieldClass();             pFieldEdit = pField as IFieldEdit;             pFieldEdit.AliasName_2 = "Distance";             pFieldEdit.Name_2 = "Distance";             pFieldEdit.Type_2 = esriFieldType.esriFieldTypeDouble ;             pFieldsEdit.set_Field(1, pFieldEdit);                       string strFCName = System.IO.Path.GetFileNameWithoutExtension(System.IO.Path.GetRandomFileName());                  IFeatureClass pFeatureClass = workspace.CreateFeatureClass(strFCName, pFieldsEdit, CLSID, null, esriFeatureType.esriFTSimple, "SHAPE", "");             return pFeatureClass;                   }         catch (Exception ex)         {             return null;         }          }


title: ecw file with mapserver

text: I have a problem with ecw file. I'm using mapserver and openlayers. I have a ecw file and I make my tileindex file with gdaltindex. I want to put in my map file, I tried many ways but no ones work. How can I show my ecw files in my map? PS: I have support for ecw files by the way.


title: GeoExplorer provokes "401 Unauthorized" request and lacks login

text: Currently I have a login problem with my deployed GeoExplorer. Everything works fine in debug mode on my local VM with the geoserver from my remote server connected. Firebug gives the following output when loading GeoExplorer.               GET http://localhost:9080/geoserver/ows?SERVICE=WMS&REQUEST=GetCapabilities&TILED=true&VERSION=1.1.1 200 OK 140ms XMLHttpRequest.js (Zeile 220)          GET http://localhost:9080/geoserver/pdf/info.json 200 OK 18ms XMLHttpRequest.js (Zeile 220)      However, once I deploy GeoExplorer and move it to my remote server the login button on the upper right doesn't appear anymore. Instead the basic authorization window "GeoServer Realm" pops up and Firebug shows the following "401 Unauthorized" message:               GET http://Server_IP/geoserver/ows?SERVICE=WMS&REQUEST=GetCapabilities&TILED=true&VERSION=1.1.1 200 OK 227ms GeoExplorer.js (Zeile 396)          GET http://Server_IP/geoserver/pdf/info.json 200 OK 75ms GeoExplorer.js (Zeile 396)          GET http://Server_IP/geoserver/rest/imports 401 Unauthorized 202ms GeoExplorer.js (Zeile 396)          "NetworkError: 401 Unauthorized - http://Server_IP/geoserver/rest/imports"      The same problem applies to the GeoExplorer that was originally installed with the OpenGeoSuite, so it rather seems to stem from the server configuration instead of the deployment process. I've tried various things in GeoServer security setup and found threads with similar symptoms, but couldn't find a solution for us yet. Maybe someone has another useful hint?! Another idea: Can I configure the deployed GeoExplorer that it requests GeoServer over localhost instead of the server IP? Could this help?


title: Java GIS framework with custom shaped markers

text: I am completely new to GIS theme and I found this board so I hope I can get good advices from experienced users. I need to find an open-source Java framework/library/api that would generate an map image with multiple coordinates markers. In particular I need to customize the shape of the markers (i.e. making the marker circle size bigger or smaller depending on external additional data). Moreover, it would be nice that the map would fit to the markers (so if I place markers only in Europe, the image would show Europe only, not the whole globe). The target display is a java webapp so I would nee to generate the image and get it in a byte form, so I can display it on the page. Does any one know an Java API that has such functionality? I looked at things like geotools or jmapviewer but they seem to lack the features I require.


title: ArcGIS - Export Dialog Options missing

text: I've just unclicked the export map options (advanced options?) and can't get them back. I'd assumed you'd be given the choice of displaying them again, but they are nowhere to be seen. I now can't choose what file type or resolution to use, which is a huge pain. http://forums.arcgis.com/threads/42377-Export-map-no-options-for-details- filename.. Of the options above, I've tried deleting normal.mxt, with no luck. There's no 'ExportDlg' folder, so I don't have that option. I'm going crazy with this. I've Arc 10.1 running on Windows 8. Any workaround / solution would be hugely appreciated. I wont have the time to uninstall. Here's a screenshot: ![enter image description here](http://i.stack.imgur.com/aeMxd.png)


title: Creating hierarchy of Feature Elements in MapForce

text: I use FME to translate SQLServer tables into gml document. In FME, I make one- one connections (every feature type produces new feature type in output). This document is not 'structured' as I want. With that gml I need to create xslt file which will transform that gml int one which is valid with my needed schema. How can I define structure in MapForce to get done this? I asked question some time ago, without knowing problem's source. This question is here.


title: PostGIS upgrade 1.5->2.1 on OSX

text: i'm struggling updating my installation of PostgreSQL/PostGIS, from 9.1/1.5 to 9.2/2.1. I'm on OSX Mountain Lion and using the very useful ports from W. Kyngesburye. I followed the procedure :   * Updated PostGIS 1.5 to 2.0   * Dumped databases   * Updated databases to postgis 2.0   * Installed PostgreSQL 9.2 and PostGIS 2.1 I'm stuck at the step of using pg_upgrade i've got an error saying : > Could not load library `$libdir/postgis-2.0` ERROR: could not access file > "$libdir/postgis-2.0": No such file or directory > > Could not load library `$libdir/rtpostgis-2.0` ERROR: could not access file > `$libdir/rtpostgis-2.0`: No such file or directory Those libraries are in `/usr/local/pgsql-9.1/lib` pg_config returns that the lib directory is /usr/local/pgsql-9.2/lib If i switch the symbolic link /usr/local/pgsql to 9.1 again, same problem. Thanks for your help, i'm very eager to use the new functions of 9.2/2.0 :-)


title: Error while executing Populate_Geometry_Columns() Postgis

text: I loaded data into postgis using easy loader(MAPINFO). Now when I run               SELECT * FROM Populate_Geometry_Columns()      I get error : > WARNING: Could not convert '"SP_GEOMETRY"' in 'vectormaps."data_bnd"' to use > typmod with srid 0, type: MultiLineString CONTEXT: PL/pgSQL function > populate_geometry_columns(boolean) line 51 at assignment   1. How do i check this error?


title: How to create dynamic legend for multiply layers with leaflet.js?

text: I want to create multiple choropleth maps which displays different statistical data with leaflet.js. So it would be nice to show different legends for each map.   I create such legends for two layers using example in this question (How to Add/Remove Legend with Leaflet Layers Control?), but I can`t figure out how to handle more layers.   Here is example of my map Dynamical legend working for the "Agricultural" and "Building" layers, but I don`t know how to add dynamical legends (LISlegend, VODlegend) to "Forest" and "Water" layers.


title: Multiple Rasters into PostGIS

text: This may have been asked before, but A) what software/libraries etc do I need to B) load multiple rasters into PostGIS? I have Python 2.7.5, QGIS 1.9 (latest Weekly Build), PostGIS (whatever the latest version is) and I've downloaded the latest batch of GDAL codecs) Any help would be greatly appreciated Thanks Steve


title: How to Buffer with Dissolve ALL on very large polyline dataset?

text: So I have a fairly large polyline network covering a decent chunk of North America, containing just shy of 900,000 table records. I want to do some very simple "GIS 101" kind of processing to it -- I want to buffer everything with the "Dissolve" parameter set to "ALL". Well, unsurprisingly, this is way more than ArcGIS is ready to take on, even w/ 6 GB RAM & with everything in a geodatabase. I know the "Dissolve" operation is more prone to running out of memory in such cases than "Buffer", so I tried to run "Buffer" without "Dissolve = ALL", but even that will not generate a usable output (the file exists but won't draw or reveal its attribute table). And, in any case, I am sure the subsequent "Dissolve" would hang. Here's what I've tried:   1. Breaking up the network into smaller pieces. Well, this seems like a non-starter, because I need to make sure I don't cut any line segments apart in doing so. Of course selecting polyline segments sequentially from the attribute table will just select random pieces of polylines all over the place, shredding the integrity of the network. I'm not sure how to break it up programmatically such that no connected pieces are broken apart.   2. Running the buffer & dissolve in QGIS, based on some hearsay that it was better at handling large datasets. Like ArcGIS, it also just hung there forever.   3. Importing everything into R with the "maptools" package and seeing if I could do it there. The data is so huge I couldn't even read it in within a reasonable time span. There's a whole lot of talk out there about this kind of situation, but it always ultimately ends with "you need to find a way to break up your large dataset." I'm not really sure how this can be done cleanly with a huge network of stream polylines, barring some crazy Python script that can tell, for each polyline segment, which ones it's connected to, and select those only, until the entire spatially discrete sub-network is selected, write that to a single feature, & then move on to the next such sub-network. This seems absurdly complex as a workaround for very simple processing on a not-even-all-that-big dataset (though I'll try it if I have to), but can anyone tell me how I can get this thing buffered?? Other software / workarounds / magic incantations / human sacrifice?? (And there'll be yet more steps after the buffer/dissolve, which would be no sweat with smaller data but will probably cause a minor earthquake if I try them with this beast). **EDIT :** After exploring this much more in-depth, I feel this question hasn't been adequately addressed by the duplicate question or the main suggestions here. Here's why:   1. The Dice tool appears to be for reducing the number of vertices in individual features? I don't have individual streamlines with >3,000 vertices, but in any case, I've used Dice with a vertices limit of 100, and still cannot run the dissolve.    2. Breaking the streams into 100 pieces and dissolving each one requires an "uber-Dissolve" at the end to combine the pieces, which will have buffers overlapping with one another. This "uber-Dissolve" experiences the same memory failure as the original approach, even though there is much less to dissolve after the 100 pieces have been processed. So this approach indeed doesn't seem to work, though for a completely different reason than what I stated earlier. It seems PostGIS has implemented an approach called "cascaded unions" that addresses this general problem, which will be fun to explore, but in the meantime I was hoping for an ArcGIS-/Python-/QGIS-based solution.


title: How to remove some features of a feature layer in the mobile cache

text: I have a Hand-Held (HH) device in which I have downloaded to the MobileCache a limited amount of Features from a FeatureLayer. Now I need to remove from the cache all features of a certain type. We were able to delete them from the cache but this places them into a table called `fs_table0_deleted_features` in the `MobileCache.db` file. We are only doing a one-way synchronize (download), but if we ever did a two-way synchronize then these features would be deleted from the server too, which is not what we want. Also, if we try to re- synchronize the now deleted features from the server they are not removed from the deleted features table and the deleted features are not placed back into the `fs_table0_current_features` table; this behavior is causing issues with our application. We also tried to use the MobileCache.RemoveData method, but this only accepts an Envelope (extent) and a list of FeatureSource. What I want to do is only remove a subset of Features that have been cached leaving the rest cached. Other solutions we have considered have been:   1. Copy in memory from the MobileCache the features we want to save, Call RemoveData, then insert into the MobileCache the features that we saved.   2. RemoveData from the MobileCache and then synchronize with the server for all the features we want   3. Delete the MobileCache and then copy a backup that was previously made to replace the now deleted MobileCache.   4. Find the Extent of all the features we want to remove and then call RemoveData with those extents. Of these, I think #1 or #3 may the better solutions, but they seem overly complex and it could be resource intensive(#1), or may have problems with stale data(#3). We don't like #2 because the load for synchronizing our devices may be heavy (thousands of devices all must synchronize in a very limited time frame) The last one, #4 could be a potential problem if there are overlapping features in which we only want to remove one and not all. Here is how we tried the delete, which is not working and isn't a valid solution.               public void PurgeCacheData(string whereClause)     {         if (!Cache.IsOpen) return;              var qFilter = new QueryFilter { WhereClause = whereClause };         var featureSources = Cache.FeatureSources;              foreach (var fSource in featureSources                  .Where(fSource => fSource.GetFeatureCount(qFilter) > 0))         {                  using (var featReader = fSource.GetDataReader(qFilter))             {                 while (featReader.Read())                 {                     featReader.Delete();                     featReader.Update();                 }             }                  // Here we tried to delete the entries in the Delete table             // but this doesn't work. Nothing happens.              using (FeatureDataTable featureDataTable =                     fSource.GetDataTable(null, EditState.Deleted, null))             {                 foreach (FeatureDataRow row in featureDataTable.Rows)                 {                     row.Delete();                 }             }              }      }


title: Python workspace

text: I wrote a simple script to rename about 15 shapefiles in a folder on our local server. When I tested it I was working on a workspace on the C:\ and it worked fine. All our data is on a local server and when I changed the pathname to (... arcpy.env.workspace = "P:\Alabama\Export") it did not work. I eventually have to use 20 or so diff paths so I don't want to have to copy data to my c drive. As I read through examples they all seem to use something on C:. I am running it in the Python window in ArcCatalog not on an IDE Sample:               import arcpy          arcpy.env.workspace = "P:\Alabama\Export"           if arcpy.Exists("Histori2.shp"):         arcpy.Rename_management("Histori2.shp", "H_BEV_CANS.shp")           if arcpy.Exists("Histori3.shp"):          arcpy.Rename_management("Histori3.shp", "H_CANS.shp")


title: How to install Openlayers plugin on Qgis 1.8?

text: After having installed Qgis 1.8 i am not able to find the Openlayer Plugin. In "the old days" I know I was supposed to check "add 3rd party repositories", but that button has disappeared? I tried adding the following repositories without luck:               http://pyqgis.org/repo/contributed     ..openlayers.org/api/2.12     ..github.com/sourcepole/qgis-openlayers-plugin.git     ..github.com/nvkelso/qgis-openlayers-plugin.git     ..plugins.qgis.org/plugins/plugins.xml     ..plugins.QGIS.org/plugins/plugins.xml     ..github.com/openlayers/openlayers      I hope someone can help me locate the OpenLayers plugin as I desperately need it.


title: How can I select "ultimate perimeter" from a group of connected lines?

text: I have a number of drawings which I converted into geodatabase feature classes. Most drawings are composed or connected 3D polylines from which I want to construct TIN surfaces and multipatch features. In order to clip the TIN I also need to have the exact footprint. How I normally do it is to manually select the perimeter polylines which I then merge into one perimeter and convert to polygon. This works but is extremely time consuming because I often have to manually select 100s of polylines without missing a single connecting line. **Can I somehow run a spatial selection to only select the polylines that form the ultimate perimeter?** Below is a typical example - i need to select only the perimeter polylines: ![enter image description here](http://i.stack.imgur.com/v68eO.jpg)


title: Developing a GIS viewer

text: I was wondering if you could give your input/opinion on what to use. I am looking into developing a small scale GIS viewer that could possibly edit the information locally (Edit would entail drawing on the map, but not data) and searching for assets on that map. Is ArcReader something that can be customized with add-ins? thats something I was looking into. What id really like is for a program to be able to view the maps with simple user tasks (drawing, searching...) that is portable. any ideas or direction on where to start? I was thinking the arc run-time (java). Any information would be helpful, thanks :).


title: What is the 'water retention capacity' tool and how is it used?

text: In the latest version of QGIS, I found a geoalgorithim in the processing toolbox that will generate 'water retention capacity.' I am familiar with DEM-based hydrologic modeling tools in QGIS, but I cannot find any information online regarding the formula for water retention capacity. Additionally, a 'plot holes' file is required along with the DEM and I don't know what that is supposed to represent. Thanks for any help.


title: Unable to view a wms layer in QGIS cpp aplication

text: I am a newbie in QGIS c++ programming, and I am trying to load a raster file from a WMS and display it. The image I want to display can be obtained from this url http://geoserveis.icc.cat/icc_mapesbase/wms/service?REQUEST=GetMap&VERSION=1.1.0&SERVICE=WMS&SRS=EPSG:23031&BBOX=290368.84,4538236.42,292203.28,4540070.86&WIDTH=520&HEIGHT=520&LAYERS=orto25c&STYLES=&FORMAT=TIF&BGCOLOR=0xFFFFFF&TRANSPARENT=TRUE&EXCEPTION=INIMAGE and I can display it without problems, if I download it first on my hard drive, and then I load it properly using a QgsRasterLayer object. However, I would like to acces this image directly from the wms server. My problem is that when I do that, nothing is displayed on the map canvas. I copy here my C++ code. Hopefully somebody can tell me what I am doing wrong. QgsDataSourceURI uriDS; uriDS.setParam("url", "http://geoserveis.icc.cat/icc_mapesbase/wms/service?"); uriDS.setParam("layers", "orto25c"); uriDS.setParam("format", "image/tif"); uriDS.setParam("srs", "epsg:23031"); uriDS.setParam("styles", ""); uriDS.setParam("bbox","290368.84,4538236.42,292203.28,4540070.86"); uriDS.setParam("width","520"); uriDS.setParam("height","520"); QString uri = uriDS.encodedUri(); // I put 'false' as the 4th parameter, because otherwise the Layer is not valid QgsRasterLayer *mypLayer = new QgsRasterLayer(uri, "layer_name", "wms", false); if (mypLayer->isValid()) { qDebug("Layer is valid"); } else { qDebug("Layer is NOT valid"); return; } connect( mypLayer, SIGNAL(repaintRequested()), mpMapCanvas, SLOT(refresh()) ); QgsMapLayerRegistry::instance()->addMapLayer(mypLayer, TRUE); QList myList; myList.append(QgsMapCanvasLayer(mypLayer, TRUE));//bool visibility mpMapCanvas->setExtent(mypLayer->extent()); mpMapCanvas->freeze(false); mpMapCanvas->setLayerSet(myList); mpMapCanvas->setVisible(true); mpMapCanvas->refresh();


title: How to put a roads layer (OpenStreetMap) on top of KML polygons?

text: I need to develop a map where KMLs are overlayed with road layers. So, KMLs are polygons, and I would need to use a map tiles like OpenStreetMaps. But the problem is how to put roads on top of KML areas? What do you suggest me? I would like to use some open source software. Thank you.


title: postGIS union multiple tables, big dataset, faster approach?

text: my goal is simple. I want to union a given set of shapefiles. In the shapefile are stored polygons. The data bounding box is germany, so the data is relatively big. I tried to use all functions implemented in QGIS. Either they failed or the precessing time is just too long to get any results. I read a lot of question throughout this forum and i set up a PostGIS Database. My Query is:                CREATE TABLE dlm_union_2 AS       SELECT       a.gid AS ID,       ST_UNION(a.geom,b.geom) AS geom        FROM sie01_f a, sie02_f b      I added the Selection of the a.gid since i have no objects returned when i export the table via pgSHAPELOADER. I want to use more than two tables in one expression and is there a way to get this faster? With only two tables my processing time is now about 800,000 ms and the progress isnt finished yet. After that i want to buffer the merged regions and would like to union it again afterwards. I have a lot of shapes that i have to union at different states of my project. Iam working on WIN 7, 4GB ram, Intel Core i5 M520. PostGIS and Postgres should be up to date, because i installed the components today. Every suggestion are highly appreciated. Thx Steffen UPDATE: I have created a buffer on a multilinestring table. Processing time was about 300,000 ms what is actually great compared to functions implemented in QGIS. My problem is that i need to union the result for further analysis. I'am using ST_UNION with a single table with polygon geometry input. Objectsize is around 300,000 polygon features. At this point my processing time is 2,800,000 ms and i dont know when it will finish. Is there any way to get such a calculation faster? After processing time 5,319,543 ms the query aborted with following exeption. This was the try to use ST_UNION with one table including 300,000 polygons.               ERROR:  GEOSUnaryUnion: std::bad_alloc          ********** Fehler **********          ERROR: GEOSUnaryUnion: std::bad_alloc     SQL Status:XX000      UPDATE: According to Mattmakesmaps i tried the way of ST_Buffer(ST_Collect(geom),0) instead of using the UNION command. With the errors iam posted above i was not able the handle the whole extend of the data. So i split the problem into smaller ones. In my case i used the state-borders in germany to divide my problem into 16 smaller ones. Mattmakesmaps solution works fine here. Recently iam try to use the ST_Simplify command to make the dimension of the polygons smaller giving me a chance to process all of the data at ones. I will keep you up to date. Thx for helping me out. Beste regards Steffen


title: Conversion of tilestache cache into mbtiles

text: I have read and I am looking for someone with documentation on how to achieve converting a cache generated by tilecache into mbtiles using mainly mbtiles. I see the mbtiles has an import function but it requires a metadata.json which tilecache does not have and also it seems the --scheme option has been depricated. I am trying to use mb-util -i tiles world.mbtiles and hopefully I would like to open the mbtiles in qgis


title: how to automate changing polygon's shape while keeping its area the same and limiting deviation from orginal centroid?

text: I have a polygon that was digitized based upon visually interpreting underlying values of rasters/vectors that represent a natural environment (soil type, surface roughness, slope, elevation, etc...). The goal of the digitization was to define a polygon shape that maximizes the average value of underlying raster/vector values underneath while keeping the polygon area relatively small. In this type of exercise there are many possible answers. I am brainstorming ways I could possibly automate tweaking the shape of the initially digitized polygon to see if there is a slightly better answer while defining some limits to the problem. Here is an example of a workflow I would like to have: 1) Begin with the initial polygon that is say, fairly circular with an area of 1 km^2 and a centroid at 0,0 2) Use a script or tool to create/try hundreds of slight permutations of the polygon's shape while having two limitations i) keeping the area the same and ii) not allowing the centroid to deviate more than say 0.5 in both x and y direction. 3) With the new polygon shape, perform an Extract by Polygon with one of the raster files and check if the average value of the cells has increased in the new polygon area compared to the original, and if so, save the new polygon to visually check later. I figure the solution could get very complex because it deals with nodes. Certainly there would need to be creation or definition of a parameter describing how rough/erratic changes in the shape are allowed to be, and also, the calculations would depend on the number of nodes in the polygon. There might also be issues with self-intersecting lines depending on the parameters chosen. Another addition to the approach would be having the changes in shape intelligently maximize the average value rather than change the shape randomly. Any other factors important to account for / think about? Do you know of some existing useful tools/scripts that could solve or help solve this problem? Any ideas are welcomed - especially ArcGIS and Open source approaches. I have a feeling that an ArcGIS solution will most likely involve programming with ArcObjects rather than any scripting. Or perhaps it's possible that a permutation problem like this one is best handled by non-GIS software - maybe something like R?


title: Where can I find North American Pipeline Data?

text: I work for an energy management company and need to couple this data with natural gas projections. Any help would be appreciated!


title: Loading / georeferencing OS rasters in QGIS

text: I have got some OS rasters but there isn't georeferencing file (TFW) with them. When I try to load them, it will obviously load in the wrong place. They are loading fine in a different GIS software (CadCorp) because the tile name is recognised (example:SJ86SW). Is that possible in QGIS? If not, while georeferencing the file can I snap to the corners of the raster to get the absolute position of these points?


title: Symbology for specific values

text: I want to apply a symbology scheme for the final layer of my model. Rules are -               All negative - Symbol A     0 - Symbol B     All positive symbol C     unknown - symbol D      I am not sure how I can define these intervals under Classification of Layer Properties. The minimum and maximum can vary since it will be used in a model. I tried to define manually (among other options - defined, standard deviation...etc.) but it worked only for that sample and failed for another example. I am using 'Apply Symbology from Layer' to display the model output finally. Thanks


title: mysql foreign data wrapper for postgres with parameter substitution

text: I have created a mysql_fdw for my postgres and when i try to run the one with a subquerry that has parameter substitution it does create a table but on querying the foreign table it throws an error:   CREATE FOREIGN TABLE overtime_2010 (   id integer,   employee_id integer,   hours integer)   SERVER mysql_svr   OPTIONS (query 'SELECT id, employee_id, hours FROM hr.overtime WHERE year = ?year;'); The sql creates the foreign table but executing a select query yields The following error: connection to the server was lost. Attempting reset: Failed. Also trying to create a local table also fails with the above error Create table local as select * from overtime_2010;


title: How do I successfully install LeCos plug in for QGIS2.0 on a Mac?

text: I have been trying without success to install the LeCos plug in on QGIS2.0 on my Mac. I originally tried it on QGIS1.8 but was unable to install the dependent packages (PIL,scipy and numpy). I updated to QGIS 2.20 and successfully installed the packages, however I am still unable to install LeCos. THis is the message that comes up: ![enter image description here](http://i.stack.imgur.com/FiQhz.jpg) The packages are definitely installed because I checked by doing this in python: ![enter image description here](http://i.stack.imgur.com/GsxJX.jpg) Has anyone else installed it successfully?


title: What is a good plotter that works well with ArcMap?

text: I will soon be setting up a small office and we will be needing a 36" plotter. I am looking for recommendations for a good plotter that works well with ArcGIS. Speed is not important as long as I can produce good maps for field work and analysis. (Budget = < $5,000 CAD)


title: How can I download images from an API using FME?

text: Using an API which enables me to retrieve XML data about Objects, I am also able to enter the ID of the object and download an image. Does anyone know how I can do this using FME?


title: Interpolation of xyz data

text: Data has been given to us in a HUGE .xyz file, which has been split into many component CSV files. When I upload the CSV file, then create the shape file by "save as", all is fine. When I go to Raster, then Interpolation, the resulting layer appears incomplete, as if it hasn't been processed. - Image below... ![enter image description here](http://i.stack.imgur.com/IXdlP.png) Question is: is it incomplete due to lack of "machine-power" (very doubtful) or data issue (doubtful, as it appears OK in the CSV) or is it a limitation of QGIS? Is there an easier way of uploading the enormous files rather than splitting them and risking a data error?


title: Oracle error when copying massive feature dataset

text: In ArcCatalog 10.1, when copying a feature dataset with roughly 3 million objects as well as a few thousand tables and relationship classes, my team is running into an > 'Underlying DBMS error [ORA-03114: not connected to ORACLE][SDE.DEFAULT]' message 9 times out of 10... I'm hoping someone has run into a similar issue. We are literally copying the feature dataset from one connection file and pasting it into another connection file/environment. This works fine with smaller sets of data. However, for this bigger dataset transfer, we are running into this problem - the process copies all objects but fails in the late stages (I think when it is copying tables). This is killing us, since the process takes about 20 hours and then bombs! X( I also get a Paste Failed Message each time: > 'Failed to paste GEODATA.NA_US_SOIL_MUPOLYGON_2013 Invalid function > arguments Underlying DBMS error [SDE.GDB_ITEMS] Underlying DBMS error > [ORA-03114: not connected to ORACLE][GEODATA.NA_US_SOIL_MUPOLYGON_x] > Underlying DBMS [ORA-03114: not connected to ORACLE] Underlying DBMS error > [ORA-03113: end-of-file on communication channel Process ID: 7804 Session > ID: 205 Serial number: 52647][GEODATA.NA_US_SOIL_MUPOLYGON_x] We've discussed this issue at length with Oracle DB team here and cannot seem to find any issue in Oracle logs. Our OS is Windows Server 2008 Std64 bit SP2, ArcSDE 10.1 SP1, and the server is a DellR710. Has anyone encountered a similar error? We're at our wits' end here. Thanks for any suggestions!


title: Display freight volume by city on a geographic map

text: I have been requested to show on a geographic map the volume (in meters cubed m3) of freight (parts) that ship to various cities (Canada/USA). The idea is to show which cities receive the most freight by volume. Presently the only mapping tool I have accessible is Microsoft Mappoint 2002. Is there another webbased tool that I could use that does a good job of showing volumes relative to each other? Mappoint has a feature to show sized pie charts - is the best type of method of showing total volume by city along with the individual contributors? Update1: I am not trying to show the flow of parts - only the size of origin and destination. For example does Toronto, Ontario have 100m3 of parts and Buffalo, New York only 30m3 of parts?


title: How does integrate work with point data sets?

text: I'm trying to combine point data into clusters and I think I've come up with a method that will work for what I want to accomplish. However, I don't understand why it works the way it does and I'm hoping someone can explain it to me. Basically I have a set of wells that are operated by different companies and are scattered throughout PA. I want to determine how clustered the wells are by operator in terms of how many clusters there are for each operator and how many wells are in each cluster (clustering at a distance of 5 miles). It appears that if I run integrate on my point data for the wells it will do this. Then I can collect event values to determine how many wells are in each cluster. What I don't understand is why some of the data is broken into more than one cluster. ![Figure](http://i.stack.imgur.com/0EwFz.png) [Image Description] The three highlighted points are the integrated points for the company represented by the light yellow well sites. When you open the attribute table for the integrated data it is clear that the group on the right is one cluster and that integrate has determined that the group on the left is two clusters split on the blue line. What I don't understand is why the points immediately above and below the line have been assigned the clusters that they have been. They are all within 5 miles (the integrate thresh hold) of the lower integrate point. And the upper ones are more than 5 miles from the upper integrate point. So why split it there? Any thoughts would be helpful. * * * To help clarify I've added a screenshot of the process that I used to integrate all the points into the clusters that I would like to use. http://imgur.com/swzxW


title: How to cancel tool execution in python

text: I've written a python function that calls a tool (Import CAD to GDB) in a loop. When I run it I have seen that it hangs sometimes, and I have to force the script to quit. I am wondering if there is a way to kill the tool execution after a certain amount of time, and move on to the next item in the loop? Here is what I have so far - I've included a timer but it still hangs:               import os, arcpy, arcgisscripting, time, sys     gp = arcgisscripting.create()          def ConvertCADtoGDB(srch_dir,gdb_dir,timeout):         for dirpath, dirnames, files in os.walk(srch_dir):             for file in files:                 f_name, f_ext = os.path.splitext(file)                 in_dgn = dirpath+'\\'+file                 out_gdb = gdb_dir+'\\'+f_name+'.gdb'                 if f_ext == '.dgn':                     if not os.path.isdir(out_gdb):                         start = time.time()                         end = start + timeout                         while time.time() < end:                             print 'creating',out_gdb                             gp.ImportCAD_conversion(in_dgn,out_gdb,'','Explode_Complex')


title: Merge multiple feature datasets while maintaining shapefiles for respective point, line, and polygon files nested ArcGIS

text: I've been doing some digging and a few days of playing around with my data in ArcInfo but I haven't been able to find a solution to merging datasets, which consist of many different point, line, and polygon shapefiles. The files are electronic navigation charts in a vector format (some examples of these are freely available: http://www.charts.noaa.gov/ENCs/ENCs.shtml). Each of these charts covers a different spatial extent or resolution and contains information relevant to mariners, fishermen etc. (e.g., depth, substrate type, buoys, lights, shipwrecks, obstructions, and so on). There is a standard for the display of these data called S-57 in which each feature is universally coded (e.g., obstructions are OBSTRN in each ENC file: http://www.s-57.com/). The coding of the attributes of each feature is also standardized. Therefore, each file you open in Arc will have identically named point, line and polygon layers for the features on the chart (although some charts may have features not found on other charts [e.g., some areas might not contain shipwrecks]). I am wanting to merge many of these ENC files, which cover portions of my region of interest, into a large regional layer. My hope is to use the respective shapefiles (and attributes) for this whole region rather than for the smaller areas the charts cover individually. The merge has to maintain all the respective point, line, and polygon shapefiles (e.g., depth, substrate type, buoys, lights, shipwrecks, obstructions) (i.e., it is basically a merge of all identically named shapefiles nested within each chart). I know the merge can be done between two respective shapefiles from 2 different ENC chart datasets (after a name change of one of them), however, I have potentially 20-30 charts to merge and there are 20+ shapefiles in each chart file. Any direction would be greatly appreciated. Thanks in advance. (I attached a screenshot of the ENC file in ArcInfo. Each chart file has these identical layers covering a different spatial extent I wish to merge). ![Screenshot](http://i.stack.imgur.com/qQ1cL.jpg)


title: Error importing XML workspace document to file geodatabase

text: Ok so I have a large SDE geodatabase that I exported to an XML Workspace Document. I tried to importing it again into a file geodatabase to test if the export worked correctly and got the following error: Import xml data failed   FDO error: -2147186428   [XML parse error: Illegal xml character   [Line: 23612914, column: 1043].] Could this be from an invalid character somewhere in my data? Is there a list of invalid characters for XML workspace documents? And is there an easy way to get to the line and column in the error message to see what it is? Also are there any file size limitations, my XML is nearly 12gb? Or could it be possible there was just an error in the exporting process? Any thoughts appreciated. Thanks


title: Debugging r.neighbors error in Processing

text: I'm having similar problems withs some GRASS algorithms in processing plugin. On the one hand i have installed Ubuntu 12.04 with Qgis 2.0.1. In this machine when I run _r.neighbors_ algorithm don't make it and log output says this:               processing.runalg("grass:r.neighbors","/home/user/Documents/Projects_GIS/Sgt/MDE-Sgt.asc",0,3,False,False,"","693360.0,698725.0,4331760.0,4336495.0",5,"/home/luismi/Documentos/Curso_Avanzado/testing/segmentacion/test34.tif")      but the curious thing is this... 11 sec later it starts another algorithm _v.voronoi_               processing.runalg("grass:v.voronoi","/home/user/.qgis2/python/plugins/processing/tests/data/points.shp",False,False,"270778.60198,270855.745301,4458921.97814,4458983.8488",-1,0.0001,0,None)      And last _log messages_ says this:               Data Source is invalid()      On the other hand I'm using OSGeoLive to compare both software: It happens the same trouble, the same two algorithm running, _r.neighbors_ , and spontaneously _v.voronoi_ , but in this case the log message is this:               OGR[3] error 4: Unable to open /usr/share/qgis/python/plugins/processing/tests/data/points.shp or /usr/share/qgis/python/plugins/processing/tests/data/points.SHP.      Seems there is something wrong inside, but I don't know what is it...


title: How to make XY to Line dynamic to access query?

text: I have connected the Access file to the ArcGIS file, and displayed it using the tool XY to line, but when I change the access file, this is not automatically adjusted in the ArcGIS template. Somehow, ArcGIS makes the XY line static, so it does not have the possibility to change. What should I do to make the produced lines 'dynamic'? Thanks a lot


title: OpenLayers Permalink getCenter() Returns strange values

text: I'm extending the openlayers permalink class to create a custom class that can create a link on the map with an href attribute that points to http://www.openstreetmap.org/edit?editor=id#map=15/11.45311/122.84363 Where: \- editor=id can be changed to other editors (potlatch2 ?) \- 15 = the current zoomlevel \- 11.xxx = lat \- 122.xxx = lon; the point is to have my new class calculate the zoom, lat and lon of the map center and pass these on to the parameters in the url mentioned above. My problem is the call to getCenter() [First line in the createParams function]:               createParams: function(center, zoom, layers) {       center = center || this.map.getCenter();            var params = {};            // If there's still no center, map is not initialized yet.       // Break out of this function, and simply return the params from the       // base link.       if (center) {         //zoom         params.zoom = zoom || this.map.getZoom();              //lon,lat         var lat = center.lat;         var lon = center.lon;              if (this.displayProjection) {           var mapPosition = OpenLayers.Projection.transform(               { x: lon, y: lat },               this.map.getProjectionObject(),               this.displayProjection);           lon = mapPosition.x;           lat = mapPosition.y;         }         params.lat = Math.round(lat * 100000) / 100000;         params.lon = Math.round(lon * 100000) / 100000;            }       return params;     }      this.map.getCenter() returns strange values like: lon=0, lat=-7.0811545516136e-10, lon=-0.0015999972820282, lat=0.0015999972820282, lon=-0.0023999959230423, lat=0.0023999959230423 which end up as lon = 0 and lat = 0 in the params.lat and params.lon variables that create the parameters for the url mentioned above. The zeros display an empty map as shown here http://www.openstreetmap.org/edit?editor=id#map=15/0/0? What could be the issue with the getCenter() call ?


title: How to speed up this simple MySQL points in the box query

text: MySQL doesn't have nearest neighbor search so. So I try to come up with a "distance" that's reasonable so the number of returned query is around 20-200. That way innodb doesn't have to compute distance for like 30k rows if we only need 20. The latitude and longitude is also stored in a myisam table called tableauxiliary Well, even after that, queries still take 13 seconds and sometimes 30 seconds.               SELECT TB.ID ,       TB.Latitude,       TB.Longitude,       111151.29341326*SQRT(pow(-6.185-TB.Latitude,2)+pow(106.773-TB.Longitude,2)*0.98839228980165) AS Distance     FROM       `tablebusiness` AS TB       join tableauxiliary as TA on TA.BusinessID=TB.ID     WHERE       MBRContains(             GeomFromText (                 'MULTIPOINT(-6.1928749092968 106.7651250907,-6.1771250907032 106.7808749093)'                 ),                 TA.Latlong)        AND MATCH (TA.FullTextSearch) AGAINST ('re*' IN BOOLEAN MODE)      AND TA.Prominent >15      ORDER BY       Distance      Number of matches are 32 by the way.               Showing rows 0 - 29 ( 32 total, Query took 16.1845 sec)      Also an even simpler query where we do not compute distance also take 16 seconds               SELECT TB.ID     FROM       `tablebusiness` AS TB       join tableauxiliary as TA on TA.BusinessID=TB.ID     WHERE       MBRContains(             GeomFromText (                 'MULTIPOINT(-6.1928749092968 106.7651250907,-6.1771250907032 106.7808749093)'                 ),                 TA.Latlong)        AND MATCH (TA.FullTextSearch) AGAINST ('res*' IN BOOLEAN MODE)         AND TA.Prominent >15       I can understand if mysql cannot do nearest neighbor search. What I am doing now is a simple search points in a box here. What am I missing? The explanation for the query is:               id  select_type     table   type    possible_keys   key     key_len     ref     rows    Extra     1   SIMPLE  TA  fulltext    PRIMARY,Prominent,LatLong,LatLong_2,FullTextSearch  FullTextSearch  0       1   Using where     1   SIMPLE  TB  eq_ref  PRIMARY     PRIMARY     752     gensurv_isikota.TA.BusinessID   1   Using index      It basically says only one row is computed for each of the process. number of data is 1.6 million, but number of data in a box is only 30


title: How to create a regular polygon grid in Openlayers?

text: I am trying to create a grid in my map and do some calculations. **I seen the following links:-** Create regular polygon grid using postgis Create regular polygon grid inside a polygon using postgis Here I understood, how to create a polygon grid using postgis spatial query. Any one, please help me to do the same concept in Openlayers(any client side code). My arguments are :   1. **Bounds value.** Example: new OpenLayers.Bounds(77.06173,28.50872,77.40918,28.74567)   2. **Distance between the points in the grid**. Example : 100 **My actual requirements is**   1. Create a grid on map.   2. Plot some POINTS.   3. Calculate the number of points in each grid.   4. Get the centroid(x,y) of the grid which have maximum number of points. I just seen this code..               var points = new OpenLayers.Layer.PointGrid({         isBaseLayer: true, dx: 15, dy: 15     });          var map = new OpenLayers.Map({         div: "map",         layers: [points],         center: new OpenLayers.LonLat(0, 0),         zoom: 2     });      Reference: Openlayers Point Grid **My Queries**   1. How can I add polygon grid layer.   2. How can I set this layer in a particular bounds/extent.   3. How can I give "select" control to this layer.   4. How can I find out Number of points on a selected grid/polygon. Note: I am using Openlayers 2.11.


title: Graduated symbols not appearing in legend in QGIS print composer

text: When I add a legend in QGIS version 2.2 print composer, the symbols of a layer containing graduated symbols does not appear. Does anyone know how to make the symbols visible? Thanks.


title: How to build a Spatialite road network from OSM data

text: I'm looking for algorithm which could help me with building my own road network (for routing) from OSM data. _I chose Spatialite DB for this, but here I have only one problem: Spatialite provides interface for build road network, but prepared road graph is saved in BLOBS so I can't get any information about it (of course, I can write queries to DB for routing facilities) but I want to draw all links on my map, so I need to get their coordinates._ Actually, I'm not sure if Spatialite could provide such API or smith like that. My question is how can I build my own road graph (network) from raw OSM data (I've already parsed XML file for highways=*) ?


title: How to shrink polygon on two sides?

text: I have a polygon which I know is a rectangle, although it may be at any orientation: ![Example polygon](http://i.stack.imgur.com/vG265.png) I want to shrink this polygon by a third of the width from two sides, so that I just get the strip in the middle - ending up with a result like the red area in this image: ![enter image description here](http://i.stack.imgur.com/Hyehc.png) I can't think how to go about doing this easily. If I try and negatively buffer the polygon then it will be shrunk on all sides. I have many tens of thousands of these polygons to do - so obviously can't do them manually. Any ideas? I have access to ArcGIS, QGIS and am relatively experienced at Python programming - if that helps.


title: R, rasterVis and RGL: how to prevent reflective surfaces?

text: I have been experimenting with digital terrain model files in ESRI ASCII grid raster format. Despite zero experience with this type of data I found it very easy to load into R, as per code below. The `plot3D` function from the `rasterVis` package defaults to a nice looking plot but the model is rather reflective. Given that it's supposed to be a landscape this shiny surface is not really appropriate. EDIT: I am looking for a more natural (if such a word can be used about a 3D rendering of a set of figures), matt, topological plot. ![shiny](http://i.stack.imgur.com/DgGcT.png) I imagine there must be ways around this but I have no experience using `raster` and would much appreciate some pointers. So my question is: how can I reduce or remove the reflectivity of the plot?               library(raster)     library(rgdal)     library(rasterVis)          foo <- raster(readGDAL("my.dtm.asc"))     plot3D(foo)


title: Why using multiversioned views?

text: I don't see the reason why multiversioned views exist. It's not possible to use joins like in normal views. And editing versions is possible by directly using the feature class itself. So when am I supposed to use them? * * * So if I understand: ![enter image description here](http://i.stack.imgur.com/FRu97.png) When I query the feature class, I query the compressed DB. And **_the only_** way to access data in the versions, is using mv views. Am I wrong?


title: More than one group for a single layer in Geoserver?

text: I have a huge number of layers in my geoserver that I have painstakingly imported and they all work correctly when viewed with GeoExplorer etc. I want to know if there is any way to create set maps which show groups of these layers. I know how to group the layers together but don't know if i can use the same layer in multiple groups. If this is entirely impossible, what is the best way to create set maps- using a new workspace/store for the layers that require duplicating? The maps will eventually be linked into a server that hopefully will have different interactive maps created from the above. Hope this makes sense, inexperienced and very much winging what I'm doing! Any help appreciated!


title: ArcMap Calculate Field..if "text" in !Field!

text: Is it possible to use Field Calculator in ArcGIS to add text to a new field based on an if in statement? Example: In an attribute table I have a some fields called !X!, !Y!, !Z! where "A (xx%)", "B (xx%)", "C (xx%)", "D (xx%)" are the data. The data can be found in any field but does not have to. I would like to use python in some form so that once I have created a New Field, I can calculate `if A in !X!` **copy** A in !X! to New Field then continue this through !Y! and !Z!. I believe outside of Field Calculator this will involve If in statements followed by arcpy.da.insertcursor() but I'm pretty stumped at the moment. I hope this is clear...thanks for any help and comments!


title: Raster format similar to ECW

text: I need an alternative to the ECW format that we currently have for a new series of Ordnance Survey maps. The ECW format looks clearer and less pixelated than the new maps in Tif format (in QGIS). Are there any alternatives formats or methods I could use to improve the appearance of the Tifs? As you can see the two map series look odd next to one another (ECW at the bottom): ![enter image description here](http://i.stack.imgur.com/QYR8Q.jpg) Thanks for your time.


title: TileMill-style based on the country the objects belong to

text: Is it possible to have different TileMill stylings in different countries? I need for instance different colours for highways based on the country the highway is in. Is there a best practice solution to this problem? I know it's possible to store the borders of all countries in PostGIS an check the intersection of an object with all countries. I could also precalulate the corresponding country for all objects like in Nominatim. Are there any better solutions?


title: Determine 'native' size of wms layer based on GetCapabilities response

text: I'm trying to understand the relation between a `<LatLongBoundingBox>` and `<BoundingBox>` in a `GetCapabilities` WMS response, given that the SRS of the bounding box is EPSG:4326. Here's an example of the WMS that I'm trying to query:               <Layer>     <!-- snip -->     <LatLonBoundingBox maxx="37.968" maxy="70.144" minx="-4.646" miny="48.833"/>     <BoundingBox SRS="EPSG:4326" maxx="0.662672" maxy="1.224240" minx="-0.081091" miny="0.852295"/>     <!-- other BoundingBoxes omitted -->     </Layer>      My questions related to this response:   * Why are the `BoundingBox` values different from the `LatLonBoundingBox` values an how is this useful? It looks like a scaling hint of some sorts.   * Is it possible to determine the 'native' resolution of the layer (in pixels) using these values?


title: Setting up Pyscripter for QGIS dev

text: I have been trying to set up PyScripter so that I can make use of code completion for QGIS python code. So far, I have been unable to get things to work.   1. PyScripter when set up as per instructions in the python cookbook is not starting. I have set up a batch file with the following contents to launch PyScripter.              SET PYTHONHOME=C:\OSGeo4W\dev\bin     PyScripter --python27 --pythonpathdll C:\OSGeo4W\dev\bin %1 %2 %3 %4 %5     PAUSE        2. The next issue is registering QGIS core in PyScripter. In the PyScripter help file this is done by just adding the names of the libraries into the Tools\Options\IDE Options..Special packages area. Some special packages such as wx and scipy are in there already but I suspect that qgis needs path specified somewhere so that correct qgis can be found. Just adding qgis to the list of special packages and then placing `from qgis import core` into a py file does not result in qgis classes and method becoming available in the code explorer.  Regarding point 1, PyScripter starts and displays a list of valid command line arguments and then quits without loading. PyScripter was installed via an installer (32 bit) and I have also tried the roaming installation based on the zip file version of PyScripter. My goal is to get PyScripter using the same python that QGIS dev is using. I have no idea why the command line options fail. I am using Win7 x64 and my QGIS install is in C:\OSGeo4W\dev. I suspect that version creep has made the existing instructions obsolete. Any help in getting PyScripter able use qgis' python and to recognise qgis' classes would be very much appreciated. Any set up instructions for an alternative IDE would also be welcome. AndrewM


title: Set custom extent in ArcGIS JSAPI 3.4/3.5

text: I am using ESRI suit 10.1 and ArcGIS JSAPI 3.4/3.5. I want to set custom extent (for India) for this map and trying to assign the coordinate in center but its not getting proper map. In previous version (3.0/3.1) we used `dojo.toJson(map.extent.toJson());` with the help of firebug to identify the extent of the specific area(x min,y min, x max,y max) then we used that extent for base map. What I tried : -I checked the API reference but I did not get any information.  -Is there any Dojo problem (might be few references gets updated?) -Tried to use lat,lon values at Center property but not getting perfect map for India In new version I am not getting this so is there any way to identify the custom extent using previous tool and use those in center property ? Any help will be great !!!


title: Integer field name changes automatically

text: Last week I had imported some shapfiles into file geodatabase. It makes sense but when I checked table of the outputs I had found that some integer field names had been changed to other names. For example original name is 2010 changed to F1. Please help.


title: Google Maps: obtaining the border of an area

text: When using the geocode API it is easy to obtain the (centre of the) coordinates of the address. For example: Rotterdam. If you look up the same location in Google Maps you also get the surface area of Rotterdam. I would like to do the same using the API, is this possible? If yes, how?


title: Where can I get neighborhood boundaries?

text: Besides Zillow's creative-commons licensed neighborhood boundaries shapefiles, what sources of neighborhood boundaries are there? Not just in the US, but all over the world. Free/open source is ideal, but paid sources could be interesting too.


title: How to migrate Opengeo Suite from localhost to live website

text: I have installed Opengeo Suite on my laptop. I was trying to follow the OpenGeo Suite user manual: http://suite.opengeo.org/opengeo- docs/webapps/index.html#webapps I implemented the sections "Creating and deploying apps with OpenGeo Webapp SDK" and "Application development with the Webapp SDK". Now I can upload maps and visualize them in my laptop through "localhost " I do not have any other computer except my personal laptop. I want to migrate the whole system from my laptop to a live website. I am beginner in web world. I searched in the internet and found the following free cloud MEGA: https://mega.co.nz/ I copied the 'OpenGeo Suite' folder and upload it on MEGA cloud. I also have a domain called 'www.worldofgis.net'. But I do not know how to make connection between my cloud storage and domain site. I would really appreciate step by step guidance from experts in this matter. Here is a snapshot from my laptop, ![enter image description here](http://i.stack.imgur.com/7Cby7.jpg)