Ryan Cooper maptastik

It is sort of tricky to run spatial queries against tables in a GeoPackage anywhere, but especially outside of QGIS. This is because we're reliant on Spatialite to carry out those queries and Spatialite is just not well supported. A real shame!

Luckily there are a few resources for helping you get up and going with installing Spatialite for working with GeoPackages in a SQL client like DB Browser.

There is a page in the DB Browser Wik on GitHub called " SpatiaLite on Windows" and it is dedicated to installing Spatialite for use with DB Browser. Note that there are special instructions for a Windows 10 issue.

Bryan McBride at Spatial Networks wrote up a nice piece on working with GeoPackages in DB Browser called "Working with Geospatial Data: An Introduction". It's got a nice walk-thru of how to use DB Browser with GeoPackage as well as some general installation instructions.

Problem

When trying to use cURL from the OSGEO4W shell, it throws the following error:

curl: (60) SSL certificate problem: unable to get local issuer certificate More details here: https://curl.haxx.se/docs/sslcerts.html

curl failed to verify the legitimacy of the server and therefore could not establish a secure connection to it. To learn more about this situation and how to fix it, please visit the web page mentioned above.

This is an example of exporting an environment so it can be used to create a new environment. It will also show how to specify a channel that should be included in the build process

First, you'll need to activate the environment you want to use.

activate old-environment

Then export then environment to environment.yml

Fade-in layers examples:

Add features from several GPX files in a directory to a GeoPackage

This is motivated by the fact that Runkeeper exports all your GPS data to individual GPX files. These GPX files contain 5 different tables:

route_points
routes
track_points
tracks
waypoints

	<!DOCTYPE html>
	<html lang="en">
	<head>
	<meta charset="UTF-8">
	<title>Map of single-part geometry data pulled directly from AGOL</title>
	</head>
	<body>


	<!-- Load require.js. Delete this if your page already loads require.js -->

	""" Generate a list of unique values for a field in a feature class

	Parameters:
	fc (str): Input feature class
	field (str): Field to find unique values for

	Returns:
	list: unique values in a field

	"""

	import arcpy
	import numpy as np

	def randomFCSample(fc, fd='sample_fd', sample_field='OBJECTID', sample_pct=10):
	"""Return a randomish sample of a feature class as a feature layer

	Parameters:
	fc (str): Input feature class to sample
	fd (str): Name of output feature layer
	sample_field (int): Numeric field to sample values from