原文发布时间：2010-04-13 18:21:39
作者：毛毛虫
在做图像分割时，需要进行对图像按国家标准进行分幅裁剪。下面是我生成用来裁剪的shape文件的过程，跟大家交流下，因为要处理的数据范围较小，数据的1: 100000的地形图行列号相同，即十位分幅图名的前四位相同。
首先，读入要进行分幅的范围（此处我用的是一个shape文件，提取这个shape文件的边界范围）。根据所要生成的标准分幅的比例尺1：1万的行列间隔，也就是纬度差是2’30”、经度差是3’45”和计算行列号的公式求出shape文件所在的范围，计算范围时保证比实际的范围要大些，且生成的整个范围的边界是标准分幅的分界的经纬度上，例如求出的最大纬度是x度，求取x度与2’30”商的整数加1再乘以2’30”（保证数据在边界内，注意单位的一致）。主要用到的转换器是范围提取转换器BoundsExtractor和表达式计算转换器ExpressionEvaluator。如下图说明：
| |
使用Tiler转换器进行分割，分割后如下图：
再次提取各个多边形的边界，根据其经纬度坐标和图幅行列号计算公式，计算出每个多边形作为标准图幅的行列号，分别命名为row_name、column_name，然后使用StringConcatenator转换器与之前的计算出的前四位名称连接起来，命名为tfmc，每个多边形即命名了按国家标准分幅1：1万的分幅名称。
如果用生成的shape文件用来裁剪栅格数据，考虑到进行重投影等操作，可以使用Scaler转换器进行放大，使图幅有重叠，设置如下：
多边形放大后的shape文件如下：
放大到一个多边形，如下图：
| |
这样就能保证用这些多边形来裁剪影像是重叠部分了。另外如果数据范围较大，数据的1: 100000的地形图行列号不同，可以用Tester进行判读来区分图幅名称的前四位。