遥感数据即得即用（Ready To Use，RTU）地理格网产品规范
| 遥感数据即得即用（Ready To Use，RTU）地理格网产品规范 作者：焦伟利 龙腾飞 何国金 王威 张兆明 王桂周 彭燕 2020年5月29日 |
 焦伟利, 龙腾飞, 何国金, 等. 遥感数据即得即用（Ready To Use，RTU）地理格网产品规范[J/OL]. 中国科学数据, 2020. (2020-05-26). DOI: 10.11922/csdata.2020.0028.zh. |
摘要&关键词
[编辑]摘要：遥感影像是一种平面栅格数据，以一定的地图投影方式反映地球表面的表象。常规的标准遥感影像分幅产品是以一定的规则按照景来分幅的，不同的卫星影像，景的分幅和编码规则不同。遥感影像的景没有完全与地球上的地理坐标对应。为了使长时序影像更容易查询、分析和管理，本文提出遥感影像即得即用（Ready To Use，RTU）地理格网产品，即将全球影像按照定义的格网系统剖分为具有统一规则和属性的数据格网。规定了遥感影像按照经纬坐标格网和直角坐标格网的分幅与编码方法和RTU地理格网产品规范，适用于不同空间分辨率的遥感影像按照不同投影地理格网的分幅与编码和RTU地理格网产品的生产。RTU地理格网产品具有使用灵活、格网编码与空间地理坐标相一致，可以按照空间维和时间维建立数据立方体，形成空间和时间上的连续产品，便于用户直接进行空间分析和应用。本规范为遥感影像共享和遥感信息整合提供以地理格网为单元的空间参照，可用于长时序、不同空间分辨率的遥感影像的综合分析和应用。
关键词：遥感数据；即得即用（Ready To Use，RTU）产品；地理格网；分幅；编码
Abstract & Keywords
[编辑]Abstract: Remote sensing image is a kind of plane grid data, which reflects the surface of the earth in a certain map projection way. Normal standard remote sensing data are divided into scenes according to certain rules. Different satellite data have different framing and coding rules. The scene of remote sensing data does not correspond to the geographical coordinates of the earth. The Ready To Use (RTU) geogrid product of remote sensing data is proposed in this paper in order to make the long-time sequence image easier to query, analyze and manage. That is, the global image is divided into data grids with resolution, corner coordinates, grid size, product type and other attributes according to the defined grid system. This specification specifies the division and encoding methods of remote sensing image according to longitude and latitude grid and Cartesian coordinate grid, and RTU geogrid product, which are applicable to different resolution remote sensing images, according to the different projection geographically dissecting and encoding, and the production of RTU geogrid products. RTU geogrid products are flexible in use, and the geographical coordinates is consistent with geocoding. Data cubes can be built according to spatial and timely dimensions to form continuous products in space and time, which is convenient for users to directly carry out spatial analysis and application. This specification provides spatial reference based on geographic grid for remote sensing image sharing and remote sensing information integration, which can be used for comprehensive analysis and application of remote sensing images with long time sequence and different spatial resolution.
Keywords: remote sensing data; RTU (Ready To Use) product; geographic grid; gridding; encoding
1 范围
[编辑]本标准规定了多种地图投影的遥感影像的地理格网的划分与编码方法和RTU地理格网产品规范。
本标准适用于遥感影像产品的地理格网分幅与编码，为遥感影像共享和遥感信息整合提供以格网为单元的空间参照，并可用于遥感影像地理格网产品生产。
2 术语与定义
[编辑]2.1 格网 grid
[编辑]由两组或多组曲线集所组成的网络，曲线集合中的曲线按某种算法相交。
注：曲线集把空间分割成格网单元[1]。
2.2 格网单元 grid cell
[编辑]构成格网系统中某级格网的基本单位[2]。
2.3 地理格网 geographic grid
[编辑]按照一定的数学规则对地球表面进行划分而形成的格网[2]。
2.4 经纬坐标格网 geographical graticule
[编辑]按一定经纬度间隔对地球表面进行划分而形成的格网[2]。
2.5 直角坐标格网 rectangular grid
[编辑]将地球表面区域按数学法则投影到平面上，按一定的纵横坐标间距和统一的坐标原点对地表区域进行划分而构成的格网[2]。
2.6 格网编码 grid encoding
[编辑]按照一定规则，赋予格网单元惟一标识代码的过程[2]。
2.7 即得即用产品ready to use product（RTU）
[编辑]按一定的标准处理的，便于用户直接分析和应用的一系列高级遥感数据产品。
2.8 数据格网 tile
[编辑]具体影像数据的剖分格网。
2.9 RTU地理格网产品RTU geogrid product
[编辑]按照一定空间基准和地理格网对遥感数据进行分幅和编码，形成具有统一规则和属性的系列遥感即得即用格网产品。这类产品可直接组成时空数据立方体，便于多源、多尺度、长时序遥感数据的高效存取和时空分析。
2.10 质量标识 quality attribute（QA）
[编辑]随数据产品一起提供的数据质量标识信息，以波段的形式存在。可用于确定产品的每个像素是否适合某种用途或应用。例如像元质量缺陷或无效像素、云、雪等标识。
3 坐标系统及地图投影
[编辑]3.1 坐标系
[编辑]采用2000国家大地坐标系（CGCS2000），比例尺小于1:10 000的影像也可以采用WGS84坐标系。
3.2 地图投影
[编辑]小比例尺大尺度（全球或全国）影像建议采用经纬度坐标。平面直角坐标投影依据尺度及地面分辨率采用不同地图投影，具体建议如下：
①全球尺度投影建议采用正弦曲线等面积伪圆柱投影（Sinusoidal），中央经线为0°。
②中国区域尺度投影建议采用阿尔伯斯等面积割圆锥投影（Albers），中央经线105°，两条标准纬线为25°和47°。
③地面分辨率优于100 m或东西跨度小于1000 km的影像建议采用UTM投影。
④极地区域建议采用兰伯特等积方位投影（Lambert Azimuthal Equal-Area）。
4 经纬坐标格网
[编辑]4.1 经纬坐标格网概述
[编辑]经纬坐标格网面向大范围（全球或全国），适用于较概略表示信息的分布和粗略定位的应用。经纬度格网采用赤道和本初子午线的交点为原点。经纬度坐标格式按照度、分、秒格式表达。经纬坐标格网按经差、纬差分级。代码由格网间隔代码、南北半球代码、纬度代码、东西半球代码、经度代码组成。图1为全球10°×10°经纬度格网分幅与编码示意图。
图1 全球10°×10°经纬度格网分幅与编码
4.2 经纬坐标格网分级
[编辑]分级规则：各层级的格网间隔为整数倍关系，同级格网单元的经差、纬差间隔相同。经纬坐标格网基本分为6级。
4.3 经纬坐标格网编码
[编辑]经纬坐标格网编码如下：
格网间隔代码＋南北半球代码＋纬度代码＋东西半球代码＋经度代码。
格网间隔代码依据表1所列的格网间隔用2位数字码表示，不足2位的在前面补0。如格网间隔1′时表示为01，格网间隔10″时，表示为10。间隔单位代码用D表示以度为单位、M表示以分为单位、S表示以秒为单位。
表1 经纬坐标格网分级
| 格网间隔 | 10°×10° | 1°×1° | 10′×10′ | 1′×1′ | 10″×10″ | 1″×1″ |
| 格网名称 | 十度 | 一度 | 十分 | 一分 | 十秒 | 一秒 |
| 格网编码 | N* E** | N** E*** | N*** E**** | N**** E***** | N***** E****** | N****** E******* |
南北半球代码用1位字母码表示，用N表示北半球，S表示南半球。
纬度代码数值和长度依据表1所列的格网间隔确定，数字码为纬度值除以格网间隔值向下取整。如格网间隔1°时，纬度代码为纬度值向下取整；格网间隔10’时，纬度代码为在1°代码基础上添加10′的数字码，10′码为纬度的分除以10后向下取整，依此类推。具体公式为：int（|纬度/格网间隔|）。
东西半球代码用1位字母码表示，用E表示东半球，W表示西半球。
经度代码比纬度代码多一位，编码方式与纬度代码相同。即int（|经度/格网间隔|）。
例：某点位于75°41′15″N，143°02′35″E，其
10°×10° 十度格网代码为：10DN7E14
1°×1° 一度格网代码为：01DN75E143
10′×10′ 十分格网代码为：10MN754E1430
1′×1′ 一分格网代码为：01MN7541E14302
10″×10″ 十秒格网代码为：10SN75411E143023
1″×1″ 一秒格网代码为：01SN754115E1430235
5 平面直角坐标格网
[编辑]5.1 平面直角坐标格网概述
[编辑]直角坐标格网与所采用的地图投影密切相关，投影相同，格网的平面坐标和编码相同，投影不同，格网的平面坐标和编码不同。
正弦曲线投影格网： 小比例尺全球尺度影像建议采用正弦曲线等面积伪圆柱投影（Sinusoidal），原点为赤道和本初子午线的交点。格网分幅与编码示意如图2所示。
图2 正弦曲线投影千公里格网示意图
阿尔伯斯投影格网： 中国区域尺度影像建议采用阿尔伯斯等面积割圆锥投影（Albers），原点为赤道和105°中央经线交点，两条标准纬线为25°和47°。格网分幅与编码示意如图3所示。
图3 阿尔伯斯投影千公里格网示意图
UTM投影格网： 地面分辨率优于100 m的影像建议采用6°分带的UTM投影。格网分幅与编码示意如图4所示。
图4 UTM百公里格网示意图
兰伯特等积方位投影格网 ：极地区域建议采用兰伯特等积方位投影（Lambert Azimuthal Equal-Area），原点为南极点或北极点。格网分幅与编码示意如图5所示。
图5 北极兰伯特等积方位投影百公里格网示意图
直角坐标格网代码主要由投影代码、南北半球代码、纵坐标格网代码、东西坐标代码（或投影带号代码）、横坐标格网代码组成。
5.2 直角坐标格网分级
[编辑]分级规则：各级格网的间隔为整数倍关系，同级格网单元在X、Y方向的间距相等。
直角坐标格网系统根据格网单元间隔分为6级，以千公里格网单元为基础，按10倍的关系细分，如表2所示。
表2 直角坐标格网系统分级
| 格网间隔（单位： m） | 1 000 000 | 100 000 | 10 000 | 1 000 | 100 | 10 | 1 |
| 格网名称 | 千公里格网 | 百公里格网 | 十公里格网 | 公里格网 | 百米格网 | 十米格网 | 米格网 |
| 格网名称代码 | A | B | C | D | F | G | H |
5.3 直角坐标格网编码方法
[编辑]平面直角坐标格网代码组成如下：
投影代码+南北坐标代码＋纵坐标格网代码＋东西坐标代码（或投影带号代码）＋横坐标格网代码。
投影代码为1位数字码。1-正弦曲线等面积伪圆柱投影（Sinusoidal），2-阿尔伯斯等面积割圆锥投影（Albers），3-UTM投影，4-北极兰伯特等积方位投影（Lambert Azimuthal Equal-Area），5-南极兰伯特等积方位投影（Lambert Azimuthal Equal-Area）。
南、北坐标代码采用1位字母码。原点以南坐标用“S”表示，原点以北坐标用“N”表示。
东、西坐标代码采用1位字母码。原点以东坐标用“E”表示，原点以西坐标用“W”表示。东西坐标代码不适用于有坐标平移的投影。
投影带号代码用于UTM投影，采用2位数字码表示。UTM投影采用6°分带，全球共分60带，投影带号代码为01-60。
纵坐标格网代码与横坐标格网代码为选用的层级格网间隔字位数值向下取整。即：int（|坐标值/间隔值|）。具体为：
千公里格网代码由坐标值千公里字位数值向下取整构成。
百公里格网代码由坐标值百公里字位数值向下取整构成。
十公里格网代码由坐标值十公里字位数值向下取整构成。
公里格网代码由坐标值一公里字位数值向下取整构成。
百米格网代码由坐标值百米字位数值向下取整构成。
十米格网代码由坐标值十米字位数值向下取整构成。
米格网代码由坐标值一米字位数值向下取整构成。
例：某点位于39°55′N, 116°30′E，其正弦曲线投影横坐标值为9960467.2 m，纵坐标值为 4420276.2 m；阿尔伯斯投影（中央经线为105°，双标准纬线为25°、47°）横坐标值为964574.1 m，纵坐标值为4346377.1 m；UTM投影6°分带带号为50，横坐标值为 457268.2 m，纵坐标值为4418627.7 m。不同投影格网代码如表3所示。
表3 直角坐标格网编码
| 格网名称（编号） | 千公里格网（A） | 百公里格网（B） | 十公里格网（C） | 公里格网（D） | 百米格网（F） | 十米格网（G） | 米格网（H） |
| 格网间隔（单位： m） | 1000000 | 100000 | 10000 | 1000 | 100 | 10 | 1 |
| 正弦曲线（1） | 1AN4E9 | 1BN44E99 | 1CN442E996 | 1DN4420E9960 | 1FN44202E99604 | 1GN442027E996046 | H1N4420276E9960467 |
| 阿尔伯斯（2） | / | 2BN43E9 | 2CN434E96 | 2DN4346E964 | 2FN43463E9645 | 2GN434637E96457 | 2HN4346377E964574 |
| UTM（3） | / | 3BN44504 | 3CN4415045 | 3DN441850457 | 3FN44186504572 | 3GN4418625045726 | 3HN441862750457268 |
例：某点位于80°30′N，116°30′E，其兰伯特等积方位投影横坐标值为948433.8 m，纵坐标值为472871.6 m。格网代码如表4所示。
例：某点位于80°30′S，116°30′W，其兰伯特等积方位投影横坐标值为-948433.8 m，纵坐标值为-472871.6 m。格网代码如表4所示。
表4 极地区域坐标格网编码
| 格网名称（编号） | 千公里格网（A） | 百公里格网（B） | 十公里格网（C） | 公里格网（D） | 百米格网（F） | 十米格网（G） | 米格网（H） |
| 格网间隔(单位：m) | 1000000 | 100000 | 10000 | 1000 | 100 | 10 | 1 |
| 北极兰伯特（4） | / | 4BN4E9 | 4CN47E94 | 4DN472E948 | 4FN4728E9484 | 4GN47287E94843 | 4HN472871E948433 |
| 南极兰伯特（5） | / | 5BS4W9 | 5CS47W94 | 5DS472W948 | 5FS4728W9484 | 5GS47287W94843 | 5HS472871W948433 |
5.4 直角坐标格网扩充
[编辑]直角坐标格网可在表2给出的直角坐标格网分级基础上按整数倍关系向小于1 m的格网单元扩展。分米格网名称代码为I。
例： 某点位于(39°55′N, 116°30′E)，其UTM投影带号为50，横坐标值为457268.2 m，纵坐标值4418627.7 m，求其UTM分米格网的代码。
分米格网代码：3IN44186277504572682
直角坐标格网亦可在表2给出的直角坐标格网分级基础上，对现有层级再细分。各层级格网单元间隔宜根据如下公式计算得到的值，见表5。
格网间隔 = 2 × 10n或5 × 10n，n∈{0，1，2，3，4，5}
表5 直角坐标格网扩充分级
| 格网名称 | 千公里区间格网 | 百公里区间格网 | 十公里区间格网 | 公里区间格网 | 百米区间格网 | 十米区间格网 |
| 格网名称代码 | A# | B# | C# | D# | F# | G# |
| 指数n取值 | 5 | 4 | 3 | 2 | 1 | 0 |
| 2 × 10n格网间隔（单位： m） | 200000 | 20000 | 2000 | 200 | 20 | 2 |
| 5 × 10n格网间隔（单位： m） | 500000 | 50000 | 5000 | 500 | 50 | 5 |
注：表中“#”为依据区间格网间隔取2或5。
扩充后的格网编码形式为：投影代码+格网名称代码+直角坐标格网代码，直角坐标格网代码公式为：int（|坐标值/间隔值|）。
例：某点位于(39°55′N，116°30′E)，其UTM投影带号为50，横坐标值为457268.2 m，纵坐标值4418627.7 m，UTM扩充的格网间隔代码见表6。
表6 UTM坐标扩充格网编码
| 格网名称 | 区间格网间隔 （单位： m） | UTM（3） |
| 千公里区间格网（A#） | 500000 | / |
| 百公里区间格网（B#） | 20000 | 3B2N2205022 |
| 十公里区间格网（C#） | 5000 | 3C5N8835091 |
| 公里区间格网（D#） | 200 | 3D2N22093502286 |
| 百米区间格网（F#） | 50 | 3F5N88372509145 |
| 十米区间格网（G#） | 2 | 3G2N220931350228634 |
6 遥感影像产品的地理格网分幅与编码
[编辑]6.1 影像格网编码
[编辑]影像格网的编码由地理格网编码、卫星标识、传感器标识、影像获取时间、影像产品版本、影像产品名称、影像波段标识7部分组成。各部分以“_”分隔，各部分说明见表7。
表7 影像格网编码各部分说明
| 名称 | 要求说明 |
| 地理格网编码（grid_id） | 符合第4、第5节要求 |
| 卫星标识（satellite_id） | 参见[3] 附录B 。例如，Landsat系列分别用L3/L4/L5/L7/L8标识 |
| 传感器标识（sensor_id） | 参见[3] 附录B 。例如，Landsat系列分别为MSS、TM、ETM、OLI、TIR |
| 获取时间（acquisition_date） | 图像获取时间，按年year (YYYY)月month (MM)日day (DD)格式 |
| 产品名称（product_name） | RTU产品名称缩略表（见表8） |
| 波段标识（band_id） | 波段号。一般为数字，质量标识波段用QA标识 |
| 版本（Version_id） | 产品生产的版本号。V后面跟2位数据表示 |
| 扩展名（extension） | 元数据与数据文件分别为xml、tif格式 |
表8 RTU产品名称缩略表
| 中文 | 英文 | 缩写 |
| 正射影像 | Digital Orthophoto Map | DOM |
| 全球一张图 | Globe Map | GlobeM |
| 全国一张图 | China Map | ChinaM |
| 区域影像图 | Regional name Image Map | 区域缩写+M |
| 星上反射率 | Top Of Atmosphere Reflectance | TOA |
| 地表反射率 | Land Surface Reflectance | LSR |
| 星上亮度温度 | Top Of Atmosphere Brightness Temperature | BT |
| 地表温度 | Land Surface Temperature | LST |
| 归一化差值植被指数 | Normalized Difference Vegetation Index； | NDVI |
| 增强植被指数 | Enhanced Vegetation Index | EVI |
| 土壤调节植被指数 | Soil Adjusted Vegetation Index | SAVI |
| 修正的土壤调节植被指数 | Modified Soil Adjusted Vegetation Index | MSAVI |
| 归一化燃烧指数 | Normalized Burnt Ratio | NBR |
| 归一化差值水体指数 | Normalized Difference Water Index | NDWI |
| 归一化差值湿度指数 | Normalized Difference Moisture Index | NDMI |
| 像元质量标识 | Pixel Quality Attribute | QA |
6.2 几种常用的遥感影像产品地理格网分幅与编码
[编辑]6.2.1 UTM投影的遥感影像产品地理格网
[编辑]影像产品格网采用UTM投影和6°分带。图6为全球UTM百公里格网分幅示意图。
图6 全球UTM百公里格网分幅示意图
对于跨带影像，跨带部分需要重投影到相应投影带，并按照该投影带的格网进行分幅。
影像可以覆盖整个格网或部分格网。部分覆盖的格网一般是位于一景影像的边缘。
地面分辨率为10-100 m的影像建议采用百公里格网分幅，地面分辨率优于10 m的影像建议采用十公里格网分幅。编码规则为：投影代码+南北坐标代码＋纵坐标格网代码＋东西坐标代码（或投影带号代码）＋横坐标格网代码。
例：某点位于39°55′N，116°30′E，其UTM投影6°分带带号为50，横坐标值为457268.2 m，纵坐标值为4418627.7 m。该点在2016年10月18日的Landsat8地表温度百公里格网产品编码为：
6.2.2 全球一张图地理格网
[编辑]全球一张图产品采用Landsat系列卫星、Sentinel系列卫星等中分辨率卫星数据制作，采用WGS84坐标系统经纬度投影（EPSG:4326）和0.00025°分辨率分块输出，输出格式为GeoTIFF。
全球一张图产品采用10°×10°经纬度格网分幅（如图7所示），编码规则为：
格网间隔代码＋南北半球代码＋纬度代码＋东西半球代码＋经度代码
图7 全球一张图产品地理格网分幅示意图
具体分幅编码示例如表9所示。
表9 全球一张图产品分幅编码示意
| 10DN8W18 | 10DN8W17 | …… | 10DN8E00 | 10DN8E01 | …… | 10DN8E17 |
| 10DN7W18 | 10DN7W17 | …… | 10DN7E00 | 10DN7E01 | …… | 10DN7E17 |
| …… | …… | …… | …… | …… | …… | …… |
| 10DN0W18 | 10DN0W17 | …… | 10DN0E00 | 10DN0E01 | …… | 10DN0E17 |
| 10DS1W18 | 10DS1W17 | …… | 10DS1W00 | 10DS1W01 | …… | 10DS1W17 |
| …… | …… | …… | …… | …… | …… | …… |
| 10DS5W18 | 10DS5W17 | …… | 10DS5W00 | 10DS5W01 | …… | 10DS5W17 |
6.2.3 全国一张图地理格网
[编辑]全国一张图采用高分一号、高分二号、资源三号系列等高分辨率卫星数据制作，采用WGS84坐标系统阿尔伯斯等面积割圆锥投影（中央经线105°，两条标准纬线为25°和47°）和2 m分辨率分块输出，输出格式为GeoTIFF。
全国一张图产品采用100公里×100公里的直角格网分幅，编码规则为：
投影代码+南北坐标代码＋纵坐标格网代码＋东西坐标代码＋横坐标格网代码
具体分幅编码示例如表10所示。
表10 全国一张图产品分幅编码示意
| 2BN60W26 | 2BN60W25 | …… | 2BN60E00 | 2BN60E01 | …… | 2BN60E23 |
| 2BN59W26 | 2BN59W25 | …… | 2BN59E00 | 2BN59E01 | …… | 2BN59E23 |
| …… | …… | …… | …… | …… | …… | …… |
| 2BN04W26 | 2BN04W25 | …… | 2BN04E00 | 2BN04E01 | …… | 2BN04E23 |
7 RTU地理格网产品
[编辑]7.1 RTU地理格网产品的构成及格式
[编辑]RTU地理格网产品由产品元数据文件、图像数据文件、缩略图文件、拇指图文件（可选）和像元质量标识文件构成。
产品元数据文件： 采用XML格式，是产品的元数据描述文件。
产品图像数据文件： 采用GeoTIFF格式，是图像实体数据。
产品缩略图文件： 采用JPEG格式，缩略图的较长边大小为1000–1200像素之间，覆盖全部图像，长宽比例与原图像一致。
产品拇指图文件（可选）： 采用JPEG格式，拇指图较长边大小为200–300像素之间，覆盖全部图像，长宽比例与原图像一致。
像元质量标识文件： 采用GeoTIFF格式，是产品的像元质量标识文件，标识图像数据质量，包括无效数据、云、云的阴影、雪、水等信息。
7.2 RTU地理格网产品的命名规则
[编辑]影像格网数据产品的名称应由产品名称和扩展名两部分组成，其表现形式如下：
影像格网的编码_起始时间_结束时间_版本号_产品名称.扩展名
产品名称要求与影像格网编码相同。
示例1：2015年全国一张图产品中格网编码为2BN59E01的文件名称如下：
2BN59E01_20150101_20151231_V01_ChinaM.XML
2BN59E01_20150101_20151231_V01_ChinaM.TIF
2BN59E01_20150101_20151231_V01_ChinaM.JPG
2BN59E01_20150101_20151231_V01_ChinaM_THUMB.JPG（可选）
2BN59E01_20150101_20151231_V01_ChinaM_PIXEL-QA.TIF
7.3 RTU地理格网产品元数据
[编辑]RTU地理格网产品元数据在用于分幅影像元数据的基础上增加有关地理格网信息。示例内容见表11。
表11 RTU地理格网产品元数据定义及示例
| 名称 | 域 | 说明 | |
| 0 | rtu_metadata | N/A | RTU级元数据 |
| 1 | tile_metadata | N/A | 格网元数据 |
| 1.1 | tile_id | 地理格网编码_卫星标识_传感器标识_获取时间_版本号_产品名称NXXXXXX_LX_XX_yyyymmdd_VVV_product | 格网标识或者格网文件名，由地理格网编码、卫星标识、传感器标识、获取时间、版本号、产品名称组成例: 3BN44504_L5_TM_20151126_V01_LSR |
| 1.2 | rtu_version | VVV | RTU产品版本号 |
| 1.3 | tile_production_date | YYYY-MM-DD | 格网生成日期 |
| 1.4 | bounding_coordinates | west, east (degrees; -180 to 180)north, south (degrees; -90 to 90) | 地理坐标范围 |
| 1.5 | spatial_reference-informmation | N/A | 有关地理格网产品空间参照信息类 |
| 1.5.1 | datum | WGS84 | 坐标系 |
| 1.5.2 | projection | UTM | 投影标识 |
| 1..5.3 | units | meters | 投影坐标单位 |
| 1.5.4 | corner_point_ x y ,location | (Variable) | 格网角点坐标（左上或右下） |
| 1.5.5 | grid_origin | corner | 格网坐标原点 |
| 1.5.6 | UTM_proj_params | N/A | UTM投影参数集 |
| 1.5.6.1 | zone | 50 | 带号 |
| 1.5.6.2 | central_meridian | 117 | 中央经线经度 |
| 1.5.6.3 | origin_latitude | 0.0 | 投影原点纬度 |
| 1.5.6.4 | false_easting | 500 000 | 东移假定值 |
| 1.5.6.5 | false_northing | 0.0 | 北移假定值 |
| 1.6 | grid_id | 3BN44504 | 格网代码由5类元素组成，即：投影代码+南北坐标代码＋纵坐标格网代码＋东西坐标代码（或投影带号代码）＋横坐标格网代码。 |
| 1.7 | cloud_cover | 6.4918 | 非填充像素的云量百分比 |
| 1.8 | cloud_shadow | 5.9551 | 非填充像素的阴影百分比 |
| 1.9 | snow_ice | 0.0148 | 非填充像素的冰雪百分比 |
| 1.10 | fill | 64.9755 | 格网中填充像素百分比 |
| 1.11 | bands | N/A | 格网的波段信息 |
| 1.11.1 | product | rtu_qa | 波段类型 |
| 1.11.2 | source | rtu scene product | 来源 |
| 1.11.3 | name | PIXEL_QA | 波段名称 |
| 1.11.4 | data_type | UINT16 | 数据类型 |
| 1.11.5 | fill_value | 1 | 填充值 |
| 1.11.6 | nsamps | 3400 | 列数 |
| 1.11.7 | nlines | 3400 | 行数 |
| 1.11.8 | short_name | QA | 简称 |
| 1.11.9 | long_name | Land surface reflectence pixel quality band | 全称 |
| 1.11.10 | file_name | 3BN44504_L5_TM_20151126_V01_LSR_QA.tif | 文件名 |
| 1.11.11 | pixel_size units , x, y | meters, 30, 30 | 像元分辨率 |
| 1.11.12 | resample_method | none | 重采样方法 |
| 1.11.13 | data_units | quality/feature classification | 数据单元描述 |
| 1.11.14 | valid_range max, min | 65535.0, 0.0 | 最大、小值 |
| 1.11.15 | bit num | 1, 2, etc. | 位数 |
| 1.11.16 | production_date | 2018-03-20T20:35:13Z | 波段生成时间 |
致 谢
[编辑]感谢中科院A类先导专项“地球大数据科学工程”（XDA190090300）、国家自然基金重点项目（61731022）、国家重点研发计划课题（2016YFA0600302）对本研究的支持。
参考文献
[编辑]
