2 lines
9.1 KiB
JavaScript
2 lines
9.1 KiB
JavaScript
|
|
define(["exports","./Transforms-a7d32f6c","./Cartesian2-40b13f31","./Check-3917d765","./when-7b580518","./EllipsoidTangentPlane-b7b148bd","./Math-87254c13","./Plane-69fa3eea"],function(a,b,O,t,P,w,y,N){"use strict";function T(a,t){this.center=O.Cartesian3.clone(P.defaultValue(a,O.Cartesian3.ZERO)),this.halfAxes=b.Matrix3.clone(P.defaultValue(t,b.Matrix3.ZERO))}T.packedLength=O.Cartesian3.packedLength+b.Matrix3.packedLength,T.pack=function(a,t,e){return e=P.defaultValue(e,0),O.Cartesian3.pack(a.center,t,e),b.Matrix3.pack(a.halfAxes,t,e+O.Cartesian3.packedLength),t},T.unpack=function(a,t,e){return t=P.defaultValue(t,0),P.defined(e)||(e=new T),O.Cartesian3.unpack(a,t,e.center),b.Matrix3.unpack(a,t+O.Cartesian3.packedLength,e.halfAxes),e};var A=new O.Cartesian3,I=new O.Cartesian3,R=new O.Cartesian3,E=new O.Cartesian3,S=new O.Cartesian3,U=new O.Cartesian3,L=new b.Matrix3,v={unitary:new b.Matrix3,diagonal:new b.Matrix3};T.fromPoints=function(a,t){if(P.defined(t)||(t=new T),!P.defined(a)||0===a.length)return t.halfAxes=b.Matrix3.ZERO,t.center=O.Cartesian3.ZERO,t;for(var e=a.length,n=O.Cartesian3.clone(a[0],A),r=1;r<e;r++)O.Cartesian3.add(n,a[r],n);var i=1/e;O.Cartesian3.multiplyByScalar(n,i,n);var s,o=0,C=0,c=0,u=0,d=0,l=0;for(r=0;r<e;r++)o+=(s=O.Cartesian3.subtract(a[r],n,I)).x*s.x,C+=s.x*s.y,c+=s.x*s.z,u+=s.y*s.y,d+=s.y*s.z,l+=s.z*s.z;C*=i,c*=i,u*=i,d*=i,l*=i;var h=L;h[0]=o*=i,h[1]=C,h[2]=c,h[3]=C,h[4]=u,h[5]=d,h[6]=c,h[7]=d,h[8]=l;var h=b.Matrix3.computeEigenDecomposition(h,v),h=b.Matrix3.clone(h.unitary,t.halfAxes),x=b.Matrix3.getColumn(h,0,E),M=b.Matrix3.getColumn(h,1,S),m=b.Matrix3.getColumn(h,2,U),f=-Number.MAX_VALUE,p=-Number.MAX_VALUE,g=-Number.MAX_VALUE,w=Number.MAX_VALUE,y=Number.MAX_VALUE,N=Number.MAX_VALUE;for(r=0;r<e;r++)s=a[r],f=Math.max(O.Cartesian3.dot(x,s),f),p=Math.max(O.Cartesian3.dot(M,s),p),g=Math.max(O.Cartesian3.dot(m,s),g),w=Math.min(O.Cartesian3.dot(x,s),w),y=Math.min(O.Cartesian3.dot(M,s),y),N=Math.min(O.Cartesian3.dot(m,s),N);x=O.Cartesian3.multiplyByScalar(x,.5*(w+f),x),M=O.Cartesian3.multiplyByScalar(M,.5*(y+p),M),m=O.Cartesian3.multiplyByScalar(m,.5*(N+g),m),h=O.Cartesian3.add(x,M,t.center);O.Cartesian3.add(h,m,h);h=R;return h.x=f-w,h.y=p-y,h.z=g-N,O.Cartesian3.multiplyByScalar(h,.5,h),b.Matrix3.multiplyByScale(t.halfAxes,h,t.halfAxes),t};var M=new O.Cartesian3,l=new O.Cartesian3;function z(a,t,e,n,r,i,s,o,C,c,u){var d=(u=!P.defined(u)?new T:u).halfAxes;b.Matrix3.setColumn(d,0,t,d),b.Matrix3.setColumn(d,1,e,d),b.Matrix3.setColumn(d,2,n,d),(e=M).x=(r+i)/2,e.y=(s+o)/2,e.z=(C+c)/2;n=l;n.x=(i-r)/2,n.y=(o-s)/2,n.z=(c-C)/2;C=u.center,e=b.Matrix3.multiplyByVector(d,e,e);return O.Cartesian3.add(a,e,C),b.Matrix3.multiplyByScale(d,n,d),u}var B=new O.Cartographic,V=new O.Cartesian3,_=new O.Cartographic,k=new O.Cartographic,W=new O.Cartographic,X=new O.Cartographic,q=new O.Cartographic,D=new O.Cartesian3,j=new O.Cartesian3,Z=new O.Cartesian3,Y=new O.Cartesian3,G=new O.Cartesian3,F=new O.Cartesian2,H=new O.Cartesian2,J=new O.Cartesian2,K=new O.Cartesian2,Q=new O.Cartesian2,$=new O.Cartesian3,aa=new O.Cartesian3,ta=new O.Cartesian3,ea=new O.Cartesian3,na=new O.Cartesian2,ra=new O.Cartesian3,ia=new O.Cartesian3,sa=new O.Cartesian3,oa=new N.Plane(O.Cartesian3.UNIT_X,0);T.fromRectangle=function(a,t,e,n,r){if(t=P.defaultValue(t,0),e=P.defaultValue(e,0),n=P.defaultValue(n,O.Ellipsoid.WGS84),a.width<=y.CesiumMath.PI){var i=O.Rectangle.center(a,B),s=n.cartographicToCartesian(i,V),o=new w.EllipsoidTangentPlane(s,n),C=o.plane,c=i.longitude,u=a.south<0&&0<a.north?0:i.latitude,d=O.Cartographic.fromRadians(c,a.north,e,_),l=O.Cartographic.fromRadians(a.west,a.north,e,k),h=O.Cartographic.fromRadians(a.west,u,e,W),x=O.Cartographic.fromRadians(a.west,a.south,e,X),M=O.Cartographic.fromRadians(c,a.south,e,q),m=n.cartographicToCartesian(d,D),f=n.cartographicToCartesian(l,j),p=n.cartographicToCartesian(h,Z),g=n.cartographicToCartesian(x,Y),s=n.cartographicToCartesian(M,G),i=o.projectPointToNearestOnPlane(m,F),u=o.projectPointToNearestOnPlane(f,H),c=o.projectPointToNearestOnPlane(p,J),d=o.projectPointToNearestOnPlane(g,K),h=o.proj
|