Daha çox

Xəritə Cəbri/Məkan Analitikində dəyişən adları olan rasterlərdə hesab?

Xəritə Cəbri/Məkan Analitikində dəyişən adları olan rasterlərdə hesab?


Dəyişən adları olan müvəqqəti rasterlərdə hesablamaların necə aparılacağını anlamağa çalışan Python təcrübəsiziyəm.

Minimum iş nümunəsi sınağım aşağıda verilmişdir, amma nəticədə bunu minlərlə daha böyük ızgaralarla etmək istərdim, buna görə də səmərəliliyi artırmaq üçün hər şeyi sona qədər müvəqqəti saxlamalı olacağımı gözləyirəm.

Aşağıdakı kod son sətirdə uğursuz olur. Səhv mesajı belədir:

RuntimeError: ERROR 000732: Giriş Raster: Dataset bin2 yoxdur və ya dəstəklənmir.

Son sətri daimi fayl adları ilə əvəz etsəm:

tmpout = Eksi ("D:/tmp/tmp" + 'il', "D:/tmp/tmp" + 'il boyu')

sonra yaxşı işləyir. Win32 -də Python 2.6.5 (r265: 79096, 19 Mar 2010, 21:48:26) ilə ArcGIS 10.0 -da IDLE çalışıram [MSC v.1500 32 bit (Intel)] win32.

arcpy.sa -dan idxal arcpy -dən idxal arcpy idxal * arcpy.env.extent = arcpy.Extent (0,0,2,2) arcpy.CreateRandomRaster_management ("D:/tmp/", "tmp1", "" , "", 1) arcpy.CreateRandomRaster_management ("D:/tmp/", "tmp2", "", "", 1) bin1 = arcpy.Raster ("D:/tmp/tmp1") bin2 = arcpy.Raster ("D:/tmp/tmp2") aralığında il üçün (2, 3, 1): illik = il - 1 binyear = "bin" + 'il' binyearold = "bin" + 'illik' çap "sonrakı ildə uğursuz olur xətt "tmpout = Minus (binyear, binyearold)

Mən düşün problem bundadırbin1bin2müvəqqəti rasters üçün göstərici fərqli olaraq simlər kimi şərh olunur, amma niyyətimi necə başa salacağımı anlaya bilmirəm. Əvəz edilməsiiki illikiləRaster (iki illik)məsələn, kömək etmədi.


ilk döngənizdə tmpX adlı rasterlər yaradırsınız. İkinci sətir yararsızdır və üçüncü sətir hər döngədə yazdığınız bir raster obyektdir. BinX əvəzinə ikinci döngənizdə tmpX adını istifadə etməlisiniz.


Yəqin ki, qarışıqam, amma 'd:/tmp' ə raster yazmaq hələ də diskə yazır? bu səbəbdən 'd:/tmp' döngəsində təyin olunmur və proses skriptlə eyni qovluqda raster tapmağa çalışdıqda bir səhvlə qarşılaşılır.

Burada qəbul edilən cavabda qeyd edildiyi kimi bir yol göstərmək əvəzinə 'in_memory/bin1' və 'in_memory/bin2' cəhd edə bilərsiniz. görünür, bu hələ də standart gdb -ə yaza bilər, amma prosesi sürətləndirə bilir.


Xəritə cəbr

Xəritə cəbrKartoqrafik modelləşdirmə olaraq da bilinən, raster məlumatlarından istifadə edərək coğrafi analiz aparmaq üçün müxtəlif vasitələrdən, operatorlardan və funksiyalardan istifadə edən sadə, lakin güclü bir cəbrdir. [1] İstənilən məlumatı və ya vahidləri çevirmək, dəyərləri birləşdirmək və ya müəyyən aralığdakı məlumatları seçmək kimi nəticələr əldə etmək üçün riyazi bir funksiyadan istifadə edərək coğrafi məlumatları idarə etmək üçün istifadə edilə bilər.

Xəritə cəbrini Dr. C. Dana Tomlin 1980 -ci illərin əvvəllərində. Tomlin coğrafi məlumatların işlənməsi üçün ibtidai operatorlar təklif etdi. Bu üsul bütün GIS əməliyyatlarını dörd sinifə bölməklə yalnız raster giriş və çıxışlarını emal edə bilər: yerli, fokus, qlobal və zonal. Bir modeli ifadə edərkən, bu model, o modelin təqib etməsi üçün sadə bir skript yaratmaqda faydalıdır. [2] Bu "skript" mürəkkəb vəzifələri yerinə yetirə biləcək bir prosedurun yerinə yetirilməsinə imkan verir.


Raster Kalkulyator alətindən istifadə

Alət informasiya qutusunda Xəritə Cəbri ifadəsi yaratmaq üçün istifadə olunan dörd əsas sahə var:

Raster Kalkulyator vasitəsi informasiya qutusu nümunəsi

Qatlar və dəyişənlər

Giriş qatları və dəyişənlər siyahısı Xəritə Cəbri ifadəsində istifadə edilə bilən girişi müəyyənləşdirir. Aləti tək başına istifadə edirsinizsə, siyahı məzmun cədvəlindəki qatları ehtiva edir. Alət ModelBuilder -də istifadə edildikdə, siyahıda məzmun cədvəlindəki təbəqələr, modeldə yaradılan nəticələr və digər model dəyişənləri olacaq.

Kalkulyator və operator düymələri

Kalkulyator düymələri ifadəyə rəqəmsal dəyərlər daxil etməyə imkan verir. Operator düymələri ifadəyə riyazi (əlavə, bölmə və sair) və məntiqi (böyük, bərabər və sair) operatorları daxil etməyə imkan verir. Bu düymələrdən hər hansı birinə basmaqla həmin nömrə və ya operator göstəricinin hazırda yerləşdiyi ifadəyə daxil ediləcək. İfadənin düzgün Xəritə Cəbri sintaksis qaydaları qaydalarına uyğun olmasını təmin etmək üçün diqqətli olun.

Alətlər

Alətlər siyahısı, Xəritə Cəbri ifadəsində istifadə edilə bilən alətlərin rahat bir seçimidir. Siyahıdakı bir aləti tıkladığınızda, alətin adı və mötərizənin [və] açılması işarənin hazırda yerləşdiyi ifadəyə yerləşdiriləcək. Alət tərəfindən tələb olunan qalan giriş daha sonra daxil edilməlidir. Bir alət ifadənin hər hansı bir yerinə yerləşdirilə bilər, ancaq etibarlı Map Cəbr sintaksisini yaradan bir mövqedə yerləşdirilməlidir.

İfadə

İfadə, icra ediləcək Map Cəbr ifadəsidir. İfadə etibarlı sintaksislə daxil edilməlidir. Xəritə Cəbrinin sintaksis qaydaları haqqında daha çox məlumat əldə etmək üçün aşağıdakı bölməyə baxın.


Sintaksis

Bu Boolean əməliyyatında istifadə ediləcək ilk giriş.

Digər parametr üçün raster göstərildiyi təqdirdə, bu parametrə giriş olaraq bir ədəd istifadə edilə bilər. Hər iki giriş üçün bir nömrə təyin edə bilmək üçün əvvəlcə hüceyrə ölçüsü və ölçüsü mühitdə təyin edilməlidir.

Bu Boolean əməliyyatında istifadə ediləcək ikinci giriş.

Digər parametr üçün raster göstərildiyi təqdirdə, bu parametrə giriş olaraq bir ədəd istifadə edilə bilər. Hər iki giriş üçün bir nömrə təyin edə bilmək üçün əvvəlcə hüceyrə ölçüsü və ölçüsü mühitdə təyin edilməlidir.

Gəlir Dəyəri

Çıxış dəyərləri 0 və ya 1 olacaq.


Təsvir, Avtomatlaşdırma və Tətbiqlər

İxtisasın son kursuna xoş gəldiniz (əlbəttə ki, capstone layihəsinə davam etməsəniz!). ArcGIS haqqında öyrəndiyiniz məlumatlardan istifadə edərək, raster hesablamalar və uyğunluq təhlili kimi texniki tapşırıqları yerinə yetirin. Bu sinifdə dörd həftə davam edən modullar ərzində GIS daxilində məkan təhlili və tətbiqləri ilə özünüzü rahat hiss edəcəksiniz: 1-ci həftə: Siz uzaqdan algılanan və peyk görüntüləri haqqında hər şeyi öyrənəcəksiniz və elektromaqnit spektri ilə tanış olacaqsınız. Bu həftənin sonunda, siz peyk görüntülərini onlayn olaraq tapa və yükləyə və iki ümumi analiz növü üçün istifadə edə bilərsiniz: NDVI və təlim keçmiş təsnifat. 2 -ci həftə: Parametrlərdən, ön şərtlərdən, dəyişənlərdən və yeni alətlər dəstindən istifadə edən böyük emal iş axınları yaratmaq üçün ModelBuilder -dən necə istifadə edəcəyinizi öyrənəcəksiniz. Həqiqətən ətraflı müzakirə etməyə vaxtımız olmayan bir neçə mövzunu da araşdıracağıq, ancaq gələcəkdə digər yollarla öyrənmək üçün iştahınızı aça bilər: coğrafi kodlaşdırma, zamana bağlı məlumatlar, məkan statistikası və ArcGIS Pro. 3 -cü həftə: Üçüncü həftədə, kəsmə doldurma vasitəsi, təpələr, mənzillər və daha çox kimi yeni, xüsusi vasitələrdən istifadə edərək rəqəmsal yüksəklik modelləri hazırlayırıq və istifadə edirik. Biz də çox əhəmiyyətli bir alqoritmdən keçirik: uyğunluq təhlili. 4 -cü həftə: Biz son həftəyə hələ ixtisaslaşmadığımız bir neçə məkan analitik vasitəsi haqqında danışmaqla başlayacağıq: Region Qrupu öz zonalarımızı quracaq, Təpə örtüyünü hamarlaşdırmaq üçün Odak Statistikası, Dəyərləri dəyişdirmək üçün Yenidən Sınıflandırın və Sıxlıq səthi yaratmaq üçün nöqtə sıxlığı. Nəhayət, GIS mövzularını təkbaşına öyrənməyə başladıqca daha çox araşdırmaq istəyə biləcəyiniz daha bir neçə şey haqqında danışaraq yekunlaşdıracağıq. Müstəqil bir kurs olaraq və ya Coğrafi İnformasiya Sistemləri (GIS) İxtisaslarının bir hissəsi olaraq Coğrafi və Ətraf Mühit Analizini götürün. Bu kursu almadan əvvəl bu ixtisas üzrə birinci, ikinci və üçüncü kursları və & quot; CBS əsasları & quot & quot; GIS Məlumat Formatları, Dizayn və Keyfiyyət & quot; və & quot; Kosmik və Ətraf Mühit Analizi & quot; Dördüncü sinfi bitirərək, İxtisas kapitonunda uğur qazanmaq üçün lazım olan bacarıqları əldə edəcəksiniz.

Получаемые навыки

Coğrafi Məlumat Sistemi (CİS), Görüntü Təhlili, Məkan Analizi, peyk görüntüləri, Gis Tətbiqləri

Рецензии

Çox yaxşı bir kurs, ancaq bu kursda qiymətləndirilən bəzi mövzular, xüsusən uzaqdan zondlama şəkillərinin işlənməsi (GIS -ə endirilməsi və idarə edilməsi) yaxşı izah edilməmişdir.

İnanılmaz bir kurs! yaxşı təşkil edilmiş, çox məlumatlı və qaynaqlar və faydalı materiallar ilə zəngin, kurs mövzularını müzakirə etmək üçün əla bir müzakirə forumu ilə.

Bu modulda, peyk və hava görüntülərinin çəkilməsi və məlumat məhsulları haqqında öyrənməzdən əvvəl, uzaqdan algılanan məlumatlara və elektromaqnit spektrinə girişdən başlayaraq, uzaqdan zondlama və peyk görüntüləri haqqında hər şeyi öyrənəcəyik. İnternetdə peyk görüntülərini necə tapmağı və yükləməyi və iki fərqli ümumi analiz növündə necə istifadə edəcəyinizi öyrənəcəksiniz: NDVI və təlim keçmiş bir təsnifat. İkinci dərsdə, Raster Kalkulyatoru və Məkan Analitikindən istifadə edərək görüntü analizini dəstəkləmək üçün bəzi əsas vasitələrdən istifadə etməyi öyrənəcəksiniz.


Məkan Analitikinə ümumi baxış

ArcGIS Spatial Analyst, hüceyrə əsaslı (raster) analiz aparmaq üçün zəngin bir alət dəsti təqdim edir. GIS məlumatlarının üç əsas növündən (raster, vektor və VÖEN), raster məlumat quruluşu məkan analizi üçün ən əhatəli modelləşdirmə mühitini təmin edir.

Hüceyrə əsaslı sistemlər, dünyanı bir ızgara quruluşuna əsaslanaraq hüceyrələr adlanan vahid vahidlərə bölür. Hər bir hüceyrə yerin müəyyən bir hissəsini təmsil edir, məsələn, bir kvadrat kilometr, hektar və ya kvadrat metr. Hər bir hüceyrəyə yerləşdiyi xüsusiyyətə və ya xüsusiyyətə uyğun bir dəyər verilir və ya yüksəklik, torpaq növü və ya yaşayış təsnifatı kimi yeri təsvir edir. Hüceyrə əsaslı bir sistemdə coğrafi mövqe bir xüsusiyyət olaraq təyin edilmir, lakin yerləşmə perspektivi olaraq bilinən saxlama quruluşuna xasdır.

Yerləşdirmə perspektivi ArcGIS Spatial Analyst -ə davamlı məlumatları (məsələn, yüksəklik, yağ konsentrasiyası və səs) daha təsirli şəkildə saxlamağa imkan verir. Davamlı məlumatlarda, hər bir yerə təyin edilmiş bir miqdar, böyüklük və ya intensivlik verilir və dəyərlər bir -birinə nisbətən mənalıdır. Yerləşdirmə perspektivi həm diskret (məsələn, torpaqdan istifadə və bitki örtüyü növü), həm də hər bir alət dəstini müşayiət edən geniş müzakirələrdə aydın olacaq məkan analizində daha çox müxtəlifliyə imkan verir.

ArcGIS Spatial Analyst ilə analiz edərkən problemləri həll etmək, raster məlumatlarını və narahatlıqlarını anlamaq üçün daha çox məlumat əldə etmək üçün aşağıdakılara baxın:


Skript sintaksisi

NetCDFConvert2DVariableInNetCDFsInArcGISTableToArcGISRasters_GeoEco (masa, inputFileField, outputRasterField, variableName, xLowerLeftCorner, yLowerLeftCorner, cellSize, nodataValue, qoymalarını, güzgü, flip, swapHemispheres, coordinateSystem, projectedCoordinateSystem, geographicTransformation, resamplingTechnique, projectedCellSize, registrationPoint, clippingRectangle, mapAlgebraExpression, buildPyramids, harada, orderby, istiqamətləri, skipExisting, basePath)

NetCDF fayllarının giriş yollarını ehtiva edən sahə.

Bu sinif hazırlandığı anda netCDF-nin iki əsas versiyası inkişaf mərhələsində idi: 3.6.2 versiyasında olan NetCDF-3 və 4.0 beta 1 versiyasında olan NetCDF-4, bu sinifin hazırkı tətbiqi işlənə bilər. NetCDF-3 faylları, lakin NetCDF-4 faylları deyil. NetCDF-4 dəstəyi yaxın gələcəkdə əlavə olunacaq.

Sıxılmış faylları dəstəklənən sıxılma formatında təqdim etsəniz, onlar avtomatik olaraq açılacaq. Fayllar arxiv formatında (məsələn .zip və ya .tar) sıxışdırılıbsa, hər bir arxivdə bir alt qovluq olmamalı olan bir fayl olmalıdır.

Yaratmaq üçün çıxış rasters olan sahə.

NetCDF faylında bir dəyişənin adı. Dəyişən iki ölçü və tam ədəd və ya üzən nöqtə məlumat tipinə malik olmalıdır.

Dəyişənin adını bilmirsinizsə, ən yaxşı təxmininizi verin. Bu adla bir dəyişən yoxdursa, bir ValueError qaldırılacaq. Səhv mesajı, giriş sənədində mövcud olan bütün 2D dəyişənlərin adlarını sadalayacaq. Bu adlardan birini seçin və yenidən cəhd edin.

NetCDF başlığını bir mətn sənədinə atmaq üçün NetCDF Başlığını Ayarla alətindən də istifadə edə bilərsiniz. Dəyişən adları "dəyişənlər:" ilə başlayan bölmənin faylını axtararaq tapa bilərsiniz. Hər bir dəyişən aşağıdakı kimi sıralanır:

harada məlumat növü dəyişənin məlumat növüdür (məsələn, "float"), ad dəyişənin adıdır və ölçülər dəyişənin ölçülərinin siyahısıdır. Məsələn, netCDF Aviso DT-MADT geostrofik cərəyanlarından alınan aşağıda, Grid_0001 və Grid_0002 iki ölçülü iki dəyişən var:

Bu vasitə netCDF dəyişənini bir ArcInfo ASCII Şəbəkəsinə çıxarır və ArcGIS ASCII -dən Raster coğrafi işləmə alətindən istifadə edərək rasterə çevirir. Bu alət və ArcGIS raster formatında uğurla ArcGIS raster formatına çevrilə bilən netCDF dəyişənlərinin növlərini məhdudlaşdıran bir neçə məhdudiyyət var.

Float və ya ikiqat məlumat növünü istifadə edən NetCDF dəyişənlərində "sonsuzluq" (INF) və ya "rəqəm deyil" (NAN) dəyərləri olmamalıdır. Bu dəyərlər aşkar edilərsə, ValueError qaldırılacaq.

ArcGIS raster formatı 32 bitlik float məlumat növünü dəstəkləyir, lakin 64 bitlik ikili məlumat tipini dəstəkləmir. İkiqat məlumat tipli bir netCDF dəyişənini təmin edərsəniz, ArcGIS ASCII-dən Raster coğrafi işləmə alətindən istifadə edərək 32 bitlik float rasterə çevriləcək. Bu vəziyyətdə alətin davranışı sənədləşdirilmir. ArcGIS 9.1 -də belə görünür:

Eksponentin -38 ilə +38 arasında dəyişdiyi dəyərlər, 32 bitlik float məlumat tipinin daha kiçik mantissa səbəbiylə müəyyən bir dəqiqlik itirilsə də, ortaya çıxan 32 bitlik float rasterində düzgün şəkildə təmsil olunur.

Eksponentin -38 -dən az olduğu dəyərlər (məsələn -39, -40 və s.) 0 -a çevrilir.

Eksponentin +38 -dən böyük olduğu dəyərlər, dəyərin işarəsindən asılı olaraq -INF və ya +INF -ə çevrilir (məsələn -5.3083635279597874e -212, ArcCatalog GUI -də -1.#INF, 2.5502286890301497e +084 görünür) kimi 1.#INF).

ArcGIS 9.1 ASCII to Raster aləti, tam ədəd rasterlərini çevirərkən bəzi qəribəliklər də nümayiş etdirir:

-128 NODATA dəyəri təyin edilmədiyi təqdirdə, 8 bitlik tam məlumatlar üçün, alət -128 dəyəri görünsə, alət 16 bitlik raster yaradacaq. 0 kimi fərqli bir NODATA dəyərinin göstərilməsi, -128 görünsə yenə də 16 bitlik bir raster verir.

Eynilə, 16 bitlik tam məlumatlar üçün, NODATA dəyəri təyin edilmədiyi təqdirdə, -32768 dəyəri ASCII faylında görünsə, alət 32 bitlik raster yaradacaq.

Daha pisi, 32 bitlik tam ədədlər üçün, NODATA dəyəri təyin edilmədiyi təqdirdə, -2147483648 dəyəri görünsə, alət bir səhv barədə məlumat verəcəkdir. Qəribə olsa belə, -2147483647 dəyəri nə olursa olsun həmişə NODATA -ya tərcümə olunur.

Bütün növ məlumatlar üçün, məlumat növü üçün mümkün olan ən kiçik dəyər olmayan bir NODATA dəyəri təyin etdiyiniz zaman alət qəribə davranışlar yaradır. Məsələn, 0-dan 255-ə qədər məlumat aralığı və 0-a NODATA dəyəri təyin edilərsə, alət imzasız 8-bitli çıxış rasterini istehsal edir. Ancaq 1 -ə NODATA dəyəri təyin olunarsa, imzalanmış 16 bitlik bir çıxış rasteri istehsal edir və ArcCatalog, Raster Dataset Xüsusiyyətləri altında NoData Dəyərinin -32768 olduğunu göstərir, lakin Identify aləti 1 dəyərinə malik hüceyrələri əslində NODATA olduğunu göstərir. Bənzər qəribə nəticələr, mümkün olan ən kiçik dəyər olmayan bir NODATA təyin edərkən digər məlumat növlərinin tam ədəd rasterləri üçün əldə edilə bilər.

Rasterin sol alt küncünün X koordinatı.

Koordinat, bu hüceyrənin ortası deyil, sol altdakı hüceyrənin küncü üçündür. Məsələn, raster bütün Yerin coğrafi proyeksiyasıdırsa, sol alt küncün koordinatı 180 dərəcə Qərb uzunluğuna uyğun olaraq -180.0 olardı.

Rasterin sol alt küncünün Y koordinatı.

Koordinat, bu hüceyrənin ortası deyil, sol altdakı hüceyrənin küncü üçündür. Məsələn, raster bütün Yerin coğrafi proyeksiyasıdırsa, sol alt küncün koordinatı 90 dərəcə Cənub eninə uyğun olaraq -90.0 olardı.

Məsələn, raster 720 sütun və 360 sətirdən ibarət olan bütün Yerin coğrafi proyeksiyasıdırsa, coğrafi dərəcənin 1/2 hissəsinə uyğun olan 0,5 ölçüsündə bir hüceyrəyə sahib olardı.

Əsas məlumat formatı hüceyrələrin kvadrat olmasını tələb edir. Hər bir ölçü üçün hüceyrə ölçüsünü təyin etmək mümkün deyil.

Bir hüceyrənin məlumat olmadığını göstərən dəyər.

Doğru olarsa, şəkil çevrilmədən əvvəl köçürülür (diaqonal ox ətrafında çevrilir). Şərq/qərb oxu sola və sağa yox, yuxarı və aşağıya doğru gedən bir görüntünü düzəltmək üçün bu seçimdən istifadə edin.

Doğru olarsa, şəkil çevrilmədən əvvəl şaquli ox ətrafında çevriləcək. Olmalı olduğu şeyin "güzgü görüntüsü" olan bir görüntünü düzəltmək üçün bu seçimdən istifadə edin.

Doğru olarsa, şəkil çevrilmədən əvvəl üfüqi ox ətrafında çevriləcək. Ters çevrilmiş bir görüntünü düzəltmək üçün bu seçimdən istifadə edin.

Doğru olarsa, görüntünün şərq və qərb yarımkürələri dəyişdiriləcək. Sakit okeanda yerləşən 0 -dan 360 -a qədər olan qlobal görüntünün istiqamətini Atlantik okeanında -180 -dən +180 -ə və ya əksinə vizaya dəyişdirmək üçün bu seçimdən istifadə edin.

Çıxış rasterini təyin etmək üçün koordinat sistemi. Bir dəyər verilmirsə, çıxış rasterinin koordinat sistemi təyin olunmamış qalacaq.

Çıxış rasterini layihələndirmək üçün yeni koordinat sistemi.

Raster yalnız orijinal proyeksiya təyin edildikdə yeni bir koordinat sisteminə proqnozlaşdırıla bilər. Orijinal koordinat sistemini təyin etmədən yeni bir koordinat sistemi göstərsəniz bir səhv ortaya çıxacaq.

ArcGIS Project Raster vasitəsi proyeksiyanı yerinə yetirmək üçün istifadə olunur. Bu alət üçün sənədlər, yeni koordinat sistemi üçün bir hüceyrə ölçüsünü də göstərməyinizi tövsiyə edir.

Müəyyən koordinat sistemləri üçün ArcGIS 9.2 Layihə Raster alətinin çox kiçik olan ixtiyari ölçüdə proqnozlaşdırılan rasteri kəsdiyini görünür. Məsələn, qlobal bir MODIS Aqua 4 km xlorofil görüntüsünü coğrafi koordinatlarda Lambert_Azimuthal_Equal_Area'ya mərkəzi meridianı -60 və mənşəyi eni -63 olan, planetin yalnız dörddə birini göstərmək üçün kəsildikdə. Bu problem, Project Raster ArcGIS istifadəçi interfeysindən interaktiv olaraq çağırıldıqda baş vermir, yalnız alət proqramla işə salındıqda baş verir (geoprosessorun ProjectRaster_management metodu). Project Raster -dən özünüz istifadə etdiyiniz zaman bunu görə bilməzsiniz, ancaq geoprosessing əməliyyatlarının bir hissəsi olaraq Project Raster -ə müraciət edən MGET vasitələrindən istifadə etdiyiniz zaman baş verə bilər.

Bu problemlə qarşılaşsanız, bu şəkildə işləyə bilərsiniz:

Birincisi, orijinal koordinat sistemində çıxış rasterini əldə etmək üçün yeni bir koordinat sistemi göstərmədən bu vasitəni işlədin.

ArcCatalogda, rasteri yeni koordinat sisteminə yansıtmaq üçün Project Raster alətindən istifadə edin. Bütün rasterin mövcud olduğunu, çox kiçik bir ölçüdə kəsilmədiyini yoxlayın.

ArcCatalogda, kataloq ağacında sağ tıklayarak, Xüsusiyyətləri seçərək və Aşağıya doğru genişlənərək proqnozlaşdırılan rasterin miqyasına baxın.

İndi, rasteri proyeksiya edən MGET alətini işə salmadan əvvəl, Ətraf mühit ayarını axtardığınız dəyərlərə qoyun. MGET alətini ArcCatalog və ya ArcMap-dən interaktiv olaraq çağırırsınızsa, alətin informasiya qutusundakı Ortamlar düyməsini basın, Ümumi Ayarları açın, Genişlik açılır menyusunu "Aşağıda göstərildiyi kimi" olaraq dəyişdirin və axtardığınız dəyərləri daxil edin. Bir coğrafi işləmə modelindən istifadə edirsinizsə, modeldəki vasitəyə sağ vurun, Dəyişən et, Ətraf mühitdən, Ümumi Ayarlar, Genişlik seçin. Bu, Extent modelinizdə MGET alətinə əlavə olaraq dəyişən olaraq yerləşdiriləcək. Genişlik dəyişənini açın, "Aşağıda göstərildiyi kimi" olaraq dəyişdirin və axtardığınız dəyərləri yazın. MGET alətini proqramlı olaraq çağırırsınızsa, geoprosessorun Genişlik xüsusiyyətini axtardığınız dəyərlərə təyin etməlisiniz. Bu və ümumiyyətlə ətraf mühit parametrləri haqqında daha çox məlumat üçün ArcGIS sənədlərinə baxın.

MGET alətini işə salın. Çıxış rasterinin ölçüsü indi uyğun ölçüdə olmalıdır.

Orijinal koordinat sistemi ilə yeni koordinat sistemi arasında çevrilmək üçün istifadə edilən çevrilmə üsulu.

Bu parametr ArcGIS 9.2 tərəfindən təqdim edilən yeni bir seçimdir. Bu parametri istifadə etmək üçün ArcGIS 9.2 olmalıdır.

Bu parametr yalnız rasterin yeni bir koordinat sisteminə yansıtılması lazım olduğunu və yeni sistemin orijinal koordinat sistemindən fərqli bir datum istifadə etdiyini və ya transformasiya tələb edən iki koordinat sistemi arasında başqa bir fərqin olduğunu göstərdiyiniz zaman lazımdır. Bir çevrilməyə ehtiyac olub olmadığını müəyyən etmək üçün aşağıdakı proseduru məsləhət görürəm:

Birincisi, orijinal koordinat sistemində çıxış rasterini əldə etmək üçün yeni bir koordinat sistemi göstərmədən bu vasitəni işlədin.

Ardından, istədiyiniz koordinat sisteminə yansıtmaq üçün çıxış rasterində ArcGIS 9.2 Project Raster alətindən istifadə edin. Bir coğrafi çevrilməyə ehtiyacınız varsa, bu vasitə sizə birini təklif edəcək. İstifadə etdiyiniz çevrilmənin tam adını yazın.

Nəhayət, bir çevrilməyə ehtiyac varsa, bu vasitəyə tam adını yazın, yenidən başladın və çıxış rasterinin istədiyiniz kimi proqnozlaşdırıldığını yoxlayın.

Orijinal rasteri yeni bir koordinat sisteminə yansıtmaq üçün istifadə ediləcək yenidən nümunə alqoritmi. ArcGIS Project Raster vasitəsi proyeksiyanı yerinə yetirmək üçün istifadə olunur və aşağıdakı dəyərləri qəbul edir:

NEAREST - ən yaxın qonşu interpolasiyası

BILINEAR - bilinear interpolasiya

Yeni bir koordinat sisteminə proyeksiya etmək üçün bu alqoritmlərdən birini göstərməlisiniz. Alqoritm seçmədən yeni bir koordinat sistemini təyin etsəniz bir səhv ortaya çıxacaq.

Proqnozlaşdırılan koordinat sisteminin hüceyrə ölçüsü. Bu parametr isteğe bağlı olsa da, ən yaxşı nəticələr əldə etmək üçün ArcGIS sənədləri yeni bir koordinat sisteminə proyeksiya edərkən hər zaman göstərməyinizi tövsiyə edir.

Piksel uyğunlaşdırılması üçün istifadə olunan x və y koordinatları (çıxış məkanında).

Bu parametr ArcGIS 9.2 tərəfindən təqdim edilən yeni bir seçimdir. Bu parametri istifadə etmək üçün ArcGIS 9.2 -yə sahib olmalısınız. Rasterin yeni bir koordinat sisteminə yansıtılması lazım olduğunu təyin etməsəniz, nəzərə alınmır.

Rasterin kəsilməli olduğu dikdörtgen.

Proqnozlaşdırılan bir koordinat sistemi göstərildisə, kəsmə proyeksiyadan sonra həyata keçirilir və yeni koordinat sistemində düzbucağın koordinatları göstərilməlidir. Heç bir proqnozlaşdırılan koordinat sistemi göstərilməyibsə, koordinatlar orijinal koordinat sistemində göstərilməlidir.

ArcGIS Clip vasitəsi klipi tamamlamaq üçün istifadə olunur. Kəsmə düzbucağı bu vasitəyə boşluqlarla ayrılmış dörd ədəddən ibarət bir simli olaraq verilməlidir. ArcGIS istifadəçi interfeysi, ArcGIS UI -dən bu vasitəni çağırarkən sətri avtomatik olaraq düzgün şəkildə formatlayır, format haqqında narahat olmağa ehtiyac yoxdur. Ancaq onu proqramatik olaraq çağırarkən, düzgün biçimlənmiş bir sətir təmin etməyə çalışın. Nömrələr SOL, ALT, SAĞ, TOP sıralanır. Məsələn, raster coğrafi koordinat sistemindədirsə, sim ilə 10 W, 15 S, 20 E və 25 N -ə kəsilə bilər:

Tam ədədlər və ya onluq ədədlər verilə bilər.

Çıxış rasterində icra etmək üçün cəbr ifadəsini xəritəyə salın.

XƏBƏRDARLIQ: ArcGIS Geoprocessing Model Builder, bu parametrin dəyərini təsadüfi və səssizcə silə bilər. Bu ArcGIS -də bir səhvdir. Saxladığınız bir modeli işə salmadan əvvəl bu vasitəni açın və parametr dəyərinin hələ də mövcud olduğunu yoxlayın.

İfadə, çevrilmiş raster proqnozlaşdırıldıqdan və kəsildikdən sonra icra olunur (əgər bu seçimlər göstərilirsə). Harita cəbrini yerinə yetirmək istədiyiniz rasteri təmsil etmək üçün hərflərə həssas olan inputRaster istifadə edin. Məsələn, rasteri tam ədəd rasterinə çevirmək və bütün hüceyrələrə 1 əlavə etmək üçün bu ifadəni istifadə edin:

İnputRaster simli hərflərə həssasdır. Xəritə cəbri ifadəsini yerinə yetirməzdən əvvəl, sim yaradılan çıxış rasterini təmsil edən müvəqqəti rasterə gedən yolla əvəz olunur. Son ifadə 4000 simvoldan az olmalıdır, əks halda ArcGIS bir səhv bildirəcək.

ArcGIS Tək Çıxış Xəritə Cəbri vasitəsi, xəritə cəbri ifadəsini yerinə yetirmək üçün istifadə olunur. Xəritə cəbrini yerinə yetirmək üçün ArcGIS Spatial Analyst uzantısı üçün lisenziyanız olmalıdır.

Xəritə cəbrinin sintaksisi çox seçici ola bilər. Bu vasitə ilə uğur qazanmağınıza kömək edəcək bəzi məsləhətlər:

Bu vasitəni istifadə etməzdən əvvəl, ArcGIS Tək Çıxış Xəritə Cəbri vasitəsi ilə xəritə cəbr ifadənizi qurun və sınayın. Sonra ifadəni bu vasitəyə yapışdırın və Tək Çıxış Xəritə Cəbrində istifadə etdiyiniz test dəyərindən çox inputRaster dəyişənini istifadə etmək üçün redaktə edin.

İfadənizi birbaşa bu vasitədə inkişaf etdirirsinizsə, çox sadə bir ifadə ilə başlayın. Düzgün işlədiyini yoxlayın, bir az əlavə edin və yenidən yoxlayın. Tam ifadəni qurana qədər bu prosesi təkrarlayın.

Həmişə boşluqlardan istifadə edərək riyazi operatorları raster yollarından ayırın. Yuxarıdakı nümunədə / operatoru hər iki tərəfdə bir boşluq ehtiva edir. Bu modeli izləyin. Bəzi hallarda, ArcGIS, boşluqlardan istifadə edən raster yollarını operatorlardan ayırmayan raster cəbr ifadələrini işlədə bilməyəcək. Bildirilən səhv mesajı ümumiyyətlə bunun problem olduğunu göstərmir və onu izləmək çox sinir bozucu ola bilər.

Doğru olarsa, çıxış rasteri üçün piramidalar qurulacaq ki, bu da ArcGIS istifadəçi interfeysində ekran sürətini artıracaq. Bu, dönüşümdən sonrakı emalda edilən son addımdır.

SQL WHERE işlənəcək satır alt hissəsini ifadə edən bir ifadə. Bu parametr təmin edilmirsə, bütün sətirlər işlənəcək. Bu parametr təmin edilsə də, əsas verilənlər bazası WHERE bəndlərini dəstəkləmirsə, bir səhv qaldırılacaq.

Bu ifadənin dəqiq sintaksisi əsas verilənlər bazasından asılıdır. ESRI, aşağıdakı sintaksisdən istifadə edərək sahələrə istinad etməyi tövsiyə edir:

ArcInfo əhatə dairələrini, şəkil fayllarını, INFO cədvəllərini və ya dBASE cədvəllərini (.dbf faylları) sorğu -sual edirsinizsə, sahə adlarını SQL ifadəsinə ikiqat dırnaqlarla əlavə edin: "MY_FIELD".

Microsoft Access cədvəllərini və ya şəxsi məlumat bazası cədvəllərini sorğulayırsınızsa, sahə adlarını kvadrat mötərizədə əlavə edin: [MY_FIELD].

ArcSDE coğrafi verilənlər bazası cədvəllərini, ArcIMS xüsusiyyət sinifini və ya ArcIMS görüntü xidməti alt qatını sorğulayırsınızsa, sahə adlarını daxil etməyin: MY_FIELD.

Satırları sıralamaq üçün istifadə ediləcək sahələr (yəni SQL SELECT ifadəsinin ORDER BY bəndində göstərilən sütunlar). Heç bir sahə təmin edilmədikdə, satırlar əsas verilənlər bazası tərəfindən təyin olunan standart qaydada sıralanacaq. Bu parametr təmin edilsə də, bu kompüter ArcGIS 9.2 və ya daha yeni bir versiyada çalışmırsa və ya əsas verilənlər bazası ORDER BY bəndlərini dəstəkləmirsə, bir səhv ortaya çıxacaq.

ORDER BY sahələrini göstərməyə əlavə olaraq, hər bir sahə üçün çeşidləmə istiqamətini də göstərməlisiniz.

ORDER BY sahələri üçün çeşid istiqamətlərini təyin edən 'Artan' və ya 'Azalan' sətirlərin siyahısı. Bu parametr təmin edilsə də, bu kompüter ArcGIS 9.2 və ya daha yeni bir versiyada çalışmırsa və ya əsas verilənlər bazası ORDER BY bəndlərini dəstəkləmirsə, bir səhv ortaya çıxacaq.

Doğru olarsa, artıq mövcud olan rasters üçün işləmə atlanacaq.

Nisbi yollardan asılı olmaq üçün əsas yol.

Bir əsas yol təmin edilərsə, girişləri (və bu vasitənin çıxışları varsa, çıxışları) göstərən sahələrdən əldə edilən nisbi yollara əlavə olunacaq. Bir əsas yol verilmirsə, cədvəli ehtiva edən iş yerinə əvvəlcədən əlavə olunacaq.


10.1 -də təqdim olunan BandArithmetic funksiyası, raster verilənlər bazasının bantlarında arifmetik əməliyyat həyata keçirir. Daha çox məlumat üçün Band Aritmetik funksiyasına baxın.

BandArithmetic funksiyası üçün arqumentlər aşağıdakılardır:

  • Normallaşdırılmış Fərq Bitki İndeksi
    • NDVI = ((NIR - Qırmızı)/(NIR + Qırmızı))
    • BandIndexes: "NIR Red", məsələn, "4 3"
    • İstinad: Pinty, B. və Verstraete, M. M. 1992, "GEMI: peyklərdən qlobal bitki örtüyünü izləmək üçün qeyri-xətti indeks", Bitki Ekologiyası, Cild. 101, 15-20.
    • SAVI = ((NIR - Qırmızı) / (NIR + Qırmızı + L)) x (1 + L)
    • burada L yaşıl bitki örtüyünün miqdarını təmsil edir, məsələn, 0,5
    • "BandIndexes": "NIR Red L", məsələn, "4 3 0.33"
    • İstinad: Huete, A. R., 1988, "Torpağa uyğunlaşdırılmış bitki indeksi (SAVI)," Uzaqdan Ətraf Mühitin Algılanması, Cild 25, 295-309.
    • TSAVI = (s (NIR-s*Red-a))/(a*NIR+Red-a*s+X*(1+s^2))
    • burada a = torpaq xəttinin kəsilməsi, s = torpaq xəttinin yamacı, X = torpağın səs -küyünü minimuma endirmək üçün təyin edilmiş bir tənzimləmə faktoru
    • "BandIndexes": "NIR Red s a X", məsələn, "4 3 0.33 0.50 1.50"
    • İstinad: Baret, F. və G. Guyot, 1991, "LAI və APAR qiymətləndirilməsi üçün bitki indekslərinin potensialları və məhdudiyyətləri", Ətrafın Uzaqdan Algılanması, Cild. 35, 161-173.
    • GEMI = eta*(1-0.25*eta)-((Qırmızı-0.125)/(1-Qırmızı))
    • burada eta = (2*(NIR^2-Qırmızı^2)+1,5*NIR+0,5*Qırmızı)/(NIR+Qırmızı+0,5)
    • "BandIndexes": "NIR Red", məsələn, "4 3"
    • İstinad: Pinty, B. və Verstraete, M. M. 1992, "GEMI: peyklərdən qlobal bitki örtüyünü izləmək üçün qeyri-xətti indeks", Bitki Ekologiyası, Cild. 101, 15-20,
    • GVITM = -0.2848*Band1-0.2435*Band2-0.5436*Band3+0.7243*Band4+0.0840*Band5-1.1800*Band7
    • "BandIndexes": "Band1 Band2 Band3 Band4 Band5 Band7", məsələn, "1 2 3 4 5 6"
    • İstinad: Todd, S. W., R. M. Hoffer və D. G. Milchunas, 1998, "Spektral indekslərdən istifadə edərək otlaqlı və otlaqsız otlaqlarda biokütlə qiymətləndirilməsi", International Journal of Remote Sensing, Vol. 19, No 3, 427-438.
    • MSAVI2 = (1/2)*(2 (NIR+1) -sqrt ((2*NIR+1)^2-8 (NIR-Qırmızı)))
    • BandIndexes: "NIR Red", məsələn, "4 3"
    • İstinad: Qi, J. et al., 1994, "Dəyişdirilmiş torpaq bitki örtüyü düzəliş indeksi," Uzaqdan Ətraf Mühitin Algılanması, Cild. 48, No 2, 119-126.
    • PVI = (NIR-a*Qırmızı-b)/(sqrt (1+a^2))
    • burada a = torpaq xəttinin yamacı, b = torpaq xəttinin qradiyenti
    • BandIndexes: "NIR Red a b", məsələn, "4 3 0.3 0.5"
    • Bant 1 = (Bant5 / Bant6) x 100
    • Bant 2 = (Bant5 / Bant1) x 100
    • Bant 3 = (Bant3 / Bant4) x (Bant5 / Bant4) x 100
    • "BandIndexes": "Band1 Band2 Band3 Band4 Band5 Band6", məsələn, "1 2 3 4 5 6"
    • VARI = (Yaşıl - Qırmızı)/(Yaşıl + Qırmızı - Mavi)
    • BandIndexes: "Qırmızı Yaşıl Mavi", məsələn, "3 2 1"
    • İstinad: Anatoli A. Gitelson, Yoram J. Kaufman, Robert Stark və Don Rundquist, 2002, "Bitki fraksiyasının uzaqdan qiymətləndirilməsi üçün yeni alqoritmlər," Uzaqdan Ətraf Mühitin Algılanması, Cild. 80, 76-87.
    • Tək bantlı çıxış yaratmaq üçün tək sətirli cəbr formulundan istifadə edin. Dəstəklənən operatorlar -, +, /, *və unary -dir.
    • Bantları müəyyən etmək üçün bant nömrəsini B və ya b hərflərindən istifadə edin. Məsələn: "BandIndexes": "(B1 + B2) / (B3 * B5)"

    Misal


    Videoya baxın: Azərbaycan bayraqlı harita Azərbaycan bayraxlı xəritə