Ini ialah arahan r.in.xyzgrass yang boleh dijalankan dalam penyedia pengehosan percuma OnWorks menggunakan salah satu daripada berbilang stesen kerja dalam talian percuma kami seperti Ubuntu Online, Fedora Online, emulator dalam talian Windows atau emulator dalam talian MAC OS
JADUAL:
NAMA
r.in.xyz - Mencipta peta raster daripada himpunan banyak koordinat menggunakan univariate
perangkaan.
Kata-kata berkaitan
raster, import, penukaran, pengagregatan, ASCII, LIDAR
SINOPSIS
r.in.xyz
r.in.xyz - membantu
r.in.xyz [-sgi] input=nama output=nama [kaedah=rentetan] [jenis=rentetan]
[pemisah=watak] [x=integer] [y=integer] [z=integer] [skip=integer]
[zrange=min, maks] [zscale=terapung] [lajur_nilai=integer] [vrange=min, maks]
[vscale=terapung] [peratus=integer] [pth=integer] [mengurangkan=terapung] [--menimpa]
[--membantu] [--kata kerja] [--tenang] [--ui]
Bendera:
-s
Imbas fail data untuk tahap kemudian keluar
-g
Dalam mod imbasan, cetak menggunakan gaya skrip shell
-i
Abaikan garis putus
--tulis ganti
Benarkan fail output untuk menulis ganti fail sedia ada
- membantu
Cetak ringkasan penggunaan
--verbose
Keluaran modul verbose
--senyap
Keluaran modul senyap
--ui
Paksa melancarkan dialog GUI
Parameter:
input=nama [diperlukan]
Fail ASCII yang mengandungi data input (atau "-" untuk dibaca daripada stdin)
output=nama [diperlukan]
Nama untuk peta raster output
kaedah=rentetan
Statistik untuk digunakan untuk nilai raster
Pilihan: n, min maksima, pelbagai, jumlah, bermaksud, stddev, kelainan, coeff_var, median,
persentil, kecondongan, trimmean
default: bermakna
jenis=rentetan
Jenis storan untuk peta raster terhasil
Pilihan: SEL, FCELL, DCELL
default: FCELL
pemisah=watak
Pemisah medan
Aksara khas: paip, koma, ruang, tab, baris baharu
default: paip
x=integer
Nombor lajur koordinat x dalam fail input (lajur pertama ialah 1)
default: 1
y=integer
Nombor lajur koordinat y dalam fail input
default: 2
z=integer
Bilangan lajur nilai data dalam fail input
Jika lajur nilai berasingan diberikan, pilihan ini merujuk kepada lajur koordinat z ke
ditapis oleh pilihan zrange
default: 3
skip=integer
Bilangan baris pengepala untuk dilangkau di bahagian atas fail input
default: 0
zrange=min, maks
Julat penapis untuk data z (min, maks)
zscale=terapung
Skala untuk digunakan pada data z
default: 1.0
lajur_nilai=integer
Bilangan lajur ganti nilai data dalam fail input
Jika tidak diberikan (atau ditetapkan kepada 0) data lajur-z digunakan
default: 0
vrange=min, maks
Julat penapis untuk data lajur nilai ganti (min, maks)
vscale=terapung
Skala untuk digunakan pada data lajur nilai ganti
default: 1.0
peratus=integer
Peratusan peta untuk disimpan dalam ingatan
Pilihan: 1-100
default: 100
pth=integer
pth persentil nilai
Pilihan: 1-100
mengurangkan=terapung
Buang peratus terkecil dan peratus daripada pemerhatian terbesar
Pilihan: 0-50
DESCRIPTION
. r.in.xyz modul akan memuatkan dan mengosongkan data ASCII x,y,z ke dalam peta raster baharu.
Pengguna boleh memilih daripada pelbagai kaedah statistik dalam mencipta raster baharu.
Data grid yang disediakan sebagai aliran mata x,y,z juga boleh diimport.
Sila ambil perhatian bahawa keluasan dan peleraian wilayah semasa digunakan untuk import. Ia adalah
oleh itu disyorkan untuk menggunakan terlebih dahulu -s bendera untuk mendapatkan keluasan mata input
diimport, kemudian laraskan rantau semasa dengan sewajarnya, dan hanya kemudian teruskan dengan
import sebenar.
r.in.xyz direka untuk memproses set data awan titik besar-besaran, contohnya LIDAR mentah atau
data swath sonar imbasan sisi. Ia telah diuji dengan set data sebesar puluhan bilion
mata (705GB dalam satu fail).
Statistik yang tersedia untuk mengisi raster ialah:
n bilangan titik dalam sel
minit nilai minimum mata dalam sel
maks nilai maksimum titik dalam sel
pelbagai julat titik dalam sel
jumlah jumlah titik dalam sel
bermakna nilai purata titik dalam sel
stddev sisihan piawai titik dalam sel
varians varians titik dalam sel
coeff_var pekali varians titik dalam sel
median nilai median titik dalam sel
persentil pth persentil titik dalam sel
skewness kecondongan titik dalam sel
trimmean min yang dipotong dalam sel
· varian dan terbitan menggunakan penganggar pincang (n). [tertakluk kepada perubahan]
· Pekali of varians diberikan dalam peratusan dan ditakrifkan sebagai (stddev/min)*100.
Ia juga mungkin untuk bin dan menyimpan lajur data lain (cth. backscatter) sambil
menapis dan menskala kedua-dua nilai lajur data dan julat z secara serentak.
NOTA
Berjalur data
Jika data diketahui berada pada grid biasa r.in.xyz boleh membina semula peta dengan sempurna sebagai
selagi berhati-hati diambil untuk menyediakan wilayah dengan betul dan peta asli data itu
unjuran digunakan. Kaedah biasa akan melibatkan penentuan resolusi grid sama ada
dengan memeriksa dokumentasi berkaitan data atau dengan mengkaji fail teks. Imbasan seterusnya
data dengan r.in.xyz's -s (Atau -g) bendera untuk mencari sempadan data input. GRASS menggunakan
konvensyen raster pusat sel di mana titik data berada di dalam pusat sel, sebagai
bertentangan dengan konvensyen nod grid. Oleh itu anda perlu mengembangkan rantau ini dengan
setengah sel ke semua arah melebihi apa yang ditemui imbasan dalam fail. Selepas wilayah
sempadan dan resolusi ditetapkan dengan betul dengan g.wilayah, jalankan r.in.xyz menggunakan n kaedah dan
sahkan bahawa n=1 di semua tempat. r.univar boleh tolong. Sebaik sahaja anda yakin bahawa rantau ini
betul-betul sepadan dengan data yang terus dijalankan r.in.xyz menggunakan salah satu bermaksud, min maks, Atau
median kaedah. Dengan n=1 keseluruhannya, hasilnya hendaklah sama tanpa mengira yang mana
kaedah tersebut digunakan.
Memori penggunaan
Sementara input fail boleh sewenang-wenangnya besar, r.in.xyz akan menggunakan sejumlah besar sistem
memori untuk kawasan raster yang besar (10000x10000). Jika modul enggan mula merungut
bahawa tidak ada memori yang mencukupi, gunakan peratus parameter untuk menjalankan modul dalam beberapa
pas. Selain menggunakan format peta yang kurang tepat (CELL [integer] atau FCELL [terapung
titik]) akan menggunakan kurang memori daripada DCELL [titik terapung berketepatan ganda] output peta
Kaedah seperti n, min maksima, jumlah juga akan menggunakan kurang memori, manakala stddev, kelainan, and
coeff_var akan menggunakan lebih banyak. Fungsi agregat median, persentil, skewness and dipotong
bermakna akan menggunakan lebih banyak memori dan mungkin tidak sesuai untuk digunakan dengan saiz yang besar
fail input.
Peta lalai jenis=FCELL bertujuan sebagai kompromi antara mengekalkan ketepatan data dan
mengehadkan penggunaan sumber sistem. Jika membaca data daripada aliran stdin, program
hanya boleh berjalan menggunakan satu hantaran.
Menetapkan rantau Sempadan and resolusi
Anda boleh menggunakan -s imbas bendera untuk mencari sejauh mana data input (dan dengan itu ketumpatan titik)
sebelum melaksanakan import penuh. guna g.wilayah untuk melaraskan sempadan rantau agar sepadan. The
-g bendera gaya shell mencetak takat yang sesuai sebagai parameter untuk g.wilayah. Sesuai
resolusi boleh didapati dengan membahagikan bilangan titik input dengan kawasan yang diliputi. cth
wc -l inputfile.txt
g.wilayah -p
# mata_setiap_sel = n_titik / (baris * kol)
g.wilayah -e
# lokasi UTM:
# mata_per_sq_m = n_mata / (ns_extent * ew_extent)
# Lokasi Lat/Lon:
# mata_per_sq_m = n_mata / (ns_extent * ew_extent*cos(lat) * (1852*60)^2)
Jika anda hanya berniat untuk menginterpolasi data dengan r.to.vect and v.surf.pertama, maka ada
sedikit perkara untuk menetapkan resolusi rantau dengan sangat baik sehingga anda hanya menangkap satu titik data
setiap sel -- anda juga boleh menggunakan "v.in.ascii -zbt" secara langsung.
Penapisan
Mata yang jatuh di luar kawasan semasa akan dilangkau. Ini termasuk mata jatuh
sebenarnya di kawasan selatan yang terikat. (untuk menangkap mereka melaraskan rantau dengan "g.region
s=s-0.000001"; lihat g.wilayah)
Baris kosong dan baris ulasan bermula dengan simbol cincang (#) akan dilangkau.
. zrange parameter boleh digunakan untuk menapis data input mengikut tahap menegak. Contoh
kegunaan mungkin termasuk menyediakan berbilang bahagian raster untuk digabungkan menjadi raster 3D
tatasusunan dengan r.to.rast3, atau untuk menapis outlier pada rupa bumi yang agak rata.
Dalam pelbagai rupa bumi pengguna mungkin mendapati itu minit peta menghasilkan penapis bunyi yang baik seperti kebanyakannya
Bunyi LIDAR adalah dari pukulan pramatang. The minit peta juga mungkin berguna untuk mencari asas
topografi dalam persekitaran berhutan atau bandar jika sel terlebih sampel.
Pengguna boleh menggunakan gabungan r.in.xyz output peta untuk mencipta penapis tersuai. cth penggunaan
r.mapcalc untuk mencipta peta min-(2*stddev). [Dalam contoh ini pengguna mungkin ingin memasukkan a
penapis sempadan bawah masuk r.mapcalc untuk mengeluarkan titik yang sangat berubah-ubah (kecil n) atau lari
r.jiran untuk melicinkan peta stddev sebelum digunakan selanjutnya.]
Ganti nilai ruang
. lajur_nilai parameter boleh digunakan dalam kes khusus apabila anda ingin menapis mengikut
julat-z tetapi tong dan simpan data lajur lain. Contohnya jika anda ingin melihat
nilai taburan belakang antara ketinggian 1000 dan 1500 meter. Ini amat berguna
semasa menggunakan r.in.xyz untuk menyediakan kepingan kedalaman untuk raster 3D — the zrange pilihan
mentakrifkan kepingan kedalaman tetapi nilai data yang disimpan dalam voxel menerangkan tambahan
dimensi. Seperti lajur z, julat penapisan dan faktor penskalaan boleh digunakan.
Unjuran semula
Jika peta raster ingin ditayangkan semula, mungkin lebih sesuai untuk mengunjurkan semula input
mata dengan m.proj or cs2cs sebelum berlari r.in.xyz.
Interpolasi ke dalam a DEM
Kebolehan topografi enjin vektor memperkenalkan overhed memori terhingga setiap vektor
titik yang biasanya akan mengehadkan peta vektor kepada kira-kira 3 juta mata (~ 1750^2
sel). Jika anda mahu lebih banyak, gunakan r.to.vect -b bendera untuk melangkau topologi bangunan. Tanpa
topologi, bagaimanapun, semua yang anda boleh lakukan dengan peta vektor dipaparkan dengan d.vect and
interpolasi dengan v.surf.pertama. Jalankan r.univar pada peta raster anda untuk menyemak bilangan
sel bukan NULL dan laraskan sempadan dan/atau resolusi seperti yang diperlukan sebelum meneruskan.
Arahan biasa untuk mencipta DEM menggunakan kesesuaian spline teratur:
r.univar lidar_min
r.to.vect -z type=point in=lidar_min out=lidar_min_pt
v.surf.rst in=lidar_min_pt elev=lidar_min.rst
Mengimport of x,y, rentetan data
r.in.xyz menjangkakan nilai berangka sebagai lajur z. Untuk melakukan kiraan kejadian
operasi walaupun pada data x,y dengan atribut bukan angka, data boleh diimport menggunakan
sama ada koordinat x atau y sebagai lajur z palsu untuk kaedah=n (kira bilangan mata setiap
sel grid), nilai z diabaikan juga.
CONTOH
Import kumpulan data Jockey's Ridge, NC, LIDAR (fail termampat "lidaratm2.txt.gz") dan
memprosesnya menjadi DEM yang bersih:
# imbas dan tetapkan sempadan rantau
r.in.xyz -s separator="," in=lidaratm2.txt out=test
g.wilayah n=35.969493 s=35.949693 e=-75.620999 w=-75.639999
g.region res=0:00:00.075 -a
# cipta peta "n" yang mengandungi kiraan mata setiap sel untuk memeriksa ketumpatan
r.in.xyz in=lidaratm2.txt out=lidar_n separator="," method=n zrange=-2,50
# ketumpatan titik semak [rho = n_sum / (baris*kol)]
r.univar lidar_n | jumlah grep
# buat peta "min" (ketinggian ditapis untuk hits pramatang)
r.in.xyz in=lidaratm2.txt out=lidar_min separator="," method=min zrange=-2,50
# tetapkan kawasan pengiraan kepada kawasan yang diminati
g.region n=35:57:56.25N s=35:57:13.575N w=75:38:23.7W e=75:37:15.675W
# semak bilangan sel bukan nol (cuba dan simpan di bawah beberapa juta)
r.univar lidar_min | grep '^n:'
# tukar kepada mata
r.to.vect -z type=point in=lidar_min out=lidar_min_pt
# interpolasi menggunakan padanan spline teratur
v.surf.rst in=lidar_min_pt elev=lidar_min.rst
# tetapkan skala warna kepada sesuatu yang menarik
r.colors lidar_min.peraturan pertama=bcyr -n -e
# sediakan versi berskala 1:1:1 untuk visualisasi NVIZ (untuk input lat/lon)
r.mapcalc "lidar_min.rst_scaled = lidar_min.rst / (1852*60)"
r.colors lidar_min.rst_scaled rule=bcyr -n -e
SEMUA
· Sokongan untuk keluaran peta berbilang daripada satu larian.
kaedah=rentetan[,rentetan,...] output=nama[,nama,...]
Ini boleh dikendalikan dengan mudah oleh skrip pembalut, dengan faedah tambahan daripadanya
sangat mudah untuk disejajarkan dengan cara itu.
· Tambah dua bendera baharu untuk sokongan untuk input binari langsung daripada libLAS untuk data LIDAR
dan mbio MB-System untuk data batimetri berbilang rasuk.
ambil perhatian: Lihat yang baru r.in.lidar modul untuk ini.
DIKENALI ISU
· n jumlah peta boleh menjadi lebih sedikit daripada `wc -l` dengan cth peratus=10 atau kurang.
Punca tidak diketahui.
· n peta peratus=100 dan peratus=xx peta berbeza sedikit (titik akan jatuh di atas/bawah
garis pembahagian)
Siasat dengan "r.mapcalc diff = bin_n.100 - bin_n.33" dsb.
Punca tidak diketahui.
· "nan" boleh bocor ke dalam coeff_var peta.
Punca tidak diketahui. Penyelesaian yang mungkin: "r.null setnull=nan"
Jika anda menghadapi sebarang masalah (atau penyelesaian!) sila hubungi Pasukan Pembangunan GRASS.
Gunakan r.in.xyzgrass dalam talian menggunakan perkhidmatan onworks.net