RichardScottOZ

raster from multiple sources, disparate grid specs

This FlatGeoBuf file contains two polygon footprints, one is a northern hemisphere sea-ice concentration GeoTIFF image, the other southern (the values are scaled from 0-100% and stored with a palette, but we're just using the raw integer value here)

https://github.com/mdsumner/gti/raw/refs/heads/main/inst/extdata/tmerc_seaice.fgb

The source GeoTIFFs are in NSIDC south Polar Stereographic, and north Polar Stereographic as is standard. We use GTI format to register these as disparate data sources within a latent raster frame for rendering at a later time. We could store all the daily sea ice files this way, with an index on field 'date', so we could coalesce to a particular date (or set of dates, for examples like Sentinel scenes).

	import logging
	import subprocess as sp
	import tempfile
	import threading

	import geopandas as gpd
	import numpy as np
	import rioxarray
	import xarray as xr
	from rasterio.enums import Resampling

	#!/bin/sh

	log() {
	format="$1"; shift
	# shellcheck disable=SC2059
	>&2 printf "$format\n" "$@"
	}

	usage() {
	>&2 cat <<EOF

	# coding=utf-8
	# Copyright 2023 The HuggingFace Inc. team. All rights reserved.
	#
	# Licensed under the Apache License, Version 2.0 (the "License");
	# you may not use this file except in compliance with the License.
	# You may obtain a copy of the License at
	#
	# http://www.apache.org/licenses/LICENSE-2.0
	#
	# Unless required by applicable law or agreed to in writing, software

	# Copyright (c) 2023 Derrick W. Turk / terminus, LLC

	# Permission is hereby granted, free of charge, to any person obtaining a copy
	# of this software and associated documentation files (the "Software"), to deal
	# in the Software without restriction, including without limitation the rights
	# to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
	# copies of the Software, and to permit persons to whom the Software is
	# furnished to do so, subject to the following conditions:

	# The above copyright notice and this permission notice shall be included in all

	#!/bin/bash
	# see https://github.com/microsoft/RoadDetections

	# 1) get the data
	wget https://usaminedroads.blob.core.windows.net/road-detections/Oceania-Full.zip

	# 2) extract the records for Australia, second column of tsv only
	zgrep AUS Oceania-Full.zip \| cut -f2 > AUS.geojson

	# 3) convert to GPKG, works better in QGIS

	#!/bin/bash

	api_call() {
	curl --no-progress-meter https://api.openai.com/v1/completions \
	-H "Content-Type: application/json" \
	-H "Authorization: Bearer your-api-key" \
	-d "{
	\"model\": \"text-davinci-003\",
	\"prompt\": \"$@\",
	\"temperature\": 0,

	# main
	llama-index
	langchain