texts:
I have been trying for a long time now to make this simple basemap work. However, there is something with the format of the data basemap expects, which I cannot figure out. Can anybody help me understand what am I doing wrong?

from mpl_toolkits.basemap import Basemap
import matplotlib.pyplot as plt
from osgeo import gdal
from numpy import linspace
from numpy import meshgrid
import numpy as np

lats  = [-17.96, -2.88, -3.12, -14.01, -9.28, -15.59, -15.69]
longs = [-38.69, -40.13, -40.04, -52.21, -37.89, -47.62, -41.43]
prod  = [ 24.6, 7.5, 8.5, 68.8, 6.8, 6.6, 31.0]

lats = np.asarray(lats)
longs = np.asarray(longs)
prod = np.asarray(prod)

map = Basemap(projection='tmerc',
          lat_0=0, lon_0=3,
          llcrnrlon=1.819757266426611,
          llcrnrlat=41.583851612359275,
          urcrnrlon=1.841589961763497,
          urcrnrlat=41.598674173123)

x = linspace(0, map.urcrnrx, prod.shape[1])
y = linspace(0, map.urcrnry, prod.shape[0])

xx, yy = meshgrid(x, y)

map.pcolormesh(xx, yy, data)

plt.show()


texts: Project only required shapefiles in folder

texts:
I have 8 shapefiles in a folder with shapefile 1 and 8 in the same projection and the others (shapefile 2-7) if different projections.

Is there a way to use python to use shapefile 1 as the projection template to project shapefiles 2-7 in a new folder but also to only copy shapefiles 1 and 8 in the output folder without reprojection them but also not using hard values to keep the script flexible?

I have been trying to use batch project, but am having trouble excluding shapefiles 1 and 8, or only including shapefiles 2-7. Should I be using the Describe.SpatialReference to filter out?


texts: Generating URL using USGS Site-Test-Tool to add as Layer in ArcGIS Online?

texts:
I am trying to use http://waterservices.usgs.gov/rest/Site-Test-Tool.html to generate a Rest URL to add as a Layer in ArcGIS Online.

I generated url (http://waterservices.usgs.gov/nwis/site/?format=ge&stateCd=va&outputDataTypeCd=iv) using USGS Site Web Service URL Generation Tool but unable to add layer in ArcGIS Online using this rest url.


texts: Finding a route with desired cost not minimum cost

texts:
The goal of routing is usually finding the route with minimum cost. That cost can be distance, time, calories or some other value but still our desire is to find the route that costs least.

But what if I want to find a route with desired cost? Not the fastest route but the route that will take around an hour (let's assume I have an hour to walk from point A to B and I want to see as much on my way there as I can). Or that I want a route that will be around 2 kilometres long.

Using a standard routing algorithm and rejecting the proposed routes until one with desired cost is found will probably lead to something like "Head straight from A to B and then circle around B a few times".

Is there any previous work on this? Does anyone have any ideas?


texts:
I have multiple large polygon datasets stored in PostgreSQL 12 and will publish them using GeoServer. I want to try to use geohash to possibly speed up the performance of serving the polygon layers. I made a geohash index on one of the larger polygon datasets in PostgreSQL and it seems to have slightly improved performance but that just might be a hopeful placebo effect.

Anyone know of any way to use geohash to speed up serving large polygon data?


texts:
I need to perform hydrological analyses and map design purposes (hillshade creation) on a region far north. It is too northern to be covered by the SRTM DEM, therefore I am using the ASTER DEM. However, in some mountainous parts, there are missing pixels, typically occuring in the space of peaks or on the other hand in monotone environments.

The GTOPO30 DEM however, is without missing pixels, even though at a comparatively coarse resolution.

What is the best workaround to resolve the problem of missing pixels? Interpolating between the edges of the ASTER, or perhaps filling these parts with a resampled GTOPO30?

(A partly similiar sounding question was asked here, but without a resolution to the announced problem of the above question.)


texts: Benefits of relationship classes and how they can be accessed in c#

texts:
At the moment I've a personal geodatabase with two imported feature classes from shape files and an own created table. This table has columns for the IDs of the elements in the two feature classes. With ArcCatalog I've created two relationship classes between each sfeature class and the table. Both are 1:1 relationsships. Within ArcMap the user can open the table and get access to more detailed informations by using this relationsships (menu -> Related Tables).

Background: One feature class holds features consisting of Polylines, the other features consisting of Points. With the Table I can create a logical relationship between one Point and one Polyline and can add more informations in this table which are only valid for this special relation.

Now I've to work with this table in C# and I've to access the other feature classes too. It's no problem to read the IDs out of the table and start a QueryFilter for the feature class where I need the informations from.

Sure, this is one possible way but there must exist a more efficient way where I can use the relationship classes for faster operations.

Does anybody know how I can query the relationship classes to get all informations at the same time without querying two different tables "manually"?

Thanks a lot!

Edit:

To be more specific. If the table has the two fields ID_EdgeTo and ID_EdgeFrom and one field has a relationship to the feature class EdgesTo and the other to EdgesFrom.

When I get through the rows and see a combination of ID_EdgeTo and ID_EdgeFrom I want to query the relationship classes to get more information for the ID's in this specific row.

Its more about whether the relationship classes are retrievable and how this has to be done. The table in the database has the form like a "manually created" relationship class with relationship attributes


texts: Reverse polygons (left-hand to right-hand rule)

texts:
I have a geojson file that is not displaying correctly, and I believe it has to do with the file containing multipolygons not following the right-hand rule. What tool can to help me fix this problem?

I've tried importing and exporting with QGIS, and converting back and forth between formats with org2org, but that didn't take me anywhere.


texts: Use multiple connections/databases in a query layer?

texts:
I have access to two databases, each with a separate connection. The first connection is an .sde connection to a SDE 10.3.1 Oracle database, and the second connection is a .odc OLE DB connection to a third party software Oracle database. Is it possible to use both databases/connections in a query layer? I have some fairly complex non-spatial SQL queries I'd like to do using both databases; the SQL includes subqueries, group by, etc.

My guess is that the answer is no; it looks like you can only select one connection when creating a query layer. I've worked through a few alternatives, but they all seem to be lacking in one way or another:

Create a query layer on a table from a single connection/database, then create an in-map join to a table from the other database/connection using the layer properties in ArcMap. This is better than nothing, but is severely limited. I can't do any of the SQL querying that I need to do such as subqueries, group by, etc.
Use the Make Query Table geoprocessing tool, which sems to be able to reference multiple connections/databases. This has similar limitations to #1. You can enter an expression, but I believe the expression is the equivalent to the WHERE clause in SQL, which is limited. I can't do subqueries, group by, etc. here either.
I don't believe it's possible to create a view in one database that would reference another database.


texts: Missing OSM Points after importing osm.pbf into PostGIS DB

texts:
Can you tell me how to import OSM data into PostGIS database without lose any data?

I am using map extracts from https://mapzen.com . At first I download Raw OpenStreetMap datasets in PBF, but after importing it by osm2pgsql, I do not have all points. For example when some street have numbers from 1 to 20, its 20 points, but after import I have only 4. But when I download Datasets split by geometry type: lines, points, or polygons (SHP) and import pure shp file into PostGIS, I have all points with that street (1 - 20).

I need this data for 2 things.

geocoding - i need compare text value of point (city, street, street number) with user entered text values in routing
generate base map and style it in geoserver

Is possible set up this tool to import all points?

Or is some other way for windows?


texts: Google Earth Engine: How to import data from Google Cloud Storage

texts:
I have finished the data process in Google Colab and upload these results to Google Cloud Storage (GCS), and I want to import these files from GCS to Google Earth Engine (GEE) via the python API provided by GEE.
The files in GCS looks like this:

While the code in Colab makes sure that they are exist:

However, when I write the import code:

outputAssetID = 'users/cjj/modis_process'
!earthengine upload image --asset_id={outputAssetID} {outputImageFile} {outputJsonFile}

The message I got is:

Running command using Cloud API.  Set --no-use_cloud_api to go back to using the API

Manifests for TfRecord ingestion must have exactly one tileset with exactly one source.

Am I using this API in the wrong way?


texts: Trim symbology label decimal places when using Arcade in ArcGIS Pro?

texts:
I have an issue trying to symbolise a layer with graduated symbols using a custom expression:

I am limited to using arcade/there are no other languages in the
drop down which I don't believe should be the case.
An arcade custom expression, even $feature['FieldName'], is producing correct values/band but the corresponding labels have 6 decimal places appended where using the field alone (no custom expression) is not. No casting and/or round, trim or replace etc will fix this.

I am using ArcGIS Pro 2.0.1.

I can manually edit the labels but any changes I make to symbology will bring back the decimal places.


texts:
I have an OpenLayers map.

On to this I have loaded a GeoJSON file.

When I pan the map the GeoJSON data moves relative to the basemap.

Surely this isn't right! Can anyone shed some light no why this is happening?


texts: LecoS stalling - table not being generated

texts:
I'm trying to collect data along highway segments in northern Alberta. I took vector data from ABMI (http://www.abmi.ca/home.html) and rasterized this to work within LecoS. I had success in early attempts to extract spatial metrics using the LecoS vector overlay function. The road segments are wildlife mortality (animal vehicle collision) hot spots and cold spots estimated using KDE+ software. I selected the center of each segment to create a 500 m circular buffer to sample the surrounding landscape. I have other raster layers that I'm hoping to use in my summary analysis and will be expanding the sampling to different buffer sizes. The following figure shows the vector overlay sampling scheme - a 500 m radial buffer from the median of each segment.

Here's a video example of LecoS (https://www.youtube.com/watch?v=mDpowxRzgOY) if anyone doesn't know about this plugin.

Everything was working in early trials with LecoS. I saved my workspace and the next day found that I could no longer extract LecoS metrics. I went through the entire process of rasterizing, projecting, and clipping the data to see if going from scratch would work. I tried uninstalling and re-installing LecoS, but this has not solved my problem. I clipped the raster in an attempt to make the analysis space smaller. I've kept the analysis simple by selecting only one metric at a time.

I'm using QGIS V. 2.18.21. LecoS prompts a response in the bottom left of my workspace:

  "Land cover calculated for feature 58 out of 58 (0 impossible)"

However, nothing happens after this - it stalls. I left it running overnight - still no table is generated. It shouldn't be a matter of waiting, because it took only a few minutes to get a table in my earlier attempts. I've tried saving to a file, adding table to workspace, and save in the attribute table - none of these options are working.

I've uploaded my raster and the buffer into Dropbox to see if I someone can try to execute and recreate the error. The vector data came from ABMI - so it is projected in Alberta 10-TM (Forest) EPSG:3400, NAD83.

https://www.dropbox.com/sh/uuzv8iw6ek09afd/AABwAGkBDBXLH3sfawm5ItxLa?dl=0

UPDATE: I clicked on the option to append the data to the Attribute Table. This worked on one trial occasion, but the excel file is blank. I had to re-open the vector file to see the appended data in the attribute table.

The process is slow because I have to create a unique vector layer to correspond to each raster layer. Further, I can only do a few metrics at a time. This workaround process creates a whole bunch of sub-files that I have piece together. I'm including a graphic overview to help contextualize the issue.

There is the option to batch process by selecting through the toolbox, but the problem still exists that a separate vector needs to be created for each type of analysis (e.g., Roads-500m.shp, Agriculture-500m.shp or Roads-1000m.shp, Agriculture-1000m.shp). It wouldn't be too bad if I could save the results as a csv - because I can amalgamate in R, but the LecoS software is not creating the csv files. I can only get the data appended to the attribute table (sometimes - it stalls often). I am working with a laptop Intel Core i7 @2.60 GHz 16GB RAM - GeForce GTX 960M. Might be a processing issue? Any suggestions?

I've posted this question on the gitHub for the software as well - but thought it would also be useful to post the question here in case someone else has had this experience and could provide an answer or suggestions to make this work more effectively.


texts:
I'm pretty new to QGIS (2.14.3) and am half way between singing its praises and hitting my head against the desk. I've been having a particularly vexing problem recently with a map I've created.

In the image (Fig 1) below you will see that I'm mapping out rivers and riparian areas. I'm projecting into WGS84/Pseudo Mercator.

(1) What I can't understand is that it appears like the OpenLayers image in the background has shifted when the jpg file is being made. In fact if you look closely, neither the full map nor the overview map is in the correct (or even the SAME) position. However, in the print composer (Fig 2) the image and layers line up perfectly.

(2) The buffer I've created (yellow and red line) should look like it does in the overview map however if you look in the full map where I've zoomed into a particular section, you'll see that it suddenly cuts across from north-west to south-east on the right-hand side of the map. Nothing I do has corrected this, including remaking the layers. This problem seems to persist in print composer, QGIS and the jpeg files. [Edit] There is also mean to be a bridge displayed in red but that doesn't seem to show in the print composer or JPEG either.

Fig 1: Jpeg file

Fig 2: Print composer view

Am I missing something? Or is this just a freak event? I'm very close to just taking a screenshot of the Print Composer and using that at the moment.


texts: Converting dem and geotiff to collada model with texture?

texts:
I have a digital elevation model (dem) file with corresponding orthophoto in geotiff format.

What I need to do is to create a 3D model of the terrain with textures (preferably in collada file)?

I would like an opensource or freeware solution, but I will consider commercial software if need be.


texts:
I am using ArcMap 10.6.1.

I am attempting to label parcels with information from two different fields, but I would like to format those labels separately for the sake of visual hierarchy.

I'm familiar with this expression: [Field1] & vbNewLine & [Field2] which shows the labels I want, in the order I want, but I do not know how to adjust the color/size/font of those labels individually.

The other option is to use the Maplex Lable engine to create classes for each field. This is good, as I can format the text the way I would like but then placement suffers, as I can't seem to recreate the centered, two-line look that the code above produces. For example, setting internal placement rules for the two classes works for small parcels, as they're still grouped together, but falls apart for larger parcels, where the labels get separated.

There are two paths to solve this and I would welcome suggestions for either one: code suggestions for formatting within the expression, so everything can be done in one go, or tips for being more fine-grained with label placement in the label engine (for example, if I recall correctly QGIS has the ability to place labels by a specified offset on a specified axis...this would solve the problem but I'm not sure how to do something similar in ArcMap).


texts: Levelling an existing site in DSM for viewshed analysis in GrassGIS

texts:
I'm doing a multi-point viewshed analysis for an urban site in GrassGIS based on a LIDAR DSM raster maps, to test for a series of proposed buildings. The site has existing building on them.
Is there way of "levelling" the site to the surrounding ground level itself, so that the calculated height would be that of "proposed building above ground", not proposed building above existing buildings would be used?
In my ignorance I have tried to export the DSM to PNG, then then painted the site in in Photoshop, with the colour of the surrounding terrain, hoping that I could re-import it and use it as a basis for viewshed. On import, the PNG is broken into r. g. b . channels, not giving me one map.
I have access to Autocad,so I was hoping there would be a way of drawing up the site boundary as a vector, and using that as a boundary to redefine the level inside the site.
Or any other approaches inside of GrassGIS that achieve the same result.


texts: Pitch Black Bstack_raster.tif image after Stacking Landsat 8 bands using SCP plugin

texts:
I encountered a disruptive problem while following a step-by-step guide of the Land cover supervised classification from a document accessed from https://arset.gsfc.nasa.gov/sites/default/files/land/webinars/Advanced_Land_Classification/supervised_exercise_v2.pdf

Chronology:

I downloaded a Landsat 8 image from the Glovis website – [Path 43, row 33] for date 22nd September, 2015.

I unzipped the folder it was saved in to extract the contents using 7-zip folder.

1st Purpose was to convert the image to Top of Atmosphere (TOA) reflectance

In the Pre-processing icon of the Semi-automatic Classification(SCP) panel, in the part ‘Directory containing Landsat bands’ I put the folder containing the 11 bands of my downloaded Landsat 8 image.

In the ‘Select MTL file’ I manually input my .MTL text file  that is the same directory folder of my Landsat 8 bands (it didn’t appear automatically as the documents says in pg.2, I don’t know why).

My metadata appeared anyway, for all 11 bands.
I clicked the run icon at the bottom right to run the process.

It really took some time to run in my Windows 7 Starter 32 bit Dell Laptop; using Qgis 2.18.11 ‘Las Palmas’. Anyway, this is where the problems started.

I then proceeded to click on the ‘Band set’ icon located at the top of the SCP toolbar.(pg.3 of the exercise).

Bands 2-7 appeared, in tandem with the exercise.
I clicked on ‘Create raster of band set (stack bands) then I clicked on the ‘Run’ icon at bottom right.

This process didn’t take long as the other one of converting to TOA reflectance.

However, the result, which is supposed to be a Bstack_raster.tif image was the astonishing bit. I got a Bstack_raster.tif image alright, but the issue was that it was pitch black.

The Problem:

My Bstack_raster.tif image was pitch black. So I right-clicked it and went to its Properties section – to ‘Style’ in order to display it in a usable color combination.

I put band 5 in the red band, band 4 in the green band and band 3 in the blue band.
My contrast enhancement was set to ‘Stretch to MinMax’.

I then clicked ‘Load’ (in the section under Load min/max values) and then clicked ‘Ok’ but my Bstack_raster.tif image was still pitch black.

From my personal overview, the reason is because all my ‘min’ values in the red, green and blue band sections read ‘nan’. Yet, I have followed each step diligently in the exercise.

Troubleshooting efforts:

I have tried to backtrack all my steps above to see any errors but none identified so far.

I have repeated the same process thrice with the same results.

I have followed all the steps of pg.1-4 of the exercise without omission.

Problem Statement:

Please advice to convert the pitch black
Bstack_raster.tif image to a usable color combination as shown in pg.4 of the exercise. A converted image of TOA is detrimental for my further uses.


texts: Matching/snapping two shapefiles using node and snapping tools

texts:
I am trying to edit my shapefiles through QGIS 2.6.1 Brighton. I used the node tools to snap the necessary nodes to eliminate spaces between two shapefiles. However, while doing it I noticed that an extra line came out. Please see image below - the line that appeared when I'm editing is pointed out by the right and left arrows. I am also using the snapping tools while doing this.

How do I remove these two lines?


texts: How to create a shapefile consisting of 2 lines using Python OGR?

texts:
How would you create a shapefile consisting of 2 lines using Python OGR?

First line is no problem, but how do you reuse MyLine and Feature objects to create the second line?

Ultimately this will be in a loop with thousands of lines created...

from osgeo import ogr
driver = ogr.GetDriverByName('ESRI Shapefile')
datasource = driver.CreateDataSource('c:/temp/testlines.shp')
layer = datasource.CreateLayer('layerName',geom_type=ogr.wkbLineString)

#create first line 0,0 to 10,0:
myLine = ogr.Geometry(type=ogr.wkbLineString)
myLine.AddPoint_2D(0,0)
myLine.AddPoint_2D(10,0)
feature.SetGeometryDirectly(myLine)
layer.CreateFeature(feature)

#create next line 0,0 to 0,10?:
#how to you reuse the myLine, feature objects?
#myLine.DeletePoints()#something like this?
myLine.AddPoint_2D(0,0)
myLine.AddPoint_2D(0,10)
#feature.NesGeometry() #something lik this?
feature.SetGeometry(myLine)
layer.CreateFeature(feature)


texts: Coordinates transformation from MGI 31276 to WGS84

texts:
I have a shapefile in MGI 31276/Balkan zone 6, and it doesn't line up with my base map.

I try to make transformation from MGI 31276 to WGS 1984 Web Mercator Auxiliary Sphere, but it's impossible, and and I don't know why.


texts: ArcGIS: Goode Homolosine Projection and Inconsistent Extents

texts:
I recently posted about a problem in ArcGIS using the Goode Homolosine (Land) projection with the Project Raster tool HERE.

The problem was related to a bug with that tool, and a workaround solution was proposed.

I have since tried to re-project different data to the Goode Homolosine projection, but the workaround solution does not work for these data. I find the proposed solution to work for some data but not others.

Here is an example of the proposed solution that produces an Inconsistent Extent warning:

I downloaded this Crop Suitability Index raster data.
I open these data in ArcMap and set the data frame's coordinate system to projected coordinate systems, world, Goode's Homolosine (Land).
In the table of contents I right click the data and select Data, Export Data.
When exporting I set the Spatial Reference to the Data Frame, and save as a .tif file.
When I open the saved .tif file in a blank map I receive an Inconsistent Extent warning message.

I do not understand why this workaround solution isn't a general solution, i.e., why this solution worked with the Malaria Ecology Index raster data I originally posted about, but not these Crop Suitability Index raster data?


texts:
Our data is of pipelines and nodes wherein for each feature there is a line and a point node as end points.

So while copying parallel lines, we need to customize the OOTB tool so as to copy specific number of lines and create a node point at end of each line.

In general tool is working fine except for the scenario where it is a 'V' shaped feature i.e any angled feature.

I have attached the data before using our customised copy parallel and after using the tool.

I am working on ArcGIS 9.3.1 ESRI product.


texts: Learning .NET for ArcGIS (mainly Desktop) use

texts:
I have basic programming skills (I program simple vba python and lisp scripts and had gone through a java course a long time ago), but it seems to me (judging by ESRI publications) that there is no escaping the need to learn .NET for future relevance in the GIS job market (even as an analyst).
I've gone through a basic c# course and can program (regrettably) quite poor programs, but I have no idea how to link the c# knowledge to actually programming for ArcGIS in order to progress.

I know the differences aren't big, but does ArcGIS use .NET or C#?

I know the best way to further my experience in .NET is to program more and more but are there any good .NET programming books or courses I can learn techniques from?

In your opinion, what is the best way to learn to program in .NET for ArcGIS? Are there any useful books or websites?

I seek your comments regarding ArcGIS 10 (and will be saving them for later) but also for ArcGIS 9.3.1.


texts: Can't get QGIS plugin translations (i18n) to work!

texts:
It's hard to be specific, but I'm having problems setting up translation files for the plug-in I'm building. The plug-in works fine, but I also want to enable translation of the GUI.

I've tried to Google a lot, but any instructions I find are poor in details. It also seems some other plug-ins in the repo don't get translated either, so it may be a bigger problem?

I'm running Ubuntu and QGIS 2.10. Plug-in is generated from the Plugin Builder, and modified in a text editor.
Ts files are built using the "makefile" and edited in the "i18n" folder. Ts/qm files get deployed correctly, but will not translate the GUI! The main python file have text to be translated on the format "self.tr('text to translate')", which calls the "QCoreApplication.translate" function so it seems to me it all should work...

I've tried locale settings in QGIS for a number of common languages as well with very mixed results. Not sure I'm alone with this problem.

Anybody know of a GOOD instruction to get translation of plug-ins to work?


texts: Create polygons from point dataset, where each polygon contains 3 points from the dataset

texts:
I'm using qgis Brighton

I have a dataset of 1500 points. I want to create polygons, where each polygon contains n=3 points from the dataset.

The aim of this is to have clusters of 3 points, which can be used for analysis, qua I can't show information on point level.

I have tried using heatmaps, but I'm not sure if this is the right tool for this analysis.

Here's a picture of my dataset.

Now, I want to create clusters of 3 points like this:

(Sorry for the amateur cutting)

I want to this, so I can make analysis of the clusters (with 3 points) e.g. energy demand, house characteristics, socioeconomics and so on.


texts:
I used ArcGIS 10.3.1 to create shapefile with a text field with 5 length. When I entered a unicode text with 5 characters, ArcGIS lost the last 3 chars (saved only 2 beginning chars).

What's the problem?


texts: Processing Toolbox, avoid using empty layer

texts: I have a process (on the processing Toolbox) I need to reproject and merge two layers (result of previous operations). Sometimes, one of these layers will be empty, in that case it is not necessary to merge (indeed the process crashes if the layer is empty). Is there a way to specify that if the results are empty, don't merge it (but the process still continues)?


texts:
I have built a Script Tool (from a modelbuilder model) in ArcGIS 10.2 which works perfectly fine in ArcGIS.

However, I am getting the following error message when I try to publish this script tool into a geoprocessing service on ArcGIS Server.

I am using in_memory workspace for outputs. Not sure if that is causing the problem.

Error Message -

Code Sample -

> # 5 & 6 - District Query & Final Query
>     if District == "Citywide":
>         gp.AddMessage (("Selected Citywide Stolen Autos..." ))
>         Final_Query = Query + Select_the_Date_Range + " AND " + SQL_Query_Refining_the_Data
>
>     else:
>         Select_the_District = "(DC_DISTRICT = " + District + ")"
>         Final_Query = Query + Select_the_Date_Range + " AND " +  Select_the_District + " AND " + SQL_Query_Refining_the_Data
>
>     arcpy.env.workspace = r"in_memory"
>     gp.AddMessage("")
>     gp.AddMessage("---------------------------------------------------------------------------------------------")
>     gp.AddMessage (("Querying the PIIN table for dates from " + Date_FROM + " to " + Date_TO + " & District = " + District ))
>
>     # Process 1 : Querying the PIIN Table and exporting the selection to a in memory FC
>     arcpy.MakeQueryLayer_management(InputDatabase, "StolenRec_ActiveStol_FC", Final_Query)
>     arcpy.CopyFeatures_management("StolenRec_ActiveStol_FC", StolenRec_ActiveStol_FC_3)
>
>
>     gp.AddMessage (("Making XY Event Layer from Recovered Location..." ))
>     # Process 2: Making XY Event Layer
>     arcpy.TableToTable_conversion("StolenRec_ActiveStol_FC", "in_memory" , "StolenRec_ActiveStol_FC_2", "")
>     arcpy.MakeXYEventLayer_management("StolenRec_ActiveStol_FC_2", "RECOVERED_X_COORD", "RECOVERED_Y_COORD", "Recovered_ELayer",
> "PROJCS['NAD_1983_StatePlane_Pennsylvania_South_FIPS_3702_Feet',GEOGCS['GCS_North_American_1983',DATUM['D_North_American_1983',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Lambert_Conformal_Conic'],PARAMETER['False_Easting',1968500.0],PARAMETER['False_Northing',0.0],PARAMETER['Central_Meridian',-77.75],PARAMETER['Standard_Parallel_1',39.93333333333333],PARAMETER['Standard_Parallel_2',40.96666666666667],PARAMETER['Latitude_Of_Origin',39.33333333333334],UNIT['Foot_US',0.3048006096012192]];-119214200
> -96198500 3048.00609601219;-100000 10000;-100000 10000;3.28083333333333E-03;0.001;0.001;IsHighPrecision", "")
>
>     gp.AddMessage (("Converting XY Event Layer to FC..." ))
>     # Process 3: Converting XY Event Layer to a in memory Point Feature Class
>     #arcpy.FeatureClassToFeatureClass_conversion(Recovered_ELayer, in_memory, "Recovered_FC_Layer","")
>     arcpy.CopyFeatures_management("Recovered_ELayer", Recovered_FC_Layer_2)


texts: Is it possible to install multiple instance of ArcGIS Portal on the same machine?

texts:
For sandboxing purposes, I'd like to be able to test something on a clean instance of Portal, without taking down the portal instance that my team members are using.

Is it possible to create another instance of Portal on the same machine, or am I going to have to spin up a new machine?


texts: Assign geolocation to var before sent to server

texts:
Im a javascript learner and have some trouble with my script i hope someone has time to help me.
In short I want the script to get position of user, save it to a variable and send the position of an object to the server.

I have 3 different objects in my actual script. I want to assign a position to the object before sent to server. I know this can be solved more elegant than having 3 functions (getPositionobject1, 2, 3) that all contain the same if statement. But this is a good starting point.

This script run OK, but sendToServer get called before geometry/position has been assign. Any advice of how i can arrange my code to make sure sendToServer() is called after geometry is assigned?

Here is a fiddle
https://jsfiddle.net/mw7cgwwz/22/

function getPositionobject1() {
    if (navigator.geolocation) {
          var timeoutVal = 10 * 1000 * 1000;
          navigator.geolocation.getCurrentPosition(
            displayPosition,
            displayError,
            { enableHighAccuracy: true, timeout: timeoutVal, maximumAge: 0 }
          );
          sendObjectOneToServer(); //This get called before geometry is defined
}

function displayPosition(position) {
  alert("Latitude: " + position.coords.latitude + ", Longitude: " + position.coords.longitude)
  geometry = (position.coords.latitude + "," + position.coords.longitude).toString()
}

sendToServer() {
    var postData =
        'stuff' + geometry
}


texts: Raster Calculator produces nan results for some multipliers

texts:
I am trying to convert a DEM file from meters to feet using a fairly simple formula.  The range of elevations in the source file is 551 - 1250.

When I multiply by 3.28 I get nan results.
When I multiply by 1.1 I get the expected results (although obviously, this is not a useful conversion to feet from meters)

I didn't have this problem with the essen version of qgis, but I recently upgraded to 2.18.


texts: How do you add this type of north arrow to ArcMap?


texts: NMEA-2000 Message 129540 ("GNSS Sats in View"): SNR field is integer or decimal?

texts:
If you're familiar with the older NMEA 0183, there was a "GSV" record for reporting information about GPS satellites in view, including their SNR (signal-to-noise ratio, also known as carrier-to-noise-density ratio, C/N0), though only as an integer value (between 0 and 99, nominally in decibels, often in dB-Hz). a,b,c

Did this situation change in the newer NMEA 2000?  There is a similar message 129540 about "GNSS Sats in View", but the free N2K Field Descriptions does not specify the fields data types:

  GNSS information on current satellites in view tagged by sequence ID.
  Information includes PRN, elevation, azimuth, SNR, defines the number
  of satellites; defines the satellite number and the information.

I found an official white paper which gives an example for a similar SNR field but for Loran stations instead of GPS satellites, which says:

  Data type: int16;  Range: +/-327.5 dB;  Resolution: 0.01 dB

Does anyone with knowledge of the specification confirm that what applies to message 130054 ("Loran-C Signal Data") is also applicable to message 129540 ("GNSS Sats in View")? Searching further I found some software documentation suggesting that it does:docs

  PGN: 129540 - GNSS Sats in View (...)   Field #8: SNR
                    Bits: 16
                    Units: dB
                    Resolution: 0.01
                    Signed: false

So now the final question is: do real receivers make use of the extra resolution available?  Could anyone with access to a GPS/GNSS unit supporting NMEA 2000 provide a text sample of that message output, please? I heard that some of latest Garmins and Hemisphere receivers work that way.

I have a scientific application that requires SNR values more precise than 1 dB, see Scientific Utility of the Signal-to-Noise Ratio (SNR) Reported by GPS Receivers.


texts: Enabling Enterprise Geodatabase in BigSQL Postgres 9.5.x install and ArcGIS Server/Desktop?

texts:
I have followed Esri's ArcGIS 10.5 documentation multiple times over: http://desktop.arcgis.com/en/arcmap/10.5/manage-data/gdbs-in-postgresql/setup-geodatabase-postgresql-windows.htm

Every time I try to Create a Spatial type for the database utilizing the ST_Geometry.dll that is for pgsql 9.5 placed in the appropriate associated directory, I get the following error:

"You must copy the latest ST_GEOMETRY library to the PostgreSQL software location.  Refer to the ArcGIS help topics for more details.  Connected RDBMS instance is not setup for Esri Spatial type configuration.  Failed to execute (CreateSpatialType)."

I get the same error if I try the enable Geodatabase tool as well.

I am pretty certain there is an issue with esri's ST_Geometry.dll and the BigSQL installation of Postgres, either the type of PostgreSQL installation or because the geodatabase tools are looking for a Postgres install in Program Files and Big SQL installs directly to C:.

Has anyone experienced this and/or know of a solution?  I would really prefer to use the BigSQL implementation if possible, but it doesn't seem like it is going to work.  In the mean time, I'm going to test with different installs/install locations.


texts: Save image within bounding box programmatically

texts:
I want to save an image from QGIS but the save image functionality feels limited.

Is it possible to use the python command line to save an image at a specific zoomed in scale and to set a bounding box over the region I want to print?

For example, I want to save an image at scale 1:1000 000 within specific coordinates as a bounding box.


texts: Seeking geoprocessing tool to find all DTED tiles that intersect an Area of Interest (AOI)?

texts:
Does a geoprocessing tool exist that can find all DTED tiles that intersect an AOI?

I have a area of interest represented by a polygon and DTED tiles (some of which intersect with the AOI).

Is there a tool that I can use to find which tiles intersect?


texts: Do the Sentinel 1 Ground Range Detected imagery available in GEE as ee.ImageCollection("COPERNICUS/S1_GRD") require no correction at all and are ready to use? The 5 pre processing steps mentioned are already done? If not, does GEE purposes a way to perform these 5 corrections?


texts: Analyzing MODIS MCD12Q1 level 3 data in ENVI?

texts: How can i see the classes of the product MODIS MCD12Q1 and analyze them in ENVI exelis? Seadas open them but ENVI not.


texts:
I am trying to calculate mean slope in given polygons. I first interpolate my (x,y,z) data to create a raster file and use the terrain analysis tool to calculate slope.
Then I use the zonal statitistics to obtain slope within each polygon (how-to-calculate-raster-statistics-for-polygons).

I have two issues.

The results are appended to the polygon shapefile. Is there a way to extract the results as a text file ?
Apparently I get only one value (i.e. one mean, one max, one median, etc) while I have 6 polygons. My polygons overlay but I get this problem even if they do not. How can I have separate results from each polygon ?


texts: QgsFileDownloader: download fails every time when run in a plugin

texts:
I'm developing a QGIS plugin that has the functionality to download raster files from a server. As the final step, the user has to name the file and select a file path on their system. After that, a download to said location using QgsFileDownloader commences. The downloader should simply download the file when called to the given location and then delete itself automatically. The following script works perfectly when run in QGIS's own editor:

#file_name houses the path selected by user
file_name, filter = (QFileDialog.getSaveFileName(iface.mainWindow(),
                        "Save the layer", filter='Virtual file (*vrt)'))
url = QUrl("http://dummypath/to/server/file.tif")

# functions triggered when download either succeeds or fails
def downloadFailed(error):
    print(error)
def downloadSucceeded():
    print('success')

downloader = QgsFileDownloader(url, file_name, delayStart=True)

# connecting signals to functions
downloader.downloadError.connect(downloadFailed)
downloader.downloadCompleted.connect(downloadSucceeded)

downloader.startDownload()

However, when I implement that script to my plugin, it errors every time. The code works identically except that it's called from another function and that successful download triggers another function instead of just a print statement. Here's what I mean:

def getData(self):
    # things happen...
    # ...
    self.saveData(url)

def downloadFailed(self, error):
    print(error)

def downloadSucceeded(self):
    # calls another function
    self.addLayerToQgis(self.file_name)

def saveData(self, url, file_format):
    self.file_name, file_filter = QFileDialog.getSaveFileName(self.dlg,
                                "Save the layer", filter='Selected format (*'+file_format+')')

    qt_url = QUrl(url)

    downloader = QgsFileDownloader(qt_url, self.file_name)
    downloader.downloadError.connect(self.downloadFailed)
    downloader.downloadCompleted.connect(self.downloadSucceeded)
    downloader.startDownload()

This always ends in error with the following statement:

['Download failed: Operation canceled']

I don't quite understand why this happens. If I stop the code from proceeding after startDownload call, it works as intended. I mean this:

def saveData(self, url, file_format):
    # code...
    downloader.startDownload()
    raise NotImplementedError

Now the file gets downloaded without a hitch. I need to do things after the download, however. Does anyone know why the file downloader behaves like this?


texts: update cursor datetime if start date greater end date

texts:
I have 2 date fields in my feature class. If the start date is greater than the end date, I want to add 24 hours to the end date. I am getting a SetValue error.

cursor = arcpy.UpdateCursor("PolyFc")
for row in cursor:
    if row.getValue("START_DATE_VAL") > row.getValue("END_DATE_VAL"):
        row.setValue("END_DATE_VAL", row.getValue("END_DATE_VAL") + datetime.timedelta(hours=24))
    cursor.updateRow(row)
del cursor


texts: Automatically get an adequate projection system based on a bounding box

texts:
The website http://projectionwizard.org/ offers a very easy way to get an appropriate projection system based on a bounding box and choosing whether projection criterion should be equal-area, equi-distant or a comprimise.

It is based on this publication that defines some principles on how to select an adequate projection system:

Jenny, B., Šavrič, B., Arnold, N. D., Marston, B. E. and Preppernau, C. A. (2017). A guide to selecting map projections for world and hemisphere maps. In: M. Lapaine and E. L. Usery (eds), Choosing a Map Projection, Lecture Notes in Geoinformation and Cartography (pp. 213–228). Berlin, Heidelberg, New York: Springer. Doi: 10.1007/978-3-319-51835-0_9

I would like to know if there is any similar tool for R (or Python) that can achieve the same. So e.g. based on a shapefile in WGS84 find an adequate projection system for the shape with equalarea as the main criterion.

Relevance: to calculate areas in R without having to manually find a projection system. Currently I only know of the function area in the raster package. However, this function does not use a projected coordinate system. From it's description:

Compute the approximate surface area of cells in an unprojected (longitude/latitude) Raster object. It is an approximation because area is computed as the height (latitudial span) of a cell (which is constant among all cells) times the width (longitudinal span) in the (latitudinal) middle of a cell. The width is smaller at the poleward side than at the equator-ward side of a cell. This variation is greatest near the poles and the values are thus not very precise for very high latitudes.


texts:
I have multistory building plans and each floor located the same place.

I need to move each floors side by side.

for example, 15 floors need to move side to side as like 5 (Column) X 3 (Row)

Any suggestions


texts: OpenLayers show interaction style without hovering

texts:
I have been making a simple Modify interaction in OpenLayers, for example

var dragInteraction = new ol.interaction.Modify({
    features: new ol.Collection([feature]),
    style: dragStyle
});

Can I force the interaction style to stay visible without having to hover over it? I am adjusting the anchor for each marker, which makes it difficult to know where to hover the cursor to actually drag the marker.


texts: Google Earth Polygon surface area to follow the curvature of the earth model

texts:
Question:
When drawing a polygon in Google Earth with latitude, longitude and elevation values at each node, how do you ensure that the surface area of the polygon also represents the desired elevation in relation to the earth's surface (interpolated between the nodes)? And, not have the surface area tangent to the nodes.

At the moment, when I draw a polygon over the earth at 800000 m absolute, it looks like this:

I would like the surface to be the same height as the perimeter of the polygon or to have the surface interpolated between the nodes. So, for example, if I had an inclined plane stretching over the ground where one end is 50 m elevation absolute and the other end is 500 m elevation absolute, the area enclosed in the polygon would represent the interpolated z value. At the moment, it is clearly not the case as can be seen in the picture.

I suspect it would be a setting in Google Earth or there is something in the KML code that can be manipulated..


texts: Saga Catchment Area not giving expected (useful) output

texts:
I am identifying a catchment network from a DEM in QGIS3 (v3.0.3) (on a mac  v10.13.4), however the step to identify the Catchment Area is not working.

I have an input DEM, with layer CRS set to UTM55S.  I am using Method 0 (Deterministic).  The output layer for the Catchment Area has no identifiable channel networks.

I am using a hydrologically corrected DEM (provided by Geosciences Australia), clipped to the area of interest.  I have saved the output raster in geotiff format and am working from this (could this be the problem?).

I am using SAGA processing toolbox.

I first run Fill Sinks, which gives an expected output - I check the difference between the original and Filled DEM and nothing stands out.  I then run the Catchment Area using the default options and don't see what I expected.

Running the Channel Network algorithm is not giving me what I would expect - the network is entirely disconnected.  I tried a few different initiation thresholds (greater than option - 0.1 -> 1e-6). These are subbasins of major rivers, I would expect to see a dominant stream network.

I have read a lot of different responses around this topic, but mostly they deal with wrong initiation thresholds (so blank output if too high), wrong projections (but UTM is suggested as appropriate), or CRS not being consistent (I have set CRS through Set Layer CRS to ensure this isn't the problem).


texts: Deleting 'garbage' nodes after polygon dissolve

texts:
I work with QGIS 2.10.1 on Windows 7 (x64). I have dissolved polygons from an attribute value and the polygons geometries being not perfectly joined together, I get some 'garbage' nodes I need to delete. I have tried to edit the layer and using nodes tool, selecting and deleting the nodes, but it seems impossible to do (nodes are not deleted at all!). Any idea on how to process this kind of clean-up?

EDIT : you can download the shapefile here for testing.


texts: MapServer & Raster picking / Elevation data

texts:
I have some GeoTIFFs with elevation data, displayed in OpenLayers or Leaflet via WMS from MapServer. I need to be able to request (from the client) the elevation at some point clicked by the user.

Does MapServer provide this functionality? If yes, how? If no, do you know a way to do so?


texts: Image collection containing multiple relative orbit numbers in Google Earth Engine

texts:
I want to select images acquired from two different relative orbits in Google Earth Engine. What i tried so far is:

var collection_orbit = collection
    .filter(ee.Filter.eq('relativeOrbitNumber_start', 54))
    .filter(ee.Filter.eq('relativeOrbitNumber_start', 113));

which throws the error "Image Collection Empty" while using only one of the lines above works well.

How can i add images from two different relative orbits to my image collection?


texts: Books or other Resources on Business Benefits from GIS?

texts:
I believe Top unsung business benefits for GIS is an excellent source of info.

My question expands on that for any industry books related to this subject.  I know there have already been several postings for recommended books but I am trying to find books that explicitly explain the value of GIS to a company.

In short, I am looking for books that will clearly explain the business benefits (financial and otherwise) that GIS can provide for a company.

Are there are any resources on the web for this information?


texts: I have a table of attributes with polygons, from the area data I sort from largest to smallest, and when using the command @row_number, he writes the order in which the polygons were created. I want to know how to change the classification of features to the size of the area and then make the count from largest to smallest?


texts: I downloaded a new shapefile and is displayed without any value in the attribute table, do you know what could be the reason?


texts: Raster of multiple polygon associating values in R?


texts: Filtering zip code shapefile by state/other zip code shapefiles?

texts:
I am looking to map one state based on zip codes (I have associated quantitative data that I am planning to join in QGIS).

I see that the Census TIGER shapefiles have an option for zip codes, but I don't see any associated state-based column that would allow for me to select a single state.

Can anyone could confirm whether it has state information, or has another source I can refer to, or does anyone have a spreadsheet that separately links zip code to state?


texts: How to see only those incidents that cannot reach facilities in ArcGIS Closest Facility analysis

texts:
I have 373 'incidents' and need to map their routes to their closest hospitals ('facilities') using a flood zone as a 'polygon barrier'

ArcGIS can only route 333 of the 373 because several 'incidents' are within the flood zone and the rest cannot reach their 'facility' because all possible routes have been cut off by the flood zone.

ArcGIS differentiates between those 'incidents' that are within the flood zone, but treats all other 'incidents', even those that cannot be routed, as the same (see image). When I check the Attributes there is no info to differentiate them.

How can I differentiate between those 'incidents' that can reach their 'facility' and those that cannot reach their 'facility' because no route is available?


texts:
I'm trying to use Geoserver Parameter SQL views, but I'm having trouble getting geoserver to recognize the parameters (and I'm not able to follow the documentation - it's not very clear to me). Here's the code I'm working with so far, does anyone have any suggestions?

SELECT *
FROM Table1 INNER JOIN Table2
ON Table1.ID = Table2.ID
WHERE Table1.Date = %Date%

I'm trying to get geoserver to recognize each "Date" as a parameter - there are about 25 total dates, and I want the user to be able to select a date (or multiple dates) for display. The dates are in yyyy-mm-dd format.

Thanks for your help.


texts:
I just installed QGIS Madeira 3.4 on my WINDOWS Vista PC (64bits).
But I can't work with it, since I keep getting the same error:

And just after, this one:

Should I uninstall it?
What could I do in order to work with QGIS?


texts: How to set up proxy for Geoserver GetFeatureInfo call?

texts:
I have problems with cross domain, I call GetFeatureInfo in Leaflet API. I was looking for solution and saw that proxy could solve that problem. Can you help me how to set it up? I have geoserver run on Tomcat7 (ubuntu).

Thanks.


texts:
I am building a report in QGIS that handels a layer with quite a large attribute table (around 200 attributes). But most of the attributes are boolean values and I only need (want to display) 'true' values in my report.

Is there a way, possibly with expression builder, to filter for boolean 'true' values?

EDIT:

This is a section of the attribute table, and i only want to display the attributes with 'true', but that differs from case to case, so i have to add all of them to the report, and want to filter them while generating the report.


texts: Making ModelBuilder tool for updating multiple fields from selected data?

texts:
Is there a way I can make a tool (in ModelBuilder) that will automatically set fields to a certain value?

For example:

I manually select a bunch of polygons in ArcMap, run a query on the selected polygons, and multiple fields are updated with the desired value.


texts: How to avoid manually typing coordinate information in CreateCircle function to select points in Mapinfo

texts:
I have a bunch of data points, each of them have a coordinate information (longitude, latitude).

I am using CreateCircle( x, y, radius) function to select points within a distance to the chosen point.

But I need to manually update the x, y (coordinate) information for each selection, which is impossible since I have thousands of points.

Does anyone know if there is a way to update the coordinate information in the CreateCircle function? Not manually typing it everytime?


texts: ArcMap Near Tool doesn't find nearest feature for some input points on the line

texts:
Using ArcMap's near tool, I have 'projected' the vertices of a trail onto a corresponding line (river), in order to find the closest point on the river to  each trail vertex. However, as you can see in the screenshot, not all the trail points (from the green line) have found corresponding near points on the river... I tried first with no search radius and again with a search radius of 2000 m (the distance between the river and the trail in the area is under 300 m).

How can I find all the corresponding near points along the river?


texts:
I would like to get WMO weather stations observations daily reports from the web in file or web-service form.
According to WMO Resolution 40 -info- this data should be availiable in text form with SYNOP codes.

Does anyone know a place to look for this data?

I've looked through NCEP and NASA sites, and didn't find the dataset itself.

I've found only http://www.ogimet.com site, that provides data in html and txt form, but you have to specify a WMO station index or a region - so no service for global coverage (all WMO stations). I suppose this site gets information from the resource, that I'm looking for.


texts: Highlighting using distinct color where features overlap in QGIS?

texts:
I have two species on my map and right now they are different colors. However, what I want to happen is when they overlap in the same spot, the color is a distinct difference from the other two colors. Only thing I've thought of is to change the transparency so the colors mix together. I'm using the stable version of GQIS.

Do you know a better way?


texts: I'm new to Qgis. Can I use similar items in the SAME layer, with each distinguished by color? For example, different crops in one layer (corn, wheat, barley...) while, a different layer would distinguish among different animals (cattle, pigs, poultry...)


texts: Copying Symbology between layers with different numbers of features

texts:
I am working on copying a layers symbology to a number of other layers which have the same geometry, table structure etc, etc but not the same number of features.

I have a state wide layer and I have written some code that selects out the different regions within the state layer based on a definition query and exports the features to shapefiles and then adds these shapefiles to ArcMap.

I am now trying to copy the symbology from the state wide layer to these other layers. However when I copy the symbology it copies the symbology that's set for the entire state.

In VBA, is there a way to remove these extra values so I end up with the same symbology as the state layer but that is specific to the layers new extent?


texts:
I need to calculate the mean slope of roads between nodes with QGIS. I have a DEM and a vector (PostGIS) with the roads, nodes are additionally displayed as point shape and also contain an elevation field derived from DEM:

I am looking for the mean slope between those nodes which are connected with a line.

As there are many roads to calculate (see the picture below), I need a workflow which processes all calculations at once.Either output as raster (with slope values on the roads and NoData for others) or shapefiles containing the the slope are useful for my purpose.

How can I do this?

I have a ArcGIS solution in mind (3D Analyst / interpolate shape), but I'd prefer a QGIS / open source solution!

I tried the Profile Tool Plugin, which only processes one line at once. Further, I was thinking of using a slope raster and failed because of the direction of the slope (contour line streets).

I was looking for other approaches, but it didn't work or I didn't understand them, e.g.:

calculate slope for a road network
Computing the slope of a road using DEM
Slope, in which direction?
Calculating slope direction in QGIS?
Determining gradient of road segment
Calculating average slope from line using QGIS?
http://youarealegend.blogspot.ch/2012/07/computing-roads-slopes-with-qgis.html


texts: Capturing accurate GPS points using mobile devices

texts:
I'm  currently working in Vietnam as a climbing guide, but am putting my GIS skills to use by doing a side project for the company involving the mapping of climbing routes. Along with another employee who is an app developer, our goal is produce a digital guidebook for the area.

For this project to work for me, I need to be able to take fairly accurate GPS points, as there are features just a few meters apart from each other that need to be distinguishable, but obviously survey grade equipment is not and option for me here. I've been using the app GPS Averaging to take several hundred points for each location.

Despite taking these averages, I'm having very little luck achieving accuracy in my output points, and I suspect my precision is low as well because the in app reported error gets very low (0.8m) after about 20 collected points, then goes up to about 3m after 200-300 collected points. However when I import these points, the actual errors appear to be well beyond 3 meters.

What can I do to improve GPS point accuracy from a mobile device?

I've currently only worked on this with android devices, but have access to ios devices too if they provide superior GPS chips.


texts: Change sde connection parameters with ArcPy

texts:
I know you can create a database connection with Arcpy using the CreateDatabaseConnection_management function.  I am wondering if anyone knows a way in which to change the parameters of that file.

For instance, say you want to use the same .sde connection file, however you want to log into a different database, or connect with a different user.

Essentially is there a way to access the properties of it as if you just right click on the connection in catalog and select connection properties.

I'd rather have my script modify the settings than delete that sde file and recreate it with the same name but different parameters.


texts: Nearest point to 30 million points using R

texts:
I have two point datasets. One with 30 million records and the other with 1.5 million. For each of the 30 million points, I would like to calculate the nearest point(s) from the 1.5 million dataset to them.
I have the following, very simple, code that would do the job if I had a spare 6 months. Is there a quick and efficient alternative way to do this?
library(data.table)

DT1 <- data.table(DT1_id = 40000001:70000000, DT1_x =sample(450000:550000, size = 30000000,replace = TRUE),DT1_y =sample(250000:350000, size = 30000000,replace = TRUE) )
DT2 <- data.table(DT2_id = 8000001:9500000, DT2_x =sample(450000:550000, size = 1500000,replace = TRUE),DT2_y =sample(250000:350000, size = 1500000,replace = TRUE) )

DT4 <- as.data.table(NULL)

for (n in 1:nrow(DT1)){
DT3 <- cbind(DT1[rep(n,nrow(DT2))],DT2)
DT3[, DIST := sqrt((DT1_x - DT2_x)^2 + (DT1_y - DT2_y)^2)]
DT4 <- rbind(DT4,DT3[ , .SD[DIST==min(DIST)], by = DT1_id])}

EDIT
Revised DT2 example dataset that more closely resembles my real data:
DT2 <- data.table(DT2_x =sample(450000:550000, size = 1200000,replace = TRUE),DT2_y =sample(250000:350000, size = 1200000,replace = TRUE))
DT2 <- rbind(DT2,DT2[sample(.N, 400000)])
DT2 <- data.table(cbind(DT2_id = sprintf("%s%0*d", "UID_", 7, 1:nrow(DT2)),DT2))


texts: Update point corrdinates in the attribute table after using using edit to move points in ArcGIS

texts: I have a shapefile of point locations that I edited after the original file was created, however the the attribute table was not updated with the new coordinates. So I am trying to figure out how to update the attribute table with the new point location coordinates. I am using ArcGIS Desktop 10.6.1.


texts: OpenLayers and TinyOWS WFS won't load - 'InvalidParameterValue'

texts:
I've got a TinyOWS on localhost that works fine (both through browser GetFeature and QGIS). But my OpenLayers code fails. I've checked here:
Display a WFS layer with OpenLayers and here : How to add a simple WFS layer from GeoServer to OpenLayers map?
and other similar questions to no avail. Can anyone help?

Here is the code:

    wfs = new OpenLayers.Layer.Vector("Editable Features", {
    strategies : [new OpenLayers.Strategy.Fixed(), saveStrategy],
    projection: new OpenLayers.Projection("EPSG:4326"),
    protocol: new OpenLayers.Protocol.WFS({
        version: "1.1.0",
        srsName: "EPSG:4326",
        url: "http://localhost/cgi-bin/tinyows",
        featurePrefix: "tows",
        featureNS :  "http://www.tinyows.org/",
        featureType: "trad",
        geometryName: "wkb_geometry",
        schema: "http://127.0.0.1/cgi-bin/tinyows?service=WFS&version=1.1.0&request=DescribeFeatureType&Typename=tows:trad"
    })
});

map.addLayers([gphy, wfs]);

Here is the xml according to Firebug:

<wfs:GetFeature xmlns:wfs="http://www.opengis.net/wfs" service="WFS" version="1.1.0" 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">

Here is the post response according to firebug:

<?xml version='1.0' encoding='UTF-8'?>
<ows:ExceptionReport
  xmlns='http://www.opengis.net/ows'
  xmlns:ows='http://www.opengis.net/ows'
  xmlns:xsi='http://www.w3.org/2001/XMLSchema-instance'
  xsi:schemaLocation='http://www.opengis.net/ows http://schemas.opengis.net/ows/1.0.0/owsExceptionReport.xsd'
  version='1.1.0' language='en'>
<ows:Exception exceptionCode='InvalidParameterValue' locator='request'>
<ows:ExceptionText>XML request isn't valid</ows:ExceptionText>
</ows:Exception>
</ows:ExceptionReport>

Here is the TinyOWS config file:

<tinyows online_resource="http://127.0.0.1/cgi-bin/tinyows"
     schema_dir="/usr/local/tinyows/schema/">
<pg host="127.0.0.1" user="postgres" password="****" dbname="****" port="5432"/>
<metadata name="TinyOWS Server"
       title="TinyOWS Server - Demo Service" />
<layer retrievable="1"
    writable="1"
    ns_prefix="tows"
    ns_uri="http://www.tinyows.org/"
    name="trad"
    title="Träd" />
</tinyows>


texts: Calculating percentage overlap in graphical modeller

texts:
Can anyone help me calculate the percentage overlap of a polygon with multipoint attributes (I.e fragmented pieces of polygon overlapping) in graphical modeller?

I need to calculate the area of my original polygon, then the area of the overlapping polygon (all the parts/attributes), then add A final filed to calculate the percentage overlap.

In an ideal world I would also like to count the sum total over overlapping parts as another total.

The final end product should end up on a single line with the attributes from the original polygon, the percentage area overlapped and the number of overlapping attributes.


texts:
Using QGIS 3.0 I try to identify features by user-form (.ui). The opening form-window is very small and hides nearly every content of the form. It has to be expanded every time manually by cursor to the full form size to see all contents of the form.
What settings have to be changed that the form opens in full-size ?


texts: What projection does ST_Area(geom, true) use?

texts:
ST_Area documentation states that "area is determined on a spheroid with units in square meters" when using geography as argument, or when setting the 2nd parameter as "true".
Usually when calculating areas we can reproject a geometry to an "Albers Equal Area", and these projections may have different source latitude and longitude depending on the region we are working. For instance +proj=aea +lat_1=-2 +lat_2=-22 +lat_0=-12 +lon_0=-54 +x_0=5000000 +y_0=10000000 +ellps=GRS80 +units=m +no_def.
On the other hand ST_Area does not state any source lat or lon or skews. Are those values dynamic? Or are they "borrowed" from the projection of the geometry? (in case, for srid 4326 would be lat=0 lon=0 x=0 y=0)


texts: Order of geometries in GeometryCollection

texts:
Does the order of geometries/features in GeometryCollection/FeatureCollection matter?
There is a point in GeoJSON RFC (last point) about order but I could not get fully:

  The word "Collection" in "FeatureCollection" and
        "GeometryCollection" does not have any significance for the
        semantics of array members.  The "features" and "geometries"
        members, respectively, of these objects are standard ordered JSON
        arrays, not unordered sets.

Are following two geometries to be considered equal or not

{
  "type":"GeometryCollection",
  "geometries":[
    {"type":"Point","coordinates":[1.0,1.0]},
    {"type":"Point","coordinates":[2.0,2.0]}
  ]
}

{
  "type":"GeometryCollection",
  "geometries":[
    {"type":"Point","coordinates":[2.0,2.0]},
    {"type":"Point","coordinates":[1.0,1.0]}
  ]
}

If collection should be ordered list then what geometry operation on collection may give different result in case geometries' order is reversed?


texts: book/tutorial on using rpy2 in python for kriging?

texts: I am proficient in Python, but have no knowledge of R. Can you suggest a good tutorial in using rpy2 or any other R library to perform Block/Cokriging "Kriging"


texts:
I have a raster in one coordinate system (SC State Plane NAD83 NSRS2007) and need to project it to another (SC State Plane NAD83 HARN). When I run Project Raster, it fails silently: the dialog box output says "Succeeded," but no output is produced. Also, clicking "Close" on the dialog box brings the tool dialog back up, as it would if a tool failed to complete.

Without an error message, I am struggling to understand what the cause of the problem is, and how to circumvent or resolve it.


texts:
A polyline layer's (green) length is being calculated via $length. Result: ~100 meters (EPSG 31466).

A buffer of 5meters is created and via Voronoi Polygons (and several other steps) a centerline is created (red). Result: ~ 200 meters (EPSG 31466).

That simply makes no sense at all if you look at the screenshot. It rather should be arround 110 meters, estimated visually. I tried merging all segments but the result stays about the same.

The background behind this situation is a sum length completeness comparism of two road networks. Because of many withstanding end pieces I want to calculate a factor I can use to divide a total network length with in order to correct the centerline length.


texts:
I am facing the issue while I am trying to get all the features back in 'Attribute Table'.
What I am doing is.

Selected the layer from Map
Performs my Edits on few features (changing filed values)

(After I perform my changes I take

While I again want to get all my features back I am getting error.

Exception from HRESULT: 0x8004151C

I tried with below code :

In both the case it's giving error.

///////1st//////////

IFeatureClass pFc = ilayer.FeatureClass;
IFeatureCursor pFCursor;
IFeature pF;

queryFilterlast.WhereClause = "SHAPE.LEN IS NOT NULL";
pFCursor = pFc.Search(queryFilterlast, false);
pF = pFCursor.NextFeature();

while (pF != null)
{
   pF = pFCursor.NextFeature();
}

/////////2nd//////////////

//initialize queryfilter
                queryFilterlast.SubFields = "*";
                queryFilterlast.WhereClause = "SHAPE.LEN IS NOT NULL";

                //Using a query filter to search a feature class:
                featureselectionlast.SelectFeatures(queryFilterlast, esriSelectionResultEnum.esriSelectionResultNew, false);

                //output features assign into featureset
                featuresetlast = featureselectionlast.SelectionSet;

                if (featuresetlast.Count != 0)
                {
                    featuresetlast.Search(null, false, out icursor);
                    //featureCursor = icursor as IFeatureCursor;
                    IFeatureCursor featureCursorL = icursor as IFeatureCursor;

                    //Get the output features
                    IFeature feature = featureCursorL.NextFeature();

                    //Loop through output features
                    while (feature != null)
                    {
                        feature = featureCursorL.NextFeature();
                    }
                    //Release the cursor
                    System.Runtime.InteropServices.Marshal.ReleaseComObject(featureCursorL);
                    featureCursorL = null;
                }


texts: How to determine Kc value and apply it to r.sun?

texts:
As outlined in this link ( r.sun delievering unrealistic Values ) to get realistic values of kilowatt / meter-2, you must multiply the results by a "clear sky factor (Kc)". The link comments recommend this website, https://eosweb.larc.nasa.gov/cgi-bin/sse/subset.cgi?email=skip@larc.nasa.gov, but there is no option to determining Kc. I have tried to follow the directions on determining the ratio of overcast situation to clear sky but neither are options on the NASA page. Also, I don't know how to multiply the output raster by whatever the Kc value might be. Can anyone assist?

From the original answer:

  The r.sun mode 2 computes just the clear-sky irridiation. To get the
  "real-sky" irridiation you have to multiply the result with a clearsky
  Factor (Kc). Kc is defined as the ratio of messured Irridation on
  Ground under overcast situation and the clearsky amount messured via
  satellite (both for horizontal orientatian) and has values from 0-1.0
  . Data to calculate the Kc-Factor can be requested here:
  https://eosweb.larc.nasa.gov/cgi-bin/sse/subset.cgi?email=skip@larc.nasa.gov
  After multiplying the results with Kc the results seem to be
  realistic.


texts: What kind of spatial index does arcgis use?

texts:
My first master's research topic was about spatial indexes.  I learned about R-Trees and R+-tree, R*-trees and other various but very similar topics.  This may not be published information, but what spatial index type does ESRI use?  As a secondary question can this type be changed and how?  There are some situations where a different type of index would boost performance.

In Oracle's case it can use the ST_Geometry spatial data type which can be indexed to either an R-tree or a quad tree or both: http://docs.oracle.com/html/A88805_01/sdo_intr.htm

Other databases use binary storage mechanisms which are non-spatial by default and do not seem to use the spatial indexing available.  In SQL 2005 and lower spatial datatype didn't exist, but ESRI was still able to use SQL databases.  According to this case study SQL 2012: http://www.microsoft.com/casestudies/Case_Study_Detail.aspx?CaseStudyID=710000000547 it uses the automated spatial performance type, but doesn't list what kind of R-Tree it is.

Bonus points for reading material answering this question.

Thanks.


texts: Is it possible to re-point operational layers ArcGIS Mobile projects v.3.0

texts:
We have an ArcGIS mobile project with two operational layers for data collection. The layers source from a GIS server on the cloud, running on top of an SDE. The database has a feature dataset, and the editable feature classes within. The editable feature classes have a lot of fields, and the mobile project has all the layers' and their fields' settings.

We want to deploy projects for a number of different regions, thus've we built a separate feature dataset and its child feature classes for each region, and published the GIS services for the mobile devices to interact with. We want to replicate the mobile project and have one project for each region.

Question: is it possible to make a copy of a mobile project and re-point the operational layers to a different service and still keep the settings (fields are the same as all feature classes have the same schema), to avoid having to configure the layers by manually doing all the ticks again? Or, is there a way to have the settings in a text file to read from? We're using Mobile Project Center 3.0.

Thanks


texts:
I've been trying to make an animated path using d3, JSON, and a list of coordinates.

This is the result in JSFiddle:

http://jsfiddle.net/5m35J/92/

As you can see, there is a red dot in the corner, but the animation is stopped. I do not know if the geographic information is correctly loaded.

$(document).ready(function () {

var json, map,
path, vis, xy,
duration, offset, origin_x, origin_y, len, group, circle;

json = {
    "type": "FeatureCollection",
        "features": [{
        "type": "Feature",
            "geometry": {
            "type": "LineString",
                "coordinates": [
                [-100.934121, 21.159216],
                [-100.807746, 20.527743],
                    [-101.427837, 20.951976],
                    [-101.257851, 21.018937],
                    [-101.189444, 19.768333],

            ]
        }
    }]
};

// Draw a GeoJSON line on the map:

map = $('#map');
xy = d3.geo.mercator().scale(12000).center([21.159216, -100.934121,]).translate([630700, 401100]);
path = d3.geo.path().projection(xy);

vis = d3.select("#map")
    .append("svg:svg")
    .attr("width", 960)
    .attr("height", 600);

vis.append("svg:g")
    .attr("class", "route")
    .selectAll("path")
    .data(json.features)
    .enter()
    .append("svg:path")
    .attr("d", path)
    .attr("fill-opacity", 0.5)
    .attr("fill", "#fff")
    .attr("stroke", "#333");

// Draw a red circle on the map:

group = vis.append("svg:g");

var targetPath = d3.selectAll('.route')[0][0];
var pathNode = d3.select(targetPath).selectAll('path').node();
var pathLength = pathNode.getTotalLength();

circle = group.append("circle")
    .attr({
    r: 10,
    fill: '#f33',
    transform: function () {
        var p = pathNode.getPointAtLength(0)
        return "translate(" + [p.x, p.y] + ")";
    }
});

// Animate the circle:

duration = 10000;
circle.transition()
    .duration(duration)
    .ease("linear")
    .attrTween("transform", function (d, i) {
    return function (t) {
        var p = pathNode.getPointAtLength(pathLength*t);
        return "translate(" + [p.x, p.y] + ")";
    }
});

});


texts: QGIS: Displaying total length of lines in print composer legend

texts:
In the print composer (also in the layer selector) I can choose to automatically display the amount of objects in that category in square brackets. That is useful for point features/layers. But in the case of line layers I would like to show the total length of all lines in that category (sum) in just the same fashion.

It is possible to get that info using the GroupStats plugin. But I am looking for a way to automatically keep my legend updated with that info. Just as neat and uncomplicated as with the square brackets for the total count of objects.


texts: Failed to add a TIFF file as a raster layer in QGIS

texts:
I am new to QGIS. The system of my computer is 32-bit operating system, with x64-based processor. My QGIS version is 3.12.1. The TIFF file that I am using is 322 MB large. But every time when I add this file as a raster layer, my QGIS cannot load it. The status bar at the bottom is always loading. The task manager of my computer shows that the QGIS software is using 99% of CPU. I don't know how to solve this problem. I have tried some samples of TIFF files on internet, and they works.

Actually, I successfully opened the TIFF file (the 322MB one) once. It was when I first downloaded it on a website. But later, I could not open it via QGIS anymore.


texts: Obtaining flow accumulation from negative elevations using ArcGIS Desktop?

texts:
I have bathymetric data extracted from INFOMAR ArcGIS grids.

Based on the data I am wondering if the fill, flow direction and flow accumulation tools work with negative elevations or if I have to use the raster calculator to change them to values lying above the sea level.


texts:
I've grouped some vector layers (from SHP files) in GeoServer (they are EPSG:2039) and created cache for them (GeoWebCache lets you choose EPSG:4326 OR EPSG:900913)
Once I've fired my App it didn't seem to read the cached tiles from the WMS (although they were created in the hard drive folder)

Could it be that because of the map projection the tiles are not being read ? (In the hard drive they are stored under EPSG_4326 folder)

This is how I init my map :

function init() {
        // if this is just a coverage or a group of them, disable a few items,
        // and default to jpeg format
        format = 'image/png';

        var bounds = new OpenLayers.Bounds(
            155000, 600000,
            246000, 745000
        );
        var options = {
            maxExtent: bounds,
            maxResolution: 566.40625,
            projection: "EPSG:2039",
            units: 'm'
        };
        map = new OpenLayers.Map('map', options);

        // setup tiled layer
        tiled = new OpenLayers.Layer.WMS(
            "background", "http://localhost:8080/geoserver/MyWorkspace/wms",
            {
                LAYERS: 'background',
                format: format,
                tiled: true,
                tilesOrigin: map.maxExtent.left + ',' + map.maxExtent.bottom
            },
            {
                buffer: 0,
                displayOutsideMaxExtent: true,
                isBaseLayer: true,
                yx: { 'EPSG:2039': false }
            }
        );

        map.addLayer(tiled);

        // build up all controls
        map.addControl(new OpenLayers.Control.LayerSwitcher());
        map.addControl(new OpenLayers.Control.PanZoomBar());
        map.addControl(new OpenLayers.Control.Navigation());
        map.addControl(new OpenLayers.Control.Scale());
        map.addControl(new OpenLayers.Control.MousePosition());
        map.setCenter(new OpenLayers.LonLat(200300,745000),7);
    }


texts:
Can someone help me understand what the last 2 lines of this code does:

import arcpy
arcpy.env.workspace = "c:/esripress/python/data/exercise07"
arcpy.env.overwriteOutput = True

copy = arcpy.CopyFeatures_management("airports.shp","Results/airports.shp")
fc = "Results/airports.shp"

cursor = arcpy.da.UpdateCursor(fc, ["STATE"], ' "STATE" <> \'AK\'')
for row in cursor:
    row[0] = "AK"
    cursor.updateRow(row)
del row
del cursor

I understand that the loop function goes through each record that does not have a value of 'AK'  and gives that record a value of "AK".  But what I don't understand is what the del row and del cursor are meant to do.


texts: ArcGIS Javascript API FeatureLayer not available immediately upon show click event

texts:
I am adding a FeatureLayer in my application as follows:

featureLayer00 = new FeatureLayer("http://webadaptorserviceurl/MapServer/0",{
          mode: FeatureLayer.MODE_ONDEMAND,
      id: "0",
          outFields: ["*"],
          infoTemplate:template
        });

Then I toggle this layer with show/hide in a common.js file as follows:

...

if (visible.indexOf('0') > -1) {
            featureLayer00.show();
        }
        else {
            featureLayer00.hide();
        }
 filteronAttribute();

...
It works fine until I get to filteronAttribute().  In the common.js file my filteronAttribute() loops through the layer and shows/hides each stat layer point as follows:

...
statLayer = featureLayer00.graphics;
dojo.forEach(statLayer, function (g) {
     if (sza.indexOf('all') > -1) {
         g.show();
     } else if (sza.indexOf('100') > -1 && g.attributes.hosp_beds < 100) {
         g.show();
}); ...

If the layer is on and I call filteronAttribute(), all works as expected.  However, if the layer is off and I turn it on, the statLayer in filteronAttribute() is empty, even though it was just set to show.

Is there something that has to happen to the featureLayer after it is switched from hide to show before it is available in other functions?

Everything works fine if the layer is visible and I filter. However, if the filter is set and the layer is switched from show to hide, the statLayer in filteronAttribute is always empty, even though it is called after the show method on the layer.  Any ideas?


texts: Connecting points with lines and saving starting and ending points in new layer

texts:
I have 3 points in one layer as shown in the first image. I am trying to connect a building points with junction point by 2 lines. For demo, I have manually drawn a lines in image 2. Here, I am trying to achieve two things:

Connect points (buildings) with Junction by lines.
Save the data of start point and end point of connecting lines in seperate  table (layer) as shown in image 2 (rightside).

Overall, there will be 3 features in points table and 2 features in line table.
What I have tried and search?

I have tried points to path using Menu Processing / Toolbox / points to path. It only connect one building with junction point with no attributes (startpoint,endpoint) saved in the table.
I have also searched for the exisiting answers 1, 2, 3, 4 but I couldn,t able to solve the problem.


texts: Why does geoprocessing service tool Table to Excel always put output in scratch folder and not folder specified?

texts:
I have the output set to a registered folder that arcgisserver has permissions to access, it works fine when I run it locally but when I run it through the web it always sends my output file to a scratch workspace instead of the specified folder.

Does anyone else have this problem?


texts:
I'm trying to use an expression to label a layer by using the "name" field but only if it has certain entries from the "type" field. The expression I've tried is:

NAME IF  "TYPE" IN  ('RIV-MAJ', 'LAKE-PER')

Can someone suggest how I would go about this?


texts:
I am trying to figure out how to use GDAL (ogr2ogr) to export spatial data from SQLServer to KML. The standard KML driver works, but it is not very robust. LIBKML has so much promise (layers, styling, etc), but I cannot find any examples of how to make use of this added functionality.

For example, this works:

ogr2ogr.exe -f LIBKML z.kml "MSSQL:server=box;database=main;trusted_conneciton=yes;tables=dbo.z" -where "state='DE'"

It exports a POINT from my table. I would like to customize labels, icons (color, scale and image) and balloon text. I would also like to create folders/layers.

Ideally, I'd love to use columns like "icon_scale" and "layer_name" in my source tables to control this customization ...

Can you provide an example?


texts: How to calculate accuracy of latitude/longitude

texts:
I'm updating a Windows Mobile application to capture location information using the built-in GPS hardware. The Microsoft API provides a method GPSGetPosition that returns the following struct (ref:GPS_POSITION):

typedef struct _GPS_POSITION {
    DWORD dwVersion;
    DWORD dwSize;
    DWORD dwValidFields;
    DWORD dwFlags;
    SYSTEMTIME stUTCTime;
    double dblLatitude;
    double dblLongitude;
    float  flSpeed;
    float  flHeading;
    double dblMagneticVariation;
    float  flAltitudeWRTSeaLevel;
    float  flAltitudeWRTEllipsoid;
    GPS_FIX_QUALITY     FixQuality;
    GPS_FIX_TYPE        FixType;
    GPS_FIX_SELECTION   SelectionType;
    float flPositionDilutionOfPrecision;
    float flHorizontalDilutionOfPrecision;
    float flVerticalDilutionOfPrecision;
    DWORD dwSatelliteCount;
    DWORD rgdwSatellitesUsedPRNs[GPS_MAX_SATELLITES];
    DWORD dwSatellitesInView;
    DWORD rgdwSatellitesInViewPRNs[GPS_MAX_SATELLITES];
    DWORD rgdwSatellitesInViewElevation[GPS_MAX_SATELLITES];
    DWORD rgdwSatellitesInViewAzimuth[GPS_MAX_SATELLITES];
    DWORD rgdwSatellitesInViewSignalToNoiseRatio[GPS_MAX_SATELLITES];
} GPS_POSITION, *PGPS_POSITION;

So I can determine the latitude and longitude but I would like to be able to calculate the accuracy of the result. The application runs on a Trimble Juno and the GPS Controller application on the device provides an accuracy value in metres. Is this something I can determine from the values in the struct or is this something specific to the hardware that would require a native function call? Is the accuracy value based on heuristics or is it a common calculation across different hardware?


texts: How can I cluster all points in a postgis table together based on their distance from each other?

texts:
I have a PostGIS (in postgres 9.3) table which has Points. I want to figure out the "clusters" in the data. If 2 points are < X distance apart, they should be added to a cluster, and this should be done recursively. In essense, a cluster is a subset of the points, where each point is within X distance of at least one other point in the cluster.

However I cannot think of any way to do this with a PostgreSQL query. Is there any single query to do this, or do I need to write a programme that will externally query the data, and build up the clusters, and continue to query for who knows how many steps?


texts: Looking for good looking maps with nice colors for presentations

texts: I'm currently working on a project where we need some good looking maps with country boundaries displaying the world. We do not need the details of the maps, but the colors are important. It will we used in a Silverlight application that also should be available in offline mode, meaning no server connection. Can anybody help me finding some map resource we could use?


texts:
I'm trying to make a Potential Distribution map for an animal species that, we think, should be dependent on edaphic variables. The problem is the only maps I could find with edaphic information are shapefiles with the info stored in their attribute tables.
Climatic data is easy because worldclim let you download the variables as .bil + .hdr, files MAXENT uses to generate the map.

I have no idea how to convert the shapefiles into, for example, individual variable files that MAXENT can use.

Of course, I'm a complete noob and have zero experience with any GIS. I installed QGIS here and I have access to ArcGIS 10.2.1 in my University.

Any ideas?


texts: AGOL Web AppBuilder - Group Filter Widget with multiple layers

texts:
I am creating a Dashboard web application for displaying utility projects using the AGOL Web AppBuilder interface. I would like to set up a Group Filter widget to filter based on a common attribute value amongst many layers. Is this possible, and if so, how can I accomplish this?

Background:
I have a number of layers in the associated webmap corresponding to their respective feature and utility type.

For example:

Sewer Facilities Project (Point),
Sewer Main Projects (Line)
Water Facilities Project
Water Main Projects
etc.

All of these layers are joined to the same table in my hosted service .mxd. This table is set with a domain value, so all layers are sharing the same field, and the same values for this field:

"Project in design"
"Project in Construction"
"Project Completed"

I would like to have a filter that displays a dropdown to select and only display projects of the selected phase -- when "Project in Design" is selected, only features with this phase will be displayed.

I currently have a filter with a drop-down selections will all layers added to it. I can select any of the phases, but once any filter is applied, all layers disappear.


texts: Buffer Selected Line Segments Using ArcGIS 10 or Python

texts: I haven't  quite found a solution for buffering selected line segments of a road network without extporting the selected line segments to a new file and then buffering.  I want to avoid creating a bunch of new datasets for this task.  Either environment will do (ArcGIS application or Python).  Any suggestions?


texts: Tool to count faulty point-curve connections and number of missing connections

texts: is there any tool in arcgis/arcmap that shows faulty point-curve connection? For example two roads should be at "+" intersection, but they are not? And second- number of missing connections due to undershoots in huge dataset?


texts:
I am trying to create the field which contains its symbol label, is there any tool  to do this?

I use ArcMap10.3


texts: Euclidean distance between a point shapefile and a line or polygon boundary

texts:
I want to calculate the Euclidean distance between (2000) points scattered on an island and the shoreline.
I have the points, a polygon of the island and also a polyline file of the coast in separate shapefiles. All files are georeferenced and all points are marked with x and y coordinates.

Is there a way to easily do this, preferably in a for loop so I get all the distances with a single code.

I already read that the package mapstat is useful, but I am somehow not able to install it.

Any tips will be welcome!


texts: Best site location based on distance and supply

texts:
The question I am trying to answer is 'where best to build a saw mill based on road distances from supply and wood volume from each supply'.

The data I have is wood supply location/volume as point features and a road network.

I have tried the 'mean center'tool with a weight field added which has been great as a starting point but I now need the center point to move along the road network based on road distance and wood supply.

I am using ArcMap 10.3 and I have a Network Analyst licence.


texts:
I have set up the new QGIS Web Client 2 (QWC2) (https://github.com/qgis/qwc2-demo-app) including the demo app and got it running quite easily.

I have created a QGIS Server WMS and also can include it in other QGIS Projects as WMS. So the WMS itself is working. The project's CRS is set to EPSG::25833.

In the qwc2-demo-app-master -> themesConfig.json I have replaced the first theme with my own WMS and specified the EPSG according to my project (EPSG::25833).

I "yarn start" the WebClient and it shows "webpack: Comlied successfully".
The WebClient shows up on localhost:8081 and also my own project shows up, since it is set to default. It is also displayed in themes and legends, the features however are NOT displayed on the map.

I'm stuck at this point and tried to change the QGIS Server CRS to EPSG::4326... but my features never show up on the map.

I cannot find much about this relatively new WebClient development.

I reproduced the exact same procedure on another machine and it works.  However, as soon as I add another Layer to my QGIS Project (the serving WMS layer) is is not showing any feature any more at all... clueless!


texts: Performance is downgrading with multiple feature datasets design for ArcGIS?

texts:
My company develops a spatial planning solution since 2009. Initially it was an ArcMap application to create and maintain spatial plans (competing with AutoCAD folks) but recently we developed "publish them to web" functionality.

Our SDE database design is based on "templating". There is a feature dataset that acts as a template for a new spatial plan. Once a new spatial plan is created we clone that feature dataset into a new feature dataset that will act as a new spatial plan. There are about 10 feature classes (polygons, lines, points, annotations) and 6 relations in the plan feature dataset.

Thus we end up with a separate feature dataset for each spatial plan. This design also helps us to solve juridical issues who can see a spatial plan and when: we set different privileges on plan feature dataset as plan goes through its status cycle (from being designed to official).

During last two years our solution became popular in target market. And here is the problem: customers with 1000 spatial plans say they want to enter 2000 more in one year. This is not a problem for our marketing and sales departments :)

But for development the problem is here: 1000 feature datasets in SDE database (SQL server or Oracle) are hell slow! Say, it takes 5 minutes to open them in ArcCatalog... Upgrading geodatabase to our next version (topologies, representations)? Leave it for the night...

People from my company did not get answer to this question during this year summit.

Do you have experience with such designs? How did you solve the performance issues?


texts: How to get any changes happened in the current SDE version?

texts:
How to get any change done on the map in the current version of sde. using arcObjects
e.g. if i had a line on the map. and in the current sde version i changed the length of that line. how to get this particular line.


texts: Compute a bounding box of a set of coordinates

texts:
Sorry if this sounds naive but I have no background in GIS and have trouble finding relevant information.

I want to zoom/move a map to fit a set of POI. So: Given a set of coordinates (WGS 84 reference frame) I want to find the bounding box (without any libraries).

Is this just plain 2D rectangle math on the latitude/longitude? For example, I found this and several like it. Does this work? How does the curvature of the earth fit into this?

There is also GMSCoordinateBounds (which I'm not allowed to link to for some reason) in the Google Maps iOS SDK that does this. Do they also just do 2D math?


texts:
I have a shapefile with real-time kinematic (RTK) measurements as a line feature, which more or less corresponds with visible features in a raster layer (shaped as lines). These measurements however are not accurate enough, they show about a 1-3 m offset with the features in different raster layers. I'm looking for a way to automate the correction of the offset.

What I've come up with so far is converting the features in the raster image to vector, then measuring the average distance between the two vector files and use this as a measure for offset, then correct this.

The problems I have now are: how to properly measure the offset between two not similar vectors? The RTK is a line feature and the shapefile created from the raster images is one big polygon which looks like two lines, the initial line feature should be in the middle of those in the end. And, after successfully measuring the offset, how do I correct the offset?

Right now I'm using QGIS, but solutions using R or Python are also more then welcome. Let me know if anything in the question is unclear or incomplete!


texts:
I have a raster elevation map that has two tiers of known elevation, and I am trying to produce a raster that will smoothly contour the area between the two tiers with a gradient corresponding to the slopes between the known elevations. Basically, the interior boundary of the raster (and everything inside it) has an elevation of, say, 600m, and the outer boundary is at 624m. The black areas are zero values that will be snipped away and are irrelevant. I would like to generate contours between these two elevation values in a raster format to visualize the slopes between the two elevations. The final raster will be subtracted from a DEM to visualize an excavation accurately.

Please see the attached image of the raster example to see what I am working with. I am flexible with my methodology but I am trying to keep my workflow 100% within QGIS as I am sure there must be a solution here. I have been looking at many of the functions included in GRASS but have not been successful in getting them to generate contours based on this data.


texts: Loading group layer in ArcGIS Map Service?

texts:
My ArcGIS JavaScript application is using an ArcGIS Server Map Service that has a bunch of layers. It also has a group layer. The structure is:

Service
    group_layer_0
        layer_1
        layer_2
    layer_3
    layer_4

After the service being loaded, it not only loads group_layer_0 (which shows layer_1 & layer_2), but also loads layer_1 and layer_2 as individual layers. In other words, layer_1 (or layer_2) is loaded twice. When I turned off group_layer_0, layer_1 is still visible because of the visibility of individual layer_1.

how should I do to avoid this? which means, when I turned off group_layer_0, it would turn off all its sub_layers; but when i turned off layer_1, it would only affect layer_1; besides, layer_1 should be loaded once, instead in current case being loaded as part of group layer and also as an individual layer.

I hope I explained this clearly. and thanks for any help!

Update:

e.g.,:

If want to turn off layer_1, i did:

var visibleLayers = mapServiceLayer.visibleLayers;
var newVisibleLayers = removeLayerIDFromLayerList(visibleLayers, 1);
//remove layer_1's layer ID from current visibleLayers;
//removeLayerIDFromLayerList is a developer-defined function

mapServiceLayer.setVisibleLayers(newVisibleLayers);
//update visible layers based on new list

In the above example, I turned off layer_1, but the new visible layer list still include group_layer_0, since there is no manipulation on this group layer's visibility. and the outcome is layer_1 is still visible. which doesn't make any sense. layer_1 should be off, even though it's part of a group layer.

but this is not the case, at least in my application.


texts: Update driving way speed in open trip planner and how to [easily] obtain the needed ids?

texts:
I'm trying to update the driving way speed limits in OTP. I have seen how opentraffic data graph updater does it, and I have created a custom one that reads from a very simple Json (just the 4 mandatory fields, wayId, startNodeId, endNodeId and average traffic speed ).

I have access to road devices that calculate the average speed in a given point which, in order to give it to OTP, I need to associate with a way, start node and end nodes ids.

To do that I have tried using QGis to load the OSM file data, now I can select the way, and know the wayId. But how can I view its nodes? Is there a simpler way to do this?


texts: Read ESRI shapefile (polygon or polyline) in R Environment

texts: How can I read ESRI polygon/polyline shapefile in R environment using rgdal or other library? can I get a syntax for it?