Suite

Comment trouver l'index d'emplacement de cellule d'un raster à l'aide d'informations Lat long?


J'ai un raster pour lequel je veux trouver l'emplacement de la cellule (Ligne, Col) si j'ai les informations Lat-long.

Je sais comment obtenir la valeur d'une cellule particulière en utilisant lat long.

extraire(ras,SpatialPoints(cbind(-82.8,35.2)))

Mais j'ai besoin de la colonne de ligne de cette cellule plutôt que de sa valeur. J'utilise le package raster dans R.

Raster : https://www.dropbox.com/s/8nhfirxr2hm3l4v/fdr_fb.tif?dl=0


# Ouvrir le fichier tif ras <- raster("fdr_fb.tif") # La fonction "extract" peut renvoyer le numéro de cellule, ici 8766 extract(ras,SpatialPoints(cbind(-82.8,35.2)), cellnumbers=TRUE) # cells fdr_fb # 8766 2 # Nous stockons le numéro de cellule id.cell <- extract(ras,SpatialPoints(cbind(-82.8,35.2)), cellnumbers=TRUE)[1] # Nous extrayons les indices de ligne et de col avec "rowColFromCell" rowColFromCell(ras , id.cell) # ligne col # [1,] 75 34 # Vérification ras[75, 34] # 2