Questions tagged [gis]

A Geographic Information System (GIS) captures, stores, analyzes, manages, and presents data that are linked to location(s). The GIS Stack Exchange should be considered for questions with this tag.

According to Wikipedia:

A geographic information system (GIS), geographical information system, or geospatial information system is any system that captures, stores, analyzes, manages, and presents data that are linked to location(s). In the simplest terms, GIS is the merging of cartography, statistical analysis, and database technology.

There is a Stack Exchange site devoted to Geographic Information Systems.

6555 questions
2
votes
0 answers

Convert GeoTiff image to Python Shapely polygon geometries representing 3 dimensional triangles

I am working on implementing the example at Quantized-Mesh-Tile on my own data sets which are in .Tif (Digital Elevation Models) format. According to the module, it accept input data in one of this fomats:- A list of shapely polygon geometries…
Umar Yusuf
  • 926
  • 2
  • 13
  • 30
2
votes
1 answer

OSError: could not find or load spatialindex_c.dll

I have installed the whl files for fiona, shapely, gdal, rtree and pyproj. Geopandas is also installed via pip install. But I get the following error when importing it: import geopandas Traceback (most recent call last): File…
IndigoChild
  • 842
  • 3
  • 11
  • 29
2
votes
2 answers

perpendicular geo distance from a point to a line using boost geometry

I want to get perpendicular distance from a point (t) to a line segment (p, q). The perpendicular may not intersect the line [p, q]. In that case I want to extend the line (p, q) hypothetically and then draw the perpendicular to get the distance. p,…
Neel Basu
  • 12,638
  • 12
  • 82
  • 146
2
votes
2 answers

OpenSkiMap tile not loading

Trying to add the OpenSkiMap to my map, but I can't seem to get it to load. Any suggestions? Define tile layer variables: var osm = L.tileLayer('https://{s}.tile.openstreetmap.fr/osmfr/{z}/{x}/{y}.png'); skiTrails =…
Vanludvig
  • 96
  • 1
  • 7
2
votes
1 answer

Discrete legend breaks in ggplot2

I'd like to break the legend into categories rather than having a continuous range of colours. Could someone kindly help me for the specific example I am using here? Below is my current trial with colour breaks at 40, 60 and 80. Thank you very…
Cecile
  • 527
  • 5
  • 22
2
votes
2 answers

Earth Engine getThumbURL blurry image

I'm trying to find a way to specify a lat and long and retrieve a close up image. The code below allows me to enter a lat and long but the image is very blurry. Is there a simple way to get a higher resolution image? My main issue is specifying the…
hmad
  • 33
  • 6
2
votes
2 answers

Convert Longitude, latitude to Pixel values using GDAL

I have a satellite GeoTIFF Image and a corresponding OSM file with only the highways. I want to convert the longitude latitude value in the OSM file to the pixels and want to highlight highway on the satellite image. I have tried several methods…
user1703942
  • 317
  • 3
  • 15
2
votes
1 answer

3D interactive surface plot with spatial data

I would like to create an interactive 3D surface plot of depths in a lake, ideally using the plotly or rgl libraries. I have extracted my data from a SpatialLinesDataFrame of contour lines in Gauss-Krueger/EPSG:31468 CRS, i.e. metric units. Now each…
M.Teich
  • 575
  • 5
  • 22
2
votes
0 answers

Using raster or rgdal packages to load MrSID raster

Is it possible to open a MrSID file using the raster or rgdal packages? Specifically I'm looking for something of the form raster(myraster.sid).
2
votes
1 answer

Land use/type information from here maps

Given a region in a map, is it possible to get data allowing me to say something like 40% is water 10% is forrest, 30% is city and 20% is open land... as an example. I know of other systems (Mapbox for one) which lets you get this land use…
JoeTaicoon
  • 1,383
  • 1
  • 12
  • 28
2
votes
3 answers

Replace elements in Numpy array by value and location

I am working on a program which will create contour data out of numpy arrays, and trying to avoid calls to matplotlib. I have an array of length L which contains NxN arrays of booleans. I want to convert this into an LxNxN array where, for example,…
D. Lef
  • 239
  • 3
  • 13
2
votes
1 answer

GIS: Merging MultiPolygons in Python

I have a geojson file of the counties in the British Isles, and I'm trying to replace the individual counties of London with a single merged county, then save the result as geojson. (The London counties can be identified because their TYPE_2…
duhaime
  • 25,611
  • 17
  • 169
  • 224
2
votes
1 answer

Installing GDAL with Anaconda

I just did a fresh install of Anaconda3 and am trying to add GDAL/OGR. I've tried installing with both conda install gdal and conda install -c conda-forge gdal. In both cases, the install completes successfully, but when I import gdal, it raises the…
triphook
  • 2,915
  • 3
  • 25
  • 34
2
votes
1 answer

Snapping a draggable point to a vector tile set line in Mapbox GL JS using turf.js

I'm trying to snap a draggable point to a vector tileset that consists of lines, but I'm not sure if it's possible using a Mapbox vector tileset. It's essentially the equivalent of this point snapping example that uses turf.js…
Mkov
  • 199
  • 1
  • 1
  • 12
2
votes
1 answer

GeoDjango. Admin not showing map location for point = models.PointField() in created instances

I am trying to create member instances in a very simple way to practice GeoDjango. I am able to create the instances but they are not showing correctly in the admin map (See images Below) below are my models.py from django.contrib.gis.db import…
Samir Tendulkar
  • 1,151
  • 3
  • 19
  • 47