30_geemap学习笔记 | 影像属性统计

原文：30 image pops stats

import ee
import geemap
geemap.set_proxy(port=10809)

Map = geemap.Map()
Map

centroid = ee.Geometry.Point(117.1439, 30.7538)

landsat = ee.ImageCollection('LANDSAT/LC08/C01/T1_SR') \
                     .filterBounds(centroid)\
                     .first()
landsat_vis = {
    'min': 0,
    'max': 3000,
    'bands': ['B5', 'B4', 'B3']
}
Map.centerObject(centroid, zoom=8)
Map.addLayer(landsat, landsat_vis, 'Landsat-8')

qsshp = "G:/learnpy/shp/潜山市.shp"
qs_ee = geemap.shp_to_ee(qsshp)

centroid = ee.Geometry.Point(117.1439, 30.7538)

landsat = ee.ImageCollection('LANDSAT/LC08/C01/T1_SR') \
                     .filterBounds(centroid) \
                     .first()
landsat = landsat.clip(qs_ee)

landsat_vis = {
    'min': 0,
    'max': 3000,
    'bands': ['B5', 'B4', 'B3']
}
Map.centerObject(centroid, zoom=8)
Map.addLayer(landsat, landsat_vis, 'Landsat-8')

worldcover = ee.ImageCollection('ESA/WorldCover/v100').first()
worldcover = worldcover.clip(qs_ee)
worldcover_vis = {
    'bands': ['Map']
}

Map.addLayer(worldcover, worldcover_vis, 'WorldCover')

landsat.propertyNames().getInfo()

# ['IMAGE_QUALITY_TIRS',
#  'CLOUD_COVER',
#  'system:id',
#  'EARTH_SUN_DISTANCE',
#  'LANDSAT_ID',
#  'system:footprint',
#  'system:version',
#  'CLOUD_COVER_LAND',
#  'GEOMETRIC_RMSE_MODEL',
#  'SR_APP_VERSION',
#  'SATELLITE',
#  'SOLAR_AZIMUTH_ANGLE',
#  'IMAGE_QUALITY_OLI',
#  'WRS_PATH',
#  'system:time_start',
#  'SENSING_TIME',
#  'ESPA_VERSION',
#  'SOLAR_ZENITH_ANGLE',
#  'WRS_ROW',
#  'GEOMETRIC_RMSE_MODEL_Y',
#  'LEVEL1_PRODUCTION_DATE',
#  'GEOMETRIC_RMSE_MODEL_X',
#  'system:asset_size',
#  'PIXEL_QA_VERSION',
#  'system:index',
#  'system:bands',
#  'system:band_names']

landsat.get('CLOUD_COVER').getInfo()

# 32.32

landsat.get('system:time_start').getInfo()

# 1364438877120

ee.Date(landsat.get('system:time_start')).format('YYYY-MM-dd').getInfo()

# '2013-03-28'

一次统计影像所有属性

landsat_props = geemap.image_props(landsat)
landsat_props.getInfo()

# {'CLOUD_COVER': 32.32,
#  'CLOUD_COVER_LAND': 32.32,
#  'EARTH_SUN_DISTANCE': 0.998104,
#  'ESPA_VERSION': '2_23_0_1a',
#  'GEOMETRIC_RMSE_MODEL': 7.379,
#  'GEOMETRIC_RMSE_MODEL_X': 4.944,
#  'GEOMETRIC_RMSE_MODEL_Y': 5.477,
#  'IMAGE_DATE': '2013-03-28',
#  'IMAGE_QUALITY_OLI': 9,
#  'IMAGE_QUALITY_TIRS': 9,
#  'LANDSAT_ID': 'LC08_L1TP_121039_20130328_20170505_01_T1',
#  'LEVEL1_PRODUCTION_DATE': 1493996942000,
#  'NOMINAL_SCALE': 30,
#  'PIXEL_QA_VERSION': 'generate_pixel_qa_1.6.0',
#  'SATELLITE': 'LANDSAT_8',
#  'SENSING_TIME': '2013-03-28T02:47:57.1239510Z',
#  'SOLAR_AZIMUTH_ANGLE': 136.924896,
#  'SOLAR_ZENITH_ANGLE': 34.939472,
#  'SR_APP_VERSION': 'LaSRC_1.3.0',
#  'WRS_PATH': 121,
#  'WRS_ROW': 39,
#  'system:asset_size': '547.006567 MB',
#  'system:band_names': ['B1',
#   'B2',
#   'B3',
#   'B4',
#   'B5',
#   'B6',
#   'B7',
#   'B10',
#   'B11',
#   'sr_aerosol',
#   'pixel_qa',
#   'radsat_qa'],
#  'system:id': 'LANDSAT/LC08/C01/T1_SR/LC08_121039_20130328',
#  'system:index': 'LC08_121039_20130328',
#  'system:time_end': '2013-03-28 02:47:57',
#  'system:time_start': '2013-03-28 02:47:57',
#  'system:version': 1522736810313046}

landsat_props.get('IMAGE_DATE').getInfo()

# '2013-03-28'

landsat_stats = geemap.image_stats(landsat, region=qs_ee)
landsat_stats.getInfo()

# {'max': {'B1': 2179,
#   'B10': 2987,
#   'B11': 2953,
#   'B2': 2336,
#   'B3': 2966,
#   'B4': 3378,
#   'B5': 4424,
#   'B6': 14613,
#   'B7': 15096,
#   'pixel_qa': 992,
#   'radsat_qa': 0,
#   'sr_aerosol': 228},
#  'mean': {'B1': 384.7877815853934,
#   'B10': 2858.190347764515,
#   'B11': 2833.479935497505,
#   'B2': 491.83200772308857,
#   'B3': 774.8468440178282,
#   'B4': 756.4534018825602,
#   'B5': 2078.2957270111347,
#   'B6': 1505.1009926042339,
#   'B7': 992.8861386863474,
#   'pixel_qa': 354.75008996351823,
#   'radsat_qa': 0,
#   'sr_aerosol': 180.83915643256015},
#  'min': {'B1': -227,
#   'B10': 2776,
#   'B11': 2751,
#   'B2': -26,
#   'B3': 209,
#   'B4': 180,
#   'B5': 436,
#   'B6': 307,
#   'B7': 205,
#   'pixel_qa': 322,
#   'radsat_qa': 0,
#   'sr_aerosol': 8},
#  'std': {'B1': 178.94111660085977,
#   'B10': 23.56345942116287,
#   'B11': 23.343755876476433,
#   'B2': 166.83367465733448,
#   'B3': 186.02055491332948,
#   'B4': 220.85221724840423,
#   'B5': 384.64481093713005,
#   'B6': 360.5429524594576,
#   'B7': 306.7369108807791,
#   'pixel_qa': 126.80182301907443,
#   'radsat_qa': 0,
#   'sr_aerosol': 58.2271929462684},
#  'sum': {'B1': 720889072.8745103,
#   'B10': 5354739127.654904,
#   'B11': 5308444866.137251,
#   'B2': 921433416.0941176,
#   'B3': 1451653741.9725494,
#   'B4': 1417194146.1058817,
#   'B5': 3893628518.117649,
#   'B6': 2819764324.819607,
#   'B7': 1860144220.3764708,
#   'pixel_qa': 664614303.5058821,
#   'radsat_qa': 0,
#   'sr_aerosol': 338797067.0039214}}

worldcover_stats = geemap.image_stats(worldcover, region=qs_ee)
worldcover_stats.getInfo()

# {'max': {'Map': 90},
#  'mean': {'Map': 22.27667768841708},
#  'min': {'Map': 10},
#  'std': {'Map': 17.44834559670808},
#  'sum': {'Map': 509914715.41176456}}

worldcover_props = geemap.image_props(worldcover)
worldcover_props.getInfo()

# {'IMAGE_DATE': '2020-01-01',
#  'Map_class_names': ['Trees',
#   'Shrubland',
#   'Grassland',
#   'Cropland',
#   'Built-up',
#   'Barren / sparse vegetation',
#   'Snow and ice',
#   'Open water',
#   'Herbaceous wetland',
#   'Mangroves',
#   'Moss and lichen'],
#  'Map_class_palette': ['006400',
#   'ffbb22',
#   'ffff4c',
#   'f096ff',
#   'fa0000',
#   'b4b4b4',
#   'f0f0f0',
#   '0064c8',
#   '0096a0',
#   '00cf75',
#   'fae6a0'],
#  'Map_class_values': [10, 20, 30, 40, 50, 60, 70, 80, 90, 95, 100],
#  'NODATA_value': '0',
#  'NOMINAL_SCALE': 9.276624232772797,
#  'copyright': 'ESA WorldCover project 2020 / Contains modified Copernicus Sentinel data (2020) processed by ESA WorldCover consortium',
#  'creator': 'ESA WorldCover project',
#  'doi': 'NA',
#  'institution': 'VITO NV',
#  'orbit_type': '',
#  'platform': 'Sentinel-2 and Sentinel-1',
#  'processing_level': 'L3',
#  'processing_mode': 'offline',
#  'product_version': 'V1.0.0',
#  'production_crs': 'EPSG:4326',
#  'production_grid': '3x3 degree tiling grid',
#  'references': 'https://esa-worldcover.org',
#  'sensor': '',
#  'source': 'Derived from EO satellite imagery',
#  'system:asset_size': '109355.097577 MB',
#  'system:band_names': ['Map'],
#  'system:id': 'ESA/WorldCover/v100/2020',
#  'system:index': '2020',
#  'system:time_end': '2021-01-01 00:00:00',
#  'system:time_start': '2020-01-01 00:00:00',
#  'system:version': 1634824277829985}