bagaimana saya bisa mengubah shapefile ke batas lat dan lon?
12
saya memiliki shapefile subdivisi negara dan saya ingin mengekstraksi array dari lat dan lon yang terikat untuk setiap divisi .. apakah mungkin untuk dilakukan?
Salah satu cara yang mungkin untuk melanjutkan menggunakan SAGA GIS http://www.saga-gis.org
Setelah membuka shapefile Anda, jalankan 3 modul ini: 1. Modul \ Bentuk \ Alat \ Dapatkan Bentuk
Modul \ Bentuk \ Alat \ Poin \ Poin dari garis [bertentangan dengan apa yang disarankan oleh namanya, Anda juga dapat menggunakan ini untuk mendapatkan poin dari poligon]
Modul \ Bentuk \ Alat \ Poin \ Tambahkan Koordinat ke poin Ini akan memberi Anda tabel yang berisi koordinat x dan y dari 4 sudut kotak pembatas file poligon Anda.
Di arcgis, ini kode python. hasilnya adalah daftar minx, miny, maxx, maxy, minM, maxM, minZ, maxZ (
import arcpy
for feat in arcpy.SearchCursor(r"c:\data\f.gdb\counties"):print feat.Shape.extent
-2.6685272725154649.4265363633626-2.5284818181812149.5079454546192NaNNaNNaNNaN-10.46333636378251.4455454544593-6.0130545458304555.3799909091533NaNNaNNaNNaN-4.7777818182761454.0555454544593-4.3534727268846854.4100000000002NaNNaNNaNNaN
Ini versi R, menggunakan contoh data dari paket rgdal:
library(rgdal)
dsn <- system.file("vectors/ps_cant_31.MIF",package="rgdal")[1]
d <- readOGR(dsn = dsn, layer="ps_cant_31")## transform if this is not longlatif(is.projected(d)) d <- spTransform(d, CRS("+proj=longlat +ellps=WGS84"))for(i in1:nrow(d)){print(bbox(d[i,]))}
Saya menggunakan fiona dan dibentuk untuk tugas semacam itu:
import fiona
from shapely.geometry import shape
with fiona.open(r'd:\Projects\_00_Data\_USstates\fe_2007_us_state00.shp','r')as features:for i, feat in enumerate(features):
geom = shape(feat['geometry'])
name = feat['properties']['NAME00']print','.join((name,)+ tuple([str(i)for i in geom.bounds]))
Outputnya terlihat seperti ini, yaitu koordinat lat / lon: -124.72583900000001,45.544321, -116.915989,49.002494 -82.626182,37.202467, -77.71951899999999999, -Umumnya-produk-sebaginya-memiliki-memilikiprodukuntukmemilikiprodukuntukmemilikiprodukuntukmemilikiprodukyangdimasukkandimulai-denganpeluang-peluang-kuat-kekuatan-kekuatan-saat-sebuahanda-dengan-satu-satu-pilihan-satu-pilihan-satu-pilihan-lebihbaik-memiliki-denganlebihbaik-kekuatan-lebihdaripada-adalahmembuat-lebihdaripada-memiliki-kekuatan-memiliki-kekuatanmemilikimemilikiprodukuntukmemiliki lebih banyakproduk—memiliki—————————————————Memilih-memiliki selulebihnya?
Jawaban:
Menggunakan modul ogr Python dari OSGEO, contoh ini akan memberi Anda sebuah tuple berisi coords yang menentukan sebuah amplop untuk setiap fitur.
sumber
Salah satu cara yang mungkin untuk melanjutkan menggunakan SAGA GIS http://www.saga-gis.org Setelah membuka shapefile Anda, jalankan 3 modul ini: 1. Modul \ Bentuk \ Alat \ Dapatkan Bentuk
Modul \ Bentuk \ Alat \ Poin \ Poin dari garis [bertentangan dengan apa yang disarankan oleh namanya, Anda juga dapat menggunakan ini untuk mendapatkan poin dari poligon]
Modul \ Bentuk \ Alat \ Poin \ Tambahkan Koordinat ke poin Ini akan memberi Anda tabel yang berisi koordinat x dan y dari 4 sudut kotak pembatas file poligon Anda.
sumber
Di arcgis, ini kode python. hasilnya adalah daftar minx, miny, maxx, maxy, minM, maxM, minZ, maxZ (
sumber
Ini versi R, menggunakan contoh data dari paket rgdal:
sumber
Saya menggunakan fiona dan dibentuk untuk tugas semacam itu:
sumber