precision mediump float; uniform float time; uniform vec2 resolution; float PI = 3.1415926535; mat2 rot(float r){ return mat2(cos(r),sin(r),-sin(r),cos(r)); } vec2 pmod(vec2 p,float n){ float np = acos(-1.)*2.0/n; float r = atan(p.x,p.y)-0.5*np; r = mod(r,np)-0.5*np; return length(p)*vec2(cos(r),sin(r)); } float random (in vec2 st) { return fract(sin(dot(st.xy, vec2(12.9898,78.233))) * 43758.5453123); } float noise (in vec2 st) { vec2 i = floor(st); vec2 f = fract(st); float a = random(i); float b = random(i + vec2(1.0, 0.0)); float c = random(i + vec2(0.0, 1.0)); float d = random(i + vec2(1.0, 1.0)); vec2 u = f*f*(3.0-2.0*f); return mix(a, b, u.x) + (c - a)* u.y * (1.0 - u.x) + (d - b) * u.x * u.y; } #define OCTAVES 6 float fbm (in vec2 st) { float value = 0.0; float amplitude = .5; float frequency = 0.; for (int i = 0; i < OCTAVES; i++) { value += amplitude * noise(st); st *= 2.; amplitude *= .5; } return value; } vec4 ground(vec3 p){ float d = p.y+0.7*fbm(p.xz)+16.0*noise(p.xz*0.08-vec2(0.0,0.1))-0.15; return vec4(0.,0.,0.,d); } // from:https://qiita.com/aa_debdeb/items/bffe5b7a33f5bf65d25b vec4 qmul(vec4 a, vec4 b) { return vec4( a.x * b.x - a.y * b.y - a.z * b.z - a.w * b.w, a.x * b.y + a.y * b.x - a.z * b.w + a.w * b.z, a.x * b.z + a.y * b.w + a.z * b.x - a.w * b.y, a.x * b.w - a.y * b.z + a.z * b.y + a.w * b.x ); } #define ITERATIONS 8 vec2 Julia(vec4 p, vec4 c) { vec4 z = p; vec4 dz = vec4(1.0, 0.0, 0.0, 0.0); vec4 pz, pdz; float r = 0.0, dr = 1.0; for (int i = 0; i < ITERATIONS; i++) { pz = z; z = qmul(pz, pz) + c; pdz = dz; dz = 2.0 * qmul(pz, pdz); r = length(z); dr = length(dz); if (r > 4.0) break; } return vec2(0.5 * log(r) * r / dr,r); } float cube(vec3 p,vec3 s){ vec3 q = abs(p); vec3 m = max(s-q,0.0); return length(abs(max(q-s,0.0)))-min(min(m.z,m.y),m.x); } vec4 tower(vec3 p){ float s = 0.11; p -= vec3(-50.7,-6.,-3.); p.xy *= rot(1.57); float ci = 0.0; float d = cube(p,vec3(15.,7.,7.)); vec2 dj =Julia(vec4(p*s,0.0),vec4(-0.9,0.2,0.2,0.0))/s; d = dj.x; ci = dj.y; vec3 mat = vec3(clamp(log(ci*0.5),0.0,1.0)); mat = vec3(1.0); float ks = 1.2; p = mod(p,ks)-0.5*ks; float cd = cube(p,vec3(0.4)); d = max(d,-cd); return vec4(mat,d); } float smoothMin(float d1, float d2, float k){ float h = exp(-k * d1) + exp(-k * d2); return -log(h) / k; } vec4 airplay(vec3 p){ p.x +=8.; p.z -= 2.2; p.y += 0.7; p.xz *= rot(0.9); float s = 1.3; p *= s; vec3 pp = p; p.y += 0.35; float d1 = cube(p+vec3(0.0,1.,0.0),vec3(3.7,1.6,2.9)); for(int i = 0;i<4;i++){ p = abs(p)-vec3(0.6,-0.13,0.4); p.xz *= rot(1.5); p.xy *= rot(1.4); } float d = cube(p,vec3(0.13,0.9,1.9)); d = max(d,-d1); p = pp; float dk1 = cube(p+vec3(0.0,1.4,0.0),vec3(3.7,1.6,2.9)); p.y -= 0.5; p.x = abs(p.x)-0.4; p.xy = pmod(p.xy,6.); for(int i = 0;i<4;i++){ p = abs(p)-0.4*vec3(0.2,0.,0.2); p.xz *= rot(0.1); p.xy *= rot(1.1); } float d2 = cube(p,vec3(0.6,0.4,1.9)); d2 = max(d2,-dk1); d = min(d,d2); p = pp; p.xz = abs(p.xz)-0.9; p.yz *= rot(1.); p.xz *= rot(-0.4); float d3 = cube(p,vec3(0.05,0.05,1.5)); d = min(d,d3); p = pp; p.x = abs(p.x)-0.8; p.z -= 1.3; p.y -= 0.9; p.xy *= rot(0.48); p.xz *= rot(0.3); p.yz *= rot(0.2); float d4 = cube(p,vec3(0.04,0.18,0.7)); d = max(d,-d4); vec3 mat = vec3(1.0); mat = mix(0.3*vec3(0.5,0.9,0.9),mat,step(0.01,d4)); p = pp; p.z += 2.; p.y -= 1.2; float d6 = cube(p+vec3(0.,0.,1.),vec3(0.2,2.,1.)); p.yz *= rot(1.57*0.5); float d5 = cube(p,vec3(0.02,0.9,0.9)); d = min(d,max(max(d5,-d6),-dk1)); return vec4(mat,d/s); } vec4 map(vec3 p){ vec4 g = ground(p); vec4 ap = airplay(p); vec4 tw = tower(p); return min(min(g,tw),ap); } vec4 volume(vec3 p){ p -= vec3(-43.7,3.2,-9.5); p.yz *= rot(0.5); p.xy *= rot(0.2); p.xz *= rot(2.3); float sx = 7.9; float sc = 2.8*pow((1.0-(p.x+sx)/(2.0*sx)),2.0); p.y += sin(time+p.x); float d1 = cube(p,vec3(sx,1.*sc,1.*sc))-2.0*fbm(p.xy*1.2-time*vec2(-2.,1.)); return vec4(-0.6*vec3(1.0),d1); } vec4 volumeFm(vec3 p){ p -= vec3(-43.7,3.2,-9.5); p.yz *= rot(0.5); p.xy *= rot(0.2); p.xz *= rot(2.3); float sx = 4.9; float sc = 1.4*(1.0-(p.x+sx)/(2.0*sx)); p.y += sin(time+p.x); float d1 = cube(p,vec3(sx,1.*sc,1.*sc))-2.0*fbm(p.xy*1.1-time*vec2(-2.,1.)); return vec4(clamp((1.0-sc*0.9),0.0,1.0)*vec3(1.0,0.2,0.0),d1); } vec3 getNormal(vec3 p) { const float h = 0.001; const vec2 k = vec2(1,-1); return normalize( k.xyy*map( p + k.xyy*h ).w + k.yyx*map( p + k.yyx*h ).w + k.yxy*map( p + k.yxy*h ).w+ k.xxx*map( p + k.xxx*h ).w ); } float softray(vec3 ro, vec3 rd , float hn) { float t = 0.01; float jt = 0.0; float res = 1.; for (int i = 0; i < 20; ++i) { jt = map(ro+rd*t).w; res = min(res,jt*hn/t); t = t+ clamp(0.02,10.,jt); } return clamp(res,0.0,1.0); } float ao(vec3 p,vec3 n,float len,float power){ float oss =0.0; for(int i =0;i<3;i++){ float d = map(p+n*len/3.0*float(i+1)).w; oss += (len-d)*power; power *=0.5; } return clamp(1.-oss,0.0,1.0); } vec3 LD = normalize(vec3(1.,1.,1.)); vec3 groundColor(vec3 p,vec3 n,vec3 col,vec4 sla){ float hl = max(dot(n, LD),0.0); hl = hl * 0.5 + 0.5; hl = hl * hl; vec3 gcol = vec3(1.0); gcol = mix(vec3(103./255.,66./255.,41./255.),gcol,smoothstep(0.8,0.9,n.y)); float sh = clamp((sla.x+0.1),0.0,1.0); return vec3(hl)*gcol*sh*sla.y; } vec3 skycolor(vec3 p){ vec3 col = vec3(0.1,0.1,0.8); float f = fbm(p.xz*0.1+21.3); col = mix(col,vec3(1.0),smoothstep(0.4,0.7,f)); return col; } vec3 towercolor(vec3 p,vec3 n,vec3 col,vec4 sla){ float hl = max(dot(n, LD),0.0); hl = hl * 0.5 + 0.5; hl = hl * hl; float sh = clamp((sla.x+0.1),0.0,1.0); return vec3(hl)*col*sh*sla.y; } float cube2D(vec2 p,vec2 s){ vec2 q = abs(p); vec2 m = max(s-q,0.0); return length(abs(max(q-s,0.0)))-min(m.x,m.y); } vec3 airplaycolor(vec3 p,vec3 n,vec3 col,vec3 rd,vec4 sla){ p.x +=8.; p.z -= 2.2; p.y += 0.7; p.xz *= rot(0.9); for(int i = 0; i<4;i++){ p = abs(p)-0.5; p.xy *= rot(1.2); p.yz *= rot(1.0); } float ds = abs(cube2D(p.xz,vec2(0.5,0.2))); vec3 em = vec3(0.0); if(ds<0.02&&col.r>0.9){ col = vec3(0.5,0.5,0.5); em = 1.0*vec3(0.9,0.1,0.1)*pow(abs(sin(time)),2.0); } float NdotL = max(dot(n, LD),0.0); vec3 R = normalize( - LD + 2.0 *n * NdotL ); float spec = pow(max(0., dot(R, -rd)), 10.0); float sh = clamp((sla.x+0.1),0.0,1.0); return (NdotL*col+spec)*sh*sla.y+em; } //https://www.shadertoy.com/view/lsKSWR float vig(vec2 uv) { uv *= 1.0 - uv.yx; float vig = uv.x*uv.y * (12.0+7.*noise(vec2(time*10.))); vig = pow(vig, 0.45); return vig; } void main() { vec2 uv = gl_FragCoord.xy / resolution.xy; vec2 uv2 = uv; uv = (uv-0.5)*2.; uv.y *= resolution.y/resolution.x; vec2 p = uv; float rkt = 3.7+0.1*mod(time*0.2,24.0*0.2); vec3 rk = vec3(-50.7+30.*sin(rkt),8.3,-3.+30.*cos(rkt)); vec3 mov = (mod(time,24.0)-8.0)*0.5*vec3(0.4,0.0,0.4); vec3 ro = mix(vec3(1.8,1.3,0.)+mov,rk,step(12.,mod(time,24.))); vec3 ta = vec3(1.0,1.2,1.3); ta += mov*step(mod(time,24.),12.); ta += vec3(-50.7,-1.,-3.); vec3 fo = normalize(ta-ro); vec3 le = cross(fo,vec3(0.,1.,0.)); vec3 up = cross(le,fo); float fov= 2.2; vec3 rd = normalize(vec3(fo*fov+up*p.y+le*uv.x)); float t = 0.002; float vd; float vdF; float d = 0.0; vec4 red; vec4 revd; vec4 revdF; vec3 accol = vec3(0.0); vec3 col = vec3(0.0); float BSC = 0.4; float ACOL = 0.4; float bac = BSC; for(int i = 0;i<159;i++){ vec4 red = map(ro+rd*t); d = red.w; revd = volume(ro+rd*t); vd = max(revd.w,bac); revdF = volumeFm(ro+rd*t); vdF = max(revdF.w,bac); float jst = step(0.2,d); accol += ACOL*revd.xyz*exp(-vd)*jst; accol += ACOL*revdF.xyz*exp(-vd); d = min(min(d,vd),vdF); t +=0.7*d; } vec3 spos = ro+rd*t; vec3 normal = getNormal(spos); vec3 ref = reflect(rd,normal); float sh = softray(spos,LD,3.3); float ao = ao(spos,normal,0.25,1.); vec4 sla = vec4(sh,ao,0.0,0.0); if(ground(spos).w<0.01){ col = groundColor(spos,normal,red.xyz,sla); }else if(tower(spos).w<0.01){ col = towercolor(spos,normal,tower(spos).xyz,sla); } else if(airplay(spos).w<0.01){ col = airplaycolor(spos,normal,airplay(spos).xyz,rd,sla); } vec3 fogcol = vec3(0.5,0.5,0.5); float far = 350.; float near =0.0; float fogk = clamp((far-t-300.*fbm(time+spos.xz*0.1+0.01*spos.y))/(far-near),0.0,1.0); col = mix(fogcol,col,fogk); col += 0.1*accol; col *= vig(uv2); gl_FragColor = vec4(col,1.0); }