texts: Attribute Assistant for 10.0 does not have AutoNumber, is there a workaround for ArcMap 10.0

texts:
I work for a water utility division and we are using ArcMap 10.0 with attribute assistant.  In attribute assistant 10.1 and up there is a value method called AutoNumber which will auto increment and populate a field with the next highest integer.  Unfortunately this option does not seem to be available in 10.0 which is my problem.  I have located a script at

http://support.esri.com/en/knowledgebase/techarticles/detail/38517

which will do this but I would like to find a better solution or one that would integrate with attribute assistant 10.0 . Is this possible and has anyone found a solution to this problem?


texts: Landsat: Searching for and downloading cloud-free images for a small area

texts: My study area is quite small (10 x 4 km). I would like to get all cloud free (say < 10% cloud cover) landsat images over that area, for a given time period (say, over 10 years). I cannot use the 10% cloud cover cutoff in EarthExplorer because it is for the whole scene. That is, I might miss a lot of images where it was mostly cloudy over the scene but was cloud free over my study area. So is there a good way to search for and download all suitable landsat images? I have my study area demarcated as a shapefile.


texts: How to add labels to points in elevation plot in R?

texts:
I am a complete noob at using GIS in R and have been trying many packages and functions but get stuck each time. What I want is to generate an elevation map for a location in the Himalayas with additional geographic points that are labeled, and possibly rivers.

What I have so far is the following code that retrieves elevation data and does some basic plotting:

    library(elevatr)
    library(raster)

    #set  projection
    prj4<-"+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs"

    # set plot bounds
    bounds <- data.frame(x=c(85.8,87.4),y=c(26.6,28.2))

    # retrieve elevation raster
    elev<-get_elev_raster(locations = bounds,
                          prj = prj4,
                          z=9,
                          clip="locations")

    plot(elev)

    # Add two dots
    dotcoords<-data.frame(x = c(86.3,87.1), y = c(27.2,27.8))
    rownames(dotcoords)<-c("locA","locB")
    coordinates(dotcoords) <- c("x","y")
    proj4string(dotcoords) <- CRS("+proj=longlat +datum=WGS84")
    pts_dotcoords <- spTransform(dotcoords, CRS(prj4))

    plot(pts_dotcoords, add=T, pch=16)

This generates the following plot:
My question is how I can add labels to the dots showing the associated names ("locA" and "locB")?

This is probably basic stuff but I really need some advice on what approach to take to get me going in the right direction (would it be better to use ggplot2 or ggmap, for example?).


texts: Group by Location within N Meters in PostGis

texts:
In PostGIS, given several Location objects, each of which looks something like this:

{
        :id => 1,
  :latitude => 32.4196,
  :ll_point => "(-10111.1877176122, -5384086.44077608, 3419435.33326163)",
 :longitude => -90.1076,
:updated_at => nil,
   :user_id => 1

}

what would the most optimal way of generating groups of locations, with the end goal of being able to say that User1 shared locations X, Y, and Z with User4. 'Sharing a location' is defined as having been within ten meters of each other.

Is there a way to do something like a group by location within X meters?


texts: Manifold 'Link drawing to query' - can I do this in QGIS?

texts:
Manifold GIS allows the creation of dynamic drawings linked to queries.

Unfortunately I do not have access to Manifold as my employer refuses to pay for it, so am currently investigating alternatives.

Suppose I would like to buffer an alignment linestring stored in a database (SQL Server 2008), sample points at intervals around the buffer and then assign elevations to those points based on a raster:

With Manifold I would simply formulate this as a single nested query
and link a drawing (vector layer) to it to contain the output, then if the
alignment changes I just refresh the query and the linked drawing updates accordingly
With QGIS I can achieve the same result via a multiple step manual procedure (generating superfluous shapefiles along the way) each time the alignment changes, but this is nowhere near adequate for my requirements

Is there a more elegant way of combining functions in QGIS? If not, is this the sort of thing I might be able to achieve server side using PostGIS or similar? (I am completely ignorant of the workings of PostGIS).


texts: Integrating NA (network analysis/routing) service in Web AppBuilder?

texts: We had a requirement for routing in our application. As the arcgis online routing service is credit based we have published our own network routing service in arcgis server 10.6 using the road layer we had. It is working in ArcMap by using only the starting and ending point coordinates (not by start/end point location names).  how do we integrate it in our WAB based application as the "Directions" widget is  not supporting it.


texts:
I have just downloaded QGIS 3.4.15-Madeira on my Mac (Mojave 10.14.6) which local settings include French.
The interface language is however American English.

I went to settings and unticked "override system locale" and restarted QGIS, but it does not solve the issue. English remain the default interface language.


texts: Using OpenLayers 3 Feature onclick event?

texts:
I looked for a way to fire onclick event and I found Adding event handler to feature in OpenLayers 3? at Stack Overflow.  I undestand from this question is that there is no event onclick to features. I want to understand why is it? And maybe there is upcoming to the API and they have added a click events.

the way they offered to do it is like that:

map.on("click", function(e) {
    map.forEachFeatureAtPixel(e.pixel, function (feature, layer) {
        //do something
    }
}

It means i need to fire a click on the whole div map, then look for my specificPixel that the feature contains and then do whatever i want. It sounds primitive and will take too long for a high amount of features.

Is there another way to solve this problem?

Secondly, as i see there is no that much events on OpenLayers 3 api .
Why is it?


texts:
I'm looking for a tool, like Curves and Lines

that would work on polygon layer. At the moment I'm reshaping them feature by feature. Some notes about the data

We use Oracle Enterprise Geodatabase (That's all I know, we outsource servers and DB management)
Data originates from CAD (I was not there for the process, it was done before I came to work here, and it was outsourced)
Each parcel boundary is covered by line features, and each line break is defined by a point (I use the tool mentioned above for fixing the lines)

I have ArcGIS 10.3.1

I've tried Smooth (arcTool) and other tools provided, but since the output is a separate shapefile, it's not useful for me, I need to maintain each parcel attributes and no deleting/adding parcels is preferred because other uses are constantly working on the areas I'm fixing (other operations).

The image shows how the parcel's curve is continuous vertices rather than being a true  curve.

It doesn't have to be a tool, If there's an more secure and reliable way to fix this problem.


texts:
I'm trying to use the pyrealsense library from Intel in QGIS in Win10

The current version of mine is QGIS 3.4, which comes with Python 3.7

But the pyrealsense supports pip install on 3.6 and 2.7

I tried installing QGIS2 but can't get easy_install and then get pip.

And also QGIS3 is more preferred

I tried a lot of other methods but still currently not one working found

I would like to know if there is either:

Downgrade the Python version of QGIS to python3.6

Or let QGIS console open the system Python


texts: Do "free" ArcGIS Online service credits expire or accumulate?

texts:
Do the "free" ArcGIS Online service credits expire or accumulate? For example, if my license comes with 100 credits, and I have used 50 at time of renewal, will after renewal I have 150 or 100?

I know that additionally-purchased credits are good for 24 months. Unclear as to the life of the "free" credits.


texts:
There seem to now be many methods for working with an Excel spreadsheet in QGIS. Some allow read-only, and some require one Excel format or another. Some also have historically been unreliable (e.g. where trivial changes to the Excel sheet break the QGIS project).

What are the main different methods for working DIRECTLY with Excel at QGIS version 2.18?

Clearly there are more convoluted methods - Excel data can be exported and then imported to many other formats (e.g. csv), but I'm asking about more direct methods of making a direct link with the original Excel file (so changes are seen in QGIS). Often other people wish to maintain live data in Excel - and GIS professionals need (if being efficient) to maintain live links to that live data.

I (think I) know of the following (but have questions about each). I suspect that the best answer will simply list the basic facts about each of the various methods - these and others.

Save using older .xls format. Bring into QGIS using a virtual xls table in a Spatialite database. This I think results in a read only table but it is possible to define that the first row has column headings.
Save using .xlsx format. Drag into QGIS or open using 'add vector layer'. No easy way to define if top row is column labels.
Use the 'Spreadsheet Layers' plugin. It is possible to define a header row, and x/y data. The plugin works by creating a VRT file which sits in parallel to the original file. Presumably this VRT file could be moved (if the file reference in this VRT file continues to point correctly at the original file). The VRT file records details about the spreadsheet (data types etc).

In all cases the header row data (i.e. field names) in Excel needs to be handled carefully - odd characters (etc) can break things.

I realise that this question doesn't ask for a simple 'how to' answer, but working with Excel spreadsheets is trivial in some software so I believe that some clear comments about the facts about the various methods are justified. The information above is incomplete.


texts: GDAL C++ API: How to create PNG or JPEG from scratch

texts:
I'm new to GIS and GDAL. My question probably is very basic, but I couldn't find answer. May be I don't understand GDAL ideology.

I need to create raster images from scratch, for example, JPEG or PNG. Their drivers don't support Create function - only CreateCopy. What is the common technique of new files creation in this case?

In principle, I can try to create Tiff because its driver suports Create(). Next, I can use CreateCopy() for PNG or JPEG using this Tiff. But such method looks indirect and unnatural for me. Also I suppose that this procedure can be too memory hungry if rasters are large.

I dealt with some image libraries before, they usually provide direct and simple way of bitmaps creation. Could somebody show me right direction for GDAL?


texts: QGIS 2.4.0 marker line interval does not work properly when line is bending

texts:
I am trying to style a line using marker line option. When I am using styling the line in the style window everything looks nice as you can see on the print screen bellow.

I am using three line features (units are mm):

simple line: dash pattern 7 2
marker line: interval 9
marker line: interval 9, offset 5

However the result in QGIS looks horrible as shown on the next picture. Could somebody tell me how to avoid this? It looks like QGIS is not properly listening to the dash pattern and interval settings.


texts: Record IDs of features loaded on request in GeoServer

texts:
I'm switching over to GeoServer, and we need to be able to record the number of features/polygons each user loads for licensing reasons.

I've got the Monitoring extension installed and the Auditing set up, but I can only see the ResponseLength.

Is there any way to be able to record the certain properties of a result from each request for a download of WFS data, whether by default or adding some kind of post-processing functionality? It'd need to be this way as opposed to just a count of the number of features to ensure that we're not counting each feature twice.

I'm fetching the data from a PostGIS database and returning to OpenLayers v4


texts: Why shapefile loaded in canvas always appear in random color?

texts:
I created the simple frame .UI file, compiled with pyuic4, and tried the ShapeViewer.py to load the shapefile "world_borders.shp". However, the map is always loaded in different random colors every time I rerun the python program.
Why is this so? I tried loading the map in QGIS program & it loaded consistently with the same colors.

#! /usr/bin/env python

from PyQt4.QtCore import *
from PyQt4.QtGui import *
from qgis.core import *
from qgis.gui import *
import sys, os

# Import our GUI.
from shapeviewer_gui import Ui_MainWindow

class ShapeViewer(QMainWindow, Ui_MainWindow):
  def __init__(self):
    QMainWindow.__init__(self)

    # Required by Qt4 to initialize the UI.
    self.setupUi(self)
    self.layers = []

    # Set the title for the app.
    self.setWindowTitle("ShapeViewer")

    # Create the map canvas.
    self.canvas = QgsMapCanvas()
    self.canvas.useImageToRender(False)

    # Set the background color to white.
    self.canvas.setCanvasColor(Qt.white)
    self.canvas.enableAntiAliasing(True)
    self.canvas.show()

    # Lay our widgets out in the main window using a vertical box layout.
    self.layout = QVBoxLayout(self.frame)
    self.layout.addWidget(self.canvas)

    actionZoomIn = QAction(QString("Zoom in"), self)
    actionZoomOut = QAction(QString("Zoom out"), self)
    actionPan = QAction(QString("Pan"), self)

    actionZoomIn.setCheckable(True)
    actionZoomOut.setCheckable(True)
    actionPan.setCheckable(True)

    self.toolbar = self.addToolBar("Canvas actions")
    self.toolbar.addAction(actionZoomIn)
    self.toolbar.addAction(actionZoomOut)
    self.toolbar.addAction(actionPan)

    # Create the map tools.
    self.toolPan = QgsMapToolPan(self.canvas)
    self.toolPan.setAction(actionPan)
    self.toolZoomIn = QgsMapToolZoom(self.canvas, False) # false = in
    self.toolZoomIn.setAction(actionZoomIn)
    self.toolZoomOut = QgsMapToolZoom(self.canvas, True) # true = out
    self.toolZoomOut.setAction(actionZoomOut)

    self.pan()

    # Layout is set - open a layer.
    # Add an OGR layer to the map.
    self.Add_Map_Layer("qgis_data/world_borders.shp")

    self.connect(actionZoomIn, SIGNAL("triggered()"), self.zoomIn)
    self.connect(actionZoomOut, SIGNAL("triggered()"), self.zoomOut)
    self.connect(actionPan, SIGNAL("triggered()"), self.pan)

  def zoomIn(self):
    self.canvas.setMapTool(self.toolZoomIn)

  def zoomOut(self):
    self.canvas.setMapTool(self.toolZoomOut)

  def pan(self):
    self.canvas.setMapTool(self.toolPan)

  def Add_Map_Layer(self, mapFileName):
    # Add the layer
    if mapFileName == "":
      layer = QgsVectorLayer()
    else:
      layer = QgsVectorLayer(os.path.abspath(mapFileName), mapFileName[mapFileName.rfind("/")+1:-4], "ogr")

    if not layer.isValid():
      print "Layer failed to load! :("
      return

    # Add layer to the registry
    QgsMapLayerRegistry.instance().addMapLayer(layer);

    # Set extent to the extent of our layer
    self.canvas.setExtent(layer.extent())

    # Set up the map canvas layer set
    self.layers.append( QgsMapCanvasLayer(layer) )
    self.canvas.setLayerSet(self.layers)

def main(argv):
  # create Qt application
  app = QApplication(argv)

  # Environment variable QGISHOME must be set to the install directory
  # before running the application
  qgis_prefix = os.getenv("QGISHOME")
  if qgis_prefix == None:
    qgis_prefix = "/usr"

  # Initialize qgis libraries
  QgsApplication.setPrefixPath(qgis_prefix, True)
  QgsApplication.initQgis()

  # create main window
  wnd = ShapeViewer()
  # Move the app window to upper left
  wnd.move(100,100)
  wnd.show()

  # run!
  retval = app.exec_()

  # exit
  QgsApplication.exitQgis()
  sys.exit(retval)

if __name__ == "__main__":
  main(sys.argv)


texts: Is there any map SDK or tool that can show the navigation that passes through the north pole in a flat way? Basically If I use a Mercator projection map it is only flat and cannot show the north part of as how it is seen in a globe. I need something that will show a plane passing through the north pole in a 2D map.


texts: Getting data manually from CARTO.js four dataview?

texts:
I have a series of CARTO.js four dataviews, but I want to manually get the data instead of listening to the 'dataChanged' event (when multiple datasets are created, this end in multiple errors). I tried using the getData

this.widgetFormula = new carto.dataview.Category(this.dataset, 'class',
                    {operation: carto.operation.SUM, operationColumn: 'pavementlengthinfeet'});

client.addDataview(this.widgetFormula);

var xx = this.widgetFormula.getData();

but this returns an empty carto.dataview.categoryData object (a JSON object with the regular keys, but values empty). The online documentation doesn't provide an example to use the getData. Any idea?


texts: Building an offline tile server for ESRI World Imagery

texts:
I have a JavaScript and Leaflet project which serves users OSM and ESRI maps offline. I have 2 map layers. I successfully built the OSM part with switch2osm tutorial. But I couldn't find the similar solution with ESRI World Imagery. I want to convert below code
var esri = L.tileLayer('http://server.arcgisonline.com/ArcGIS/rest/services/World_Imagery/MapServer/tile/{z}/{y}/{x}');

to:
var esri = L.tileLayer('http://127.0.0.1/MapServer/tile/{z}/{y}/{x}');


texts: Using ArcMap Topology to make admin boundaries match up?

texts:
I have two separate datasets ('Wards' and 'Census Output Areas'), whose boundaries don't quite match up (see below - wards in red, COAs in grey).

Is there a way of automatically adjusting the Ward boundaries so that they align with the Census Output Areas?  I know that I could use the 'Align Edge Tool' on the Topology toolbar to fix individual errors, but it would take a long time to do them all. After reading about topology rules (here), I think the most appropriate rule is 'area boundary must be covered by boundary of', but I'm not quite sure how I would use it and it seems there are no fixes for this anyway..  Any suggestions?


texts: Rule-based styling road crossing in QGIS

texts:
Styling/symbolizing (rule-based option) a roads vector layer with intersections of streets with different styles I have incorrect (and annoying) overlaps of the different symbolization styles. is it possible fix the problem?

I also add the symbolization level window


texts: Measuring percentage overlap between home ranges and another shape file in R

texts:
I have tracking data from multiple animals and I have created kernel density estimates of their home ranges. I want to find out the proportion of each home range area that falls within protected areas. I'm using the adehabitatHR package to calculate the home ranges and the wdpar package to download the protected areas from the World Database on Protected Areas. But I'm not sure how to link the two to calculate the proportion of overlap. Here's an example of the code I have so far:

library(amt)
library(wdpar)
library(rgdal)
library(adehabitatHR)

data(deer)
mini_deer <- deer[1:100, ]
mini_deer$id <- ifelse(mini_deer$x_ < 4313000, "A", "B")

mini_deer <-
  mk_track(
    mini_deer,
    .x = x_,
    .y = y_,
    .t = t_,
    id = id,
    crs = CRS(
      "+init=epsg:3035 +proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs"
    )
)

# Keep only the id column and the x and y coordinates
mydata.sp <- mini_deer[, c("id", "x_", "y_")]
coordinates(mydata.sp) <- c("x_", "y_")

# kde
kernel.ref <- kernelUD(mydata.sp, h = "href")  # href = the reference bandwidth
mydata.kernel.poly <- getverticeshr(kernel.ref, percent = 95, unout = "m2")
print(mydata.kernel.poly)  # returns the area of each polygon

# get the protected area data
ger_raw_data <- wdpa_fetch("Germany", wait = TRUE)
pa_data <- st_transform(ger_raw_data, 3035)

# scene missing
# proportion of overlap for A and B


texts:
I have a mini project to do which is to assign coordinates to each tiff image and later covert these images to ECW format.  Im using Qgis 216 Nobedo. I have a txt file with the names of pictures and coordinates +pictures themselves.
Is there a way I could assign coordinates to each picture on QGIS?


texts:
My aim is to multiply all x coordinates by 2 and all y coordinates by 3 of an R sf object. How can I achieve this?

When running the following code the polygons are duplicated:

library(sf)

# Create Polygon: https://cran.r-project.org/web/packages/sf/vignettes/sf1.html
p1 <- rbind(c(1000,1000), c(1000,2000), c(2000,2000), c(1000,1000))
p2 <- rbind(c(500,500), c(500,100), c(1000,1000), c(500,500))
pol <-st_polygon(list(p1, p2))
pol_geom = st_sfc(pol, crs = 4326)

# Trying to multiply x coordinate by 2 and y coordinates by 3
pol_geom_new <- st_geometry(pol_geom) * c(2,3)
st_bbox(pol_geom)
st_bbox(pol_geom_new)
pol_geom; pol_geom_new


texts: Identifying which lines cross building polygons using ArcMap?

texts: I have the following lines generated from coordinates with XY to Line. I would like to select and delete the lines which cross and are contained within the buildings. I need to do this to find the distance between buildings. I have tried using the select by attribute tool with no luck so far. Also due to the high number of lines Arcmap 10.3 tends to crash a lot.


texts:
My question is a follow up to Convert a raster to GeoJSON polygons in PostGIS.

I'm trying to accomplish querying point data from the database, and then turn that into a raster for calculating (average of a field), and finally get a GeoJSON output.

My query is as follows:

SELECT row_to_json(fc)
     FROM
     ( SELECT 'FeatureCollection' As type, array_to_json(array_agg(f)) As features
         FROM (SELECT 'Feature' As type,
                ST_AsGeoJSON((px).geom)::json As geometry
                , row_to_json((SELECT props FROM (SELECT (px).id as Id, (px).particles as Particles) as props )) As properties
               FROM
                (SELECT ST_PixelAsPolygons(
                    ST_SetValue(
                        ST_SetValue(
                            ST_AddBand(
                                ST_AsRaster(sd.point,5, 5),
                            '8BUI'::text, 1, 0),
                        2, 2, 10),
                    1, 1, NULL))
                FROM sensordata as sd) px
                 ) f
    ) fc;

This results in an error:

pg_query(): Query failed: ERROR: column px.geom does not exist LINE 5: ST_AsGeoJSON((px).geom)::json As geometry

I think it's something obvious, but I just can't seem to figure out what's up.

After attempting the solution proposed by Tilt I get a new kind of error message:

ERROR: function st_asgeojson(record) does not exist LINE 5: ST_AsGeoJSON(px.geom)::json As geometry ^ HINT: No function matches the given name and argument types. You might need to add explicit type casts.

Here's the final working query.

WITH pixelrecords AS (
        SELECT ST_PixelAsPolygons(
            ST_SetValue(
            ST_SetValue(
                ST_AddBand(
                ST_AsRaster(data.point,5, 5),
                '8BUI'::text, 1, 0),
            2, 2, 10),
            1, 1, NULL)
    ) px, data.id as Id, data.particles as Particles
            FROM sensordata as data
        )
        ,features AS (
        SELECT 'Feature' As type,
            ST_AsGeoJSON((px).geom)::json As geometry
            --This will likely go wrong since there is no id or Particles in the px record, only val, geom, x and y
            --You can better get the values from the original sensordata records and transfer them along with the polygon
            ,row_to_json((SELECT props FROM (SELECT (Id), (Particles)) as props )) As properties
               FROM
            pixelrecords px
        )

        SELECT
            'FeatureCollection' As type,
            array_to_json(array_agg(f)) As features
        FROM features f;


texts: Downloading Large AGOL hosted data with attachments

texts:
We have a large data collection project with attachments, and have gotten to the point where the data can no longer be downloaded as a file geodatabase without timing out (it can be downloaded as a shapefile).

Is there a method other than the below or downloading it in separate chunks to get a large dataset from AGOL as a file geodatabase?

Trying these workarounds:

Unable to export a large hosted feature service with attachments to a file geodatabase from a REST endpoint

Unable to export a feature layer with attachments from ArcGIS Online for Organizations to a file geodatabase

Download attachments from a hosted feature service

Results in the following:

arcgisscripting.ExecuteError: ERROR 999999: Error executing function.The object has been deleted and is no longer valid.ERROR: code:504, Your request has timed out., Proxy server got bad address from remote server (verify the server is running).Failed to execute (FeatureClassToFeatureClass).

OR


texts: I'm heading out to the Arctic Circle (Spitsbergen) I have a Sokkia DGPS and Magnet on a tablet device (connects to the DGPS and logs the data), they connect via blue tooth etc. Can I bring the data straight into QGIS bypassing Magnet?


texts: Draw heatmap only over water layer using Leaflet?

texts:
I need to draw a heatmap, but only over water because I want to represent water heat level. I am using Leaflet. How can I do it?

I found an only-water layer, maybe it would be usefult for developing a solution: http://openlayers.org/en/v3.2.1/examples/tile-vector.html?q=tile

Edit:

An option would be:
 1. draw water + 2. draw noise +  3. draw land.
However, I can not find a way to draw only the land-side of a map.

Another option would be, before painting, discarding the points that are not water. Would it be possible in the client side?


texts: Extracting LAS/LAZ info in a performant way with PDAL?

texts:
I can do this:
import json
import subprocess as sp

laz_file = somefile.laz
r = (sp.run(['pdal', 'info', laz_file], stderr=sp.PIPE, stdout=sp.PIPE))
json_info = json.loads(r.stdout.decode())

and parse through the JSON, get the info I need (EPSG code and min/max x/y).  It does "work."
However calling this workflow on a laz file, especially when some of mine are over 1 gb, takes significantly longer than a las file.  As I am iterating over hundreds of laz files (which change daily and sometimes hourly) and intending the output to be used as a dynamic map tracker for our holdings.
How can I extract the EPSG code and the bounding box in a more performant way?


texts: GeoServer generates no tiles in Format EPSG:4326, preview is completly blank

texts:
I am currently trying to setup a geoserver. I need to request it in MATLAB. I tried this TMS request
http://localhost:8080/geoserver/gwc/service/tms/1.0.0/layer/0/0/0.jpeg

I always get this error message:
Coverage [minx,miny,maxx,maxy] is [1, 0, 0, 0, 0], index [x,y,z] is [0, 0, 0]

I looked in the Cached Layer Preview. When I used the native projection it worked and showed me the tiles in the browser, using EPSG4326 nothing shows up only the interface and a blank map. I tried to switch from JAIEXT to JAI which didn't helped. Furthermore, the native projection is a wrong one, though I do not know how to change it.
I am using an ImageMosaic made out of many air images (GeoTIFF).The original projection is : EPSG:25832 while the server claims it's native projection is EPSG:404000. I tried to specified EPSG:25832 as the correct projection but it didn't helped either. I tried entering the boundary boxes manual as the calculated values were a little bit off.


texts: Combining two euclidean distance rasters in Raster Calculator

texts: I have two euclidean disance rasters and would like to combine them for use in a multi criteria evaluation. Could anyone advise me on the syntax to use within Raster Calculator (ArcGIS 10.3.1)? There are not any areas of NoData.


texts: How to get the geometry type of a WFS vector layer in OpenLayers 3?

texts:
I've been trying to get the geometry type of a vector layer OpenLayers 3 with no success. I know that the geometry type can be obtained from a feature, but in my case I have some layers without features (empty layers to start drawing, I need to choose the propper draw interaction).

I've found nothing in the docs, and looking at the ol.layer.Vector and the ol.source.Vector objects in the Chrome console, seems to be some obscure methods and properties like 'A', 'B', 'Fa', 'U', etc. I've been exploring them but I couldn't find wich one has the geometry type, if there is.

I must be missing something... Any clue?

EDIT

And I've realised I'm missing an important thing. The source of the layer is always an WFS layer. I was thinking that the vector layer object was bound to a geometry type (like in most GIS software), and is not, it can contain mixed geometry types. And source object maybe can contain anything too, so that's why there is nothing like a 'getGeometryType' method in ol3? Should I query directly the WFS server?


texts: Accessing the Google Earth Pro image dataset for exporting a GeoTIFF

texts:
I would like to extract with Earth Engine, Google Earth Pro GeoTIFF images in UTM projection using opposing corner coordinates.  In Google Earth Pro, the ~6” imagery I assume to be already processed raw satellite imagery that it is now in an ImageCollection. Can this ImageCollection be accessed directly where opposing corners are provided in the earth engine coder and then exported as a GeoTIFF? What would be the name of the ImageCollection?  Is there example code for doing this?
Is there example code for doing this?


texts: Fill in spaces between polygon features, voronoi style

texts:
I'm looking for a way to merge polygons across the empty space between them with a Voronoi-type result. There is a very similar question here, but the approaches outlined in the comments did not produce the result I wanted (it created too many extra polygons, and all were very jagged). So, I did it manually as you can see below. I'd like to know a way to do this non-manually as I have even more complicated polygons that need this treatment.

What I have (before merge):

What I want (after merge):

Edit: There's a buffering question here that may be useful to me so long as I clip the resulting buffer's world extent, but so far it has no answer. Also in my case whenever I try buffering this shapefile, the rectangular ocean boundaries wind up getting buffered as well.


texts:
Is there a way to select a point from my map's current extent? I want to use ArcObjects to tell me if there is at least one point visible in my extent. I'd like to do this entirely in VBA, and without selecting any features or using any mouse-click events. (Selection can be done programmatically, just not manually.)

An demonstration would show that I have a point-type shapefile active and positioned at the top of my Table of Contents. I run my code which then displays a message box showing the longitude and latitude of the first point that is found within the current extent. Any point can be found and it's coordinates returned, just as long as that point is within the given extent.


texts: Is it possible to use different symbols instead of colours in one shapefile in QGIS?

texts:
I want to show 50 bird species on a map and give each species a different symbol (with a different colour) so they can be visually distinguished in the map.

Is this possible?


texts: Clipping image with Basic (ArcView) level license of ArcGIS Desktop?

texts:
Is there a way to trim a geo-referenced air photo to an irregularly shaped boundary?

I tried the clip tool shown in my attachment, but it didn't work.  It looks as though the output extent field requires a rectangularly shaped polygon as the clipping feature.

Can I even do this in ArcView?

I've never had a request to do anything like this before.

I do use Photoshop and other similar software packages.  I could just go in to an environment like this & start painting everything outside of my boundary in white, but then I would have to come back in and geo-reference the image again, and it just strikes me that there must be a faster and easier way (even though I only have ArcView 9.3).


texts: QGIS write in attribut table while mouse clicking

texts:
I want to draw a track in QGIS with mouse-click to use it inside the print composer(atlas generation).

The best way to do this for me is too click on the screen and the python script in the background set an autoincrement on my first column in the attribute table of my layer.

I already use it with the "add feature" function in QGIS, but it's not really comfortable.

For each click I make, a new row should be created. Does someone have an idea how to program this task with PyQgis?


texts:
Batch Renaming Relationship Classes

I have consulted the post above, but found that the solution is out of date and no longer works. It suggests two options to initially get the list of relationship classes:

Option 1

import arcpy
import os

desc = arcpy.Describe(r'C:\Users\<username>\Documents\ArcGIS\Default.gdb')
print desc
for child in desc.children:
    if child.datatype == "RelationshipClass":
        print child

Option 2

import arcpy, os
workspace = r'C:\Users\OWNER\Documents\ArcGIS\Default.gdb'

# Generate the list of relationship classes
rc_list = [c.name for c in arcpy.Describe(workspace).children if c.datatype == "RelationshipClass"]
print rc_list

Both options returned "[]" which is an empty array. This means these are no longer viable options since I have literally 150 relationship classes to rename and none were returned.

Eventually I want to rename my relationship classes as such:

import arcpy, os
from arcpy import env

# Set workspace
env.workspace = "C:/Users/XXX/Desktop/test.gdb"

#get list of relationship classes
rc_list = [c.name for c in arcpy.Describe(env.workspace).children if c.datatype == "RelationshipClass"]

#remove all underscores
for rc in rc_list:
    out_data = rc.replace("_","")
    arcpy.Rename_management(os.path.join(env.workspace,rc), out_data)

#get the list of relationship classes again
rc_list = [c.name for c in arcpy.Describe(env.workspace).children if c.datatype == "RelationshipClass"]
for rc in rc_list:
    #
    if "reading" in rc OR "action" in rc:
        #add underscores between name and ATTACHREL
        if "ATTACH" in rc:
            out_data = rc.replace("ATTACHREL","TblGPS__ATTACHREL")
            arcpy.Rename_management(os.path.join(env.workspace,rc), out_data)
            print (rc)

        #add underscores between name and rel
        if "rel" in rc:
            out_data = rc.replace("rel","TblGPS__rel")
            arcpy.Rename_management(os.path.join(env.workspace,rc), out_data)
            print (rc)
    else:
        #add underscores between name and ATTACHREL
        if "ATTACH" in rc:
            out_data = rc.replace("ATTACHREL","PtGPS__ATTACHREL")
            arcpy.Rename_management(os.path.join(env.workspace,rc), out_data)
            print (rc)

        #add underscores between name and rel
        if "rel" in rc:
            out_data = rc.replace("rel","PtGPS__rel")
            arcpy.Rename_management(os.path.join(env.workspace,rc), out_data)
            print (rc)

How can I extract the relationship classes from the geodatabase so I can rename them in bulk? I have them in a file geodatabase right now but eventually I will have them in SDE.


texts:
I'm pretty new to GIS and I've learned about both the geographic coordinate system, which defines a location on a 3d model of the world and the projected coordinate system to map this to a 2d surface and I thought I had a decent grasp of it.
Other questions like Is WGS84 a Coordinate system or projection system? explain that EPSG:4326 / WGS84 is the first (coordinate system) and not the latter (projection).
And that's why I'm confused about this page in the openlayers example: https://openlayers.org/en/latest/examples/reprojection-wgs84.html
Here they explicitily talk about the WGS84 as a projection:
This example shows client-side reprojection of OpenStreetMap in WGS84.
...
new View({
    projection: 'EPSG:4326',
    center: [0, 0],
    zoom: 2,
  }))

So I'm a bit confused by this, can someone explain what is meant here?


texts: How to set schema parameter in osm2pgsql calling

texts:
How I can to set schema name parameter when call osm2pgsql?
For example:

osm2pgsql --create -s -C 2500 -G  -S D:\default.style -d osm_db -o input -U osm_user -W -H localhost -r pbf D:\osm.pbf

I set the -o parameter but information was loaded in public...
How do I set a target schema?


texts:
Is it possible to change visibility or opacity of google base map ?
I am adding it like in simple example, and adding data layer as overlay and want to see only vector markers. Is it possible?

var map;
function initialize() {
  var mapOptions = {
    zoom: 8,
    center: new google.maps.LatLng(-34.397, 150.644)
  };
  map = new google.maps.Map(document.getElementById('map-canvas'),
      mapOptions);
}

google.maps.event.addDomListener(window, 'load', initialize);


texts: Error when zooming to layer in PyQGIS application. On-the-fly transformation the issue?

texts:
I am trying to set the extent of my window to a layer on my canvas. This has worked smoothly until I started using Web Mercator maps instead of WGS84 maps. My features are getting their geolocation from an external source which feeds in decimal degrees, so I am using on-the-fly transformation to get it projected properly on Web Mercator. The problem is that in my code, when I set the extent to my layer, it appears to try to do it in decimal degrees, putting my window in the middle of the Atlantic Ocean.

My extent code is very simple:

self.mapCanvas.setExtent(self.myLayer.extent())

My canvas CRS code is straightforward as well:

self.mapCanvas = QgsMapCanvas()
self.mapCanvas.useImageToRender(False)
self.mapCanvas.setCanvasColor(Qt.white)
self.mapCanvas.setCrsTransformEnabled(True)
self.mapCanvas.setDestinationCrs(QgsCoordinateReferenceSystem(3857, QgsCoordinateReferenceSystem.EpsgCrsId))
self.mapCanvas.show()

Is anything not right?


texts:
In QGIS 2.18.7 I want to change the projection of my shapefiles into a form that preserves distances for the sake of analysis. But when I do this from WGS84 to "WGS84 UTM zone 48N" or "south asia equidistant conic" I get very small measurements of distance for known features. For example 1km real world features are about 6cm using the measuring tool.

As suggested in other posts I have tried changing the project projection and resaving shapefiles in new projections but I still keep getting the same tiny feature representations. I have also tried turning 'on the fly' CRS transformations off

How can I preserve the real world distances on my map?


texts:
We are generating cache from external WMS. In some zoom levels (17, 18 and up) all map is like chessboard, some tiles are ok but some tiles are just blank. When zooming in or out, different tiles are blank. We have tried to reseed, add the layer again, also add external WMS to Store again. We are out of ideas.

Map looks like this. Ok area surrounded by blank area.

No errors in logs. Tile size is 256x256. We generate tiles for EPSG:900913 jpeg. External WMS is ArcGIS. When we explore the map with QGIS, then external WMS seems to be ok.

Has anyone had the same issue?


texts: Keeping polygon's field when generating random points within overlapping polygons

texts:
I have a shapefile with a lot of polygons, representing events occuring for a given year. Every area can be covered several times, for every year it occured on the same coverage. Each polygon is characterized by a certain quantity (ie. number of hosues builded).

I generate random points within thoses polygons to show the amount of this quantity (with the toolbox). By doing this, I'm loosing all the fields of the polygon, and am not able to categorize the points according to the year they're referring to. It's not possible to perform a spatial join because of the overlapping polygons.

The answer to this question (How to create random points in overlapping polygons?) won't work because I'm not using fixed points generation.

So how can I generate random points within polygons and keep all the fields of the polygons the points are referring to?


texts: Calculate area with OpenLayers4 and Turf.js

texts:
I've been working with OL4 for a few months now but this is the first time I try to calculate an area and after some research looks like Turf.js seems like the way to go.

First, allow me to add all the code which, IMO, is relevant here, projection-wise:

This is where I add all the features to the map:

proj4.defs("EPSG:3795", "+proj=lcc +lat_1=23 +lat_2=21.7 +lat_0=22.35 +lon_0=-81 +x_0=500000 +y_0=280296.016 +datum=NAD27 +units=m +no_defs");

var parser = new ol.format.WMSGetFeatureInfo();
var ajaxCall = $.ajax({
    url: url,
    success: function(response) {
        results = parser.readFeatures(response);
        $.each(results, function(i, v) {
            var feature = v;
            feature.setGeometry(v.getGeometry().clone().transform('EPSG:3795', 'EPSG:3857'));
            vector.getSource().addFeature(feature);
        });
    }
});

This is my view definition:

var view = new ol.View({
        center: ol.extent.getCenter(
            ol.proj.transformExtent(extent, 'EPSG:4326', 'EPSG:3857')
        ),
        zoom: 11.7
});

I have two layers on the map and I am trying to calculate the area of the intersection between two features, each belonging to one of the layers. This is how I am calculating the area:

        var g1 = feature.getGeometry(); //<-- 'feature' comes from openlayers' event 'forEachFeatureAtPixel'

        var gc1 = g1.getCoordinates();
        var gt1 = g1.getType();
        var poly1 = getTurfHelper(gt1, gc1);
        $.each(__map.getLayerByName('layer2').getSource().getFeatures(), function(i, f) {
            var gt2 = f.getGeometry().getType();
            var gc2 = f.getGeometry().getCoordinates();
            var poly2 = getTurfHelper(gt2, gc2);
            var intersection = turf.intersect(poly1, poly2);
            if (intersection !== undefined && intersection !== null) {
                console.log(turf.area(intersection));
            }
        });

function getTurfHelper(type, coordinates) {
    if (type === 'Polygon') {
        return turf.polygon(coordinates);
    } else if (type === 'MultiLineString') {
        return turf.multiLineString(coordinates);
    } else if (type === 'MultiPolygon') {
        return turf.multiPolygon(coordinates);
    } else if (type === 'Point') {
        return turf.point(coordinates);
    } else if (type === 'LineString') {
        return turf.lineString(coordinates);
    } else if (type === 'Feature') {
        var geometry = {
            "type": "Feature",
            "coordinates": coordinates
        };
        return turf.feature(geometry);
    }
    return undefined;
}

Now, the problem at hand:

20450142473549.008 is what I am getting as an area. According to turf's docs, turf.area's result is in square meters and in my test data the intersection should be somewhere around 8000m^2 which is why I think I'm doing something really wrong. Maybe it's a projection issue or maybe just not using turf properly. Any help appreciated.


texts: How to export settings for a network dataset

texts: Is there a way to export/save the settings for a network dataset? Setting up the connectivity groups, attributes, evaluators etc can take quite a while, even for a small dataset. If the network dataset needs to be recreated multiple times for whatever reason, it quickly becomes a painful task.


texts:
Does satellite imagery have a unique azimuth for each of its spectral bands?

Is it right to say that if a sensor has 8 spectral bands, then it may have 8 different azimuth angle?


texts: Does Earth Engine store data that is in Google Maps?

texts:
I noticed that Google Maps has a bicycle layer that can be turned on.  Extracting this polyline layer and bringing it into ArcGIS Online would be useful.  However, I don't know how to get that out of Google Maps.
Is Earth Engine associated with Google Maps and, if so, can someone point me in the right direction to finding that data?


texts: Simplest way to export a data set out of OpenStreetMap's Overpass API

texts:
I would like to create a web map showing just one type of point of interest. They are made available via the Overpass API and count up to only few thousands for the whole world.

I don't expect them to change at all, they will at most slowly increase in number (they represent some very expensive equipment).

According to all the above, I thought to export all of them manually, store them as a JSON file that gets shipped to the clients when they connect for the first time. In some months I will redo the procedure and update the data.

I was able to implement this scenario using a very small subset of data. Now I would like to export the whole data set. How should I do?


texts:
I have come across a set of poorly generated rasters which I need to align and regularize. I am open to any quick and dirty way to do this, but prefer QGIS since I have been using it for all my operations thus far.

The rasters contain positive ranges of floating point values, have a NODATA value defined, but also have many cells containing nan values. I would like to convert all nan values to the currently defined NODATA value. To achieve this, I tried many combinations of using  the following:

gdal raster calculator (using numpy where and isnan functions)
raster calculator
GRASS r.null tool

However I can't seem to get a workflow going. The gdal raster calc with numpy functions was close, but unfortunately theres no way of referring to the currently defined NODATA value by variable (as far as I know).

I would like to avoid using Python tools like rasterio since rasters should not be delivered in this format anyway (no use scripting this).

Does anyone have a way of achieving this?


texts: How to calculate the area of classified pixels in Google Earth Engine?

texts:
I'm using the Python API of GEE in a Jupyter notebook.
I'd like to perform a surface evaluation of classified pixels.
I have a ee_map with one band called gfc that have values in [1,2 ... 19, 30, 40, 50, 51] characterising the gain or loss of treecover.
I try to evaluate the surface of each class in hectares:
#define the pixel resolution
    res = 30

    hist = ee_map.reduceRegion(**{
      'reducer': ee.Reducer.autoHistogram(),
      'geometry': ee.FeatureCollection(assetId).geometry(),
      'scale': res,
      'maxPixels': 1e12
    })

    hist = pd.DataFrame(hist.getInfo()['gfc'])

    #add column name
    hist.columns= ['code', 'pixels']

    #won't work
    #hist['area'] = hist['pixels']*res*res/10000

Problem : due to projection issues I don't know the surface of each pixel.
The documentation suggest to use ee_map.multiply(ee.Image.pixelArea()) but it will only work if the map have binaries values.
Is there a way to combine it with ee.Reducer.autoHistogram() or am I force each of them in a loop?


texts: Verifying an answer in a GIS book (D8 flow direction)

texts:
I am currently reading this book and I found the answer a bit confusing. Would be nice if someone could elaborate on the accuracy of the following sentence:

The values 64 and 16 have the highest frequency, which means the direction of the surface flow is north and west.

Also, if there is a better reading in regards to this topic.

The book: https://books.google.com/books?id=FURADwAAQBAJ&pg=PA240&dq=The+values+64+and+16+have+the+highest+frequency,+which+means+the+direction+of+the+surface+flow+is+north+and+west.+Step&hl=en&sa=X&ved=0ahUKEwibpMCSo4_aAhWBdt8KHWvrDrQQ6AEIJzAA#v=onepage&q=The%20values%2064%20and%2016%20have%20the%20highest%20frequency%2C%20which%20means%20the%20direction%20of%20the%20surface%20flow%20is%20north%20and%20west.%20Step&f=false


texts:
I am using Geoserver 2.4.0 to view the UK map in my web map application by using the openlayers and wms service.

Now I am looking into view the terrain data in 3d with WebGL enabled browser and got some examples in Geoserver W3DS for the same. But there is no extension / plugin found to enable this service in my Geoserver. Hence please help me to activate W3DS in step by step.Thanks in Advance.


texts: Annotations from different and completely unrelated QGIS projects appear

texts:
Annotations from completely separate QGIS projects appear in unrelated projects, not all created or necessarily accessible by me. This question was asked around 6 months ago here, a different use case, but may make sense to merge (I am new here):
QGIS showing text annotations from other geographically adjacent projects
To elaborate a bit further, not all the annotations appear, just a few random ones. Once removed, annotations can come back - but I cannot see a consistent pattern. I don't think the annotations disappear from the originating project (which is a life saver obviously), if deleted from the offended project. We usually create a design label folder with separate point layers to attach the annotations to, but sometimes even these randomly detach themselves (user errors may explain this, of course).
The particular use case is that I work with a number of different people on NextCloud. We share a small number of project assets, but none that would necessarily explain this problem (as I don't have access to the projects with the offending annotations). In a recent project I was working on my own NextCloud instance, unrelated to my co-workers, and again the same issue surfaced ... the image below is self explanatory I think, and just for completeness the annotation belongs to a project I have no access to, have never loaded.
If it were up to me I would not use annotations!


texts: Measuring accuracy of latitude and longitude?

texts:
I have latitude and longitude as 19.0649070739746 and 73.1308670043945 respectively.

In this case both coordinates are 13 decimal places long, but sometimes I also get coordinates which are 6 decimal places long.

Do fewer decimal points affect accuracy, and what does every digit after the decimal place signify?


texts: What vertical datum for OpenDroneMap with no GCP

texts: I am processing drone imagery with OpenDroneMap using a native Ubuntu installation of OpenDroneMap.  I see the the LAS generated returns points in UTM for X and Y coordinates as stated in the documentation.    The Z values are about 18 meters below the actual surface. I cannot explain this shift based on differences in vertical datums for my area of interest.  What is the vertical datum values produced by OpenDroneMap when no ground control point file is incorporated?


texts: How to delete multiple selected polygons at once?

texts: I'm new at this, a novice at best. I'm using QGIS 1.7.3. I simply want to select a large number of polygons on screen and hit the delete key to delete them and resave the shapefile. Can this be done?


texts:
ArcGIS Server Version: 10.4.0
ArcGIS JS API: 3.25
ArcGIS Desktop: 10.2.1
Browser: IE 11

In my web application, I am consuming secured services and have a requirement to print the map contents.

Initially, the token expiry was set to default (i.e. 60 mins) and everything worked fine including printing, using the Utilities/PrintingTools GP Service.

As part of the requirement, the token expiry was increased from 60 mins to 240 mins - everything still worked fine except for printing! Looking up the esri help it was clear I had to create a custom service. I created the custom service as described in the help documents, but it still fails with the following error:

  Error executing tool. Export Web Map Task Job ID:
  jb88972125ace4b6593c8055ee779e920 : Layer "dynLyrId": Unable to
  connect to map server at
  https://server.domain:6443/arcgis/rest/services/WebViewer/MapServer/.
  Failed to execute (Export Web Map). Failed to execute (Export Web Map
  Task).

The custom service has a user connection (with saved credentials) embedded into it as required.

Upon investigating further, I found out that if I remove the token from the WebMapAsJson, it works like a charm (in ArcCatalog at least)!

I don't know how to do that as I am not creating the WebMapAsJson manually, but the esri JS API creates it while sending the request. I also tried to increase the ArcGIS Servers' shortTimeout value from 60 to 240 but even that failed.

Any ideas to make the printing work with long-validity tokens?


texts: Raster Calculator dynamic map algebra expression

texts:
Is it possible to dynamically construct a Map Algebra expression that is dependent on the output of another tool, either using model builder or arcpy? For example, say I have a model that performs a raster reclassification that has a list of raster inputs and their respective outputs, i.e.

rasterOne -> reclassification -> outputRasterOne

rasterTwo -> reclassification -> outputRasterTwo

EDIT: To clarify, the number of input rasters is not known and therefore the number and names of the output rasters are not either. Because of this, I cannot hardcode them into the map algebra expression.

I would then want the map algebra expression to be similar to:

%"outputRasterOne"% + %"outputRasterTwo"% + ...


texts: Selecting new feature in attribute table of ArcMap without using mouse?

texts:
I am using ArcGIS 10.1 for Desktop.

When editing fields of a feature in an attribute table (right-click -> open attribute table with start editing enabled), is there a way to quickly select a new record without using the mouse?

I want the record to be selected so I can see it on the map in order to know what to change the field to.

I am not creating a new record.  I have added a new column to the table and want to add data to the fields in that column.  Before I add the data to the field I want to be sure I am editing the correct record, hence the selection highlighting in both the map and the table.  I can do this by double clicking on the field with the mouse, but it seems like it would be faster to stay on the keyboard.


texts: WebGL error when switching visualization

texts:
I am getting similar WebGL error messages when switching CARTO VL visualizations:

[.WebGL-0x59c22054800]RENDER WARNING: there is no texture bound to the unit 0

[.WebGL-0x59c222e5500]RENDER WARNING: there is no texture bound to the unit 1

These errors are being returned when switching from a current visualization to a new created one as shown in this example, but also when blending the visualization with a new filter as in this other example.

I have seen the last option working fine in this block. What is wrong with my examples code?

UPDATE: it looks like the problem is related to the (GeoJSON) data source. Replicating the same examples but using Dataset or SQL works (viz, filter examples).


texts: Migrating from Mapinfo to PostGIS/QGIS and preserving styles?

texts:
We're moving all of our Mapinfo tab files into a PostGIS database, and will be using QGIS instead of Mapinfo afterwards.

I've noticed that PostGIS records the tabs' style formats in a column named 'MI_STYLE', which preserves the styling when the tables are loaded back into Mapinfo. However when the tables are loaded into QGIS for viewing, it doesn't seem to recognise the styles.

Any thoughts on how to get around this?


texts: GEE converting to feature to geometry for .filterBounds

texts:
I am trying to filter a MODIS image collection by country. I have seen a lot of examples where you can map the image collection over a function and clip to the country, but I was trying to simplify the code by filtering my image collection by the country geometry (.filterBounds) instead. I know I need to feed in a geometry to do this, so I tried to convert my country feature to a geometry, and then filter by .filterBounds. It does not return an error, but also does not filter to my country shape. Is there a a way to get this to work? If not, why is creating a function to clip an image better? thanks so much for the help!
var canadaShape = ee.FeatureCollection('USDOS/LSIB_SIMPLE/2017').filter(ee.Filter.eq('country_na', 'Canada'));
var geometry = CanadaShape.geometry();
Map.addLayer(geometry);
var terraCollection = ee.ImageCollection("MODIS/006/MOD10A1").select(['NDSI_Snow_Cover']).filterBounds(geometry).filter(ee.Filter.date('2019-11-01', '2020-05-01'));
Map.addLayer(terraCollection);


texts: Cannot add ArcGIS ModelBuilder output to table of contents

texts:
I've made a model that creates some points with Z value and a buffer, which I then do a 3D analyst IDW on, and a Cut Fill analysis.

The model input is a line along with af DEM/Raster, and some names for the output files. However I would like to show the IDW output as a layer that will be opened in the Table Of Contents and displayed as well as the CutFill, but for some reason I can't get it to work.

Answer: Adding a Make Raster Layer after the IDW output seems to do the trick as well.

Link for the toolbox with the model in it


texts:
How can I solve this problem with my SLD code for styles. When I do two scale ranges for a one lineSymbolizer to use another color or another thickness they appear twice in the legend of my viewer because I have to separate them into two rules. Can you do this in the same rule?

        <Title>Rua</Title>

          <ogc:Filter>
            <ogc:PropertyIsEqualTo>
              <ogc:PropertyName>TIPO_LOGRADOURO</ogc:PropertyName>
              <ogc:Literal>RUA</ogc:Literal>
            </ogc:PropertyIsEqualTo>
          </ogc:Filter>

      <MinScaleDenominator>17000</MinScaleDenominator>
      <MaxScaleDenominator>36000</MaxScaleDenominator>

      <LineSymbolizer>

        <Stroke>
      <CssParameter name="stroke">#f6efc0</CssParameter>
      <CssParameter name="stroke-width">0.6</CssParameter>
    </Stroke>

      </LineSymbolizer>

     </Rule>

       <Title>Rua</Title>

         <ogc:Filter>
            <ogc:PropertyIsEqualTo>
              <ogc:PropertyName>TIPO_LOGRADOURO</ogc:PropertyName>
              <ogc:Literal>RUA</ogc:Literal>
            </ogc:PropertyIsEqualTo>
          </ogc:Filter>

      <MinScaleDenominator>8500</MinScaleDenominator>
      <MaxScaleDenominator>17000</MaxScaleDenominator>

      <LineSymbolizer>

        <Stroke>
          <CssParameter name="stroke">#e098a9</CssParameter>
      <CssParameter name="stroke-width">1</CssParameter>
    </Stroke>

      </LineSymbolizer>

The legend looks like

[]


texts: Spatial Join if more than 50% of join feature overlay target feature?

texts:
My ModelBuilder in ArcMap 10.3 includes 2 Shapefiles:

Buildings footprint (polygons).
Red shapes (polygons) which includes information for each building.

The task is creating new layer that contains Red shapes (target features) with corresponding footprints (join features), if more than 50% of "Red Shape" overlay building footprint.

I tried "Spatial Join" with match options INTERSECT, HAVE_THEIR_CENTER_IN and WITHIN_A_DISTANCE with negative value in search radius, but i can't achieve the desired result. For example in attached picture; Red shape No. 4 shall be included while the shapes 1, 2 and 3 shall be excluded since they contribution very small on building footprint.

How I would do that in ModelBuilder?


texts:
I am working with a .tif raster which has an attribute table. I can successfully add a number of extra fields, however I am having problems populating those fields.

When I add the updateCursor part to my code the script runs correctly but the attribute table is empty.

What I want to do is to take the value in the "rv_q20" field and if that is equal to a number write a new given value in "rv_q20s". I tried to pass the values from "rv_q20" between "" (e.g. "5") but encounter the same problem.

I would like to avoid using the arcpy.CalculateField_management

import arcpy

megaraster = r"E:\GIS Data\Raster1.tif"

#add score fields
arcpy.AddField_management(megaraster, "rv_q20s", "SHORT")
arcpy.AddField_management(megaraster, "rv_q50s", "SHORT")
arcpy.AddField_management(megaraster, "rv_q100s", "SHORT")
arcpy.AddField_management(megaraster, "rv_q200s", "SHORT")
arcpy.AddField_management(megaraster, "rv_q500s", "SHORT")
arcpy.AddField_management(megaraster, "rv_q1500s", "SHORT")

arcpy.AddField_management(megaraster, "sw_q20s", "SHORT")
arcpy.AddField_management(megaraster, "sw_q50s", "SHORT")
arcpy.AddField_management(megaraster, "sw_q100s", "SHORT")
arcpy.AddField_management(megaraster, "sw_q200s", "SHORT")
arcpy.AddField_management(megaraster, "sw_q500s", "SHORT")
arcpy.AddField_management(megaraster, "sw_q1500s", "SHORT")

#add field descriptions
arcpy.AddField_management(megaraster, "rv_q20d", "TEXT")
arcpy.AddField_management(megaraster, "rv_q50d", "TEXT")
arcpy.AddField_management(megaraster, "rv_q100d", "TEXT")
arcpy.AddField_management(megaraster, "rv_q200d", "TEXT")
arcpy.AddField_management(megaraster, "rv_q500d", "TEXT")
arcpy.AddField_management(megaraster, "rv_q1500d", "TEXT")

arcpy.AddField_management(megaraster, "sw_q20d", "TEXT")
arcpy.AddField_management(megaraster, "sw_q50d", "TEXT")
arcpy.AddField_management(megaraster, "sw_q100d", "TEXT")
arcpy.AddField_management(megaraster, "sw_q200d", "TEXT")
arcpy.AddField_management(megaraster, "sw_q500d", "TEXT")
arcpy.AddField_management(megaraster, "sw_q1500d", "TEXT")

RFscore = "RFscore"
SWscore = "SWscore"

arcpy.AddField_management(megaraster, RFscore, "SHORT")
arcpy.AddField_management(megaraster, SWscore, "SHORT")

#add field for total score
Totscore = "TotScore"
arcpy.AddField_management(megaraster, Totscore, "SHORT")

#add field for Risk Level
RiskLevel = "RiskLevel"
arcpy.AddField_management(megaraster, RiskLevel, "TEXT")

RFscoreName = "RFRetPer"
SWscoreName = "SWRetPer"

arcpy.AddField_management(megaraster, RFscoreName, "TEXT")
arcpy.AddField_management(megaraster, SWscoreName, "TEXT")

#call the update cursor to iterate through all the rows

rows = arcpy.UpdateCursor(megaraster)

for row in rows:
    fldDepth = row.getValue("rv_q20")
    if fldDepth is 1:
        row.setValue("rv_q20s", 10)
    elif fldDepth is 2:
        row.setValue("rv_q20s", 11)
    elif fldDepth is 3:
        row.setValue("rv_q20s", 12)
    elif fldDepth is 4:
        row.setValue("rv_q20s", 13)
    elif fldDepth is 5:
        row.setValue("rv_q20s", 14)
    elif fldDepth is 6:
        row.setValue("rv_q20s", 15)
    else:
        row.setValue("rv_q20s", 0)


texts: Serving clean transparency with mbtiles in GeoServer?

texts:
I'm running across an issue with the display of some historic line maps that I'm trying to serve as a WMS with GeoServer.

I have created a mosaic of georeferenced line drawings in ArcMap, and I want to serve it with a transparent background.  I put the single mosaic .tif in TileMill (dev version) and styled it there using the raster colorizer to only show values from 0 to 150 (the black/dark greys).  This produces a perfect black line drawing with a transparent background that looks great over an aerial.  I then export to a .mbtiles file and serve that through GeoServer using the Mbtiles plugin.

The main problem is that when viewing the WMS in ArcMap, "white noise" is visible around the lines, which is not in TileMill.  I unpacked the .mbtiles file with mb-util and the lines are clean there, so it seems to come from GeoServer somewhere.

I've tried modifying the color band details in the layer settings, but am not having any luck.

I think the problem was just with the default raster style in GeoServer, so I was able to create a new style and upload it to achieve exactly what I was looking for.  All is well!

I'm still refining the SLD that I will use for this, but I'm basically using the Alpha Channel effect that is described here.  Also, you can select which band to use as described here.

Download the full SLD.xml file available in the first link for the Alpha Channel
In GeoServer, go to the Styles page, and Add New Style
Upload the document, and you'll then be able to edit it.  I ended up with something simple like:

...
<RasterSymbolizer>
  <ChannelSelection>
    <GrayChannel>
      <SourceChannelName>4</SourceChannelName>
    </GrayChannel>
  </ChannelSelection>
  <ColorMap>
    <ColorMapEntry color="#000" quantity="175" opacity="0" />
    <ColorMapEntry color="#000" quantity="255" />
  </ColorMap>
 </RasterSymbolizer>
 ...
Be sure that the correct style is associated with your raster layer.
What's really nice about this method is that you can load the WMS in ArcMap/QGIS (or any webmap) and any changes to the stylesheet in GeoServer are immediately reflected in the WMS.


texts: Remove point cloud from Potree/Threejs-scene

texts:
I'm updating the cloud.js document while Potree is running on the client. When finished updating, I'd like to update the point cloud for the client. I tried to really update the client to re-read the cloud.js document, but cannot get it to work.

Now, as a workaround I'd like to remove the old pointcloud and add a new one using Potree.addPointcloud(). The second part works, but I want to delete the prevous version of the pointcloud. Is there anyone who can point me in the right direction.

I found something with the remove() method like scene.remove(pointcloud), but it's not recognized as a command.


texts: Is there a way to load a kml raster super-overlay in QGIS?

texts: I have created a set of tiled rasters as a super-overlay in Google Earth, and it performs remarkably well. Is there some way to load a super-overlay as a layer in QGIS? I am wondering if this could be an option to manage large numbers of rasters without setting up a WMS.


texts:
I'm trying to add a column in my point shapefile that lists points as degrees, minutes, not degrees decimal.

If it matters, I'm using 2.18.  The closest I've gotten is geom_to_wkt($geometry)


texts: "show only visible features" doesn't seem to work

texts:
I'm trying (unsuccesfully) to place the attribute table next to a map containing the same in print composer. Let's say for example my map in the primary window contains 6 districts and my map in the print composer only shows 1 district. If I unselect the "show only visible features" box then my table shows all entries (e.g. from 6 districts). However when I select the "show only visible features" it shows only the title row of the attribute table: it does not show any of the data entries.

What am I doing wrong?


texts: Creating permanent tile URL in Google Earth Engine

texts:
Is there a way to create a permanent tile server URL for an image?
Earlier I had used ee.Image.getMap() to get a URL like - https://earthengine.googleapis.com/v1alpha/projects/planet-webapp/maps/74f663849d42e3e8e71ff9a118942791-d7e8513d1d7c2b8cd15147beb65f2c17/tiles/{x}/{y}/{z}.png
What's the validity of such URLs?
I would prefer a URL for long term use once generated.
I am working on an open-source project which needs to show satellite imagery and analysis for multiple years, for 140+ planting projects. I need a way to automate the imagery generation process whenever a new location is added and save the tile URLs in the backend and use it.
This is how it will look on frontend - https://planet-webapp-jrckm5v45.vercel.app/yucatan-reforestation


texts: Does spherical equilateral triangle T contain point P?

texts:
I am working on a Python script to:

compute the vertex coordinates of a geodesic sphere/icosahedron,
project the triangles onto a sphere, then
find which spherical triangle contains an arbitrary point.

I would ultimately like to provide my script with a list of sample point coordinates, and the script would indicate which triangles contain points. These characters would then be used for presence-absence data in biological studies.

I'm mostly done, but struggling with step 3. Is there a fast algorithm to identify whether a point is contained within a spherical equilateral triangle? In two dimensions for general polygons, you can count how many times a ray crosses a perimeter, but this seems challenging to do in three dimensions since the perimeter is curved.

It seems like there should be a clever way to exploit the symmetry of spherical equilateral triangles, but I'm just not seeing it. I don't work with GIS or spherical trigonometry, so just let me know if my question is improperly phrased.


texts: Topology in QGIS shapefile, GRASS and v.clean (nested polygons)

texts:
Following this: Nested polygons in QGIS overlap and don't show the real colour gradient

I applied the GRASS function v.clean -> break as suggested by @ahmadhanb and the shapefile now shows correctly all the nested basins without using the Feature blending mode = Multiply.

However I noticed that in the attribute table the v.clean function created ~2 times more variables than the original shapefile (e.g. original shapefile 250 polygons, after v.clean 500 polygons). There are then duplicated polygons that eventually create issues when doing a Join.

How can I remove the duplicated polygons created by v.clean and keep ONLY the 'good' ones that appear in the QGIS environment when selected?

I tried to apply rmdupl within v.clean but the shapefile did not change at all!


texts: How to display layer name to GetLegendGraphic request WMS service?

texts:
I have a Mapping and WMS Service that I have published via ArcServer 10.3.1.

When communicating with the WMS service via a web browser, using a GetLegendGraphic Request. I am returned a graphic of my legend, such as:

In addition to this, I need the layer name at the top of the legend. Similar to what I see in my MXD:

Is there any way to get the layer name displayed when I use the GetLegendGraphic request?


texts: Polyline insertPoint issue using ArcGIS API for JavaScript?

texts:
Using ArcGIS API for JavaScript 2.5 I am trying to draw a path on an esri map.  with circles like so, but when I call Polyline.insertPoint it doesn't show the new segment:

I keep adding points as time progresses. My implementation is to have a Graphic object with a  Polyline geometry and a SimpleLineSymbol. I call the Polyline's insertPoint method when I need to add a segment. The problem I'm coming across is that when I call insertPoint, the point point is added to the array, BUT doesn't color in the line. The interesting thing is that if I remove the graphic from the graphics layer and then add it back on, then the new line segment shows up. This would be sufficient, except I need to have circles drawn on top of the line, so when I remove+add the Polyline it is drawn on top of the circles.

So in my current implementation every time I add a line point, I have to:

Call Polyline.insertPoint
Remove then add the graphic from the graphics layer
Remove + add all the circles on the path (because they just got hidden from the above action)

The problem with the above steps is that the line and circles flash on and off during this process which is very undesirable.

I feel like I am using the graphic objects provided by esri incorrectly because it's so hard to accomplish a seemingly easy task. Do you have any suggestions on an alternative implementation that has a good probability of not having these same graphical artifacts that I'm experiencing?


texts: Creating multiple address points for each unit using arcpy

texts:
I have a point feature class of addresses and an excel spreadsheet of all the suites associated with each address. I need to create a point for each suite belonging to an address. Is there a way, using arcpy, to scan the excel, match addresses, and then explode each point to create a new point for each suite?

I would like to create a tool so that I don't have to go through each record manually but i'm not sure where to begin. I've done some searching but most of the information I have found is about deleting duplicate points not creating them. Has anyone come across this sort of problem before or know of any good tools that do this?

I am using Arc 10.2.2


texts:
I am running v10 of ArcGIS for Desktop.
I am aware of the OTB templates provided with ArcGIS.
- However, I would like to supplement these layouts with additional layouts.

I am also aware of the Map Template Gallery, but this seems to be more around data models and I am after the actual print layout templates.

In particular I am looking for Local Government style templates, but the only downloads I can find seem to be again around the data model side of things as opposed to the actual print templates.

Has anyone seen (perhaps on ArcScripts) a bundle of templates that I can add to ArcGIS Desktop so that I can switch to them in my Layout View?


texts: OpenStreetMaps routing - starting far from a road

texts:
I am writing an web application, that needs routes between some markers (A, B, C, D).
So far I used the this api with PHP, where I pass each pair of points (A and B, B and C, C and D .... and so on) and loop through each pair to get the coordinates and build polyline. But some of the points are not on the street, like the starting point in this example.

So, technically is there a way to find if any of the points is not on the road?


texts:
I have the following code:

//qf  is a query filter
var fCursor = FC.Search(qf, true);
var feat = fCursor.NextFeature();
int selCount = FC.FeatureCount(qf);
var userId = new int[selCount];
if (selCount != 0)
{
    for (int i = 0; i < selCount; i++)
    {
        userId.SetValue(feat.OID, i);

        **Marshal.ReleaseComObject(feat);**

        feat = fCursor.NextFeature();
    }
    ...
}

Does it make any difference if we release IFeature (feat) outside of the loop?


texts: How to transform/convert Coordinate Reference System with Proj4 and php

texts:
Basically my question is:

Can I combine the functions of proj4 with php in my server?

For example, I would like to know if it is possible get data in WGS84 and transform in ETRS89.

I need to get WGS84 data with php, then use proj4 to transform them to ETRS89 and finally use them in php.

If the answer is right, how can I configure php and proj4 to work together???

Thanks!!!


texts: Oracle versioned geodatabase feature class field value update on multiple versions

texts:
I have an Oracle enterprise geodatabase with multiple versions.

In one field value of one point feature class I need to change all features within all versions.

I have more than 50 versions.

Is there any solution for this?

I am using ArcGIS 10.2.1 products.


texts: QGIS - Labeling points with numbers and adding text label to legend as list to represent what the numbers are

texts:
I am new to GIS and therefore have very limited knowledge of this software.

I have a list of reserves that have been put into the system. To label these I was wondering if there is a way of making the labels into numbers and then beside the map a key of what the numbers represent.


texts: It seems there is no geoalgorithm to drop z values in the QGIS processing toolbox. Can this be done programmatically using the Python console? At the moment I have to save the layer with unchecked 'Include z-dimension' checkbox (non automatic geometry type):


texts: Export virtual fields in attribute table to Excel using QGIS

texts:
I am using QGIS.

I am trying to export my attribute table to Excel but I want the virtual fields also to be exported with their function. For example I have a field that sums up two other fields and I want in my Excel file to keep the sum function and not only the context of the number.


texts:
I need to get a feature where i clicked on:

from PyQt4.QtCore import *
from qgis.core import QGis, QgsMapLayerRegistry
from qgis.gui import QgsMapToolEmitPoint, QgsRubberBand,  QgsMapToolIdentifyFeature

polylayer =QgsMapLayerRegistry.instance().mapLayersByName('Polys')[0]#shapefile
#save current tool
prevTool = iface.mapCanvas().mapTool()
#create tool
curTool = QgsMapToolEmitPoint(iface.mapCanvas())

def on_click(ptnClick, clickedButton):
    gpoint = QgsGeometry.fromPoint(ptnClick)
    pntBuff = gpoint.buffer((iface.mapCanvas().mapUnitsPerPixel() * 2), 0)
    rect = pntBuff.boundingBox()
    if clickedButton == Qt.LeftButton:
        print ptnClick
        index = QgsSpatialIndex(polylayer.getFeatures())
        # returns array of IDs of features which intersect the rectangle
        intersect = index.intersects(rect)
        print intersect
        request = QgsFeatureRequest().setFilterFids(intersect)
        for feature in polylayer.getFeatures(request):
            print "found:", feature.id()
        for feature in polylayer.getFeatures():
            gpoly = feature.geometry()
            if gpoint.within(gpoly):
                 print "Found within", feature.id()
            if gpoly.within(gpoint):
                 print "Found within", feature.id()
            if gpoly.intersects(rect):
                print "Found intersects rect:", feature.id()
            if gpoly.intersects(gpoint):
                print "Found intersects point:", feature.id()
            if gpoly.contains(ptnClick):
                print "Found contains point:", feature.id()

def tool_changed(tool):#another tool was activated
    iface.mapCanvas().mapToolSet.disconnect(tool_changed)
    curTool.deleteLater()

curTool.canvasClicked.connect(on_click)
iface.mapCanvas().setMapTool(curTool)
iface.mapCanvas().mapToolSet.connect(tool_changed)

But I get no result.
Any idea whats wrong or missing?


texts: Quantum GIS Standalone Installer Windows folder

texts:
Which is the "Quantum GIS Standalone Installer Windows" folder?
Is it, as usual, something like "c:\program files\quantum gis"?

I cannot make a test (at this moment) in a Windows machine.


texts:
I have geotagged pictures which I have linked in QGIS to a point shapefile. In the attribute table of the shapefile, you can find the complete path to every picture, so theoretically it is clickable. In QGIS I have an action which opens the photo.

Now I would like to create a 3D map with qgis2threejs where I have these points for the photos and the photos should be openable / clickable within the browser view.

So by clicking on the object of the 3D map in my web browser I would like to have a clickable path to the photo (opens by clicking) in the attribute table. Now I have only the path which I can copy and paste.

Even better would be if I could somehow put the whole picture into the object's attribute table or a clickable thumb nail of it.

Do I need to put the file path somehow in a clickable mode (hHTML?) after I have created the 3D map?

When I open my 3D map's JavaScript file with Notepad++ I see the following code for the point layer that contains the file path to the photo (excerpt):

// Layer 1
lyr = project.addLayer(new Q3D.PointLayer({q:1,objType:"Box",type:"point",name:"161020_Site_visit"}));
lyr.a = ["filepath","filename","DateTimeOr","GPSAltitud","GPSAltit_1","GPSDateSta","GPSImgDire","GPSImgDi_1","GPSLatitud","GPSLatit_1","GPSLongitu","GPSLongi_1","GPSMapDatu","GPSSatelli","GPSTimeSta","GPSVersion","ImageDescr"];

lyr.f[0] = {d:0.267671923173,h:0.267671923173,m:0,rotateX:90,w:0.267671923173,pts:[[-51750.4998443,-122555.375632,0.178447948782]]};

lyr.f[0].a = ["D:\\Photos\\DSCN2858.JPG","DSCN2858.JPG","2016:10:20 10:44:55","undef","0",null,"undef",null,"0",null,"0",null,null,null,"00:00:00","2 3 0 0",null];

Is there no way to manually (by adding some code like Photo1) make the photos clickable like an clickable URL, so that one could just click on D:\Photos\DSCN2858.JPG to open the photo?

I found this Stack Overflow Q&A that deals with a very similar question, but I don't know how to implement it for this purpose.

The problem seems to be that web browsers (this is where you watch the 3D map) don't allow top open local files (the photos). There seems to be a way to circumvent this, but apparently it requires more coding skills to not break the JavaScript file which comes with your 3D map.


texts: Reason for better save capability of QGIS 3

texts:
I have an internal processes meeting with work where I am trying to describe the benefits to migrating to QGIS 3 from QGIS 2.

I need to explain the technicalities of how the save option works better for some of the processing tools. In this example I use the dissolve tool:

QGIS 2.18.14

QGIS 3.2.3

I know that QGIS 3 is so much better because you don't have to painfully navigate through your file system and that you have access to the Quick Access Panel, but I just don't know what has changed technically for this to be possible.

Basically I need the technical jargon to explain to them why?


texts:
I have three rasters that have been reclassified into 5 classes each. I would like to use decision tree analysis, using the classes in the 3 rasters to come up with a final raster that has 5 classes, developed by using the classes in the 3 rasters as criteria.

How can I go about this?


texts: Grid creation with specific rows and columns

texts: I have a shapefile for a property boundary and i want to create a grid over this boundary. However I want to be able to specify the number of Rows and Columns, I would also like each individual grid to be able to hold information. I have tried using "Create vector index plugin", and vector grid but it has its problems. Is there any other plugin or application I can use to achieve this. I am using QGIS 2.18.11.


texts:
I am still a newbie at the Flex Viewer; I have been using the Out of the Box Web Application for Server 10.0.

I tried to put a simple service as an additional basemap (I also tried it as an additional operational layer) and now the viewer will not fully open up.

After I refresh the page multiple times it finally comes up but doesn't work properly.

The map service I have added is not cached but I don't cache it for the Out of the box web application.

The map service is relatively small geographically, and currently has no large feature classes in it or imagery. Would caching the map service change it that much??

I'd like to add imagery into but I want to get it working without it first. I would like to replace the web mapping application with the Flex viewer as there are more options.

The data server and the web server are on the same machine. The map service is also using UNC paths as I am creating the original map on a different machine but if that is huge problem that could be adjusted.

Trouble is the machine that has ArcGIS Server on it does not have Desktop and vice versa. Kind of limited with a single use license.


texts: How to add point features from longitude and latitude coordinates on QGIS

texts: Will it be possible to add a set of points from a table of longitude and latitude coordinates as show in the pic attached? Or will I have to convert them to X Y coordinates and then add it as a .CSV layer? if this is the case how do I go about the conversion.


texts: Measuring Average Distance Polygons Extend Past Line

texts: I need to find the average distance of a set of polygons that extend beyond a line for risk management purposes.  There are multiple polygons representing potential spill zones and different pipeline segments associated with these areas.  I am attempting to find an average distance from the furthest part of the polygon to the line.  This information does not need to follow any elevation contours or other geographic features.  I simply need a way to find the average so I can create a meaningful buffer for some data I need to clip.  Can someone suggest a method for doing so?


texts: How to Import ESRI Geodatabase format .gdb into PostGIS

texts:
I ran into a problem when loading a ESRI Geodatabase format .gdb into PostGIS. I have data 2.5GB GDB file.
I followed some tutorials on the internet, but it seems that doesn't work out.

I created a database "SampleNY"
I executed this from the console:

ogr2ogr -f "PostgreSQL" PG:"dbname=SampleNY user=postgres" NYPluto/Pluto.gdb

But nothing happens, I didn't got neither an error nor a successful operation.
Did I miss any steps?


texts: Seamlessly merging USGS DEM data with higher resolution survey?

texts:
I have a 10m USGS DEM for my entire map extent (~100sq miles)

Inside of that, I have specific regions of higher quality GPS point data. (From a rover using a OPUS corrected RTK base station)

How do I seamlessly merge my data/data products with the "background" USGS data? If I simply merge/replace the USGS data with mine, the edges don't align and my entire region looks "sunk in" or "raised up" in comparison to the USGS "background"

The gaps that are being filled in by the USGS DEM are up to a full square mile, so I can't just interpolate the blank area from the edges of my survey.

My first though is to "feather" the USGS DEM. Obviously at the edge where the two meet, my data is more accurate, so feathering the USGS data to meet smoothly with my survey seems like it should work, but I am unsure of the technical way to approach this or the accuracy of this method.

The data products that are being produced from this are hillshades, viewsheds, streamflow, contours, and more - and none of them seem to work well with a sharp boundary transition.


texts: Carry over attributes from a point layer to a new buffer layer using Python in QGIS

texts:
Using the code outlined in the question: How to create square buffers around points in QGIS with Python? I created a layer of square buffers from coordinate points, however I was unable to carry over all of the point attributes to the new layer.

We have the following columns in our point attributes:

BigCellID (integer),
BigCellTyp (integer),
x,
y

The following code carried over the x and y columns perfectly:

layer = iface.activeLayer()

feats = [ feat for feat in layer.getFeatures() ]

epsg = layer.crs().postgisSrid()

uri = "Polygon?crs=epsg:" + str(epsg) + "&field=id:integer&field=x:real&field=y:real&field=point_id:integer""&index=yes"

mem_layer = QgsVectorLayer(uri,
                           'square_buffer',
                           'memory')

prov = mem_layer.dataProvider()

for i, feat in enumerate(feats):
    point = feat.geometry().asPoint()
    new_feat = QgsFeature()
    new_feat.setAttributes([i, point[0], point[1], feat.id()])
    tmp_feat = feat.geometry().buffer(1000, -1).boundingBox().asWktPolygon()
    new_feat.setGeometry(QgsGeometry.fromWkt(tmp_feat))
    prov.addFeatures([new_feat])

QgsMapLayerRegistry.instance().addMapLayer(mem_layer)

However to get it to also carry over the other two columns I made the following modification to one line:

uri = "Polygon?crs=epsg:" + str(epsg) + "&field=id:integer&field=x:real&field=y:real&field=BigCelID:integer&field=BigCellTyp:integer""&index=yes"

instead of:

uri = "Polygon?crs=epsg:" + str(epsg) + "&field=id:integer&field=x:real&field=y:real&field=point_id:integer""&index=yes"

This modification did not affect the x and y columns being carried over, and did create new columns; but the column BigCellID was filled with sequentially increasing numbers from 0, while BigCellTyp was filled with 'Null'. These are not the original values.

I would like to know how I can modify this code to ensure the contents of my two columns BigCellID and BigCellTyp are carried over to the new buffer layer?


texts: Getting table of coordinates from polygon in QGIS

texts:
I would like to gather all of the GPS coordinates within a set area as a table to transfer it to Excel or sheets.
I designed a rough shape of the area with the polygon tool in the layer and now im stuck.


texts: Raster projection changes to none when loading into MapWindows after already being defined in ArcMap

texts: I am trying to load a number of raster files (.tif) into MapWindows v.4.8.8. These rasters have previously been reprojected in ArcMap to have the projection WGS_1984_UTM_Zone_60S and saved with this projection. When I load these rasters into MapWindows their projection comes up as "none" and I can't run any analysis on them as a "no spatial reference information" error pops up. I am very new to using MapWindows but I need to use it for a specific plug-in.


texts: Using R to do logistic geographically weighted regression(GWR) prediction

texts:
I'm using the R package "GWmodel" to do the logistic GWR, but this package doesn't provide the prediction function for generalized GWR(Binomial/Poisson). How can I do the Logistic GWR prediction using R? I've searched other R packages, but there were no good solutions.

Also, I tried to apply Python library “PySAL”, this library is, it seems, providing logistic GWR prediction. However, the difference of the bandwidths calculated by "PySAL" and "GWmodel" was too large(using the same dataset), and the bandwidth from "GWmodel" was more reasonable, so I'm still considering to work on R. I don't want to use GWR4 too.


texts:
I have a local Geoserver connected to a PostGIS database and a table with an attribute field called "style" which is of data type integer. On the OpenLayers preview of the Geoserver the data is shown and accessible.

The web application I load via an local Apache (XAMPP) server. The same policy origin is considered through an cgi proxy script which is linked within the script.js. with this code line "OpenLayers.ProxyHost = "/cgi-bin/proxy.cgi?url=";". In the proxy script I added localhost for the Apache server and localhost:8080 for the Geoserver.

The problem is that i want to style a WFS over the attribute style, and every style: 1, 2, 3, ... should get an own color and a own checkbox.

For now only the default rule "rule_3" works and I don't know why. Here is the part of the source code:

// Marker Clustering for WFS Layer
var strat = [
    new OpenLayers.Strategy.BBOX(),
    new OpenLayers.Strategy.AnimatedCluster({
        distance: 15,
        animationMethod: OpenLayers.Easing.Expo.easeOut,
        animationDuration: 10
    })
    // new OpenLayers.Strategy.Cluster()
];

var style = new OpenLayers.Style();

var rule_1 = new OpenLayers.Rule({
    filter : new OpenLayers.Filter.Comparison({
        type : OpenLayers.Filter.Comparison.EQUAL_TO,
        property : 'style',
        value : 1
    }),
    symbolizer : {
        fillColor : '#ffffff'
    }
});

var rule_2 = new OpenLayers.Rule({
    filter : new OpenLayers.Filter.Comparison({
        type : OpenLayers.Filter.Comparison.EQUAL_TO,
        property : 'style',
        value : 2
    }),
    symbolizer : {
        fillColor : '#ff00ff'
    }
});

var rule_3 = new OpenLayers.Rule({
    symbolizer : {
        fillColor : '#ff6600'
    }
});

style.addRules([rule_1, rule_2, rule_3]);

var wfsLayer = new OpenLayers.Layer.Vector(
    'Christmas Traditions',
    {
        strategies: strat,
        styleMap: new OpenLayers.StyleMap(style),
        protocol: new OpenLayers.Protocol.WFS({
            version: '1.0.0',
            url: 'http://localhost:8080/geoserver/wfs',
            featurePrefix: 'cite',
            featureType: [
                'op_god_jul_fts',
                'op_vrolij_kerstfeest_fts'
            ],
            featureNS: 'http://www.opengeospatial.net/cite',
            //propertyName: ['txt', 'style'],
            srsName: 'EPSG:4326',
            geometryName: 'coords'
    }),
    isBaseLayer: false,
    visibility: true,
    });

// Add layers to map object
map.addLayers([
    cloudmadeMidnightCommander,
    stamenToner,
    mapQuestAerials,
    wfsLayer
]);

Current all points have the same color of style rule_3. The dots over Belgium and Netherlands should have an other color. The table "op_god_jul_fts" should have the style of rule_1 and the table "op_vrolij_kerstfeest_fts" should have the style of rule_2.


texts: Recording all vertices on shortest path with igraph?

texts:
I am working with a river network where vegetation surveys were conducted (black dots). I would like to:

calculate the minimum distance between all points and the point highlighted by the red arrow
for every shortest path between a point and "the red-arrow point" record all points located on that path. For example, the shortest path between the point highlighted by a circle goes through the points highlighted by  rectangles. Therefore, I would to like to to record the unique ID of the "rectangle points".

I have been trying to use shp2graph and igraph to find a solution but so far I was only able to calculate pairwise distances.

This is what I managed to do so far:

library(igraph)
library(shp2graph)
library(rgdal)
library(tidyverse)

#Import rivers and points into R

rivers<-readOGR(dsn="/path/to/test_area/",layer="river")
points<-readOGR(dsn="/path/to/test_area/",layer="points_all")

#Assign the points to the river network

ptsxy<-coordinates(points)[,1:2]
ptsxy<-cbind(ptsxy[, 1],ptsxy[, 2])
res.nv<-points2network(ntdata=rivers, pointsxy = ptsxy,Detailed=TRUE, ELComputed = TRUE,approach=2,ea.prop = rep(1, 2))

#Convert to igraph object

igr2 <-nel2igraph(res.nv[[1]],res.nv[[2]], weight = res.nv[[8]])

#Calculate minimum distances between every pair of points

id_num <- as.numeric(res.nv[[3]])
mat_dist<-distances(igr2,v=id_num,to=id_num)

The result is a matrix of minimum distances between points.

Here's a link to the data I am currently using:

https://www.dropbox.com/s/zvjti3b7zp08ruy/sample_data.zip?dl=0


texts: For what tools in ArcMap is the Export Script option active?

texts:
What tools allow you to export scripts?

I am assuming that these script would be blank python scripts.  I have not found a single tool where this option is active in my installation.


texts: edge Data labelling OSM local tile server

texts:
i have OSM tile server running locally on my machine

I have used Leaflet to add markers

Now I am looking for graph script in which markers to act as nodes and polygon (more precisely an arc) to show edged data? is there any API available?


texts: Bringing State Boundary closer to land using QGIS?

texts:
I am new to QGIS.

I downloaded a data from a website and tried uploading it to the base map plugin which would be OpenStreetMap.

However, the boundary is not close to the land.

For example, like this:

Is there a way to bring the boundary closer to the land?

I have also used other base map plugins such as Google and Bing. It is still causing the same issue.


texts: How to create an oriented buffer using arcpy?

texts:
I would like to create an oriented buffer for every polygon in my shapefile using arcpy. By oriented I mean that I have two angles a1 and a2 that constraint the direction of the buffer. This is represented in the graph below:

Any ideas?


texts: Editing QGIS point layer to retrieve attributes of intersecting line layer?

texts:
I'm creating a custom editing form using QGIS 2.0 and would like to modify the mouse event when creating a point. I've gotten a start by reading Programatically check for mouse click in PyQGIS?

In fact, the code below is copied and modified directly from it.

My goal: When creating a new point, I would like to retreive feature attributes of another (line) layer that intersects with the mouse click, using a rectangle filter. The mouse click opens the form when creating the point, while automatically retrieving data from another line layer intersecting with the rectangle buffer. This is in the interest of comparing data in these two layers that are supposed to correspond.

Following the lead of the mentioned Q&A, i'm attempting to make use of the following code without creating a separate Map Canvas. To put it to use, what would I provide, if anything, as the Python Init function under Layer Properties > Fields? Currently I have mypyfile.formOpen as I've added python logic as detailed here: http://nathanw.net/2011/09/05/qgis-tips-custom-feature-forms-with-python-logic/

I may be going about this in the wrong way, so i'm open to suggestions to direct me towards creating this data collection tool.

class PointTool(QgsMapTool):
    def __init__(self, canvas):
        QgsMapTool.__init__(self, canvas)
        self.canvas = self.iface.mapCanvas()
        self.iface = iface
        self.pointTool = QgsMapToolEmitPoint(self.canvas)

    def run(self):
        self.canvas.setMapTool(self.pointTool)
        QObject.connect(self.pointTool, SIGNAL("canvasClicked(const QgsPoint, Qt::MouseButton)", self.selectFeature))

    def selectFeature(self, point):
        #resLayer previously defined
        if not resLayer or resLayer.type() != QgsMapLayer.VectorLayer():
            QMessageBox.warning(None, "No!", "Select a vectorlayer")
            return

        width = self.iface.mapCanvas().mapUnitsPerPixel()*2
        rect = QgsRectangle(point.x() - width, point.y() - width, point.x() + width, point.y() + width)
        rect = self.iface.mapCanvas().mapRenderer().mapToLayerCoordinates(layer, rect)

        resLayer.select([], rect)
        feat = QgsFeature()
        ids = []
        while resLayer.nextFeature(feat):
            ids.append(feat.id())
        #selecting the feature or features that intersect with mouse click buffer
        resLayer.setSelectedFeatures(ids)

#code continues to guide user data entry using python logic


texts: Area calculation in QGIS done two different ways?

texts:
It seems that there are at least two different ways to calculate the area of a polygon in QGIS using Python:

1) area = feature.geometry().area()
2) d = QgsDistanceArea()
   area_v2 = d.measurePolygon(feature.geometry().asPolygon()[0])

These two approaches produce very similar outcomes, but not the same. It also seems that QGIS prefers the second approach -once the geometry column is added through Vector -> Export/Add Geometry Columns. Presumingly, this is related with the projections. However, everything I deal on QGIS is on EPSG:3857 since I am using OSM data. Furthermore, the perimeter values obtained from this approach match with the first approach when the following statement is executed:

 perimeter = feature.geometry().length()

I am inclined to choose the first approach due to its simplicity, so why does QGIS rely on the second approach?


texts: Installation of QGIS 3.4.4 on Ubuntu 16.04 fails to install qgis-plugin-grass because of uninstallable grass760 dependency

texts:
I have used the instructions  here marked Tested 29/06/2018 to install QGIS 3.4.4 on Ubuntu 16.04.  While QGIS installs nicely I get an error on installing qgis-plugin-grass.  Subsequently the Grass tools listed in the processing window all throw errors.

Upon installation I get the following error...

The following packages have unmet dependencie: qgis-plugin-grass :
Depends; grass760 but it is not installable.
E: Unable to correct problems, you have held broken packages.

I have tried reinstalling using aptitude but this did not make any difference.  How do I install the grass760 dependencies and qgis-plugin-grass?