Suite

Différences de projection dans R

Différences de projection dans R


J'essaie de reprojeter certaines données climatiques que j'ai téléchargées à partir de http://gdo-dcp.ucllnl.org/downscaled_cmip_projections/dcpInterface.html#Welcomein

Les données BCSD-CMIP5 se trouvent sous l'onglet « Projections : Archives complètes ».

Les fichiers sont au format netCDF, mais je rencontre quelques problèmes. Les données raster d'origine (appelons-les rast1) lorsqu'elles sont importées ont l'étendue

xmin 235,25 xmax 293 ymin 25,125 ymax 52,875

et un crs de +proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0

Les autres objets (à la fois sp et raster) avec lesquels je travaille contiennent des valeurs négatives pour xmin et xmax

xmin : -122.8021 xmax : -75.38249 ymin : 31.78881 ymax : 47.15952

Lorsque j'essaie de reprojeter rast1 à l'aide des autres objets (rast2), l'étendue ne change pas et même si le rast1 reprojeté est dans le même crs, les étendues ne s'alignent pas.

rast_new = projetRaster(rast1, rast2)

J'ai donc dû faire les choses suivantes :

rast_new = projectExtent(rast1,rast2) rast_new = projectRaster(rast1,rast_new) values(rast_new) = values(rast1) ### Ceci est nécessaire car toutes les valeurs étaient NA

Cependant, même si rast_new est dans la même projection que rast2

+proj=aea +lat_1=20 +lat_2=60 +lat_0=40 +lon_0=-96 +x_0=0 +y_0=0 +datum=NAD83 +units=m +no_defs +ellps=GRS80+towgs84=0,0, 0

le tracé n'est pas le même. Voir l'image

Je me demande simplement ce que je peux faire pour corriger cela ? Inutile de dire que je suis relativement nouveau dans la compréhension des projections dans R (et des projections en général !)

J'ai exporté les deux rasters et visualisés dans ArcGIS, pour voir s'il s'agissait d'un problème de plot() dans R, mais les images sont les mêmes que ci-dessus, mais les projections sont les mêmes (selon ArcGIS 10.2). Même lorsque je reprojette le "mauvais" raster dans ArcGIS, la projection ne semble pas être dans Albers Equal Area. Lorsque j'importe directement depuis le fichier netCDF dans ArcGIS, la projection semble correcte.


Il semble que la première image utilise le système de référence géodésique WSG84, tandis que la deuxième image utilise GRS80, c'est pourquoi elle semble incurvée. Cela suggère que les objets n'ont pas été convertis dans le même SCR.

Pouvez-vous fournir les données climatiques, ou un échantillon, afin que moi-même ou un autre utilisateur puisse approfondir la question ?


J'ai posté sur R-SIG-GEO et j'ai obtenu cette réponse du Dr Hijmans.

http://r-sig-geo.2731867.n2.nabble.com/Raster-projection-alignment-issues-td7587564.html

Fondamentalement, vous devez ajuster manuellement l'étendue de l'objet raster. Spécifiquement,

étendue(r) <- étendue(235.25-360, 293-360, 25.125, 52.875),

ou utilisez la fonction de rotation du package raster


Voir la vidéo: First angle projection vs Third angle projection in Tamil. Planes of projection in Tamil