@maplat/transform 0.1.0 → 0.1.1

package/package.json

@@ -4,7 +4,7 @@
   "publishConfig": {
     "access": "public"
   },
-  "version": "0.1.0",
+  "version": "0.1.1",
   "description": "A JavaScript library that performs coordinate transformation between two plane coordinate systems using transformation definitions generated by Maplat.",
   "type": "module",
   "main": "./dist/maplat_transform.cjs",

package/src/index.ts

@@ -10,6 +10,9 @@ import type {
 } from "./geometry";
 import { unitCalc, transformArr } from "./geometry";
 import type { EdgeSet, EdgeSetLegacy } from "./edgeutils";
+export * from './geometry';
+export * from './triangulation';
+export * from './edgeutils';
 
 /**
  * 座標ペアの型定義。[ソース座標, ターゲット座標] の形式

package/dist/maplat_transform.cjs

@@ -1 +0,0 @@
-"use strict";var it=Object.defineProperty;var ot=(r,t,e)=>t in r?it(r,t,{enumerable:!0,configurable:!0,writable:!0,value:e}):r[t]=e;var A=(r,t,e)=>ot(r,typeof t!="symbol"?t+"":t,e);Object.defineProperty(exports,Symbol.toStringTag,{value:"Module"});const B=11102230246251565e-32,k=134217729,at=(3+8*B)*B;function $(r,t,e,o,f){let s,c,g,l,y=t[0],w=o[0],i=0,u=0;w>y==w>-y?(s=y,y=t[++i]):(s=w,w=o[++u]);let d=0;if(i<r&&u<e)for(w>y==w>-y?(c=y+s,g=s-(c-y),y=t[++i]):(c=w+s,g=s-(c-w),w=o[++u]),s=c,g!==0&&(f[d++]=g);i<r&&u<e;)w>y==w>-y?(c=s+y,l=c-s,g=s-(c-l)+(y-l),y=t[++i]):(c=s+w,l=c-s,g=s-(c-l)+(w-l),w=o[++u]),s=c,g!==0&&(f[d++]=g);for(;i<r;)c=s+y,l=c-s,g=s-(c-l)+(y-l),y=t[++i],s=c,g!==0&&(f[d++]=g);for(;u<e;)c=s+w,l=c-s,g=s-(c-l)+(w-l),w=o[++u],s=c,g!==0&&(f[d++]=g);return(s!==0||d===0)&&(f[d++]=s),d}function ft(r,t){let e=t[0];for(let o=1;o<r;o++)e+=t[o];return e}function D(r){return new Float64Array(r)}const ct=(3+16*B)*B,ut=(2+12*B)*B,ht=(9+64*B)*B*B,X=D(4),j=D(8),J=D(12),z=D(16),S=D(4);function lt(r,t,e,o,f,s,c){let g,l,y,w,i,u,d,v,h,a,n,b,M,x,E,_,m,T;const I=r-f,P=e-f,N=t-s,C=o-s;x=I*C,u=k*I,d=u-(u-I),v=I-d,u=k*C,h=u-(u-C),a=C-h,E=v*a-(x-d*h-v*h-d*a),_=N*P,u=k*N,d=u-(u-N),v=N-d,u=k*P,h=u-(u-P),a=P-h,m=v*a-(_-d*h-v*h-d*a),n=E-m,i=E-n,X[0]=E-(n+i)+(i-m),b=x+n,i=b-x,M=x-(b-i)+(n-i),n=M-_,i=M-n,X[1]=M-(n+i)+(i-_),T=b+n,i=T-b,X[2]=b-(T-i)+(n-i),X[3]=T;let U=ft(4,X),Y=ut*c;if(U>=Y||-U>=Y||(i=r-I,g=r-(I+i)+(i-f),i=e-P,y=e-(P+i)+(i-f),i=t-N,l=t-(N+i)+(i-s),i=o-C,w=o-(C+i)+(i-s),g===0&&l===0&&y===0&&w===0)||(Y=ht*c+at*Math.abs(U),U+=I*w+C*g-(N*y+P*l),U>=Y||-U>=Y))return U;x=g*C,u=k*g,d=u-(u-g),v=g-d,u=k*C,h=u-(u-C),a=C-h,E=v*a-(x-d*h-v*h-d*a),_=l*P,u=k*l,d=u-(u-l),v=l-d,u=k*P,h=u-(u-P),a=P-h,m=v*a-(_-d*h-v*h-d*a),n=E-m,i=E-n,S[0]=E-(n+i)+(i-m),b=x+n,i=b-x,M=x-(b-i)+(n-i),n=M-_,i=M-n,S[1]=M-(n+i)+(i-_),T=b+n,i=T-b,S[2]=b-(T-i)+(n-i),S[3]=T;const et=$(4,X,4,S,j);x=I*w,u=k*I,d=u-(u-I),v=I-d,u=k*w,h=u-(u-w),a=w-h,E=v*a-(x-d*h-v*h-d*a),_=N*y,u=k*N,d=u-(u-N),v=N-d,u=k*y,h=u-(u-y),a=y-h,m=v*a-(_-d*h-v*h-d*a),n=E-m,i=E-n,S[0]=E-(n+i)+(i-m),b=x+n,i=b-x,M=x-(b-i)+(n-i),n=M-_,i=M-n,S[1]=M-(n+i)+(i-_),T=b+n,i=T-b,S[2]=b-(T-i)+(n-i),S[3]=T;const st=$(et,j,4,S,J);x=g*w,u=k*g,d=u-(u-g),v=g-d,u=k*w,h=u-(u-w),a=w-h,E=v*a-(x-d*h-v*h-d*a),_=l*y,u=k*l,d=u-(u-l),v=l-d,u=k*y,h=u-(u-y),a=y-h,m=v*a-(_-d*h-v*h-d*a),n=E-m,i=E-n,S[0]=E-(n+i)+(i-m),b=x+n,i=b-x,M=x-(b-i)+(n-i),n=M-_,i=M-n,S[1]=M-(n+i)+(i-_),T=b+n,i=T-b,S[2]=b-(T-i)+(n-i),S[3]=T;const nt=$(st,J,4,S,z);return z[nt-1]}function gt(r,t,e,o,f,s){const c=(t-s)*(e-f),g=(r-f)*(o-s),l=c-g,y=Math.abs(c+g);return Math.abs(l)>=ct*y?l:-lt(r,t,e,o,f,s,y)}function yt(r,t){var e,o,f=0,s,c,g,l,y,w,i,u=r[0],d=r[1],v=t.length;for(e=0;e<v;e++){o=0;var h=t[e],a=h.length-1;if(w=h[0],w[0]!==h[a][0]&&w[1]!==h[a][1])throw new Error("First and last coordinates in a ring must be the same");for(c=w[0]-u,g=w[1]-d,o;o<a;o++){if(i=h[o+1],l=i[0]-u,y=i[1]-d,g===0&&y===0){if(l<=0&&c>=0||c<=0&&l>=0)return 0}else if(y>=0&&g<=0||y<=0&&g>=0){if(s=gt(c,l,g,y,0,0),s===0)return 0;(s>0&&y>0&&g<=0||s<0&&y<=0&&g>0)&&f++}w=i,g=y,c=l}}return f%2!==0}function p(r,t,e={}){const o={type:"Feature"};return(e.id===0||e.id)&&(o.id=e.id),e.bbox&&(o.bbox=e.bbox),o.properties=t||{},o.geometry=r,o}function W(r,t,e={}){if(!r)throw new Error("coordinates is required");if(!Array.isArray(r))throw new Error("coordinates must be an Array");if(r.length<2)throw new Error("coordinates must be at least 2 numbers long");if(!K(r[0])||!K(r[1]))throw new Error("coordinates must contain numbers");return p({type:"Point",coordinates:r},t,e)}function wt(r,t,e={}){for(const f of r){if(f.length<4)throw new Error("Each LinearRing of a Polygon must have 4 or more Positions.");if(f[f.length-1].length!==f[0].length)throw new Error("First and last Position are not equivalent.");for(let s=0;s<f[f.length-1].length;s++)if(f[f.length-1][s]!==f[0][s])throw new Error("First and last Position are not equivalent.")}return p({type:"Polygon",coordinates:r},t,e)}function F(r,t={}){const e={type:"FeatureCollection"};return t.id&&(e.id=t.id),t.bbox&&(e.bbox=t.bbox),e.features=r,e}function K(r){return!isNaN(r)&&r!==null&&!Array.isArray(r)}function dt(r){if(!r)throw new Error("coord is required");if(!Array.isArray(r)){if(r.type==="Feature"&&r.geometry!==null&&r.geometry.type==="Point")return[...r.geometry.coordinates];if(r.type==="Point")return[...r.coordinates]}if(Array.isArray(r)&&r.length>=2&&!Array.isArray(r[0])&&!Array.isArray(r[1]))return[...r];throw new Error("coord must be GeoJSON Point or an Array of numbers")}function Q(r){if(Array.isArray(r))return r;if(r.type==="Feature"){if(r.geometry!==null)return r.geometry.coordinates}else if(r.coordinates)return r.coordinates;throw new Error("coords must be GeoJSON Feature, Geometry Object or an Array")}function bt(r){return r.type==="Feature"?r.geometry:r}function _t(r,t,e={}){if(!r)throw new Error("point is required");if(!t)throw new Error("polygon is required");const o=dt(r),f=bt(t),s=f.type,c=t.bbox;let g=f.coordinates;if(c&&vt(o,c)===!1)return!1;s==="Polygon"&&(g=[g]);let l=!1;for(var y=0;y<g.length;++y){const w=yt(o,g[y]);if(w===0)return!e.ignoreBoundary;w&&(l=!0)}return l}function vt(r,t){return t[0]<=r[0]&&t[1]<=r[1]&&t[2]>=r[0]&&t[3]>=r[1]}var q=_t;function At(r){const t=[0,1,2,0].map(o=>r[o][0][0]),e={a:{geom:r[0][0][1],index:r[0][1]},b:{geom:r[1][0][1],index:r[1][1]},c:{geom:r[2][0][1],index:r[2][1]}};return wt([t],e)}function V(r,t,e,o,f,s=!1,c){const g=r.map(l=>{(!c||c<2.00703)&&(l=tt(l));const y=isFinite(l)?t[l]:l==="c"?o:l==="b0"?f[0]:l==="b1"?f[1]:l==="b2"?f[2]:l==="b3"?f[3]:function(){const w=l.match(/e(\d+)/);if(w){const i=parseInt(w[1]);return e[i]}throw"Bad index value for indexesToTri"}();return s?[[y[1],y[0]],l]:[[y[0],y[1]],l]});return At(g)}function tt(r){return typeof r=="number"?r:r.replace(/^(c|e|b)(?:ent|dgeNode|box)(\d+)?$/,"$1$2")}function Mt(r,t){return Array.isArray(r[0])?r:r.map(e=>[e.illstNodes,e.mercNodes,e.startEnd])}function H(r,t){for(let e=0;e<t.features.length;e++)if(q(r,t.features[e]))return t.features[e]}function rt(r,t,e){const o=t.geometry.coordinates[0][0],f=t.geometry.coordinates[0][1],s=t.geometry.coordinates[0][2],c=r.geometry.coordinates,g=t.properties.a.geom,l=t.properties.b.geom,y=t.properties.c.geom,w=[f[0]-o[0],f[1]-o[1]],i=[s[0]-o[0],s[1]-o[1]],u=[c[0]-o[0],c[1]-o[1]],d=[l[0]-g[0],l[1]-g[1]],v=[y[0]-g[0],y[1]-g[1]];let h=(i[1]*u[0]-i[0]*u[1])/(w[0]*i[1]-w[1]*i[0]),a=(w[0]*u[1]-w[1]*u[0])/(w[0]*i[1]-w[1]*i[0]);if(e){const n=e[t.properties.a.index],b=e[t.properties.b.index],M=e[t.properties.c.index];let x;if(h<0||a<0||1-h-a<0){const E=h/(h+a),_=a/(h+a);x=h/b/(E/b+_/M),a=a/M/(E/b+_/M)}else x=h/b/(h/b+a/M+(1-h-a)/n),a=a/M/(h/b+a/M+(1-h-a)/n);h=x}return[h*d[0]+a*v[0]+g[0],h*d[1]+a*v[1]+g[1]]}function xt(r,t,e,o){const f=r.geometry.coordinates,s=e.geometry.coordinates,c=Math.atan2(f[0]-s[0],f[1]-s[1]),g=Et(c,t[0]),l=t[1][g];return rt(r,l.features[0],o)}function mt(r,t,e,o,f,s,c,g){let l;if(c&&(l=H(r,F([c]))),!l){if(e){const y=r.geometry.coordinates,w=e.gridNum,i=e.xOrigin,u=e.yOrigin,d=e.xUnit,v=e.yUnit,h=e.gridCache,a=R(y[0],i,d,w),n=R(y[1],u,v,w),b=h[a]?h[a][n]?h[a][n]:[]:[];t=F(b.map(M=>t.features[M]))}l=H(r,t)}return g&&g(l),l?rt(r,l,s):xt(r,o,f,s)}function R(r,t,e,o){let f=Math.floor((r-t)/e);return f>=o&&(f=o-1),f}function Et(r,t){let e=Z(r-t[0]),o=Math.PI*2,f;for(let s=0;s<t.length;s++){const c=(s+1)%t.length,g=Z(r-t[c]),l=Math.min(Math.abs(e),Math.abs(g));e*g<=0&&l<o&&(o=l,f=s),e=g}return f}function Z(r,t=!1){const e=t?function(o){return!(o>=0&&o<Math.PI*2)}:function(o){return!(o>-1*Math.PI&&o<=Math.PI)};for(;e(r);)r=r+2*Math.PI*(r>0?-1:1);return r}const L=2.00703,O=class O{constructor(){A(this,"points",[]);A(this,"pointsWeightBuffer");A(this,"strict_status");A(this,"vertices_params");A(this,"centroid");A(this,"edgeNodes");A(this,"edges");A(this,"tins");A(this,"kinks");A(this,"yaxisMode",O.YAXIS_INVERT);A(this,"strictMode",O.MODE_AUTO);A(this,"vertexMode",O.VERTEX_PLAIN);A(this,"bounds");A(this,"boundsPolygon");A(this,"wh");A(this,"xy");A(this,"indexedTins");A(this,"stateFull",!1);A(this,"stateTriangle");A(this,"stateBackward")}setCompiled(t){if(t.version||!t.tins&&t.points&&t.tins_points){this.points=t.points,this.pointsWeightBuffer=!t.version||t.version<2.00703?["forw","bakw"].reduce((o,f)=>{const s=t.weight_buffer[f];return s&&(o[f]=Object.keys(s).reduce((c,g)=>{const l=tt(g);return c[l]=s[g],c},{})),o},{}):t.weight_buffer,t.strict_status?this.strict_status=t.strict_status:t.kinks_points?this.strict_status=O.STATUS_ERROR:t.tins_points.length==2?this.strict_status=O.STATUS_LOOSE:this.strict_status=O.STATUS_STRICT,this.vertices_params={forw:[t.vertices_params[0]],bakw:[t.vertices_params[1]]},this.vertices_params.forw[1]=[0,1,2,3].map(o=>{const f=(o+1)%4,s=V(["c",`b${o}`,`b${f}`],t.points,t.edgeNodes||[],t.centroid_point,t.vertices_points,!1,L);return F([s])}),this.vertices_params.bakw[1]=[0,1,2,3].map(o=>{const f=(o+1)%4,s=V(["c",`b${o}`,`b${f}`],t.points,t.edgeNodes||[],t.centroid_point,t.vertices_points,!0,L);return F([s])}),this.centroid={forw:W(t.centroid_point[0],{target:{geom:t.centroid_point[1],index:"c"}}),bakw:W(t.centroid_point[1],{target:{geom:t.centroid_point[0],index:"c"}})},this.edges=Mt(t.edges||[]),this.edgeNodes=t.edgeNodes||[];const e=t.tins_points.length==1?0:1;this.tins={forw:F(t.tins_points[0].map(o=>V(o,t.points,t.edgeNodes||[],t.centroid_point,t.vertices_points,!1,t.version))),bakw:F(t.tins_points[e].map(o=>V(o,t.points,t.edgeNodes||[],t.centroid_point,t.vertices_points,!0,t.version)))},this.addIndexedTin(),t.kinks_points&&(this.kinks={bakw:F(t.kinks_points.map(o=>W(o)))}),t.yaxisMode?this.yaxisMode=t.yaxisMode:this.yaxisMode=O.YAXIS_INVERT,t.vertexMode&&(this.vertexMode=t.vertexMode),t.strictMode&&(this.strictMode=t.strictMode),t.bounds?(this.bounds=t.bounds,this.boundsPolygon=t.boundsPolygon,this.xy=t.xy,this.wh=t.wh):(this.xy=[0,0],t.wh&&(this.wh=t.wh),this.bounds=void 0,this.boundsPolygon=void 0)}else{t=JSON.parse(JSON.stringify(t).replace('"cent"','"c"').replace(/"bbox(\d+)"/g,'"b$1"')),this.tins=t.tins,this.addIndexedTin(),this.strict_status=t.strict_status,this.pointsWeightBuffer=t.weight_buffer,this.vertices_params=t.vertices_params,this.centroid=t.centroid,this.kinks=t.kinks;const e=[];for(let o=0;o<this.tins.forw.features.length;o++){const f=this.tins.forw.features[o];["a","b","c"].map((s,c)=>{const g=f.geometry.coordinates[0][c],l=f.properties[s].geom,y=f.properties[s].index;e[y]=[g,l]})}this.points=e}return{tins:this.tins,strict_status:this.strict_status,weight_buffer:this.pointsWeightBuffer,vertices_params:this.vertices_params,centroid:this.centroid,kinks:this.kinks}}addIndexedTin(){const t=this.tins,e=t.forw,o=t.bakw,f=Math.ceil(Math.sqrt(e.features.length));if(f<3){this.indexedTins=void 0;return}let s=[],c=[];const g=e.features.map(h=>{let a=[];return Q(h)[0].map(n=>{s.length===0?s=[Array.from(n),Array.from(n)]:(n[0]<s[0][0]&&(s[0][0]=n[0]),n[0]>s[1][0]&&(s[1][0]=n[0]),n[1]<s[0][1]&&(s[0][1]=n[1]),n[1]>s[1][1]&&(s[1][1]=n[1])),a.length===0?a=[Array.from(n),Array.from(n)]:(n[0]<a[0][0]&&(a[0][0]=n[0]),n[0]>a[1][0]&&(a[1][0]=n[0]),n[1]<a[0][1]&&(a[0][1]=n[1]),n[1]>a[1][1]&&(a[1][1]=n[1]))}),a}),l=(s[1][0]-s[0][0])/f,y=(s[1][1]-s[0][1])/f,w=g.reduce((h,a,n)=>{const b=R(a[0][0],s[0][0],l,f),M=R(a[1][0],s[0][0],l,f),x=R(a[0][1],s[0][1],y,f),E=R(a[1][1],s[0][1],y,f);for(let _=b;_<=M;_++){h[_]||(h[_]=[]);for(let m=x;m<=E;m++)h[_][m]||(h[_][m]=[]),h[_][m].push(n)}return h},[]),i=o.features.map(h=>{let a=[];return Q(h)[0].map(n=>{c.length===0?c=[Array.from(n),Array.from(n)]:(n[0]<c[0][0]&&(c[0][0]=n[0]),n[0]>c[1][0]&&(c[1][0]=n[0]),n[1]<c[0][1]&&(c[0][1]=n[1]),n[1]>c[1][1]&&(c[1][1]=n[1])),a.length===0?a=[Array.from(n),Array.from(n)]:(n[0]<a[0][0]&&(a[0][0]=n[0]),n[0]>a[1][0]&&(a[1][0]=n[0]),n[1]<a[0][1]&&(a[0][1]=n[1]),n[1]>a[1][1]&&(a[1][1]=n[1]))}),a}),u=(c[1][0]-c[0][0])/f,d=(c[1][1]-c[0][1])/f,v=i.reduce((h,a,n)=>{const b=R(a[0][0],c[0][0],u,f),M=R(a[1][0],c[0][0],u,f),x=R(a[0][1],c[0][1],d,f),E=R(a[1][1],c[0][1],d,f);for(let _=b;_<=M;_++){h[_]||(h[_]=[]);for(let m=x;m<=E;m++)h[_][m]||(h[_][m]=[]),h[_][m].push(n)}return h},[]);this.indexedTins={forw:{gridNum:f,xOrigin:s[0][0],yOrigin:s[0][1],xUnit:l,yUnit:y,gridCache:w},bakw:{gridNum:f,xOrigin:c[0][0],yOrigin:c[0][1],xUnit:u,yUnit:d,gridCache:v}}}transform(t,e,o){if(e&&this.strict_status==O.STATUS_ERROR)throw'Backward transform is not allowed if strict_status == "strict_error"';this.yaxisMode==O.YAXIS_FOLLOW&&e&&(t=[t[0],-1*t[1]]);const f=W(t);if(this.bounds&&!e&&!o&&!q(f,this.boundsPolygon))return!1;const s=e?this.tins.bakw:this.tins.forw,c=e?this.indexedTins.bakw:this.indexedTins.forw,g=e?this.vertices_params.bakw:this.vertices_params.forw,l=e?this.centroid.bakw:this.centroid.forw,y=e?this.pointsWeightBuffer.bakw:this.pointsWeightBuffer.forw;let w,i;this.stateFull&&(this.stateBackward==e?w=this.stateTriangle:(this.stateBackward=e,this.stateTriangle=void 0),i=d=>{this.stateTriangle=d});let u=mt(f,s,c,g,l,y,w,i);if(this.bounds&&e&&!o){const d=W(u);if(!q(d,this.boundsPolygon))return!1}else this.yaxisMode==O.YAXIS_FOLLOW&&!e&&(u=[u[0],-1*u[1]]);return u}};A(O,"VERTEX_PLAIN","plain"),A(O,"VERTEX_BIRDEYE","birdeye"),A(O,"MODE_STRICT","strict"),A(O,"MODE_AUTO","auto"),A(O,"MODE_LOOSE","loose"),A(O,"STATUS_STRICT","strict"),A(O,"STATUS_ERROR","strict_error"),A(O,"STATUS_LOOSE","loose"),A(O,"YAXIS_FOLLOW","follow"),A(O,"YAXIS_INVERT","invert");let G=O;exports.Transform=G;exports.format_version=L;

package/dist/maplat_transform.js

@@ -1,545 +0,0 @@
-var it = Object.defineProperty;
-var ot = (r, t, e) => t in r ? it(r, t, { enumerable: !0, configurable: !0, writable: !0, value: e }) : r[t] = e;
-var A = (r, t, e) => ot(r, typeof t != "symbol" ? t + "" : t, e);
-const B = 11102230246251565e-32, T = 134217729, at = (3 + 8 * B) * B;
-function $(r, t, e, o, f) {
-  let s, c, g, l, y = t[0], w = o[0], i = 0, u = 0;
-  w > y == w > -y ? (s = y, y = t[++i]) : (s = w, w = o[++u]);
-  let d = 0;
-  if (i < r && u < e)
-    for (w > y == w > -y ? (c = y + s, g = s - (c - y), y = t[++i]) : (c = w + s, g = s - (c - w), w = o[++u]), s = c, g !== 0 && (f[d++] = g); i < r && u < e; )
-      w > y == w > -y ? (c = s + y, l = c - s, g = s - (c - l) + (y - l), y = t[++i]) : (c = s + w, l = c - s, g = s - (c - l) + (w - l), w = o[++u]), s = c, g !== 0 && (f[d++] = g);
-  for (; i < r; )
-    c = s + y, l = c - s, g = s - (c - l) + (y - l), y = t[++i], s = c, g !== 0 && (f[d++] = g);
-  for (; u < e; )
-    c = s + w, l = c - s, g = s - (c - l) + (w - l), w = o[++u], s = c, g !== 0 && (f[d++] = g);
-  return (s !== 0 || d === 0) && (f[d++] = s), d;
-}
-function ft(r, t) {
-  let e = t[0];
-  for (let o = 1; o < r; o++) e += t[o];
-  return e;
-}
-function D(r) {
-  return new Float64Array(r);
-}
-const ct = (3 + 16 * B) * B, ut = (2 + 12 * B) * B, ht = (9 + 64 * B) * B * B, X = D(4), L = D(8), G = D(12), j = D(16), I = D(4);
-function lt(r, t, e, o, f, s, c) {
-  let g, l, y, w, i, u, d, v, h, a, n, b, x, M, m, _, E, k;
-  const S = r - f, P = e - f, N = t - s, C = o - s;
-  M = S * C, u = T * S, d = u - (u - S), v = S - d, u = T * C, h = u - (u - C), a = C - h, m = v * a - (M - d * h - v * h - d * a), _ = N * P, u = T * N, d = u - (u - N), v = N - d, u = T * P, h = u - (u - P), a = P - h, E = v * a - (_ - d * h - v * h - d * a), n = m - E, i = m - n, X[0] = m - (n + i) + (i - E), b = M + n, i = b - M, x = M - (b - i) + (n - i), n = x - _, i = x - n, X[1] = x - (n + i) + (i - _), k = b + n, i = k - b, X[2] = b - (k - i) + (n - i), X[3] = k;
-  let U = ft(4, X), Y = ut * c;
-  if (U >= Y || -U >= Y || (i = r - S, g = r - (S + i) + (i - f), i = e - P, y = e - (P + i) + (i - f), i = t - N, l = t - (N + i) + (i - s), i = o - C, w = o - (C + i) + (i - s), g === 0 && l === 0 && y === 0 && w === 0) || (Y = ht * c + at * Math.abs(U), U += S * w + C * g - (N * y + P * l), U >= Y || -U >= Y)) return U;
-  M = g * C, u = T * g, d = u - (u - g), v = g - d, u = T * C, h = u - (u - C), a = C - h, m = v * a - (M - d * h - v * h - d * a), _ = l * P, u = T * l, d = u - (u - l), v = l - d, u = T * P, h = u - (u - P), a = P - h, E = v * a - (_ - d * h - v * h - d * a), n = m - E, i = m - n, I[0] = m - (n + i) + (i - E), b = M + n, i = b - M, x = M - (b - i) + (n - i), n = x - _, i = x - n, I[1] = x - (n + i) + (i - _), k = b + n, i = k - b, I[2] = b - (k - i) + (n - i), I[3] = k;
-  const et = $(4, X, 4, I, L);
-  M = S * w, u = T * S, d = u - (u - S), v = S - d, u = T * w, h = u - (u - w), a = w - h, m = v * a - (M - d * h - v * h - d * a), _ = N * y, u = T * N, d = u - (u - N), v = N - d, u = T * y, h = u - (u - y), a = y - h, E = v * a - (_ - d * h - v * h - d * a), n = m - E, i = m - n, I[0] = m - (n + i) + (i - E), b = M + n, i = b - M, x = M - (b - i) + (n - i), n = x - _, i = x - n, I[1] = x - (n + i) + (i - _), k = b + n, i = k - b, I[2] = b - (k - i) + (n - i), I[3] = k;
-  const st = $(et, L, 4, I, G);
-  M = g * w, u = T * g, d = u - (u - g), v = g - d, u = T * w, h = u - (u - w), a = w - h, m = v * a - (M - d * h - v * h - d * a), _ = l * y, u = T * l, d = u - (u - l), v = l - d, u = T * y, h = u - (u - y), a = y - h, E = v * a - (_ - d * h - v * h - d * a), n = m - E, i = m - n, I[0] = m - (n + i) + (i - E), b = M + n, i = b - M, x = M - (b - i) + (n - i), n = x - _, i = x - n, I[1] = x - (n + i) + (i - _), k = b + n, i = k - b, I[2] = b - (k - i) + (n - i), I[3] = k;
-  const nt = $(st, G, 4, I, j);
-  return j[nt - 1];
-}
-function gt(r, t, e, o, f, s) {
-  const c = (t - s) * (e - f), g = (r - f) * (o - s), l = c - g, y = Math.abs(c + g);
-  return Math.abs(l) >= ct * y ? l : -lt(r, t, e, o, f, s, y);
-}
-function yt(r, t) {
-  var e, o, f = 0, s, c, g, l, y, w, i, u = r[0], d = r[1], v = t.length;
-  for (e = 0; e < v; e++) {
-    o = 0;
-    var h = t[e], a = h.length - 1;
-    if (w = h[0], w[0] !== h[a][0] && w[1] !== h[a][1])
-      throw new Error("First and last coordinates in a ring must be the same");
-    for (c = w[0] - u, g = w[1] - d, o; o < a; o++) {
-      if (i = h[o + 1], l = i[0] - u, y = i[1] - d, g === 0 && y === 0) {
-        if (l <= 0 && c >= 0 || c <= 0 && l >= 0)
-          return 0;
-      } else if (y >= 0 && g <= 0 || y <= 0 && g >= 0) {
-        if (s = gt(c, l, g, y, 0, 0), s === 0)
-          return 0;
-        (s > 0 && y > 0 && g <= 0 || s < 0 && y <= 0 && g > 0) && f++;
-      }
-      w = i, g = y, c = l;
-    }
-  }
-  return f % 2 !== 0;
-}
-function p(r, t, e = {}) {
-  const o = { type: "Feature" };
-  return (e.id === 0 || e.id) && (o.id = e.id), e.bbox && (o.bbox = e.bbox), o.properties = t || {}, o.geometry = r, o;
-}
-function W(r, t, e = {}) {
-  if (!r)
-    throw new Error("coordinates is required");
-  if (!Array.isArray(r))
-    throw new Error("coordinates must be an Array");
-  if (r.length < 2)
-    throw new Error("coordinates must be at least 2 numbers long");
-  if (!J(r[0]) || !J(r[1]))
-    throw new Error("coordinates must contain numbers");
-  return p({
-    type: "Point",
-    coordinates: r
-  }, t, e);
-}
-function wt(r, t, e = {}) {
-  for (const f of r) {
-    if (f.length < 4)
-      throw new Error(
-        "Each LinearRing of a Polygon must have 4 or more Positions."
-      );
-    if (f[f.length - 1].length !== f[0].length)
-      throw new Error("First and last Position are not equivalent.");
-    for (let s = 0; s < f[f.length - 1].length; s++)
-      if (f[f.length - 1][s] !== f[0][s])
-        throw new Error("First and last Position are not equivalent.");
-  }
-  return p({
-    type: "Polygon",
-    coordinates: r
-  }, t, e);
-}
-function F(r, t = {}) {
-  const e = { type: "FeatureCollection" };
-  return t.id && (e.id = t.id), t.bbox && (e.bbox = t.bbox), e.features = r, e;
-}
-function J(r) {
-  return !isNaN(r) && r !== null && !Array.isArray(r);
-}
-function dt(r) {
-  if (!r)
-    throw new Error("coord is required");
-  if (!Array.isArray(r)) {
-    if (r.type === "Feature" && r.geometry !== null && r.geometry.type === "Point")
-      return [...r.geometry.coordinates];
-    if (r.type === "Point")
-      return [...r.coordinates];
-  }
-  if (Array.isArray(r) && r.length >= 2 && !Array.isArray(r[0]) && !Array.isArray(r[1]))
-    return [...r];
-  throw new Error("coord must be GeoJSON Point or an Array of numbers");
-}
-function z(r) {
-  if (Array.isArray(r))
-    return r;
-  if (r.type === "Feature") {
-    if (r.geometry !== null)
-      return r.geometry.coordinates;
-  } else if (r.coordinates)
-    return r.coordinates;
-  throw new Error(
-    "coords must be GeoJSON Feature, Geometry Object or an Array"
-  );
-}
-function bt(r) {
-  return r.type === "Feature" ? r.geometry : r;
-}
-function _t(r, t, e = {}) {
-  if (!r)
-    throw new Error("point is required");
-  if (!t)
-    throw new Error("polygon is required");
-  const o = dt(r), f = bt(t), s = f.type, c = t.bbox;
-  let g = f.coordinates;
-  if (c && vt(o, c) === !1)
-    return !1;
-  s === "Polygon" && (g = [g]);
-  let l = !1;
-  for (var y = 0; y < g.length; ++y) {
-    const w = yt(o, g[y]);
-    if (w === 0) return !e.ignoreBoundary;
-    w && (l = !0);
-  }
-  return l;
-}
-function vt(r, t) {
-  return t[0] <= r[0] && t[1] <= r[1] && t[2] >= r[0] && t[3] >= r[1];
-}
-var q = _t;
-function At(r) {
-  const t = [0, 1, 2, 0].map((o) => r[o][0][0]), e = {
-    a: { geom: r[0][0][1], index: r[0][1] },
-    b: { geom: r[1][0][1], index: r[1][1] },
-    c: { geom: r[2][0][1], index: r[2][1] }
-  };
-  return wt([t], e);
-}
-function V(r, t, e, o, f, s = !1, c) {
-  const g = r.map(
-    (l) => {
-      (!c || c < 2.00703) && (l = tt(l));
-      const y = isFinite(l) ? t[l] : l === "c" ? o : l === "b0" ? f[0] : l === "b1" ? f[1] : l === "b2" ? f[2] : l === "b3" ? f[3] : function() {
-        const w = l.match(/e(\d+)/);
-        if (w) {
-          const i = parseInt(w[1]);
-          return e[i];
-        }
-        throw "Bad index value for indexesToTri";
-      }();
-      return s ? [[y[1], y[0]], l] : [[y[0], y[1]], l];
-    }
-  );
-  return At(g);
-}
-function tt(r) {
-  return typeof r == "number" ? r : r.replace(/^(c|e|b)(?:ent|dgeNode|box)(\d+)?$/, "$1$2");
-}
-function xt(r, t) {
-  return Array.isArray(r[0]) ? r : r.map((e) => [
-    e.illstNodes,
-    e.mercNodes,
-    e.startEnd
-  ]);
-}
-function K(r, t) {
-  for (let e = 0; e < t.features.length; e++)
-    if (q(r, t.features[e]))
-      return t.features[e];
-}
-function rt(r, t, e) {
-  const o = t.geometry.coordinates[0][0], f = t.geometry.coordinates[0][1], s = t.geometry.coordinates[0][2], c = r.geometry.coordinates, g = t.properties.a.geom, l = t.properties.b.geom, y = t.properties.c.geom, w = [f[0] - o[0], f[1] - o[1]], i = [s[0] - o[0], s[1] - o[1]], u = [c[0] - o[0], c[1] - o[1]], d = [l[0] - g[0], l[1] - g[1]], v = [y[0] - g[0], y[1] - g[1]];
-  let h = (i[1] * u[0] - i[0] * u[1]) / (w[0] * i[1] - w[1] * i[0]), a = (w[0] * u[1] - w[1] * u[0]) / (w[0] * i[1] - w[1] * i[0]);
-  if (e) {
-    const n = e[t.properties.a.index], b = e[t.properties.b.index], x = e[t.properties.c.index];
-    let M;
-    if (h < 0 || a < 0 || 1 - h - a < 0) {
-      const m = h / (h + a), _ = a / (h + a);
-      M = h / b / (m / b + _ / x), a = a / x / (m / b + _ / x);
-    } else
-      M = h / b / (h / b + a / x + (1 - h - a) / n), a = a / x / (h / b + a / x + (1 - h - a) / n);
-    h = M;
-  }
-  return [
-    h * d[0] + a * v[0] + g[0],
-    h * d[1] + a * v[1] + g[1]
-  ];
-}
-function Mt(r, t, e, o) {
-  const f = r.geometry.coordinates, s = e.geometry.coordinates, c = Math.atan2(f[0] - s[0], f[1] - s[1]), g = mt(c, t[0]), l = t[1][g];
-  return rt(r, l.features[0], o);
-}
-function Et(r, t, e, o, f, s, c, g) {
-  let l;
-  if (c && (l = K(r, F([c]))), !l) {
-    if (e) {
-      const y = r.geometry.coordinates, w = e.gridNum, i = e.xOrigin, u = e.yOrigin, d = e.xUnit, v = e.yUnit, h = e.gridCache, a = R(y[0], i, d, w), n = R(y[1], u, v, w), b = h[a] ? h[a][n] ? h[a][n] : [] : [];
-      t = F(b.map((x) => t.features[x]));
-    }
-    l = K(r, t);
-  }
-  return g && g(l), l ? rt(r, l, s) : Mt(r, o, f, s);
-}
-function R(r, t, e, o) {
-  let f = Math.floor((r - t) / e);
-  return f >= o && (f = o - 1), f;
-}
-function mt(r, t) {
-  let e = Q(r - t[0]), o = Math.PI * 2, f;
-  for (let s = 0; s < t.length; s++) {
-    const c = (s + 1) % t.length, g = Q(r - t[c]), l = Math.min(Math.abs(e), Math.abs(g));
-    e * g <= 0 && l < o && (o = l, f = s), e = g;
-  }
-  return f;
-}
-function Q(r, t = !1) {
-  const e = t ? function(o) {
-    return !(o >= 0 && o < Math.PI * 2);
-  } : function(o) {
-    return !(o > -1 * Math.PI && o <= Math.PI);
-  };
-  for (; e(r); )
-    r = r + 2 * Math.PI * (r > 0 ? -1 : 1);
-  return r;
-}
-const H = 2.00703, O = class O {
-  constructor() {
-    A(this, "points", []);
-    A(this, "pointsWeightBuffer");
-    A(this, "strict_status");
-    A(this, "vertices_params");
-    A(this, "centroid");
-    A(this, "edgeNodes");
-    A(this, "edges");
-    A(this, "tins");
-    A(this, "kinks");
-    A(this, "yaxisMode", O.YAXIS_INVERT);
-    A(this, "strictMode", O.MODE_AUTO);
-    A(this, "vertexMode", O.VERTEX_PLAIN);
-    A(this, "bounds");
-    A(this, "boundsPolygon");
-    A(this, "wh");
-    A(this, "xy");
-    A(this, "indexedTins");
-    A(this, "stateFull", !1);
-    A(this, "stateTriangle");
-    A(this, "stateBackward");
-  }
-  /**
-   * コンパイルされた設定を適用します
-   * 
-   * @param compiled - コンパイルされた設定オブジェクト
-   * @returns 変換に必要な主要なオブジェクトのセット
-   * 
-   * 以下の処理を行います：
-   * 1. バージョンに応じた設定の解釈
-   * 2. 各種パラメータの復元
-   * 3. TINネットワークの再構築
-   * 4. インデックスの作成
-   */
-  setCompiled(t) {
-    if (t.version || !t.tins && t.points && t.tins_points) {
-      this.points = t.points, this.pointsWeightBuffer = !t.version || t.version < 2.00703 ? ["forw", "bakw"].reduce((o, f) => {
-        const s = t.weight_buffer[f];
-        return s && (o[f] = Object.keys(s).reduce((c, g) => {
-          const l = tt(g);
-          return c[l] = s[g], c;
-        }, {})), o;
-      }, {}) : t.weight_buffer, t.strict_status ? this.strict_status = t.strict_status : t.kinks_points ? this.strict_status = O.STATUS_ERROR : t.tins_points.length == 2 ? this.strict_status = O.STATUS_LOOSE : this.strict_status = O.STATUS_STRICT, this.vertices_params = {
-        forw: [t.vertices_params[0]],
-        bakw: [t.vertices_params[1]]
-      }, this.vertices_params.forw[1] = [0, 1, 2, 3].map((o) => {
-        const f = (o + 1) % 4, s = V(
-          ["c", `b${o}`, `b${f}`],
-          t.points,
-          t.edgeNodes || [],
-          t.centroid_point,
-          t.vertices_points,
-          !1,
-          H
-        );
-        return F([s]);
-      }), this.vertices_params.bakw[1] = [0, 1, 2, 3].map((o) => {
-        const f = (o + 1) % 4, s = V(
-          ["c", `b${o}`, `b${f}`],
-          t.points,
-          t.edgeNodes || [],
-          t.centroid_point,
-          t.vertices_points,
-          !0,
-          H
-        );
-        return F([s]);
-      }), this.centroid = {
-        forw: W(t.centroid_point[0], {
-          target: {
-            geom: t.centroid_point[1],
-            index: "c"
-          }
-        }),
-        bakw: W(t.centroid_point[1], {
-          target: {
-            geom: t.centroid_point[0],
-            index: "c"
-          }
-        })
-      }, this.edges = xt(t.edges || []), this.edgeNodes = t.edgeNodes || [];
-      const e = t.tins_points.length == 1 ? 0 : 1;
-      this.tins = {
-        forw: F(
-          t.tins_points[0].map(
-            (o) => V(
-              o,
-              t.points,
-              t.edgeNodes || [],
-              t.centroid_point,
-              t.vertices_points,
-              !1,
-              t.version
-            )
-          )
-        ),
-        bakw: F(
-          t.tins_points[e].map(
-            (o) => V(
-              o,
-              t.points,
-              t.edgeNodes || [],
-              t.centroid_point,
-              t.vertices_points,
-              !0,
-              t.version
-            )
-          )
-        )
-      }, this.addIndexedTin(), t.kinks_points && (this.kinks = {
-        bakw: F(
-          t.kinks_points.map((o) => W(o))
-        )
-      }), t.yaxisMode ? this.yaxisMode = t.yaxisMode : this.yaxisMode = O.YAXIS_INVERT, t.vertexMode && (this.vertexMode = t.vertexMode), t.strictMode && (this.strictMode = t.strictMode), t.bounds ? (this.bounds = t.bounds, this.boundsPolygon = t.boundsPolygon, this.xy = t.xy, this.wh = t.wh) : (this.xy = [0, 0], t.wh && (this.wh = t.wh), this.bounds = void 0, this.boundsPolygon = void 0);
-    } else {
-      t = JSON.parse(
-        JSON.stringify(t).replace('"cent"', '"c"').replace(/"bbox(\d+)"/g, '"b$1"')
-      ), this.tins = t.tins, this.addIndexedTin(), this.strict_status = t.strict_status, this.pointsWeightBuffer = t.weight_buffer, this.vertices_params = t.vertices_params, this.centroid = t.centroid, this.kinks = t.kinks;
-      const e = [];
-      for (let o = 0; o < this.tins.forw.features.length; o++) {
-        const f = this.tins.forw.features[o];
-        ["a", "b", "c"].map((s, c) => {
-          const g = f.geometry.coordinates[0][c], l = f.properties[s].geom, y = f.properties[s].index;
-          e[y] = [g, l];
-        });
-      }
-      this.points = e;
-    }
-    return {
-      tins: this.tins,
-      strict_status: this.strict_status,
-      weight_buffer: this.pointsWeightBuffer,
-      vertices_params: this.vertices_params,
-      centroid: this.centroid,
-      kinks: this.kinks
-    };
-  }
-  /**
-   * TINネットワークのインデックスを作成します
-   * 
-   * インデックスは変換処理を高速化するために使用されます。
-   * グリッド形式のインデックスを作成し、各グリッドに
-   * 含まれる三角形を記録します。
-   */
-  addIndexedTin() {
-    const t = this.tins, e = t.forw, o = t.bakw, f = Math.ceil(Math.sqrt(e.features.length));
-    if (f < 3) {
-      this.indexedTins = void 0;
-      return;
-    }
-    let s = [], c = [];
-    const g = e.features.map((h) => {
-      let a = [];
-      return z(h)[0].map((n) => {
-        s.length === 0 ? s = [Array.from(n), Array.from(n)] : (n[0] < s[0][0] && (s[0][0] = n[0]), n[0] > s[1][0] && (s[1][0] = n[0]), n[1] < s[0][1] && (s[0][1] = n[1]), n[1] > s[1][1] && (s[1][1] = n[1])), a.length === 0 ? a = [Array.from(n), Array.from(n)] : (n[0] < a[0][0] && (a[0][0] = n[0]), n[0] > a[1][0] && (a[1][0] = n[0]), n[1] < a[0][1] && (a[0][1] = n[1]), n[1] > a[1][1] && (a[1][1] = n[1]));
-      }), a;
-    }), l = (s[1][0] - s[0][0]) / f, y = (s[1][1] - s[0][1]) / f, w = g.reduce(
-      (h, a, n) => {
-        const b = R(
-          a[0][0],
-          s[0][0],
-          l,
-          f
-        ), x = R(
-          a[1][0],
-          s[0][0],
-          l,
-          f
-        ), M = R(
-          a[0][1],
-          s[0][1],
-          y,
-          f
-        ), m = R(
-          a[1][1],
-          s[0][1],
-          y,
-          f
-        );
-        for (let _ = b; _ <= x; _++) {
-          h[_] || (h[_] = []);
-          for (let E = M; E <= m; E++)
-            h[_][E] || (h[_][E] = []), h[_][E].push(n);
-        }
-        return h;
-      },
-      []
-    ), i = o.features.map((h) => {
-      let a = [];
-      return z(h)[0].map((n) => {
-        c.length === 0 ? c = [Array.from(n), Array.from(n)] : (n[0] < c[0][0] && (c[0][0] = n[0]), n[0] > c[1][0] && (c[1][0] = n[0]), n[1] < c[0][1] && (c[0][1] = n[1]), n[1] > c[1][1] && (c[1][1] = n[1])), a.length === 0 ? a = [Array.from(n), Array.from(n)] : (n[0] < a[0][0] && (a[0][0] = n[0]), n[0] > a[1][0] && (a[1][0] = n[0]), n[1] < a[0][1] && (a[0][1] = n[1]), n[1] > a[1][1] && (a[1][1] = n[1]));
-      }), a;
-    }), u = (c[1][0] - c[0][0]) / f, d = (c[1][1] - c[0][1]) / f, v = i.reduce(
-      (h, a, n) => {
-        const b = R(
-          a[0][0],
-          c[0][0],
-          u,
-          f
-        ), x = R(
-          a[1][0],
-          c[0][0],
-          u,
-          f
-        ), M = R(
-          a[0][1],
-          c[0][1],
-          d,
-          f
-        ), m = R(
-          a[1][1],
-          c[0][1],
-          d,
-          f
-        );
-        for (let _ = b; _ <= x; _++) {
-          h[_] || (h[_] = []);
-          for (let E = M; E <= m; E++)
-            h[_][E] || (h[_][E] = []), h[_][E].push(n);
-        }
-        return h;
-      },
-      []
-    );
-    this.indexedTins = {
-      forw: {
-        gridNum: f,
-        xOrigin: s[0][0],
-        yOrigin: s[0][1],
-        xUnit: l,
-        yUnit: y,
-        gridCache: w
-      },
-      bakw: {
-        gridNum: f,
-        xOrigin: c[0][0],
-        yOrigin: c[0][1],
-        xUnit: u,
-        yUnit: d,
-        gridCache: v
-      }
-    };
-  }
-  /**
-   * 座標変換を実行します
-   * 
-   * @param apoint - 変換する座標
-   * @param backward - 逆方向の変換かどうか
-   * @param ignoreBounds - 境界チェックを無視するかどうか
-   * @returns 変換後の座標、または境界外の場合はfalse
-   * 
-   * @throws {Error} 逆方向変換が許可されていない状態での逆変換時
-   */
-  transform(t, e, o) {
-    if (e && this.strict_status == O.STATUS_ERROR)
-      throw 'Backward transform is not allowed if strict_status == "strict_error"';
-    this.yaxisMode == O.YAXIS_FOLLOW && e && (t = [t[0], -1 * t[1]]);
-    const f = W(t);
-    if (this.bounds && !e && !o && !q(f, this.boundsPolygon))
-      return !1;
-    const s = e ? this.tins.bakw : this.tins.forw, c = e ? this.indexedTins.bakw : this.indexedTins.forw, g = e ? this.vertices_params.bakw : this.vertices_params.forw, l = e ? this.centroid.bakw : this.centroid.forw, y = e ? this.pointsWeightBuffer.bakw : this.pointsWeightBuffer.forw;
-    let w, i;
-    this.stateFull && (this.stateBackward == e ? w = this.stateTriangle : (this.stateBackward = e, this.stateTriangle = void 0), i = (d) => {
-      this.stateTriangle = d;
-    });
-    let u = Et(
-      f,
-      s,
-      c,
-      g,
-      l,
-      y,
-      w,
-      i
-    );
-    if (this.bounds && e && !o) {
-      const d = W(u);
-      if (!q(d, this.boundsPolygon)) return !1;
-    } else this.yaxisMode == O.YAXIS_FOLLOW && !e && (u = [u[0], -1 * u[1]]);
-    return u;
-  }
-};
-/**
- * 各種モードの定数定義
- * すべてreadonlyで、型安全性を確保
- */
-A(O, "VERTEX_PLAIN", "plain"), A(O, "VERTEX_BIRDEYE", "birdeye"), A(O, "MODE_STRICT", "strict"), A(O, "MODE_AUTO", "auto"), A(O, "MODE_LOOSE", "loose"), A(O, "STATUS_STRICT", "strict"), A(O, "STATUS_ERROR", "strict_error"), A(O, "STATUS_LOOSE", "loose"), A(O, "YAXIS_FOLLOW", "follow"), A(O, "YAXIS_INVERT", "invert");
-let Z = O;
-export {
-  Z as Transform,
-  H as format_version
-};