Conversion de fichiers ascii windspeed en rasters

11

J'ai essayé de convertir des fichiers de vitesse du vent en rasters sans succès. Il a été discuté dans un forum Arcgis que j'ai suivi. Cela impliquait de remplacer les points-virgules par des espaces et d'ajouter un en-tête.

J'espérais inclure l'un des fichiers .asc d'origine avec une version .csv que j'ai «nettoyée». Malheureusement, je ne peux pas voir comment faire cela sans copier et coller le texte. La première ligne du fichier .asc d'origine est:

Conversion of speed45i.dat   Wind m/s  Wind Speed @ 45m N. Ire 
(100,  499); 0.0; 0.0; 0.0; 0.0; 0.0; 0.0; 0.0; 0.0; 0.0; 0.0; 0.0; 0.0; 0.0; 0.0; 0.0; 0.0; 0.0; 0.0; 0.0; 0.0; 0.0; 0.0; 0.0; 

J'ai changé cela en:

ncols 700
nrows 1300
x||corner 0
y||corner 0
cellsize 1000                                                                                               "
"( 0 1299) 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 ;"

Lorsque j'utilise le plug-in de format de conversion de conversion raster, j'ai reçu le message d'erreur:

Le processus n'a pas pu démarrer. Soit le programme appelé est manquant, soit vous ne disposez pas des autorisations suffisantes pour appeler le programme.

Toutes les suggestions reçues avec reconnaissance.

tedwalsh
la source
vous pouvez toujours essayer de mettre un exemple dans DropBox et partager un lien, si vous voulez une aide plus spécifique avec l'ensemble de données particulier.
RyanKDalton

Réponses:

14

Je connais très bien ces données. Ce sont les données de vitesse du vent NOABL légèrement notoires. Vous êtes également sur la bonne voie pour les convertir en raster ASCII. L'en-tête que j'ai créé pour eux (il y a de nombreuses années) était le suivant:

ncols 700
nrows 1300
xllcorner 0
yllcorner 0
cellsize 1000
nodata_value -999
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 

Il semble que vous ayez utilisé des symboles de tuyau ('|') au lieu de 'l' pour xll et yll dans l'en-tête (qui signifie 'en bas à gauche'). Perdez également les chiffres entre parenthèses au début de chaque ligne ainsi que les virgules inversées. Rendez-le simplement délimité par des espaces. Les chiffres entre parenthèses (comme je suis sûr que vous avez calculé) indiquent les coordonnées de départ de 100 kilomètres de la ligne de données (qui est à une résolution de 1 km - il y a donc une centaine de valeurs après chaque ensemble de parenthèses). Ces nombres sont non seulement redondants dans un raster ESRI ASCII, mais ils le «cassent» également.

CAVEAT:
utilisez ces données avec prudence. Ils sont très anciens et ont été interpolés à partir d'une toute petite poignée de stations météorologiques réparties à travers le Royaume-Uni. Ils ont également été interpolés sans tenir compte de la rugosité de la surface et de la plupart des autres facteurs qui influencent la vitesse du vent au niveau local. Ils sont OK à titre indicatif mais sont connus pour surestimer la vitesse du vent dans de nombreux endroits. Par exemple, je trouve régulièrement que la vitesse du vent NOABL 45m est souvent comparable aux vitesses du vent mesurées par un anémomètre à des hauteurs de 80m, mais là encore, puisque NOABL n'a utilisé aucune donnée des îles du Nord, j'ai trouvé NOABL sous-estimer tristement les vitesses du vent dans Shetland. Le plus au nord de l'une des rares stations rencontrées (seulement 10 je pense!) Dans les données NOABL était Wick.

Tous les avertissements mis à part cependant, j'ai trouvé que lorsqu'ils sont utilisés à bon escient, en comprenant la surestimation (et avec une pincée d'expérience), j'ai trouvé que les résultats des calculs des rendements des parcs éoliens utilisant ces données étaient comparables aux calculs utilisant des données plus à jour et des données rigoureusement collectées. Cependant, ne prenez pas cela comme un évangile car chaque site est différent et la vitesse du vent réelle est grandement affectée par la rugosité, la portée ouverte, la topographie pour ne nommer que quelques facteurs.

Si la vitesse du vent est essentielle à votre projet, utilisez-la certainement comme point de départ, mais ne pariez pas sur la ferme!

EDIT
Pour supprimer les nombres entre parenthèses, vous pouvez écrire un script Python pour analyser les données ou une alternative simple consiste à:

  • ouvrir le fichier dans un éditeur de texte, supprimer la première ligne (" Conversion de vitesse ... "), remplacer globalement tous les espaces par rien et ré-enregistrer avec l'extension '.csv'
  • Ouvrez ensuite dans un package de feuille de calcul en définissant le délimiteur comme un point-virgule (dans Excel et OpenOffice Calc, vous serez invité à ouvrir le fichier). Sélectionnez maintenant la première colonne qui contient tous les numéros entre crochets et supprimez-la.
  • Insérez six lignes au tout début et ajoutez votre en-tête ("ncols 700", etc.).
  • Enfin, enregistrez avec l'extension '.asc' mais assurez-vous que votre délimiteur est un espace de sauvegarde.
MappaGnosis
la source
Sylvester, Merci beaucoup. Problème lors de la suppression des numéros entre parenthèses. Les premières lignes doivent-elles apparaître comme ci-dessous? nrows 1300 xllcorner 0 yllcorner 0 cellsize 1000 nodata_value -999 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0 0,0
Tedwalsh
@tedwalsh: c'est à peu près tout. J'ai modifié ma réponse ci-dessus pour ajouter la plupart de la première ligne à l'exemple de «code». J'ai également ajouté des instructions pour retirer facilement les supports. Si vous rencontrez toujours des difficultés, envoyez-moi un MP avec votre adresse e-mail (trouvez la section Contacts de mon site Web) et je trierai quelque chose pour vous.
MappaGnosis
1
L'ensemble de données NOABL était assez soigné pour l'époque, mais il y a près de 20 ans. Je me souviens aussi avoir travaillé sur / avec. Il est préférable de ne pas l'utiliser pour de vrais travaux à des hauteurs d'éoliennes typiques maintenant.
scruss
11

La valeur NoData est manquante dans votre fichier ascii et vous avez x || et y || au lieu de xll et yll. Je suppose que votre NoData est -999.

essaye ça:

NCOLS 700
NROWS 1300
XLLCORNER 0
YLLCORNER 0
CELLSIZE 1000
NODATA_VALUE -999

Votre ascii ouvert dans un éditeur de texte devrait ressembler à ceci:

NCOLS 700
NROWS 1300
XLLCORNER 0
YLLCORNER 0
CELLSIZE 1000
NODATA_VALUE -999
0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 
Gago-Silva
la source
2
En termes de format, cela semble correct, mais je pense que ce peut être une erreur importante de confondre zéro avec une valeur NoData. Un zéro est nécessaire pour indiquer une vitesse du vent nulle, tandis qu'une valeur NoData indiquerait l'absence d'informations sur la vitesse du vent.
whuber
Je supposais que NoData était 0, car @tedwalsh n'a pas donné cette information.
Gago-Silva