Skip to content

Instantly share code, notes, and snippets.

View hobu's full-sized avatar

Howard Butler hobu

434 followers · 164 following
View GitHub Profile
@hobu
hobu / gist:62326dd6a98be59ead75316af1078faf
Created October 9, 2019 20:33
COMPD_CS["GCS_NZGD_2000", PROJCS["NZGD2000 / New Zealand Transverse Mercator 2000", GEOGCS["NZGD2000", DATUM["New_Zealand_Geodetic_Datum_2000", SPHEROID["GRS 1980",6378137,298.257222101004, AUTHORITY["EPSG","7019"]], AUTHORITY["EPSG","6167"]], PRIMEM["Greenwich",0], UNIT["degree",0.0174532925199433, AUTHORITY["EPSG","9122"]], AUTHORITY["EPSG","4167"]], PROJECTION["Transverse_Mercator"], PARAMETER["latitude_of_origin",0], PARAMETER["central_meridian",173], PARAMETER["scale_factor",0.9996], PARAMETER["false_easting",1600000], PARAMETER["false_northing",10000000], UNIT["metre",1, AUTHORITY["EPSG","9001"]], AXIS["Easting",EAST], AXIS["Northing",NORTH], AUTHORITY["EPSG","2193"]], VERT_CS["unknown", VERT_DATUM["unknown",2005], UNIT["metre",1, AUTHORITY["EPSG","9001"]], AXIS["Up",UP]]]"
@hobu
hobu / driver
Created October 8, 2019 19:42
GDALDriverManager* hMgr = GetGDALDriverManager();
std::unique_ptr<GDALDriver> hDriver = std::unique_ptr<GDALDriver>(hMgr->GetDriverByName(aoDrivers[0].c_str()));
if (!hDriver)
{
CPLError(CE_Failure, CPLE_AppDefined,
"Cannot get driver for %s", aoDrivers[0].c_str());
return CE_Failure;
}
/* create output raster */
// GDALDatasetH hDstDS = nullptr;
@hobu
hobu / gist:2fba34643999f1add9e93ca895f209d2
Created September 30, 2019 19:15
unsigned int openFlags = GA_ReadOnly | GDAL_OF_VECTOR | GDAL_OF_VERBOSE_ERROR;
GDALDataset* ds = (GDALDataset*) GDALOpenEx( dsString.c_str(),
openFlags,
NULL, NULL, NULL );
@hobu
hobu / gist:e12b6adc2186b5598ce850f3237af1b3
Created September 27, 2019 15:13
void*
GEOSGeom_getUserData_r(GEOSContextHandle_t extHandle, const Geometry* g)
{
assert(0 != g);
if(nullptr == extHandle) {
return 0;
}
GEOSContextHandleInternal_t* handle = 0;
@hobu
hobu / gist:cf9977c55dfdedc2951be8ee13dab423
Created September 20, 2019 13:17
"AUTOSRS_V_NAME","COMPOUND_WKT"
"NAVD88 - Geoid03 (Feet)","COMPD_CS[""unknown"",PROJCS[""NAD83 / Pennsylvania South"",GEOGCS[""NAD83"",DATUM[""North_American_Datum_1983"",SPHEROID[""GRS 1980"",6378137,298.2572221010042,AUTHORITY[""EPSG"",""7019""]],AUTHORITY[""EPSG"",""6269""]],PRIMEM[""Greenwich"",0],UNIT[""degree"",0.0174532925199433],AUTHORITY[""EPSG"",""4269""]],PROJECTION[""Lambert_Conformal_Conic_2SP""],PARAMETER[""standard_parallel_1"",39.93333333333333],PARAMETER[""standard_parallel_2"",40.96666666666667],PARAMETER[""latitude_of_origin"",39.33333333333334],PARAMETER[""central_meridian"",-77.75],PARAMETER[""false_easting"",1968500],PARAMETER[""false_northing"",0],UNIT[""US survey foot"",0.3048006096012192,AUTHORITY[""EPSG"",""9003""]],AUTHORITY[""EPSG"",""32129""]],VERT_CS[""NAVD88 - Geoid03 (Feet)"",VERT_DATUM[""unknown"",2005],UNIT[""US survey foot"",0.3048006096012192,AUTHORITY[""EPSG"",""9003""]],AXIS[""Up"",UP]]]"
"NAVD88 - Geoid03 (Meters)","COMPD_CS[""unknown"",PROJCS[""WGS 84 / UTM zone 11N"",G
@hobu
hobu / gist:464735dd74144b1e097a74f867f96e48
Created August 29, 2019 21:50
In file included from ../tutorial.cpp:1:
In file included from /Users/hobu/miniconda3/envs/pdal-python-test/include/pdal/PointView.hpp:37:
In file included from /Users/hobu/miniconda3/envs/pdal-python-test/include/pdal/DimDetail.hpp:37:
In file included from /Users/hobu/miniconda3/envs/pdal-python-test/include/pdal/Dimension.hpp:9:
In file included from /Users/hobu/miniconda3/envs/pdal-python-test/include/pdal/DimUtil.hpp:40:
/Users/hobu/miniconda3/envs/pdal-python-test/include/pdal/util/Utils.hpp:572:21: error: expected expression
auto pred = [tChar](char c){ return(c == tChar); };
^
/Users/hobu/miniconda3/envs/pdal-python-test/include/pdal/util/Utils.hpp:588:21: error: expected expression
auto pred = [tChar](char c){ return(c == tChar); };
@hobu
hobu / gist:919c7b36046a2f256325f7afe0eab88d
Created August 22, 2019 21:37
Package libcxx conflicts for:
numpy=1.14 -> libcxx[version='>=4.0.1']
zstd=1.4.0 -> libcxx[version='>=4.0.1']
nitro -> libcxx[version='>=4.0.1']
geos=3.7.2 -> libcxx[version='>=4.0.1']
libgdal=3.0 -> libcxx[version='>=4.0.1']
sqlite=3 -> libcxx[version='>=4.0.1']
tiledb=1.6.0 -> libcxx[version='>=4.0.1']
laszip -> libcxx[version='>=4.0.1']
python=3.6 -> libcxx[version='>=4.0.1']
@hobu
hobu / gist:a0d709976c32c4ad95bcf93f7b12c1ef
Created August 22, 2019 21:16
Package python conflicts for:
libgdal=3.0 -> libkml[version='>=1.3.0,<1.4.0a0'] -> boost==1.66.0 -> numpy[version='>=1.8'] -> python[version='>=2.7,<2.8.0a0|>=3.5,<3.6.0a0|>=3.6,<3.7.0a0|>=3.7,<3.8.0a0']
geotiff=1.5.1 -> libtiff[version='>=4.0.9,<5.0a0'] -> zstd[version='>=1.4.0,<1.4.1.0a0'] -> lz4 -> python[version='>=2.7,<2.8.0a0|>=3.5,<3.6.0a0|>=3.6,<3.7.0a0|>=3.7,<3.8.0a0']
postgresql -> python[version='2.7.*|3.4.*|3.5.*|3.6.*|>=2.7,<2.8.0a0|>=3.5,<3.6.0a0|>=3.6,<3.7.0a0']
geotiff=1.5.1 -> libtiff[version='>=4.0.9,<5.0a0'] -> zstd[version='>=1.4.0,<1.4.1.0a0'] -> lz4 -> setuptools -> python[version='2.7.*|3.5.*|3.6.*']
libkml=1.3 -> boost==1.66.0 -> numpy[version='>=1.8'] -> python[version='>=2.7,<2.8.0a0|>=3.5,<3.6.0a0|>=3.6,<3.7.0a0|>=3.7,<3.8.0a0']
libkml=1.3 -> boost==1.66.0 -> python[version='2.7.*|3.5.*|3.6.*']
libgdal=3.0 -> libkml[version='>=1.3.0,<1.4.0a0'] -> boost==1.66.0 -> python[version='2.7.*|3.5.*|3.6.*']
numpy=1.14 -> python[version='2.7.*|3.4.*|3.5.*|3.6.*|>=2.7,<2.8.0a0|>=3.5,<3.6.0a0|>
@hobu
hobu / gist:7746a89ad7abf8fa82785717b784165a
Created July 24, 2019 22:28
UnsatisfiableError: The following specifications were found to be incompatible with each other:
- conda-forge/osx-64::boost-cpp==1.70.0=h75728bb_2 -> libboost[version='<0']
- conda-forge/osx-64::cfitsio==3.470=h389770f_1 -> libcurl[version='>=7.64.1,<8.0a0'] -> krb5[version='>=1.16.3,<1.17.0a0'] -> libedit[version='>=3.1.20170329,<4.0a0']
- conda-forge/osx-64::curl==7.65.3=h22ea746_0 -> krb5[version='>=1.16.3,<1.17.0a0'] -> libedit[version='>=3.1.20170329,<4.0a0']
- conda-forge/osx-64::krb5==1.16.3=hcfa6398_1001 -> libedit[version='>=3.1.20170329,<4.0a0']
- conda-forge/osx-64::libcurl==7.65.3=h16faf7d_0 -> krb5[version='>=1.16.3,<1.17.0a0'] -> libedit[version='>=3.1.20170329,<4.0a0']
- conda-forge/osx-64::libedit==3.1.20170329=hcfe32e1_1001
- conda-forge/osx-64::libkml==1.3.0=hed7d534_1010 -> boost-cpp[version='>=1.70.0,<1.70.1.0a0'] -> libboost[version='<0']
- conda-forge/osx-64::libpq==11.4=h756f0eb_3 -> krb5[version='>=1.16.3,<1.17.0a0'] -> libedit[version='>=3.1.20170329,<4.0a0']
@hobu
hobu / gist:6e52149e108c37ae6d6da0203b2320e7
Created June 10, 2019 14:22
aspect-slope.sh
# adaptation of https://kingsgeocomputation.org/2016/03/16/aspect-slope-maps-in-qgis/ in
# gdal_calc
cat > color-ramp.txt <<EOL
# QGIS Generated Color Map Export File
10,181,181,181,255,Custom color map entry
12,159,192,133,255,Custom color map entry
14,157,221,94,255,Custom color map entry
16,154,251,12,255,Custom color map entry
20,181,181,181,255,Custom color map entry