Mere

Masseprojektion af rasterfiler i en loop [Linux]

Masseprojektion af rasterfiler i en loop [Linux]


Jeg har en mappe med mange.ascrasterfiler, der hver ledsages af en tilsvarende.prjfil. Filens basenavn er den samme for begge filer. Fremskrivningen kan dog variere mellem dem alle (de er i øjeblikket i forskellige UTM-projicering).
Hvad jeg ønsker at opnå er at standardisere projiceringen af ​​alle de individuelle rasterfiler. Min nuværende tilgang er dennebash-manuskript

#! / bin / bash echo Starter konvertering til WGS84 for f i * .asc; gør ekko Konvertering af "$ f"; ekko Til "$ {f% .asc} _wgs84.tif"; gdalwarp "$ f" "$ {f% .asc} _wgs84.tif" -t_srs "+ proj = longlat + ellps = WGS84"; færdig ekko Konvertering afsluttet

Imidlertid mislykkes scriptet for hver fil, der returneres:

"Ugyldig anmodning til pj_param, dødelig"

Så jeg vil gerne se en masse reprojicerede .tif (eller ethvert andet format) filer, der alle har den samme projektion.

Mens jeg skriver dette, arbejder jeg også på enpythonscript, som dog ikke synes at være så ligetil sombash.


opkald:gdalinfo grid2000.ascDet ser ud til at fejle selv i slutningen?

Driver: AAIGrid / Arc / Info ASCII Grid Files: grid2000.asc grid2000.prj Størrelsen er 16, 16 Koordinatsystemet er: PROJCS ["WGS_1984_UTM_Zone_30N", GEOGCS ["GCS_WGS_1984", DATUM ["WGS_1984", SPHEROID ["WGS_84" 6378137,298.257223563]], PRIMEM ["Greenwich", 0], UNIT ["Grad", 0,017453292519943295]], PROJEKTION ["Transverse_Mercator"], PARAMETER ["latitude_of_origin", 0], PARAMETER ["central_meridian", , PARAMETER ["scale_factor", 0.9996], PARAMETER ["false_easting", 500000], PARAMETER ["false_northing", 0], UNIT ["Meter", 1]] Origin = (565146.201264380943030,680158.854243048001081) Pixelstørrelse = (30,0000000000 -30.000000000000000) ugyldig anmodning til pj_param, dødelig

call: gdalsrsinfo xyz.prj # (Bemærk at xyz.prj kan være en anden fil)

PROJ.4: '+ proj = utm + zone = 17 + syd + datum = WGS84 + enheder = m + ingen_defs' OGC WKT: PROJCS ["WGS_1984_UTM_Zone_17S", GEOGCS ["GCS_WGS_1984", DATUM ["WGS_1984", SPHEROID [" WGS_84 ", 6378137,298.257223563]], PRIMEM [" Greenwich ", 0], UNIT [" Grad ", 0,017453292519943295]], PROJEKTION [" Transverse_Mercator "], PARAMETER [" latitude_of_origin ", 0], PARAMETER ["_ -81], PARAMETER ["scale_factor", 0.9996], PARAMETER ["false_easting", 500000], PARAMETER ["false_northing", 10000000], UNIT ["Meter", 1]]

Inden du kører dit script i mit Debian-system:

Det kørte perfekt uden nogen fejl (og uden ændringer hvor som helst):

Hvert billede blev tilstrækkeligt projiceret til WGS84.

Hvad er resultatet af en enkelt projektion med kommandoen ved bash-konsollen (uden script)? Er det med succes?


Se videoen: Grafički, zvučni i video formati datoteka