worldwind的技术手册.docVIP

worldwind技术手册 目录①首先从from the Global Land Coverage Facility下载Landsat7的全球卫片 Remarks: The scenes come as a set of grayscale GeoTiff files for each row/path/date scene ②然后使用波段1 (Blue),2 (Green),3 (Red)合成RGB的GeoTiff image。 Remarks: 由于我不不想编程实现，我使用了PixelSense，但是PixelSense 提供了一个 linear clip contrast stretch，因此必须 replace the image data after PixelSense outputted the RGB GeoTiff。 ③再对RGB Landsat GeoTiff使用颜色增强算法 Remarks: The basic idea for the algorithm I used is to divide the image into a grid with each grid square being n by n pixels, where n is about 25 for a 30meter/pixel Landsat scene. For each square, find the least-squares for both the high and low, and then save this information in order to later apply a linear stretch of the pixels in that square, weighting the R, G, B bands individually to the desired tastes. I found that generally, the source Landsat scenes have too much blue and red, and I weight the stretch to compensate for this for a more pleasing image. Next, use a cubic convolution to smooth out the stretch information in the grid squares. This is only applied to the information in the grid squares, not the actual pixel data. I did this 3 times for increased smoothness, but it seems that some areas still needed more smoothing. Next, stretch the pixels in each square according to the new grid stretch information. I tried to smooth the differences between the squares again by using bi-linear averaging (I think thats what to call it) for the grid squares that surround the current square being processed. ④ 重新将图像从UTM 投影到大地投影（Geodetic projection）。 Remarks: I used GDAL for this. Its free, fast, and gives decent results. ⑤镶嵌和重复将Landsat影像切割成规则的大地网格方块（Geodetic grid squares）。 Remarks: I used ERMapper for this, with the ESG Utilities add-on. 因为Landsat图像之间互相交迭；因此在切割瓦片前必须将它们缝合。I used grid squares that were 2.25 x 2.25 degrees, with a resulting image size of 8152 x 8152 pixels, which roughly equates to 30 m/pixel at the equator. ⑥建立图像瓦片式金字塔 Remarks: I wrote a custom prog