Most Popular

1500 questions
23
votes
1 answer

QGIS 3.0 point cluster renderer, is it possible to control rendering order (z-order)?

I've been having a play with the new point cluster renderer in QGIS 3.0.2 and it's great - a lot easier (and faster) than having to do the clustering yourself. To show what I mean, here's an example of road accidents in Edinburgh using STATS19 data…
Steven Kay
  • 20,405
  • 5
  • 31
  • 82
23
votes
4 answers

Using OSGeo4W shell with Python 3

I would like to use OSGeo4W shell with Python 3 but when typing python3 I get the following error: Fatal Python error: Py_Initialize: unable to load the file system codec File "C:\OSGEO4~1\apps\Python27\lib\encodings\__init__.py", line 123 …
bennos
  • 2,069
  • 1
  • 15
  • 29
23
votes
2 answers

equivalent of sp package's point in polygon / overing using sf

I am migrating code from sp package to the newer sf package. My previous code I had a polygon SpatialDataFrame (censimentoMap) and a SpatialPointDataFrame (indirizzi.sp) and I got the polygon cell id ("Cell110") for each point laying within with the…
Giorgio Spedicato
  • 503
  • 2
  • 4
  • 9
23
votes
7 answers

GIS Software Choice for a small university research centre

I need to choose a GIS system for a small university research centre. We are handling a broad range of data, (for example, numerical tidal analyses, weather data, poverty, isolated economic activity, skills availability and renewable energy resource…
Alan Owen
  • 239
  • 1
  • 3
23
votes
3 answers

How to calculate the angle at which two lines intersect in PostGIS?

I want to calculate the angle between two lines where they intersect in PostGIS. The starting point for angle calculations in PostGIS seems to be ST_Azimuth - but that takes points as input. My first thought was to take the endpoints of the…
mvexel
  • 1,680
  • 13
  • 30
23
votes
4 answers

Layer stacking rasters in QGIS?

I have some multi-spectral data, where the individual bands are separate Raster files. How do I layer stack these rasters so that I get one raster with 3 bands instead of 3 different rasters? This is quite easy to do in ERDAS and ArcGIS, but I…
Devdatta Tengshe
  • 41,311
  • 35
  • 139
  • 263
23
votes
6 answers

Converting shapefile to CSV including attributes and geometry

I have a shapefile with 60k+ entries, all of which are polygons with corresponding attributes (acreage totals, landowner names, tax ID #s, etc.). What I ultimately need is a CSV file with all of these attributes and their corresponding geometry (in…
Rikk
  • 231
  • 1
  • 2
  • 3
23
votes
1 answer

Displaying numbers of single band raster layer in QGIS without additional vector layer?

I wonder if there is a raster symbolizer in QGIS with a label option, that renders the cell values as labels of single band raster layer, a scalar field. The standard raster symbolizer dialog in QGIS has no label option. So I use two data sources,…
huckfinn
  • 3,583
  • 17
  • 35
23
votes
7 answers

Alternatives to Shapefiles as open-source, cross platform dataset types

I am working on software that is very ESRI oriented, but a future version will likely not be using ESRI software. It uses Shapefiles and Geodatabases. I'm planning on getting all of my data to Shapefiles in anticipation for future versions of the…
Tanner
  • 7,375
  • 14
  • 62
  • 109
23
votes
12 answers

Discontinuous or Interrupted map projections

Do any of the current open source GIS softwares support discontinous projection systems e.g. sinusoidal, dymaxion, butterfly, etc? If not, are there any affordable commercial ones that do (affordable meaning < $150)? Apparently it was a bad idea to…
matt wilkie
  • 28,176
  • 35
  • 147
  • 280
23
votes
2 answers

Can qgis_core work with external python distribution such as Anaconda?

I've built a Python application that has its own GUI. I am trying to use qgis.core/PyQGIS as map display only. Specifically, I have anaconda 2.3 and its Python distribution (2.7.11). It seems that this scheme can only work if the application I am…
For Comment
  • 1,115
  • 7
  • 21
23
votes
3 answers

Smoothing/interpolating raster in Python using GDAL?

I am developing in Python and using GDAL from OSGEO to manipulate and interact with rasters and shapefiles. I want to take a shapefile that has point features and interpolate it into a surface raster. Right now I am using the method…
Doug
  • 617
  • 1
  • 6
  • 10
23
votes
5 answers

Android devices and apps with support for Galileo positioning?

While the European's Galileo positioning system is only expected to be fully deployed by 2020, the first positioning made using already deployed Galileo satellites was made back in 2013. Also, hardware devices supporting Galileo are already…
augustin
  • 369
  • 1
  • 2
  • 10
23
votes
5 answers

Obtaining extent of each polygon in shapefile using ArcPy with ArcMap

In ArcGIS 10 and Python I want to get the extent (xmax, ymax, xmin, ymin) info of each of the polygons in a shapefile. I can get the extent of the whole shapefile using file=r"D:\SCRATCH\ARCGIS\100k_trc_tiles_TVM.shp" desc=arcpy.Describe(file) print…
GeorgeC
  • 8,228
  • 7
  • 52
  • 136
23
votes
4 answers

Cut/Clip/Split a polygon with another polygon from same layer

To not have overlapping polygons/features, I normally use one or more polygons to clip/shape other polygons. This is often much faster than tracing when adding polygons etc. I cannot find this functionality in QGIS. The highlighted polygon is…
Johan Fast
  • 319
  • 1
  • 2
  • 5