Suite

Rogner le raster du problème vectoriel


J'ai donc un fichier de formes d'un vecteur pour les lignes de la ville de Philadelphie et une image raster du comté de Philadelphie que j'ai reçue des charmantes personnes de l'USDA:NRCS. J'ai converti le raster d'un fichier sid en TIFF (à l'aide de gdal_translate), puis j'ai reprojeté le raster TIFF nouvellement converti en 3857 (à l'aide de gdalwarp).

Ce nouveau raster a désormais la même projection native que le vecteur des limites de la ville : EPSG:3857. Je charge donc à la fois le fichier de formes pour les limites de la ville et ce nouveau tif raster que j'ai préparé dans QGIS et je lance cet outil de découpage utile (Raster-> Extraction-> Clipper) et je m'occupe de définir le raster comme source et en sélectionnant le vecteur comme calque de masque et appuyez sur "OK". #c'était simple. L'ordinateur commence à faire son travail et je suis excité par l'arrivée de mon nouveau raster.

Ding ! ordinateur terminé… seuls South Philly et une infime partie de Central Philly sont extraits - West, North et East Philly semblent s'être retirés de la fête et tout est noir.

Hm… c'est bizarre, peut-être que je peux utiliser la ligne de commande GDAL pour forcer cela, parce que vous savez, si c'est fait sur la ligne de commande, je suis un vrai hacker et ce n'est pas comme si QGIS le faisait déjà en arrière-plan, mais je m'égare.

J'ai donc commandé ces mofos (en utilisant à nouveau gdalwarp !) 0… 10… 100 ! DING ! D'accord, d'accord, d'accord, alors que je me dirige moi-même vers le répertoire où j'ai sorti le TIFF tronqué. Ouvrez cette ventouse et je suis avec exactement la même sortie de Southside et Central qu'il ne reste plus qu'à représenter.

J'ai joint deux images, la première montre un polygone de couleur rose que je veux découper du raster, et la seconde montre la sortie les deux fois.

Photo 1 :

Photo 2:

FICHIERS SOURCE

Shapefile utilisé (en zip)

TIFF utilisé (en zip)

METTRE À JOUR: J'ai remarqué que les zones tronquées sont en fait des régions en dehors de la ville et des coupes juste au nord de la ville. (Notez que le sud de philly n'est plus traversé par la rivière Schuylkill). J'ai essayé d'en faire un polygone pour pouvoir télécharger une image pour plus de clarté, mais je ne pouvais pas le comprendre. Veuillez simplement reproduire le découpage de trame avec les fichiers fournis ci-dessus et vous pouvez également voir de quoi je parle.

Cela m'amène à croire qu'il doit y avoir un problème de projection, cependant, après avoir gdalinfo-ed le raster et consulté le .prj du fichier city_limits, j'ai confirmé qu'ils ont la même projection susmentionnée. #sayhwatt


Il semble y avoir une différence petite mais significative entre les projections qui fait trébucher gdalwarp.

En coursgdalsrsinfo philly_3857.tif > city_limits.prjet en retirer les premières lignes a fait l'affaire pour moi.

En comparaison:

philly_3857.tif :

PROJCS["WGS 84 / Pseudo-Mercator", GEOGCS["WGS 84", DATUM["WGS_1984", SPHEROID["WGS 84",6378137,298.257223563, AUTORITE["EPSG","7030"]], AUTORITE[" EPSG","6326"]], PRIMEM["Greenwich",0], UNIT["degré",0.0174532925199433], AUTORITE["EPSG","4326"]], PROJECTION["Mercator_1SP"], PARAMETRE["central_meridian ",0], PARAMETER["scale_factor",1], PARAMETER["false_easting",0], PARAMETER["false_northing",0], UNIT["metre",1, AUTHORITY["EPSG","9001"] ], EXTENSION["PROJ4","+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m [email protected] +wktext +no_defs"], AUTORITÉ["EPSG","3857"]]

city_limits.prj (l'original) :

PROJCS["WGS_84_Pseudo_Mercator", GEOGCS["GCS_WGS_1984", DATUM["D_WGS_1984", SPHEROID["WGS_1984",6378137,298.257223563]], PRIMEM["Greenwich",0], UNIT["Degree",0.017453292519943295 ["Mercator"], PARAMETRE["central_meridian",0], PARAMETER["false_easting",0], PARAMETER["false_northing",0], UNIT["Meter",1], PARAMETER["standard_parallel_1",0.0] ]

Pour un utilisateur normal, ils semblent être les mêmes. Mais si vous supprimez les fichiers .q** du zip, le fichier de formes est chargé dans QGIS en tant queEPSG : 54004, et ce sera ne pas aligner avec le raster. gdalwarp ne regarde pas dans les fichiers .q** et produit donc des déchets.


Voir la vidéo: Map Algebra Functionality and Advanced Raster Calculator Calculations (Octobre 2021).