斯潘8
斯潘4
Start:教程:在FME中读取和转换卫星图像|下一步:光栅计算和光栅选项板
This is the third article in the教程:在FME中读取和转换卫星图像.This article will run you through the process of creating an RGB image from single band raster images collected by Sentinel-2 with FME.还演示了一个由FME云提供的Web服务,它非常适合预览来自陆地卫星和哨兵来源的图像。我们还将展示如何利用AWS托管的Mapzen全球数字高程模型(DEM),快速轻松地创建覆盖卫星图像的三维曲面。
Complete workspace (without source data):geotiff2obj.fmwt
1) Using FME,我们可以构建一个为我们进行云/亚搏在线数据过滤的工作流。Go to Landsat/Sentinel selector page:http://demos.fmeserver.com/landsat-sentinel-selector/.使用该工具可以从您选择的位置查找和下载图像。在本练习中,仅选择Sentinel-2提供者。
First select location.
尝试选择有一些绿色植被的图像(否则NDVI将不可见)。
Then select an image with green vegetation.
2)下载波段2,3,4 (BGR) and 8 (NIR) AND their corresponding previews of the same bands for the same tile.在下图中,这包括下载常规的高分辨率b02.jp2,b03.jp2型,B04.jp2,B08.jp2,从预览文件夹,the previews of B02.jp2,b03.jp2型,B04.jp2,b08.jp2.
查找要下载的文件:上面的工作流程应该可亚搏在线以帮助您找到我们要下载的Sentinel-2图像。Here you can also see the Tile name from the online folder/URL path.在这种情况下,磁贴名称为“10uev”。
3)使用图块名称创建文件夹,把所有的高分辨率波段文件放进去。The tile name is located at the top of the Sentinel public image browser,在上面的图像中,10>u>ev导致文件夹名为10uev。
4)在tile文件夹中创建一个文件夹,call it "Preview",将所有预览放入该文件夹(我们将在构建工作区时使用预览进行测试)。
5)同时,在下载页面上,打开JSON元数据文件,查看云和数据覆盖率。
6)使用jpeg 2000(geojp2/gmljp2)格式,使用FME数据检查器检查一些下载的文件。请注意像素值和解释(灰色16)。您还将注意到这些图像大部分是黑色的。If you look closely,however,you may notice that there are come areas which appear a little lighter.本文的下一部分将对如何“照亮”这些进行说明。
NOTE:我们还可以使用AWS阅读器上的Sentinel-2来收集这些图像,而不是本手册的过程。however,这样做可以更深入地了解图像是如何在线存储的。
现在我们已经将图像下载到您的机器上,我们可以构建一个工作流程,将波段组合在一起并“照亮”图像,使其看起来更像我们应该期待的。亚搏在线
1) Open Workbench,add a JPEG2000 reader,读取所有下载的预览文件。确保功能类型是由文件名构成的。预览允许在创建和测试阶段进行非常快速的处理。对于最终图像,只应调整源文件的路径。
2)当四种功能类型出现在画布上时(带区2,3,4,8)添加说明它们代表什么的注释(波段2-蓝色,Band 3 - Green,Band 4 - Red,波段8-红外)-这样我们总能看到,我们目前使用的波段。
3)让我们做个实验,我们该怎么处理乐队?我们可以先添加一个Rasterbandcombiner。将所有RGB波段连接到它,运行并检查结果,can you tell what's wrong?
4)在组合之前,让我们先用RasterbandInterpretation胁迫器更改每个波段的解释(蓝色到蓝色16,绿色到绿色16,红色到红色16)运行并检查,你能告诉我这次怎么了吗?
5) We can change the brightness with RasterExpressionEvaluator.在RasterBandCombiner之后添加RasterExpressionEvaluator。
Create a common multiplier published parameter,MULTIPLIER,set by default to 20:
Multiplier parameter settings.
为每个波段添加以下表达式:
解释 | Expression |
Red8 | A[0]*$(乘数)/255 |
Green8 | A[1]*$(MULTIPLIER)/255 |
Blue8 | A[2]*$(MULTIPLIER)/255 |
RASTERExpressionEvaluator公司:这里我们增加了像素的亮度,并将每个波段的解释从16位更改为8位。The MULTIPLIER determines how much brighter to make the pixels.We are dividing by 255 to get the pixel values into the 8-bit range.8-bit color bands can have pixel values up to 255 (or 2^8-1) whereas 16-bit color bands can have values up to 65535 (or 2^16-1).
6)重新运行,检查,更好?也许不是。
7) Add an AttributeCreator after each of the RasterBandInterpretationCoercer's and create a numeric sort attribute,for example _sort with a value of 1 for Red,绿色2个,蓝色3个。把它们都放进分梳机前的分梳机。Make sure Red has the lowest value,Green - second lowest,蓝色-最高。
分拣机参数
8)运行并检查。现在输出应该看起来不错。尝试调整乘数以获得更好的结果。真正的色彩平衡需要更多,本课程不涉及。
9) What if we want to make an infrared image?Connect Infrared to red stream,从红到绿,从绿到蓝,disable the old connections and re-run.在这里,我们将所有的数据在可见光谱中移动得更高。If you have picked an image with a lot of vegetation you should see quite a ‘red' image as trees reflect a lot of infrared light.
InfraRed satellite image
10)添加geotiff writer,将“光栅文件定义”设置为“自动”,并输入光栅文件名的光栅图块。运行翻译以生成本地RGB图像-我们可以使用它进行进一步处理或保留当前工作流。亚搏在线
11)(可选)切换到高分辨率光栅(通过删除“预览”部分调整文件路径)以生成高分辨率图像。For further experiments,确保切换回预览。
Note for coordinate system information,必须使用高分辨率光栅图像。
在本练习中,我们将利用一个定制的变换器,该变换器可以下拉并剪切mapzen的全局DEM的一个子集。您可以阅读有关服务的更多信息在这里.然后,我们将使用DEM创建一个三维曲面,并将先前创建的彩色卫星图像作为外观添加到曲面中。
1)在RasteExpressionEvaluator之后,将MapzenawsdemDownloader从FME Hub添加到现有工作区。重要的是,您使用的图像有一个坐标系,否则工作区将失败。请注意,高分辨率图像标记有坐标系信息,而预览没有。
2)设置一些低分辨率(缩放级别8)并运行和检查。
3)将缩放级别增加到10。您可能需要增加“要下载的最大瓷砖数量”的数量。将输出DEM保存为geotiff格式(或者,如果您的连接良好且不关心流量,只需使用变压器)。
4) Add a TINGenerator,并将mapzenawsdem输出端口连接到点/线输入端口。将表面公差设置为DEM光栅分辨率(暗示: look at the zoom level in the custom transformer).
5) Add a Scaler to exaggerate the vertical component of the TIN Surface.Scale 'z' by an appropriate value (depends on location) - for BC/California values 2.5-4 work well,for flat areas - 10 or even more.
6)增加外观设计,send the RGB raster (from the RasterExpressionEvaluator output port) to Appearance input port,三角网表面到几何输入端口。Under Texture Coordinate Generation Parameters,将纹理映射类型设置为“从俯视图”
外观参数
7)运行。在DI中检查结果
8) (Optional) Switch to hi-res to get the best quality final product in DI
9)添加obj writer,保存曲面。
10)(可选)在Meshlab中检查结果。看起来meshlab不喜欢大于4096*4096的纹理-必要时在FME中应用纹理之前重新取样-如果没有,可以免费下载meshlab在这里.
最终工作区应如下所示:
?2019安全亚搏在线软件公司|Legal