Daha çox

GDAL istifadə edərək bir proyeksiya faylını necə yarada bilərəm

GDAL istifadə edərək bir proyeksiya faylını necə yarada bilərəm


Aşağıdakı funksiya ilə GDAL istifadə edərək bir polyline faylından bir tampon əldə etməyə çalışıram

def createBuffer (inputfn, outputBufferfn, bufferDist): inputds = ogr.Open (inputfn) inputlyr = inputds.GetLayer () shpdriver = ogr.GetDriverByName ('ESRI Shapefile') os.path.exists (outputBufferSDete): shpdata outputBufferfn) outputBufferds = shpdriver.CreateDataSource (outputBufferfn) bufferlyr = outputBufferds.CreateLayer (outputBufferfn, geom_type = ogr.wkbPolygon) featureDefn = bufferlyr.GetLayerDefer ()) ) outFeature = ogr.Feature (featureDefn) outFeature.SetGeometry (geomBuffer) bufferlyr.CreateFeature (outFeature)

Yaxşı işləyir və .bdf .shp .shx faylları yaradır

Yaxşı. Amma eyni qovluğa (funksiya tərəfindən yaradılmış) projectioin faylını daxil etmək istərdim .prj

Aşağıdakı funksiyaya sahibəm.

def construye_proyeccion (nombre_arhivo): spatialRef = osr.SpatialReference () spatialRef.ImportFromEPSG (4326) spatialRef.MorphToESRI () file = open (' % s.prj' % fayl adı, 'w') fileRef.wite (spatialRef = osr. ) file.close ()

Tamponu və onun .prj faylını eyni qovluqda yaratmaq üçün nə edə bilərəm?


İki şəkil faylının proyeksiyasının eyni olduğunu düşündüyüm iki funksiyanı birləşdirə bilərsiniz.

1)spatialRef = inputlyr.GetSpatialRef ()orijinal şəkil faylının proyeksiyasını verir
2)prj = os.path.splitext (outputBufferfn) [0] + ".prj"tampon şəkil faylının prj faylının adını verir

outputBufferfn = "a_shapefile.shp" os.path.splitext (outputBufferfn) [0] + ".prj" 'a_shapefile.prj'

Yeni

outputBufferfn = "/users/shared/a_shapefile.shp" (filepath, fayl adı) = os.path.split (outputBufferfn) os.path.join (filepath, os.path.splitext (fayl adı) [0] + ".prj" ) '/users/shared/a_shapefile.prj'

Sonra:

def createBuffer (inputfn, outputBufferfn, bufferDist):… # spatialRef = inputlyr.GetSpatialRef () prj = os.path.splitext (outputBufferfn) [0] + ".prj" file = open (prj, 'w') fayl .write (spatialRef.ExportToWkt ()) file.close ()…

Shapely və Fiona (GDAL/OGR kitabxanasının başqa bir Python sarğısı) ilə ən asan və avtomatikdir (prj faylı):

şekli.geometrinin idxal şəklindən, giriş olaraq fiona.open ("a_shapefile.shp") ilə xəritələndirmə: # yalnız şəkil faylının həndəsəsini Polygon input.schema ['geometry'] = "Polygon" # yazın fiona.open ilə tampon shapefile ('buffer_a_shapefile.shp', 'w', 'ESRI Shapefile', input.schema, input.crs) kimi çıxış üçün: girişdə feat üçün: # GeoJSONdan düzgün geometriyə geom = forma (feat [' həndəsə]]) # GeoJSON feat ['geometry'] üçün # formalı tampon = mapping (geom.buffer (5)) output.write (feat)

Bəzi izahatlar:

1) orijinal şəkil faylının sxemi (giriş.schema, JSON/Python lüğət)

{'geometry': 'Point', 'properties': OrderedDict ([(u'name ',' str: 254 ')])}

2) orijinal şəkil faylının crs

input.crs {'init': u'epsg: 4326 '}

3) tamponlanmış şəkil faylının sxemi:

input.schema ['geometry'] = "Polygon" input.schema {'geometry': 'Polygon', 'properties': OrderedDict ([(u'name ',' str: 254 ')])}

Müxtəlif Dəstək GDAL Raster Formatları

Proqnozlar dəstəyi ( *.prj faylı varsa oxuyun) olduqca məhduddur. Əlavə nümunə .prj faylları, saxlayıcıya, [email protected] ünvanına göndərilə bilər.

Şəbəkə oxunması üçün NODATA dəyəri də mövcud olduqda qorunur. İlk tarama xəttində ondalık dəyərlər olmadıqda, ızgaralar imzalı 16 bit tamsayı kimi qəbul edilir, bu halda görüntü 32 bit üzən nöqtə kimi qəbul edilir.

AIG GIO sürücüləri Arc/Info Binary Grid formatında da mövcuddur.

Qeyd: kimi tətbiq olunur gdal/frmts/aaigrid/aaigriddataset.cpp.


GDAL - Coğrafi İnformasiya Sistemlərindən istifadə edərək bir proyeksiya faylını necə yarada bilərəm

Chris Goodell tərəfindən yazılmışdır | 14 may 2015

RAS Mapper, 2D sahələri əvvəlcədən işlətməyə, nəticələrini xəritələməyə və arxa plan şəkillərini idarə etməyə imkan verən HEC-RAS-dakı bir pəncərədir. HEC-RAS-ın gələcək versiyalarında RAS Mapper, HEC-RAS modelləşdirmə təcrübəmizdə getdikcə daha çox yayılacaq. Anladığım kimi, nəticədə RAS Mapper və həndəsə redaktoru pəncərəsi simvolizə olunan "Əsas RAS Pəncərəsi" ni əvəz edərək HEC-RAS üçün ön interfeys yaratmaq üçün birləşə bilər.

Amma bu fərziyyədir. HEC-RAS 5.0-ı çox sadə bir modeldən başqa bir şey üçün istifadə etmək istəyirsinizsə, həqiqətən RAS Mapper-dən istifadə etməyiniz və uzantı ilə ... Projeksiyon Faylının olması fərziyyə deyil. Layihə faylı xüsusi bir coğrafi koordinat sistemini və proyeksiyasını təyin edir və bizim HEC-RAS modelləşdiriciləri üçün bir qədər yeni bir şeydir (əslində yarandığı gündən RAS Mapper-in bir hissəsi idi, lakin indi Version 5.0-də 2D modelləşdirmə və veb-görüntü ilə birlikdə) , RAS Mapper HEC-RAS modelləşdirməsinin ayrılmaz hissəsinə çevrilir). Layihə faylı, RAS istifadəçiləri olaraq əvvəllər çox şey bilmədiyimiz coğrafi xəritələşdirmə (yəni GIS) haqqında bir şey bilməyimizi tələb edir. Ancaq RAS Mapper -dən istifadə etmək istəyirsinizsə, və edəcəksən, proyeksiya faylının nə olduğunu və onu necə əldə edəcəyinizi başa düşməlisiniz. Onsuz RAS Mapper çox yararsızdır. RAS layihəmizə bir proyeksiya faylı əlavə etmək, layihənin coğrafi istinadını təyin edir. Projeksiyon faylları *.prj uzantısına malikdir. Bunu HEC-RAS layihə faylı ilə qarışdırmamaq üçün diqqətli olun *.prj. Hər ikisi də HEC-RAS layihə qovluğunda yerləşir, ancaq proyeksiyanızı təyin etmək üçün RAS Mapper-də yalnız düzgün formatlanmış bir proyeksiya faylı işləyəcək.

Layihə faylı, həqiqətən müəyyən bir formatda açar sözləri olan sadə bir mətn faylıdır. Əslində, "Tanınmış Mətn" formatında və ya WKT ilə yazılmış tək bir simdir. Sadə, oxumaq asandır və tərəfindən yaradılmışdır GIS Konsorsiumunu açın. İçəridə bir proyeksiya faylının necə göründüyü budur. Diqqət yetirin ki, bütün CAPS tərəfindən müəyyən edilmiş bəzi açar sözlər var, ardınca mötərizədə olan açar sözlə əlaqəli bəzi məlumatlar var []. Mən nəyin nə ilə getdiyini görməyi asanlaşdırmaq üçün onu rənglə kodladım. Bənövşəyi rəng iyerarxiyadakı ən yüksək nizamı ifadə edir, sonra mavi, sonra yaşıl, sonra qırmızı. Başqa sözlə, qırmızı açar söz yaşıl söz üçün "uşaq", yaşıl mavi üçün uşaq, mavi bənövşəyi uşaqdır.

Aşağıdakı hər bir açar söz və mötərizəli məlumatlar "bənd" adlanır. Birinci və əsas bənd, PROJCS, Proqnozlaşdırılan Koordinat Sistemini ifadə edir. Proqnozlaşdırılan koordinat sistemi aşağıdakı yarımbəndlərdən ibarətdir:

1. Enlem və uzunluq dərəcələrinə əsaslanan və üfüqi istinad datumu (DATUM) və uzunluq ölçüləri üçün istinad meridianını (PRIMEM) ehtiva edən Coğrafi Koordinat Sistemi (GEOGCS). DATUM da yuxarıdakı nümunədə 1866 -cı il Clarke Ellipsoid olan yerin formasının (SPHEROID) təsvirini ehtiva edir. Vahidlər (UNIT) burada yalnız GEOGCS, bu halda dərəcələr üçün çıxarılır.
2. Coğrafi koordinatlardan (lat/uzun) proqnozlaşdırılan koordinatlara proyeksiyası olan proyeksiya (PROJEKSİYA). Üç ölçülü sferoidin (Yerin) iki ölçülü görüntü mühitinə necə yansıtıldığı mahiyyətcə budur. Yuxarıdakı nümunədə, Transvers Merkator istifadə edən seçilir Universal Transvers Merkator (UTM) koordinat sistemi.

3. Müxtəlif proyeksiya parametr dəyərləri (PARAMETER). Parametr kotirovkalarda etiketlənir, ardınca isə dəyəri verilir.
4. Proqnozlaşdırılan koordinat sistemi (UNIT) üçün vahidlər. Burada sayğaclar çevrilmə əmsalı “1” olan xətti vahid kimi istifadə olunur. Dönüşüm əmsalı təsvir olunan vahidləri sayğaclara çevirir. "FOOT_US" istifadə edilərsə, dönüşüm əmsalı 0.30480060960121924 olardı.

Layihə sənədinizdə bəzi əlavə müddəalar ola bilər, ancaq yuxarıda sadalananlar tipik görünür. Təsvirləri ilə WKT formatında istifadə olunan bütün açar sözləri GeoAPI -də tapa bilərsiniz burada.
Öz proyeksiya sənədlərinizi yaza bilərsiniz və GeospatialPython.com bir üsul təqdim edir (başqa nümunələr var, yalnız Google). Bununla birlikdə, artıq tərtib edilmiş bir proyeksiya faylını tapmaq və ondan istifadə etmək daha asan və daha praktikdir. Artıq coğrafi istinadlı bir HEC-RAS layihəniz varsa, həndəsə komponentlərinizi (axın mərkəz xətti, xscutlines, axın xətləri və s.) Kompüterinizdə harada saxlandığını tapın və istifadə edin.

Coğrafi istinad HEC-RAS layihənizi yaratmaq üçün istifadə olunan GIS sənədləriniz yoxdursa, ən azı üç fərqli yerdə proyeksiya faylları tapa bilərsiniz: ArcGIS 10.0 və ya daha əvvəl, spatialreference.org və EPSG Projection Database. Başqalarını tanıyırsınızsa, aşağıda şərh yazın!

Bu mənbələrdən hər hansı birini istifadə edərkən, doğru proyeksiya faylını seçdiyinizə əmin olun. Doğru olub olmadığını anlayacaqsınız veb görüntülərini RAS Mapper -ə gətirmək və hər şeyin məkan baxımından düzgün olduğunu yoxlamaq. ArcGIS (10.1 -dən əvvəlki versiyalar) 5000 -dən çox coğrafi, proqnozlaşdırılan və şaquli koordinat sistemini ehtiva edən bir Koordinat Sistemləri qovluğunu ehtiva edir. Təəssüf ki, ArcGIS -in daha yeni versiyaları bu qovluqla birlikdə gəlmir. ArcGIS 10.1 və ya daha yeni bir versiyanız varsa və ya ümumiyyətlə ArcGISiniz yoxdursa, böyük bir məkan istinad sistemləri məlumat bazasına daxil ola bilərsiniz. http://spatialreference.org/

Spatialreference.org saytına daxil olsanız, yalnız EPSG, IAU2000 və ya spatialreference.org istinadlarından seçim etdiyinizə əmin olun. ESRI istinadlarında yalnız PROJCS bəndini deyil, proyeksiya faylında GEOGCS maddəsi var. RAS Mapper -də işləməyəcəklər. İstədiyiniz istinadı tapdıqdan sonra üzərinə vurun və sonra mövcud formatların siyahısından ".PRJ Faylını" seçin. Bir proyeksiya faylı kompüterinizə yüklənəcək və RAS Mapper -də istifadə etməyə hazırsınız. İstinad üçün açar sözləri axtarmağa imkan verən rahat bir axtarış qutusu var. Məsələn, layihəniz Havaydadırsa və üfüqi datumun NAD83 olduğunu bilirsinizsə, açar sözlər daxil edə bilərsiniz: Hawaii NAD83, Axtarış düyməsini basın və aşağıdakı məkan istinadlarının siyahısını görəcəksiniz:

Geospatialpython.org tərəfindən GoogleCode'da yerləşdirilən EPSG Projeksiyon Veritabanı da onlayn olaraq mövcuddur. Burada asanlıqla kopyalanıb öz proyeksiya sənədinizə yapışdırıla bilən düzgün WKT formatında bir çox proyeksiya faylının mətn faylını tapa bilərsiniz.

Bir proyeksiya faylını seçib RAS Mapper-ə təyin etdikdən sonra, düzgün olduğunu iki dəfə yoxlayın veb görüntü əlavə etmək. Hər şey sıralanırsa, getməyiniz yaxşıdır. Aşağıdakı şəkil, RAS Mapper -də səhv bir proyeksiya faylı təyin edilmiş Muncie verilənlər bazasını göstərir.

Aydındır ki, bu nümunədəki axınlarımız və kəsişmələrimiz aid olduqları Muncie Indianada düzgün uyğunlaşdırılmamışdır. Yanlış proyeksiya faylı bizi Kanadanın Alberta şəhərinə gətirdi!

Düzgün proyeksiya faylını yenidən təyin etmək, modeli düzgün məkan istinadına qaytarır.

HEC-RAS-da istifadə edə biləcəyimiz başqa proyeksiya faylları haqqında hər kəs bilirmi? Əgər belədirsə, aşağıda şərh yazın.

Şərhlər

Anonim

Blogunuzdakı proqnozlar haqqında yazınızı gördüm və modelin koordinat sistemini HEC-RAS-da indi çevirə biləcəyinizi qeyd etmək yaxşı ola biləcəyini düşündüm.

Chris Goodell

Cam Ackerman

HEC-RAS, PRJ-in esri WKT formatını dəstəkləyir. Digər WKT formatlarından istifadə YMMV -dir.

Cam Ackerman

HEC-RAS hazırda esri 's WKT formatını dəstəkləyir. Digər WKT formatları YMMV -dir.

Chris Goodell

Başınızı qaldırdığınız üçün təşəkkürlər. İnanıram ki, spatialreference.org saytındakı bütün prj faylları, 447 ESRI istinadları istisna olmaqla (istehza ilə) HEC-RAS tərəfindən istifadə olunan formata uyğun gəlir (ən azı yoxladıqlarım). Onların yalnız GEOGCS maddəsi var, digərləri yoxdur.

Anonim

RasMapper -də öz WMS linkimi yükləyə bilərəmmi?

Çünki mənim modelim S-JTSK / Krovak East North və WGS-84-də bacgkround / WMS (ortofoto, peyk görüntüsü) yaradılmışdır. Bir bacgkround yüklədiyim zaman modelim yerindən çıxdı- bu səhv proqnozlaşdırıla bilər.

Chris Goodell

Modelinizin yaradıldığı proyeksiya və RAS Mapper -də yüklədiyiniz proyeksiya eyni olmalıdır.

Robert S.

Kris fikrimi təsdiqlədi:
Digər koordinat proqnozlarını WGS84 olaraq işləyərkən WMS görüntü serverlərindən (RasMapper- Tools) istifadə edə bilməzsiniz. Bu gün RasMapper digər WMS -ə qoşulmağı dəstəkləmir.
Statik görüntü (coğrafi istinad şəkli) istifadə edərək bir fon həll etməlisiniz.

Ümid edirəm ki, yaxşı yazıram, Chris!

Chris Goodell

Robert. Sualınızı daha yaxşı başa düşdüyümü düşünürəm (Kameronun YÖK -də köməyi ilə). HEC-RAS layihə kataloquna girsəniz, GDAL adlı bir qovluq görəcəksiniz. Orada Veb Xəritə Xidmətləri adlı bir qovluq görürsünüz. Orada HEC-RAS-da mövcud olan hər bir veb görüntü xidməti üçün bir xml faylını görəcəksiniz. Daha çox veb görüntü xidməti əlavə etmək istəyirsinizsə, öz xml -nizi yarada və bu qovluğa yerləşdirə bilərsiniz. Daha sonra HEC-RAS-da Veb Görüntü Pəncərəsində görünməlidir. Bunu öz veb görüntülərinizlə edə bilirsinizsə, xml -i mənə göndərin ki, başqaları ilə paylaşaq.

Mike Chilson

ArcGIS 10.1+ ilə sadəcə sevdiyiniz siyahıya bir proyeksiya əlavə edin və bir proyeksiya faylı yaradacaq və yerləşdirəcək:
C: İstifadəçilər & İstifadəçi Adı & quotDataRoamingESRIDesktop10.2ArcMapCoordinate Sistemləri

Gerrit K.

Bu mesajın hər hansı bir sirli hissəsinə görə çox üzr istəyirəm, ancaq Macar EOV -un (Vahid Milli Projeksiyon) RAS Mapper -də niyə düzgün yansıtılmadığını araşdırdıqdan bir neçə saat sonra, bunu düzgün etmək üçün bütün məlumatların mövcud olduğunu aşkar etdim. RAS Mapper -in istifadə etdiyi GDAL kitabxanası.

Səhv proyeksiya əvvəlcə məlumat mənbəsindəki bir parametr dəyərindəki bir səhvdən qaynaqlana bilər, lakin bu çoxdan düzəldilmişdir və düzəliş GDAL kitabxana faylının bir hissəsidir 'datum_shift.csv '.

Bu faylda sadalanan 5 (!) Datum shift parametrləri var (& quot; SEQ_KEY & quot 350 & quot; SEQ_KEY & quot 354 -ə qədər) və bunların hər hansı biri HD72 -dən WGS 84 -ə daha cari tətbiqdən daha dəqiq keçmək olar.

Beləliklə, GDAL kitabxanası yalnız qismən tətbiq olunurmu? Yoxsa RAS Mapper -in GDAL kitabxana faylında və#39datum_shift.csv ' -da təyin olunan üstünlük verilən keçidini tətbiq etməsinə mane olan bir səhv varmı?

Chris Goodell

Gerrit- məlumat üçün təşəkkürlər. Problemin nə ola biləcəyindən əmin deyiləm. Yaxşı bir cavab ala biləcəyimizi görmək üçün başqalarına sualınıza işarə edəcəyəm.

Salam Chris, bu yazı üçün təşəkkürlər. Həndəsəmin WMS çayından/ərazisindən kənara çıxması ilə bağlı problemim var (aşağıdakı linkə baxın). .Prj faylını əl ilə düzəltməklə bir həll yolu tapdım, amma bu 2D modelin proyeksiyasına təsir edəcək. sel xəritələri?

Chris Goodell

Salam Mike. Yaxşı sualdır. Bunun necə işləyəcəyindən əmin deyiləm. Bağırsaq instinktim yaxşı işləməli olduğunu söyləyir, amma düşünürəm ki, bu, dəqiq bilmək üçün sınamalı olduğunuz şeylərdən biridir. Necə getdiyini mənə bildir.
Kris

Naməlum

Layihə faylı yalnız RAS-Mapper və ya həndəsə üçün əhəmiyyətlidirmi? İşlədiyi güman edilən bir RAS modelini miras aldım, hətta son iş nəticələrini göndərdilər. Bir həndəsə səhvi səbəbiylə işlədə bilmirəm –, bəzi çayların əlaqəli kəsişmələrinə düzgün düzülməmiş görünür. Bəzi qəribələr ərazi yüksəkliyinin altına düşür. Düşünürəm ki, bu problemlər proyeksiya faylının olmaması səbəbindən ola bilər.

Chris Goodell

Həndəsə quruluşundakı səhvlər kimi səslənir. Proqnozlaşdırma ilə heç bir əlaqəsi olmamalıdır.

Colin Holdgrafer

RAS Mapper-in bir UTM Zone 15 NAD83 (metr) proyeksiya faylı istifadə edərək xəritəni çox gözəl göstərdiyi bir layihəm var, amma Options- & gtConvert Project Units-dən istifadə edərək layihəni metrdən ayağa çevirdikdə xəritə artıq işləmirdi. Yeni layihə üçün proyeksiya faylımı necə dəyişə bilərəm ki, layihəni lazımi yerdə göstərsin? Yoxsa vahid çevrilməsinə əsaslanan yeni bir proyeksiya faylını yükləmək daha məntiqli olardı? Və belədirsə, bu yeni üfüqi məlumat nə olardı?

Chris Goodell

Doğru vahid sisteminə malik yeni bir proyeksiya faylı əldə etməlisiniz (və ya özünüz yaratmalısınız).

Colin Holdgrafer

İş ondadır ki, proyeksiya hələ də UTM Zone 15 NAD83 (metr) dir, ancaq indi ayaqlara çevrildi.

Bir kəsişmənin çayla kəsişdiyi müəyyən bir nöqtədə vahid çevrilməsindən sonra koordinatlar necə dəyişdi:

(Original) Metrlər: 451305.5383, 4625115.6106
(Dönüşümdən sonra) Ayaqlar: 1480661.2147, 15174263.8183

Dönüştürmədən əvvəl istifadə etdiyim forma proyeksiya faylı:

Layihəmin vahid çevrilməsindən sonra Arktik Okeanının ortasında bitməməsi üçün bunu necə dəyişdirmək lazımdır? UNIT [& quotAyaq & quot, 3.28084]] demək üçün UNIT [& quotMetre & quot, 1]] dəyişdirməyə çalışdım, amma bu problemi həll etmədi.

Ben Snyder

Mən uzun müddətdir ki, RAS -ın etdiyi sirli üfüqi məlumat dəyişməsini həll etməyə çalışıram. CA State Plane Zone 1, US Feet, EPSG 2225 istifadə edərək Civil 3D-dən bir səth ixrac edirəm. Ən son işləmə testim üçün Civil 3D tif ilə eyni proyeksiya məlumatlarına malik RAS üçün bir proyeksiya faylı hazırladım. Yeni bir ərazi quranda & quot; Eyni Layihə & quot; yazır. Sonra həndəsə redaktorunu açıb bir işarənin koordinatlarına baxdığımda, şərq istiqamətdə təxminən 12 fut, şimalda isə bir neçə metr uzaqlaşır. TİF -i CİS -ə gətirdiyim zaman dəyişiklik YOXDUR, ona görə də Mülki 3D ixracatında problem kimi görünmür. Mən məğlub oldum! Hər hansı bir fikir? Problem yaradır, çünki bütün RAS nəticələrinin rasterləri dəyişdirilir və sıralanmır.

Ben Snyder

Əslində, bu işlədi! Civil 3D ilə işləyən RAS istifadəçiləri, sirr dəyişikliyinin baş verməməsi üçün SRS -in & quot; Layihə Eyni & quot; olmasını təmin etməli ola bilərlər. Ancaq tam olaraq eyni olduqları müddətcə getmələri yaxşıdır.

Anonim

Salam, köməyinizə ehtiyacım var. Prj faylını: & quothttp: //spatialreference.org/&quot -dan işləyə bilmirəm. & QuotEPSG: 22174: POSGAR 98 / Argentina 4 & quot üçün ehtiyacım var
Səhv tapa bilmirəm. Çox təşəkkür edirəm.

Anonim

HEC-RAS layihəsi aşağıdakı kimi qurulur:

• Vahidlər: ABŞ adətidir
• LiDAR verilənlər bazası
o Dövlət təyyarəsi proyeksiyası
o Yatay vahidlər = beynəlxalq ayaqlar
• RAS Mapper -də istifadə olunan proyeksiya
o Dövlət təyyarəsi proyeksiyası
o Yatay vahidlər = metr

Hava görüntüləri düzgün sıralanır. Ölçmə vasitəsi RAS Mapper -də istifadə edildikdə, ayaqla məsafə verir, ancaq ölçmənin metr olması lazım olduğunu göstərir. Məsələn, ölçü aləti ilə 10 fut ölçən bir məsafə əslində 10/0.30458 = 32.8 futdur. Bu model düzgün nəticələr verəcəkmi, yoxsa axının yönəldildiyi mənzərə əsasən 3.28 (1/0.3048) faktoru ilə azaldı?

Chris Goodell

Düşünürəm ki, bu sadəcə bir səhvdir və sən yaxşısan. Ancaq ikiqat yoxlamağa əmin olun.

Bu yazıda vahid haqqında danışdınız.

HECRAS ilə vahid haqqında oxşar sualım var.

Budur mənim vəziyyətim. Bu su hövzəsində (Souris çayı havzasında) ABŞ tərəfdaşları ilə işləyirik. HECRAS -ı inkişaf etdirmək və kalibr etməkdən məsuldurlar, bütün məlumatlarda və hesablamalarda ingilis dilindən, yəni ayaqlardan istifadə edirlər. Sonra bu modeli Kanadada istifadə edirik, lakin SI -ni axın, yüksəklik, məsafə və s. Baxımından istifadə edirik.

Sadəcə vahidi SI -ə çevirərək bu modeli istifadə edə bilərəm, yoxsa bu modeli istifadə etmək üçün düzgün bir yol deyil. Kanadada bu modeli istifadə etməyin düzgün yolu, giriş axınını İngilis vahidi olaraq hazırlamaq, modeli işə salmaq və seçdiyimiz vahidlərin nəticələrini oxumaqdır. Düşünürəm ki, bu fikir məntiqli və CivilGEO mütəxəssisləri ilə təsdiqləndi. Amma Chris hələ də fikrinizi bilmək istəyirəm.

Chris Goodell

Salam Qərb. Teorik olaraq, əsas RAS pəncərəsindəki Layihə Bölmələrini Dönüştür seçimindən istifadə edərək modelinizi başqa bir vahid sisteminə çevirə bilərsiniz. Ancaq bunu edərkən çox ehtiyatlı olardım, çünki RAS -ın bu dönüşüm proqramında bəzi səhvlər olduğu bilinirdi. Ən azından, dönüşümdən sonra hər şeyin düzgün edildiyini yoxlamaq üçün yoxlamalısınız. Lazım gələrsə proyeksiya faylınızı yeniləməyi unutmayın. RAS bunu avtomatik olaraq etməyəcək.


Məlumatların yenidən layihələndirilməsi və yenidən nümunələndirilməsi¶

Əvvəlki hissədə raster məlumatlarının necə saxlanılacağı göstərilmişdir. Bununla birlikdə, bir çox hallarda bu məlumatların yenidən işlənməsinə və yenidən nümunələnməsinə ehtiyac var. Praqmatik bir həll gdalwarp istifadə edər və bunu qabıqda edər. Bir tərəfdən, bu əlverişlidir, amma bəzən bu işi ara bir addım olaraq yerinə yetirməlisiniz və faylların yaradılması və silinməsi yorucu və səhvlərə meyllidir. İdeal olaraq, sizin üçün proyeksiyanı yerinə yetirəcək bir python funksiyanız olardı. GDAL bunu təyin edərək buna imkan verir yaddaşdakı raster faylları. Bunlar normal GDAL məlumat dəstləridir, lakin bunlar fayl sistemində, yalnız kompüterin yaddaşında yoxdur. Tez aralıq hesablamalar üçün əlverişli bir “scratchpad ”. GDAL eyni zamanda bir funksiyanı da təqdim edir. gdal.ProjectImage gdalwarp qabiliyyətlərinin çoxunu ortaya qoyur. Təklifi yerinə yetirmək üçün bu iki hiyləni birləşdirəcəyik. Nümunə olaraq, əvvəlki hissədə qeyd olunan İngilis Adaları üçün NDVI məlumatlarının İngiltərə üçün uyğun bir proyeksiya olan Ordnance Survey National Grid -də təkrarlanmasının lazım olduğu vəziyyətə baxacağıq.

Əsas komplikasiya ehtiyacdan irəli gəlir gdal.ProjectImage GDAL məlumat dəstləri üzərində işləmək. Əvvəlki bölmədə, NDVI məlumatlarını bir GeoTIFF faylına saxladıq, bu da bizə başlanğıc məlumatları verir. Hələ də çıxış verilənlər bazası yaratmalıyıq. Bu o deməkdir ki, proyeksiya edilməzdən əvvəl çıxış verilənlər bazasının geotransformasiyasını və ölçüsünü təyin etməliyik. Bu, orijinal verilənlər bazasının həcmi haqqında məlumat toplamaq, təyinat proyeksiyasına yansıtmaq və oradan piksel və geotransformasiya parametrlərinin sayını hesablamaqdır. Bu, bu vəzifəni yerinə yetirən (çox şərh olunan) bir funksiyadır:

Funksiya, edə biləcəyiniz bir GDAL yaddaşlı fayl obyekti qaytarır ReadAsArray və s. reproject_dataset diskə yazmır. Bununla birlikdə, məlumatların bir nüsxəsini çıxarmaqla yaddaşdakı rasteri GDAL tərəfindən dəstəklənən istənilən formatda saxlaya bilərik. Bu sözün həqiqi mənasında iki kod xətti alır.

Proqramın əsas hissəsini genişləndiririk ki, (i) yenidən redaktənin nəticəsini GeoTIFF faylı olaraq saxlayaq, (ii) əldə edilən məlumat faylını oxuyaq və (iii) planı hazırlayaq:


.QPJ Fayl Uzantısı

QPJ faylı, açıq mənbə, platformalararası GIS tətbiqi olan QGIS tərəfindən yaradılan bir GIS xəritəsi proyeksiya faylıdır. Yerin 3D təsvirini 2D müstəviyə çevirmək üçün istifadə olunan koordinat istinad sistemi (CRS) və xəritə proyeksiyasını ehtiva edir. QPJ faylları digər GIS sistemləri tərəfindən istifadə olunan PRJ xəritəsi proyeksiya fayllarına bənzəyir, lakin QPJ sənədlərində QGIS-ə məxsus məlumatlar var.

Müxtəlif həvəskar və peşəkar coğrafiyaşünarlar coğrafi məlumatları xəritələmək, redaktə etmək və dərc etmək üçün QGIS -dən istifadə edirlər. QGIS -in bir hissəsi olaraq coğrafiyaçılar coğrafi xüsusiyyətləri və sərhədləri ifadə etmək üçün xəritələnmiş nöqtələr, xətlər və çoxbucaqlardan istifadə edən şəkil faylları yaratmalıdırlar. QGIS şəkil faylları əslində biri QPJ faylı olan bir neçə komponentli fayllardan ibarətdir.

Bir QPJ faylı, QGIS -in Yerin əyri, ellipsoidal səthini 2D müstəviyə çevirmək üçün istifadə etdiyi bir koordinat istinad sistemini və xəritə proyeksiya funksiyasını təyin edir, belə ki, coğrafiyaçılar bu müstəviyə 2D fiqurlar (və beləliklə də xəritədə verilən məlumatları) çəkə bilərlər. Fərqli QGIS xəritələri, coğrafiyanın öyrəndiyi yerdən və coğrafiyanın yerinə yetirmək istədiyi təhlildən asılı olaraq Yer səthini 2D müstəviyə çevirmək üçün fərqli CRS və xəritə proqnozlarından istifadə edir.

QGIS xəritəsi proyeksiyalarına yalnız QGIS tərəfindən istifadə olunan QXIS-ə aid xüsusi parametrlər daxildir (AXIS, EXTENSION və AUTHORITY parametrləri kimi). Bu səbəbdən, bir coğrafiyaçı QGIS daxilində yeni bir şəkil faylı yaratdıqda, proqram iki xəritə proyeksiya faylı yaradır: QPJ faylı (QGIS tərəfindən istifadə olunur) və PRJ faylı (ötürmə zamanı ArcGIS kimi digər GIS proqramları tərəfindən istifadə edilə bilər). məlumatlar).


Tam müəyyən etmək üçün əlavə parametrlər tələb edən proqnozlar üçün, Xüsusiyyətlər cədvəli hər tətbiq olunan parametr üçün dəyərləri daxil etməyə imkan verir. Boz rəngdə göstərilən parametrlər düzəliş edilmir, çünki onlar adlandırılmış proyeksiya tərifinin bir hissəsidir.

Adlandırılmış proyeksiyanın parametrlərini dəyişdirmək üçün uyğun olan ümumi proyeksiya metodunu seçdi və bütün uyğun atributlar kopyalanacaq. Məsələn, UTM  zone tərifindən başlayaraq xüsusi bir proyeksiya yaratmaq üçün UTM  zone seçin, sonra UTM   -dan Transverse Mercator -a proyeksiya adını dəyişin. The


SAGA GIS

Gdal modulundan istifadə edərək Aster məlumatlarını (Geotiff olaraq) yüklədim
lavabo modulu, istifadə edildiyindən daha çox su tutma sahəsinin paralel modulu istifadə edildiyindən daha çoxdur
kanal şəbəkə modulu, ancaq kanal modulu heç bir şey göstərmir
boş ekran.

Geotiff layihəsi etməliyəm, yoxsa lazım deyil.

Güman edirəm ki, "Başlama Şəbəkəsi" və ızgara kimi su tutma sahəsi ızgarasından istifadə etmisiniz
ağlabatan nəticələr göstərir? "Başlatma Eşiği" parametrini düzəltmisinizmi?

Məlumatlarınızın proyeksiyası hansıdır?

Sahə sahəsini 1000000 həddi ilə başlanğıc ızgarası olaraq istifadə etdim. və
məlumatların proyeksiyası EPSG 4326 -da.

Ancaq mürəkkəb analizdən istifadə etsəm, nəticə verir, yəni kanal və hövzə
vektor faylı amma hövzəm kiçikdir və 200 kv. km ərazini tapmalıyam
o .

Başlanğıc həddini siçan ilə tutma sahəsinin üzərinə "tapdınız"
grid? Başlama həddinizdə bir problem olduğundan şübhələnirəm
qəbulu və ya toplama sahəsi ızgarası ilə.

BTW, mürəkkəb analiz modulu "Kanal Şəbəkəsi və Drenajdan istifadə edir
Hövzələr "modulu.

. və yalnız proyeksiyanıza baxdım, bu lat/lon - daha yaxşı bir təkrar
sahə hesablamalarını əldə etmək üçün düzbucaqlı bir koordinat sisteminə
sağ

Proj4 dialoq ızgarasından istifadə edərək ızgaranı yenidən yansıtdım. amma yoxlayanda
proqnozlaşdırılan şəbəkə təsviri, müəyyən edilməmiş bir koordinat sistemini göstərir
dəyişmiş hüceyrə dəyərini göstərir. yaxşıdı ya bir iş görmüşəm
çünki Geo tiff fayl proyeksiyam EPSG 4326 göstərirdi

Təsvir yeni təyin edilmiş koordinat sistemini göstərməlidir, yox
"müəyyən edilməmiş".

Enlem və Boylam koordinat sistemi ilə işiniz (Coğrafi WGS
84 EPSG: 4326), ölçü vahidiniz dərəcədir. 1 dərəcədir

111 km. Belə ki
bu koordinat sistemindəki pikseliniz 1 qövs saniyə olmalıdır (

30 m
proqnozlaşdırılan koordinat sistemi). Beləliklə, kanal çıxarmaq üçün eşik dəyəri
şəbəkə kiçik olmalıdır (1-0.001). Volkerin dediyi kimi, baxmalısınız
Catchment Area grid və izləmək istədiyiniz yerdən dəyərin nə olduğunu görün
kanal şəbəkəsi (Catchment Area vizualizasiyası üçün Rainbow istifadə edirəm
logaritmik miqyaslı və logaritmik uzanma faktoru ilə əvvəlcədən təyin edilmiş rəng ölçüsü
10000).

SAGA -nın proyeksiyanızı tanıdığını yoxlamalısınız (Coğrafi, görməlisiniz
Obyekt Xüsusiyyətləri altındakı Təsvir sekmesinde - sağ pəncərə). Yoxdursa, sən
Projeksiyon - Proj istifadə edərək proyeksiyanı qura bilərsiniz. 4 / Koordinat Referansını təyin edin
Sistem. Sonra, Koordinat Dönüşümünü (Şəbəkə) istifadə edərək yenidən layihələndirə bilərsiniz.
Bölgənizdən asılı olaraq yerli bir CRS seçə bilərsiniz və ya UTM -dən istifadə edə bilərsiniz
qlobal əhatə dairəsi.

İstinad olaraq, kanal şəbəkəsini müəyyən etmək istədiyiniz zaman bunu etməlisiniz
addımları izləyin:

doldurma lavaboları (Arazi Analizi - Ön İşləmə / Doldurma Lavaboları)
(Planchon/Darboux)), bu adı istifadə edərək DEM əldə edirsiniz
Giriş olaraq DEM aşağıdakı addımları yerinə yetirir

tutma sahəsini hesablayın (Arazi Analizi - Hidrologiya / Tutma Sahəsi
(Paralel)) kanal üçün ərəfənizi əldə etmək üçün Catchment Area ızgarasına baxın
mülayim bölgələr üçün başlanğıc 10 000 - 30 000 kvadrat aralığından istifadə edirəm
kanal izləmək üçün metr (11-33 piksel)

SAGA 2.0.8 istifadə edirsinizsə bu addım isteğe bağlıdır, ancaq istifadə edirsinizsə tələb olunur
2.0.7 və Terrain Analysis - Kanallar / D8 Flow Analizini işləməlisiniz
axın istiqaməti cədvəli almalısınız

kanalları izləmək və Terror Analizi - Kanallar /
Kanal Şəbəkəsi, DEM -ni Yüksəklik olaraq, Başlama Şəbəkəsi Catchment -dən istifadə edir
Sahə, Başlama Eşiği olaraq 10000-30000 (və ya dəyəriniz) olaraq təyin edə bilərsiniz
minimum seqment uzunluğu piksellərlə (kanallarınızda bir neçə piksel varsa
bu dəyərdən kiçik olsa izlənilməz).

Demək olar ki, yuxarıda təsvir olunan eyni problemlə qarşılaşmışam və edə bilmirəm
Tövsiyə olunanlara əməl edərkən boş məlumat dəstləri xəritələrində keçmiş kimi görünür
addımlar.

Məqsədim 30 metrlik ASTER DEM istifadə edərək "Havza hövzəsi təbəqəsi" yaratmaqdır
plitələr. Problemimin təsviri demək olar ki, yuxarıdakı Fella 41 ilə eynidir:

"Gdal modulundan istifadə edərək Aster məlumatlarını (Geotiff olaraq) yüklədim
lavabo modulu, istifadə edildiyindən daha çox su tutma sahəsinin paralel modulu istifadə edildiyindən daha çoxdur
kanal şəbəkə modulu, ancaq kanal modulu heç bir şey göstərmir
boş bir ekran ".

Etdiyimi təkrarlamaq üçün istifadəçidə göstərilən hər iki addımı izlədim
su hövzələri haqqında 2 bələdçisi və bu yazıda "Mihai" tərəfindən edilən təkliflər.

Arazi Analizi - Ön İşləmə / Doldurma Lavaboları (Planchon / Darboux) kimi
MiHai tərəfindən təklif edilmişdir).

Ardından Arazi Analizini - Hidrologiya / Tutma Sahəsini (Paralel) və hesablayıram
uyğun bir xəritə əldə edin.

Daha sonra vəziyyəti yaxşılaşdırmaq üçün səhifə 140 -da istifadəçi təlimatında göstərilən addımları izləyirəm
toplama sahəsi vizual:

Bir az təcrübədən sonra "Logaritmik" ə 100000 daxil olduğunu görürəm
Stretch Factor dəyər sahəsi daha yaxşı vizual şərh verir. İndi bacarıram
toplama sahələrinin yaradılmasının nəticələri haqqında daha yaxşı təəssürat yaratmaq üçün "

  1. Ancaq toplama sahələri yarandıqda boş Kanal Şəbəkəsi və
    Drenaj hövzələri oradan irəliləyə bilməz.

"Mürəkkəb analizdən istifadə etsəm, nəticə verir, yəni kanal və hövzə
vektor faylı "Kanal Şəbəkəsi və Drenaj Hövzələri ızgaraları, şəkil faylları və
digər ızgaraların hamısı gözəl yaradılmışdır. Ancaq yaratmağa gedəndə
Su hövzələri boş nəticələr alıram.

Volkerin "BTW, mürəkkəb analiz modulu" Kanalını istifadə edir
Şəbəkə və Drenaj Hövzələri "modulu." Kompleks analizdən istifadə edərkən bunları
məlumat dəstləri yaradılır, amma istifadəçi təlimatında və ya MiHai -də göstərilən addımları təqib etsəm
yuxarıda boş Kanal Şəbəkəsi və Drenaj Hövzələri alıram.

Su yaratmaq üçün mürəkkəb analizdən əldə edilən məlumat dəstlərindən istifadə etsəm
hövzələr - təəssüf ki, onlar da boşdur

Fella41 -dən kiçik bir fərq, Grid & gt tools & gt Merging -dən istifadə etməyimdir
İki ASTER plitəsini birləşdirmək üçün modul. Sonra bir kvadrat istifadə edərək DEM -i kəsirəm
maraq dairəsi üçün çoxbucaqlı. Nəticənin proyeksiyasını təyin edirəm
Projeksiyon = Proj istifadə edərək DEM birləşdirildi və kəsildi. 4 & gt Koordinat sistemini təyin edin
orijinal ASTER DEM plitələri ilə eynidir. Yəni

Proyektiv Coğrafi Koordinat Sistemi:

Fərqli eşikləri sınamaqla çox sınaq və səhv edərək nə olduğunu öyrənə bilmirəm
səhv gedir.

Həmçinin istinad edilən eşik dəyərlərinin funksiyasını başa düşmürəm
Volker, Fella41 və Mihai. Ancaq onları yüksək səviyyəyə qaldırmağımdan asılı olmayaraq
(1.000.000) və ya aşağı (0.1) Hələ də boş ekranlar ala bilmirəm.

Su Hövzələrinin tamamlanmasına kömək etmək üçün hər cür kömək çox böyük olardı
təqdir olunur.

Sahənizdəki ızgara yaxşı görünürsə, bu proyeksiya olmalıdır. Kimi
Coğrafi koordinatlarla işləyərkən hər şey çirkin ola bilər.
UTM və ya yerli CRS -ə yenidən müraciət edərdim.

Eşik aşağıdakı kimi işləyir:

  • ümumi tutma sahəsi, Arazi Analizi ilə əldə etdiyiniz şeydir - Hidrologiya / Tutma Sahəsi (Paralel)
  • total catchment area represent the area which generate flow to a certain pixel, and in your geographic data will have values given by the number of pixels multiplied by the area of an pixel
  • the threshold for channel network extraction is the value of the catchment area of a pixel from which water is assumed to flow in a channel
    -in your geographic data, although you obtain catchment area, the values are given by the following computation: in decimal degrees

In general, there is need to transform from geographic to metric coordinates,
because geographic data is used only for representing global and regional
data, and not for computing geomorphometric variables.

By changing the project it worked straight away.

For the benefit of others I did the following :

Using Projection - proj.4 > Coordinate Transformation (Grid), I converted
the merged ASTER DEM to WGS84 UTM 56S

Using Projection - proj.4 > Coordinate Transformation (Shapes), I also
converted my bounding area polygon shap to clip the DEM with to WGS84 UTM 56S

Followed the steps in the manual and your variation above to :

a) Terrain Analysis > Preprocessing > Sink removal (or Fill Sinks
(Planchon/Darboux))

b) Terrain Analysis > Hydrology > Catchment Area (Parallel)

c) Terrain Analysis > Channels > Channel Network

d) Terrain Analysis > Channels > Water Basins

a) Terrain Analysis - Compound Analyses > Basic Terrain Analysis and

b) Terrain Analysis > Channels > Water Basins

So thanks very much. This was very helpful.

The reference in my last post should say

"By changing the projection it worked straight away"

Now I can create these outstanding drainage datasets my next question is :

How can I improve the resolution of the DEM and thus the granularity and
smoothness of the channel network (ie not too many straight lines between
junctions)?

Using the 30 M ASTER DEM, (I assume best reprojected to UTM) and using the
contours module to enerate 1-5 metre contours can I then create a new higher
resolution DEM and rerun the drainage modules, or is there a better way to do
this.

If creating a new DEM from higher resolution contours :

a) is there a way to create smoother curvature contours that don't have
straight edges and pointy corners?

b) can anyone point me to any resources / tutorials to show how do create a
dem from contours ?

Well, the DEM can be improved in several ways.

I do not use ASTER GDEM only, because it has a lot of artefacts. Only where
SRTM is not available (at polar latitudes). If there is SRTM cover on your
area I would use an averaging of the two DEMs.

Also you can filter the DEM. You have the function Grid - Filter /
Multidirection Lee Filter. I use with great succes as estimated noise
(absolute and relative) 1 and the method noise variance given relative to mean
standard deviation.

There is another filter, which works very fine, especially when you have
lowlands (where smooth more) and highlands (where smooth less), but it must be
applied on maximum 2000x2000 ArcASCII grids (using cmd):

When you compute catchment area, the method used will influence the result.
You probably used Deterministic 8. You can try with Multiple Flow Direction
(MFD) or even Multiple Triangular Flow Direction.

I wouldn't use contours, because they not necessary increase the resolution
(numerically yes, but there will not be more features). I rather use Grid -
Splie Interpolation / Multilevel B-Spline Interpolation (from Grid) or a
simple Grid - Tools / Resampling with B-Splines.

a) Shapes - Grid / Contours Lines from Grid

b) you need to use an external application for smoothing lines (GRASS
v.generalize maybe)

c) Shapes - Points / Convert Lines to Points

d) Grid - Spline Interpolation / Multilevel B-Spline Interpolation (11 to 14
levels, 11 give smoother results)

Thanks for your quick response and helpful suggestions.

I tried the resampling module and changing cell size from 31 (which I assume
equates to 30 metre pixels from the ASTER DEM) down to 5 (assume 5 metres)
using BSpline interpretation.

I then created 5 metre contours which were very smooth and curved - perfect.
(However I wIll look also at your Grid - Spline Interpolation / Multilevel
B-Spline Interpolation (11 to 14 levels, 11 give smoother results)
suggestions)

I then went back and filled sinks, created the catchment area (Parallel),
created channel networks and compound analysis using default settings.

Tt this resolution the Wetness Index is quite beautiful and perfect also for
my needs.

For some reason Network Channels were created in great detail but Channels
were not created at all.

I have been fiddling around trying to generate the Chanels (ie the streams)
and trying to reduce the density of Network Channels so get a useful visual
both at full extent, (about 15,000 acres) and also in a small area (ie 150
acres).

I will try your various suggestions above and see how they work. I suspect it
has to do with the various thresholds which I am beginning to understand
better thank you.

Any further suggestions on guestimate thresholds for tidying up the channels
and channel network would be appreciated if they are easy to speculate on, but
in the meantime I will try som more trial and error to see if I can tighten
this up. I tried loading these datasets from the 30 m analysis but the main
channels dont line up well enough.

As I have experienced on all fronts to date, the SAGA GIS community assistance
has been most helpful.

Hopefully by documenting my experiences will help others who are in the same
boat as me in coming up to speed.

Thanks again for your assistance.

You can control the density of channels by setting also a minimum segment
length (in pixels). So if the threshold for extraction is 10000 acres, but the
length of the channel extracted by this threshold does not exceed 1000 pixels
(so with 5 m resolution will be 5000 m), then that channel will not be traced.

This way you can reduce parallel and very dense channels.

When you say that Network Channels were traced but not channels, do you refer
to the D8 Flow analysis function. That function trace all the channels, but
later when you trace them using the threshold value you will get the needed
Channel Network. You also need the Flow Direction grid as input (at least in
2.0.7).

That's making sense - will try experimenting with that.

Re Channels not showing I followed two approaches.

Firstly the steps to creating catchment areas then network channels.

Secondly doing compound analysis. If this multi step process uses the D8 Flow
analysis function it was done. If not then I didn't use it.

I will have a look at this and your other suggestions.

I now largely get it. Thank you again for your helpful comments and
təkliflər.

By using your tips in this post and Chapter 5 - Generating a Watershed Basins
Layer, in Volume 2 of Vern Cimmery's User Guide for Saga (version 2.05) I have
been able to generate some extremely useful drainage datasets.

For the benefit of others -

Starting with importing the ASTER 30m DEM using File >GDAL > import raster
(being the best resolution I can find for this area), Merging the DEMs and
resampled the grid from 30 m to 5 m, I can generate excellent realistic
looking and smooth curvature contours. It also sets me up for more realistic
looking channels at small resolution (ie 100 acres versus 10,000 acres)

I ran several different Channel Networks scenarios and saved them each time
with a description of the thresholds and segment lengths. I adjusted the
initation grid, initiation threshold to include 0, 1000, 100,000 and
1,1000,000. For various of these outputs I also adjusted minimal segment
length in pixels from 10 to 500.

Running these scenarios gave me a good understanding of how Network Channels
yaranır. By zooming in so I could see and count the number of pixels
that make up a channel network segment I was able to adjust the minimal
segment length that worked best for producing tidy datasets (ie to reduce the
number / clutter of parallel channels within a few picels of each other) at
each of the different initiation thresholds. I also discovered the ruler at
the bottom of the map represents the number of pixels. This was helpful as a
guide to adjusting minimal segment length rather than counting individual
number of pixels that made up a segment.

I was then able to use the different Channel Networks thresholds (ie 0, 1000,
100,000, 1,000,000) to create different size waterbasins, also adjusting for
the minimum number of pixels which make up a basin which was only necessary
the small initiation thresholds (ie 0, 1000).

After all this I then ran Terrain Analysis > Compound Analyses > Basic Terrain
Analysis module to generate a multitude of other useful datasets, in
particular wetness index.

I rediscovered the useful "Channels" layer that I reported before was not
being generated was in fact generated by the Terrain Analysis > Compound
Analyses > Basic Terrain Analysis module

I havent yet worked out how to generate the D8 Flow Analysis function.

Of the datasets generated from Basic Terrain Analysis I dont yet understand
the role, function or application of the the following layers :

Anyway this has been a very useful exerise not only to generat the datasets I
was trying to generate but also other datasets I didnt know were possible and
very importantly has significantly increased my cognition and understanding of
how to us Saga GIS tools.

So thanks very much again MiHai for all of your help, Volker for your helpful
comments in this and other posts, Vern Cimmery for the manual and also Kim
Cimmery for all your help in bringing me up to speed with Saga GIS generally
and in relation to pattern classification to identify Tree Crowns.


#Styling Polygons

As I mentioned earlier, a benefit of SVG is that we can style with CSS: we can color constituent countries by style rules that assign the fill property. However, we first need to give each country its own path element, rather than sharing one. Without distinct path elements, we have no way to assign distinct colors.

Within the uk.json TopoJSON file, the Admin-0 map subunits are represented as a feature collection. By pulling out the features array, we can compute a data join and create a path element for each feature:

We can also compute the "class" attribute as a function of data, so that in addition to the subunit class, each country has a three-letter class corresponding to its ISO-3166 alpha-3 country code. These allow us to apply separate fill styles for each country:

The styles above also hide Ireland completely, but we’ll bring it back in the next step when we display boundary lines. The map so far:


3 Cavab 3

From the content of the error message, it would appear that you're trying to run the algorithm on a field which contains null values, so that's the main problem.

Açın Field Calculator and check the box to Mövcud sahəni yeniləyin and select the the field you're creating your points from, which I'm assuming is the "total" field in your data.

Enter the following expression:

This will insert 0's in place of all the nulls, which will allow you to run the Random Points tool now.

I've done this, but I still get the same error message. I wonder if my data is not setup to run a dot map? I don't know how to articulate things as well as I should - I'm still relatively new to this.

In Excel, re-arrange the table so that you have only one row for each Geography value. Create a column for each language, and put the number of speakers of that language in that column. Eg, convert this

I did this manually, but I think you can do it automatically using something called a pivot table, as mentioned here and here.

Save the new table as a CSV. Add the new table to your QGIS project.

Join the new table to the neighborhood layer. Now you have a one-to-one relationship between table rows and neighborhood geometries.

Run Random points inside polygons once for each language category. Be sure to name each new layer with the name of the language.

There must be something wrong with how I'm joining the tables together. I recreated my table like you instructed. I've loaded both the shape file and the .CSV in QGIS. When I go to join, I select the shape file and join them. Then when I open the Attribute Table for my data set, the numbers appear correct. Nothing is Null. But when I run the Random Points tool it gives me the same error.

There are so many different screenshots at the link you posted I'm having a hard time finding the error message. Can you please add the text of the error message to your original question? Use the edit link under the question.

The error message is now in the original question.

The input data you have has blank cells as one of the values. QGIS uses the field you specified to add points up to the number in the field. To do this it converts the field into an integer, but it cannot convert a blank cell into an integer.

The easiest way for you to fix this is to put 0's into the blank cells of your source data.


FWTools: Open Source GIS Binary Kit for Windows and Linux

The FWTools kits also aims to track the latest development versions of the packages included as opposed to official releases. While this may mean the packages are less stable, it is intended to give folks a chance to use the latest and greatest. FWTools releases also are a means by which I make recent development version bug fixes available to a wider audience than would be prepared to build them from the source.

With FWTools releases, I also endeavor to build in as many optional components as possible. Thus, I include support for ECW, JPEG2000, HDF and other file formats that require extra libraries.

Linux FWTools releases are intended to be distribution and packaging system agnostic. They should install on pretty much any x86 style Linux system released within the last few years.


Videoya baxın: GDALOGR - Beginners Tutorial