-
Notifications
You must be signed in to change notification settings - Fork 26
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
Making Protest Code Snippets, Artworks and T-shirts #2
Comments
Once you have your artwork, it is easy to submit it to Redbubble for making and selling it automatically. |
some posters, calendars like this artist: http://www.redbubble.com/people/cem-?ref=artist_title_name |
A "not quite code" transcription of the forward algo is available at https://www.wikidata.org/wiki/Q29043602. I guess that can be something to start with? |
Condensed the forward ops, 37lines×56cols (CC0 as usual): (EDIT: just debugged it, so it's all 58col for now) /**
* MiniPRCoords. Public Domain.
* WGS → GCJ → BD
* (35, 105) → (34.999093, 105.003286) → (35.005398, 105.009667)
* (34, 106) → (33.998424, 106.003994) → (34.004217, 106.010579)
* https://github.com/Artoria2e5/PRCoords */
let {sqrt, abs, sin, cos, hypot, atan2, PI} = Math
function wgs_to_gcj({lat, lon}) {
/* 常数 / Constants: */
let A = 6378245 // Krasovsky 1940/SK-42: a
let _F = 1 / 298.3 // Krasovsky 1940/SK-42: f
let EE = 2*_F - _F**2 // e^2 = 2*f - f^2
/* 偏移 / Shifts: */
// noise origin: (35, 105); default shift: <300, -100>
let y = lat - 35, x = lon - 105
let spi = ((z) => sin(z*PI))
let Δ_N = -100+2*x+3*y+0.2*y*y+0.1*x*y+0.2*sqrt(abs(x))
+20/3*(2*spi(6*x)+2*spi(2*x)+2*spi(y)
+4*spi(y/3)+16*spi(y/12)+32*spi(y/30))
let Δ_E = 300+1*x+2*y+0.1*x*x+0.1*x*y+0.1*sqrt(abs(x))
+20/3*(2*spi(6*x)+2*spi(2*x)+2*spi(x)
+4*spi(x/3)+15*spi(x/12)+30*spi(x/30))
/* 经纬 / Lat-lon packing: */
let lat_r = lat * PI / 180
let common = 1 - EE * sin(lat_r)**2
// meter → °
let lat1 = (PI / 180) * (A * (1 - EE)) * common**1.5
let lon1 = (PI / 180) * (A * cos(lat_r)) / common**0.5
return {lat: lat + Δ_N/lat1, lon: lon + Δ_E/lon1}
}
let B_F = 3000/180*PI
let B_DLAT = 0.0060, B_DLON = 0.0065
function gcj_to_bd({lat, lon}) {
let r = hypot(lat, lon) + 2e-5 * sin(B_F * lat)
let θ = atan2(lat, lon) + 3e-6 * cos(B_F * lon)
return {lat: r*sin(θ) + B_DLAT, lon: r*cos(θ) + B_DLON}
}
df = (a, b) => ({lat: a.lat - b.lat, lon: a.lon - b.lon})
function gcj_to_wgs(a) {
return df(a, df(wgs_to_gcj(a), a))
}
function bd_to_gcj(a) {
let {lat, lon} = df(a, {lat: B_DLAT, lon: B_DLON})
let r = hypot(lat, lon) - 2e-5 * sin(B_F * lat)
let θ = atan2(lat, lon) - 3e-6 * cos(B_F * lon)
return {lat: r*sin(θ), lon: r*cos(θ)}
}
itery = (fwd, rev) => function(bad) {
let curr = rev(bad)
let prev = bad
let diff = {lat: Infinity, lon: Infinity}
while (abs(diff.lat) > 1e-5 || abs(diff.lon) > 1e-5) {
diff = df(fwd(curr), bad)
prev = curr
curr = df(curr, diff)
}
return curr
} I guess I will get it slightly syntax highlighted. Function names and comments should be emphasized so they look like section titles. I don't quite feel like expressing the gcj and iterative backward ones in JS. Some formalae should do. |
Fits on an A4 at unknown size too. The attachment sure looks ugly, but hey you got something to give to friends now! |
Ayy. I… uploaded this thing to… teespring. You can get a tote bag, a mug, or a sticker. https://teespring.com/miniprcoords-tote-v1 |
i still dont think the current code fits well on a t-shirt and i have no good idea how it will. the manual version minifies to ~1100 characters with variable renaming and other uglyfied things, which is around twice the length of the perl thing. so.. 11 lines of garbage or 37 lines of loose small text? |
🤔 |
I was thinking about some possible ways to eliminate the huge chunk of code on its nonlinear transformation by, e.g. extracting constants or making a lookup-table, but it seems the idea doesn't work well, not to mention the understandability of the code should be a focus. Is using a language with higher-level builtin math operations (Mathematica? Maple? Sage? Julia? Octave?) can help simplifying it? 🤔 For example, the Matlab code looks considerably shorter. https://github.com/Artoria2e5/PRCoords/blob/master/matlab/PRCoords.m Avoid coding high-level transformation in the code should help I guess, if the external function performs a well-known operation, it shouldn't be much of a problem. |
I… don't think the MATLAB code is shorter than the js. and OH FUCK THE CAIJUN ITERATION SHOULD BE A DO WHILE |
It would be helpful to make some protest code snippets, artworks, T-shirts that are easy to express and reproduce the core algorithm of the transformation of coordinates, i.e something like RSA in CrytoWar, DeCSS in DVD, or Free Speech Flag in Blu-ray. It also makes the act of activism more interesting to participate.
Perl program:
Favicon:
Idea: embed the icon in the favicon of PRCoords
DeCSS Haiku:
The text was updated successfully, but these errors were encountered: