All Questions

1
vote
0answers
23 views

Output layer doesn't get created

Why isn't my heat layer valid and I don't get a .tif in this location (this location exists 100%)? Code below: output = "D:/QGIS_Project_New/heat_map_last.tif" heat = QgsRasterLayer(output, "Heatmap")...
-2
votes
2answers
23 views

Seeking data for spatial interpolation?

Can someone point me to a Kriging dataset? If it is a natural data, it would be great! (Such as rainfall, porosity, soil-density, atmosphere-concentration, etc. ), temperature observations.
0
votes
0answers
18 views

SDO Oracle Table not recognized as Geometry in QGIS

I can connect to the Oracle Database with QGIS but spatial table/layers are not recognized as such. I have read the following in the QGIS documentation: Normally, an Oracle Spatial layer is ...
1
vote
0answers
9 views

Dissolve Line Strings to one string - “No Shapely geometry can be created from null value”

I have a layer with 7000 LineStrings, all represent bus routes: It seems that all the bus routes are contiguous, but I'm not sure of it. I'm trying to use python shapely to dissolve all those ...
0
votes
0answers
11 views

Convert DWG to Shapefile Using FME and R

I converted a DWG file to an ESRI Shapefile using FME Workbench. Now I wish to create an elevation raster. I'm not sure how to extract the elevation info. library(sp) library(raster) # https://www....
2
votes
1answer
15 views

Raster calculator - set default output path - QGIS 3.4.6

Raster calculator: How can I set the default output folder? By default it's set to "C:\Program Files\QGIS 3.4\bin" I'm using QGIS 3.4.6
0
votes
0answers
16 views

Error in FillSinks with Python

I'm trying to use FillSinks from ArcHydroTools with Python My code is: >>> import arcpy >>> import ArcHydroTools >>> fill = arcpy.GetParameterAsText(0) + '\\fill' >>...
0
votes
1answer
25 views

Open the popups of all markers that are visible on the map with a zoom > 10

I need to open all the markers that are visible on the map at Zoom > 10. I am using leaflet.markercluster plugin. Init map: initMap() { this.map = L.map('map', { center: [55.55, 37.61], ...
1
vote
0answers
24 views

Changing the field colour/design on a QGIS 3 QT attribute form?

I have problems with changing the physical appearance of a custom attribute form (QGIS 3, Windows 7). A custom attribute form can be designed, for use in QGIS, using QT / QT Designer. In QT Designer ...
0
votes
0answers
8 views

org.geotools.data.wfs.internal.WFSException [on hold]

I upgraded recently to geotools 21 since am working with java11 and I'm getting now the following error when doing anything on my FeatureCollection which wasn't there using the older version of ...
0
votes
0answers
13 views

Problem with Fill Sinks - ArcHydroTools

I'm getting use FillSinks of ArchHydroTools with Python but some problems are happen My code is: Import arcpy Import ArcHydroTools ArcHydroTools.FillSinks('projectraster.tif','fil.lyr', "10", "...
-1
votes
0answers
11 views

How to setup geofencing in MySQL? [on hold]

Right now I am working with PostgreSQL database and want to work with the MySQL database to setup geofencing. How can I do that?
3
votes
1answer
33 views

Cheapest path through polygons and lines from point to point

I want to find the cheapest path or the least cost path from one point to another point. The points are part of a vector map and do not lie on a network. All features being part of the map have ...
1
vote
0answers
13 views

Raster correlation: Error (arg6) in raster.modified.ttest (R)

I am performing pairwise spatial correlation of 6 rasters using the raster.modified.ttest() function in R, as per this very useful post. I have consistently had the following error with my data; ...
2
votes
1answer
26 views

Central deployment of startup.py in QGIS

In QGIS 3 on Win10 we use a startup.py script for setting up e.g. database connections and a couple of other things, residing in C:\Users\me\AppData\Roaming\QGIS\QGIS3 What we are after is to ...
0
votes
0answers
16 views

Modis Reprojection Tool (MRT) software

I am trying to use MODIS (Aqua) images but I want it to be geo-referenced using the MODIS Reprojection Tool. The problem is USGS has recently removed it from their website. Does anybody here have a ...
1
vote
0answers
9 views

Creating Scale field based on polyline length (ArcMap)(Python)

I'm playing around with data driven pages, and have a polyline feature class that I'm using as my index. Now, because the lines are all of variable lengths I can't just use the default settings as ...
0
votes
0answers
11 views

How to use PostGIS to get the exact polygon outline of a pointcloud patch (pcpatch) [on hold]

How can I use PostGIS or PostgreSQL to get the exact polygon outline of a pointcloud patch (pcpatch)?
4
votes
2answers
31 views

Distance to a ridge as function of angle (Python)

No success googling this for a while, but likely I just don't know correct terminology. Also, tell me if you see another exchange more suited for this. Given: A boolean 2-dim array that stores if a ...
0
votes
0answers
6 views

Repoint enterprise geodatabase connection using ArcPy in ArcGIS Pro

I have an ArcPy script that changes the connection parameters for every item in a layer file from a SQL Server connection to a PostGIS connection. This script works perfectly in ArcGIS Python version ...
0
votes
0answers
22 views

Leaflet markers at same geo location (ClusterGroup)

I have a map created with leaflet with different markers. Sometimes the markers have the same (1:1) geo position. I have already tested the great addon clustergroup, but I need the function to be ...
1
vote
0answers
28 views

Download data on click via Leaflet webmap [on hold]

I am working on a webmap which will allow users to click on a portion of the globe and download associated river data. Here is my code so far. <!DOCTYPE html> <html> <head> <...
2
votes
1answer
39 views

Get OGR driver Name for QGIS VectorLayer

Is it possible in the Python script to get the name of the OGR driver (for example, "Esri Shapefile" or "MapInfo File") for the QGIS vector layer? __ QGIS 2.18
0
votes
1answer
22 views

is it possible to show all labels of WMS layer at any extent without SLD ol5?

I am showing WMS layer that shows label on it by defining SLD in GeoServer. The problem is that is is not showing all labels at particular extent. But when I zoom in, all labels are displayed. My ...
0
votes
0answers
11 views

ArcGIS is importing orginal numpy while VirtualEnv is activated

I could activate ArcGIS's python command windows to use VirtualEnv by "execute execfile(r'C:\Users\jason\myenv\Scripts\activate_this.py', {'file': r'C:\Users\jason\myenv\Scripts\activate_this.py'})" ...
0
votes
0answers
19 views

Print elevation of different coordinates from .hgt file

I am trying iterate through the data(elevation) from .hgt file. I have attached the codes which i used to get the elevation from the file. from gmalthgtparser import HgtParser hgt_file='N09E076....
0
votes
0answers
9 views

Parameters of the Muskingum-Cunge Model in HEC-HMS

I'm trying to simulate the flow of a river in HEC-HMS and for the routing model, I have chosen to use Muskingum-Cunge. In preparing the model, I used HEC-GeoHMS. When I got to check the parameters ...
1
vote
0answers
11 views

How to save user drawn layers into files(.shp or .tiff)? [on hold]

I am developing a wpf gis application using arcgis runtime sdk. In my application user can define FeatureCollectionTable schema and draw features(such as polyline,polygon,point and etc) on map. Now, I ...
0
votes
0answers
5 views

Get Feature Service syncCapabilities with ArcPy

I'm working on an ArcPy script to test layers. I have either an arcpy mapping module layer object, or a feature service URL as a text string. I wish to get the feature service's capabilities, in ...
0
votes
0answers
10 views

Using output of ArcGIS Euclidean distance tool?

How do I use the output of the ArcGIS Desktop Euclidean distance tool? It gives raster output and I want to generate a table enlisting the distance between various centroids from one centroid.
-1
votes
1answer
17 views

Linear model in R - Multiplication Expression [on hold]

I have 3 numerical variables A, B and C. I am trying to create a linear model capable of predicting A. The expression that I am using is the product of B*C in order to predict A; however, when looking ...
0
votes
1answer
10 views

Using Ground Control Points in Agisoft Photoscan?

I'm fairly new to the Agisoft Photoscan program, but I am trying to build several models using aerial imagery and highly accurate GCPs, with coordinates collected in the field. I align my photos and ...
2
votes
1answer
26 views

Transform CRS of Shapefile in R

I started with an AutoCAD DWG file and used FME Workbench to extract elevation data into a shapefile. read_sf yields: > shapename Simple feature collection with 41620 features and 1 field (with ...
2
votes
2answers
140 views

Apply Line Pattern Fill to LineStrings?

I am trying to print out a network map using black/white hatching patterns, similar to Black/white hatching pattern fills in QGis but for Lines. Is there some trick I can use?
-1
votes
1answer
30 views

Taking data from a CSV file to create points on a shapefile

I am trying to import data from a CSV file and write a script that will generate points on a new shapefile. I am running into issues for when I try to update the InsertCursor. I get an error of "not ...
0
votes
0answers
11 views

ArcGIS Online Custom Label Expressions

I'm not familiar with Arcade, but I'm trying to do a simple expression. I need to write an if-else label expression that will label addresses by their address number, unless there is a address unit ...
0
votes
0answers
11 views

Jason-3 Data Mapped in ArcGIS

Has anyone figured out an easy way to convert Jason-3 data in NetCDF format to a GIS format (raster, shapefile, etc.) so it can be viewed in ArcGIS? Jason-3 data is available for free in NetCDF format,...
0
votes
0answers
11 views

XML to TIFF raster and zonal statistics in QGIS

I was sent several XML files of a flood model's peak extent that works fine in ArcGIS, I need to use zonal statistics, and I don't have that available in my workplace's ArcGIS subscription, so I need ...
1
vote
0answers
11 views

Emit map click using Esri 4x JS API

I am trying to emit a map click from the MapView object using its emit method. I have successfully done this in the 3x API from the Map object. This same method, however, is not working for me in the ...
1
vote
0answers
10 views

QGIS quit unexpectedly

Running QGIS 2.18 on Mac OS 10.14. Python 2.7 and GDAL 2.1 installed. Trying to clip raster by polygon. (raster--> extraction --> clipper) Get error message. ("qgis 2.18 clipper The process failed ...
0
votes
1answer
25 views

How can I export multiple, consistently sized raster files?

I'm working with a computer that uses multiple ASCII raster files and ran into an issue with the file alignments. As you can see in ArcMap, the rasters are being aligned correctly. However, upon ...
0
votes
1answer
18 views

Adjusting symbology in layer file?

I have a layer file for a geology layer containing attributes for "Age", "Rock Type" and "L_Code". The layer has symbology defined using the "L_Code" attribute. But I want the same symbology applied ...
1
vote
0answers
9 views

Cannot find Grass7 folder in QGIS 2.18 content folder

I have freshly re-installed QGIS 2.18.28 in the prescribed order. I was under the impression Grass7 in included in one of the install packages. But when trying to point QGIS to the grass7 folder under ...
0
votes
1answer
29 views

Is saving rasters to VRT bad convention?

I'm working with a lot of rasters (dozens of thousands) and would like to store them for a while. They will be used intermittently, but never all at once! The size of VRT (~2 kb) is such a positive ...
1
vote
1answer
42 views

Raster Volume Returning Wrong Volume

I have a clipped surface using Washington State Plane North (feet) added to a QGIS project using the same CRS. I am trying to calculate the area under this surface (a dirt pile) using SAGA Raster ...
0
votes
0answers
6 views

VeloxRaster_crop returns numeric

Curious as to why the following only returns a numeric vector with row and column information rather than a cropped raster: ras <- velox("dir.tif") cty <- st_read("cty.shp") ras <- ras$crop(...
1
vote
1answer
36 views

How to split vector layer in PostgreSQL/PostGIS

QGIS has this nice tool to Split a layer into multiple layers based on one attribute and create a new layer for each unique value within that attribute. Is this also possible using PostgreSQL/PostGIS?...
-1
votes
1answer
20 views

Working with Raster data/DEM (mosaic to new raster, viewshed) in ArcGIS Desktop?

I'm trying to work with raster data from this source: http://www.soest.hawaii.edu/coasts/data/maui/index.html I downloaded the DEMs from different islands and imported them into ArcGIS. But I can't ...
2
votes
1answer
42 views

How to use QgsRasterInterface.histogram?

This is really a question of how to better understand QGIS API documentation (C++ or Python) but I'll use this example: How do I use QgsRasterInterface.histogram? I was told in this thread I can use ...
0
votes
0answers
8 views

Open NITF containing multiple images using Python GDAL

I have a file in NITF 2.1 format, which contains multiple images within it (in ENVI speak "image segments"). But with Python's GDAL, opening the file only gives me access to the first image. How can I ...

15 30 50 per page