[World Wind学习]22.相机高度和瓦片等级计算

在这里我们看到判断Lod的级别主要有三个条件：
 * 1、相机视角范围，视角范围越大，所包含的tileSize就越大
 * 2、相机与瓦片距离，距离越远，所包含的tileSize也就越大
 * 3、相机视锥与瓦片是否相交
 相对应我们可以把视角剔除方法理解成以下三步：
 * 1、根据视角范围，画个大圈，把大圈里的大瓦片全加进来
 * 2、根据相机与瓦片的距离进行细化瓦片，如果太远的瓦片仍然保持很大，而近处的瓦片需要进一步更新，计算子瓦片。
 * 3、每计算一个瓦片都需要判断它是否在视锥里，否则剔除身后看不到的瓦片
 在理解以上三个条件的前提下，需要理解的是视角、距离到底与tileSize有什么关系。
 对于不同级别影像或地形，ww中设置了相应的大小，根据一定的比例系数和经验参数寻找与距离和视角的关系。即通过调整参数可控制在怎样的距离下可以看到怎样级别的地形或影像。例如在地形中设置为3.0 2.9等常数，在影像中设置为TileDrawDistance和TileDrawSpread 等常数。只是反映了一定的经验参数而已。

http://blog.csdn.net/jk276993857/article/details/5937941

1.采用dstile制作WW瓦片，需要传人参数Level0的大小和建立金字塔的等级，这是一种自顶向下的方法。顶是指金字塔的顶部。

　　创建瓦片的通用命令行格式如下：

　　dstile.exe tile --lztsd tile_size --wwcache --overviews output_directory georeferenced_image_file(s)

　　参数如下：

　　tile_size - 0图层瓦片大小，十进制数，能被180正常（最好为180/(2^N)）。

　　output_directory - 瓦片金字塔的输出路径。（如何没有盘符只要文件夹名，是在FWTool安装目录下）

　　georeferenced_image_file(s) - 影像数据的完整路径（如果影像数据在FWTool安装目录下的文件夹下，可以使用相对路径）.一系列的多重影像可以自动组合成单一的瓦片金字塔。（ A list of multiple image files can be supplied and will be automatically combined into a single set of tiles.）

举例：确定金字塔的顶部大小比如2.25度，确定了0级的瓦片数量160*160，根据地球的周长，确定了每个瓦片的长度范围，对于了256pixels的瓦片，则每个像素的分辨率可知。
　　确定了level0的tileSize，也就间接确定了Level0显示的相机（CameraBase）高度（_altitude），通过viewRange可视范围、视域体和相机中心与瓦片球面夹角确定。

　　制作WW瓦片时，确定了levlZeroTileSizeDegrees后，瓦片分多少级合适应该由影像的最高精度决定。 同样的相机高度，对于不同的瓦片数据集的levelZeroTileSizeDegrees是不一样的，所以在同样的高度存在多个瓦片等级，没有一个统一的值。

CameraBase的Update()方法中计算可视范围：

// Old view range (used in quadtile logic)
   double factor = (this._altitude) / this._worldRadius;
   if(factor > 1)
    viewRange = Angle.FromRadians(Math.PI);
   else
    viewRange = Angle.FromRadians(Math.Abs(Math.Asin((this._altitude) / this._worldRadius))*2);

　　在QuadTile和SurfaceTile的Update()方法中，一项任务就是：初始化可见瓦片的瓦片；将视点中心部分瓦片细化，外围不细化。首先初始化瓦片，

if (!isInitialized)
    {
        if (DrawArgs.Camera.ViewRange * 0.5f < Angle.FromDegrees(QuadTileSet.TileDrawDistance * tileSize)&&
            MathEngine.SphericalDistance(CenterLatitude, CenterLongitude, DrawArgs.Camera.Latitude, DrawArgs.Camera.Longitude) < Angle.FromDegrees(QuadTileSet.TileDrawSpread * tileSize * 1.25f) && DrawArgs.Camera.ViewFrustum.Intersects(BoundingBox)
           )
           Initialize();
     }

　　接着细化，对于超出可视范围的剔除：

if (isInitialized && World.Settings.VerticalExaggeration != verticalExaggeration || m_CurrentOpacity != QuadTileSet.Opacity ||
                    QuadTileSet.RenderStruts != renderStruts)
                {
                    CreateTileMesh();
                }

                if (isInitialized)
                {
                    if (DrawArgs.Camera.ViewRange < Angle.FromDegrees(QuadTileSet.TileDrawDistance * tileSize)
    && MathEngine.SphericalDistance(CenterLatitude, CenterLongitude,
                            DrawArgs.Camera.Latitude, DrawArgs.Camera.Longitude) < Angle.FromDegrees(QuadTileSet.TileDrawSpread * tileSize)
    && DrawArgs.Camera.ViewFrustum.Intersects(BoundingBox)
                        )
                    {
                        if (northEastChild == null || northWestChild == null || southEastChild == null || southWestChild == null)
                        {
                            ComputeChildren(drawArgs);
                        }

                        if (northEastChild != null)
                        {
                            northEastChild.Update(drawArgs);
                        }

                        if (northWestChild != null)
                        {
                            northWestChild.Update(drawArgs);
                        }

                        if (southEastChild != null)
                        {
                            southEastChild.Update(drawArgs);
                        }

                        if (southWestChild != null)
                        {
                            southWestChild.Update(drawArgs);
                        }
                    }
                    else
                    {
                        if (northWestChild != null)
                        {
                            northWestChild.Dispose();
                            northWestChild = null;
                        }

                        if (northEastChild != null)
                        {
                            northEastChild.Dispose();
                            northEastChild = null;
                        }

                        if (southEastChild != null)
                        {
                            southEastChild.Dispose();
                            southEastChild = null;
                        }

                        if (southWestChild != null)
                        {
                            southWestChild.Dispose();
                            southWestChild = null;
                        }
                    }
                }

                if (isInitialized)
                {
                    if (DrawArgs.Camera.ViewRange / 2 > Angle.FromDegrees(QuadTileSet.TileDrawDistance * tileSize * 1.5f)
                            || MathEngine.SphericalDistance(CenterLatitude, CenterLongitude, DrawArgs.Camera.Latitude, DrawArgs.Camera.Longitude) > Angle.FromDegrees(QuadTileSet.TileDrawSpread * tileSize * 1.5f))
                    {
                        if (Level != 0 || (Level == 0 && !QuadTileSet.AlwaysRenderBaseTiles))
                            this.Dispose();
                    }
                }

2.Globe Mapper中默认采用的方法是根据最高精度建立ww瓦片金字塔。另外一个选项就是传人Level0参数。

　　WW适用SurfaceTile、QuadTile，VirtualEarth Tile，可以加载多个不同的瓦片集合QuadTileSet。而不同的QuadTileSet有不同的levelZeroTileSizeDegrees。

　　BingMap的瓦片体系只是一种固定的瓦片层级。由于Bing地图是正射地图，所以根据屏幕像素可以计算比例尺。但是放在球上显示，需要经过纹理映射，世界坐标转换到屏幕坐标，因此与比例尺对应似乎存在问题。

  

3.WW中几个参数理解尝试：

　　如图：AB为相机海拔高度，令AB=BC做B点切线，作CE//AB交地球与E，作EG垂直AB于G，则EG=BC=AB，若认为弧BE长度等于EG，则α的弧度值=viewRange/2。事实上α的弧度值是大于viewRange/2的。

　　β为相机的FOV/2，β=22.5度。α对应弧BE，包括多少个瓦片？假设β只能看到 1个（对称分布）1/2个瓦片，则弧BE大概包括tan(PI/4)/tan(PI/8)=1/tan(PI/8)=1/0.392699=2.414个。

　　所以QuadTileSet.TileDrawDistance=3.5 ，每个Level等级加载大约5或者（3+3）个瓦片。QuadTileSet.TileDrawSpread=2.9。

　　初始化时候的范围α<QuadTileSet.TileDrawDistance*tileSize，QuadTileSet.TileDrawSpread*1.25=3.625。比细化的范围大，

　　细化范围 2α<QuadTileSet.TileDrawDistance*tileSize，细化了（1+1）或者3个瓦片；

　　超过某个范围是剔除，剔除的临界范围更大些α>QuadTileSet.TileDrawDistance*tileSize*1.5=5.25*tileSize，QuadTileSet.TileDrawSpread*1.5=4.35。

　　回到原来的问题，通过高度如何确定显示范围的金字塔等级Level。首先需要明确当前显示的范围内金字塔等级有多个，我要中心的那个。

　　我想大概中心瓦片大小sizeTileCenter=AB*tan(PI/8)/R，再根据levelZeroTileSizeDegrees求出Level。

--------------------------------------------------------------------------------------------------------------------------------

WW支持改变视域体，相当于改变了Frustrm。
R*sin（α）=Altitude
临界条件：ViewRange=2*α=2*3.5*TileSize；

地球半径R(m)	6378137
赤道周长 =2*PI*R(m)	40075016.69
Level 0瓦片 Tile0（度）	2.25	36
瓦片纹理大小（像素）	512	512
Level 0 Tile赤道分辨率（m） =2*PI*R/360*Tile0/512	489.196981	7827.152

　　以上不知道分析的对不对！