gdal_grid - Online sa Cloud

PROGRAMA:

NAME

gdal_grid - Lumilikha ang gdal_grid ng regular na grid mula sa nakakalat na data

SINOPSIS

gdal_grid [-ot {Byte/Int16/UInt16/UInt32/Int32/Float32/Float64/
CInt16/CInt32/CFloat32/CFloat64}]
[-of format] [-co "NAME=VALUE"]
[-zfield field_name] [-z_increase increase_value] [-z_multiply multiply_value]
[-a_srs srs_def] [-spat xmin ymin xmax ymax]
[-clipsrc |WKT|datasource|spat_extent]
[-clipsrcsql sql_statement] [-clipsrclayer layer]
[-clipsrcwhere expression]
[-l layername]* [-where expression] [-sql select_statement]
[-txe xmin xmax] [-tye ymin ymax] [-outsize xsize ysize]
[-a algorithm[:parameter1=value1]*] [-q]

DESCRIPTION

Ang program na ito ay lumilikha ng regular na grid (raster) mula sa nakakalat na data na nabasa mula sa OGR
pinanggalingan ng Datos. Ang data ng pag-input ay isasama upang punan ang mga grid node ng mga halaga, maaari kang pumili
mula sa iba't ibang paraan ng interpolation.

Simula sa GDAL 1.10, posibleng itakda ang GDAL_NUM_THREADS opsyon sa pagsasaayos
para iparallelize ang processing. Ang halagang tutukuyin ay ang bilang ng mga thread ng manggagawa, o
ALL_CPUS para gamitin ang lahat ng mga core/CPU ng computer.

-hindi na uri:
Para ang mga output band ay nasa ipinahiwatig na uri ng data.

-ng format:
Piliin ang format ng output. Ang default ay GeoTIFF (GTiff). Gamitin ang pangalan ng maikling format.

-txe xmin Xmax:
Itakda ang georeferenced X na mga lawak ng output file na gagawin.

-tye ymin ymax:
Itakda ang georeferenced Y na mga lawak ng output file na gagawin.

-malaki xsize ysize:
Itakda ang laki ng output file sa mga pixel at linya.

-a_srs srs_def:
I-override ang projection para sa output file. Ang srs_def maaaring alinman sa karaniwan
GDAL/OGR form, kumpletong WKT, PROJ.4, EPSG:n o isang file na naglalaman ng WKT.

-zfield field_name:
Tinutukoy ang field ng attribute sa mga feature na gagamitin para makakuha ng Z value. Ito
Ino-override ng value ang Z value na nabasa mula sa feature geometry record (natural, kung mayroon kang Z
halaga sa geometry, kung hindi, wala kang pagpipilian at dapat tukuyin ang isang pangalan ng field
naglalaman ng halaga ng Z).

-z_pagtaas pagtaas_halaga:
Pagdaragdag sa field ng attribute sa mga feature na gagamitin para makakuha ng Z value mula sa. Ang
Ang karagdagan ay dapat na kapareho ng yunit ng halaga ng Z. Ang halaga ng resulta ay magiging Z value + Z
dagdagan ang halaga. Ang default na halaga ay 0.

-z_multiply multiply_value:
Ito ay multiplication ratio para sa Z field. Ito ay maaaring gamitin para sa paglipat mula sa hal. paa hanggang
metro o mula sa taas hanggang sa malalim. Ang halaga ng resulta ay magiging (halaga ng Z + pagtaas ng Z
value) * Z multiply value. Ang default na halaga ay 1.

-a [algorithm[:parameter1=value1][:parameter2=value2]...]:
Itakda ang interpolation algorithm o pangalan ng sukatan ng data at (opsyonal) ang mga parameter nito.
Tingnan INTERPOLATION ALGORITHMS at DATA METRICS mga seksyon para sa karagdagang talakayan ng
Magagamit na Mga Pagpipilian.

- dumura xmin ymin Xmax ymax:
Nagdaragdag ng spatial na filter upang piliin lamang ang mga tampok na nasa loob ng kahon ng hangganan
inilarawan ni (xmin, ymin) - (xmax, ymax).

-clipsrc [xmin ymin Xmax ymax]|WKT|datasource|spat_extent:
Nagdaragdag ng spatial na filter upang piliin lamang ang mga feature na nasa loob ng tinukoy na hangganan
box (ipinahayag sa source SRS), WKT geometry (POLYGON o MULTIPOLYGON), mula sa isang
datasource o sa spatial na lawak ng - dumura opsyon kung gagamitin mo ang spat_extent
keyword. Kapag tumukoy ng datasource, karaniwang gusto mong gamitin ito
kumbinasyon ng -clipsrclayer, -cliprcwhere or -clipsrcsql mga pagpipilian.

-clipsrcsql sql_statement:
Piliin ang ninanais na mga geometry gamit ang isang SQL query sa halip.

-clipsrclayer layername:
Piliin ang pinangalanang layer mula sa source clip datasource.

-cliprcwhere pagpapahayag:
Limitahan ang mga gustong geometries batay sa query ng katangian.

-l layername:
Isinasaad ang (mga) layer mula sa datasource na gagamitin para sa mga feature ng pag-input. May
matukoy nang maraming beses, ngunit hindi bababa sa isang pangalan ng layer o a -sql ang pagpipilian ay dapat na
tinukoy.

-saan pagpapahayag:
Isang opsyonal na SQL WHERE style query expression na ilalapat upang pumili ng mga feature
proseso mula sa (mga) input layer.

-sql piliin ang_pahayag:
Isang SQL statement na susuriin laban sa datasource para makagawa ng virtual na layer ng
mga tampok na ipoproseso.

-co 'NAME=VALUE':
Nagpapasa ng opsyon sa paggawa sa driver ng format ng output. Maramihan -co ang mga pagpipilian ay maaaring
nakalista. Tingnan ang dokumentasyong tukoy sa format para sa mga opsyon sa legal na paggawa para sa bawat format.

-q:
Pigilan ang progress monitor at iba pang non-error na output.

src_datasource:
Sinusuportahan ng anumang OGR ang nababasang datasource.

dst_filename:
Sinusuportahan ng GDAL ang output file.

INTERPOLATION ALGORITHMS

Mayroong bilang ng mga interpolation algorithm na mapagpipilian.

invdist
Baligtad na distansya sa isang kapangyarihan. Ito ay default na algorithm. Mayroon itong mga sumusunod na parameter:

kapangyarihan:
Lakas ng pagtimbang (default 2.0).

makinis:
Smoothing parameter (default 0.0).

radius1:
Ang unang radius (X axis kung ang anggulo ng pag-ikot ay 0) ng ellipse ng paghahanap. Itakda ang parameter na ito
sa zero upang magamit ang buong point array. Ang default ay 0.0.

radius2:
Ang pangalawang radius (Y axis kung ang anggulo ng pag-ikot ay 0) ng ellipse ng paghahanap. Itakda ito
parameter sa zero upang magamit ang buong point array. Ang default ay 0.0.

anggulo:
Anggulo ng pag-ikot ng ellipse ng paghahanap sa mga degree (counter clockwise, default 0.0).

max_points:
Maximum na bilang ng mga data point na gagamitin. Huwag maghanap ng mas maraming puntos kaysa sa numerong ito.
Ito ay ginagamit lamang kung ang search ellipse ay nakatakda (parehong radii ay hindi zero). Zero ang ibig sabihin niyan
lahat ng nahanap na punto ay dapat gamitin. Ang default ay 0.

min_points:
Minimum na bilang ng mga data point na gagamitin. Kung mas kaunting halaga ng mga puntos ang natagpuan ang grid node
itinuturing na walang laman at pupunan ng NODATA marker. Ito ay ginagamit lamang kung naghahanap
ellipse ay nakatakda (parehong radii ay hindi zero). Ang default ay 0.

nodata:
NODATA marker upang punan ang mga walang laman na puntos (default 0.0).

karaniwan
Moving average na algorithm. Mayroon itong mga sumusunod na parameter:

radius1:
Ang unang radius (X axis kung ang anggulo ng pag-ikot ay 0) ng ellipse ng paghahanap. Itakda ang parameter na ito
sa zero upang magamit ang buong point array. Ang default ay 0.0.

radius2:
Ang pangalawang radius (Y axis kung ang anggulo ng pag-ikot ay 0) ng ellipse ng paghahanap. Itakda ito
parameter sa zero upang magamit ang buong point array. Ang default ay 0.0.

anggulo:
Anggulo ng pag-ikot ng ellipse ng paghahanap sa mga degree (counter clockwise, default 0.0).

min_points:
Minimum na bilang ng mga data point na gagamitin. Kung mas kaunting halaga ng mga puntos ang natagpuan ang grid node
itinuturing na walang laman at pupunan ng NODATA marker. Ang default ay 0.

nodata:
NODATA marker upang punan ang mga walang laman na puntos (default 0.0).

Tandaan, na mahalagang itakda ang ellipse ng paghahanap para sa moving average na paraan. Isa itong bintana
na i-a-average kapag nag-compute ng mga halaga ng mga grid node.

pinakamalapit na
Pinakamalapit na algorithm ng kapitbahay. Mayroon itong mga sumusunod na parameter:

radius1:
Ang unang radius (X axis kung ang anggulo ng pag-ikot ay 0) ng ellipse ng paghahanap. Itakda ang parameter na ito
sa zero upang magamit ang buong point array. Ang default ay 0.0.

radius2:
Ang pangalawang radius (Y axis kung ang anggulo ng pag-ikot ay 0) ng ellipse ng paghahanap. Itakda ito
parameter sa zero upang magamit ang buong point array. Ang default ay 0.0.

anggulo:
Anggulo ng pag-ikot ng ellipse ng paghahanap sa mga degree (counter clockwise, default 0.0).

nodata:
NODATA marker upang punan ang mga walang laman na puntos (default 0.0).

DATA METRICS

Bukod sa interpolation functionality gdal_grid ay maaaring gamitin upang kalkulahin ang ilang sukatan ng data
gamit ang tinukoy na window at output grid geometry. Ang mga sukatan na ito ay:

pinakamaliit:
Pinakamababang halaga na makikita sa grid node na ellipse ng paghahanap.

pinakamataas:
Maximum na halaga na natagpuan sa grid node search ellipse.

saklaw:
Isang pagkakaiba sa pagitan ng minimum at maximum na mga halaga na makikita sa grid node search ellipse.

bilangin:
Isang bilang ng mga data point na matatagpuan sa grid node search ellipse.

average_distance:
Isang average na distansya sa pagitan ng grid node (gitna ng search ellipse) at lahat ng
ang mga data point na matatagpuan sa grid node search ellipse.

average_distance_pts:
Isang average na distansya sa pagitan ng mga data point na makikita sa grid node search ellipse. Ang
ang distansya sa pagitan ng bawat pares ng mga puntos sa loob ng ellipse ay kinakalkula at average ng lahat
itinakda ang mga distansya bilang halaga ng grid node.

Ang lahat ng mga sukatan ay may parehong hanay ng mga opsyon:

radius1:
Ang unang radius (X axis kung ang anggulo ng pag-ikot ay 0) ng ellipse ng paghahanap. Itakda ang parameter na ito
sa zero upang magamit ang buong point array. Ang default ay 0.0.

radius2:
Ang pangalawang radius (Y axis kung ang anggulo ng pag-ikot ay 0) ng ellipse ng paghahanap. Itakda ito
parameter sa zero upang magamit ang buong point array. Ang default ay 0.0.

anggulo:
Anggulo ng pag-ikot ng ellipse ng paghahanap sa mga degree (counter clockwise, default 0.0).

min_points:
Minimum na bilang ng mga data point na gagamitin. Kung mas kaunting halaga ng mga puntos ang natagpuan ang grid node
itinuturing na walang laman at pupunan ng NODATA marker. Ito ay ginagamit lamang kung naghahanap
ellipse ay nakatakda (parehong radii ay hindi zero). Ang default ay 0.

nodata:
NODATA marker upang punan ang mga walang laman na puntos (default 0.0).

Pagbabasa PARAGRAPH HIWALAY Mga halaga

Kadalasan mayroon kang isang text file na may listahan ng mga halaga ng XYZ na pinaghihiwalay ng kuwit upang gumana (kaya
tinatawag na CSV file). Madali mong magagamit ang ganoong uri ng data source sa gdal_grid. Lahat ng kailangan mo
ay gumawa ng virtual dataset header (VRT) para sa iyong CSV file at gamitin ito bilang input datasource
para gdal_grid. Makakakita ka ng mga detalye sa format ng VRT sa pahina ng paglalarawan ng Virtual Format.

Narito ang isang maliit na halimbawa. Hayaan kaming magkaroon ng isang CSV file na tinatawag dem.csv Na naglalaman ng

Silangan, Hilaga, Elevation
86943.4,891957,139.13
87124.3,892075,135.01
86962.4,892321,182.04
87077.6,891995,135.01

Para sa data sa itaas ay gagawa kami dem.vrt header na may sumusunod na nilalaman:



dem.csv
wkbPoint




Tinutukoy ng paglalarawang ito ang tinatawag na 2.5D geometry na may tatlong coordinate X, Y at Z. Z
ang halaga ay gagamitin para sa interpolation. Ngayon ay maaari mong gamitin dem.vrt kasama ang lahat ng OGR program (start
sa oginfo upang subukan na ang lahat ay gumagana nang maayos). Maglalaman ang datasource ng isang layer
tinatawag 'dem' puno ng mga tampok na punto na binuo mula sa mga halaga sa CSV file. Gamit ito
diskarteng maaari mong pangasiwaan ang mga CSV file na may higit sa tatlong column, lumipat ng column, atbp.

Kung ang iyong CSV file ay hindi naglalaman ng mga header ng column, maaari itong pangasiwaan sa mga sumusunod
paraan:



Ang pahina ng paglalarawan ng Comma Separated Value ay naglalaman ng mga detalye sa format na CSV na sinusuportahan ng
GDAL/OGR.

Halimbawa

Ang sumusunod ay lilikha ng raster TIFF file mula sa VRT datasource na inilarawan sa Pagbabasa PARAGRAPH
HIWALAY Mga halaga seksyon gamit ang kabaligtaran na distansya sa isang paraan ng kapangyarihan. Mga halaga sa
mababasa ang interpolate mula sa Z value ng geometry record.

gdal_grid -a invdist:power=2.0:smoothing=1.0 -txe 85000 89000 -tye 894000 890000 -outsize 400 400 -of GTiff -ot Float64 -l dem dem.vrt dem.tiff

Ang susunod na utos ay gumagawa ng parehong bagay tulad ng nauna, ngunit nagbabasa ng mga halaga upang i-interpolate
mula sa field ng katangian na tinukoy sa -zfield opsyon sa halip na geometry record. Kaya sa
sa kasong ito, ang mga coordinate ng X at Y ay kinukuha mula sa geometry at ang Z ay kinuha mula sa
'Elevation' patlang. Ang GDAL_NUM_THREADS ay nakatakda rin upang iparallelize ang pag-compute.

gdal_grid -zfield "Elevation" -a invdist:power=2.0:smoothing=1.0 -txe 85000 89000 -tye 894000 890000 -outsize 400 400 -of GTiff -ot Float64 -l dem_CPDALffNUM dem.

MGA AUTHORS

Andrey Kiselev [protektado ng email]

Gumamit ng gdal_grid online gamit ang mga serbisyo ng onworks.net