Python – How to Determine Ground Type from Given Coordinates

coordinatespythonsoil

I would like to know what type of soil (ground) is a given coordinate (lat&lon). For instance, if a coordinate is rural soil, a road, a building, etc. Is it possible? I've tried:

import requests
url = "https://rest.isric.org/soilgrids/v2.0/properties/query?lon=34.5&lat=1.299"
r = requests.get(url)
print(r.json())

However, the output JSON is not clear the type of soil.

Best Answer

I guess you're looking for landcover classifications...

Some good global open-access datasets are:

Related Solutions

[GIS] Given a set of coordinates, how can I calculate the minimum bounds

Hope i understand the question correctly...

We can solve the problem for longitude and latitude separately, so i will take your example with the longitudes: -76, -135 and 164.

First i would order them:

-135, -76, 164

Then i would add the left most coordinate to the right again: -135 + 360 = 225

-135, -76, 164, 225

Now we can calculate the gaps between the coordinates:

-135 (59) -76 (240) 164 (61) 225
             .......

The biggest gap (240) must be the boundary of the minimum bounding box, the part that does not belong to the box. The dotted line is the biggest part of the circle we can spare out. In our example that means, the boundary box starts with 164, includes -135 and ends with -76.

[GIS] How to get the coordinates inside the polygon using python

A point inside a polygon

canvas = qgis.utils.iface.mapCanvas()
aLayer = canvas.layer(0)  # or code to select the layer of interest
for f in aLayer.getFeatures():
    geom = f.geometry()  
    p = geom.pointOnSurface()
    print p.asPoint()

Get Inner Rings

canvas = qgis.utils.iface.mapCanvas()
aLayer = canvas.layer(0)  # or code to select the layer of interest
for f in aLayer.getFeatures():
    geom = f.geometry()
    if len(geom.asPolygon()) > 1 : # we have inner rings
      print geom.asPolygon()[1:]