第七色在线视频,2021少妇久久久久久久久久,亚洲欧洲精品成人久久av18,亚洲国产精品特色大片观看完整版,孙宇晨将参加特朗普的晚宴

為了賬號(hào)安全,請(qǐng)及時(shí)綁定郵箱和手機(jī)立即綁定
已解決430363個(gè)問(wèn)題,去搜搜看,總會(huì)有你想問(wèn)的

在 Google Earth Engine 中的 imageCollection 上對(duì)每個(gè)圖像執(zhí)行

在 Google Earth Engine 中的 imageCollection 上對(duì)每個(gè)圖像執(zhí)行

qq_花開(kāi)花謝_0 2022-10-08 17:06:01
我需要對(duì)圖像集合執(zhí)行每個(gè)圖像的 PCA。然后,我只想保留主組件軸 1,并將其作為一個(gè)帶添加到我的圖像集中的每個(gè)圖像中。最終,我想導(dǎo)出一個(gè) .csv 文件,其中行標(biāo)題處的 GPS 采樣位置和圖像 ID 作為列標(biāo)題,平均主成分軸 1 作為值。這樣做背后的想法是,我想要一個(gè)代理(光譜異質(zhì)性)用于 R 中的進(jìn)一步統(tǒng)計(jì)分析。這是我到目前為止的代碼://Create an test image to extract information to be used during PCAvar testImage =ee.Image('LANDSAT/LC08/C01/T1_SR/LC08_168080_20130407').select(['B2', 'B3', 'B4', 'B5', 'B6', 'B7'],        ['Blue', 'Green', 'Red', 'NIR', 'SWIR1', 'SWIR2']);// Define variables for PCAvar region = Extent;var scale = testImage.projection().nominalScale();var bandNames = testImage.bandNames();Map.centerObject(region);// Function for performing PCAfunction doPCA(image){  // This code is from https://code.earthengine.google.com/7249153a8a0f5c79eaf562ed45a7adadvar meanDict = image.reduceRegion({    reducer: ee.Reducer.mean(),    geometry: region,    scale: scale,    maxPixels: 1e9});var means = ee.Image.constant(meanDict.values(bandNames));var centered = image.subtract(means);// This helper function returns a list of new band names.var getNewBandNames = function(prefix) {  var seq = ee.List.sequence(1, bandNames.length());  return seq.map(function(b) {    return ee.String(prefix).cat(ee.Number(b).int());  });};// [START principal_components]var getPrincipalComponents = function(centered, scale, region) {  var arrays = centered.toArray();  var covar = arrays.reduceRegion({    reducer: ee.Reducer.centeredCovariance(),    geometry: region,    scale: scale,    maxPixels: 1e9  });Extent是我的 ROI,而LandsatCol是經(jīng)過(guò)預(yù)處理的圖像集。嘗試將 PCA 映射到圖像集合(代碼的倒數(shù)第二行)時(shí),此處的代碼會(huì)產(chǎn)生錯(cuò)誤。錯(cuò)誤顯示:“數(shù)組:需要參數(shù)‘值’”。關(guān)于如何處理這個(gè)問(wèn)題的任何建議?以及如何在圖像集合上添加主分量軸 1 作為每個(gè)圖像的帶?
查看完整描述

1 回答

?
智慧大石

TA貢獻(xiàn)1946條經(jīng)驗(yàn) 獲得超3個(gè)贊

我想到了。錯(cuò)誤“Array: Parameter 'values' is required”與稀疏矩陣有關(guān),稀疏矩陣是執(zhí)行 PCA 的過(guò)濾、裁剪和指定區(qū)域的產(chǎn)物。Earth Engine 無(wú)法處理稀疏矩陣。


這是工作代碼。LandsatCol是我的預(yù)處理圖像集。


// Display AOI

var point = ee.Geometry.Point([30.2261, -29.458])

Map.centerObject(point,10);


// Prepairing imagery for PCA

var Preped = LandsatCol.map(function(image){

  var orig = image;

  var region = image.geometry();

  var scale = 30;

  var bandNames = ['Blue', 'Green', 'Red', 'NIR', 'SWIR1', 'SWIR2'];

  var meanDict = image.reduceRegion({

    reducer: ee.Reducer.mean(),

    geometry: region,

    scale: scale,

    maxPixels: 1e9

  });

  var means = ee.Image.constant(meanDict.values(bandNames));

  var centered = image.subtract(means);

  var getNewBandNames = function(prefix) {

  var seq = ee.List.sequence(1, 6);

  return seq.map(function(b) {

    return ee.String(prefix).cat(ee.Number(b).int());

    });

  };


  // PCA function

  var getPrincipalComponents = function(centered, scale, region) {

    var arrays = centered.toArray();

    var covar = arrays.reduceRegion({

      reducer: ee.Reducer.centeredCovariance(),

      geometry: region,

      scale: scale,

      maxPixels: 1e9

    });

    var covarArray = ee.Array(covar.get('array'));

    var eigens = covarArray.eigen();

    var eigenValues = eigens.slice(1, 0, 1);

    var eigenVectors = eigens.slice(1, 1);

    var arrayImage = arrays.toArray(1);

    var principalComponents = ee.Image(eigenVectors).matrixMultiply(arrayImage);

    var sdImage = ee.Image(eigenValues.sqrt())

    .arrayProject([0]).arrayFlatten([getNewBandNames('sd')]);

    return principalComponents.arrayProject([0])

    .arrayFlatten([getNewBandNames('pc')])

    .divide(sdImage);

    };


  var pcImage = getPrincipalComponents(centered, scale, region);

  return ee.Image(image.addBands(pcImage));

});

print("PCA imagery: ",Preped);


查看完整回答
反對(duì) 回復(fù) 2022-10-08
  • 1 回答
  • 0 關(guān)注
  • 172 瀏覽
慕課專欄
更多

添加回答

舉報(bào)

0/150
提交
取消
微信客服

購(gòu)課補(bǔ)貼
聯(lián)系客服咨詢優(yōu)惠詳情

幫助反饋 APP下載

慕課網(wǎng)APP
您的移動(dòng)學(xué)習(xí)伙伴

公眾號(hào)

掃描二維碼
關(guān)注慕課網(wǎng)微信公眾號(hào)