第8页
第4排
你好,
我有来自Postgis数据库和光栅文件的多边形。我想使用裁剪器获取部分光栅文件,该文件位于用作裁剪器的多边形内。
这是剪子:
It should be a 10x10 km box.
请查看最小范围:520000.0,616万.0
这是剪报(光栅文件)
现在,我对clipper输出(内部)感到非常惊讶:
我不知道为什么。3出现在小数点后?
你知道我做错了什么吗?
我的期望是从剪子的边界开始,正好是6160000.0。
你没有做错什么,但是当剪切光栅时,FME会返回整个像素(这很有意义,否则剪切将涉及重新采样)。
@红色地理谢谢你的回答,但我有一个问题:
我还使用了其他工具,如这里介绍的:
https://gis.stackexchange.com/questions/200407/merging-gdal-tiles-and-crop-via-bounding-box
我使用的gdalwarp命令是:
gdalwarp.exe-of gtiff-co compress=lzw-cutline clipper.kml-crop-to-cutline tifffilestobeclippedandmerged.vrt result.t
以这种方式生成的文件的gdal信息是:
Origin = (540000.000000000116415,6170000.000000001862645) Pixel Size = (0.699986000280060,-0.699986000280581) Metadata: AREA_OR_POINT=Area Image Structure Metadata: COMPRESSION=LZW INTERLEAVE=PIXEL ________________________.000,6170000.000)(9d38'9.67“E,55D40'27.01“N)左下(540000.000,6160000.000)(9d38'4.44“E,55d35'3.56“N)右上(550000.000,6170000.000)(9D47'42.04“E,55D40'23.68“N)右下(550000.000,6160000.000)(9D47'35.50“E,55d35'0.23“N)中心(545000.000,6165000.000)(9D42'52.91“E,55d37'43.71"N)
而对于使用上述FME处理的光栅,它是:
原点=(540000.0000039999553346169999.80003999998186993)像素大小=(0.7000000000000,-0.699951000979980)元数据:区域或点=区域图像结构元数据:压缩=LZW交错=像素角坐标:左上(540000.000,6169999.800)(9d38'9.67“E,55D40'27.01“N)左下(540000.000,6160000.300)(9d38'4.44“E,55d35' 3.57"N)Upper Right ( 550000.200,6169999.800)(9D47'42.05“E,55D40'23.67“N)右下(550000.200,6160000.300)(9D47'35.51“E,55d35'0.24“N)中心(545000.100,6165000.050)(9D42'52.92“E,55D37'43.72“N)
所以我想知道在FME中我能做些什么来获得和使用GDAL相同的结果?
从我收集的gdalwarp不仅剪辑,而且重新采样图像。您可以通过检查原始光栅上的信息(特别是像素大小)来确认这一点。
基本上,有两种方法可以处理剪裁光栅,我们在这里看到:
FME方法的好处是,每个生成的像素仍然与原始像素匹配。用GDAL的方式,我不太确定。
如果你想让fme以gdal的方式运行,你必须计算出当前的像素大小,然后重新对光栅取样,以便在剪切之前,它的新像素与剪辑多边形完全对齐。说实话,我现在真的说不出该怎么做:)(但我有时差反应)
?2019安全亚搏在线软件公司|法律