Python – Fixing Google Earth Engine EEException for Type Mismatch

geemapgoogle-earth-enginepython

This code selects Sentinel-2 data and de-clouds it using the S2_CLOUD_PROBABILITY data set. Then, cloud-free Sentinel-2 images and NDVI were synthesized.The problem is that the images cannot be displayed and I get the next error message:

EEException: Filter.eq, argument 'leftField': Invalid type.
Expected type: String.
Actual type: Dictionary<String>.
Actual value: {rightField=system:index, leftField=system:index}

I do all this in Python Jupyter notebook and using the package 'geemap'. This is my code:

river_shp='D:/paper1 folder/data/shp/Changjiang/CJ.shp'

river = geemap.shp_to_ee(river_shp)

S2_Cloud = ee.ImageCollection("COPERNICUS/S2_CLOUD_PROBABILITY")

S2 =ee.ImageCollection("COPERNICUS/S2")

roi = river

styling = {'color':'black', 'fillColor':'00000000'}

Map=geemap.Map()

Map.centerObject(roi,6)


def CloudProbability(img,thread):
  prob = img.select('probability')
  return img.updateMask(prob.lte(thread))

def MergeImages(primary, secondary):
  join = ee.Join.inner()
  filter = ee.Filter.equals({
    'leftField': 'system:index',
    'rightField': 'system:index'
  })
  Col = join.apply(primary, secondary, filter)
  def func_ach(image):
    img1 = ee.Image(image.get("primary"))
    img2 = ee.Image(image.get("secondary"))
    return ee.Image.cat(img1,img2)
  Col = Col.map(func_ach)
  return ee.ImageCollection(Col)

def main():
  startDate = "2021-5-1"
  endDate = "2021-10-1"
  S2Img1 = S2.filterDate(startDate, endDate) \
                  .filterBounds(roi)
  S2Img2 = S2_Cloud.filterDate(startDate, endDate) \
                        .filterBounds(roi)
  S2Imgs = MergeImages(S2Img1, S2Img2)
  def func_irn(image):
    return CloudProbability(image, 30)
  S2Imgs = S2Imgs.map(func_irn)

  def func_tjh (img):
    return img.normalizedDifference(['B8','B4'])
  ndvi_max = S2Imgs.map(func_tjh).max().clip(roi)
  Map.addLayer(S2Imgs,{'bands':['B8', 'B4', 'B3'], 'min':0, 'max':3000},'S2Imgs')
  Map.addLayer(ndvi_max, {'min':0, 'max':0.9}, "ndvi")
main()
Map.addLayer(roi.style(styling),{},'CJ')

Best Answer

This is how to use ee.Filter.equals in python: filter = ee.Filter.equals(leftField= 'system:index',rightField= 'system:index')

Related Solutions

Google Earth Engine – Fix Layer Error: Element.get, Argument ‘Property’: Invalid Type

Convert the filter ImageCollection to a List of images.

When the get method is applied to an ImageCollection it expects the name (string) of a property to get; when it is applied to a List, it expects an integer index value and returns the respective image from the List.

filter = filter.toList(filter.size()); // convert ImageCollection to List

See line 97: https://code.earthengine.google.com/e4f2efb32d16a45e918a5ff60ab78b03

Note, however, that converting collections to lists can negatively affect performance. Avoid using .toList() if you can use alternative methods. For instance, use .filter() to select and image of interest.

Google Earth Engine – Resolving Image.subtract Error: Invalid Image Type

This error appears because you are mixing in indices_terra_2 function, images with ImageCollections in mndisi calculation. Following code fixes such error by using list of images, instead ImageCollection, for retrieving adequate index in LST images.

var Jawa = ee.FeatureCollection("users/mohfaisal181099/KP/shp_jawa");
Map.centerObject(Jawa,7);

Map.addLayer(Jawa);

//Define Research Time
  var start = ee.Date('2020-01-01');
  var end = ee.Date('2020-12-31');

//Time Variable
var ndays = ee.Number(end.difference(start,'day')).round();
var nmonths = ee.Number(end.difference(start,'month')).round();
var nyears = ee.Number(end.difference(start,'year')).round();
print('ndays', ndays);
print('nmonths', nmonths);
print('nyears', nyears);

//Convert Kelvin to Celcius
function celcius(img){
    return img.multiply(0.02).subtract(273.15)
              .copyProperties(img,['system:time_start']);
}
  
//LST Terra MODIS 
var terra_kelvin = ee.ImageCollection('MODIS/006/MOD11A1')
                       .filterBounds(Jawa)
                       .select('LST_Day_1km')
                       .filterDate(start,end);

var terra_LST = terra_kelvin.map(celcius);
print("terra_LST", terra_LST);

var terra_LST_lst = terra_LST.toList(terra_LST.size());

var modisProjection = terra_LST.first().projection();
print('MODIS Projection:', modisProjection.nominalScale());


//Downsetting Function
function reduceres(image){
  return image.reduceResolution({
    reducer : ee.Reducer.mean(),
    maxPixels : 1024
  })
  .reproject({
    crs : modisProjection
  });
}

//Define Dataset
//Terra MODIS Surface Reflectance
var terra_SR = ee.ImageCollection('MODIS/006/MOD09GA')
                     .filterBounds(Jawa)
                     .filterDate(start,end)
                     .map(reduceres);

print("terra_SR", terra_SR);

//Function Indices Terra MODIS
var indices_terra_1 = function(a){
  //MNDWI
  var mndwi = a.expression(
                  '(Green - SWIR)/(Green + SWIR)',
                   {'Green': a.select('sur_refl_b04'),
                    'SWIR': a.select('sur_refl_b02')
                   }).rename('MNDWI');
    return a.addBands(mndwi)
            .copyProperties(a,['system:time_start']);
    };
  
var terra_mndwi = terra_SR.map(indices_terra_1);
var tst = terra_mndwi.first().projection();

print("terra_mndwi", terra_mndwi);

print('MNDWI Projection:', tst.nominalScale());

//MNDISI
var terra_mndwi_lst = terra_mndwi.toList(terra_mndwi.size());

function indices_terra_2(aa){
  
  var idx = terra_mndwi_lst.indexOf(aa);
  
  var mndisi = ee.Image(aa).expression(
                  '(LST - ((MNDWI + NIR + SWIR)/3))/(LST + ((MNDWI + NIR + SWIR)/3))',
                   {'LST': ee.Image(ee.List(terra_LST_lst).get(idx)).select('LST_Day_1km'),
                    'MNDWI': ee.Image(aa).select('MNDWI'),
                    'NIR' : ee.Image(aa).select('sur_refl_b02'),
                    'SWIR': ee.Image(aa).select('sur_refl_b06')
                   }).rename('MNDISI');
    
    return ee.Image(aa).addBands(mndisi)
             .copyProperties(aa,['system:time_start']);

    }
    
var terra_mndisi = terra_mndwi_lst.map(indices_terra_2);

print("terra_mndisi", terra_mndisi);

var tst2 = ee.Image(terra_mndisi.get(0)).projection();
print('MNDISI Projection:', tst2.nominalScale());

After running above code in GEE code editor, it can be observed in following image that MNDWI and MNDISI bands were added to Image Collection list and nominalScale of first image was printed as expected.

Best Answer

Related Solutions

Google Earth Engine – Fix Layer Error: Element.get, Argument ‘Property’: Invalid Type

Google Earth Engine – Resolving Image.subtract Error: Invalid Image Type

Related Question