title: Problem with WFS filters (BBOX, INTERSECTS)

text: I'm trying to get a feature/features from a small area using the BBOX (or INTERSECTS) filter in my WFS query. I should be expecting only a small set of features or a single feature and yet the query returns quite a handsome bunch of them. I have GeoServer version 2.2.2 with an Oracle data store. I have disabled the 'loose bbox' selection so that shouldn't be the cause of the problem. Here's the query with a BBOX filter:               <wfs:GetFeature     xmlns:wfs="http://www.opengis.net/wfs"     service="WFS"     version="1.1.0"     outputFormat="json"     xsi:schemaLocation="http://www.opengis.net/wfs     http://schemas.opengis.net/wfs/1.1.0/wfs.xsd"     xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance">       <wfs:Query typeName="LiVi:LIIKENNE_ELEMENTTI" srsName="EPSG:3067" xmlns:LiVi="http://172.17.14.211:8080/LiVi">         <ogc:Filter xmlns:ogc="http://www.opengis.net/ogc">           <ogc:BBOX>             <ogc:PropertyName>GEOMETRY</ogc:PropertyName>             <gml:Envelope xmlns:gml="http://www.opengis.net/gml" srsName="EPSG:3067">               <gml:lowerCorner>316600 6838986</gml:lowerCorner>               <gml:upperCorner>327696 6844298</gml:upperCorner>             </gml:Envelope>           </ogc:BBOX>         </ogc:Filter>       </wfs:Query>     </wfs:GetFeature>      The problem persists with the INTERSECTS filter. The size of the polygon is about 10*10 meters.               <wfs:GetFeature     xmlns:wfs="http://www.opengis.net/wfs"     service="WFS"     version="1.1.0"     outputFormat="json"     xsi:schemaLocation="http://www.opengis.net/wfs     http://schemas.opengis.net/wfs/1.1.0/wfs.xsd"     xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance">       <wfs:Query typeName="LiVi:LIIKENNE_ELEMENTTI" srsName="EPSG:3067" xmlns:LiVi="http://172.17.14.211:8080/LiVi">         <ogc:Filter xmlns:ogc="http://www.opengis.net/ogc">           <ogc:Intersects>             <ogc:PropertyName>GEOMETRY</ogc:PropertyName>             <gml:Polygon xmlns:gml="http://www.opengis.net/gml" srsName="EPSG:3067">               <gml:exterior>                 <gml:LinearRing>                   <gml:posList>308082.07106781186 6833724.928932188 308082.07106781186 6833739.071067812 308067.92893218814 6833739.071067812 308067.92893218814 6833724.928932188 308082.07106781186 6833724.928932188</gml:posList>                 </gml:LinearRing>               </gml:exterior>             </gml:Polygon>           </ogc:Intersects>         </ogc:Filter>       </wfs:Query>     </wfs:GetFeature>      **Edit:** Here's another request (this time a GET request):               http://172.17.14.211:8080/geoserver/wfs?request=GetFeature&version=1.1.0&srsName=EPSG:3067&typeName=LiVi:LIIKENNE_ELEMENTTI&BBOX=316600,6838986,316605,6838991,EPSG:3067      Even though the BBOX area is quite small, I get 1263 features back, most of which don't even fall inside the bounds specified in the request. For instance, from one of the features in the response:               <gml:boundedBy>       <gml:Envelope srsDimension="2" srsName="http://www.opengis.net/gml/srs/epsg.xml#3067">         <gml:lowerCorner>317629.69841038 6841957.478078741</gml:lowerCorner>         <gml:upperCorner>317902.64972173725 6841987.000123474</gml:upperCorner>       </gml:Envelope>     </gml:boundedBy>


title: How to open AWS ports for inbound traffic

text: I am configuring an instance of Geoportal Server on an Amazon EC2 instance. I have ports 8080 open for tomcat, but need to open ports for Apache Directory Studio. How do I view ports which are accessible for inbound traffic?


title: QGIS Print Composer with ESPG:32615

text: I have imported a PostGIS layer as the only object into my project and set the project and layer CRS to EPSG:32615 (UTM Zone 15N). My layer appears just fine in the work area. When I try to create a print composer, I cannot get the map to appear in the print composer, even though I click the `Set to map canvas extent` button in the composer. The print composer only shows a blank, white area. If I set the project CRS to EPSG:4326 I can add the map to the print composer. Is there a setting I need to change to adjust to CRS of the print composer? **EDIT** Following some advice, I created a layer and manually added a polygon. I noticed the issue pops up when I try to place the polygon near the origin. If the polygon is hundreds of thousands of meters from the origin, when I add the map to the print composer, things work fine. If I place the polygon around the origin (positive and negative x and y values), the print composer seems not to like it. With the linked file, I noticed that if I changed the scale of `Print Composer 3` to 32 or less a zoomed in version of the polygon appears. Anything higher than 32 gives me a blank page. Test File: https://drive.google.com/folderview?id=0B-SpdXvKI1JgX1lGbnFDMWFxRHM&usp=sharing This makes me think I've configured something incorrectly. Can anyone spot it? Sample WKT:               "MULTIPOLYGON Z (((11.0959798314622 6.46039688549521 0,11.885913117052 5.85981372972084 0,11.6611675157088 5.56421098165167 0,10.8708850668735 6.16452289633918 0,11.0959798314622 6.46039688549521 0)))"     "MULTIPOLYGON Z (((11.885913117052 5.85981372972084 0,11.0959798314622 6.46039688549521 0,11.3224760481292 6.75811300197837 0,12.1123181897999 6.15759914256564 0,11.885913117052 5.85981372972084 0)))"     "MULTIPOLYGON Z (((11.5466005766636 7.05271167221563 0,12.3363525288016 6.45226638354819 0,12.1123181897999 6.15759914256564 0,11.3224760481292 6.75811300197837 0,11.5466005766636 7.05271167221563 0)))"     "MULTIPOLYGON Z (((12.3363525288016 6.45226638354819 0,11.5466005766636 7.05271167221563 0,11.7764387175175 7.35482055386346 0,12.5660981809189 6.75444558401706 0,12.3363525288016 6.45226638354819 0)))"     "MULTIPOLYGON Z (((12.0073425546073 7.65833023129213 0,12.7984300880318 7.06002642651748 0,12.5660981809189 6.75444558401706 0,11.7764387175175 7.35482055386346 0,12.0073425546073 7.65833023129213 0)))"     "MULTIPOLYGON Z (((12.7984300880318 7.06002642651748 0,12.0073425546073 7.65833023129213 0,12.232087222072 7.95374403812074 0,13.0236497939066 7.35625275343114 0,12.7984300880318 7.06002642651748 0)))"     "MULTIPOLYGON Z (((12.4618565878502 8.25576251892767 0,13.2512402335275 7.65559725220707 0,13.0236497939066 7.35625275343114 0,12.232087222072 7.95374403812074 0,12.4618565878502 8.25576251892767 0)))"     "MULTIPOLYGON Z (((13.2512402335275 7.65559725220707 0,12.4618565878502 8.25576251892767 0,12.6819923681615 8.54511821138774 0,13.470326977221 7.94375703426402 0,13.2512402335275 7.65559725220707 0)))"     "MULTIPOLYGON Z (((12.9075850363348 8.84164744281931 0,13.6967905731785 8.24161942134774 0,13.470326977221 7.94375703426402 0,12.6819923681615 8.54511821138774 0,12.9075850363348 8.84164744281931 0)))"     "MULTIPOLYGON Z (((14.157865342396 8.84806051755353 0,13.9250314628096 8.5418194421421 0,13.1383077374417 9.14491982213477 0,13.3700411704104 9.44952075323603 0,14.157865342396 8.84806051755353 0)))"      Project CRS: ![enter image description here](http://i.stack.imgur.com/uaSQt.png) Layer CRS: ![enter image description here](http://i.stack.imgur.com/kP8yR.png)


title: Is it possible to export MapJ object into TAB file in MapXtreme?

text: I am trying to open a TAB file having geometry information and bind/join it with JDBC Table having attribute information. Now I have to export the joined table along with map/geometry information from source TAB file and export to new TAB file.


title: ESRI Geoportal Linux install error

text: I'm following the How to Set Up Esri Geoportal Server on Linux white paper. Going everything OK until I have to run ./grants_linuxpg.sh (page 10) `./grants_linuxpg.sh: 20: ./grants_linuxpg.sh: Syntax error: "(" unexpected` Did someone already had this error? What do I wrong? I just can't believe that there is a syntax failure in the original script downloaded from the geoportal sourceforge website?! Thanks for helping!


title: How to get screen coordinate from geographic one?

text: I need to know, how can I convert some geographic coordinates to the screen one. I want to remind you about desktop coordinates system, how it look like: ![enter image description here](http://i.stack.imgur.com/ftgP8.png) I've developed soft, which is loading map fragments into matrix ( 10x10 ) in the fullscreen mode. The map images are loading perfectly. But I have a task to draw some real elements on map, which do have longitude, latitude values. So... I first of all think to convert geographic data to the web mecator values by using such a formulas:                   internal Models.PointGeospatial ConvertFromMercatorToSpherical(double mercatorX, double mercatorY)         {             if (Math.Abs(mercatorX) < 180 && Math.Abs(mercatorY) < 90) return null;             if ((Math.Abs(mercatorX) > 20037508.3427892) || (Math.Abs(mercatorY) > 20037508.3427892)) return null;                  double subitemFirst = mercatorX / 6378137.0;             double subitemSecond = subitemFirst * 57.295779513082323;             double subitemThird = Math.Floor((double)((subitemSecond + 180.0) / 360.0));             double subitemFourth = subitemSecond - (subitemThird * 360.0);             double subitemFifth = 1.5707963267948966 - (2.0 * Math.Atan(Math.Exp((-1.0 * mercatorY) / 6378137.0)));                  var model = new Models.PointGeospatial()             {                 Longitude = subitemFourth,                 Latitude = subitemFifth * 57.295779513082323             };                  return model;         }              internal Models.PointMercator ConvertFromSphericalToMercator(double longitude, double latitude)         {             if ((Math.Abs(longitude) > 180 || Math.Abs(latitude) > 90)) return null;                  double subitemFirst = longitude * 0.017453292519943295;             double x = 6378137.0 * subitemFirst;             double subitemSecond = latitude * 0.017453292519943295;                  var model = new Models.PointMercator()             {                 X = 6378137.0 * subitemFirst,                 Y = 3189068.5 * Math.Log((1.0 + Math.Sin(subitemSecond)) / (1.0 - Math.Sin(subitemSecond)))             };                  return model;         }      Then I'm getting such values after conversion: > * * * > > Building Id: 302 > > Geographic data: { Longitude: 37.5303305475 ; Latitude: 55.7481343179 } > > Mercator data: { X: 4177857.2858509403; Y: 7508438.839224923 } Really... I don't even imagine how could I place correctly the item on map just only knowing its longitude, latitude values and its mercator equivalent with having 640x480 opened window for e.g., on which I must place some pushpins... The one detail I've foud, that first three digits from mercator values are not changing for the one city (and I understand they must not ), so the task to set a correct place for the item in the map with exact zoom level is rather diffcult for me and I don't understand how to solve this problem. I'm listening for you suggestions, thanks!


title: Garmin GPS 78 Interface

text: Does anyone know how to interface a Garmin GPS 78 Marine GPS with their computer to manage waypoints (upload, download, save, etc.) Previously I was using the MN DNR Garmin tool, but it no longer recognizes the GPS. Any solutions are welcome Thanks


title: Clip GAPAnalysis and save as a new file output strange looking colors

text: I download USGS GAP Analysis and use a polygon to clip raster from that GAP Analysis. I am wondering why there is a strange looking color when you clipped the raster. If you zoom in, you can see myraid colors, but when zoom out you see some white spotty areas. I can't insert 2 saved image for you to review it, but here is the link you can see it . Gap Analysis example


title: postgresql 8.4 and postgis 1.5 error: type modifier is not allowed for type 'geometry'

text: I would like to generate a buffer (within a limit of 3km ) where the user clicked point. In pgAdmin, select my database and open a query editor and create 2 PostGIS functions (UTM Zone Function and ST_Buffer_Meters Function) Used following SQL command to create another column in to my table. my table is “schools” with lat and lon fields. ALTER TABLE schools ADD COLUMN geom_buffer geometry(POLYGON,4326); but appear error: type modifier is not allowed for type 'geometry' Why is that..??/ Pls. help me anyone. ..


title: Find extent of data inside many polygons

text: This is similar to Calculate polygon area inside extent I have a few different polygon datasets which describe things like soil type, bushfire risk etc. For each property boundary we want to see how much of each soil type etc there is. I tried Tabulate Interesection ![enter image description here](http://i.stack.imgur.com/z2OcD.jpg) but it gives me percentages over 100% ![enter image description here](http://i.stack.imgur.com/sk8xp.jpg) . The 3 in blue are the same property. Any suggestions?


title: Input field not showing up in Distance Matrix In QGIS 2.0

text: I tried to perform a Distance Matrix analysis of building centroids. I have a two building centroid layers. The attributes in both the tables are `name, osm_id, perimeter and area`. In the Distance Matrix Tool, I select this layer as input field then I have to select the unique field ID. My unique field ID is osm_id, but it does not show up in the list. Does someone know why that could be? I get all the other attributes... Thank you!


title: How to set Window id in Mapbasic?

text: Last couple of days i am working with creating thematic map using MapBasic. I am facing some difficulties while doing scripting and compilation.i am getting error due to Dynamic Window Id. does any code or solution to set window id? current window id i got while doing thematic map manually using mapinfo tools. if i used that window id in to my code it working, but only thing is that i have to open that window. if i closed that window and run my code with that window id it showing me incorrect window id. anybody can tell me..how to solve this issue thanks and regards


title: edit lock layer

text: I have a question. I want to make a protected processing imported layers from the sqlite database. Is there in QGIS 2.0 a function or a plugin or do I have to insert directly into the data? Thank you for your help


title: Overlapping of two layers in MySQL spacial extensions

text: I have been into GIS for almost a week now. Right now, I am learning MySQL spacial extensions. I have a question regarding the functions MySQL provides: I have a table that consists of LineString, say roads, and another table that contains Points, say all the fuel stations. Is there any function that will help me combine the two? Or other way i can think of putting the question is, I have a layer of roads and another of all fuel stations. I want to overlap one layer on other and create a new layer. Is this possible using MySQL? If not, which GIS system should I be looking at to achieve the above task? Thank you. EDIT: I am doing a term paper and one of the tasks assigned to me was to do a feasibility study of of certain paper that was written a few years ago(guessing it to be very old :-P). There the author had claimed to have written an algorithm for overlay operation. However not much details have been provided by the author on what the contents of layers were. So I was wondering if MySQL or any other vendors have provided any features to do the overlap operation now. Went through MySQL spatial extensions document, couldn't find anything there.


title: Is there any statistical method to proof if the method is suitable or not

text: I have recently done a phenological study using threshold method from MODIS NDVI images. But I'm suffered that I don't have any situ observation (ground truth) data to do validation. Is there any statistical method to validate that this method is suitable for remote sensing based phenological study without ground truth data. Phenology is the study of periodic plant life cycle events (onset of green-up, maturity, offset of green-up etc) and how these are influenced by seasonal and interannual variations in climate, as well as habitat factors. Thanks in advance


title: Calculate the values of range, sill, nugget in ordinary kriging?

text: I am working on temperature data, and I want to interpolate the data to any location that I specify. My programming language is asp.net4 C#. I don't know how can I determine the sill, range, nugget for experimental variogram or theoretical variogram. I make the following code to calculate these variables in experimental variogram:               public void Calculate_EmpericalVari_nugget_sill_range()     {     double H=0,VaiogramValue;     int N = KnownPoints.Count;     for (int j = 0; j < N; j++)//calculate the emperical variogram     {     for (int i = 0; i < N; i++)     {     H = AucledianDistance(KnownPoints[i].X, KnownPoints[i].Y,      KnownPoints[j].X,     KnownPoints[j].Y);     VaiogramValue = 0.5 * Math.Pow((KnownPoints[i].Z - KnownPoints[j].Z), 2);     Emper_Variogram[i][j] = VaiogramValue;     if (H==0) nugget = VaiogramValue;     if (sill < VaiogramValue)     {     sill = VaiogramValue;     range = H;     }     }     }         }      I work according to the following that sill is the maximum value and range is the distance when I get the sill and nugget is the value at distance equal zero? Is this right? And do they remain constants when I calculate the theoretical semivariogram?


title: why is there a difference between using EPSG and using basic code

text: I have the OSgeo4w installed on my PC and I am trying to upload ITM(israel transverse mercator) landmarks to google earth.   To do this I am using the ` cs2cs` command.   I noticed that i get a deviation using the next code:   ` cs2cs +proj=longlat +ellps=WGS84 +datum=WGS84 +to +proj=tmerc +lat_0=31.73439361111111 +lon_0=35.20451694444445 +k=1.000007 +x_0=219529.584 +y_0=626907.39 +ellps=GRS80 +units=m`   I tried the next code:   `cs2cs +init=epsg:4326 +to +init=EPSG:2039`   this time my deviation decreased significantly.   why does this happen?   Another thing is, when i use the EPSG code I get a third number (that is always zero using the first code).   The third number, after the lon and lat is around 20, it cannot be the elevation so what does it represent?


title: Creating numpy.array with variable number of fields to test arcpy.da.ExtendTable performance?

text: There was an **ArcPy Café** blog posting entitled **Adding Fields: Performance Tips** which advocated: > Two approaches to help increase performance when adding numerous fields to a > table or feature class. ... > `2`. Use the data access and NumPy modules. The data access module function > named ExtendTable() joins the contents of a NumPy structured array to a > table based on a common attribute field.   >  This is the faster approach, however, the types of fields you can add using > numpy are limited. There is no support for adding blobs, raster, and date > fields. In addition, the field alias can not be defined or altered. I have done some performance testing, using **ArcGIS 10.2 for Desktop** , of this method versus multiple (three and nine) Add Fields on file geodatabase feature classes of 16 to 1,000,000 polygons. My initial findings suggest that any gains for nine and especially three fields are lost as the number of polygons increases. Consequently, I suspect that this method should only be advocated for many fields and not so many features. Below is my current test code, and its current output but these are presented purely for interest/background. **What I would like know is how to make the number of fields placed into the NumPy Array become configurable?** I suspect that there is a Python technique unknown to me but my searches have not yet turned it up so I will try to explain ... If I have a variable `numFields = 10` can I use something like `for i in range(numFields)` to expand `'TEST_INTEGER'+str(i)` to become ...               narray = numpy.array([],                          numpy.dtype([('_ID', numpy.int),                                       ('TEST_INTEGER'+str(0), numpy.int),                                       ('TEST_INTEGER'+str(1), numpy.int),                                       ('TEST_INTEGER'+str(2), numpy.int),                                       ('TEST_INTEGER'+str(3), numpy.int),                                       ('TEST_INTEGER'+str(4), numpy.int),                                       ('TEST_INTEGER'+str(5), numpy.int),                                       ('TEST_INTEGER'+str(6), numpy.int),                                       ('TEST_INTEGER'+str(7), numpy.int),                                       ('TEST_INTEGER'+str(8), numpy.int),                                       ('TEST_INTEGER'+str(9), numpy.int),                                       ]))      * * *               import arcpy,numpy,time          fc = r"C:\polygeo\Projects\test.gdb\testFishnet"     fc2 = r"C:\polygeo\Projects\test.gdb\testFishnet2"     cellWidthHeightList = ["0.25","0.1","0.025","0.01","0.0025","0.001"]          for cellWidthHeight in cellWidthHeightList:         numCells = (1 / float(cellWidthHeight)) ** 2         print "Creating fishnet of {0} polygons".format(str(int(numCells)))         if not arcpy.Exists(r"C:\polygeo\Projects\test.gdb"):             arcpy.CreateFileGDB_management(r"C:\polygeo\Projects","test.gdb")         if arcpy.Exists(r"C:\polygeo\Projects\test.gdb\testFishnet"):             arcpy.Delete_management(r"C:\polygeo\Projects\test.gdb\testFishnet")             arcpy.CreateFishnet_management("C:/polygeo/Projects/test.gdb/testFishnet",                                            "0 0","0 1",cellWidthHeight,cellWidthHeight,                                            "#","#","1 1","NO_LABELS","#","POLYGON")         if arcpy.Exists(r"C:\polygeo\Projects\test.gdb\testFishnet2"):             arcpy.Delete_management(r"C:\polygeo\Projects\test.gdb\testFishnet2")             arcpy.Copy_management(r"C:\polygeo\Projects\test.gdb\testFishnet",                                   r"C:\polygeo\Projects\test.gdb\testFishnet2",                                   "FeatureClass")              start = time.clock()         arcpy.AddField_management(fc,"TEST_INTEGER","LONG","#","#","#","#",                                   "NULLABLE","NON_REQUIRED","#")         arcpy.AddField_management(fc,"TEST_TEXT","TEXT","#","#","100","#",                                   "NULLABLE","NON_REQUIRED","#")         arcpy.AddField_management(fc,"TEST_FLOAT","DOUBLE","#","#","#","#",                                   "NULLABLE","NON_REQUIRED","#")         arcpy.AddField_management(fc,"TEST_INTEGER2","LONG","#","#","#","#",                                   "NULLABLE","NON_REQUIRED","#")         arcpy.AddField_management(fc,"TEST_TEXT2","TEXT","#","#","100","#",                                   "NULLABLE","NON_REQUIRED","#")         arcpy.AddField_management(fc,"TEST_FLOAT2","DOUBLE","#","#","#","#",                                   "NULLABLE","NON_REQUIRED","#")         arcpy.AddField_management(fc,"TEST_INTEGER3","LONG","#","#","#","#",                                   "NULLABLE","NON_REQUIRED","#")         arcpy.AddField_management(fc,"TEST_TEXT3","TEXT","#","#","100","#",                                   "NULLABLE","NON_REQUIRED","#")         arcpy.AddField_management(fc,"TEST_FLOAT3","DOUBLE","#","#","#","#",                                   "NULLABLE","NON_REQUIRED","#")         elapsed = (time.clock() - start)         print " - AddField took {0} seconds to add NINE fields".format(elapsed)              start = time.clock()         narray = numpy.array([],                              numpy.dtype([('_ID', numpy.int),                                           ('TEST_INTEGER', numpy.int),                                           ('TEST_TEXT', '|S100'),                                           ('TEST_FLOAT', numpy.float),                                           ('TEST_INTEGER2', numpy.int),                                           ('TEST_TEXT2', '|S100'),                                           ('TEST_FLOAT2', numpy.float),                                           ('TEST_INTEGER3', numpy.int),                                           ('TEST_TEXT3', '|S100'),                                           ('TEST_FLOAT3', numpy.float),                                           ]))         print " - numpy.array() took {0} seconds".format((time.clock() - start))         arcpy.da.ExtendTable(fc2, "OID@", narray, "_ID")         print " - arcpy.da.ExtendTable() took {0} seconds".format((time.clock() - start))         elapsed2 = (time.clock() - start)         print " - NumPy and ExtendTable took {0} seconds to add NINE fields".format(elapsed2)         print " - NumPyExtendTable:AddField ratio = {0}".format(str(float(elapsed2) / float(elapsed)))      * * *               >>> ================================ RESTART ================================     >>>      Creating fishnet of 16 polygons      - AddField took 6.15530941159 seconds to add NINE fields      - numpy.array() took 8.66656530336e-05 seconds      - arcpy.da.ExtendTable() took 0.251475596777 seconds      - NumPy and ExtendTable took 0.25444199483 seconds to add NINE fields      - NumPyExtendTable:AddField ratio = 0.0413369950748     Creating fishnet of 100 polygons      - AddField took 7.01486277938 seconds to add NINE fields      - numpy.array() took 6.49992397754e-05 seconds      - arcpy.da.ExtendTable() took 0.523825072221 seconds      - NumPy and ExtendTable took 0.554671962901 seconds to add NINE fields      - NumPyExtendTable:AddField ratio = 0.0790709640866     Creating fishnet of 1600 polygons      - AddField took 6.93555420404 seconds to add NINE fields      - numpy.array() took 2.96487760387e-05 seconds      - arcpy.da.ExtendTable() took 0.297160559526 seconds      - NumPy and ExtendTable took 0.334817165881 seconds to add NINE fields      - NumPyExtendTable:AddField ratio = 0.0482754738887     Creating fishnet of 10000 polygons      - AddField took 6.60505587654 seconds to add NINE fields      - numpy.array() took 2.77482134763e-05 seconds      - arcpy.da.ExtendTable() took 0.612328569257 seconds      - NumPy and ExtendTable took 0.672550554964 seconds to add NINE fields      - NumPyExtendTable:AddField ratio = 0.101823598094     Creating fishnet of 160000 polygons      - AddField took 9.54403527444 seconds to add NINE fields      - numpy.array() took 3.57305762435e-05 seconds      - arcpy.da.ExtendTable() took 5.36002166641 seconds      - NumPy and ExtendTable took 5.40552797628 seconds to add NINE fields      - NumPyExtendTable:AddField ratio = 0.566377619198     Creating fishnet of 1000000 polygons      - AddField took 28.7716610917 seconds to add NINE fields      - numpy.array() took 3.61106887681e-05 seconds      - arcpy.da.ExtendTable() took 34.4609012468 seconds      - NumPy and ExtendTable took 34.4646480158 seconds to add NINE fields      - NumPyExtendTable:AddField ratio = 1.19786785706     >>>


title: How to use OpenStreetMap in Java?

text: I have a requirement of GIS to be implemented using Java. I have read the developer guide for OpenStreetMap but there is insufficient information for its setup technically. I am completely new to GIS. Any help regarding Technical Setup of OpenStreetMap as well as information about using them will be very much appreciated. Thanks in advance.


title: Translate data from MapServer to Qgis server

text: Might you help me in this question: I've got operational Mapserver with shape- and map-files on it (*.shp *.map). How can I show them on Qgis server using Qgis web-client or Lizmap? This is my sequence of actions:   1. Open Qgis desktop and add new WMS layer from MapServer (in main windowI can see all of my shape-files with styles, described in map-files).   2. Save Qgis project (*.qgs)   3. Install Qgis project to Qgis server (I make copy of it to cgi-bin directory with qgis_mapserv.fcgi script. Altered way is install it by Lizmap plugin)   4. Try to open this project using Qgis desktop (in new project by utilizing new WMS layer), using Qgis web client or Lizmap. And I see no images in each case. I think that problem may be in Qgis server, which can't open data from MapServer. Have any ideas how can it be fixed?


title: GIS QUERY on point and buffer

text: I have a point table and the buffer geometry is created around the point. I have defined the if the point is within the buffer_distance then its value is =1 and if it is outside the buffer its value is =0, and when it is on buffer_polygon then =2. Now when ever i change the location of point then it should check the following conditions and the flag value should change according to the location of point new position. I am changing the position of the point in QGIS. My table structure is like this:-               GID(PK) Flag(INT)   Buffer_Distance THE_Geom    Buffer_GEOM     1   1   100 point(X,Y Location) Buffer Created around the Point      2   2   200 point(X,Y Location) Buffer Created around the Point      3   2   200 point(X,Y Location) Buffer Created around the Point      4   3   300 point(X,Y Location) Buffer Created around the Point      5   3   300 point(X,Y Location) Buffer Created around the Point      6   1   100 point(X,Y Location) Buffer Created around the Point


title: How to upload a shapefile to geoserver using cURL and replace file?

text: In the sequence of this: How to upload a shapefile to geoserver using cURL? which was properly solved, I now need to be able to replace the existing file, cause my process is to be repeated every day. Any clues? Thanks


title: Advanced GMaps directions

text: I try to create an application for taxis. The application will be GMap based. On the map will be marked the region in which the taxi driver is working in. The region will be a polygon designed by the application owner. Also the client charge will be based on kilometer. The end user can define the departure and destination points on the map. In case that both locations are inside the region all are fine. In case the destination point is outside the region, the client must be charged in double rate for the kilometers outside the region. What I like to ask is the following:   * Is there a way to know what is the distance from the point inside the region until the region border?    * And what is the distance from the point outside the region until the region border? To be more specific, I attach an image here: ![enter image description here](http://i.stack.imgur.com/ksI6G.jpg) What I need is to get the distance from departure location to the border of the region (the yellow circle) and from yellow circle to the destination location. I know that this is very complex, but I don't even know if is possible. Can somebody to help me with that problem ?


title: How to redistribute population data using Python in ArcGIS?

text: I am trying to redistribute population data from census centroids among to housing density.Would you please help me to do that? I do have census area centroids (points) and the boundaries belong to each centroid. Point data has population value. Also I have a raster set housing density. I would like to redistribute population data (points) in their own site-boundaries among housing density. Output data needs to be a raster data 100×100 meter pixels. I am familiar with FORTRAN and VB but I am not a programmer. I have access to ArcGIS 10.1 however I have never used python. Your help and advice will be most appreciated. Best regards Behnam


title: How to perform trilateration using 3 lat/lon points without distances?

text: I read this thread Trilateration using 3 latitude and longitude points, and 3 distances explaining about trilateration using 3 latitude and longitude with measured distance. The problem is I want to find out an unknown target location with only 3 knowns latitude and longitude co-ordinates **without knowing each distance point**. For example, I have 3 following longitudes and latitudes: > lat="-6.28530175" lon="106.9004975375"   >  lat="-6.28955287" lon="106.89573839"   >  lat="-6.28388865789474" lon="106.908087643421"   > Thank you very much :) * * * Based on @D.E.Wright answer I am googling to find out more and I get this following forum http://www.coderanch.com/t/453432/Programming/implement-cell- triangulation-mobile-phones, it gives 3 steps to triangulation. The problem I couldn't understand what he means about "Calculate distance from the first tower based on speed which gives a radius value".      And with @Kirk Kuykendall suggestion I decide to put my 3 lat/long to google map and I get this following image : ![Google Maps](http://i.stack.imgur.com/fKtHQ.png) Where :   A, B, C : Cell Towers location   A : lat **-6.28955287** , lon **106.89573839** with cell tower Radius **6000m**   B : lat **-6.28530175** , lon **106.9004975375** with cell tower Radius **6000m** , and   C : lat **-6.28388865789474** , lon **106.908087643421** with cell tower Radius **6000m**      Now, with these data can I doing a triangulation to get my real position? And how is the mathematical method? I've try to modify phyton code to remove the distance variable from @wwnick in this thread Trilateration using 3 latitude and longitude points, and 3 distances but until now I couldn't find the answer.


title: Large scale geoprocessing through ESRI REST interface?

text: A typical limitation of web based processing is that urls are limited to 2 thousand characters for a GET request. From my understanding of how the ESRI Javascript API works, JSONP is utilized to interact with ArcGIS Server Geoprocessing services via the REST API, requiring the use of a GET request? This seems to pose a problem if I'm attempting to send large quantities of user defined runtime data (> 100 coordinate pairs and several other user parameter strings) to be worked upon in a geoprocessing service via the standard ESRI Javascript API. What's the best way to tackle this problem? I'm thinking about writing a proxy service that will receive GET requests, but will issue POST requests to the REST API.


title: Is it possible to display two categorical symbolizes for two different attributes in one layer alongside each other in ArcMAP 10.1?

text: I have a road layer that I would like to display attribute_1 = "speed" on top of the line segment and attribute_2 = "traffic_volume" below the line segment. For each attribute I would like to setup speed bins 0-30, 30-40, 50-60, 60+ and for volume bins: 0-200, 200-1000, 1000-3000, 5000-10000, 10000+ As a rough Schematic, red represents speed, and black represents volumes and the road should run between those two symbols: ![Displaying the two symbologies for one road layer](http://i.stack.imgur.com/CMDwC.jpg) **Question: I know how to display each attribute as symbol; however, I'm not sure what's the best way to display both symbologies so that they don't overlap** Note that I'm open to using ArcPY


title: Is postgres 9.1.9 supported by ArcGis 10.1 Service Pack 1?

text: I want to install Postgres and Postgis. Now we use ArcGis 10.1 So I install Postgres 9.0.5 But there is not exist Postgis for this version. I try install postgis 2.0.3 But I fail(it is only 32 bit application ) http://download.osgeo.org/postgis/windows/pg90/ So I'd like to install postgres 9.1.3 But I find only 9.1.9 version :( I'd like to install Postgis because we will use geoserver.


title: Uninstall ERDAS ECW/JP2 SDK Version 5.0 (Ubuntu)

text: I am using Ubuntu Linux as main OS when it comes to geoprocessing but I am not yet an expert using it. I recently needed QGIS to have ECW support, so I followed the instructions of Luke in the following question and got it working. ECW for QGIS 1.8 in Ubuntu 12.04 Now I discovered that I keep getting following error message when doing python geoprocessing which involves any GDAL features ERROR 1: /usr/lib/gdalplugins/1.10/gdal_ECW_JP2ECW1.so: undefined symbol: _ZTI14GDALPamDataset It doesn't affect the outcome of the script or crashes it, but it is annoying to get the terminal spammed with this error message during the usage of the script. Is there a solution for that? If not, could anybody guide me through the uninstall process? Note to myself: temporary stuff always in a VM ;)


title: Service published but cache generation failed on ArcGIS Server

text: I am trying to publish a map service containing a cached raster layer on our ArcGIS Server Enterpise 10.1. I analyze the map service before publishing and no error. The error I get when publishing is: > The service has been published successfully. The cache for the service has > been enabled but the cache generation failed." See attached image ![enter image description here](http://i.stack.imgur.com/UdFM6.png) Further details:   * the GIS Server and the folder are on the same machine, with Windows OS.    * we're talking about a virtual machine in the cloud.    * the raster layer is created from a mosaic raster in .tif format (size= 1,5 GB).   * I use ArcGIS Desktop 10.1   * I am trying to create a 5 scale cache.  PS: Found the error on esri forum, but it's for ArcGIS Online servers: http://forums.arcgis.com/threads/60488-Error-Publishing-a-Map-Service-for- ArcMap-to-ArcGIS-Online Thank you in advance for your help.


title: What is the best way to host a GeoDatabase on an ArcGISServer?

text: I have a simple Geodatabase in access format. I'd really like to host this database on an ArcGIS server through one of their services. The idea is that my users would be able to connect to this Geodatabase to look at both their features, and their databases, and then relate them. I do not want to allow editing of the database, just viewing purposes. If at all possible, I'd really like to stay away from using ArcGIS Engine Runtime on the client side. What might be the best service I can run? And does anyone have a good example of how to connect to that service to extract the database information? Thank you!


title: How can I create an adapter between the ArcGIS REST interface and WMS or TMS?

text: I'm trying to use imagery available via the ArcGIS REST tile cache interface, however the tool I'm using, JOSM, only supports WMS or TMS. In addition the data from the ArcGIS server is in a "special" projection, Virginia Lambert Conformal Conic, which doesn't have an ESRI or EPSG wkid but does have WKT provided by the server. Is there a way using free or open source tools to convert the tiles available from the server to WMS/TMS? I know this would involve downloading multiple tiles, then mosaicking and reprojecting them to create tiles suitable for WMS/TMS. I've come across MapServer and GeoServer, which seem very capable, but I wouldn't know where to begin, so any pointers would be helpful. The server I'm interested in is here: http://gismaps.virginia.gov/arcgis2/rest/services/VBMP2006_2007/MapServer


title: 'Import arcpy' error (RuntimeError: NotInitialized) when using IDLE

text: ArcGIS Version: ArcGIS 10, Python 2.6.5 Machine: Windows 7 64bit (Intel) I'm receiving the following error when importing arcpy using the the IDLE (gui and cmd prompt) and when using pythonwin. I had removed python26 and python25 from my system and have done a fresh install of ArcGIS 10 and python. Therefore am calling IDLE from here: "C:\Python26\ArcGIS10.0\pythonw.exe" The Error: **IDLE 2.6.5 > > > import arcpy Traceback (most recent call last): File "", line 1, in > import arcpy File "C:\Program Files > (x86)\ArcGIS\Desktop10.0\arcpy\arcpy__init__.py", line 17, in from > geoprocessing import gp File "C:\Program Files > (x86)\ArcGIS\Desktop10.0\arcpy\arcpy\geoprocessing__init__.py", line 14, in > from _base import * File "C:\Program Files > (x86)\ArcGIS\Desktop10.0\arcpy\arcpy\geoprocessing_base.py", line 568, in > env = GPEnvironments(gp) File "C:\Program Files > (x86)\ArcGIS\Desktop10.0\arcpy\arcpy\geoprocessing_base.py", line 565, in > GPEnvironments return GPEnvironment(geoprocessor) File "C:\Program Files > (x86)\ArcGIS\Desktop10.0\arcpy\arcpy\geoprocessing_base.py", line 521, in > **init** self._refresh() File "C:\Program Files > (x86)\ArcGIS\Desktop10.0\arcpy\arcpy\geoprocessing_base.py", line 523, in > _refresh envset = (set(env for env in self._gp.listEnvironments())) > RuntimeError: NotInitialized** EsriUK and Esri tech support have recommended fresh reinstall and registry edits = but this has not worked. I'm waiting on them for a solution. I've investigated whether it was a module arcpy path issue. There are a couple of links on the Esri forums and on the ArcGIS help regarding paths. Also, a possible related thread on gis.stackexchange. I followed the advise from 'kimo' on the Esri forum thread. But still the import error.   1. My system variable (that I've edited) is: PYTHONPATH = "c:\python26";"C:\Python26\ArcGIS10.0";"C:\Python26\ArcGIS10.0\Lib\site-packages";"C:\Program Files (x86)\ArcGIS\Desktop10.0\arcpy"   2. I have checked I have the contents of "C:\Python26\ArcGIS10.0\Lib\site-packages\Desktop10.pth". The contents are: C:\Program Files (x86)\ArcGIS\Desktop10.0\bin C:\Program Files (x86)\ArcGIS\Desktop10.0\arcpy C:\Program Files (x86)\ArcGIS\Desktop10.0\ArcToolbox\Scripts However, the traceback message shows that a .py script is found in the folder: C:\Program Files (x86)\ArcGIS\Desktop10.0\arcpy\arcpy. Therefore, suggests the arcpy module is found by the interpreter but then not initialized(?). Arcpy can be called from the Python Window within ArcMap and ArcCatalog okay. It is only when called from a Python interpreter that the 'import arcpy' fails. Therefore, arcpy is initalized okay by ArcMap/ArcCatalog, but not IDLE - how can this be? I previously had an install of python25, then python26 then after that upgraded to ArcGIS 10. Therefore, there was quite a lot of legacy python on my machine. However, I've uninstalled all and registry cleaned using ccleaner. Can anyone suggest another line of enquiry as I'm not experienced with python enough to solve initialize/runtime errors? I'd like to avoid reinstalling the OS too! Cheers Sam


title: Remote sensing software options for creating structure footprints from aerials

text: As shown in screen shot here, just wondering what other software options for creating polygons based on structure footprints are available. Something cheaper than ENVI but still versatile would be ideal. Typically need a swath of about 1000 feet on either side of polyline ranging from 1 to 1000 miles. It has been about a year since a similar question was asked and thought there might be some new technology around now. same screen shot as in above link attached here: ![enter image description here](http://i.stack.imgur.com/188CH.jpg)


title: Why to use Add Attachments (Data Management) and not Hyperlink field?

text: Why should I use Add Attachments and not a hyperlinked field? What are the advantages in case of a lot of layers with a lot of attachments at the GDB? What are the best way to manage attachments for a large database that contain a lot pf layers with attachments (JPG, PDFs, DWGs) in term of the paths (pre stage before publishing to ArcGIS Server)?


title: How to convert a track of points into walking directions?

text: Imagine you have a list of points (latitude + longitude) you have previously recorded with a GPS receiver. Now, you need to convert it to a textual representation of the route you followed. Example:   1. Start at Redmond Avenue 12   2. Walk 200 m   3. Turn left to Another St.   4. Walk 100 m   5. Stop at Another St. 10   6. End It doesn't really need to be very detailed. The first idea that comes to my mind is to use an external service (google, yahoo...) to reverse geocode each point. Then apply a simple algorithm to write the route when I detect a change in the street name. Any tip about what is the best algorithm to use or any existing (open) implementation is welcome. It does not seem to be trivial.


title: Best open source operating system for open source GIS?

text: I'm wanting to get into open source GIS, and it looks like PostGIS on Postgres is a good DBMS to go with, and there are plenty of posts on desktop editors. What I haven't come across (and forgive me if I just didn't find it) is info on what open source _operating systems_ folks are using with open source GIS. Ubuntu, some other form of Linux, etc. Thoughts and experiences appreciated.


title: What kind of datasets are necessary for geomarketing applications?

text: If you would have to create a data base for geomarketing applications yourself, what kind of datasets would you purchase/collect? The first goal would be to find new store locations. Another would be to assess the efficiency of advertising efforts in different areas.


title: Master Thesis- topic?

text: It comes the time for me to decide what will be the topic of my Master thesis. I want to write something related to GIS. Maybe you can give me some interested topics (My idea was to do a basic geoportal, but now i think i can't manage to do that). Please write your suggestions. THANK'S A LOT !!!


title: Using Python and open source modules, how can I perform Raster Calculator functions?

text: I would like to perform mathematical functions similar to those you would use Raster Calculator for in ArcMap, however I'd like to do it in Python without using the Arcpy module. I have the osgeo module, and it seems like it should be possible to do however I'm having a hard time finding any useful syntax. Any help would be appreciated. An example function would be to multiply every cell value in an input raster by 2. Thanks, -Scott


title: Is there a source for leaf-off aerial imagery?

text: I am looking for leaf-off imagery for Oklahoma, Kansas and Texas. 1m leaf-off imagery would be ideal, although any substitute would be helpful. I need fairly high resolution imagery to distinguish conifer from deciduous trees-- Landsat would likely be too coarse. How can I obtain leaf-off imagery for KS, OK, TX?


title: Multiple labels in QGIS

text: Is it possible to show 2 or more labels for 1 layer? For example - 1 label will show the name of a school and the 2nd label with show the number of students? I wanT to avoid having 2 layers - one for each label. Thanks


title: Real time positioning of vector features using database

text: I am a newbie in GIS and want to implement something like a vehicle tracker wherein I will get a series of coordinates in a database. I have no idea how to go about it, I mean what should my process consist of? I have found something about GeoJSON, but since I'm using GWT should I be really using it? I guess I can directly interact with the Java data structure. Apart from that, how to initiate the loop for animation? I have seen rotate and move methods in openlayers, but would they suffice in case of a complex motion like vehicular movement.


title: Keep Getting Black OSM Tiles in Openlayers

text: When trying to view OSM tiles using the external storage example of OpenLayers, I keep getting black tiles when I zoom in and out. I don't know what the issue is but every time I zoom in/out there are black tiles that appear. Please help me.


title: Problems while installing Mapserver on Ubuntu

text: I followed instructions at http://www.mapserver.org/installation/unix.html but it didn't work well. My system is ubuntu 14.04 and I'm a total rookie when it comes to unix/compiling/working with the terminal **First Question:** _Is there a functioning script which can be executed in the terminal?_ ( I don't know c - compiling well and don't know how to get all required packages) All I did so far is set up the Apache2. When I browse localhost with Firefox, I get the "Apache2 Ubuntu Default Page". **Second Question:** I want to make a mapserver and a mapclient as WMS-server which is based on geotiff and shapefiles. _Where can I find information to get things started?_ ( I would take ANY similar solution for any other mapserverclient, as long it is understandable for a 10/10 rookie) cheers


title: Openlayers Tasmania example

text: I'm trying to implement an OpenLayers example   http://openlayers.org/dev/examples/getfeatureinfo-control.html using GeoServer that is installed on the same machine as my browser. The map of Tasmania loads fine as well as the roads and lakes layers, but when I left click on the map to display the layer information, I get the spinning wheel of death. I've downloaded the OpenLayers-2.13/OpenLayers.debug.js file for debugging but not sure where to begin. Any suggestions are appreciated.


title: What could cause issues migrating a Python plugin to another Windows version?

text: I am developing a plugin in my laptop and the plugin runs well. So, I installed the same version in other computer and the plugin doesn't run well. Sometimes I develop the code in laptop, sometimes in the other computer. The same language (python), the same version. The windows is different, one is vista and the other is windows 7. The difference can be here? Can someaone explain me why?


title: Using National Agriculture Statistics Service (NASS) county level crop data to estimate crop acreage in hydrologic regions?

text: National Agriculture Statistics Service (NASS) provides county level crop (Corn, Soybean, Sorghum, etc) data. I want to calculate crop acres on each subwatershed from the county based data. I'm using ArcGIS 9.3. I created a join in ArcGIS with the spreadsheet I acquired from NASS. The county level crop data is a shapefile. Watershed (subwatershed) boundaries is also in a shapefile format. What is the best way to deal with this situation?


title: How to analyze census data and display on map (shape file) in ArcGIS

text: first time on this part of Stack Exchange world. I'm hoping you guys can give me an idea on how to start. I basically have all this shape files representing a country, states of that country or towns etc... For each shape file there is associated data with them. So the country shape file has 100 rows and maps to .csv files with 100 rows, same goes for towns, states etc. Basically it gives me a flexibility at how much level I can experiment with. Now, let's say I chose the states data and I wanted to find whether folks with education are more likely to get jobs than those who don't, the first and easy part was to join the shape file to its .csv file that contains data about employment. After that I simply got confused.. the csv files contains raw data and how to use them confuses me, an example (fake data):               SID STYPE   SNAME        MALEP  FEMALEP EMPLOYED  UMEMPLOYED EDUCATION     S1  STA     Alabama      4345   1003    100000    2000       1000     S2  STA     Alaska       1456   8343    8888      4000     S3  STA     Arizona      90000  1247    1000     S4  STA     Arkansas     18831  1213    0      Uh, so how do I use the data above? I don't think a **select by attributes** will do me any good as it shows just a blue border, and I don't know how I could possible query that? It doesn't give me any specific data, it is just an aggregate. So I think my best option is to use either **categories** or **quantities** symbology. However, I'm not sure how to get meaningful data. If I was to use **quantities** (I think **quantities** is the way to go), do I simply normalize people who are employed against people having education? How would I display the output.. density dots or gradient symbology/color? As I said I just need a good shave, thanks! I will be over in the room if you are feeling helpful.. I have been waiting forever for someone to show up :(.


title: Problem with specific arrivals and departure time in Google transit

text: I am developing a transit app and using Google transit GTFS structure. So for table stops_times.txt i have to specify the arrivals and departure times of bus. But in my case buses doesn't have specific times, they keep on moving all the day on a route (say from 0900 to 1800 hrs). How i can tackle this problem ?


title: Installing GeoNetwork on Debian with Tomcat?

text: How can I install GeoNetwork on a Debian server that already has tomcat7. Could you give me help using the guide, please? I followed the steps but I get the following error: > Error Raised exception while starting appl handler. Skipped. Exception > java.lang.IllegalArgumentException: Repository uses database URL of > 'jdbc:h2:~/.geonetwork/geonetwork;LOCK_TIMEOUT=20000;DB_CLOSE_ON_EXIT=FALSE;MVCC=TRUE' > which does not match current database URL > 'jdbc:h2:~/.geonetwork/geonetwork;LOCK_TIMEOUT=20000;MVCC=TRUE'. Modify the > svn property geonet:dburl on the root of the subversion repository at > /var/lib/tomcat7/webapps/geonetwork/WEB-INF/data/data/metadata_subversion or > specify a different subversion repository > > Message Repository uses database URL of > 'jdbc:h2:~/.geonetwork/geonetwork;LOCK_TIMEOUT=20000;DB_CLOSE_ON_EXIT=FALSE;MVCC=TRUE' > which does not match current database URL > 'jdbc:h2:~/.geonetwork/geonetwork;LOCK_TIMEOUT=20000;MVCC=TRUE'. Modify the > svn property geonet:dburl on the root of the subversion repository at > /var/lib/tomcat7/webapps/geonetwork/WEB-INF/data/data/metadata_subversion or > specify a different subversion repository > > Handler org.fao.geonet.Geonetwork > > Stack java.lang.IllegalArgumentException: Repository uses database URL of > 'jdbc:h2:~/.geonetwork/geonetwork;LOCK_TIMEOUT=20000;DB_CLOSE_ON_EXIT=FALSE;MVCC=TRUE' > which does not match current database URL > 'jdbc:h2:~/.geonetwork/geonetwork;LOCK_TIMEOUT=20000;MVCC=TRUE'. Modify the > svn property geonet:dburl on the root of the subversion repository at > /var/lib/tomcat7/webapps/geonetwork/WEB-INF/data/data/metadata_subversion or > specify a different subversion repository at > org.fao.geonet.kernel.SvnManager.(SvnManager.java:179) at > org.fao.geonet.Geonetwork.start(Geonetwork.java:345) at > jeeves.server.JeevesEngine.initAppHandler(JeevesEngine.java:600) at > jeeves.server.JeevesEngine.init(JeevesEngine.java:206) at > jeeves.server.sources.http.JeevesServlet.init(JeevesServlet.java:91) at > javax.servlet.GenericServlet.init(GenericServlet.java:160) at > org.apache.catalina.core.StandardWrapper.initServlet(StandardWrapper.java:1266) > at > org.apache.catalina.core.StandardWrapper.loadServlet(StandardWrapper.java:1185) > at org.apache.catalina.core.StandardWrapper.load(StandardWrapper.java:1080) > at > org.apache.catalina.core.StandardContext.loadOnStartup(StandardContext.java:5026) > at > org.apache.catalina.core.StandardContext.startInternal(StandardContext.java:5313) > at org.apache.catalina.util.LifecycleBase.start(LifecycleBase.java:150) at > org.apache.catalina.core.ContainerBase.addChildInternal(ContainerBase.java:899) > at org.apache.catalina.core.ContainerBase.addChild(ContainerBase.java:875) > at org.apache.catalina.core.StandardHost.addChild(StandardHost.java:618) at > org.apache.catalina.startup.HostConfig.deployWAR(HostConfig.java:963) at > org.apache.catalina.startup.HostConfig$DeployWar.run(HostConfig.java:1600) > at java.util.concurrent.Executors$RunnableAdapter.call(Executors.java:471) > at java.util.concurrent.FutureTask$Sync.innerRun(FutureTask.java:334) at > java.util.concurrent.FutureTask.run(FutureTask.java:166) at > java.util.concurrent.ThreadPoolExecutor.runWorker(ThreadPoolExecutor.java:1146) > at > java.util.concurrent.ThreadPoolExecutor$Worker.run(ThreadPoolExecutor.java:615) > at java.lang.Thread.run(Thread.java:701)


title: 3D analyst tool - add surface information problem with polygons

text: When I'm using 3D analyst tools - functional surface - add surface information to add a z-value to polygons I have created in Arcmap, error 9999 always occurs. I'm using Lidar-data as input feature class. It works when I try it with points instead of polygons. Sorry if I formulated anything incorrectly. I'm just a GIS-beginner. Does anyone have any idea about what I could be doing wrong?


title: What qualifications should an auditor have if they plan on auditing a GIS?

text: Background summary: Various systems are being audited and GIS is the next to be audited. The initial kick-off for the audit was a meeting between me and a single auditor. Initially, I thought it was going to be the typical security audit that is done for all RDBMS, which is fine. That's good. But then she said they would be executing a series of tests against our data (deletes, updates, inserts). The geodatabase is versioned, so naturally I'm concerned. When I received the meeting minutes from the auditor, I had to re-write it because the vocabulary was all wrong (lots of GIS words used in the wrong way/wrong order). I'm not concerned about failing the audit--I'm worried about the auditors really mucking things up because they have no GIS experience. Does anyone have experience with GIS audits and can provide some guidelines as to what qualifications the auditor (or audit team) should have? Do internal audit departments typically have GIS professionals on staff? I added the esri-geodatabase and arcgis tags to this question because of the versioning--I might consider dropping those once a read answers.


title: How to create a circle with ArcEngine?

text: I want to create a circle with IPoint P and radius r in ArcEngine and then I need to cast it to an IGeometry. Does anyone know how to do this? ESRI ArcEngine 10 VS2010 C#


title: OpenLayers request tiles out of range

text: OpenLayers map is requesting tiles out of range and all the map is pink. If I inspect single tile, the address is next: > /Tiles/15/ **16384** /16384.png But the folder with generated tiles does not contain such values: The range is only: /Tiles/15/ **16099** \- /Tiles/15/ **16119** The tiles are generated by Maperitive. Is it a Maperative problem or I should set additional options to OpenLayers.Map?               map = new OpenLayers.Map("demoMap");     var cLayer2 = new OpenLayers.Layer.OSM("New Layer",      "http://localhost:15157/Map/Tiles/${z}/${x}/${y}.png", { numZoomLevels: 19 });     map.addLayer(cLayer2);     map.zoomToMaxExtent();


title: Determing, whether WMS service supports CQL filtering

text: Is it possible to query, whether WMS service supports CQL filtering (like Geoserver does) or not? I didn't find any hints from GetCapabilities query result.


title: ArcGIS Online Topo Maps

text: ArcView 10.1 does not have this kind of topo map to "Add Data". Where do I find this type of topo map, preferably free. ArcGIS online does not seem to offer this type of topo, unless I am not seeing it on my side. ![enter image description here](http://i.stack.imgur.com/2uxyy.jpg)


title: Identify sample doesn't work

text: I'm trying to follow the Identify snippet from the ESRI ArcEngine SDK http://help.arcgis.com/en/sdk/10.0/arcobjects_net/componenthelp/index.html#/Do_Identify_Snippet/00490000000n000000/ The third line of the tutorial throws a COM exception               public void DoIdentify(ESRI.ArcGIS.Carto.IActiveView activeView, System.Int32 x, System.Int32 y)     {       if(activeView == null) return;       ESRI.ArcGIS.Carto.IMap map = activeView.FocusMap;       ESRI.ArcGIS.CartoUI.IIdentifyDialog identifyDialog = new ESRI.ArcGIS.CartoUI.IdentifyDialogClass();  //THIS LINE THROWS THE EXCEPTION      {"Retrieving the COM class factory for component with CLSID {F8AAB776-8B1E-11D3-9F7A-00C04F6BC886} failed due to the following error: 80040111 ClassFactory cannot supply requested class (Exception from HRESULT: 0x80040111 (CLASS_E_CLASSNOTAVAILABLE))."} I expect to run into some trouble using the snippets, but I can't even call a constructor... : ( Any thoughts are appreciated... specs: ESRI ArcEngine 10 VS 2010 C# **Edit:** Maybe I don't have access to the CartoUI library. I also tried this:               AdvancedDrawingDialogClass t = new ESRI.ArcGIS.CartoUI.AdvancedDrawingDialogClass();       But I got a similar error. **Edit 2** I found out how to do the identification so I'm going to post it for the community.               for (int i = 0; i < mapControl.LayerCount; i++)     {       ILayer layer = mapControl.get_Layer(i);       if (!(layer is IDynamicLayer))       {         ESRI.ArcGIS.Carto.IIdentify identify1 = layer as IIdentify;         IArray array = identify1.Identify(geometry);         if (array != null)         {           object obj = array.get_Element(0);           IFeatureIdentifyObj fobj = obj as IFeatureIdentifyObj;           IRowIdentifyObject irow = fobj as IRowIdentifyObject;           IFeature feature = irow.Row as IFeature;           IGeometry geometry1 = feature.Shape;          }       }     }


title: Setting up PostGIS: Postgres can't see the PostGIS extensions?

text: I've just installed Postgres 9.3.4 and PostGIS 2.1 on OSX Mavericks using the very latest versions of the KyngChaos libraries. This is a new machine, and I deliberately haven't installed Postgres using Homebrew or Macports. However, when I try to create spatial extensions on a database, I can't, because Postgres claims not to be able to see the extension file:               :~ anna$ psql -d land -c "CREATE EXTENSION postgis;"     ERROR:  could not open extension control file      "/usr/local/Cellar/postgresql/9.3.4/share/postgresql/extension/postgis.control":     No such file or directory      It looks as though the file is located somewhere else:               :~ anna$ mdfind postgis.control     /usr/local/pgsql-9.3/share/extension/postgis.control      I'm definitely using the expected version of Postgres, although I don't know it's actually the Kyngchaos version:               :~ anna$ psql --version     psql (PostgreSQL) 9.3.4          :~ anna$ which psql     /usr/local/pgsql-9.3/bin/psql      It looks like perhaps there are multiple versions of Postgres floating around on my system, and I'm using the wrong one?               :~ anna$ mdfind -name "postgres" | grep -G "postgres$"     /usr/local/pgsql-9.3/bin/postgres     /usr/local/var/postgres     /usr/local/Cellar/postgresql/9.3.4/bin/postgres      Can anyone advise me on what to do next? * * * **FURTHER UPDATE** Here's the output of `pg_config`:               BINDIR = /usr/local/pgsql-9.3/bin     DOCDIR = /usr/local/pgsql-9.3/share/doc     INCLUDEDIR = /usr/local/pgsql-9.3/include     PKGINCLUDEDIR = /usr/local/pgsql-9.3/include     INCLUDEDIR-SERVER = /usr/local/pgsql-9.3/include/server     LIBDIR = /usr/local/pgsql-9.3/lib     PKGLIBDIR = /usr/local/pgsql-9.3/lib     LOCALEDIR =     MANDIR = /usr/local/pgsql-9.3/man     SHAREDIR = /usr/local/pgsql-9.3/share     SYSCONFDIR = /usr/local/pgsql-9.3/etc     PGXS = /usr/local/pgsql-9.3/lib/pgxs/src/makefiles/pgxs.mk     CONFIGURE = '--with-openssl' '--with-pam' '--with-krb5' '--with-gssapi' '--with-ldap' '--enable-thread-safety' '--with-bonjour' '--with-python' '--without-perl' '--enable-nls' '--with-libxml' 'CC=clang' 'CXX=clang++' 'CFLAGS=-Os -arch x86_64 -isysroot /Developer/SDKs/MacOSX10.7.sdk -D_FILE_OFFSET_BITS=64' 'LD=clang' 'LDLFAGS=-arch x86_64 -isysroot /Developer/SDKs/MacOSX10.7.sdk'     CC = gcc     CPPFLAGS = -I/usr/include/libxml2     CFLAGS = -Os -arch x86_64 -D_FILE_OFFSET_BITS=64 -Wall -Wmissing-prototypes -Wpointer-arith -Wdeclaration-after-statement -Wendif-labels -fno-strict-aliasing -fwrapv     CFLAGS_SL =     LDFLAGS = -arch x86_64 -Wl,-dead_strip_dylibs     LDFLAGS_SL =     LIBS = -lintl -liconv -lxml2 -lpam -lssl -lcrypto -lgssapi_krb5 -ledit -lz     VERSION = PostgreSQL 9.3.4


title: How to make a copy of a temp raster in QGIS?

text: I have a temp raster (saga raster) file. How can I make a copy of this file (like "save as")? Thank you!


title: WMS Image Services Do Not Display

text: We recently upgraded to Image Server 10.0.1 and in the process our WMS Image Services quit displaying. The connection to them works in ArcCatalog, but it never displays. If I preview the service it gives coordinates but no image. I will admit that we had this problem when we originally created them and I figured out the problem...but that was over a year ago and I can't remember that far back. The XML can be viewed here: https://maps.sjcounty.net/ArcGIS/services/2009_SJCounty_1ft_IS/ImageServer/WMSServer?request=GetCapabilities&service=WMS Any ideas?


title: Extent is zooming map to 0,0

text: When I use this code (from the example):                       ESRI.ArcGIS.Client.Geometry.Envelope myEnvelope = new ESRI.ArcGIS.Client.Geometry.Envelope();             myEnvelope.XMin = -130;             myEnvelope.YMin = 5;             myEnvelope.XMax = -70;             myEnvelope.YMax = 65;             mapControl.Extent = myEnvelope;      The map is centered over that ever-popular and fascinating point off the coast of Africa: 0,0. What am I doing wrong?f


title: How can I turn QGIS into a simple hiking planner?

text: I'd like to try and turn QGIS into a simple hiking planner. I have found and worked out how to import OS raster mapping at the correct scale and how to add freehand vector lines to the mapping, which is a start. I'd now like to take things a bit further and make QGIS do something more useful. I realise there are commercial products available that do a similar thing, but they cost money! I have started to play around with extra attributes to the vector lines I am adding in my vector 'routes' layer, such as length in Km, length in miles etc. However, I'm not sure how to do this effectively. The length of the line in Km is already encapsulated in each line in the (Derived) section of the id. I'd like to grab that and put it in its own, easier to view column. I'd also like to automatically have QGIS convert this to miles in another column. Finally, if possible, I'd like QGIS to calculate and somehow display the length and bearing of each leg (ie the distance between each node), which I suspect might be trickier. Is any of this possible?


title: Get Raster Catalog Names from a File Geodatabase Programmatically

text: I'm trying to get only raster catalog entries from a file geodatabase using the `arcgisscripting` geoprocessing API. So far I am able to do this using the `ListDatasets` function using the `ALL` argument. The problem is that it also returns raster datasets as well and I suspect there are other types it would return if they were present in the file geodatabase. I've tried the `Raster` argument but it returns the raster datasets and omits the raster catalogs. I guess I could filter the items returned from `gp.ListDatasets("", "Raster")` from the items returned from `gp.ListDatasets("", "ALL")` leaving me with what I'm looking for, however there might be some other type I don't know of that would be left in the list if it were present. Here is the code I'm starting with:               import arcgisscripting     gp = arcgisscripting.create(9.3)     gp.workspace = r"C:\FileGeodatabase"     datasets = gp.ListDatasets("", "ALL")     for dataset in dataset:         #Only want raster catalogs here...      I'm looking for a way to get **only** raster catalogs. What approach would ensure I'm working only with raster catalogs?


title: How to get maximum value where two shapefiles overlap using ArcGIS for Desktop 10.1?

text: I'm trying to merge 2 datasets (both shapefiles, but I can convert them to raster layers if necessary) where the values are compared, and whichever is higher becomes the new value for the layer. How would I got about this? They aren't quite adjacent to one another, because there's a fairly large overlap in the middle. I'm using ArcMap 10.1


title: Problem import shapefile in PostGIS database

text: I am having a Polygon shape-file. When I try to import this shape-file in PostGIS some of the attributes in Geom Field is found empty. But in attribute table when I select that polygon It is visible in QGIS. What is mistake I have committed?


title: How to drill down through multiple polygons/points in QGIS?

text: I have a layer that has multiple polygons sitting on top of each other, is there a way to easily scroll through those polygons that overlap each other? Thanks


title: Which Windows OS is supported?

text: Please can you tell me which Windows OS Quantum Desktop supports. i.e. Windows 7, XP, Vista etc.


title: How to load a NetCDF-file into QGIS using Python?

text: I have a CF-conform NetCDF-file which can be loaded into Quantum GIS (Version 1.8) by using the "Add Raster Layer"-button and subsequently choosing the file and the raster layers of interest. This works without problems. However, when trying to load the same file by following these steps from the python cookbook, the layer remains invalid. I tried it using the console:               >>> fileName = "/path/to/file/myfile.nc"     >>> fileInfo = QtCore.QFileInfo(fileName)     >>> baseName = fileInfo.baseName()     >>> rlayer = QgsRasterLayer(fileName, baseName)     >>> rlayer.isValid()     False      When trying to add the layer to the map registry, it doesn't show up:               >>> QgsMapLayerRegistry.instance().addMapLayer(rlayer)     <qgis.core.QgsRasterLayer object at 0x3715f80>      I tried the same procedure with a jpeg with world file and it works. The cookbook mentions that GDAL might not support some formats. I don't know how to check this in QGIS, but since it works manually, I don't think that this is the problem. Maybe, I have to specify the raster layer of interest somewhere in the code? Any help is highly appreciated.


title: ArcGIS Server 9.3.1 - Creating a MXD and Exporting a PDF From It

text: I've been tasked with developing a geoprocessing service in ArcGIS Server 9.3.1 that in addition to providing the resultant feature sets, it needs to produce a pdf in which the layers can be toggled on & off. Upgrading to ArcGIS Server 10.0 is not going to happen (even with my begging & pleading). The most simplistic workflow possible is outlined below. Step 1 is easy, I am already doing this with the python geoprocessor - steps 2 and 3 are where it gets difficult. I know creating a mxd and exporting to pdf are **NOT** possible with the python gp, so what are my options? I have been looking at some SOE examples, but they seem far to complex for my needs.   1. Geoprocessing service does some magical processing and creates a geodatabase with result feature classes.   2. Create a mxd document containing result feature classes and apply symbology from layer files   3. Export pdf from mxd - store on server somewhere. Just to clarify, the mxd is only temporary to create the pdf. No services will need to come from it and it can be destroyed after pdf creation. This workflow be initiated from a Flex app.


title: How do I rename rasters in an mxd by using my csv?

text: I have the following data in a csv:                   RasterName  First name         Raster 1    Aaron         Raster 1    Paul         Raster 2    Bryan         Raster 3    Cranston         Raster 3    Vince      I have added these rasters to my mxd as many times as necessary (twice for Raster 1, once for Raster 2 and twice for Raster 3), so that my TOC looks like the `RasterName` column. Now I need to rename the raster in my map document to the appropriate `First name`. I read the csv into csvdict using reader and when I use the following:               csvdict = []     for row in reader:             csvdict.append(row)          for layer in arcpy.mapping.ListLayers(mxd):           for i,j in csvdict:             if layer.name == i:                 layer.name = j      I get this:                   Aaron         Aaron         Bryan         Cranston         Cranston      When I really want:                   Aaron         Paul         Bryan         Cranston         Vince      I'm sure I'm missing something simple, but I haven't been able to figure it out.


title: How to generate a thumbnail of an openlayers map?

text: I would like to allow the user to generate a thumbnail image of what the user is seeing in the OpenLayers map canvas so the user may print, save, share, email, or whatever. **Any ideas of how i would go about doing this?** I'm looking for solutions that generate an image. I don't want to save the innerHTML, feature data or anything other than an image. An image is much easier to share and cannot later be modified. Examples: ![enter image description here](http://i.stack.imgur.com/v2nb8.png) or this will also work: ![enter image description here](http://i.stack.imgur.com/KTHOy.png)


title: Distance between points and polygons in QGIS / PostGIS

text: I need to work out the distance between a table of polygons and a table of points. The table of polygons has different site boundaries, and the point table has points relating to each of the site boundaries. I need to work out the distance between each point and the site it relates to. What is the best way to do this? Is this best done with a PostGIS query or is there a way to do it in QGIS (I notice there is distance matrix tool in QGIS but it doesn't support polygons). Ideally I would like a way to do this in QGIS as its something my users do quite regularly.


title: ArcScene: creating 3D polygon features from a DEM

text: I have a DEM with a lot of topographic variance and a separate layer of polygon features (w/o Z values)throughout. When I use 3D Analyst to convert features(polygon) to 3D, only a portion of the polygons are overlaid. Many of the polygons seem to "be inside the DEM" or underneath, specifically the polygons that appear to be on a hillside. I've tried several extrusion methods, but can not find a clear visualization of the data. Any recommendations are welcomed.


title: How to Set the Vertical Coordinate system when creating 3D Shapefile?

text: I have data of wells in CSV format, which I'm converting to a 3D shapefile, using ogr2ogr. The latitude and longitude are in WGS84, while the depth is from the WGS84 geoid. I have the OSGeo4W installation of ogr2ogr. The process completes successfully, and Points with XYZ are created; However I can't figure out how to set the Vertical Coordinate system when creating this shapefile. My CSV file is like this:               id, name, Lat, Long, depth     1, A, 12, 35, -10     2, B, 13, 40, -15     3, C, 15, 45, -20      My vrt file is like this:               <OGRVRTDataSource>         <OGRVRTLayer name="data">             <SrcDataSource>data.csv</SrcDataSource>             <GeometryType>wkbPoint25D</GeometryType>             <LayerSRS>WGS84</LayerSRS>             <GeometryField encoding="PointFromColumns" x="Long" y="Lat" z="depth"/>         </OGRVRTLayer>     </OGRVRTDataSource>      I'm using the following command to create the shapefile: `ogr2ogr -f "ESRI Shapefile" input.shp data.vrt` How to Set the Vertical Coordinate system when creating 3D Shapefile with ogr2ogr?


title: Enabling SAGA in QGIS

text: A similiar question has been asked recently, however I have some specific queries and process I have tried which are probably better in its own thread for clarity. Hope this is ok. See link to recent question. How to make GRASS and SAGA tools capable to use in QGIS？. The process I used is as follows: 1) Installed QGIS using OSGeo4W. Installed Dufour 2.1 under Program files. Apps directory shows saga_gui.exe (2.0.8). Grass plugin working ok. Install looks ok. 2) The SAGA instructions mention a configuration dialogue without being a bit more accurate. https://www.qgis.org/en/docs/user_manual/processing/3rdParty.html. I assume this is in QGIS Settings \ Options \System variables, and that you should be able to put in the path to the .exe? If so I cannot seem to get this to work. Or maybe it's somewhere else. 3) Also mentioned is that SAGA 2.1 is needed. I've downloaded this and, as it appears to be a "standalone programme", put the whole thing into the Dufour folder rather than reinstalling everything with the aim of configuring the path to the .exe from QGIS. Is this correct, as usually installs put files all over the place? From my quick look at the program SAGA seems to look pretty good. However there appears to be only basic file support (dxf, txt, etc) and no ability to directly load MapInfo files which is what I was interested in. Clarity on running from QGIS appreciated. Thanks.


title: Create base maps from vector data

text: I may be trying to attack this problem incorrectly, so I'm sure I'll be corrected. I currently have quite a few shapefiles containing all of the major highways, streets, railroads, counties, etc. of Texas. I need to combine all of these files so that they can be used as a base map and then pre-render the tiles that will be served out to the client. I haven't found a way to do this and was wondering if there is a way that this can be done or a different way to go about it. The reason that we are using the shapefiles as a base map is because we have to have historical maps. In other words, if a new road is laid and subsequently added to the maps, we will have to keep our old maps and create new tiles for the changes so that if the user has to see why a specific road wasn't taken in the past they can see that at that time the road didn't exist. **UPDATE: A bit more context.** The shapefiles will be loaded into SQL Server 2012 for performing queries, and thus it will be running on a Windows box. I understand that I can use the SQL Server plugin for GeoServer to create a store from the database and then layers from the tables, but the problem with that is there will be several layers used to make up a coherent map, i.e., streets, major highways, counties, etc. In that scenario the client will have to connect to at least 5 layers just to get the map on the page. I don't know enough about GeoServer to know if it will process all of those layers and then send out tiles containing the content for all specified layers or if it is going to send tiles for each layer laying them on top of each other. What I want to do is simply take all of the shapefiles and essentially merge them, applying styles, and then generate the tiles from the resulting map so there would only be one set of tiles being served. Additional Notes: \- Application is being developed in ASP.NET MVC using C#.


title: Software/scripts to calculate Frontal Area Index (λF)?

text: Maybe you know some software/scripts to calculate Frontal Area Index (λF)? This **article** (page 2) describes that index. For now I have not found any software or script that can calculate that (even though google gives large number of hits on "frontal area index arcgis"). In some articles I found, there is a mention that ArcGIS can do it (probably yes, but at this point I cannot create script to do it).


title: Calculating percentiles of a raster (TRMM DATA)

text: I have a set of raster (floating point) i.e A daily rainfall intensity data for a year, now i need to select and group these raster in to various percentile categories, Above 90th percentile, 50-90, 10-50, below 10. how to achieve this using ArcGIS 10.2 ?


title: autocad-map of real estates with own coordinate system and information to be converted to use it in QGIS

text: I have a question regarding CAD and QGIS, but unfortunately I don't have a clue of CAD. Here it goes: Real estates data has been mapped. There are several reference points surrounding the area. So there is an own coordinate system being used. The coordinates of the corners of each polygon are known. Each polygon holds different information regarding size, owner etc. The task is to have this finished mapping converted to be used in QGIS. The person wants me to have it correctly transferred into UTM-PCS. The outcome is an editable map of the real estates and each polygon is given it's equivalent attribut table. Any hints???


title: How to get mouseposition (openlayers) in decimal degrees?

text: I want to display the **mouseposition** in **Degree Decimal** rather than in **meters**. ![enter image description here](http://i.stack.imgur.com/tXBN7.png) See my code:                   <script type="text/javascript">         var map, ia_wms;         function init(){             map = new OpenLayers.Map('map', {zoomDuration: 1,projection: 'EPSG:3857'});         var esr_all = new OpenLayers.Layer.WMS("ESR_Boundary","http://localhost/cgi-bin/mapserv?map=/var/www/mstest/test_sr.map",{'layers':"ESR_outline",'transparent':true,'format':'image/png'}                ,{                 isBaseLayer: false,opacity: 100                 }                 );         var mslayer = new OpenLayers.Layer.WMS( "MapServer Layer",                                               "http://localhost/cgi-bin/mapserv?map=/var/www/osm-demo/mapserver-utils-svn/osm-google_test.map&",                                               {'layers': "land0 borders0 places0 land1 borders1 places1 land2 borders2 places2 land3 borders3 places3 land4 landuse4 transport_areas4 waterarea4 borders4 places4 land5 landuse5 transport_areas5 waterarea5 roads5 borders5 places5 land6 landuse6 transport_areas6 waterarea6 waterways6 roads6 borders6 places6 land7 landuse7 transport_areas7 waterarea7 waterways7 roads7 borders7 places7 land8 landuse8 transport_areas8 waterarea8 waterways8 railways8 roads8 borders8 places8 land9 landuse9 transport_areas9 waterarea9 waterways9 railways9 roads9 borders9 places9 land10 landuse10 transport_areas10 waterarea10 waterways10 railways10 roads10 aeroways10 borders10 places10 land11 landuse11 transport_areas11 waterarea11 waterways11 railways11 roads11 aeroways11 borders11 places11 land12 landuse12 transport_areas12 waterarea12 waterways12 railways12 roads12 aeroways12 borders12 places12 land13 landuse13 transport_areas13 waterarea13 waterways13 railways13 roads13 aeroways13 borders13 places13 land14 landuse14 transport_areas14 waterarea14 waterways14 railways14 roads14 aeroways14 borders14 places14 land15 landuse15 transport_areas15 waterarea15 waterways15 railways15 roads15 aeroways15 borders15 places15 land16 landuse16 transport_areas16 waterarea16 waterways16 railways16 roads16 aeroways16 borders16 places16 land17 landuse17 transport_areas17 waterarea17 waterways17 railways17 roads17 aeroways17 borders17 places17 land18 landuse18 transport_areas18 waterarea18 waterways18 railways18 roads18 aeroways18 borders18 places18"}                                                );             map.addLayers([esr_all, mslayer]);         map.addControl(new OpenLayers.Control.MousePosition(         {                     prefix: 'coordinates: ',                     separator: ' | ',                     //numDigits: 2,                     //emptyString: 'Mouse is not over map.'                 })                     );               map.addControl(new OpenLayers.Control.LayerSwitcher());         map.zoomToExtent(new OpenLayers.Bounds(69.00,21.00,70.00,22.00));         map.setCenter(new OpenLayers.LonLat(69.57,21.55).transform(         new OpenLayers.Projection('EPSG:4326', 'EPSG:3857'),         map.getProjectionObject()         ), 14);         }     </script>      In mapfile I have given: > >     WEB >     IMAGEPATH '/opt/fgs/tmp/ms_tmp/' >     IMAGEURL  '/ms_tmp/' >     TEMPLATE  '/var/www/mstest/templatesample.html' >         METADATA >             "wms_enable_request"    "*" >             #"ows_enable_request" "*" >             "wms_srs" "EPSG:4326 EPSG:900913 EPSG:3857 EPSG:2154 > EPSG:310642901 EPSG:4171 EPSG:310024802 EPSG:310915814 EPSG:310486805 > EPSG:310702807 EPSG:310700806 EPSG:310547809 EPSG:310706808 EPSG:310642810 > EPSG:310642801 EPSG:310642812 EPSG:310032811 EPSG:310642813 EPSG:2986" >             "labelcache_map_edge_buffer" "-10" >             "wms_title" "osm" >         END >      END >   > > DEBUG 1 CONFIG "MS_ERRORFILE" "stderr" CONFIG "PROJ_LIB" "/var/www/osm- > demo/mapserver-utils-svn" PROJECTION "init=epsg:900913" END I tried to change the init=epsg to 4326 , But few layers(osm data from postgres db) ended up not drawing. How to tackle this?


title: Terminology problem - Raster vs Grid vs Lattice

text: I have a question about GIS terminology - what is the difference between grid and raster and grid and lattice - in theory?? I thought, that raster is regular network and value is in the cell and the grid have its value in node. In this case i don´t understand what meaning have a lattice? Is the same as a grid? Or grid is raster with value in one of the corner?? What is more confusing are some pictures, which describes difference between grid and lattice - http://www.geog.ubc.ca/courses/geob370/notes/interpolation/grid_vs_lattice.png or http://www.innovativegis.com/basis/mapanalysis/Topic18/Topic18_files/image005.jpg In case of Matlab example: http://www.mathworks.com/help/releases/R2013b/matlab/ref/meshgrid_ex.gif grid looks like lattice! Please can anybody helps? Thaks a lot.


title: (ms4w) oracle spatilal problem of drawing compound polygon for Oracle 10g

text: When i draw a compound polygon SDO_GEOMETRY( 2003, -- two-dimensional polygon NULL, NULL, SDO_ELEM_INFO_ARRAY(1,1005,2, 1,2,1, 5,2,2), -- compound polygon SDO_ORDINATE_ARRAY(6,10, 10,1, 14,10, 10,14, 6,10) ) by using msw4 accompanied with mapserver 5.4 or the previous version, they are successful to create an image. But when i change to msw4 3.0.6 with mapserver 6.0.3 (libmap.dll of Oracle 11g), i always get the problem : msDrawShape(): General error message. Only polygon shapes can be drawn using a POLYGON layer definition. I think it comes from the file libmap.dll, maybe libmap.dll for Orl 11g doesn't work welll with Orc 10g Can you suggest me what to do ? Iam using Oracle client 10.2.0.3, and the underlying database is Oracle db 10g (10.2.0.4)


title: QGIS and usage of ExportToKML()

text: I am trying to dump to a csv file the geometry (and attributes) of a shapefile, including a KML dump (using a python script as described in shp-to- csv) but I am unable to get the correct coordinates (while exporting the shapefile as a whole, using save as kml provides a correct dump):               # Write attributes and kml out to csv     for feat in lyr:         attributes = feat.items()         geom       = feat.GetGeometryRef()         attributes['CAMM']        = camm         attributes['kmlgeometry'] = geom.ExportToKML()         attributes['centroid']    = geom.Centroid()         attributes['kmlcentroid']    = geom.Centroid().ExportToKML()         csvwriter.writerow(attributes)         print(attributes)         break      which results, e.g., in coordinates such as               <Point><coordinates>175.56640841253102,5096344.3355032187</coordinates></Point>      instead of something close to 11,46 ... Do I need to perform explicitly some transformation before invoking the kml export function? Thanks,


title: How to implement subtype functionality in Oracle Spatial?

text: I am planning to shift my entire database from arcgis to oracle. I will not be using ArcSDE services also. I need subtitutes for ESRI's functionalities like domain, subtype, relationships in oracle. I have some how worked out two them i.e. domains and relationships, but no alternatives for subtypes have striked me.


title: AML routines and python

text: I´m trying to use some AML files in ArcMap 10, but I have noticed that this option is not available for this version. Does anyone know or tried to "translate" the aml sequences to another language for example python? or another application like QGIS or SAGA. Regards.


title: free way to author, share ESRI gis data

text: I hope someone can help. THE ISSUE: I work for a company where 90% of the personnel work remotely. The GIS department acts as a support tool authoring content for a variety of projects. There is a need to create a workflow where employees can access basic base maps to view authored content on base maps without having to put in a request to the GIS department. The hope is that someone can RDP to a computer or login via a website to access data and to capture basic elements of a geographic location (lat/long, geographic extent, etc) while viewing and printing authored content. CURRENT WORKFLOW: Employees put in a request via email to the GIS department. This workflows creates a backlog. QUESTION: Is there a free solution that can be utilized out of the box with little or no required knowledge of programming code to set up the tool? thanks Jay


title: How to create a polyline feature from a point layer

text: I would like to create a route line based on a point layer. The existing tool does not meet the customer requirement. But using the code snippet below, the start point and end point are always the same. It means they point to the same instance? Please help. Thanks.               rows = arcpy.SearchCursor(pntLyr, whereStr, "", "", "STOPA A")      start=True     for row in rows:               if (start==True):             start_feat=row             start=False         else:             end_feat=row             # create the line based on the start_feat.shape and end_feat.shape             # Before I reassign the start_feat, start_feat and end_feat already point to the same feature now at this point             start_feat=row


How do I convert from points to raster making the point the lower left corner of the cell? I'm using ArcGIS 10.1. The default for the "Convert Point to Raster" tool is to make the point the centroid of the resulting raster cell.


title: QGIS 2.2 - expression string builder to ignore NULL values

text: Simply put, I would like to create a rule-based filter whereby you add 2 columns together and you get a certain colour as a result. However, some fields in **Column_1** may be NULL and those exact same fields may have values in **Column_2** and vice-versa. Is there an expression to add the 2 columns whilst ignoring the NULL values or atleast treating them as zero? The following is the logic I am trying to achieve: 1 + 2 = 3 1 + NULL = 1 NULL + 2 = 2 Additional info: The data from the columns are JOINED from other layers, I do not want to create any new columns as I am wanting the colours of layers to dynamically change when the user adds or removes the JOINED layers. Thank you.


title: QGIS autotrace plugin - when will it be ready?

text: Is there any news on when the MapInfo-style auto-trace plugin that Lutra Consulting has been developing, will be put up to the QGIS Plugins Repository? It was supposed to be this month, according to comments on the YouTube video for it. See http://www.youtube.com/watch?v=sKRu_MJcbmo


title: How to design a road with a certain slope

text: Is there a way that you want to connect two different points of altitude with a constant slope of 7% or another, to draw a line in the 3D model?


title: ArCGIS Server 10/10.1 And Sql Server 2008

text: Is there any compulsion that SQL server 2008 and ArcGIS server 10/10.1 should be installed on same server machine and is it possible to install SQL server 2008 and ArcGis server on separate server machines?If is it possible then what are the complications and performance issue involved in above scenario.?


title: ArcGIS cached map service sometimes pausing

text: I have published an ArcGIS service based on a MXD I have generated the cached tiles on a shared NAS. When I access the Map Service with the build in js admin viewer (http://SERVER/arcgis/rest/service/....) The map display quickly but, every now and then no tiles are coming for maybe 5 seconds, and then they all come quickly. Why could that be? Any tip / ideas welcome? How to diagnose the problem?


title: problem after merge shapefiles, can't show it on OpenLayers with GeoServer

text: I have 5 .shp files, after merge it with geomerge, I can't use it on OpenLayers with GeoServer. it's only show white layer and when i click it, it show error message: Internal error occurred ShapeType changed illegally from Arc to Point . what should I do?


title: Reference for IDockableWindow?

text: Kind of a newbie with ArcObjects .NET. Trying to open a dockable window with the following code:               public class CAD_Test_01_button : ESRI.ArcGIS.Desktop.AddIns.Button     {         public CAD_Test_01_button()         {         }              protected override void OnClick()         {                  UID dockWinID = new UIDClass();             dockWinID.Value = ThisAddIn.IDs.DockableWindow1;                  // Use GetDockableWindow directly as we want the client IDockableWindow not the internal class             IDockableWindow dockWindow = ArcMap.DockableWindowManager.GetDockableWindow(dockWinID);             dockWindow.Show(true);              }              protected override void OnUpdate()         {         }     }      But I get the errors: Error: The type or namespace name 'IDockableWindow' could not be found (are you missing a using directive or an assembly reference?) Error: The name 'ArcMap' does not exist in the current context What references could I be missing? Thx, Rich


title: Error When Trying to Run Slope Analysis Script

text: I have written a script which is supposed to loop through all the rasters in a folder and perform a slope analysis on each raster and save the results in another folder. I keep getting the following error: > < type 'exceptions.RuntimeError'>: ERROR 000875: Output raster: > C:\Slope_Outputs\burd_25m_dem_slope's workspace is an invalid output > workspace. Here's my code:               class LicenseError(Exception):         pass          # Set desktop license used to ArcView     #     import arcview     import arcpy     from arcpy import env          try:         if arcpy.CheckExtension("Spatial") == "Available":             arcpy.CheckOutExtension("Spatial")         else:             # Raise a custom exception             raise LicenseError          except LicenseError:         arcpy.AddMessage("Spatial Analyst license is unavailable")     except:         print arcpy.GetMessages(2)          from arcpy.sa import *          originLocation = arcpy.GetParameterAsText(0)     slopeMeasurement = arcpy.GetParameterAsText(1)     destinationLocation = arcpy.GetParameterAsText(2)          arcpy.env.workspace = originLocation          rasterList = arcpy.ListRasters("*")          for raster in rasterList:         finalDestination = destinationLocation+"\\"+raster+"_slope"         arcpy.AddMessage(finalDestination)         outSlope=arcpy.sa.Slope(raster, slopeMeasurement)         arcpy.AddMessage(outSlope)         outSlope.save(finalDestination)      I'm not really sure why I'm getting this error. Any ideas? Thanks.


title: Why unable to create interactive FeatureSet features for ModelBuilder tool?

text: I am attempting to set up an "interactive clip" tool that will allow the user to draw clip features on the fly and use those as the clip features in the Clip tool in the analysis toolbox. I have added the Clip tool to ModelBuilder and exposed the clip features as a model parameter. I also changed the type to a FeatureSet, which I understand is supposed to allow the user to interactively enter the features when the tool is opened within ArcMap. http://help.arcgis.com/en/arcgisdesktop/10.0/help/index.html#//002w0000004t000000 However, when I open the tool in ArcMap, the Clip Features input does not allow the interactive drawing of features, but rather only allows the user to select an already existing feature class. ArcGIS help and the above link indicate that under this setup the FeatureSet should by default point to an in_memory location where the user can temporarily draw features. Here's how I have it set up in ModelBuilder: ![tool interface](http://i.imgur.com/XfG3FCC.jpg) Here's the resulting tool interface: ![model builder](http://i.imgur.com/Kd2qxpy.jpg) What am I missing here? Any ideas why this is not working as expected?


title: Problem to transform events from one route to another

text: I have a road feature class layer containing parts like 1_10, 1_20. 2_10, 2_20... etc. 1_, 2_ is later used to create the route and _10, _20 is indicating the direction (the road feature classes were digitized in the same direction - they told me at least!). In this feature classe I calculate the "from" and "to" values for each segment. Ater that I create a route with arcpy.CreateRoutes_lr that seems to work great. Now for this route, linear events are collected over time. After some time there are changes in the road feature class. So my arcpy script now recalculates the new lengths of the road segments and their "from" and "to" values, then creates a new route. This all seems to work fine BUT I am not able to transform the old events from the old route to the new route with arcpy.TransformRouteEvents_lr. I get a lot of event artifacts (like 0.5m long) and even events from lets say route 7 are "mapped" to route 10, or Event X is on route 10 and 11... However I think the most errors are like ths (the events are basing on the original segments):   * Event1 on route1 from 0 to 20   * Event2 on route1 from 20 to 30   * Event3 on route1 from 0 to 25 ??!!! Another problem is that my routes consists of loops which I think is an additional challenge. I also tried a workaround as creating a new feature class with the old routes and the old events and then with arcpy.LocateFeaturesAlongRoutes_lr to map it on the new routes but that does not seem to work either (actually I have quite similar results as with the transform function). I dont want to map old events on a changed/new routes (since that makes no sence for pavement). And this seems to work at least. The problems are events on routes that have not changed, but maybe the from/to values of the events must be recalculated. I hope this wall of text makes sense somehow and unfortunately I can not post pictures yet. Any help is really appreciated! Edit: after doing some research I found out that all features that are not a single line create errors so for example Y-or T-junctions in a route, or looping routes.