Lebih

Kesalahan saat Memproyeksikan Raster menggunakan ArcPy?

Kesalahan saat Memproyeksikan Raster menggunakan ArcPy?


Menggunakan ArcGIS 10.3 Saya menerima 'ERROR 999999' yang ditakuti ketika mencoba memproyeksikan geodatabase yang penuh dengan raster dari WGS84_UTM_zone_10N ke Nad83_BC_Environment_Albers

import arcpy InRasterGDB = r"C:UsersPoirierMDocumentsGISDataLandsatSSTImagery.gdb" OutRasterGDB = r"C:UsersPoirierMDocumentsGISDataLandsatSSTImageryAlbers.gdb" OutCoordSystem = "3005 CellSize = "30 30" GeoTransform = "WGS_1984_(ITRF00)_To_NAD_1983" InCoordSystem = "32610" arcpy.env.workspace = InRasterGDB untuk nama f di arcpy.ListRasters(): print fname print OutRasterGDB + "" + fname arcpy.ProjectRaster_management (InRasterGDB + "" + fname, OutRasterGDB + "" + fname, OutCoordSystem, ResampleType, CellSize, GeoTransform, "", InCoordSystem)

Ini seharusnya berfungsi, tetapi mengembalikan kesalahan berikut saat dijalankan di Jendela Python di ArcCatalog:

Runtime error Traceback (panggilan terakhir terakhir): File "", baris 16, di  File "c:program files (x86)arcgisdesktop10.3arcpyarcpymanagement.py", baris 9009, di ProjectRaster, naikkan e ExecuteError: ERROR 999999: Error execution function. Gagal mengeksekusi (ProjectRaster).

Saya dapat melakukan satu raster secara manual baik melalui antarmuka alat atau dalam pembuat model. Adakah yang bisa menyarankan mengapa itu tidak berfungsi sebagai skrip Python?


Saya punya skrip yang berfungsi sekarang - lihat jawaban (MPoirier) saya di bawah ini.


Saya percaya masalahnya mungkin dalam definisi sistem koordinat Anda - ref. halaman bantuan Project Raster.

Sistem koordinat tempat raster input akan diproyeksikan. Nilai default diatur berdasarkan pengaturan lingkungan Sistem Koordinat Output.

Nilai yang valid untuk parameter ini adalah

  • File dengan ekstensi ".prj" (file prj yang dikirimkan bersama ArcGIS dapat ditemukan di "C:Program FilesArcGISCoordinate Systems").
  • Kelas fitur yang ada, kumpulan data fitur, katalog raster (pada dasarnya apa saja dengan sistem koordinat).
  • Representasi string dari sistem koordinat. String panjang ini dapat dibuat dengan menambahkan variabel sistem koordinat ke ModelBuilder, menyetel nilai variabel sesuai keinginan, lalu mengekspor model ke skrip Python.

Anda memberinya wkid secara langsung, yang tidak terdaftar sebagai opsi.

Tapi, kamu bisa membuat referensi spasial melalui wkid menggunakan properti SpatialReference.

OutCoordSystem = arcpy.SpatialReference(3005) InCoordSystem = arcpy.SpatialReference(32610)

Berikut ini skrip yang berfungsi. Saya mulai dari kode pembuat model untuk memproyeksikan ulang satu raster, dan secara bertahap mengubah input ke pengujian variabel di setiap langkah. Tampaknya satu-satunya perbedaan antara skrip yang berfungsi di bawah ini dan skrip saya yang tidak berfungsi sebelumnya adalah bahwa variabel didefinisikan di dalam badan for loop.

Saya juga menambahkan beberapa komentar untuk memperjelas.

##Script untuk mengulang melalui file geodatabase untuk menemukan raster dan memproyeksikannya kembali dari WGS_1984_UTM_zone_10N ke NAD_1983_BC_Environment_Albers ## Matthew Poirier, Terakhir diubah 18 Februari 2015 # Impor modul arcpy impor arcpy #Define workspace (lokasi geodatabase) InRasterC: = r" PoirierMDocumentsGISDataLandsatSSTImagery.gdb" arcpy.env.workspace = InRasterGDB #Loop melalui daftar semua raster di ruang kerja untuk Ras di arcpy.ListRasters(): # Variabel lokal untuk input (opsional) dan output: InRas = "C:UsersPoirierMDocumentsGISDataLandsatSSTImagery.gdb" + Ras OutRas = "C:UsersPoirierMDocumentsGISDataLandsatSSTImageryAlbers.gdb" + Ras # Tentukan string proyeksi - bisa juga menggunakan string bentuk panjang yang disalin dari skrip pembuat model (diuji) atau file .prj yang ada InProj = "PROJCS['WGS_1984_UTM_zone_10N',GEOGCS['GCS_WGS_1984',DATUM['D_WGS_1984',SPHEROID[' WGS_1984',6378137.0,298.257223563]],PRIMEM['Greenwich',0.0],UNIT['Gelar',0.0174532925199433 ]],PROJECTION['Transverse_Mercator'],PARAMETER['false_easting',500000.0],PARAMETER['false_northing',0.0],PARAMETER['central_meridian',-123.0],PARAMETER['scale_factor',0.9996],PARAMETER[' latitude_of_origin',0.0],UNIT['Meter',1.0]]" OutProj = "PROJCS['NAD_1983_BC_Environment_Albers',GEOGCS['GCS_North_American_1983',DATUM['D_North_American_1983',SPHEROID[',GRS_1983] ['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Albers'],PARAMETER['False_Easting',1000000.0],PARAMETER['False_Northing',0.0],PARAMETER['Central_Meridian',- 126.0],PARAMETER['Standard_Parallel_1',50.0],PARAMETER['Standard_Parallel_2',58.5],PARAMETER['Latitude_Of_Origin',45.0],UNIT['Meter',1.0]]" # Tentukan metode transformasi datum geografis yang akan digunakan GeoTran = "WGS_1984_(ITRF00)_To_NAD_1983" # Tentukan jenis resampling (bilinear atau kubik untuk data kontinu, default ke terdekat) ReSample = "BILINEAR" # Tentukan ukuran sel output (default ke raster input Cel lSize = "30 30" # Proses: Project Raster - perhatikan tanda kutip kosong adalah pengganti untuk parameter opsional 'registration_point' arcpy.ProjectRaster_management(InRas, OutRas, OutProj, ReSample, CellSize, GeoTran, "", InProj)

Coba lewati ukuran sel sebagai satu dimensi.

Suka:

Ukuran Sel = "30"

Jika tidak, coba uji skrip tanpa parameter opsional apa pun dan lihat apa yang terjadi.


Menghubungkan ikan ke beberapa metrik kompleksitas struktural terumbu karang menggunakan teknologi tiga dimensi

Kompleksitas struktural sangat mempengaruhi keanekaragaman hayati dan produktivitas ekosistem. Pada terumbu karang, kompleksitas struktural biasanya diukur menggunakan metrik skala tunggal dan kecil ('kekasaran') yang mewakili beberapa atribut spasial yang dieksploitasi secara berbeda oleh spesies, sehingga membatasi pemahaman lengkap tentang bagaimana ikan berasosiasi dengan struktur terumbu. Kami menggunakan pendekatan baru untuk membandingkan hubungan antara ikan dan komponen kompleksitas terumbu yang sebelumnya tidak tersedia, dan membandingkan hasilnya dengan indeks rugositas tradisional. Penelitian ini berfokus pada ikan damsel untuk mengeksplorasi hubungan antara ikan dan struktur terumbu. Tiga spesies teritorial, dengan kebiasaan trofik yang kontras dan penggunaan struktur terumbu yang diharapkan, diperiksa untuk menyimpulkan mekanisme spesifik spesies potensial yang terkait dengan bagaimana kompleksitas memengaruhi pemilihan habitat. Rekonstruksi terumbu tiga dimensi dari fotogrametri mengukur metrik kualitas habitat berikut: 1) paparan visual terhadap predator dan pesaing, 2) kepadatan tempat perlindungan predator dan 3) ketersediaan makanan terkait substrat. Metrik ini menjelaskan distribusi spesies lebih baik daripada ukuran rugositas tradisional, dan setiap spesies merespons komponen kompleksitas yang berbeda. Mengingat bahwa efek kritis dari degradasi terumbu adalah hilangnya struktur, mengadopsi teknologi tiga dimensi berpotensi menawarkan alat baru untuk memahami asosiasi spesies-habitat dan membantu memperkirakan bagaimana ikan akan terpengaruh oleh perataan terumbu.


Menabrak beranda oleh Komunitas ♦ 2 hari yang lalu

Pertanyaan ini memiliki jawaban yang mungkin baik atau buruk sistem telah menandainya aktif sehingga mereka dapat ditinjau.

Tidak yakin kami dapat menguji ini tanpa setidaknya satu sampel HDF - dari mana Anda mendapatkannya? Bisakah Anda berbagi satu?

Contoh File: drive.google.com/file/d/1ZQqF8IHds8uxHqevxTee4JQmZ8I3rOpU/…

– Ritesh Porwal
28 Februari pukul 10:44

Gunakan gdal_translate dengan parameter -tr

saya sudah menggunakan parameter -tr tetapi tidak berhasil.

– Ritesh Porwal
28 Februari pukul 11:58

Saya mencoba mengekstrak subdataset dari MODIS LST HDF menggunakan kode berikut.

Tetapi raster Output (sds4.tif) memiliki ukuran sel yang berbeda [929.378147, 907]. [926.6254331, 926.6254331] adalah ukuran sel asli dari file HDF.


Ytdjtkui

Sebuah eksperimen sosial. Apa hal terburuk yang bisa terjadi?

Acara yang dikenal oleh pecandu sejarah

Apakah roda pendarat yang lebih tinggi buruk untuk pesawat, khususnya pesawat besar?

Lingkungan kerja yang tidak bersahabat setelah whistle-blowing pada rekan kerja dan bos kita. Apa yang saya lakukan?

Balas ‘no posisi’ selagi lowongan masih ada (‘HiWi’ posisi di Jerman)

Bagaimana cara memeriksa peserta di acara?

Apakah boleh / apakah masuk akal bagi pemain lain untuk bergabung dengan permainan lari Munchkin?

Bagaimana jika seseorang berinvestasi dalam aplikasi saya?

Mengajarkan integral tak tentu yang memerlukan selubung khusus

Mengapa bagian landasan peluncuran Pesawat Ulang-alik ini tampak mengambang di udara?

Bisakah saya membuat dinding setinggi 7 kaki × 5 kaki dengan mantra Ilusi Kecil?

perintah node saat mendefinisikan koordinat di TikZ

Apakah sah untuk melakukan diskriminasi karena obat yang digunakan untuk mengobati kondisi medis?

Mensimulasikan probabilitas 1 dari 2^N dengan kurang dari N bit acak

Menginstal PowerShell pada Kali OS 32-bit gagal

Mengapa semua pintu di Ferenginar (dunia asal Ferengi) jauh lebih pendek daripada rata-rata Ferengi?

Apakah ada situs perbandingan Kesepakatan Brexit yang Tidak Berpihak?

Bisakah addon jahat mengakses riwayat internet dan semacamnya di chrome/firefox?

Giant Toughroad SLR 2 sejauh 200 mil dalam dua hari, apakah akan berhasil?

Sinyal I2C dan daya pada jarak jauh (kabel 10 meter)

Apakah tak terhingga dapat diamati secara matematis?

Bisakah saya mengandalkan file repositori GitHub ini?

Apa istilah ketika dua orang bernyanyi dalam harmoni, tetapi mereka tidak menyanyikan nada yang sama?

R - Bagaimana memproyeksikan ulang koordinat sedemikian rupa sehingga dihitung dalam Meter?

Bagaimana mengaktifkan transformasi proyeksi di QGIS2 Vektor Shapefile tidak akan diplot dengan benar di ruang yang sama Memproyeksikan ulang Google Maps (latlong) ke Irish Grid TM65 menggunakan spTransform di R? Perbedaan besar antara hasil untuk transformasi datum 7-param Terjemahan ringan dengan transformasi koordinat dari Gauss-Kruger ke WGS 84. Kemungkinan Alasan?Ekspor/cetak QGis: ketidakcocokan data OSM pada OSM-layerCRS Masalah proyeksi mencegah penggabungan spasial di RNumerical Vertex Edit Error: lintang atau bujur melebihi batasMasalah menggunakan 'masker' saat memplot KDE - RUnit dari proyeksi CRS tertentu?

Saya memiliki satu set poin dari pelacakan hewan yang berbeda dan saya ingin mencoba dan memperkirakan wilayah jelajah mereka (KED, MCP). Saya berhasil menyelesaikan perhitungan tetapi output unit tampak aneh: sepertinya unit dalam derajat .

CRS dari Koordinat adalah EPSG:4326 (Luas: -180.00, -90.00, 180.00, 90.00
Proj4: +proj=longlat +datum=WGS84 +no_defs) .

Ini adalah bagaimana saya membaca data di R.

Judul data saya pada dasarnya adalah:

untuk menghitung MCP yang saya gunakan

Adakah yang bisa membantu saya untuk memproyeksikan ulang ini ke beberapa CRS lain sehingga perhitungan luas akhir dalam meter? (masih menemukan kakiku dengan R)

Sudahkah Anda membaca data ke dalam R? Bisakah Anda menunjukkan kepada kami caranya? Ada paket untuk data spasial di R dan kecuali Anda memberi tahu kami apa yang Anda gunakan, tidak mudah bagi kami untuk membantu. Misalnya jika Anda memiliki data poligon sf, maka st_area(d) akan mengembalikan luas dalam meter persegi meskipun data dalam derajat lat-long. Harap edit pertanyaan Anda untuk memberikan beberapa informasi lebih lanjut.

terima kasih. Tolong saya mengedit posting awal untuk menambahkan beberapa informasi lebih lanjut.

Saya memiliki satu set poin dari pelacakan hewan yang berbeda dan saya ingin mencoba dan memperkirakan wilayah jelajah mereka (KED, MCP). Saya berhasil menyelesaikan perhitungan tetapi output unit tampak aneh: sepertinya unit dalam derajat .

CRS dari Koordinat adalah EPSG:4326 (Luas: -180.00, -90.00, 180.00, 90.00
Proj4: +proj=longlat +datum=WGS84 +no_defs) .

Ini adalah bagaimana saya membaca data di R.

Judul data saya pada dasarnya adalah:

untuk menghitung MCP yang saya gunakan

Adakah yang bisa membantu saya untuk memproyeksikan ulang ini ke beberapa CRS lain sehingga perhitungan luas akhir dalam meter? (masih menemukan kakiku dengan R)

Sudahkah Anda membaca data ke dalam R? Bisakah Anda menunjukkan kepada kami caranya? Ada paket untuk data spasial di R dan kecuali Anda memberi tahu kami apa yang Anda gunakan, tidak mudah bagi kami untuk membantu. Misalnya jika Anda memiliki data poligon sf, maka st_area(d) akan mengembalikan luas dalam meter persegi meskipun data dalam derajat lat-long. Harap edit pertanyaan Anda untuk memberikan beberapa informasi lebih lanjut.

terima kasih. Tolong saya mengedit posting awal untuk menambahkan beberapa informasi lebih lanjut.

Saya memiliki satu set poin dari pelacakan hewan yang berbeda dan saya ingin mencoba dan memperkirakan wilayah jelajah mereka (KED, MCP). Saya berhasil menyelesaikan perhitungan tetapi output unit tampak aneh: sepertinya unit dalam derajat .

CRS dari Koordinat adalah EPSG:4326 (Luas: -180.00, -90.00, 180.00, 90.00
Proj4: +proj=longlat +datum=WGS84 +no_defs) .

Ini adalah bagaimana saya membaca data di R.

Judul data saya pada dasarnya adalah:

untuk menghitung MCP yang saya gunakan

Adakah yang bisa membantu saya untuk memproyeksikan ulang ini ke beberapa CRS lain sehingga perhitungan luas akhir dalam meter? (masih menemukan kakiku dengan R)

Saya memiliki satu set poin dari melacak hewan yang berbeda dan saya ingin mencoba dan memperkirakan wilayah jelajah mereka (KED, MCP). Saya berhasil menyelesaikan perhitungan tetapi output unit tampak aneh: sepertinya unit dalam derajat .

CRS dari Koordinat adalah EPSG:4326 (Luas: -180.00, -90.00, 180.00, 90.00
Proj4: +proj=longlat +datum=WGS84 +no_defs) .

Ini adalah bagaimana saya membaca data di R.

Judul data saya pada dasarnya adalah:

untuk menghitung MCP yang saya gunakan

Adakah yang bisa membantu saya untuk memproyeksikan ulang ini ke beberapa CRS lain sehingga perhitungan luas akhir dalam meter? (masih menemukan kakiku dengan R)

Sudahkah Anda membaca data ke dalam R? Bisakah Anda menunjukkan kepada kami caranya? Ada paket untuk data spasial di R dan kecuali Anda memberi tahu kami apa yang Anda gunakan, tidak mudah bagi kami untuk membantu. Misalnya jika Anda memiliki data poligon sf, maka st_area(d) akan mengembalikan luas dalam meter persegi meskipun data dalam derajat lat-long. Harap edit pertanyaan Anda untuk memberikan beberapa informasi lebih lanjut.

terima kasih. Tolong saya mengedit posting awal untuk menambahkan beberapa informasi lebih lanjut.

Sudahkah Anda membaca data ke dalam R? Bisakah Anda menunjukkan kepada kami caranya? Ada paket untuk data spasial di R dan kecuali Anda memberi tahu kami apa yang Anda gunakan, tidak mudah bagi kami untuk membantu. Misalnya jika Anda memiliki data poligon sf, maka st_area(d) akan mengembalikan luas dalam meter persegi meskipun data dalam derajat lat-long. Harap edit pertanyaan Anda untuk memberikan beberapa informasi lebih lanjut.

terima kasih. Tolong saya mengedit posting awal untuk menambahkan beberapa informasi lebih lanjut.

Sudahkah Anda membaca data ke dalam R? Bisakah Anda menunjukkan kepada kami caranya? Ada paket untuk data spasial di R dan kecuali Anda memberi tahu kami apa yang Anda gunakan, tidak mudah bagi kami untuk membantu. Misalnya jika Anda memiliki data poligon sf, maka st_area(d) akan mengembalikan luas dalam meter persegi meskipun data dalam derajat lat-long. Harap edit pertanyaan Anda untuk memberikan beberapa informasi lebih lanjut.

Sudahkah Anda membaca data ke dalam R? Bisakah Anda menunjukkan kepada kami caranya? Ada paket untuk data spasial di R dan kecuali Anda memberi tahu kami apa yang Anda gunakan, tidak mudah bagi kami untuk membantu. Misalnya jika Anda memiliki data poligon sf, maka st_area(d) akan mengembalikan luas dalam meter persegi meskipun data dalam derajat lat-long. Harap edit pertanyaan Anda untuk memberikan beberapa informasi lebih lanjut.


Ytdjtkui

Mengapa komputer on-board diizinkan untuk mengubah kontrol tanpa memberi tahu pilot?

Mengajarkan integral tak tentu yang memerlukan selubung khusus

Bisakah Saya Mengambil Alamat Email dari BCC?

Memetakan daftar menjadi plot fase

Bagaimana cara memverifikasi apakah g adalah generator untuk p?

Menemukan semua interval yang cocok dengan predikat dalam vektor

Bagaimana Bumi seorang diri mampu menciptakan 3 dari 4 dewa kekacauan?

Irreduksibilitas polinomial sederhana simple

Penerapan Prinsip Tanggung Jawab Tunggal

Apakah saya memerlukan visa masuk ganda untuk perjalanan Inggris -> Swedia -> Inggris?

Kutukan bekerja dengan berteriak - Bagaimana cara menghindari kerusakan tambahan?

Lingkungan kerja yang tidak bersahabat setelah whistle-blowing pada rekan kerja dan bos kita. Apa yang saya lakukan?

Pada titik mana karakter mendapatkan kembali semua Dadu Hit mereka?

Menggabungkan pegangan cabang fitur paralel secara efisien di SFDX

Bagaimana saya bisa menggunakan tanda panah di bash Prompt saya?

Apakah area gambar CRT adalah jajaran genjang (bukan persegi panjang yang sebenarnya)?

Mengapa Kant, Hegel, dan Adorno meninggalkan beberapa kata dan frasa dalam alfabet Yunani?

Apakah panjang goresan Kanji yang tepat penting?

Bagaimana saya bisa melewati dokumen teknis yang sangat panjang dan sangat kering, tetapi juga sangat berguna saat mempelajari alat baru?

Haruskah tesis PhD saya diajukan dengan nama resmi saya?

Bisakah penipuan kriminal ada tanpa kerusakan?

Apakah ada alasan untuk tidak memakan makanan yang dijatuhkan di permukaan bulan?

Cara untuk mempercepat pengguna menerapkan RK4

Konversi file GeoTIFF ke NetCDF tidak membuat file NetCDF tetapi database saya menganggapnya demikian

Konversi file NetCDF ke GeotiffProgram/ bahasa untuk memproses file netcdf “cepat”Menyandikan proyeksi azimut yang berjarak sama dalam konversi NetCDFcara menyiapkan tiff untuk membuat file netcdfAdakah yang tahu cara mengatasi masalah GDAL dengan raster penerjemahan batch di QGIS?Tentukan geotransformasi yang akan dikonversi aset impor mesin netCDF ke geotiffgoogle earth: masalah dengan proyeksiMengonversi netCDF ke file geotiff?Mengonversi file GeoTiff ke netCDFMenggabungkan NetCDF yang memiliki dimensi identik

Saya menggunakan miniconda dan proyek Open Data Cube. Saya mencoba untuk mencerna beberapa data Landsat_5 ke dalam datacube, di mana geotiff dikonversi ke file netCDF. Menggunakan kode datacube saya menggunakan alat menelan:

Itu membuat file NetCDF dan terlihat benar, namun ketika menjalankan kode notebook jupyter saya di mana saya memuat dataset landsat dengan:

Saya menemukan kesalahan saat mencari LS5. 6_48. nc file yang tidak ada. Saya memiliki LS5. 6_49. nc , 6_50, dan seterusnya, tetapi bukan 6_48. Itu terjadi untuk banyak file dan bukan hanya satu. Saya pikir mungkin ini kesalahan saat mengonversi geotiff ke NetCDF tetapi tidak yakin bagaimana cara memperbaikinya. Saya sudah mencoba memastikan perpustakaan saya diperbarui dan semuanya terlihat bagus. Kesalahan adalah sebagai berikut:

Ada saran? Saya telah mengurangi gambar saya menjadi 15 dari 8000 dan telah sepenuhnya menjatuhkan datacube dan membuat yang baru hanya dengan 15 itu dan masih menghadapi masalah yang sama.

Saya menggunakan miniconda dan proyek Open Data Cube. Saya mencoba untuk mencerna beberapa data Landsat_5 ke dalam datacube, di mana geotiff dikonversi ke file netCDF. Menggunakan kode datacube saya menggunakan alat menelan:

Itu membuat file NetCDF dan terlihat benar, namun ketika menjalankan kode notebook jupyter saya di mana saya memuat dataset landsat dengan:

Saya menemukan kesalahan saat mencari LS5. 6_48. nc file yang tidak ada. Saya memiliki LS5. 6_49. nc , 6_50, dan seterusnya, tetapi bukan 6_48. Itu terjadi untuk banyak file dan bukan hanya satu. Saya pikir mungkin ini kesalahan saat mengonversi geotiff ke NetCDF tetapi tidak yakin bagaimana cara memperbaikinya. Saya sudah mencoba memastikan perpustakaan saya diperbarui dan semuanya terlihat bagus. Kesalahan adalah sebagai berikut:

Ada saran? Saya telah mengurangi gambar saya menjadi 15 dari 8000 dan telah sepenuhnya menjatuhkan datacube dan membuat yang baru hanya dengan 15 itu dan masih menghadapi masalah yang sama.

Saya menggunakan miniconda dan proyek Open Data Cube. Saya mencoba untuk mencerna beberapa data Landsat_5 ke dalam datacube, di mana geotiff dikonversi ke file netCDF. Menggunakan kode datacube saya menggunakan alat menelan:

Itu membuat file NetCDF dan terlihat benar, namun ketika menjalankan kode notebook jupyter saya di mana saya memuat dataset landsat dengan:

Saya menemukan kesalahan saat mencari LS5. 6_48. nc file yang tidak ada. Saya memiliki LS5. 6_49. nc , 6_50, dan seterusnya, tetapi bukan 6_48. Itu terjadi untuk banyak file dan bukan hanya satu. Saya pikir mungkin ini kesalahan saat mengonversi geotiff ke NetCDF tetapi tidak yakin bagaimana cara memperbaikinya. Saya sudah mencoba memastikan perpustakaan saya diperbarui dan semuanya terlihat bagus. Kesalahan adalah sebagai berikut:

Ada saran? Saya telah mengurangi gambar saya menjadi 15 dari 8000 dan telah sepenuhnya menjatuhkan datacube dan membuat yang baru hanya dengan 15 itu dan masih menghadapi masalah yang sama.

Saya menggunakan miniconda dan proyek Open Data Cube. Saya mencoba untuk mencerna beberapa data Landsat_5 ke dalam datacube, di mana geotiff dikonversi ke file netCDF. Menggunakan kode datacube saya menggunakan alat menelan:

Itu membuat file NetCDF dan terlihat benar, namun ketika menjalankan kode notebook jupyter saya di mana saya memuat dataset landsat dengan:

Saya menemukan kesalahan saat mencari LS5. 6_48. nc file yang tidak ada. Saya memiliki LS5. 6_49. nc , 6_50, dan seterusnya, tetapi bukan 6_48. Itu terjadi untuk banyak file dan bukan hanya satu. Saya pikir mungkin ini kesalahan saat mengonversi geotiff ke NetCDF tetapi tidak yakin bagaimana cara memperbaikinya. Saya sudah mencoba memastikan perpustakaan saya diperbarui dan semuanya terlihat bagus. Kesalahan adalah sebagai berikut:

Ada saran? Saya telah mengurangi gambar saya menjadi 15 dari 8000 dan telah sepenuhnya menjatuhkan datacube dan membuat yang baru hanya dengan 15 itu dan masih menghadapi masalah yang sama.


1 Jawaban 1

Saya telah diberi jawaban oleh Esri untuk ini dan mencantumkannya di sini kalau-kalau ada orang yang menemukan masalah yang sama.

ArcGIS Pro mengambil sampel dari data untuk menemukan nilai min-max. Sampel ini secara default diatur ke 10.000.

Jika pengguna mengubah ukuran sampel tersebut melalui menu overflow (hamburger) pada tab Symbology ke angka yang lebih tinggi, itu akan membaca nilai min-max dengan benar. Dalam kasus saya, saya harus mengatur ukuran sampel ke 163000 karena begitu sedikit baris dari keseluruhan yang lebih tinggi dari 661.

Alasan ia bekerja di QGIS tanpa ini tampaknya karena QGIS tidak menggunakan ukuran sampel, ia hanya menghitung min-max data dari seluruh bingkai data secara default.

Sangat menarik. Saya tidak menyadari bahwa ArcGIS Pro menggunakan sampel data untuk menghitung nilai min-max untuk simbologi.


Bdfyjk

Saya mencoba mengubah file GRIB menjadi GeoTIFF untuk digunakan di ArcGIS, tetapi saya kesulitan menampilkan gambar di lokasi yang benar. Saya telah dapat membuat GeoTIFF, menggunakan GDAL dengan Python, yang menunjukkan data tetapi tidak muncul di lokasi yang benar saat dibawa ke ArcGIS. Gambar yang dihasilkan di bawah ini.

Data yang saya kerjakan dapat diunduh dari: https://gimms.gsfc.nasa.gov/SMOS/SMAP/L05/

Saya mencoba memproyeksikan data ke WGS84 Web Mercator (Auxiliary Sphere), EPSG: 3857

Catatan: Saya telah mencoba memasukkan data melalui ArcMap dengan membuat Raster Mosaic yang seharusnya dapat bekerja dengan data GRIB, tetapi saya tidak beruntung. Saya juga sudah mencoba menggunakan alat Project Raster, tetapi ArcGIS tidak menyukai proyeksi default yang berasal dari file GRIB dan memberikan kesalahan.

Saya tidak terlalu berpengalaman dalam menggunakan GDAL atau GDAL dengan Python.

Saya mencoba mengubah file GRIB menjadi GeoTIFF untuk digunakan di ArcGIS, tetapi saya kesulitan menampilkan gambar di lokasi yang benar. Saya telah dapat membuat GeoTIFF, menggunakan GDAL dengan Python, yang menunjukkan data tetapi tidak muncul di lokasi yang benar saat dibawa ke ArcGIS. Gambar yang dihasilkan di bawah ini.

Data yang saya kerjakan dapat diunduh dari: https://gimms.gsfc.nasa.gov/SMOS/SMAP/L05/

Saya mencoba memproyeksikan data ke WGS84 Web Mercator (Auxiliary Sphere), EPSG: 3857

Catatan: Saya telah mencoba memasukkan data melalui ArcMap dengan membuat Raster Mosaic yang seharusnya dapat bekerja dengan data GRIB, tetapi saya tidak beruntung. Saya juga sudah mencoba menggunakan alat Project Raster, tetapi ArcGIS tidak menyukai proyeksi default yang berasal dari file GRIB dan memberikan kesalahan.

Saya tidak terlalu berpengalaman dalam menggunakan GDAL atau GDAL dengan Python.

Saya mencoba mengubah file GRIB menjadi GeoTIFF untuk digunakan di ArcGIS, tetapi saya kesulitan menampilkan gambar di lokasi yang benar. Saya telah dapat membuat GeoTIFF, menggunakan GDAL dengan Python, yang menunjukkan data tetapi tidak muncul di lokasi yang benar saat dibawa ke ArcGIS. Gambar yang dihasilkan di bawah ini.

Data yang saya kerjakan dapat diunduh dari: https://gimms.gsfc.nasa.gov/SMOS/SMAP/L05/

Saya mencoba memproyeksikan data ke WGS84 Web Mercator (Auxiliary Sphere), EPSG: 3857

Catatan: Saya telah mencoba memasukkan data melalui ArcMap dengan membuat Raster Mosaic yang seharusnya dapat bekerja dengan data GRIB, tetapi saya tidak beruntung. Saya juga sudah mencoba menggunakan alat Project Raster, tetapi ArcGIS tidak menyukai proyeksi default yang berasal dari file GRIB dan memberikan kesalahan.

Saya tidak terlalu berpengalaman dalam menggunakan GDAL atau GDAL dengan Python.

Saya mencoba mengubah file GRIB menjadi GeoTIFF untuk digunakan di ArcGIS, tetapi saya kesulitan menampilkan gambar di lokasi yang benar. Saya telah dapat membuat GeoTIFF, menggunakan GDAL dengan Python, yang menunjukkan data tetapi tidak muncul di lokasi yang benar saat dibawa ke ArcGIS. Gambar yang dihasilkan di bawah ini.

Data yang saya kerjakan dapat diunduh dari: https://gimms.gsfc.nasa.gov/SMOS/SMAP/L05/

Saya mencoba memproyeksikan data ke WGS84 Web Mercator (Auxiliary Sphere), EPSG: 3857

Catatan: Saya telah mencoba memasukkan data melalui ArcMap dengan membuat Mosaik Raster yang seharusnya dapat bekerja dengan data GRIB, tetapi saya tidak beruntung. Saya juga sudah mencoba menggunakan alat Project Raster, tetapi ArcGIS tidak menyukai proyeksi default yang berasal dari file GRIB dan memberikan kesalahan.


1 Jawaban 1

Karena data Anda melintasi dua zona UTM 36N dan 37N, dan dengan asumsi Anda memiliki dua gambar, Anda memiliki tiga opsi:

  • Jika kedua gambar terpisah dan setiap gambar terletak di zona yang berbeda, maka proyeksikan ulang masing-masing dengan zona UTM yang benar.
  • Jika dua gambar digabungkan bersama (gambar mosaik), maka pilih salah satu zona seperti UTM 37N, karena sebagian besar Yordania terletak di zona UTM 37N dan hanya area Aqaba yang terletak di UTM 36N, yang dalam hal ini tidak akan deformasi besar.
  • Opsi terakhir adalah Anda dapat menggunakan sistem referensi koordinasi (CRS) EPSG:3066 Jordan TM yang merupakan CRS UTM khusus yang sesuai dengan Yordania dalam satu zona.

Pilihannya terserah anda mau menggunakan UTM mengikuti standar CRS atau CRS lokal yaitu Jordan TM CRS.


Vjylku

Mengapa asam asetat (pKa = 4,76) lebih kuat dari asam karbonat (pKa = 6,36)?

Apakah Artikel Ini Tentang Kemungkinan Ilmu Sampah Alam Semesta Cermin?

Apakah mungkin memiliki planeswalker yang sama dari edisi yang berbeda di dek Komandan?

Mengapa Bernie Sanders selalu disebut "sosialis"?

Produk Minimum yang Layak untuk game RTS?

Apakah ada budaya selebriti di dunia akademis dan haruskah kita mencegahnya?

Bagaimana cara menghadapi bawahan yang berkinerja buruk?

Ruang disk penuh selama penyisipan, apa yang terjadi?

Bagaimana cara menambahkan "rasa bawang" yang kuat ke biryani (dalam gaya restoran)?

Lari dari konflik

Serial fantasi dengan dewa yang bersembunyi di depan mata sebagai pelayan pria mantan saudara raja dan wanita bajak laut mayat hidup yang menjalankan rumah bordil

Mengapa Ylvis menggunakan "go" alih-alih "say" dalam frasa seperti "Anjing pergi 'guk'"?

Bagaimana jika saya ketinggalan koneksi dan tidak punya uang untuk memesan penerbangan berikutnya?

Apakah pengurangan batas kecepatan akibat polusi juga berlaku untuk mobil listrik di Prancis?

Pilih semua kolom kecuali geometri menggunakan lapisan virtual

Apakah mungkin untuk secara metodis menemukan total cara membaca frasa yang diberikan membuat tumpukan?

Apakah semua kabel listrik dibuat sama?

Bagaimana senjata bisa dilawan dengan pertarungan jarak dekat tanpa kemampuan mentah atau penjelasan luar biasa?

Apakah simbol persentase konstan?

Apakah pesan terakhir Opportunity ke Earth "Baterai saya hampir habis dan hari mulai gelap"?

Kepadatan titik dari tabel XY di ArcGIS Pro

Akankah ArcGIS Pro memiliki ModelBuilder (lebih baik)?Menggunakan Iterator di ArcGIS Pro ModelBuilder?Mengekstrusi Tinggi di ModelBuilder dari ArcGIS Pro?Memperbarui peta ArcGIS Online dari ArcGIS ProArcGIS Pro ModelBuilder menghentikan eksekusi pada warningSpacing di tabel atribut ArcGIS Pro?Mengekstrak data baru dengan ArcGIS Pro ModelBuilder?ArcGIS Pro tidak memiliki catatan dalam tabel atribut dari SQL Server?Menggabungkan statistik zona dari beberapa raster menjadi satu tabel ArcGIS ProMenyembunyikan batas tampilan kepadatan titik di ArcGIS Pro?

Saya memiliki tabel CSV 1.1gb poin besar di Lat/Long WGS84 yang saya ingin menghasilkan raster kepadatan titik. Saya mengambil langkah-langkah berikut tetapi prosesnya sangat lambat. Saya menyambut saran tentang cara meningkatkan ini atau lebih memahami di mana sebenarnya hambatan saya.

Di pembuat model saya memiliki langkah-langkah berikut:


  1. Tabel XY ke Poin - output ke ruang kerja in_memory
  2. Proyek - proyeksikan ulang keluaran in_memory dari LatLong WGS84 ke WGS84 Web
    Mercator (Auxiliary Sphere) - di sini saya tidak dapat menyimpan hasilnya
    in_memory (tidak yakin mengapa) jadi saya menulisnya ke gdb
  3. Point Density - hitung raster kepadatan titik saya dari output
    2.

Saya menjalankan Windows Server 2019 dengan Intel Xeon Platinum 8175M CPU @ 2.5GHz (2 prosesor), RAM 760GB, OS 64-bit yang saya buat dengan instans AWS dengan hard disk 500GB.

Saat ini pembuat model saya berada di langkah 1 pada 1% dan sudah ada di sana untuk sementara waktu. Penggunaan CPU saya 1% dan Memori juga 1%

Bagaimana saya bisa mempercepat proses saya?

Saya memiliki tabel CSV 1.1gb poin besar di Lat/Long WGS84 yang saya ingin menghasilkan raster kepadatan titik. Saya mengambil langkah-langkah berikut tetapi prosesnya sangat lambat. Saya menyambut saran tentang cara meningkatkan ini atau lebih memahami di mana sebenarnya hambatan saya.

Di pembuat model saya memiliki langkah-langkah berikut:


  1. Tabel XY ke Poin - output ke ruang kerja in_memory
  2. Proyek - proyeksikan ulang keluaran in_memory dari LatLong WGS84 ke WGS84 Web
    Mercator (Auxiliary Sphere) - di sini saya tidak dapat menyimpan hasilnya
    in_memory (tidak yakin mengapa) jadi saya menulisnya ke gdb
  3. Point Density - hitung raster kepadatan titik saya dari output
    2.

Saya menjalankan Windows Server 2019 dengan Intel Xeon Platinum 8175M CPU @ 2.5GHz (2 prosesor), RAM 760GB, OS 64-bit yang saya buat dengan instans AWS dengan hard disk 500GB.

Saat ini pembuat model saya berada di langkah 1 pada 1% dan sudah ada di sana untuk sementara waktu. Penggunaan CPU saya 1% dan Memori juga 1%

Bagaimana saya bisa mempercepat proses saya?

Saya memiliki tabel CSV 1.1gb poin besar di Lat/Long WGS84 yang saya ingin menghasilkan raster kepadatan titik. Saya mengambil langkah-langkah berikut tetapi prosesnya sangat lambat. Saya menyambut saran tentang cara meningkatkan ini atau lebih memahami di mana sebenarnya hambatan saya.

Di pembuat model saya memiliki langkah-langkah berikut:


  1. Tabel XY ke Poin - output ke ruang kerja in_memory
  2. Proyek - proyeksikan ulang keluaran in_memory dari LatLong WGS84 ke WGS84 Web
    Mercator (Auxiliary Sphere) - di sini saya tidak dapat menyimpan hasilnya
    in_memory (tidak yakin mengapa) jadi saya menulisnya ke gdb
  3. Point Density - hitung raster kepadatan titik saya dari output
    2.

Saya menjalankan Windows Server 2019 dengan Intel Xeon Platinum 8175M CPU @ 2.5GHz (2 prosesor), RAM 760GB, OS 64-bit yang saya buat dengan instans AWS dengan hard disk 500GB.

Saat ini pembuat model saya berada di langkah 1 pada 1% dan sudah ada di sana untuk sementara waktu. Penggunaan CPU saya 1% dan Memori juga 1%

Bagaimana saya bisa mempercepat proses saya?

Saya memiliki tabel poin CSV 1.1gb besar di Lat/Long WGS84 yang saya ingin menghasilkan raster kepadatan titik. Saya mengambil langkah-langkah berikut tetapi prosesnya sangat lambat. Saya menyambut saran tentang cara meningkatkan ini atau lebih memahami di mana sebenarnya hambatan saya.

Di pembuat model saya memiliki langkah-langkah berikut:


  1. Tabel XY ke Poin - output ke ruang kerja in_memory
  2. Proyek - proyeksikan ulang keluaran in_memory dari LatLong WGS84 ke WGS84 Web
    Mercator (Auxiliary Sphere) - di sini saya tidak dapat menyimpan hasilnya
    in_memory (tidak yakin mengapa) jadi saya menulisnya ke gdb
  3. Point Density - hitung raster kepadatan titik saya dari output
    2.

Saya menjalankan Windows Server 2019 dengan Intel Xeon Platinum 8175M CPU @2.5GHz (2 prosesor), RAM 760GB, OS 64-bit yang saya buat dengan instans AWS dengan hard disk 500GB.

Saat ini pembuat model saya berada di langkah 1 pada 1% dan sudah ada di sana untuk sementara waktu. Penggunaan CPU saya 1% dan Memori juga 1%


2 Jawaban 2

Following @metasequoia's suggestion, I'm posting the workaround I found for this issue. This doesn't actually answer the question since it relies on rasterio instead of GDAL, but it might help other people.

This minimal example worked for me (with rasterio 1.0.1):

Thanks, that helped me out when looking for a solution to read tifs in google cloud storage from compute nodes (without data duplication).

– metasequoia
Dec 4 '18 at 0:17

Rasterio added support for gcs urls in version 1.0.15. If you're working with a recent version you can now use the following:

This is convenient, as you no longer need to wrap rasterio requests in the env.


Tonton videonya: GIS Tools: ArcMapArcGIS Tool for counting inside features. points