Vairāk

Veidot punktu shapefile no platuma / garuma, izmantojot GeoTools?

Veidot punktu shapefile no platuma / garuma, izmantojot GeoTools?


Esmu jauns ģeotīklos, un man ir mazliet grūti uzzīmēt punktu, izmantojot platuma un garuma vērtības. vai kāds mani var pamudināt, lai atrisinātu manu problēmu.

GeometryFactory geometryFactory = JTSFactoryFinder.getGeometryFactory (): SimpleFeatureBuilder featureBuilder = jauns SimpleFeatureBuilder (TYPE); Punkta punkts = geometryFactory.createPoint (jauna koordināta (garums, platums)); featureBuilder.add (punkts);

taču šis kods nedarbojas pareizi. vai kāds varētu paskaidrot, kāpēc ...?


Kodā jums trūkst definēt objekta tipu: GeoTools (un ne tikai) vispirms jums būs jādefinē datu struktūra, ti Funkcijas tips (~ shēma). Tam jums ir jāizmantoSimpleFeatureTypeBuilderkas palīdzēs izveidot objekta tipu, pēc tam jūs varēsiet pareizi izveidotSimpleFeatureBuilderpamatojoties uz iepriekš definēto iezīmju tipu (VEIDSarguments jūsu kodā)

Izveidojiet vienkāršās funkcijas veidu (datu struktūra):

SimpleFeatureTypeBuilder b = jauns SimpleFeatureTypeBuilder (); // iestatiet nosaukumu b.setName ("MyFeatureType"); // pievienot ģeometrijas rekvizītu b.setCRS (DefaultGeographicCRS.WGS84); // iestatiet crs first b.add ("location", Point.class); // pēc tam pievienojiet ģeometriju // izveidojiet tipu final SimpleFeatureType TYPE = b.buildFeatureType ();

Izveidojiet objektu (-us), izmantojot iepriekš definēto objekta tipu

SimpleFeatureBuilder featureBuilder = jauns SimpleFeatureBuilder (TIPS); GeometryFactory geometryFactory = JTSFactoryFinder.getGeometryFactory (); Punkta punkts = geometryFactory.createPoint (jauna koordināta (garums, platums)); featureBuilder.add (punkts); Funkcija SimpleFeature = featureBuilder.buildFeature ("fid.1"); // izveidojiet 1. funkciju

Nākamajā solī varat izveidot objektu kolekciju, izmantojot iepriekš definēto (-ās) iezīmi (-es):

DefaultFeatureCollection featureCollection = jauns DefaultFeatureCollection ("iekšējais", TYPE); featureCollection.add (funkcija); // Pievienot iezīmi 1 // Pievienot iezīmi 2, 3 utt

Es iesaku jums sekot piemēriem no vietnes http://docs.geotools.org/latest/userguide/library/main. Tas arī palīdzēs saprast, kā ģeotelpas dati tiek strukturēti Geotools

Plānošanas daļā iesaku iepazīties ar ģeotehnikas rīku renderēšanas moduli: http://docs.geotools.org/latest/userguide/library/render/map.html.


Pienākumi

Berjēnas apgabala ĢIS un zemes aprakstīšanas departamentam ir divas galvenās atbildības jomas:

Ģeogrāfiskās informācijas sistēmas

Visaptverošu ģeogrāfiskās informācijas sistēmu (ĢIS) izstrāde, kas nodrošina vērtīgus instrumentus efektīvākai un efektīvākai piekļuvei, sasaistei, analīzei un informācijas uzturēšanai novadam un tā iedzīvotājiem.

Zemes apraksts

Īpašuma aprakstu, nodokļu maksātāju vārdu, adrešu ierakstu, īpašumu karšu, kā arī aerofotogrāfijas uzturēšana visām 39 vietējām vērtēšanas vienībām novadā.


Konvertējot punktu koordinātas


Objekts break.df ir tikai pamatkabelis (datu rāmis), nevis telpisks sf punkti objekts. Lai šos datus pārveidotu telpiskos punktos, mums tie būs “jāatrod” kartē. Par laimi, mūsu datiem ir informācija par garumu (x-koordinātu) un platumu (y-koordinātu), kur garums ir jūsu X koordināta un aptver austrumus / rietumus, bet platums ir jūsu Y koordināta un aptver ziemeļus / dienvidus, tāpēc mēs jau zinām, kur punkti atrodas uz Zemes, un jautājums ir tikai par to ievietošanu tur.

Lai izveidotu punktu, mēs izmantosim funkciju st_as_sf () sf objekts no break.df.

Arguments coords = norāda katra punkta (X, Y) koordinātas - vai garumu, platumu. Datu ietvarā break.df, šie mainīgie ir ērti nosaukti X un un tiek norādīti st_as_sf (), izmantojot argumentu coords kā coords = c (& quotX & quot, & quotY & quot).

Arguments crs = norāda koordinātu atskaites sistēmu (CRS). DRS attiecas uz veidu, kā telpiskie dati, kas attēlo zemes virsmu (kas ir apaļa un trīsdimensiju), tiek saplacināti, lai jūs varētu tos “uzzīmēt” uz divdimensiju virsmas. Šī koncepcija ir tik svarīga, ka tai ir jābūt savai sadaļai.

Koordinātu atskaites sistēma


Mēs neiziesim cauri īsts DRS uzgriežņi un bultskrūves, kuras jūs varat izlasīt Geocomputation with R (GWR) 2.4. un 6. nodaļā, taču mēs to pietiekami pārdzīvosim, lai jūs varētu tikt galā ar lielāko daļu ar DRS saistīto telpisko datu apkarošanas uzdevumu.

DRS satur divus galvenos komponentus: ģeogrāfisko koordinātu sistēmu (GCS) un projicēto koordinātu sistēmu (PCS). GCS var tālāk sadalīt divās daļās: elipse un atsauces punkts. Elipse ir Zemes formas modelis. Datums nosaka šī modeļa koordinātu sistēmu - sākuma punktu un asis. Jums ir nepieciešami šie divi pamatkomponenti, lai novietotu punktus uz Zemes trīsdimensiju virsmas. Padomājiet par to kā par mēģinājumu izveidot globusu (elipse) un izdomāt, kur novietot punktus uz šī globusa (nulles punkts).

Pēc tam PCS pārvērš šos punktus no globusa divdimensiju telpā. Mums tas jādara, jo, veidojot plakanu papīru vai ekrānā redzamas kartes, nevis globus (tas ir sava veida grūti pārnest pasauli apkārt, kad atrodaties pilsētā).

Izmantojot funkciju st_crs (), varat izsaukt telpisko datu kopas DRS. Piemēram, skaitīšanas trakta datu DRS ir

DRS aprakstīšanai ir divi veidi: epsg kods un proj4string. Epsg ir saīsnes veids, kā definēt visus DRS komponentus ar vienu numuru. Šeit varat meklēt visus iespējamos epsg.

Turpretī proj4string izsaka visu. Projicēto koordinātu sistēmu jeb PCS definē ar + proj =, šajā gadījumā longlat, jo mums ir X, Y vai garuma / platuma koordinātas. Kā es zināju, kuru nulles punktu un elipsi precizēt? Maģija. Faktiski pastāv kopīgi lietotu DRS / elipses kopums, no kuriem viens ir NAD83 / GRS80. Pārbaudiet izdales materiālu Canvas 6. nedēļas mapē, lai atrastu citas izplatītas kombinācijas.

Šajā diskusijā par DRS jums vajadzētu atcelt trīs svarīgas lietas. Pirmkārt, kā aprakstīts iepriekš, DRS ir divas sastāvdaļas, GCS un PCS, kas jums vienmēr jānorāda. Tas nozīmē, ka, ieviešot telpiskos datus, jums jāzina, kas ir GCS un PCS. Piemēram, jums saka vietne, kurā lejupielādējāt datus. Vai arī persona, kas jums sniedz datus, stāsta jums. Varbūt burvju feja to čukstēja tev ausī. Lai vai kā, jums vienkārši jāzina. Ja norādīsiet nepareizu DRS, radīsies dažas problēmas. Visiem datiem, kas ir (X, Y) formā, ir projicēta garuma / platuma koordinātu sistēma. Šajā gadījumā vienkārši norādiet DRS, kuru izmantojām iepriekš minētajā kodā, un jums vajadzētu būt labi.

Otrkārt, visiem jūsu pašreizējās R sesijas telpiskajiem datiem jābūt vienādam DRS. Ja tas nenotiek, jūs nevarēsit pārklāt objektus kartē vai veikt nevienu no daudzslāņu telpiskajām operācijām, kuras mums bija jāveic 5. laboratorijā. DRS mēs varam pārbaudīt, izmantojot funkciju st_crs (). Pārbaudīsim, vai sf. piesaista un pārtraukums.sf ir tāda pati DRS

Punktu pārprojektēšana


Trešā svarīgā lieta ir tā, ka DRS jābūt saderīgai ar funkcijām, kuras plānojat izmantot savos datos. Svarīga saderības problēma ir tā, ka jūsu telpiskajiem datiem ir jābūt atbilstošam attāluma metrikai, piemēram, metriem vai jūdzēm. Attālums garuma / platuma grādos ir decimālgrādos, ko (1) lielākā daļa cilvēku nesaprot un (2) nav saderīga ar vairākām telpiskām funkcijām.

Mēs varam uzzināt telpisko datu kopas vienības, izmantojot funkciju st_crs () un izsaucot vienības šādi

Pārprojektēsim DRS uz kaut ko tādu, kas aizņem metrus. Populāra skaitītāju bāzes projicēto koordinātu sistēma ir universālā šķērsvirziena Mercator (UTM). UTM atdala ASV atsevišķās zonās, un Ziemeļkalifornija atrodas 10. zonā, kā parādīts zemāk redzamajā attēlā.

Pārprojektēsim abus sf. piesaista un pārtraukumi.sf uz UTM 10. zonas projicēto koordinātu sistēmu. Izmantojiet + proj = utm kā PCS, NAD83 kā nulles punktu un GRS80 kā elipsi (populāras izvēles iespējas ASV projekcijai / nulles punktam / elipsei). Ikreiz, kad izmantojat UTM, jums arī jānorāda zona, ko mēs darām, izmantojot + zona = 10. Lai pārprojektētu, izmantojiet funkciju st_transform () šādi

“M” apzīmē metrus. Lieliski. Tagad kartēsim ielaušanās.

Vēl viena svarīga problēma, ar kuru varat saskarties, ir tāda, ka shapefile vai jebkura no avota lejupielādētā telpisko datu kopa nesatur DRS (neprojektēta vai nezināma). Šajā gadījumā izmantojiet funkciju st_set_crs (), lai iestatītu DRS. Sīkāku informāciju skatiet GWR 6.1.

  1. DRS vienmēr jānosaka jebkurai telpisko datu kopai, kuru ievedat R.
  2. Ja plānojat strādāt ar vairākām telpisko datu kopām vienā un tajā pašā projektā, pārliecinieties, vai tiem ir viens un tas pats DRS.
  3. Pārliecinieties, vai DRS ir piemērots tiem telpisko analīžu veidiem, kurus plānojat veikt.

Ievērojot šos principus, jums vajadzētu būt iespējai pārvarēt lielāko daļu ar DRS saistīto jautājumu. Ja jūs iestrēgstat, izlasiet GWR Ch. 2.4 un 6.

/>
Šis darbs ir licencēts saskaņā ar starptautisko Creative Commons Attribution-NonCommercial 4.0 licenci.


Ģeogrāfiskās informācijas sistēmu (GIS) datu zinātnes paātrināšana ar RAPIDS cuSpatial un GPU

Ģeogrāfiskās informācijas sistēmu (ĢIS) ekonomiskā ietekme kopš 2013. gada ir dubultojusies. Pilsoniskie dati, informācija par pandēmiju un autonoma informācija par transportlīdzekļiem rada plašas datu bāzes, kurām nepieciešama izpratne. Mūsu mērķis NVIDIA ir sniegt GPU paātrinātus risinājumus visbiežāk izmantotajiem ĢIS datu zinātnes aspektiem. Pēdējos deviņos ...


ESRI formas failu sagatavošana pārveidošanai tiešsaistē

ESRI shapefile ir vektoru datu glabāšanas formāts, kas sastāv no trim failiem ar vienādu faila nosaukuma prefiksu:

  • shp: ģeometrija tiek glabāta kā virsotņu saraksts
  • shx: ģeometrijas indekss, lai nodrošinātu ātru, uz indeksu balstītu meklēšanu
  • dbf: dBase fails, ko izmanto atribūtu informācijas glabāšanai

Shapefiles var būt papildu faili, piemēram, .prj, kurā glabājas koordinātu sistēmas informācija (to izmanto ArcGIS), vai .cpg, kas norāda kodu lapu izmantojamo rakstzīmju kopas identificēšanai. Pēc definīcijas vektors ir uz Dekarta balstīta (t.i. X, Y) datu struktūra, ko izmanto telpisko datu glabāšanai. Lielākā daļa ģeogrāfiskās informācijas sistēmu, tostarp ESRI ArcGIS, ģeotelpiskos datus glabā vektoru formātā. Citi izplatītie telpisko datu formāti ietver rastru un TIN. Lai iegūtu vairāk informācijas, izlasiet ESRI Shapefile tehnisko aprakstu: ESRI Baltā grāmata, 1988. gada jūlijs

ESRI prj un WKT faili

Shapefilei parasti pievieno .prj failu, kas satur labi zināmu teksta (WKT) virkni, kurā glabājas koordinātu atskaites sistēmas informācija. Labi zināms teksts (WKT) ir teksta iezīmēšanas valoda vektoru ģeometrijas objektu attēlošanai kartē. Formātu sākotnēji definēja Atvērtais ģeotelpiskais konsorcijs (OGC) un aprakstīja viņu vienkāršās funkcijas piekļuves specifikācijā, kurā ESRI bija līdzdalībnieks. Pašreizējā WKT standarta versija tika publicēta 2019. gada 13. augustā.

Ne visos .prj failos ir vienādi parametri. Starp ESRI WKT failiem un OGC WKT failiem ir arī sintaktiskas atšķirības. Apmeklējiet EPSG.io, lai apskatītu dažādu WKT failu piemērus. Šeit ir EPSG WKT failu piemēri: 25832:

Labi pazīstama teksta kā HTML piemērs

OGC WKT piemērs

ESRI WKT piemērs

Datuma transformācijas

Datuma transformācijas ir nepieciešamas, ja transformācija no vienas koordinātu sistēmas uz citu ietver ģeogrāfisko koordinātu sistēmu. Datuma pārveidojumi jāveic piesardzīgi, jo tie var izraisīt ievērojamas datu nobīdes par 1 metru vai vairāk. Ir vairākas atšķirīgas metodes, lai veiktu nulles punkta pārveidošanu - Bursa-Volfs ir viena no metodēm.

WKT parametrs TOWGS84 tiek izmantots, lai tuvinātu transformāciju no horizontālā punkta uz WGS84 nulles punktu. Sākotnējā WKT vienkāršo funkciju specifikācijā TOWGS84 nav norādīts kā derīgs atslēgvārds. ESRI neatbalsta parametru TOWGS84 WKT. Tā vietā, lai izmantotu fiksētu transformāciju, ESRI programmatūra lūdz lietotāju pēc vajadzības izvēlēties piemērotu transformācijas metodi, tāpat kā QGIS.

Pašreizējai WKT versijai nav atgriezeniskas savietojamības ar TOWGS84. Ģeodēzisko datu WKT apraksti, kas rakstīti OGC 01-009 specifikācijā, nav lasāmi, ieviešot pašreizējo WKT specifikāciju, ja ir paredzēts papildu TOWGS84 objekts.

Koordinātu atskaites sistēmas identifikācija hale »savienot

hale »connect telpisko datu apstrādei izmanto atvērtā koda Java ĢIS rīku komplektu GeoTools. Savukārt GeoTools satur EPSG datu bāzes kopiju. Kad platformā tiek augšupielādēts shapefile, hale »connect izmanto GeoTools bibliotēku, lai mēģinātu saskaņot .prj failā esošo WKT ar DRS definīcijām EPSG datu bāzē. Ja tiek atrasta atbilstība, dati pēc tam tiek apstrādāti, izmantojot saskaņoto EPSG CRS kodu.

Dažos gadījumos GeoTools nevar noteikt DRS definīciju, pamatojoties uz WKT ietverto informāciju. Lai atbalstītu GeoTools pareizas atbilstības atrašanā, lietotājs var pievienot parametru AUTHORITY kā pēdējo parametru WKT failam. Pēc tam GeoTools izmanto lietotāja ievadīto EPSG kodu AUTHORITY parametrā, lai formas failam piemērotu EPSG datu bāzē definēto DRS definīciju. Šīs definīcijas satur arī Bursa-Wolf parametrus, kas nepieciešami nulles punkta pārveidošanai, no DRS definīcijām EPSG datu bāzē.

Manuāli pievienojot parametru AUTHORITY ESRI prj failam

Varat manuāli rediģēt ESRI .prj failu, lai iekļautu papildu parametru AUTHORITY. Parametrs failā jāpievieno kā pēdējais beigu parametrs, izmantojot šādu modeli: AUTHORITY ["EPSG", "25832"]. Lai gan parametra AUTHORITY izmantošana pašreizējā WKT versijā ir oficiāli novecojusi, atpakaļsaderības nolūkos tā izmantošana joprojām tiek pieņemta.

Pašreizējā WKT versijā objekts AUTHORITY ir aizstāts ar ID objektu. Identifikācijas objekts nav tik šauri definēts kā AUTHORITY bija iepriekšējās versijās.

Asis uzsist

Galvenais iemesls, kāpēc lietotāji saskaras ar asu pagriešanas problēmām, ir tāds, ka formas failā ģeometriju asu secība ir garums, platums / austrumu virziens, ziemeļdaļa, bet WKT definīcija .prj failā apraksta DRS, kurai ir platums, garums / Northing, austrumu ass kārtība. Piemēram, WGS 84 = EPSG: 4326.

Tas ir problemātiski, ja nepastāv EPSG definīcija DRS ar formas secībā izmantoto asu secību, jo parametram AUTHORITY nevar nodrošināt atbilstošu kodu. Šādos gadījumos ir nepieciešams izmantot slēdzi “Flip koordinātu secība”, lai pēc datu kopas publicēšanas mainītu asu secību.

Galu galā vienīgais veids, kā apstrādāt asu secību, ir analizēt shapefile ģeometriju asu secību un noteikt, vai asu secība atbilst DRS definīcijai .prj failā, un rīkoties atbilstoši, ti, vai nu pievienot parametram AUTHORITY .prj failu vai izmantojiet pārslēgu “Pārvērst koordinātu secību” Hale »savienojumā.

& laquo Profila izveide, izmantojot hale »connect Profile Management Tool (PMT)


Kartes displeju izveidošana ar datiem projicēto koordinātu atskaites sistēmā

Šis piemērs parāda, kā projektētajā koordinātu atskaites sistēmā importēt un attēlot ģeogrāfiskos datus, kas satur koordinātas.

Šis piemērs jo īpaši parāda, kā to izdarīt

Importēt īpašas rastra un vektoru datu kopas

Izveidojiet karšu displejus datu vizualizēšanai

Parādīt vairākas datu kopas kartes displejā

Vienā kartes attēlojumā parādiet vairākas datu kopas ar koordinātām ģeogrāfiskās un projicētās koordinātu atskaites sistēmās

1. piemērs: Rastra datu importēšana projicēto koordinātu atskaites sistēmā

Ģeogrāfiskos rastra datus, kas projektētajā koordinātu atskaites sistēmā satur koordinātas, var saglabāt dažādos formātos, ieskaitot standarta failu formātus, piemēram, GeoTIFF, Spatial Data Transfer Standard (SDTS), NetCDF, HDF4 vai HDF5. Šis piemērs ilustrē datu importēšanu no GeoTIFF faila. Faila dati satur koordinātas projektētās kartes koordinātu atskaites sistēmā Masačūsetsas štata lidostas kontinentālās zonas koordinātu sistēma .

Attēla koordinātas GeoTIFF failā boston.tif atrodas projicētā koordinātu atskaites sistēmā. To var noteikt, izmantojot ģeotiffinfo funkciju, un pārbaudiet lauka PCS un Projection vērtības.

Koordinātu garuma vienību nosaka informācijas struktūras lauks UOMLength.

Lai importētu attēlu un telpiskās atsauces objektu, izmantojiet readgeoraster.

2. piemērs: Rastra datu parādīšana projicēto koordinātu atskaites sistēmā

Attēlu var parādīt uz parastām MATLAB asīm, izmantojot mapshow, kas parāda attēlu un asu robežas nosaka atsauces objekta R noteiktajām robežām. Kā minēts iepriekš, koordinātas atrodas ASV apsekojuma pēdās un ir saistītas ar izcelsmi uz kartes dienvidrietumiem, tāpēc skaitļi ir lieli. Zonas ietvaros koordinātas vienmēr ir pozitīvas.

3. piemērs: vektoru datu importēšana projicēto koordinātu atskaites sistēmā

Ģeogrāfisko vektoru datus, kas projektētās koordinātu atskaites sistēmā satur koordinātas, var uzglabāt formas failos. Šis piemērs parāda, kā vektoru datus importēt projicētā koordinātu atskaites sistēmā no shapefile, boston_roads.shp.

Importējiet vektora līnijas datus no faila boston_roads.shp.

Lai iegūtu informāciju par projicēto koordinātu atskaites sistēmu, vispirms atgrieziet informāciju par shapefile kā struktūru. Pēc tam vaicājiet laukam CoordinateReferenceSystem.

4. piemērs: Parādīt vektoru un rastra datus projicēto koordinātu atskaites sistēmā

Vektora un rastra dati šajā piemērā ir vienā projektētajā koordinātu atskaites sistēmā. Tomēr vektoru dati ir metru garuma vienībās, bet rastra dati - apsekojuma pēdas garuma vienībās. Pārveidojiet rastra datus metru garuma vienībās un parādiet datus uz tām pašām asīm.

Konvertējiet rastra attēla koordinātas no ASV apsekojuma pēdas vienībām uz metru.

Parādiet rastra attēlu un vektoru datus, izmantojot mapshow.

5. piemērs: Datu parādīšana gan ģeogrāfiskās, gan projicētās koordinātu atskaites sistēmās

Jums var būt ģeogrāfiski dati, kuru koordinātas ir platumā un garumā, un citi dati, kuru koordinātas atrodas projicētajā koordinātu atskaites sistēmā. Šīs datu kopas varat parādīt tajā pašā kartes displejā. Šis piemērs parāda, kā attēlot datus ģeogrāfisko koordinātu atskaites sistēmā (platums un garums) ar datiem projicētajā kartes koordinātu atskaites sistēmā (Masačūsetsas štata plaknes kontinentālās zonas koordinātu sistēmā).

Lasiet rastra attēlu ar pasaules failu, kura koordinātas ir platumā un garumā. Izmantojiet imread, lai lasītu attēlu, un worldfileread, lai lasītu pasaules failu un izveidotu telpiskas atsauces objektu.

Lai tajā pašā kartes displejā parādītu pārskata attēlu un GeoTIFF attēlu, jums jāizveido kartes displejs, izmantojot Mapping Toolbox & # 8482 projekcijas struktūru, kas satur projekcijas informāciju par datiem projektētajā koordinātu atskaites sistēmā, Masačūsetsas štata lidostas kontinentālās zonas koordinātu sistēma . Lai šajā sistēmā izveidotu kartes attēlojumu, varat izmantot projekcijas informāciju, kas atrodas GeoTIFF failā. Izmantojiet funkciju geotiff2mstruct, lai izveidotu Mapping Toolbox & # 8482 projekcijas struktūru no GeoTIFF informācijas struktūras satura. Funkcija geotiff2mstruct atgriež projekciju metru vienībās. Izmantojiet projekcijas struktūru, lai definētu projekcijas parametrus kartes attēlojumam.

Izmantojiet Bostonas pārskata attēla platuma un garuma robežas.

Izmantojot kartes projekcijas struktūrā saglabāto projekcijas informāciju, izveidojiet kartes displeju un iestatiet kartes platuma un garuma robežas. Parādiet ģeogrāfiskos datus kartes asīs. geoshow projicē platuma un garuma koordinātas.

Tā kā GeoTIFF attēla koordinātas atrodas projicētā koordinātu atskaites sistēmā, izmantojiet mapshow, lai displejā pārklātu detalizētāku Bostonas attēlu. Sarkanā krāsā uzzīmējiet Bostonas attēla robežas.

Tuviniet GeoTIFF attēla ģeogrāfisko reģionu, nosakot asu robežas Bostonas attēla robežām un pievienojiet nelielu buferi. Ņemiet vērā, ka bufera lielums (delta) ir izteikts metros.

Jūs varat pārklāt ceļa vektorus uz kartes displeja. Izmantojiet simbolu specifikāciju, lai katrai ceļa klasei piešķirtu savu krāsu.

Varat arī pārklāt datus no GPS, kas saglabāts GPX failā. Importējiet ģeogrāfisko vektoru datus no faila boston_placenames.gpx, kas iekļauts programmatūrā Mapping Toolbox & # 8482. Fails satur ģeogrāfisko punktu pazīmju platuma un garuma koordinātas Bostonā, Masačūsetsā, ASV. Izmantojiet gpxread, lai lasītu GPX failu un atgrieztu ģeopunkta vektoru.

Pārklājiet vietvārdus uz kartes un palieliniet marķiera izmēru, nomainiet marķierus uz apļiem un iestatiet to malas un sejas krāsas dzeltenā krāsā.

Datu kopas informācija

Failos boston.tif un boston_ovr.jpg ir iekļauti materiāli, kuru autortiesības ir aizsargātas ar GeoEye, visas tiesības aizsargātas. GeoEye tika apvienots ar DigitalGlobe korporāciju 2013. gada 29. janvārī. Lai iegūtu papildinformāciju par datu kopām, izmantojiet komandas tipa boston.txt un tipa boston_ovr.txt.

Faili boston_roads.shp un boston_placenames.gpx ir no Ģeogrāfiskās informācijas biroja (MassGIS), Masačūsetsas Sadraudzības, Tehnoloģiju un drošības dienestu izpildbiroja. Lai iegūtu papildinformāciju par datu kopām, izmantojiet komandas tipa boston_roads.txt un ierakstiet boston_placenames_gpx.txt.


Kartes displeju izveidošana ar datiem projicēto koordinātu atskaites sistēmā

Šis piemērs parāda, kā projektētajā koordinātu atskaites sistēmā importēt un attēlot ģeogrāfiskos datus, kas satur koordinātas.

Šis piemērs jo īpaši parāda, kā to izdarīt

Importēt īpašas rastra un vektoru datu kopas

Izveidojiet karšu displejus datu vizualizēšanai

Parādīt vairākas datu kopas kartes displejā

Vienā kartes attēlojumā parādiet vairākas datu kopas ar koordinātām ģeogrāfiskās un projicētās koordinātu atskaites sistēmās

1. piemērs: Rastra datu importēšana projicēto koordinātu atskaites sistēmā

Ģeogrāfiskos rastra datus, kas projektētajā koordinātu atskaites sistēmā satur koordinātas, var saglabāt dažādos formātos, ieskaitot standarta failu formātus, piemēram, GeoTIFF, Spatial Data Transfer Standard (SDTS), NetCDF, HDF4 vai HDF5. Šis piemērs ilustrē datu importēšanu no GeoTIFF faila. Faila dati satur koordinātas projektētās kartes koordinātu atskaites sistēmā Masačūsetsas štata lidostas kontinentālās zonas koordinātu sistēma .

Attēla koordinātas GeoTIFF failā boston.tif atrodas projicētā koordinātu atskaites sistēmā. To var noteikt, izmantojot ģeotiffinfo funkciju, un pārbaudiet lauka PCS un Projection vērtības.

Koordinātu garuma vienību nosaka informācijas struktūras lauks UOMLength.

Lai importētu attēlu un telpiskās atsauces objektu, izmantojiet readgeoraster.

2. piemērs: Rastra datu parādīšana projicēto koordinātu atskaites sistēmā

Attēlu var parādīt uz parastām MATLAB asīm, izmantojot mapshow, kas parāda attēlu un asu robežas nosaka atsauces objekta R noteiktajām robežām. Kā minēts iepriekš, koordinātas atrodas ASV apsekojuma pēdās un ir saistītas ar izcelsmi uz kartes dienvidrietumiem, tāpēc skaitļi ir lieli. Zonas ietvaros koordinātas vienmēr ir pozitīvas.

3. piemērs: vektoru datu importēšana projicēto koordinātu atskaites sistēmā

Ģeogrāfisko vektoru datus, kas projektētās koordinātu atskaites sistēmā satur koordinātas, var uzglabāt formas failos. Šis piemērs parāda, kā vektoru datus importēt projicētā koordinātu atskaites sistēmā no shapefile, boston_roads.shp.

Importējiet vektora līnijas datus no faila boston_roads.shp.

Lai iegūtu informāciju par projicēto koordinātu atskaites sistēmu, vispirms atgrieziet informāciju par shapefile kā struktūru. Pēc tam vaicājiet laukam CoordinateReferenceSystem.

4. piemērs: Parādīt vektoru un rastra datus projicēto koordinātu atskaites sistēmā

Vektora un rastra dati šajā piemērā ir vienā projektētajā koordinātu atskaites sistēmā. Tomēr vektoru dati ir metru garuma vienībās, bet rastra dati - apsekojuma pēdas garuma vienībās. Pārveidojiet rastra datus metru garuma vienībās un parādiet datus uz tām pašām asīm.

Konvertējiet rastra attēla koordinātas no ASV apsekojuma pēdas vienībām uz metru.

Parādiet rastra attēlu un vektoru datus, izmantojot mapshow.

5. piemērs: Datu parādīšana gan ģeogrāfiskās, gan projicētās koordinātu atskaites sistēmās

Jums var būt ģeogrāfiski dati, kuru koordinātas ir platumā un garumā, un citi dati, kuru koordinātas atrodas projicētajā koordinātu atskaites sistēmā. Šīs datu kopas varat parādīt tajā pašā kartes displejā. Šis piemērs parāda, kā attēlot datus ģeogrāfisko koordinātu atskaites sistēmā (platums un garums) ar datiem projicētajā kartes koordinātu atskaites sistēmā (Masačūsetsas štata plaknes kontinentālās zonas koordinātu sistēmā).

Lasiet rastra attēlu ar pasaules failu, kura koordinātas ir platumā un garumā. Izmantojiet imread, lai lasītu attēlu, un worldfileread, lai lasītu pasaules failu un izveidotu telpiskas atsauces objektu.

Lai parādītu pārskata attēlu un GeoTIFF attēlu vienā un tajā pašā kartes displejā, jums jāizveido kartes displejs, izmantojot Mapping Toolbox & # 8482 projekcijas struktūru, kas satur projekcijas informāciju par datiem projektētajā koordinātu atskaites sistēmā, Masačūsetsas štata lidostas kontinentālās zonas koordinātu sistēma . Lai šajā sistēmā izveidotu kartes attēlojumu, varat izmantot projekcijas informāciju, kas atrodas GeoTIFF failā. Izmantojiet funkciju geotiff2mstruct, lai izveidotu Mapping Toolbox & # 8482 projekcijas struktūru no GeoTIFF informācijas struktūras satura. Funkcija geotiff2mstruct atgriež projekciju metru vienībās. Izmantojiet projekcijas struktūru, lai definētu projekcijas parametrus kartes attēlojumam.

Izmantojiet Bostonas pārskata attēla platuma un garuma robežas.

Izmantojot kartes projekcijas struktūrā saglabāto projekcijas informāciju, izveidojiet kartes displeju un iestatiet kartes platuma un garuma robežas. Parādiet ģeogrāfiskos datus kartes asīs. geoshow projicē platuma un garuma koordinātas.

Tā kā GeoTIFF attēla koordinātas atrodas projicētā koordinātu atskaites sistēmā, izmantojiet mapshow, lai displejā pārklātu detalizētāku Bostonas attēlu. Sarkanā krāsā uzzīmējiet Bostonas attēla robežas.

Tuviniet GeoTIFF attēla ģeogrāfisko reģionu, nosakot asu robežas Bostonas attēla robežām un pievienojiet nelielu buferi. Ņemiet vērā, ka bufera lielums (delta) ir izteikts metros.

Jūs varat pārklāt ceļa vektorus uz kartes displeja. Izmantojiet simbolu specifikāciju, lai katrai ceļa klasei piešķirtu savu krāsu.

Varat arī pārklāt datus no GPS, kas saglabāts GPX failā. Importējiet ģeogrāfisko vektoru datus no faila boston_placenames.gpx, kas iekļauts programmatūrā Mapping Toolbox & # 8482. Fails satur ģeogrāfisko punktu pazīmju platuma un garuma koordinātas Bostonā, Masačūsetsā, ASV. Izmantojiet gpxread, lai lasītu GPX failu un atgrieztu ģeopunkta vektoru.

Pārklājiet vietvārdus uz kartes un palieliniet marķiera izmēru, nomainiet marķierus uz apļiem un iestatiet to malas un sejas krāsas dzeltenā krāsā.


Sižetu veidošana

R ir nepieciešami vairāki paketes, lai lasītu Shapefiles un pārveidotu koordinātas, lai tās atbilstu jūsu zemes gabalu kartes projekcijām, vai otrādi. Šīs paketes ietvēra: rgdal, rgeos, sp un ggplot2. Lai sasniegtu līdzīgus rezultātus, var izmantot dažādas citas paketes.

Pirmkārt, shapefile tiek importēts R, izmantojot “readOGR” metodi, un pēc tam projekcija jāpārveido tā, lai tā atbilstu platuma un garuma koordinātu formātam, kāds man ir manā kūrētajā noziegumu datu kopā. Tas izrādījās visgrūtākais šī projekta posms, un, lai to izstrādātu, bija nepieciešams liels daudzums Google.

Kad tas ir izdarīts, no šejienes ir viegli doties, it īpaši, ja esat iepazinies ar ggplot2. Es to paņēmu pa solim un vispirms gribēju redzēt, kā iznāca kartes gabals. Tas tiek darīts, izmantojot geg_polygon metodi, kas pieejama ar ggplot2, un pārveidotos shapefile datus no augšas:

Nākamais solis bija pievienot visu noziegumu vietu slāni un nošķirt nozieguma veidus:

Teorētiski tas bija šī projekta mērķis, jo no iepriekš minētā sižeta ir skaidrs, ka tas nav ideāls, jo tas ir pārāk aizņemts un grūti interpretējams.

Es nolēmu apskatīt tikai viena veida noziegumus konkrētajā mēnesī. Redzot, ka man patīk braukt ar velosipēdu un es nekad negribētu, lai mans kalnu velosipēds tiktu nozagts, es vēlētos uzzināt, vai man vajadzētu būt īpaši modram, novietojot savu velosipēdu noteiktās vietās, tāpēc es apgriezu datu kopu, lai apskatītu velosipēdu zādzības mēnesī Decembris.

Tas radīja vizuāli pievilcīgāku sižetu, ko varēja interpretēt. Kā jau bija paredzēts, noziegumu gadījumi vairāk koncentrējas uz pilsētas centru.

Negatīvais ir tas, ka gūtā ieskata apjoms ir ierobežots ar viena veida noziegumiem uz vienu mēnesi gadā, savukārt ir vēl 13 noziegumu veidi un pārējie gada mēneši.


Vai konvertēt ģeokodus platumā / garumā NAD1983 uz UTM_1983_Zone_18N?

Es vēlos redzēt ģeokodēšanas rezultātu atšķirību apmēru, ko veikuši ArcGis un Texas A & ampM Geocoding Services. Man ir 500 ģeokodētas adreses, izmantojot Texas A & ampM ģeokodēšanas pakalpojumus Latitude / Longitude NAD1983. Tagad es vēlos to konvertēt uz (UTM_1983_Zone_18N). Lai es varētu aprēķināt attālumu starp diviem UTM, izmantojot Excel. Otrās ģeokodētās adreses mapē (UTM_1983_Zone_18N) tika izveidotas ar ArcGis 10.2.

Vai ir kāda ideja, kā veikt šo pārveidošanu?

by DanPatterson_Re noguris

Es joprojām esmu neizpratnē par to, kas jums ir. Ja koordinātu sistēmu īpašības ir pareizas (GCS NAD83 jeb decimāldaļu fails ar NAD83 nulles punktu), esat gatavs to projektēt, kā aprakstīts. Ja pārbaudījāt tās koordinātu sistēmu un tā nav zināma, tad izmantojiet rīku Definēt projekciju Arctoolbox, lai to definētu kā GCS NAD83 formātu, pēc tam izmantojiet rīku Projekts tajā pašā rīkjoslā. Ja tas ir kāds cits faila formāts, jums tas būs jāprecizē


Kartes displeju izveidošana ar datiem projicēto koordinātu atskaites sistēmā

Šis piemērs parāda, kā projektētajā koordinātu atskaites sistēmā importēt un attēlot ģeogrāfiskos datus, kas satur koordinātas.

Šis piemērs jo īpaši parāda, kā to izdarīt

Importēt īpašas rastra un vektoru datu kopas

Izveidojiet karšu displejus datu vizualizēšanai

Parādīt vairākas datu kopas kartes displejā

Vienā kartes attēlojumā parādiet vairākas datu kopas ar koordinātām ģeogrāfiskās un projicētās koordinātu atskaites sistēmās

1. piemērs: Rastra datu importēšana projicēto koordinātu atskaites sistēmā

Ģeogrāfiskos rastra datus, kas satur koordinātas projicētā koordinātu atskaites sistēmā, var saglabāt dažādos formātos, tostarp standarta failu formātos, piemēram, GeoTIFF, Spatial Data Transfer Standard (SDTS), NetCDF, HDF4 vai HDF5. Šis piemērs ilustrē datu importēšanu no GeoTIFF faila. Faila dati satur koordinātas projektētās kartes koordinātu atskaites sistēmā Massachusetts State Plane Mainland Zone coordinate system .

The coordinates of the image in the GeoTIFF file, boston.tif , are in a projected coordinate reference system. You can determine that by using the geotiffinfo function and examine the PCS and Projection field values.

The length unit of the coordinates are defined by the UOMLength field in the info structure.

To import the image and the spatial referencing object, use readgeoraster .

Example 2: Display Raster Data in Projected Coordinate Reference System

You can display the image on a regular MATLAB axes using mapshow , which displays the image and sets the axes limits to the limits defined by the referencing object, R . The coordinates, as mentioned above, are in US survey foot and are relative to an origin to the southwest of the map, which is why the numbers are large. The coordinates are always positive within the zone.

Example 3: Import Vector Data in Projected Coordinate Reference System

Geographic vector data that contains coordinates in a projected coordinate reference system can be stored in shapefiles. This example illustrates how to import vector data in a projected coordinate reference system from the shapefile, boston_roads.shp .

Import vector line data from the boston_roads.shp file.

To get information about the projected coordinate reference system, first return information about the shapefile as a structure. Then, query the CoordinateReferenceSystem field.

Example 4: Display Vector and Raster Data in Projected Coordinate Reference System

The vector and raster data in this example are in the same projected coordinate reference system. However, the vector data is in length units of meter, while the raster data is in length unit of survey foot. Convert the raster data to length units of meter and display the data on the same axes.

Convert the coordinates of the raster image from units of US survey foot to meter.

Display the raster image and vector data using mapshow .

Example 5: Display Data in both Geographic and Projected Coordinate Reference Systems

You may have geographic data whose coordinates are in latitude and longitude and other data whose coordinates are in a projected coordinate reference system. You can display these data sets in the same map display. This example illustrates how to display data in a geographic coordinate reference system (latitude and longitude) with data in a projected map coordinate reference system (Massachusetts State Plane Mainland Zone coordinate system).

Read a raster image with a worldfile whose coordinates are in latitude and longitude. Use imread to read the image and worldfileread to read the worldfile and construct a spatial referencing object.

To display the overview image and the GeoTIFF image in the same map display, you need to create a map display using a Mapping Toolbox™ projection structure containing the projection information for the data in the projected coordinate reference system, Massachusetts State Plane Mainland Zone coordinate system . To make a map display in this system, you can use the projection information contained in the GeoTIFF file. Use the geotiff2mstruct function to construct a Mapping Toolbox™ projection structure, from the contents of the GeoTIFF information structure. The geotiff2mstruct function returns a projection in units of meters. Use the projection structure to define the projection parameters for the map display.

Use the latitude and longitude limits of the Boston overview image.

Create a map display by using the projection information stored in the map projection structure and set the map latitude and longitude limits. Display the geographic data in the map axes. geoshow projects the latitude and longitude coordinates.

Since the coordinates of the GeoTIFF image are in a projected coordinate reference system, use mapshow to overlay the more detailed Boston image onto the display. Plot the boundaries of the Boston image in red.

Zoom to the geographic region of the GeoTIFF image by setting the axes limits to the limits of the Boston image and add a small buffer. Note that the buffer size ( delta ) is expressed in meters.

You can overlay the road vectors onto the map display. Use a symbol specification to give each class of road its own color.

You can also overlay data from a GPS stored in a GPX file. Import point geographic vector data from the boston_placenames.gpx file included with the Mapping Toolbox™ software. The file contains latitude and longitude coordinates of geographic point features in part of Boston, Massachusetts, USA. Use gpxread to read the GPX file and return a geopoint vector .

Overlay the placenames onto the map and increase the marker size, change the markers to circles and set their edge and face colors to yellow.

Data Set Information

The files boston.tif and boston_ovr.jpg include materials copyrighted by GeoEye, all rights reserved. GeoEye was merged into the DigitalGlobe corporation on January 29th, 2013. For more information about the data sets, use the commands type boston.txt and type boston_ovr.txt .

The files boston_roads.shp and boston_placenames.gpx are from the Bureau of Geographic Information (MassGIS), Commonwealth of Massachusetts, Executive Office of Technology and Security Services. For more information about the data sets, use the commands type boston_roads.txt and type boston_placenames_gpx.txt .


Skatīties video: How to make Attribute for shapefiles in ArcGIS - ArcGIS for Beginners