Himawari8
针对葵花8(Himawari8)数据的显示和专题产品合成算法集。
# AOT
葵花8数据的气溶胶光学厚度产品算法,返回气溶胶Image。
函数 | 返回值 |
---|---|
pie.Algorithm.Himawari8.AOT(image) | Image |
参数 | 类型 | 说明 |
---|---|---|
image | Image | 要进行计算的葵花8影像数据 |
# 示例
import pie
#选择葵花8数据
Map = pie.Map()
filter = pie.Filter.calendarRange(11, 12,"hour")
imageCollection = pie.ImageCollection("H08/AHI-L1-G").filter(filter)
image = imageCollection.first()
#计算AOT
AOT = pie.Himawari8.AOT(image).select("B1")
#默认分级
visParam= {
'min': 0.05,
'max': 1.5,
'opacity':1,
'palette': '00ff7f,ffff00,de9037,dd104e,b0287c,9c3741'
}
Map.setCenter(108.92, 34.54, 4)
Map.addLayer(AOT,visParam, "H8_AOT")
Map
# cloudMask
葵花8数据的云掩膜算法,返回云掩膜Image,云值为1,其余为0。
函数 | 返回值 |
---|---|
cloudMask(image) | Image |
参数 | 类型 | 说明 |
---|---|---|
image | Image | 要进行计算的葵花8影像数据 |
# 示例
#选择葵花8数据
Map = pie.Map()
image_H8 = pie.ImageCollection("H08/AHI-L1-G").first()
#计算云掩膜
cloudMask = pie.Himawari8.cloudMask(image_H8)
#选择前3个波段做掩膜运算
image_rm_cloud = image_H8.select(["B1","B2","B3"]).updateMask(cloudMask)
visParam = {
'min': 0,
'max': 3000,
'bands': ["B3","B2","B1"]
}
#添加到地图
Map.setCenter(109, 33, 3)
Map.addLayer(image_rm_cloud.select(["B1","B2","B3"]),visParam, "H8_cloludMask")
Map
# legend
葵花8数据的RGB合成图例,目前支持气团模式、夜间云微物理模式、白天云微物理模式、白天对流风暴模式等。
函数 | 返回值 |
---|---|
legend(type) | Legend |
参数 | 类型 | 说明 |
---|---|---|
type | String | 合成类型字符标识,主要包括: Airmass:气团模式 Night_Microphysics:夜间云微物理模式 Day_Microphysics:白天云微物理模式 Day_Convective_Storms:白天对流风暴模式 |
# 示例
Map=pie.Map()
legend1 = pie.Himawari8.legend(type="Airmass")
Map.addLegend(legend1)
Map
# PM25
葵花8数据的PM2.5产品算法,返回PM2.5 Image。
函数 | 返回值 |
---|---|
PM25(image) | Image |
参数 | 类型 | 说明 |
---|---|---|
image | Image | 要进行计算的葵花8影像数据 |
# 示例
# H8_PM25
import pie
Map = pie.Map()
imageCollection = pie.ImageCollection("H08/AHI-L1-G")
image = imageCollection.first()
#计算pm2.5
PM25 = pie.Himawari8.PM25(image).select("B1")
#分级
visParam = {
'opacity':1,
'classify':'35,75,115,150,250,500',
'palette': '00ff7f,ffff00,de9037,dd104e,b0287c,9c3741'
}
Map.setCenter(108.92, 34.54, 3)
Map.addLayer(PM25,visParam, "H8_PM25")
Map
# PM100
葵花8数据的PM10产品算法,返回PM10 Image。
函数 | 返回值 |
---|---|
pie.Algorithm.Himawari8.PM100(image) | Image |
参数 | 类型 | 说明 |
---|---|---|
image | Image | 要进行计算的葵花8影像数据 |
# 示例
#选择葵花8数据
Map = pie.Map()
filter = pie.Filter.calendarRange(11, 12,"hour")
imageCollection = pie.ImageCollection("H08/AHI-L1-G").filter(filter)
image = imageCollection.first()
#计算pm10
PM100 = pie.Himawari8.PM100(image).select("B1")
#默认分级
visParam = {
'min': 32,
'max': 38,
'opacity':1,
'palette': '00ff7f,ffff00,de9037,dd104e,b0287c,9c3741'
}
Map.setCenter(108.92, 34.54, 3)
Map.addLayer(PM100,visParam, "H8_PM10")
Map
# rgb
葵花8数据的多模式的RGB合成算法,返回合成后的影像Image,目前支持气团模式、夜间云微物理模式、白天云微物理模式、白天对流风暴模式、沙尘模式、对流云模式、强风暴模式、山火灰模式、积雪雾模式、自然色模式、真彩色模式等。
函数 | 返回值 |
---|---|
rgb(image,type) | Image |
参数 | 类型 | 说明 |
---|---|---|
image | Image | 要进行计算的葵花8影像数据 |
type | String | 合成类型字符标识,主要包括: Airmass:气团模式 Night_Microphysics:夜间云微物理模式 Day_Microphysics:白天云微物理模式 Day_Convective_Storms:白天对流风暴模式 Dust:沙尘模式 Clouds_Convection:对流云模式 Severe_Storm:强风暴模式 Ash:山火灰模式 Snow_Fog:积雪、雾模式 Natural:自然色模式 TrueColor:真彩色模式 |
# 示例
#影像RGB模式渲染,并加载到地图
Map = pie.Map()
imageCollection = pie.ImageCollection("H08/AHI-L1-G")
filter = pie.Filter.calendarRange(11, 12,"hour")
images_result = imageCollection.filter(filter)
imageRGB = pie.Himawari8.rgb(images_result.getAt(0),"Airmass")
Map.addLayer(imageRGB, {'min':0,'max':255}, "LayerRGB")
#设置地图显示范围
Map.setCenter(107.7002, 31.3354, 2)
Map
# sst
葵花8数据的海表温度产品提取算法,返回海温影像Image。
函数 | 返回值 |
---|---|
sst(image,type) | Image |
参数 | 类型 | 说明 |
---|---|---|
image | Image | 要进行计算的葵花8影像数据 |
# 示例
# 影像海表温度产品生成,并加载到地图
Map = pie.Map()
image = pie.ImageCollection("H08/AHI-L1-G").first()
imageSST = pie.Himawari8.sst(image)
visParam = {
'opacity':1.0,
'classify':'0,' +
'0.5,1,1.5,2,2.5,3,3.5,4,4.5,5,5.5,6,6.5,7,7.5,8,8.5,9,9.5,10,' +
'10.5,11,11.5,12,12.5,13,13.5,14,14.5,15,15.5,16,16.5,17,17.5,18,18.5,19,19.5,20,' +
'20.5,21,21.5,22,22.5,23,23.5,24,24.5,25,25.5,26,26.5,27,27.5,28,28.5,29,29.5,30,' +
'30.5,31',
'palette': '00008B,'+'000097,0000AF,0000BB,0000C7,0000DF,0000EC,0000F8,0010FF,001CFF,0028FF,0040FF,004CFF,0058FF,0070FF,007CFF,0087FF,0093FF,00ABFF,00B7FF,00C3FF,' +
'00DBFF,00E7FF,00F3FF,0cFFF8,18FFEC,24FFDF,3CFFC7,48FFBB,54FFAF,6CFF97,78FF8B,83FF80,8FFF74,A7FF5C,B3FF50,BFFF44,D7FF2C,E3FF1F,EFFF14,FFFB00,' +
'FFEF00,FFE300,FFCB00,FFBF00,FFB300,FF9B00,FF8F00,FF8300,FF7800,FF6000,FF5400,FF4800,FF3000,FF2400,FF1800,FF0000,F30000,E70000,CF0000,C30000,' +
'B70000,9F0000'
}
Map.addLayer(imageSST,visParam,"SST")
# 设置地图显示范围
Map.setCenter(112.6, 32.5, 2)
Map
有问题,我来改改 (opens new window)
上次更新: 2022/05/25, 07:37:39