Airy GLSL / Fil | Observable

Fil

Vocateur.

Workspace

Published

Edited

Jul 7, 2021

Fork of Azimuthal Equidistant GLSL

•

1 fork

GLSL projections

Rectangular Polyconic GLSL Satellite GLSL Van Der Grinten IV GLSL Van Der Grinten III GLSL Van Der Grinten II GLSL Winkel Tripel GLSL Lee Modified Stereographic GLSL Miller Oblated Stereographic GLSL Modified Stereographic GS48 GLSL Modified Stereographic GS50 GLSL Foucaut GLSL Lagrange GLSL Kavrayskiy VII GLSL Eisenlohr GLSL Eckert VI GLSL Eckert V GLSL Eckert IV GLSL Eckert III GLSL Eckert II GLSL Eckert I GLSL Bertin1953 GLSL Wagner VI GLSL Phytoplankton Briesemeister GLSL Rectilinear GLSL Azimuthal Equal-Area GLSL WebMercator to globe Cordiform GLSL Transverse Mercator GLSL Mercator GLSL Mollweide GLSL Atlantis GLSL Wiechel GLSL Wagner VII GLSL Times GLSL Wagner IV GLSL Equal Earth GLSL Conic Equal-Area GLSL Albers GLSL Equirectangular GLSL Conic Conformal GLSL Conic Equidistant GLSL Van Der Grinten GLSL Littrow GLSL Azimuthal Equidistant GLSL Armadillo GLSL Larrivée GLSL Cylindrical Equal-Area GLSL Aitoff GLSL

Airy GLSL

August GLSL Berghaus GLSL Baker Transverse GLSL Baker GLSL Boggs GLSL Bonne GLSL Hammer GLSL Bromley GLSL Bottomley GLSL Collignon GLSL Fahey GLSL Transverse Fahey GLSL Loximuthal GLSL Cylindrical Stereographic GLSL Cylindrical Equal-Area GLSL Polyconic GLSL Patterson GLSL Miller GLSL Robinson GLSL

Also listed in…

Projections

projection = d3

.geoAiry()

.radius(radius)

.clipAngle(clipAngle)

.rotate([0, -90])

.fitExtent([[2, 2], [width - 2, height - 2]], { type: "Sphere" })

import { image } from "@fil/armadillo-glsl"

glproj = {

const raw = d3.geoAiryRaw(radius * radians);

return `

float beta = ${radius * radians},

tanBeta_2 = tan(beta / 2.),

b = 2. * log(cos(beta / 2.)) / (tanBeta_2 * tanBeta_2);

float rho = sqrt(x * x + y * y);

float z = -beta / 2., z_2, cosz_2, sinz_2, tanz_2, lnsecz_2, delta = 1.;

for (int i = 0; i < 25; i++) {

if (abs(delta) > 1e-6) {

z_2 = z / 2.,

cosz_2 = cos(z_2),

sinz_2 = sin(z_2),

tanz_2 = sinz_2 / cosz_2,

lnsecz_2 = -log(abs(cosz_2));

z -= delta = (2. / tanz_2 * lnsecz_2 - b * tanz_2 - rho) / (-lnsecz_2 / (sinz_2 * sinz_2) + 1. - b / (2. * cosz_2 * cosz_2)) * (cosz_2 < 0.0 ? 0.7 : 1.);

}

float sinz = sin(z);

lambda = atan(x * sinz, rho * cos(z));

phi = asin(y * sinz / rho);

if (rho > ${raw(0, projection.clipAngle() * radians)[1]}) transparent = true;

}

height = width

import { reglCanvas, createDrawCommand, regl } with {

width,

height,

projection,

glproj,

image

} from "@fil/phytoplankton"

d3 = require("d3@7", "d3-geo-projection@3")

import { drag } from "@d3/versor-dragging"

import { slider } from "@jashkenas/inputs"

import { cos, radians } from "@fil/math"

land = d3.json(

"https://unpkg.com/visionscarto-world-atlas@0.0.6/world/110m_land.geojson"

)

One platform to build and deploy the best data apps

Experiment and prototype by building visualizations in live JavaScript notebooks. Collaborate with your team and decide which concepts to build out.

Use Observable Framework to build data apps locally. Use data loaders to build in any language or library, including Python, SQL, and R.

Seamlessly deploy to Observable. Test before you ship, use automatic deploy-on-commit, and ensure your projects are always up-to-date.

Learn more