Le portail francophone de la géomatique

Bonjour à tous,

dans le cadre de mon travail, je dois reprojeter un raster 100*100m (Corine Land Cove)r de EPSG:3035 (LAEA/ETRS89) à EPSG:2154 Lambert 93.

Pour cela, j'ai utilisé R (package raster) et en sortie le raster a l'air d'être parfait mais en zoomant et en comparant le résultat à l'original on se rend compte qu'il y a quelques déformations...surtout que le CLC est constitué d'une multitude d'ilots. Voici mon script :

clc<-raster("D:/users/.../clc2012/g100_clc12_V18_5.tif")#CLC copernicus epsg:3035##import
crs(clc)<-"+proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 +ellps=GRS80 +units=m +no_defs"#assignation de la proj epsg 3035

fr<-readOGR(dsn="D:/users/.../ZE-SHP", layer="france_corse_laea")#shape pour decouper la France
clc1<-crop(clc, ext)#decoupage de l'emprise

clc_fr<-mask(clc1, fr) #découpage selon le polygone
sr<-"+proj=lcc +lat_1=49 +lat_2=44 +lat_0=46.5 +lon_0=3 +x_0=700000 +y_0=6600000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs"
clc93<-projectRaster(clc_fr, crs=sr, method = 'ngb')#reprojection en spécifiant le nouveau crs, ici sr (lambert 93)

j'ai choisi "nearest neighbour" car ce sont des valeurs discrètes

J'ai réessayé 2 ou 3 fois, et à chaque fois les déformations sont différentes (toutes légères mais pas identiques). Ce ne sont pas les autres opérations qui produisent ces erreurs car j'ai sorti clc_fr (raster original découpé ) et après vérification; aucun problème.
Idem pour gdal Warp (sous QGIS), la sortie est aussi légèrement déformée...

Y a t-il une solution pour obtenir un raster identique à l'original? ou bien les déformation sont simplement dues aux différents systèmes de projection et on ne peut pas y échapper ?

Merci d'avance